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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.36 by pam-fi, Fri Apr 11 13:44:36 2008 UTC revision 1.44 by pam-fi, Fri Sep 2 14:49:02 2011 UTC
# Line 15  Line 15 
15  #include <TrkLevel1.h>  #include <TrkLevel1.h>
16    
17  // z-coordinate of track state-vector reference-plane  // z-coordinate of track state-vector reference-plane
18  #define ZINI 23.5    #define ZINI 23.5 ///< z-coordinate of track state-vector reference-plane.
19  // (mechanical) z-coordinate of the tracker planes  // (mechanical) z-coordinate of the tracker planes
20  #define ZTRK6 -22.22  #define ZTRK6 -22.22
21  #define ZTRK5 -13.31  #define ZTRK5 -13.31
# Line 67  class Trajectory : public TObject{ Line 67  class Trajectory : public TObject{
67      Trajectory();      Trajectory();
68      Trajectory(int n);      Trajectory(int n);
69      Trajectory(int n, float* pz);      Trajectory(int n, float* pz);
70      ~Trajectory(){Delete();};      ~Trajectory(){Delete();}
71      void Dump();      void Dump();
72      void Delete();      void Delete();
73    
74      int DoTrack2(float* al);      int DoTrack(float* al, float zini);
75      float GetLength(){float l=0; for(int i=0; i<npoint;i++)l=l+tl[i]; return l;};      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        float GetLength(){float l=0; for(int i=0; i<npoint;i++)l=l+tl[i]; return l;}
81      float GetLength(int,int);      float GetLength(int,int);
82    
83      ClassDef(Trajectory,3);      ClassDef(Trajectory,3);
# Line 104  public: Line 109  public:
109      int   seqno;           ///<stored track sequential number      int   seqno;           ///<stored track sequential number
110      int   image;           ///<sequential number of track-image      int   image;           ///<sequential number of track-image
111                    
112      float al[5];           ///<TRACK STATE VECTOR      /*! @brief Track state vector.
113         *
114         *  This is the track state vector on reference plane defined by #ZINI.
115         *
116         *  al[0]: X coordinate [cm]
117         *  al[1]: Y coordinate [cm]
118         *  al[2]: sin theta (altitude; theta = 0 is normal incidence)
119         *  al[3]: phi (azimuth; phi = 0 is negative X axis)
120         *  al[4]: deflection (with sign) [1/GV]
121         *
122         */
123        float al[5];
124      float coval[5][5];     ///<covariance matrix      float coval[5][5];     ///<covariance matrix
125      int   xgood[6];        ///<cluster id for x-view (0 = view not included in the fit)      int   xgood[6];        ///<cluster id for x-view (0 = view not included in the fit)
126      int   ygood[6];        ///<cluster id for y-view (0 = view not included in the fit)      int   ygood[6];        ///<cluster id for y-view (0 = view not included in the fit)
# Line 137  public: Line 153  public:
153      TrkTrack();      TrkTrack();
154      TrkTrack(const TrkTrack&);      TrkTrack(const TrkTrack&);
155    
156      ~TrkTrack(){ Delete(); };      ~TrkTrack(){ Delete(); }
157                    
158      void Dump();      void Dump();
159      void Clear();      void Clear();
160      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
161      void Delete();      void Delete();
162      void Copy(TrkTrack&);      void Copy(TrkTrack&);
163  //    void Set();  //    void Set();
# Line 151  public: Line 167  public:
167      Bool_t HasImage(){return !(image==-1);} ///< Returns true if the track has an image      Bool_t HasImage(){return !(image==-1);} ///< Returns true if the track has an image
168      int DoTrack(Trajectory* t);             ///< Evaluates the trajectory in the apparatus.      int DoTrack(Trajectory* t);             ///< Evaluates the trajectory in the apparatus.
169      int DoTrack2(Trajectory* t);            ///< Evaluates the trajectory in the apparatus.      int DoTrack2(Trajectory* t);            ///< Evaluates the trajectory in the apparatus.
170      float BdL(){return 0;};                 ///< Evaluates the integral of B*dL along the track.      float BdL(){return 0;}              ///< Evaluates the integral of B*dL along the track.
171      Int_t GetNX(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)XGood(i); return n;};      Int_t GetNX(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)XGood(i); return n;}
172      Int_t GetNY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i); return n;};      Int_t GetNY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i); return n;}
173      Int_t GetNtot(){return GetNX()+GetNY();};      Int_t GetNXY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i)*XGood(i); return n;}
174        Int_t GetNtot(){return GetNX()+GetNY();}
175      Float_t GetRigidity();      Float_t GetRigidity();
176      Float_t GetDeflection();      Float_t GetDeflection();
177      Bool_t IsSaturated(int,int);      Bool_t IsSaturated(int,int);
# Line 164  public: Line 181  public:
181      Float_t GetDEDX();      Float_t GetDEDX();
182      Float_t GetDEDX(int ip);      Float_t GetDEDX(int ip);
183      Float_t GetDEDX(int ip,int iv);      Float_t GetDEDX(int ip,int iv);
184        Int_t GetLeverArmXY();
185      Int_t GetLeverArmX();      Int_t GetLeverArmX();
186      Int_t GetLeverArmY();      Int_t GetLeverArmY();
187      Float_t GetChi2X();      Float_t GetChi2X();
# Line 180  public: Line 198  public:
198      void SetGood(int *xg, int *yg);      void SetGood(int *xg, int *yg);
199      void LoadField(TString s);      void LoadField(TString s);
200      void Fit(double pfixed, int& fail, int iprint, int froml1);      void Fit(double pfixed, int& fail, int iprint, int froml1);
201      void Fit(double pfixed, int& fail, int iprint){ Fit(pfixed,fail,iprint,0); };      void Fit(double pfixed, int& fail, int iprint){ Fit(pfixed,fail,iprint,0); }
202      void FitReset();      void FitReset();
203      void SetTrackingMode(int trackmode);      void SetTrackingMode(int trackmode);
204      void SetPrecisionFactor(double fact);      void SetPrecisionFactor(double fact);
# Line 188  public: Line 206  public:
206      void SetDeltaB(int id, double db);      void SetDeltaB(int id, double db);
207    
208      Bool_t IsInsideCavity(float);      Bool_t IsInsideCavity(float);
209      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); };      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); }
210      Bool_t IsInsideAcceptance();      Bool_t IsInsideAcceptance(float);
211        Bool_t IsInsideAcceptance(){ return IsInsideAcceptance(0.); }
212        Bool_t IsInsideGFSurface(const char*,float);
213        Bool_t IsInsideGFSurface(const char* surf){ return IsInsideGFSurface(surf,0.); }
214    
215      Bool_t EvaluateClusterPositions();      Bool_t EvaluateClusterPositions();
216    
# Line 200  public: Line 221  public:
221      Int_t GetClusterY_ID(int ip);      Int_t GetClusterY_ID(int ip);
222      Int_t GetLadder(int ip);      Int_t GetLadder(int ip);
223      Int_t GetSensor(int ip);      Int_t GetSensor(int ip);
224      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; };      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; }
225      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; };      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; }
226      void ResetXGoo(int ip){ xgood[ip]=0; };      void ResetXGood(int ip){ xgood[ip]=0; }
227      void ResetYGood(int ip){ ygood[ip]=0; };      void ResetYGood(int ip){ ygood[ip]=0; }
228  /*     void SetXGood(int ip, int clid, int is); */  /*     void SetXGood(int ip, int clid, int is); */
229  /*     void SetYGood(int ip, int clid, int is); */  /*     void SetYGood(int ip, int clid, int is); */
230      void SetXGood(int ip, int clid, int il, int is, bool bad);      void SetXGood(int ip, int clid, int il, int is, bool bad);
231      void SetYGood(int ip, int clid, int il, int is, bool bad);      void SetYGood(int ip, int clid, int il, int is, bool bad);
232      void SetXGood(int ip, int clid, int il, int is){ SetXGood(ip,clid,il,is,false); };      void SetXGood(int ip, int clid, int il, int is){ SetXGood(ip,clid,il,is,false); }
233      void SetYGood(int ip, int clid, int il, int is){ SetYGood(ip,clid,il,is,false); };      void SetYGood(int ip, int clid, int il, int is){ SetYGood(ip,clid,il,is,false); }
234    
235    
236      Bool_t BadClusterX(int ip){ return IsBad(ip,0); };      Bool_t BadClusterX(int ip){ return IsBad(ip,0); }
237      Bool_t BadClusterY(int ip){ return IsBad(ip,1); };      Bool_t BadClusterY(int ip){ return IsBad(ip,1); }
238    
239      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }
240      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }
241    
242      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); }
243      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); }
244      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); }
245      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); }
246      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; }
247      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; }
248  /*     Float_t GetClusterX_oordinatePU(int ip); */  /*     Float_t GetClusterX_oordinatePU(int ip); */
249  /*     Float_t GetClusterY_CoordinatePU(int ip); */  /*     Float_t GetClusterY_CoordinatePU(int ip); */
250            
# Line 234  public: Line 255  public:
255      Int_t GetNColumns();      Int_t GetNColumns();
256    
257      Float_t GetDEDX_max(int ip, int iv);      Float_t GetDEDX_max(int ip, int iv);
258      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); }
259      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); }
260      Float_t GetDEDX_min(int ip, int iv);      Float_t GetDEDX_min(int ip, int iv);
261      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); }
262      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); }
263    
264      Float_t GetResidual_max(int ip, int iv);      Float_t GetResidual_max(int ip, int iv);
265      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); }
266      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); }
267        Float_t GetResidual_av(int ip, int iv);
268        Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); }
269        Float_t GetResidual_av(){ return GetResidual_av(-1,-1); }
270    
271      Int_t GetClusterX_Multiplicity_max();      Int_t GetClusterX_Multiplicity_max();
272      Int_t GetClusterX_Multiplicity_min();      Int_t GetClusterX_Multiplicity_min();
# Line 252  public: Line 276  public:
276      Float_t GetClusterX_Seed_min();      Float_t GetClusterX_Seed_min();
277      Float_t GetClusterY_Seed_min();      Float_t GetClusterY_Seed_min();
278    
279      TrkTrack* GetTrkTrack(){return this;};      TrkTrack* GetTrkTrack(){return this;}
280    
281      friend class TrkLevel2;      friend class TrkLevel2;
282    
# Line 278  public: Line 302  public:
302    
303      TrkSinglet();      TrkSinglet();
304      TrkSinglet(const TrkSinglet&);      TrkSinglet(const TrkSinglet&);
305      ~TrkSinglet(){Delete();};      ~TrkSinglet(){Delete();}
306    
307      void Dump();      void Dump();
308      void Clear();      void Clear();
309      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
310      void Delete(){Clear();};      void Delete(){Clear();};
311      Float_t GetSignal(){return fabs(sgnl);}      Float_t GetSignal(){return fabs(sgnl);}
312      Bool_t IsSaturated(){return (sgnl<0); };      Bool_t IsSaturated(){return (sgnl<0); }
313    
314      Bool_t IsBad()                 { return multmax<=0; };      Bool_t IsBad()                 { return multmax<=0; }
315      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); }
316      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); };      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); }
317    
318    
319      friend class TrkLevel2;      friend class TrkLevel2;
# Line 370  class TrkLevel2 : public TObject { Line 394  class TrkLevel2 : public TObject {
394    
395      TrkLevel2();      TrkLevel2();
396  //    TrkLevel2(cTrkLevel2 *);  //    TrkLevel2(cTrkLevel2 *);
397      ~TrkLevel2(){Delete();};      ~TrkLevel2(){Delete();}
398                    
399      void Clear();      void Clear();
400      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
401      void Delete();      void Delete();
402      void Set();      void Set();
403        int UnpackError(){ for(int i=0; i<12; i++)if(!StatusCheck(i,0x12))return 1; return 0;}
404            
405      int ntrk() {return Track->GetEntries();}    ///< number of stored track      int ntrk() {return Track->GetEntries();}    ///< number of stored track
406      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets
# Line 383  class TrkLevel2 : public TObject { Line 408  class TrkLevel2 : public TObject {
408    
409      void Dump();      void Dump();
410      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);
411      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          };      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          }
412      void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); };          void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); }    
413      void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     };          void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     }    
414      void GetLevel2Struct(cTrkLevel2 *) const;      void GetLevel2Struct(cTrkLevel2 *) const;
415      void LoadField(TString);      void LoadField(TString);
416      float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss)      float GetBX(float* v){return TrkParams::GetBX(v);} ///< Bx (kGauss)
417      float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)      float GetBY(float* v){return TrkParams::GetBY(v);} ///< By (kGauss)
418      float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)      float GetBZ(float* v){return TrkParams::GetBZ(v);} ///< Bz (kGauss)
419      Float_t GetZTrk(Int_t);      Float_t GetZTrk(Int_t);
420      Float_t GetXTrkLeft(){return XMAGNLOW;};      Float_t GetXTrkLeft(){return XMAGNLOW;}
421      Float_t GetXTrkRight(){return XMAGNHIGH;};      Float_t GetXTrkRight(){return XMAGNHIGH;}
422      Float_t GetYTrkLeft(){return YMAGNLOW;};      Float_t GetYTrkLeft(){return YMAGNLOW;}
423      Float_t GetYTrkRight(){return YMAGNHIGH;};      Float_t GetYTrkRight(){return YMAGNHIGH;}
424            
425      Bool_t IsMaskedVK(int,int);      Bool_t IsMaskedVK(int,int);
426      Bool_t GetVKMask(int,int);      Bool_t GetVKMask(int,int);
# Line 405  class TrkLevel2 : public TObject { Line 430  class TrkLevel2 : public TObject {
430      TrkSinglet   *GetSingletY(int);      TrkSinglet   *GetSingletY(int);
431            
432      TrkTrack     *GetStoredTrack(int i);      TrkTrack     *GetStoredTrack(int i);
433      Int_t         GetSeqNo(Int_t i)  {return (((TrkTrack *)Track->At(i))->seqno);}; ///< Returns track sequential number      Int_t         GetSeqNo(Int_t i)  {return (((TrkTrack *)Track->At(i))->seqno);} ///< Returns track sequential number
434    
435      TRefArray *GetTracks_NFitSorted();      TRefArray *GetTracks_NFitSorted();
436      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();};      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();}
437            
438      Int_t     GetNTracks();      Int_t     GetNTracks();
439      TrkTrack* GetTrack(int i);      TrkTrack* GetTrack(int i);
440      TrkTrack* GetTrackImage(int i);      TrkTrack* GetTrackImage(int i);
441            
442      TrkLevel2*    GetTrkLevel2(){return this;}      TrkLevel2*    GetTrkLevel2(){return this;}
443      TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array      TClonesArray* GetTrackArray(){return Track;}///< returns pointer to the track array
444            
445      void   StatusDump(int view);      void   StatusDump(int view);
446      Bool_t StatusCheck(int view, int flagmask);      Bool_t StatusCheck(int view, int flagmask);
Legend:
Removed from v.1.36  
changed lines
  Added in v.1.44
  ViewVC Help
Powered by ViewVC 1.1.23