X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONTrackExtrap.h;h=66ee9ed2d116be1187f5d1b5682316fd6e56c6f9;hb=ade4e6f9b75c6e778780ad9bb9dfd80200612999;hp=7faf95dfe0b1efe9f874b6a2f81dc360b8834ddc;hpb=fac70e254db1e9dc11f104c2f914eb4073bf4a3c;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONTrackExtrap.h b/MUON/AliMUONTrackExtrap.h index 7faf95dfe0b..66ee9ed2d11 100644 --- a/MUON/AliMUONTrackExtrap.h +++ b/MUON/AliMUONTrackExtrap.h @@ -15,6 +15,7 @@ ////////////////////////////////////////////////////////////// #include +#include class AliMagF; class AliMUONTrackParam; @@ -27,54 +28,100 @@ class AliMUONTrackExtrap : public TObject /// Destructor virtual ~AliMUONTrackExtrap(){}; - /// set field map - static void SetField(const AliMagF* magField) {fgkField = magField;} + static void SetField(); + + /// return kTRUE if the field is switched ON + static Bool_t IsFieldON() {return fgFieldON;} static Double_t GetImpactParamFromBendingMomentum(Double_t bendingMomentum); static Double_t GetBendingMomentumFromImpactParam(Double_t impactParam); - static void ExtrapToZ(AliMUONTrackParam *trackParam, Double_t Z); - static void ExtrapToZCov(AliMUONTrackParam* trackParam, Double_t zEnd); - static void ExtrapToStation(AliMUONTrackParam *trackParamIn, Int_t station, AliMUONTrackParam *trackParamOut); + // Linearly extrapolate track parameters + static void LinearExtrapToZ(AliMUONTrackParam* trackParam, Double_t zEnd); + + // Linearly extrapolate track parameters and covariances + static void LinearExtrapToZCov(AliMUONTrackParam* trackParam, Double_t zEnd, Bool_t updatePropagator = kFALSE); + + // Extrapolate track parameters in magnetic field + static Bool_t ExtrapToZ(AliMUONTrackParam *trackParam, Double_t zEnd); + + // Extrapolate track parameters and covariances in magnetic field + static Bool_t ExtrapToZCov(AliMUONTrackParam* trackParam, Double_t zEnd, Bool_t updatePropagator = kFALSE); + + // Extrapolate track parameters to vertex, corrected for multiple scattering and energy loss effects + // Add branson correction resolution and energy loss fluctuation to parameter covariances + static void ExtrapToVertex(AliMUONTrackParam* trackParam, + Double_t xVtx, Double_t yVtx, Double_t zVtx, + Double_t errXVtx, Double_t errYVtx); + + // Extrapolate track parameters to vertex, corrected for multiple scattering effects only + // Add branson correction resolution to parameter covariances + static void ExtrapToVertexWithoutELoss(AliMUONTrackParam* trackParam, + Double_t xVtx, Double_t yVtx, Double_t zVtx, + Double_t errXVtx, Double_t errYVtx); + + // Extrapolate track parameters to vertex, corrected for energy loss effects only + // Add dispersion due to multiple scattering and energy loss fluctuation to parameter covariances + static void ExtrapToVertexWithoutBranson(AliMUONTrackParam* trackParam, Double_t zVtx); + + // Extrapolate track parameters to vertex without multiple scattering and energy loss corrections + // Add dispersion due to multiple scattering to parameter covariances static void ExtrapToVertexUncorrected(AliMUONTrackParam* trackParam, 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 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 Double_t GetMCSAngle2(const AliMUONTrackParam& param, Double_t dZ, Double_t x0); + 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 Bool_t ExtrapOneStepRungekutta(Double_t charge, Double_t step, Double_t* vect, Double_t* vout); private: - static const AliMagF* fgkField; //!< field map + static const Double_t fgkSimpleBPosition; //!< position of the dipole + static const Double_t fgkSimpleBLength; //!< length of the dipole + static Double_t fgSimpleBValue; //!< magnetic field value at the centre + static Bool_t fgFieldON; //!< kTRUE if the field is switched ON static const Bool_t fgkUseHelix; //!< Tell whether to use Helix or not (default is Runge-Kutta) static const Int_t fgkMaxStepNumber; //!< Maximum number of steps for track extrapolation static const Double_t fgkHelixStepLength; //!< Step lenght for track extrapolation (used in Helix) 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 ExtrapToZRungekutta(AliMUONTrackParam *trackParam, Double_t Z); + static Bool_t ExtrapToZHelix(AliMUONTrackParam *trackParam, Double_t Z); + static Bool_t ExtrapToZRungekutta(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 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 void ExtrapToVertex(AliMUONTrackParam* trackParam, + Double_t xVtx, Double_t yVtx, Double_t zVtx, + Double_t errXVtx, Double_t errYVtx, + Bool_t correctForMCS, Bool_t correctForEnergyLoss); - static Double_t TotalMomentumEnergyLoss(Double_t pTotal, Double_t pathLength, Double_t rho); - static Double_t ApproximateBetheBloch(Double_t beta2); + static void AddMCSEffectInAbsorber(AliMUONTrackParam* trackParam, Double_t signedPathLength, Double_t f0, Double_t f1, Double_t f2); + static void CorrectMCSEffectInAbsorber(AliMUONTrackParam* param, + Double_t xVtx, Double_t yVtx, Double_t zVtx, + Double_t errXVtx, Double_t errYVtx, + Double_t absZBeg, Double_t pathLength, Double_t f0, Double_t f1, Double_t f2); + static void CorrectELossEffectInAbsorber(AliMUONTrackParam* param, Double_t eLoss, Double_t sigmaELoss2); + static Bool_t GetAbsorberCorrectionParam(Double_t trackXYZIn[3], Double_t trackXYZOut[3], Double_t pTotal, + Double_t &pathLength, Double_t &f0, Double_t &f1, Double_t &f2, + Double_t &meanRho, Double_t &totalELoss, Double_t &sigmaELoss2); + + static Double_t BetheBloch(Double_t pTotal, Double_t pathLength, Double_t rho, Double_t atomicZ, Double_t atomicZoverA); + static Double_t EnergyLossFluctuation2(Double_t pTotal, Double_t pathLength, Double_t rho, Double_t atomicZoverA); + + static void Cov2CovP(const TMatrixD ¶m, TMatrixD &cov); + static void CovP2Cov(const TMatrixD ¶m, TMatrixD &cov); 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); - - static void GetField(Double_t *Position, Double_t *Field); ClassDef(AliMUONTrackExtrap, 0) // Tools for track extrapolation in ALICE dimuon spectrometer };