1 #ifndef ALIITSURECOPARAM_H
2 #define ALIITSURECOPARAM_H
3 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 /* $Id: AliITSURecoParam.h 57215 2012-06-17 14:47:08Z masera $ */
8 ///////////////////////////////////////////////////////////////////////////////
10 // Class with ITS reconstruction parameters //
12 ///////////////////////////////////////////////////////////////////////////////
14 class AliITSUTrackCond;
16 #include <TObjArray.h>
17 #include "AliDetectorRecoParam.h"
19 class AliITSURecoParam : public AliDetectorRecoParam
22 enum {kClus2Tracks,kPropBack,kRefitInw,kNTrackingPhases}; // tracking phases
25 AliITSURecoParam(Int_t nLr);
26 virtual ~AliITSURecoParam();
28 static AliITSURecoParam *GetLowFluxParam();// make reco parameters for low flux env.
29 static AliITSURecoParam *GetHighFluxParam();// make reco parameters for high flux env.
30 static AliITSURecoParam *GetCosmicTestParam();// special setting for cosmic
32 Double_t GetMaxDforV0dghtrForProlongation() const {return fMaxDforV0dghtrForProlongation;}
33 Double_t GetMaxDForProlongation() const {return fMaxDForProlongation;}
34 Double_t GetMaxDZForProlongation() const {return fMaxDZForProlongation;}
35 Double_t GetMinPtForProlongation() const {return fMinPtForProlongation;}
36 Double_t GetTanLorentzAngle(Int_t lr) const;
37 Double_t GetSigmaY2(Int_t lr) const;
38 Double_t GetSigmaZ2(Int_t lr) const;
39 Bool_t GetAllowDiagonalClusterization(Int_t lr) const;
41 Double_t GetTPCITSWallRMin() const {return fTPCITSWallRMin;}
42 Double_t GetTPCITSWallRMax() const {return fTPCITSWallRMax;}
43 Double_t GetTPCITSWallZSpanH() const {return fTPCITSWallZSpanH;}
44 Double_t GetTPCITSWallMaxStep() const {return fTPCITSWallMaxStep;}
45 TObjArray* GetTrackingConditions() const {return (TObjArray*)&fTrackingConditions;}
46 Int_t GetNTrackingConditions() const {return fTrackingConditions.GetEntriesFast();}
47 AliITSUTrackCond* GetTrackingCondition(Int_t i) const {return (AliITSUTrackCond*)fTrackingConditions[i];}
48 Bool_t GetUseMatLUT(Int_t step) const {return (step<0||step>=kNTrackingPhases) ? kFALSE:fUseMatLUT[step];}
49 Int_t GetTracker() const {return fTracker;}
50 Bool_t GetSAonly() const {return fSAonly;}
52 void SetNLayers(Int_t n);
53 void SetTanLorentzAngle(Int_t lr, Double_t v);
54 void SetSigmaY2(Int_t lr, Double_t v);
55 void SetSigmaZ2(Int_t lr, Double_t v);
56 void SetAllowDiagonalClusterization(Int_t lr, Bool_t v);
58 void SetMaxDforV0dghtrForProlongation(Double_t v) {fMaxDforV0dghtrForProlongation = v;}
59 void SetMaxDForProlongation(Double_t v) {fMaxDForProlongation = v;}
60 void SetMaxDZForProlongation(Double_t v) {fMaxDZForProlongation = v;}
61 void SetMinPtForProlongation(Double_t v) {fMinPtForProlongation = v;}
63 void SetTPCITSWallRMin(double v) {fTPCITSWallRMin = v;}
64 void SetTPCITSWallRMax(double v) {fTPCITSWallRMax = v;}
65 void SetTPCITSWallZSpanH(double v) {fTPCITSWallZSpanH = v;}
66 void SetTPCITSWallMaxStep(double v) {fTPCITSWallMaxStep = v;}
68 void SetUseMatLUT(Int_t step, Bool_t v) {if (step>=0&&step<kNTrackingPhases) fUseMatLUT[step]=v;}
69 void AddTrackingCondition(AliITSUTrackCond* cond);
70 virtual void Print(Option_t *opt="") const;
71 void SetTracker(Int_t opt) {fTracker=opt;}
72 void SetSAonly(Bool_t sa=kTRUE) {fSAonly=sa;}
75 Int_t fNLayers; // number of layers
77 Double_t fMaxDforV0dghtrForProlongation; // max. rphi imp. par. cut for V0 daughter
78 Double_t fMaxDForProlongation; // max. rphi imp. par. cut
79 Double_t fMaxDZForProlongation; // max. 3D imp. par. cut
80 Double_t fMinPtForProlongation; // min. pt cut
82 Double_t fTPCITSWallRMin; // minR
83 Double_t fTPCITSWallRMax; // maxR
84 Double_t fTPCITSWallZSpanH; // half Z span
85 Double_t fTPCITSWallMaxStep; // max tracking step
87 Bool_t fUseMatLUT[kNTrackingPhases]; // allow usage of material lookup table at given tracking step
88 Bool_t* fAllowDiagonalClusterization; //[fNLayers] allow clusters of pixels with common corners only
89 Double_t* fTanLorentzAngle; //[fNLayers] Optional Lorentz angle for each layer
90 Double_t* fSigmaY2; //[fNLayers] addition to road width^2
91 Double_t* fSigmaZ2; //[fNLayers] addition to road width^2
93 TObjArray fTrackingConditions; // array of tracking conditions for different iterations
95 static const Double_t fgkMaxDforV0dghtrForProlongation; // default
96 static const Double_t fgkMaxDForProlongation; // default
97 static const Double_t fgkMaxDZForProlongation; // default
98 static const Double_t fgkMinPtForProlongation; // default
99 static const Double_t fgkTanLorentzAngle; // default
100 static const Double_t fgkSigmaRoadY; // default
101 static const Double_t fgkSigmaRoadZ; // default
103 // hardwired params for TPC-ITS border layer
104 static const Double_t fgkTPCITSWallRMin; // fiducial R min
105 static const Double_t fgkTPCITSWallRMax; // fiducial R max
106 static const Double_t fgkTPCITSWallZSpanH; // half Z span
107 static const Double_t fgkTPCITSWallMaxStep; // max tracking step
109 static const Bool_t fgkUseMatLUT[kNTrackingPhases]; // def. mat. LUT usage
110 static const Bool_t fgkAllowDiagonalClusterization; // clusters of pixels with common corners
112 Int_t fTracker; // 0 - global, 1 - cooked, 2 - CA
113 Bool_t fSAonly; // kTRUE, if ITSU standalone tracking only
115 AliITSURecoParam(const AliITSURecoParam & param);
116 AliITSURecoParam & operator=(const AliITSURecoParam ¶m);
118 ClassDef(AliITSURecoParam,6) // ITS reco parameters
121 //_____________________________________________________________________________
122 inline Double_t AliITSURecoParam::GetTanLorentzAngle(Int_t lr) const
124 // get tg of Lorentz Angle for the layer
125 return (lr<fNLayers)&&fTanLorentzAngle ? fTanLorentzAngle[lr]:0;
128 //_____________________________________________________________________________
129 inline Double_t AliITSURecoParam::GetSigmaY2(Int_t lr) const
131 // get tg of Lorentz Angle for the layer
132 return (lr<fNLayers)&&fSigmaY2 ? fSigmaY2[lr]:fgkSigmaRoadY*fgkSigmaRoadY; //0;
135 //_____________________________________________________________________________
136 inline Double_t AliITSURecoParam::GetSigmaZ2(Int_t lr) const
138 // get tg of Lorentz Angle for the layer
139 return (lr<fNLayers)&&fSigmaZ2 ? fSigmaZ2[lr]:fgkSigmaRoadZ*fgkSigmaRoadZ;//0;
142 //_____________________________________________________________________________
143 inline Bool_t AliITSURecoParam::GetAllowDiagonalClusterization(Int_t lr) const
145 // are diagonal clusters permitted
146 return (lr<fNLayers)&&fAllowDiagonalClusterization ? fAllowDiagonalClusterization[lr]:fgkAllowDiagonalClusterization;