static void ExtrapToZCov(AliMUONTrackParam* trackParam, Double_t zEnd);
static void ExtrapToStation(AliMUONTrackParam *trackParamIn, Int_t station, AliMUONTrackParam *trackParamOut);
static void ExtrapToVertexUncorrected(AliMUONTrackParam* trackParam, Double_t zVtx);
- static void ExtrapToVertex(AliMUONTrackParam *trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx);
+ static void ExtrapToVertex(AliMUONTrackParam *trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx,
+ Bool_t CorrectForMCS = kTRUE, Bool_t CorrectForEnergyLoss = kTRUE);
- static void AddMCSEffectInTrackParamCov(AliMUONTrackParam *param, Double_t dZ, Double_t x0);
+ static Double_t TotalMomentumEnergyLoss(AliMUONTrackParam* trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx);
+
+ static void AddMCSEffect(AliMUONTrackParam *param, Double_t dZ, Double_t x0);
static void ExtrapOneStepRungekutta(Double_t charge, Double_t step, Double_t* vect, Double_t* vout);
static const Double_t fgkRungeKuttaMaxResidue;//!< Maximal distance (in Z) to destination to stop the track extrapolation (used in Runge-Kutta)
// Functions
+
+ /// Not implemented
AliMUONTrackExtrap(const AliMUONTrackExtrap& trackExtrap);
+ /// Not implemented
AliMUONTrackExtrap& operator=(const AliMUONTrackExtrap& trackExtrap);
static void ExtrapToZHelix(AliMUONTrackParam *trackParam, Double_t Z);
static void ConvertTrackParamForExtrap(AliMUONTrackParam* trackParam, Double_t forwardBackward, Double_t *v3);
static void RecoverTrackParam(Double_t *v3, Double_t Charge, AliMUONTrackParam* trackParam);
- static void BransonCorrection(AliMUONTrackParam *trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx);
- static Double_t TotalMomentumEnergyLoss(Double_t thetaLimit, Double_t pTotal, Double_t theta);
- static void FieldCorrection(AliMUONTrackParam *trackParam, Double_t Z);
+ static void AddMCSEffectInAbsorber(AliMUONTrackParam* trackParam, Double_t pathLength, Double_t f0, Double_t f1, Double_t f2);
+ static void GetAbsorberCorrectionParam(Double_t trackXYZIn[3], Double_t trackXYZOut[3], Double_t &pathLength,
+ Double_t &f0, Double_t &f1, Double_t &f2, Double_t &meanRho);
+
+ static Double_t TotalMomentumEnergyLoss(Double_t pTotal, Double_t pathLength, Double_t rho);
+ static Double_t ApproximateBetheBloch(Double_t beta2);
static void ExtrapOneStepHelix(Double_t charge, Double_t step, Double_t *vect, Double_t *vout);
static void ExtrapOneStepHelix3(Double_t field, Double_t step, Double_t *vect, Double_t *vout);