diff --git a/Projects/Strasse/Analysis.cxx b/Projects/Strasse/Analysis.cxx
index 229b6e1e678bb690d1ce700c2967275058f0fbbf..6d3dabf82b2c59536d087de5223ea57b3685bc61 100644
--- a/Projects/Strasse/Analysis.cxx
+++ b/Projects/Strasse/Analysis.cxx
@@ -77,12 +77,13 @@ void Analysis::TreatEvent(){
 
     double deltaPhi = abs(Proton1.Phi()/deg-Proton2.Phi()/deg);
     double sumTheta = Proton1.Theta()/deg+Proton2.Theta()/deg;
-    double OpeningAngle = Proton1.Angle(Proton2)/deg;
+    Theta12  = Proton1.Angle(Proton2)/deg;
+
     // reject event that make no physical sense
-    if(deltaPhi<170 && sumTheta<80){
+    /*if(deltaPhi<170 && sumTheta<80){
       return;
       }
-    
+    */
     // computing minimum distance of the two lines
     TVector3 Vertex;
     TVector3 delta;
@@ -127,7 +128,7 @@ void Analysis::End(){
 void Analysis::InitOutputBranch() {
   RootOutput::getInstance()->GetTree()->Branch("Ex",&Ex,"Ex/D");
   RootOutput::getInstance()->GetTree()->Branch("ELab",&ELab,"ELab/D");
-  RootOutput::getInstance()->GetTree()->Branch("ThetaLab",&ThetaLab,"ThetaLab/D");
+  RootOutput::getInstance()->GetTree()->Branch("Theta12",&Theta12,"Theta12/D");
   RootOutput::getInstance()->GetTree()->Branch("ThetaCM",&ThetaCM,"ThetaCM/D");
   RootOutput::getInstance()->GetTree()->Branch("VertexX",&VertexX,"VertexX/D");
   RootOutput::getInstance()->GetTree()->Branch("VertexY",&VertexY,"VertexY/D");
@@ -150,7 +151,7 @@ void Analysis::InitInputBranch(){
 void Analysis::ReInitValue(){
   Ex = -1000 ;
   ELab = -1000;
-  ThetaLab = -1000;
+  Theta12 = -1000;
   ThetaCM = -1000;
   VertexX=-1000;
   VertexY=-1000;
@@ -158,7 +159,6 @@ void Analysis::ReInitValue(){
   deltaX=-1000;
   deltaY=-1000;
   deltaZ=-1000;
-
   Distance=-1000;
 }
 
diff --git a/Projects/Strasse/Analysis.h b/Projects/Strasse/Analysis.h
index c1e0c3f3801b287aed3159d887877a1ff6973a2b..c6561c0a487cfed78aba2d304e30cf8c8ce23987 100644
--- a/Projects/Strasse/Analysis.h
+++ b/Projects/Strasse/Analysis.h
@@ -53,7 +53,7 @@ class Analysis: public NPL::VAnalysis{
   private:
     double Ex;
     double ELab;
-    double ThetaLab;
+    double Theta12;
     double ThetaCM;
     double VertexX;
     double VertexY;
diff --git a/Projects/Strasse/macro/Control.cxx b/Projects/Strasse/macro/Control.cxx
index 2627945ebed43b3e76d899efd4f51adf25d2fec9..4dded1f5d16caba29ad868744e226b1960725fef 100644
--- a/Projects/Strasse/macro/Control.cxx
+++ b/Projects/Strasse/macro/Control.cxx
@@ -25,7 +25,7 @@ void Control(){
   PhysicsTree->Draw("OuterPosZ:fDetected_Position_Z[3]",cond.c_str(),"col") ; 
   cOuter->cd(4);
   PhysicsTree->Draw("OuterPosY:OuterPosX:OuterPosZ", cond.c_str(),"");            
-*/
+
   TCanvas* cVertex = new TCanvas("ControlVertex","ControlVertex",1000,1000);
   cVertex->Divide(2,2);
   cond = "VertexX!=-1000";
@@ -44,14 +44,19 @@ void Control(){
   cdelta->Divide(2,2);
   cond = "deltaX!=-1000";
   cdelta->cd(1);
-  PhysicsTree->Draw("deltaX>>dx(500,-2,2)",cond.c_str(),"col") ; 
+  PhysicsTree->Draw("deltaX>>dx(500,-0.01,0.01)",cond.c_str(),"col") ; 
   cdelta->cd(2);
-  PhysicsTree->Draw("deltaY>>dy(500,-2,2)",cond.c_str(),"col") ; 
+  PhysicsTree->Draw("deltaY>>dy(500,-0.01,0.01)",cond.c_str(),"col") ; 
   cdelta->cd(3);
-  PhysicsTree->Draw("deltaZ>>dz(5000,-20,20)",cond.c_str(),"col") ; 
+  PhysicsTree->Draw("deltaZ>>dz(500,-0.01,0.01)",cond.c_str(),"col") ; 
   cdelta->cd(4);
-  PhysicsTree->Draw("Distance>>hd(500,0,80)", cond.c_str(),"");           
+  PhysicsTree->Draw("Distance>>hd(500,0,0.01)", cond.c_str(),"");           
+*/
+  TCanvas* ctheta= new TCanvas("ControlTheta","ControlTheta",1000,1000);
+  cond = "Theta12!=-1000";
+  PhysicsTree->Draw("Theta12>>ht",cond.c_str(),"col") ; 
  
+
 /*
   TCanvas* c2 = new TCanvas("Control 2", "Control2",500,500,2000,1000);  
   c2->Divide(2,1);
diff --git a/Projects/Strasse/macro/Shift.cxx b/Projects/Strasse/macro/Shift.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..057b95256ca840b97a3b9439b79870bfe6d8b661
--- /dev/null
+++ b/Projects/Strasse/macro/Shift.cxx
@@ -0,0 +1,14 @@
+void Shift(){
+  TFile* file_ok = TFile::Open("../../Outputs/Analysis/strasse_ok.root");
+  TFile* file_shifted = TFile::Open("../../Outputs/Analysis/strasse_shifted.root");
+  TTree* ok= (TTree*) file_ok->FindObjectAny("PhysicsTree");
+  TTree* shifted= (TTree*) file_shifted->FindObjectAny("PhysicsTree");
+  TCanvas* ctheta= new TCanvas("ControlTheta","ControlTheta",1000,1000);
+  string cond = "Theta12!=-1000";
+  ok->Draw("Theta12>>ht(5000)",cond.c_str(),"") ; 
+  shifted->Draw("Theta12>>hts(5000)",cond.c_str(),"same") ; 
+  TH1* hts= (TH1*) gDirectory->FindObjectAny("hts");
+  hts->SetFillColor(kOrange-3);
+  hts->SetLineColor(kOrange-3);
+
+}