Possibility of selecting larger windows in SA tracker, in order to use the tracker...
[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.cxx
1 /**************************************************************************
2  * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
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
16 #include "AliITSRecoParam.h"
17
18 /* $Id$ */
19
20 ///////////////////////////////////////////////////////////////////////////////
21 //                                                                           //
22 // Class with ITS reconstruction parameters                                  //
23 // Origin: andrea.dainese@lnl.infn.it                                        //
24 //                                                                           //
25 ///////////////////////////////////////////////////////////////////////////////
26
27
28
29 ClassImp(AliITSRecoParam)
30
31 const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
32 const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
33 const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
34 const Double_t AliITSRecoParam::fgkriw=80.0;
35 const Double_t AliITSRecoParam::fgkdiw=0.0053;
36 const Double_t AliITSRecoParam::fgkX0iw=30.0;
37 const Double_t AliITSRecoParam::fgkrcd=61.0;
38 const Double_t AliITSRecoParam::fgkdcd=0.0053;
39 const Double_t AliITSRecoParam::fgkX0cd=30.0;
40 const Double_t AliITSRecoParam::fgkyr=12.8;
41 const Double_t AliITSRecoParam::fgkdr=0.03;
42 const Double_t AliITSRecoParam::fgkzm=0.2;
43 const Double_t AliITSRecoParam::fgkdm=0.40;
44 const Double_t AliITSRecoParam::fgkrs=50.0;
45 const Double_t AliITSRecoParam::fgkds=0.001;
46 const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
47 const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
48 const Double_t AliITSRecoParam::fgkrPipe=3.;
49 const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
50 const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
51 const Double_t AliITSRecoParam::fgkdPipe=0.0028;
52 const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
53 const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
54 const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
55 const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
56 const Double_t AliITSRecoParam::fgkX0Air=21.82;
57 const Double_t AliITSRecoParam::fgkX0Be=65.19;
58 const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
59 const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
60 const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
61 const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
62
63 //_____________________________________________________________________________
64 AliITSRecoParam::AliITSRecoParam() : TObject(),
65 fMaxSnp(1.),
66 fNSigmaYLayerForRoadY(0),
67 fNSigmaRoadY(0),
68 fNSigmaZLayerForRoadZ(0),
69 fNSigmaRoadZ(0),
70 fNSigma2RoadZC(0),
71 fNSigma2RoadYC(0),
72 fNSigma2RoadZNonC(0),
73 fNSigma2RoadYNonC(0),
74 fMaxNormChi2NonCForHypothesis(0),
75 fMaxChi2(0),
76 fMaxRoad(0),
77 fMaxChi2In(0),
78 fChi2PerCluster(0),
79 fXV(0), 
80 fYV(0),
81 fZV(0),
82 fSigmaXV(0),
83 fSigmaYV(0),
84 fSigmaZV(0),
85 fVertexCut(0),
86 fMaxDZforPrimTrk(0),
87 fMaxDZToUseConstraint(0), 
88 fMaxDforV0dghtrForProlongation(0),
89 fMaxDForProlongation(0),
90 fMaxDZForProlongation(0),
91 fMinPtForProlongation(0),
92 fAddVirtualClustersInDeadZone(kTRUE),
93 fZWindowDeadZone(0),
94 fSigmaXDeadZoneHit2(0),
95 fSigmaZDeadZoneHit2(0),
96 fXPassDeadZoneHits(0),
97 fUseTGeoInTracker(0),
98 fAllowSharedClusters(kTRUE),
99 fClusterErrorsParam(1),
100 fFindV0s(kTRUE),
101 fFactorSAWindowSizes(1.)
102 {
103   //
104   // constructor
105   //
106   SetLayersParameters();
107   SetUseTGeoInTracker(0);
108   SetAllowSharedClusters(kTRUE);
109   SetFindV0s(kTRUE);
110   SetAddVirtualClustersInDeadZone(kTRUE);
111   SetUseAmplitudeInfo(kTRUE);
112   SetClusterErrorsParam(1);
113 }
114 //_____________________________________________________________________________
115 AliITSRecoParam::~AliITSRecoParam() 
116 {
117   //
118   // destructor
119   //  
120 }
121 //_____________________________________________________________________________
122 AliITSRecoParam *AliITSRecoParam::GetHighFluxParam() 
123 {
124   //
125   // make default reconstruction  parameters for hig  flux env.
126   //
127   AliITSRecoParam *param = new AliITSRecoParam();
128
129   param->fMaxSnp = 0.95;
130
131   param->fNSigmaYLayerForRoadY = 4.;
132   param->fNSigmaRoadY = 7.5;
133   param->fNSigmaZLayerForRoadZ = 4.;
134   param->fNSigmaRoadZ = 7.5;
135
136   param->fNSigma2RoadZC = 60.; //7.75^2
137   param->fNSigma2RoadYC = 60.; //7.75^2
138   param->fNSigma2RoadZNonC = 50.; //7.07^2
139   param->fNSigma2RoadYNonC = 50.; //7.07^2
140
141   param->fMaxChi2PerCluster[0] = 11.; //7
142   param->fMaxChi2PerCluster[1] = 12.; //5
143   param->fMaxChi2PerCluster[2] = 12.; //8
144   param->fMaxChi2PerCluster[3] = 5.;  //8
145   param->fMaxChi2PerCluster[4] = 12.; //6.5
146
147   param->fMaxNormChi2NonC[0] = 7.;
148   param->fMaxNormChi2NonC[1] = 8.;
149   param->fMaxNormChi2NonC[2] = 8.;
150   param->fMaxNormChi2NonC[3] = 11.;
151   param->fMaxNormChi2NonC[4] = 14.;
152   param->fMaxNormChi2NonC[5] = 25.;
153
154   param->fMaxNormChi2C[0] = 11.;
155   param->fMaxNormChi2C[1] = 13.;
156   param->fMaxNormChi2C[2] = 15.;
157   param->fMaxNormChi2C[3] = 18.;
158   param->fMaxNormChi2C[4] = 30.;
159   param->fMaxNormChi2C[5] = 35.;
160
161   param->fMaxNormChi2NonCForHypothesis = 7.;
162   
163   param->fMaxChi2 = 35.;
164
165   param->fMaxChi2s[0] = 25.; //40   
166   param->fMaxChi2s[1] = 25.; //40   
167   param->fMaxChi2s[2] = 25.; //40   
168   param->fMaxChi2s[3] = 25.; //40   
169   param->fMaxChi2s[4] = 40.; //40   
170   param->fMaxChi2s[5] = 50.; //40
171
172   param->fMaxRoad = 6.;
173
174   // not used
175   param->fMaxChi2In = 16.;
176    
177   param->fMaxChi2sR[0] = 10.;   
178   param->fMaxChi2sR[1] = 10.;   
179   param->fMaxChi2sR[2] = 10.;   
180   param->fMaxChi2sR[3] = 10.;   
181   param->fMaxChi2sR[4] = 30.;   
182   param->fMaxChi2sR[5] = 40.;   
183
184   param->fChi2PerCluster = 9.;
185   // not used
186
187   param->fXV = 0.;
188   param->fYV = 0.;
189   param->fZV = 0.;
190   param->fSigmaXV = 0.0050;
191   param->fSigmaYV = 0.0050;
192   param->fSigmaZV = 0.0100;
193
194   param->fVertexCut = 25.;
195
196   param->fMaxDZforPrimTrk = 0.4;
197   param->fMaxDZToUseConstraint = 3.;
198
199   param->fMaxDforV0dghtrForProlongation = 30.;
200   param->fMaxDForProlongation = 10.;
201   param->fMaxDZForProlongation = 20.;
202   param->fMinPtForProlongation = 0.120;
203
204   param->fZWindowDeadZone = 2.0;
205   param->fSigmaXDeadZoneHit2 = 0.004/12.;
206   param->fSigmaZDeadZoneHit2 = 0.001/12.;
207   param->fXPassDeadZoneHits = 0.018;
208
209   
210   return param;
211 }
212 //_____________________________________________________________________________
213 AliITSRecoParam *AliITSRecoParam::GetLowFluxParam() 
214 {
215   //
216   // make default reconstruction  parameters for low  flux env.
217   //
218   return GetHighFluxParam();
219 }
220 //_____________________________________________________________________________
221 AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam() 
222 {
223   //
224   // make default reconstruction  parameters for cosmics
225   //
226   AliITSRecoParam *param = new AliITSRecoParam();
227
228   param->SetFactorSAWindowSizes(3.);
229
230   param->fMaxSnp = 0.95;
231
232   param->fNSigmaYLayerForRoadY = 4.;
233   param->fNSigmaRoadY = 7.5;
234   param->fNSigmaZLayerForRoadZ = 4.;
235   param->fNSigmaRoadZ = 7.5;
236
237   param->fNSigma2RoadZC = 60.; //7.75^2
238   param->fNSigma2RoadYC = 60.; //7.75^2
239   param->fNSigma2RoadZNonC = 50.; //7.07^2
240   param->fNSigma2RoadYNonC = 50.; //7.07^2
241
242   param->fMaxChi2PerCluster[0] = 11.; //7
243   param->fMaxChi2PerCluster[1] = 12.; //5
244   param->fMaxChi2PerCluster[2] = 12.; //8
245   param->fMaxChi2PerCluster[3] = 5.;  //8
246   param->fMaxChi2PerCluster[4] = 12.; //6.5
247
248   param->fMaxNormChi2NonC[0] = 7.;
249   param->fMaxNormChi2NonC[1] = 8.;
250   param->fMaxNormChi2NonC[2] = 8.;
251   param->fMaxNormChi2NonC[3] = 11.;
252   param->fMaxNormChi2NonC[4] = 14.;
253   param->fMaxNormChi2NonC[5] = 25.;
254
255   param->fMaxNormChi2C[0] = 11.;
256   param->fMaxNormChi2C[1] = 13.;
257   param->fMaxNormChi2C[2] = 15.;
258   param->fMaxNormChi2C[3] = 18.;
259   param->fMaxNormChi2C[4] = 30.;
260   param->fMaxNormChi2C[5] = 35.;
261
262   param->fMaxNormChi2NonCForHypothesis = 7.;
263   
264   param->fMaxChi2 = 35.;
265
266   param->fMaxChi2s[0] = 25.; //40   
267   param->fMaxChi2s[1] = 25.; //40   
268   param->fMaxChi2s[2] = 25.; //40   
269   param->fMaxChi2s[3] = 25.; //40   
270   param->fMaxChi2s[4] = 40.; //40   
271   param->fMaxChi2s[5] = 50.; //40
272
273   param->fMaxRoad = 6.;
274
275   // not used
276   param->fMaxChi2In = 16.;
277    
278   param->fMaxChi2sR[0] = 10.;   
279   param->fMaxChi2sR[1] = 10.;   
280   param->fMaxChi2sR[2] = 10.;   
281   param->fMaxChi2sR[3] = 10.;   
282   param->fMaxChi2sR[4] = 30.;   
283   param->fMaxChi2sR[5] = 40.;   
284
285   param->fChi2PerCluster = 9.;
286   // not used
287
288   param->fXV = 0.;
289   param->fYV = 0.;
290   param->fZV = 0.;
291   param->fSigmaXV = 0.0050;
292   param->fSigmaYV = 0.0050;
293   param->fSigmaZV = 0.0100;
294
295   param->fVertexCut = 25.;
296
297   param->fMaxDZforPrimTrk = 0.4;
298   param->fMaxDZToUseConstraint = 3.;
299
300   param->fMaxDforV0dghtrForProlongation = 30.;
301   param->fMaxDForProlongation = 10.;
302   param->fMaxDZForProlongation = 20.;
303   param->fMinPtForProlongation = 0.120;
304
305   param->fZWindowDeadZone = 2.0;
306   param->fSigmaXDeadZoneHit2 = 0.004/12.;
307   param->fSigmaZDeadZoneHit2 = 0.001/12.;
308   param->fXPassDeadZoneHits = 0.018;
309
310   
311   return param;
312 }
313 //_____________________________________________________________________________
314 void AliITSRecoParam::SetLayersParameters() 
315 {
316   //
317   // number of layers and layers spatial resolutions
318   //
319
320   // spatial resolutions of the detectors
321   // y: 12 12 38 38 20 20 micron
322   fSigmaY2[0]=1.44e-6;
323   fSigmaY2[1]=1.44e-6;
324   fSigmaY2[2]=1.444e-5;
325   fSigmaY2[3]=1.444e-5;
326   fSigmaY2[4]=4.0e-6;
327   fSigmaY2[5]=4.0e-6;
328   // z: 120 120 28 28 830 830 micron
329   fSigmaZ2[0]=1.44e-4;
330   fSigmaZ2[1]=1.44e-4;
331   fSigmaZ2[2]=7.84e-6;
332   fSigmaZ2[3]=7.84e-6;
333   fSigmaZ2[4]=6.889e-3;
334   fSigmaZ2[5]=6.889e-3;
335
336   return;
337 }
338