Modifications used for addendum to Dimuon TDR (JP Cussonneau):
[u/mrichter/AliRoot.git] / MUON / AliMUONSegment.cxx
index 01671570e80cc2e02eb63c71086850806c7e275b..fcf369a5750deab291b85f31e2c0e9f46b9fa261 100644 (file)
 
 /*
 $Log$
 
 /*
 $Log$
+Revision 1.4  2000/06/30 10:15:48  gosset
+Changes to EventReconstructor...:
+precision fit with multiple Coulomb scattering;
+extrapolation to vertex with Branson correction in absorber (JPC)
+
 Revision 1.3  2000/06/25 13:06:39  hristov
 Inline functions moved from *.cxx to *.h files instead of forward declarations
 
 Revision 1.3  2000/06/25 13:06:39  hristov
 Inline functions moved from *.cxx to *.h files instead of forward declarations
 
@@ -253,18 +258,27 @@ void AliMUONSegment::UpdateFromStationTrackParam(AliMUONTrackParam *TrackParam,
   // The "road" is parametrized from the old reco_muon.F
   // with 8 cm between stations.
   AliMUONTrackParam *param0;
   // The "road" is parametrized from the old reco_muon.F
   // with 8 cm between stations.
   AliMUONTrackParam *param0;
-  Double_t cReso2, sReso2;
+//   Double_t cReso2, sReso2;
   // parameters to define the widths of the searching roads in station 0,1,2
   // width = p0 + p1/ (momentum)^2
   //                  station number:        0         1          2
   // parameters to define the widths of the searching roads in station 0,1,2
   // width = p0 + p1/ (momentum)^2
   //                  station number:        0         1          2
-  static Double_t p0BendingCoor[3] =     { 6.43e-2, 1.64e-2,   0.034 };   
-  static Double_t p1BendingCoor[3] =     {    986.,    821.,    446. };  
-  static Double_t p0BendingSlope[3] =    { 3.54e-6, 3.63e-6,  3.6e-6 };  
-  static Double_t p1BendingSlope[3] =    { 4.49e-3,  4.8e-3,   0.011 };  
-  static Double_t p0NonBendingCoor[3] =  { 4.66e-2, 4.83e-2,   0.049 };   
-  static Double_t p1NonBendingCoor[3] =  {   1444.,    866.,    354. };  
-  static Double_t p0NonBendingSlope[3] = { 6.14e-4, 6.49e-4, 6.85e-4 };  
-  static Double_t p1NonBendingSlope[3] = {      0.,      0.,      0. };  
+//   static Double_t p0BendingCoor[3] =     { 6.43e-2, 1.64e-2,   0.034 };   
+//   static Double_t p1BendingCoor[3] =     {    986.,    821.,    446. };  
+//   static Double_t p0BendingSlope[3] =    { 3.54e-6, 3.63e-6,  3.6e-6 };  
+//   static Double_t p1BendingSlope[3] =    { 4.49e-3,  4.8e-3,   0.011 };  
+//   static Double_t p0NonBendingCoor[3] =  { 4.66e-2, 4.83e-2,   0.049 };   
+//   static Double_t p1NonBendingCoor[3] =  {   1444.,    866.,    354. };  
+//   static Double_t p0NonBendingSlope[3] = { 6.14e-4, 6.49e-4, 6.85e-4 };  
+//   static Double_t p1NonBendingSlope[3] = {      0.,      0.,      0. };
+  
+  static Double_t p0BendingCoor[3] =     { 6.43e-2, 6.43e-2,   6.43e-2  };   
+  static Double_t p1BendingCoor[3] =     {    986.,    986.,       986. };  
+  static Double_t p0BendingSlope[3] =    {   3.6e-6,   3.6e-6,     3.6e-6  };  
+  static Double_t p1BendingSlope[3] =    {  1.1e-2,  1.1e-2,    1.1e-2  };  
+  static Double_t p0NonBendingCoor[3] =  {   0.049,   0.049,     0.049  };   
+  static Double_t p1NonBendingCoor[3] =  {   1444.,   1444.,      1444. };  
+  static Double_t p0NonBendingSlope[3] = {   6.8e-4,   6.8e-4,     6.8e-4  };  
+  static Double_t p1NonBendingSlope[3] = {      0.,      0.,         0. };  
   param0 = &(TrackParam[0]);
 
 // OLD version
   param0 = &(TrackParam[0]);
 
 // OLD version
@@ -296,15 +310,15 @@ void AliMUONSegment::UpdateFromStationTrackParam(AliMUONTrackParam *TrackParam,
   // because they are added in the functions "NormalizedChi2WithSegment"
   // and "NormalizedChi2WithHitForRec"
   // Bending plane
   // because they are added in the functions "NormalizedChi2WithSegment"
   // and "NormalizedChi2WithHitForRec"
   // Bending plane
-  cReso2 = fBendingCoorReso2;
-  sReso2 = (2. * cReso2 )/ (Dz3*Dz3) ;
-  fBendingCoorReso2 = p0BendingCoor[Station] + p1BendingCoor[Station]*InverseMomentum*InverseMomentum - cReso2;
-  fBendingSlopeReso2 = p0BendingSlope[Station] + p1BendingSlope[Station]*InverseMomentum*InverseMomentum - sReso2;
+//   cReso2 = fBendingCoorReso2;
+//   sReso2 = (2. * cReso2 )/ (Dz3*Dz3) ;
+  fBendingCoorReso2 = p0BendingCoor[Station] + p1BendingCoor[Station]*InverseMomentum*InverseMomentum ;  // - cReso2
+  fBendingSlopeReso2 = p0BendingSlope[Station] + p1BendingSlope[Station]*InverseMomentum*InverseMomentum; //  - sReso2;
   // Non bending plane
   // Non bending plane
-  cReso2 = fNonBendingCoorReso2;
-  sReso2 =  (2. * cReso2 )/ (Dz3*Dz3) ;
-  fNonBendingCoorReso2 = p0NonBendingCoor[Station] + p1NonBendingCoor[Station]*InverseMomentum*InverseMomentum - cReso2;
-  fNonBendingSlopeReso2 = p0NonBendingSlope[Station] + p1NonBendingSlope[Station]*InverseMomentum*InverseMomentum - sReso2;
+//   cReso2 = fNonBendingCoorReso2;
+//   sReso2 =  (2. * cReso2 )/ (Dz3*Dz3) ;
+  fNonBendingCoorReso2 = p0NonBendingCoor[Station] + p1NonBendingCoor[Station]*InverseMomentum*InverseMomentum; // - cReso2;
+  fNonBendingSlopeReso2 = p0NonBendingSlope[Station] + p1NonBendingSlope[Station]*InverseMomentum*InverseMomentum; //  - sReso2;
   return;
 }
 
   return;
 }