1 #ifndef ALIMUONTRACKEXTRAP_H
2 #define ALIMUONTRACKEXTRAP_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 // Revision of includes 07/05/2004
10 /// \class AliMUONTrackExtrap
11 /// \brief Track parameters in ALICE dimuon spectrometer
13 //////////////////////////////////////////////////////////////
14 /// Tools for track extrapolation in ALICE dimuon spectrometer
15 //////////////////////////////////////////////////////////////
20 class AliMUONTrackParam;
22 class AliMUONTrackExtrap : public TObject
26 AliMUONTrackExtrap() : TObject(){};
28 virtual ~AliMUONTrackExtrap(){};
31 static void SetField(const AliMagF* magField) {fgkField = magField;}
33 static Double_t GetImpactParamFromBendingMomentum(Double_t bendingMomentum);
34 static Double_t GetBendingMomentumFromImpactParam(Double_t impactParam);
36 static void LinearExtrapToZ(AliMUONTrackParam* trackParam, Double_t zEnd);
37 static void ExtrapToZ(AliMUONTrackParam *trackParam, Double_t zEnd);
38 static void ExtrapToZCov(AliMUONTrackParam* trackParam, Double_t zEnd, Bool_t updatePropagator = kFALSE);
39 static void ExtrapToStation(AliMUONTrackParam *trackParamIn, Int_t station, AliMUONTrackParam *trackParamOut);
40 static void ExtrapToVertexUncorrected(AliMUONTrackParam* trackParam, Double_t zVtx);
41 static void ExtrapToVertex(AliMUONTrackParam *trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx,
42 Bool_t CorrectForMCS = kTRUE, Bool_t CorrectForEnergyLoss = kTRUE);
44 static Double_t TotalMomentumEnergyLoss(AliMUONTrackParam* trackParam, Double_t xVtx, Double_t yVtx, Double_t zVtx);
46 static Double_t GetMCSAngle2(const AliMUONTrackParam& param, Double_t dZ, Double_t x0);
47 static void AddMCSEffect(AliMUONTrackParam *param, Double_t dZ, Double_t x0);
49 static void ExtrapOneStepRungekutta(Double_t charge, Double_t step, Double_t* vect, Double_t* vout);
53 static const AliMagF* fgkField; //!< field map
54 static const Bool_t fgkUseHelix; //!< Tell whether to use Helix or not (default is Runge-Kutta)
55 static const Int_t fgkMaxStepNumber; //!< Maximum number of steps for track extrapolation
56 static const Double_t fgkHelixStepLength; //!< Step lenght for track extrapolation (used in Helix)
57 static const Double_t fgkRungeKuttaMaxResidue;//!< Maximal distance (in Z) to destination to stop the track extrapolation (used in Runge-Kutta)
62 AliMUONTrackExtrap(const AliMUONTrackExtrap& trackExtrap);
64 AliMUONTrackExtrap& operator=(const AliMUONTrackExtrap& trackExtrap);
66 static void ExtrapToZHelix(AliMUONTrackParam *trackParam, Double_t Z);
67 static void ExtrapToZRungekutta(AliMUONTrackParam *trackParam, Double_t Z);
69 static void ConvertTrackParamForExtrap(AliMUONTrackParam* trackParam, Double_t forwardBackward, Double_t *v3);
70 static void RecoverTrackParam(Double_t *v3, Double_t Charge, AliMUONTrackParam* trackParam);
72 static void AddMCSEffectInAbsorber(AliMUONTrackParam* trackParam, Double_t pathLength, Double_t f0, Double_t f1, Double_t f2);
73 static void GetAbsorberCorrectionParam(Double_t trackXYZIn[3], Double_t trackXYZOut[3], Double_t &pathLength,
74 Double_t &f0, Double_t &f1, Double_t &f2, Double_t &meanRho);
76 static Double_t TotalMomentumEnergyLoss(Double_t pTotal, Double_t pathLength, Double_t rho);
77 static Double_t ApproximateBetheBloch(Double_t beta2);
79 static void ExtrapOneStepHelix(Double_t charge, Double_t step, Double_t *vect, Double_t *vout);
80 static void ExtrapOneStepHelix3(Double_t field, Double_t step, Double_t *vect, Double_t *vout);
82 static void GetField(Double_t *Position, Double_t *Field);
84 ClassDef(AliMUONTrackExtrap, 0) // Tools for track extrapolation in ALICE dimuon spectrometer