Record changes.
[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.cxx
CommitLineData
44347160 1/**************************************************************************
572f41f9 2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
44347160 3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
572f41f9 16#include "AliITSRecoParam.h"
17
18/* $Id$ */
44347160 19
20///////////////////////////////////////////////////////////////////////////////
21// //
22// Class with ITS reconstruction parameters //
23// Origin: andrea.dainese@lnl.infn.it //
24// //
25///////////////////////////////////////////////////////////////////////////////
26
27
44347160 28
29ClassImp(AliITSRecoParam)
30
e50912db 31const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
32const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
33const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
34const Double_t AliITSRecoParam::fgkriw=80.0;
35const Double_t AliITSRecoParam::fgkdiw=0.0053;
36const Double_t AliITSRecoParam::fgkX0iw=30.0;
37const Double_t AliITSRecoParam::fgkrcd=61.0;
38const Double_t AliITSRecoParam::fgkdcd=0.0053;
39const Double_t AliITSRecoParam::fgkX0cd=30.0;
40const Double_t AliITSRecoParam::fgkyr=12.8;
41const Double_t AliITSRecoParam::fgkdr=0.03;
42const Double_t AliITSRecoParam::fgkzm=0.2;
43const Double_t AliITSRecoParam::fgkdm=0.40;
44const Double_t AliITSRecoParam::fgkrs=50.0;
45const Double_t AliITSRecoParam::fgkds=0.001;
46const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
47const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
48const Double_t AliITSRecoParam::fgkrPipe=3.;
49const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
50const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
51const Double_t AliITSRecoParam::fgkdPipe=0.0028;
52const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
53const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
54const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
55const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
56const Double_t AliITSRecoParam::fgkX0Air=21.82;
57const Double_t AliITSRecoParam::fgkX0Be=65.19;
58const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
59const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
60const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
61const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
44347160 62
63//_____________________________________________________________________________
572f41f9 64AliITSRecoParam::AliITSRecoParam() : TObject(),
e50912db 65fMaxSnp(1.),
572f41f9 66fNSigmaYLayerForRoadY(0),
67fNSigmaRoadY(0),
68fNSigmaZLayerForRoadZ(0),
69fNSigmaRoadZ(0),
70fNSigma2RoadZC(0),
71fNSigma2RoadYC(0),
72fNSigma2RoadZNonC(0),
73fNSigma2RoadYNonC(0),
74fMaxNormChi2NonCForHypothesis(0),
75fMaxChi2(0),
76fMaxRoad(0),
77fMaxChi2In(0),
78fChi2PerCluster(0),
79fXV(0),
80fYV(0),
81fZV(0),
82fSigmaXV(0),
83fSigmaYV(0),
84fSigmaZV(0),
85fVertexCut(0),
86fMaxDZforPrimTrk(0),
87fMaxDZToUseConstraint(0),
88fMaxDforV0dghtrForProlongation(0),
89fMaxDForProlongation(0),
90fMaxDZForProlongation(0),
91fMinPtForProlongation(0),
92fAddVirtualClustersInDeadZone(kTRUE),
93fZWindowDeadZone(0),
94fSigmaXDeadZoneHit2(0),
95fSigmaZDeadZoneHit2(0),
96fXPassDeadZoneHits(0),
e50912db 97fUseTGeoInTracker(0),
572f41f9 98fAllowSharedClusters(kTRUE),
e50912db 99fClusterErrorsParam(1),
572f41f9 100fFindV0s(kTRUE)
44347160 101{
102 //
103 // constructor
104 //
105 SetLayersParameters();
e50912db 106 SetUseTGeoInTracker(0);
afd25725 107 SetAllowSharedClusters(kTRUE);
108 SetFindV0s(kTRUE);
109 SetAddVirtualClustersInDeadZone(kTRUE);
572f41f9 110 SetUseAmplitudeInfo(kTRUE);
e50912db 111 SetClusterErrorsParam(1);
44347160 112}
113//_____________________________________________________________________________
114AliITSRecoParam::~AliITSRecoParam()
115{
116 //
117 // destructor
118 //
119}
120//_____________________________________________________________________________
121AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
122{
123 //
124 // make default reconstruction parameters for hig flux env.
125 //
126 AliITSRecoParam *param = new AliITSRecoParam();
127
128 param->fMaxSnp = 0.95;
129
130 param->fNSigmaYLayerForRoadY = 4.;
131 param->fNSigmaRoadY = 7.5;
132 param->fNSigmaZLayerForRoadZ = 4.;
133 param->fNSigmaRoadZ = 7.5;
134
135 param->fNSigma2RoadZC = 60.; //7.75^2
136 param->fNSigma2RoadYC = 60.; //7.75^2
137 param->fNSigma2RoadZNonC = 50.; //7.07^2
138 param->fNSigma2RoadYNonC = 50.; //7.07^2
139
140 param->fMaxChi2PerCluster[0] = 11.; //7
141 param->fMaxChi2PerCluster[1] = 12.; //5
142 param->fMaxChi2PerCluster[2] = 12.; //8
143 param->fMaxChi2PerCluster[3] = 5.; //8
144 param->fMaxChi2PerCluster[4] = 12.; //6.5
145
146 param->fMaxNormChi2NonC[0] = 7.;
147 param->fMaxNormChi2NonC[1] = 8.;
148 param->fMaxNormChi2NonC[2] = 8.;
149 param->fMaxNormChi2NonC[3] = 11.;
150 param->fMaxNormChi2NonC[4] = 14.;
151 param->fMaxNormChi2NonC[5] = 25.;
152
153 param->fMaxNormChi2C[0] = 11.;
154 param->fMaxNormChi2C[1] = 13.;
155 param->fMaxNormChi2C[2] = 15.;
156 param->fMaxNormChi2C[3] = 18.;
157 param->fMaxNormChi2C[4] = 30.;
158 param->fMaxNormChi2C[5] = 35.;
afd25725 159
160 param->fMaxNormChi2NonCForHypothesis = 7.;
44347160 161
162 param->fMaxChi2 = 35.;
163
164 param->fMaxChi2s[0] = 25.; //40
165 param->fMaxChi2s[1] = 25.; //40
166 param->fMaxChi2s[2] = 25.; //40
167 param->fMaxChi2s[3] = 25.; //40
168 param->fMaxChi2s[4] = 40.; //40
169 param->fMaxChi2s[5] = 50.; //40
170
171 param->fMaxRoad = 6.;
172
173 // not used
174 param->fMaxChi2In = 16.;
175
176 param->fMaxChi2sR[0] = 10.;
177 param->fMaxChi2sR[1] = 10.;
178 param->fMaxChi2sR[2] = 10.;
179 param->fMaxChi2sR[3] = 10.;
180 param->fMaxChi2sR[4] = 30.;
181 param->fMaxChi2sR[5] = 40.;
182
183 param->fChi2PerCluster = 9.;
184 // not used
185
186 param->fXV = 0.;
187 param->fYV = 0.;
188 param->fZV = 0.;
189 param->fSigmaXV = 0.0050;
190 param->fSigmaYV = 0.0050;
191 param->fSigmaZV = 0.0100;
192
193 param->fVertexCut = 25.;
afd25725 194
195 param->fMaxDZforPrimTrk = 0.4;
196 param->fMaxDZToUseConstraint = 3.;
197
198 param->fMaxDforV0dghtrForProlongation = 30.;
199 param->fMaxDForProlongation = 10.;
200 param->fMaxDZForProlongation = 20.;
201 param->fMinPtForProlongation = 0.120;
202
203 param->fZWindowDeadZone = 2.0;
204 param->fSigmaXDeadZoneHit2 = 0.004/12.;
205 param->fSigmaZDeadZoneHit2 = 0.001/12.;
206 param->fXPassDeadZoneHits = 0.018;
207
44347160 208
209 return param;
210}
211//_____________________________________________________________________________
212AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
213{
214 //
215 // make default reconstruction parameters for low flux env.
216 //
217 return GetHighFluxParam();
218}
219//_____________________________________________________________________________
220AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
221{
222 //
223 // make default reconstruction parameters for cosmics
224 //
225 return GetHighFluxParam();
226}
227//_____________________________________________________________________________
228void AliITSRecoParam::SetLayersParameters()
229{
230 //
231 // number of layers and layers spatial resolutions
232 //
233
234 // spatial resolutions of the detectors
235 // y: 12 12 38 38 20 20 micron
236 fSigmaY2[0]=1.44e-6;
237 fSigmaY2[1]=1.44e-6;
238 fSigmaY2[2]=1.444e-5;
239 fSigmaY2[3]=1.444e-5;
240 fSigmaY2[4]=4.0e-6;
241 fSigmaY2[5]=4.0e-6;
242 // z: 120 120 28 28 830 830 micron
243 fSigmaZ2[0]=1.44e-4;
244 fSigmaZ2[1]=1.44e-4;
245 fSigmaZ2[2]=7.84e-6;
246 fSigmaZ2[3]=7.84e-6;
247 fSigmaZ2[4]=6.889e-3;
248 fSigmaZ2[5]=6.889e-3;
249
250 return;
251}
252