/[PAMELA software]/DarthVader/TrackerLevel2/inc/TrkLevel2.h
ViewVC logotype

Annotation of /DarthVader/TrackerLevel2/inc/TrkLevel2.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.42 - (hide annotations) (download)
Wed Mar 11 14:19:09 2009 UTC (15 years, 9 months ago) by pam-fi
Branch: MAIN
CVS Tags: v9r00, v9r01
Changes since 1.41: +4 -1 lines
File MIME type: text/plain
new methods to define the acceptance, with tolerance.

1 mocchiut 1.1 /**
2     * \file TrkLevel2.h
3     * \author Elena Vannuccini
4     */
5     #ifndef trklevel2_h
6     #define trklevel2_h
7    
8     #include <TObject.h>
9     #include <TObjArray.h>
10     #include <TClonesArray.h>
11 pam-fi 1.7 #include <TRefArray.h>
12 pam-fi 1.8 #include <TRef.h>
13 pam-fi 1.3
14 pam-fi 1.18 #include <TrkParams.h>
15 pam-fi 1.8 #include <TrkLevel1.h>
16 mocchiut 1.1
17 pam-fi 1.2 // z-coordinate of track state-vector reference-plane
18     #define ZINI 23.5
19 pam-fi 1.5 // (mechanical) z-coordinate of the tracker planes
20 pam-fi 1.27 #define ZTRK6 -22.22
21     #define ZTRK5 -13.31
22     #define ZTRK4 -4.41
23     #define ZTRK3 4.49
24     #define ZTRK2 13.39
25     #define ZTRK1 22.29
26     // magnet cavity dimensions
27     #define ZMAGNHIGH 21.83
28     #define ZMAGNLOW -21.83
29     #define XMAGNHIGH 8.07
30     #define XMAGNLOW -8.07
31     #define YMAGNHIGH 6.57
32     #define YMAGNLOW -6.57
33 pam-fi 1.35 // tof planes
34     #define ZS11 53.74
35     #define ZS12 53.04
36     #define ZS21 23.94
37     #define ZS22 23.44
38     #define ZS31 -23.49
39     #define ZS32 -24.34
40    
41 pam-fi 1.5 // (mechanical) x/y-coordinates of magnet cavity
42 pam-fi 1.34 /* #define XTRKL -8.1 */
43     /* #define XTRKR 8.1 */
44     /* #define YTRKL -6.6 */
45     /* #define YTRKR 6.6 */
46 pam-fi 1.2
47 mocchiut 1.1 /**
48     * \brief Class to describe, by points, a particle trajectory in the apparatus.
49     *
50     * The idea is to create it by integrating the equations of motion, given the
51     * track state vector and the z coordinates where to evaluate track position.
52     */
53     // ==================================================================
54     class Trajectory : public TObject{
55     private:
56    
57     public:
58    
59     int npoint; ///< number of evaluated points along the trajectory
60 pam-fi 1.35 float* x; //[npoint]
61     float* y; //[npoint]
62     float* z; //[npoint]
63     float* thx; //[npoint]
64     float* thy; //[npoint]
65     float* tl; //[npoint]
66 mocchiut 1.1
67 pam-fi 1.2 Trajectory();
68 mocchiut 1.1 Trajectory(int n);
69     Trajectory(int n, float* pz);
70 pam-fi 1.15 ~Trajectory(){Delete();};
71 mocchiut 1.1 void Dump();
72 pam-fi 1.15 void Delete();
73 mocchiut 1.1
74 pam-fi 1.41 int DoTrack(float* al, float zini);
75     int DoTrack(float* al){ return DoTrack(al,23.5); };
76    
77     int DoTrack2(float* al, float zini);
78     int DoTrack2(float* al){ return DoTrack2(al,23.5); };
79    
80 pam-fi 1.2 float GetLength(){float l=0; for(int i=0; i<npoint;i++)l=l+tl[i]; return l;};
81     float GetLength(int,int);
82    
83 pam-fi 1.35 ClassDef(Trajectory,3);
84 mocchiut 1.1
85     };
86     /**
87     * \brief Class to describe fitted tracks.
88     *
89     * A track is defined by the measured coordinates associated to it, the
90     * track status vector, plus other quantities.
91     * A track may have an "image", due to the ambiguity in the y view.
92 pam-fi 1.24 *
93     * Cluster flags: xgood[6], ygood[6]
94     *
95     * xgood/ygood = +/- 0lsccccccc
96 pam-fi 1.35 * ccccccc ID (1-7483647) of the included cluster
97     * s sensor number (1,2 - increasing y)
98     * l ladder number (1,2,3 - increasing x)
99     * +/- does-not/does include bad strips
100     *
101 mocchiut 1.1 */
102     // ==================================================================
103     class TrkTrack : public TObject {
104    
105     private:
106    
107 pam-fi 1.32 public:
108    
109 pam-fi 1.3 int seqno; ///<stored track sequential number
110     int image; ///<sequential number of track-image
111 pam-fi 1.8
112 pam-fi 1.22 float al[5]; ///<TRACK STATE VECTOR
113 mocchiut 1.1 float coval[5][5]; ///<covariance matrix
114 pam-fi 1.31 int xgood[6]; ///<cluster id for x-view (0 = view not included in the fit)
115     int ygood[6]; ///<cluster id for y-view (0 = view not included in the fit)
116 mocchiut 1.1 float xm[6]; ///<measured x coordinates
117     float ym[6]; ///<measured y coordinates
118     float zm[6]; ///<measured z coordinates
119     float resx[6]; ///<spatial resolution on X view
120     float resy[6]; ///<spatial resolution on y view
121 pam-fi 1.24 float tailx[6]; ///<spatial resolution tail on X view
122     float taily[6]; ///<spatial resolution tail on y view
123 mocchiut 1.1 float chi2; ///<chi2
124 pam-fi 1.31 int nstep; ///<n.step
125 pam-fi 1.12 float xv[6]; ///<calculated x coordinates
126 mocchiut 1.1 float yv[6]; ///<calculated y coordinates
127     float zv[6]; ///<calculated z coordinates
128     float axv[6]; ///<calculated angles (deg) on x view
129     float ayv[6]; ///<calculated angles (deg) on y view
130 pam-fi 1.24 float dedx_x[6]; ///<dE/dx in MIP (<0 if saturated)
131     float dedx_y[6]; ///<dE/dx in MIP (<0 if saturated)
132 pam-fi 1.31 int multmaxx[6]; ///<cluster multiplicity and strip of maximum on x view
133     int multmaxy[6]; ///<cluster multiplicity and strip of maximum on y view
134     float seedx[6]; ///< seed of the cluster x
135     float seedy[6]; ///< seed of the cluster y
136     float xpu[6]; ///< x coordinate in pitch units
137     float ypu[6]; ///< y coordinate in pitch units
138 pam-fi 1.3
139 pam-fi 1.35 float xGF[14]; ///<calculated x coordinates on GF reference planes
140     float yGF[14]; ///<calculated y coordinates on GF reference planes
141    
142 mocchiut 1.1 TrkTrack();
143     TrkTrack(const TrkTrack&);
144    
145 pam-fi 1.15 ~TrkTrack(){ Delete(); };
146 pam-fi 1.10
147 mocchiut 1.1 void Dump();
148 pam-fi 1.12 void Clear();
149 pam-fi 1.15 void Clear(Option_t *option){Clear();};
150 pam-fi 1.12 void Delete();
151 pam-fi 1.15 void Copy(TrkTrack&);
152 pam-fi 1.16 // void Set();
153    
154 pam-fi 1.3 Int_t GetSeqNo(){return seqno;} ///< Returns the track sequential number
155     Int_t GetImageSeqNo(){return image;} ///< Returns the track image sequential number
156 mocchiut 1.1 Bool_t HasImage(){return !(image==-1);} ///< Returns true if the track has an image
157 pam-fi 1.35 int DoTrack(Trajectory* t); ///< Evaluates the trajectory in the apparatus.
158     int DoTrack2(Trajectory* t); ///< Evaluates the trajectory in the apparatus.
159     float BdL(){return 0;}; ///< Evaluates the integral of B*dL along the track.
160 pam-fi 1.24 Int_t GetNX(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)XGood(i); return n;};
161     Int_t GetNY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i); return n;};
162 pam-fi 1.39 Int_t GetNXY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i)*XGood(i); return n;};
163 pam-fi 1.3 Int_t GetNtot(){return GetNX()+GetNY();};
164 mocchiut 1.1 Float_t GetRigidity();
165     Float_t GetDeflection();
166 pam-fi 1.24 Bool_t IsSaturated(int,int);
167     Bool_t IsSaturated(int);
168     Bool_t IsSaturated();
169     Bool_t IsBad(int,int);
170 mocchiut 1.1 Float_t GetDEDX();
171 pam-fi 1.28 Float_t GetDEDX(int ip);
172     Float_t GetDEDX(int ip,int iv);
173 pam-fi 1.39 Int_t GetLeverArmXY();
174 pam-fi 1.24 Int_t GetLeverArmX();
175     Int_t GetLeverArmY();
176 pam-fi 1.29 Float_t GetChi2X();
177     Float_t GetChi2Y();
178     Float_t GetLnLX();
179     Float_t GetLnLY();
180 pam-fi 1.12
181 pam-fi 1.30 Float_t GetEffectiveAngle(int ip, int iv);
182    
183 pam-fi 1.12 void SetMeasure(double *xmeas, double *ymeas, double *zmeas);
184     void SetResolution(double *rx, double *ry);
185 pam-fi 1.26 void SetTail(double *tx, double *ty, double factor);
186     void SetStudentParam(int flag);
187 pam-fi 1.12 void SetGood(int *xg, int *yg);
188     void LoadField(TString s);
189 pam-fi 1.25 void Fit(double pfixed, int& fail, int iprint, int froml1);
190     void Fit(double pfixed, int& fail, int iprint){ Fit(pfixed,fail,iprint,0); };
191 pam-fi 1.12 void FitReset();
192 pam-fi 1.19 void SetTrackingMode(int trackmode);
193 pam-fi 1.21 void SetPrecisionFactor(double fact);
194     void SetStepMin(int istepmin);
195 pam-fi 1.33 void SetDeltaB(int id, double db);
196    
197 pam-fi 1.35 Bool_t IsInsideCavity(float);
198     Bool_t IsInsideCavity(){ return IsInsideCavity(0.); };
199 pam-fi 1.42 Bool_t IsInsideAcceptance(float);
200     Bool_t IsInsideAcceptance(){ return IsInsideAcceptance(0.); };
201     Bool_t IsInsideGFSurface(const char*,float);
202     Bool_t IsInsideGFSurface(const char* surf){ return IsInsideGFSurface(surf,0.); };
203 pam-fi 1.14
204 pam-fi 1.28 Bool_t EvaluateClusterPositions();
205 pam-fi 1.25
206 pam-fi 1.14 void FillMiniStruct(cMini2track&);
207     void SetFromMiniStruct(cMini2track*);
208 pam-fi 1.12
209 pam-fi 1.24 Int_t GetClusterX_ID(int ip);
210     Int_t GetClusterY_ID(int ip);
211     Int_t GetLadder(int ip);
212     Int_t GetSensor(int ip);
213     Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; };
214     Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; };
215 pam-fi 1.39 void ResetXGood(int ip){ xgood[ip]=0; };
216 pam-fi 1.25 void ResetYGood(int ip){ ygood[ip]=0; };
217 pam-fi 1.36 /* void SetXGood(int ip, int clid, int is); */
218     /* void SetYGood(int ip, int clid, int is); */
219     void SetXGood(int ip, int clid, int il, int is, bool bad);
220     void SetYGood(int ip, int clid, int il, int is, bool bad);
221     void SetXGood(int ip, int clid, int il, int is){ SetXGood(ip,clid,il,is,false); };
222     void SetYGood(int ip, int clid, int il, int is){ SetYGood(ip,clid,il,is,false); };
223    
224 pam-fi 1.24
225     Bool_t BadClusterX(int ip){ return IsBad(ip,0); };
226     Bool_t BadClusterY(int ip){ return IsBad(ip,1); };
227    
228     Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };
229     Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };
230 pam-fi 1.20
231 pam-fi 1.31 Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };
232     Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };
233     Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };
234     Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };
235     Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };
236     Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };
237 pam-fi 1.36 /* Float_t GetClusterX_oordinatePU(int ip); */
238 pam-fi 1.31 /* Float_t GetClusterY_CoordinatePU(int ip); */
239    
240 pam-fi 1.34 Float_t GetYav();
241     Float_t GetXav();
242     Float_t GetZav();
243    
244     Int_t GetNColumns();
245    
246     Float_t GetDEDX_max(int ip, int iv);
247     Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };
248     Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };
249     Float_t GetDEDX_min(int ip, int iv);
250     Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };
251     Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };
252    
253     Float_t GetResidual_max(int ip, int iv);
254     Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };
255     Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };
256 pam-fi 1.37 Float_t GetResidual_av(int ip, int iv);
257     Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); };
258     Float_t GetResidual_av(){ return GetResidual_av(-1,-1); };
259 pam-fi 1.34
260     Int_t GetClusterX_Multiplicity_max();
261     Int_t GetClusterX_Multiplicity_min();
262     Int_t GetClusterY_Multiplicity_max();
263     Int_t GetClusterY_Multiplicity_min();
264    
265     Float_t GetClusterX_Seed_min();
266     Float_t GetClusterY_Seed_min();
267 pam-fi 1.31
268 mocchiut 1.1 TrkTrack* GetTrkTrack(){return this;};
269    
270 pam-fi 1.3 friend class TrkLevel2;
271    
272 pam-fi 1.35 ClassDef(TrkTrack,5);
273 mocchiut 1.1
274     };
275     /**
276     * \brief Class to describe single clusters ("singlets").
277     *
278     * Single clusters are clusters not associated to any track.
279     */
280     class TrkSinglet : public TObject {
281    
282     private:
283 pam-fi 1.8
284 mocchiut 1.1
285     public:
286 pam-fi 1.8
287 mocchiut 1.1 int plane; ///<plane
288     float coord[2]; ///<coordinate (on sensor 1 and 2)
289 pam-fi 1.24 float sgnl; ///<cluster signal in MIP (<0 if saturated)
290 pam-fi 1.35 int multmax; ///<cluster multiplicity and strip of maximum
291 mocchiut 1.1
292     TrkSinglet();
293     TrkSinglet(const TrkSinglet&);
294 pam-fi 1.15 ~TrkSinglet(){Delete();};
295 mocchiut 1.1
296     void Dump();
297 pam-fi 1.15 void Clear();
298     void Clear(Option_t *option){Clear();};
299     void Delete(){Clear();};
300 pam-fi 1.24 Float_t GetSignal(){return fabs(sgnl);}
301     Bool_t IsSaturated(){return (sgnl<0); };
302 pam-fi 1.35
303     Bool_t IsBad() { return multmax<=0; };
304     Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };
305     Int_t GetCluster_MaxStrip() { return (Int_t)(abs(multmax)%10000); };
306    
307    
308 pam-fi 1.3 friend class TrkLevel2;
309    
310 pam-fi 1.35 ClassDef(TrkSinglet,4);
311 mocchiut 1.1
312     };
313    
314     /**
315     * \brief Class to describe tracker LEVEL2 data.
316     *
317     * A tracker events is defined by some general variables, plus the collection of all the fitted tracks and all
318     * single clusters on X and Y views.
319     * Tracks and single clusters ("singlets") are described by the classes TrkTrack and TrkSinglet respectivelly.
320     *
321     * Each track may have an "image", due to the ambiguity on the Y view, which is stored also.
322     * Thus, the number of stored tracks ( ntrk() ) differs from the number of "physical" tracks ( GetNTracks() ).
323     * Proper methods allow to sort tracks and select the physical ones ( GetTracks() ).
324 pam-fi 1.28 *
325     * The event status indicates the processing status of data from each DSP, according to the following
326     * notation:
327     *
328 pam-fi 1.35 * LSB --> 0 missing packet
329     * 1 CRC error
330     * 2 on-line software alarm (latch-up, timeout ecc...)
331     * 3 jump in the trigger counter
332     * 4 decode error
333     * 5 n.clusters > maximum number (level1 processing)
334     * 6
335     * 7
336     * 8 n.clusters > maximum value (level2 processing)
337     * 9 n.couples per plane > maximum values (vector dimention)
338     * 10 n.doublets > maximum values
339     * 11 n.triplets > maximum values
340     * 12 n.yz-clouds > maximum values
341     * 13 n.xz-clouds > maximum values
342     * 14 n.candidate-tracks > maximum values
343     * 15 n.couples per plane > maximum values (for Hough transform)
344     * MSB --> 16
345 pam-fi 1.28 *
346     *
347 pam-fi 1.39 * For all data processed before June 2007 the event status was coded according to
348     * a different rule:
349     *
350     * Status of level1 processing
351     * 0 -- OK
352     * 1 -- missing packet
353     * 2 -- 1 CRC error
354     * 3 -- 2 on-line software alarm (latch-up flags asserted or n.transmitted-words = 0)
355     * 4 -- 3 jump in the trigger counter
356     * 10 -- 4 decode error
357     * 11 -- 5 n.clusters > maximum number (for level1 processing)
358     * Status of level2 processing
359     * 21 -- 0 n.clusters > maximum value (for level2 processing)
360     * 22 -- 1 n.couples per plane > maximum values (vector dimention)
361     * 23 -- 2 n.doublets > maximum values
362     * 24 -- 3 n.triplets > maximum values
363     * 25 -- 4 n.yz-clouds > maximum values
364     * 26 -- 5 n.xz-clouds > maximum values
365     * 27 -- 6 n.candidate-tracks > maximum values
366     * 28 -- 7 n.couples per plane > maximum values (for Hough transform)
367 pam-fi 1.28 *
368     *
369 mocchiut 1.1 */
370     class TrkLevel2 : public TObject {
371    
372     private:
373 pam-fi 1.15
374 mocchiut 1.1 public:
375    
376 pam-fi 1.15 Int_t good[12]; ///< event status
377 pam-fi 1.24 UInt_t VKmask[12]; ///< Viking-chip mask
378     UInt_t VKflag[12]; ///< Viking-chip flag
379 mocchiut 1.1
380     TClonesArray *Track; ///< fitted tracks
381     TClonesArray *SingletX; ///< x singlets
382     TClonesArray *SingletY; ///< y singlets
383    
384     TrkLevel2();
385     // TrkLevel2(cTrkLevel2 *);
386 pam-fi 1.11 ~TrkLevel2(){Delete();};
387 pam-fi 1.10
388 pam-fi 1.11 void Clear();
389 pam-fi 1.15 void Clear(Option_t *option){Clear();};
390 pam-fi 1.11 void Delete();
391 pam-fi 1.16 void Set();
392 pam-fi 1.40 int UnpackError(){ for(int i=0; i<12; i++)if(!StatusCheck(i,0x12))return 1; return 0;};
393 pam-fi 1.11
394     int ntrk() {return Track->GetEntries();} ///< number of stored track
395 mocchiut 1.1 int nclsx(){return SingletX->GetEntries();} ///< number of x singlets
396     int nclsy(){return SingletY->GetEntries();} ///< number of y singlets
397    
398     void Dump();
399 pam-fi 1.11 void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);
400 pam-fi 1.18 void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL); };
401     void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); };
402     void SetFromLevel2Struct() { SetFromLevel2Struct(&level2event_); };
403 pam-fi 1.11 void GetLevel2Struct(cTrkLevel2 *) const;
404 pam-fi 1.3 void LoadField(TString);
405 pam-fi 1.25 float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss)
406     float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)
407     float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)
408 pam-fi 1.6 Float_t GetZTrk(Int_t);
409 pam-fi 1.34 Float_t GetXTrkLeft(){return XMAGNLOW;};
410     Float_t GetXTrkRight(){return XMAGNHIGH;};
411     Float_t GetYTrkLeft(){return YMAGNLOW;};
412     Float_t GetYTrkRight(){return YMAGNHIGH;};
413 pam-fi 1.6
414 pam-fi 1.24 Bool_t IsMaskedVK(int,int);
415     Bool_t GetVKMask(int,int);
416     Bool_t GetVKFlag(int,int);
417    
418 pam-fi 1.6 TrkSinglet *GetSingletX(int);
419     TrkSinglet *GetSingletY(int);
420    
421     TrkTrack *GetStoredTrack(int i);
422 pam-fi 1.3 Int_t GetSeqNo(Int_t i) {return (((TrkTrack *)Track->At(i))->seqno);}; ///< Returns track sequential number
423 pam-fi 1.24
424 pam-fi 1.11 TRefArray *GetTracks_NFitSorted();
425     TRefArray *GetTracks(){return this->GetTracks_NFitSorted();};
426    
427     Int_t GetNTracks();
428     TrkTrack* GetTrack(int i);
429 mocchiut 1.1 TrkTrack* GetTrackImage(int i);
430 pam-fi 1.11
431 pam-fi 1.3 TrkLevel2* GetTrkLevel2(){return this;}
432     TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array
433    
434 pam-fi 1.28 void StatusDump(int view);
435     Bool_t StatusCheck(int view, int flagmask);
436    
437 pam-fi 1.24 ClassDef(TrkLevel2,3);
438 mocchiut 1.1
439     };
440    
441     #endif

  ViewVC Help
Powered by ViewVC 1.1.23