]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
Possibility of angle correction in CorrectForMeanMaterial (A.Dainese)
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 5 Oct 2007 15:18:08 +0000 (15:18 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 5 Oct 2007 15:18:08 +0000 (15:18 +0000)
STEER/AliExternalTrackParam.cxx
STEER/AliExternalTrackParam.h

index 4b63368379bbf26e859aa4735bf79faec945624b..84f001e2e525debeab520184abb5bcd605cbf037 100644 (file)
@@ -201,8 +201,8 @@ Double_t AliExternalTrackParam::GetLinearD(Double_t xv,Double_t yv) const {
 }
 
 Bool_t AliExternalTrackParam::CorrectForMeanMaterial
-(Double_t xOverX0,  Double_t xTimesRho, Double_t mass, 
-Double_t (*Bethe)(Double_t)) {
+(Double_t xOverX0,  Double_t xTimesRho, Double_t mass, Bool_t anglecorr, 
+ Double_t (*Bethe)(Double_t)) {
   //------------------------------------------------------------------
   // This function corrects the track parameters for the crossed material.
   // "xOverX0"   - X/X0, the thickness in units of the radiation length.
@@ -218,6 +218,13 @@ Double_t (*Bethe)(Double_t)) {
   Double_t &fC43=fC[13];
   Double_t &fC44=fC[14];
 
+  //Apply angle correction, if requested
+  if(anglecorr) {
+    Double_t angle=TMath::Sqrt((1.+ fP3*fP3)/(1.- fP2*fP2));
+    xOverX0 *=angle;
+    xTimesRho *=angle;
+  } 
+
   Double_t p=GetP();
   Double_t p2=p*p;
   Double_t beta2=p2/(p2 + mass*mass);
index 7916b85e2037cf33dc540a369c28429f2ddd206d..e1b5af9f9f925dc635b75e7a4a1b22e65f9d26fb 100644 (file)
@@ -106,7 +106,8 @@ class AliExternalTrackParam: public AliVParticle {
   Double_t GetD(Double_t xv, Double_t yv, Double_t b) const; 
   Double_t GetLinearD(Double_t xv, Double_t yv) const; 
   Bool_t CorrectForMeanMaterial(Double_t xOverX0, Double_t xTimesRho, 
-        Double_t mass, Double_t (*f)(Double_t)=ApproximateBetheBloch);
+        Double_t mass,  Bool_t anglecorr=kFALSE,
+       Double_t (*f)(Double_t)=ApproximateBetheBloch);
   Double_t GetPredictedChi2(Double_t p[2],Double_t cov[3]) const;
 
   Double_t