/[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.37 by pam-fi, Wed Oct 22 15:17:32 2008 UTC revision 1.45 by pam-fi, Thu Feb 27 11:24:42 2014 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    
217      void FillMiniStruct(cMini2track&);      void FillMiniStruct(cMini2track&);
218      void SetFromMiniStruct(cMini2track*);      void SetFromMiniStruct(cMini2track*);
219        void FillMiniStruct(){ extern cMini2track track_; FillMiniStruct(track_); };
220        void SetFromMiniStruct(){extern cMini2track track_; SetFromMiniStruct(&track_);};
221            
222      Int_t GetClusterX_ID(int ip);      Int_t GetClusterX_ID(int ip);
223      Int_t GetClusterY_ID(int ip);      Int_t GetClusterY_ID(int ip);
224      Int_t GetLadder(int ip);      Int_t GetLadder(int ip);
225      Int_t GetSensor(int ip);      Int_t GetSensor(int ip);
226      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; };      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; }
227      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; };      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; }
228      void ResetXGoo(int ip){ xgood[ip]=0; };      void ResetXGood(int ip){ xgood[ip]=0; }
229      void ResetYGood(int ip){ ygood[ip]=0; };      void ResetYGood(int ip){ ygood[ip]=0; }
230  /*     void SetXGood(int ip, int clid, int is); */  /*     void SetXGood(int ip, int clid, int is); */
231  /*     void SetYGood(int ip, int clid, int is); */  /*     void SetYGood(int ip, int clid, int is); */
232      void SetXGood(int ip, int clid, int il, int is, bool bad);      void SetXGood(int ip, int clid, int il, int is, bool bad);
233      void SetYGood(int ip, int clid, int il, int is, bool bad);      void SetYGood(int ip, int clid, int il, int is, bool bad);
234      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); }
235      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); }
236    
237    
238      Bool_t BadClusterX(int ip){ return IsBad(ip,0); };      Bool_t BadClusterX(int ip){ return IsBad(ip,0); }
239      Bool_t BadClusterY(int ip){ return IsBad(ip,1); };      Bool_t BadClusterY(int ip){ return IsBad(ip,1); }
240    
241      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }
242      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }
243    
244      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); }
245      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); }
246      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); }
247      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); }
248      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; }
249      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; }
250  /*     Float_t GetClusterX_oordinatePU(int ip); */  /*     Float_t GetClusterX_oordinatePU(int ip); */
251  /*     Float_t GetClusterY_CoordinatePU(int ip); */  /*     Float_t GetClusterY_CoordinatePU(int ip); */
252            
# Line 234  public: Line 257  public:
257      Int_t GetNColumns();      Int_t GetNColumns();
258    
259      Float_t GetDEDX_max(int ip, int iv);      Float_t GetDEDX_max(int ip, int iv);
260      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); }
261      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); }
262      Float_t GetDEDX_min(int ip, int iv);      Float_t GetDEDX_min(int ip, int iv);
263      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); }
264      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); }
265    
266      Float_t GetResidual_max(int ip, int iv);      Float_t GetResidual_max(int ip, int iv);
267      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); }
268      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); }
269      Float_t GetResidual_av(int ip, int iv);      Float_t GetResidual_av(int ip, int iv);
270      Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); };      Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); }
271      Float_t GetResidual_av(){ return GetResidual_av(-1,-1); };      Float_t GetResidual_av(){ return GetResidual_av(-1,-1); }
272    
273      Int_t GetClusterX_Multiplicity_max();      Int_t GetClusterX_Multiplicity_max();
274      Int_t GetClusterX_Multiplicity_min();      Int_t GetClusterX_Multiplicity_min();
# Line 255  public: Line 278  public:
278      Float_t GetClusterX_Seed_min();      Float_t GetClusterX_Seed_min();
279      Float_t GetClusterY_Seed_min();      Float_t GetClusterY_Seed_min();
280    
281      TrkTrack* GetTrkTrack(){return this;};      TrkTrack* GetTrkTrack(){return this;}
282    
283      friend class TrkLevel2;      friend class TrkLevel2;
284    
# Line 281  public: Line 304  public:
304    
305      TrkSinglet();      TrkSinglet();
306      TrkSinglet(const TrkSinglet&);      TrkSinglet(const TrkSinglet&);
307      ~TrkSinglet(){Delete();};      ~TrkSinglet(){Delete();}
308    
309      void Dump();      void Dump();
310      void Clear();      void Clear();
311      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
312      void Delete(){Clear();};      void Delete(){Clear();};
313      Float_t GetSignal(){return fabs(sgnl);}      Float_t GetSignal(){return fabs(sgnl);}
314      Bool_t IsSaturated(){return (sgnl<0); };      Bool_t IsSaturated(){return (sgnl<0); }
315    
316      Bool_t IsBad()                 { return multmax<=0; };      Bool_t IsBad()                 { return multmax<=0; }
317      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); }
318      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); };      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); }
319    
320    
321      friend class TrkLevel2;      friend class TrkLevel2;
# Line 373  class TrkLevel2 : public TObject { Line 396  class TrkLevel2 : public TObject {
396    
397      TrkLevel2();      TrkLevel2();
398  //    TrkLevel2(cTrkLevel2 *);  //    TrkLevel2(cTrkLevel2 *);
399      ~TrkLevel2(){Delete();};      ~TrkLevel2(){Delete();}
400                    
401      void Clear();      void Clear();
402      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
403      void Delete();      void Delete();
404      void Set();      void Set();
405      int UnpackError(){ for(int i=0; i<12; i++)if(StatusCheck(i,0x12))return 1; return 0;};      int UnpackError(){ for(int i=0; i<12; i++)if(!StatusCheck(i,0x12))return 1; return 0;}
406            
407      int ntrk() {return Track->GetEntries();}    ///< number of stored track      int ntrk() {return Track->GetEntries();}    ///< number of stored track
408      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets
# Line 387  class TrkLevel2 : public TObject { Line 410  class TrkLevel2 : public TObject {
410    
411      void Dump();      void Dump();
412      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);
413      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          };      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          }
414      void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); };          void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); }    
415      void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     };          void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     }    
416      void GetLevel2Struct(cTrkLevel2 *) const;      void GetLevel2Struct(cTrkLevel2 *) const;
417      void LoadField(TString);      void LoadField(TString);
418      float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss)      float GetBX(float* v){return TrkParams::GetBX(v);} ///< Bx (kGauss)
419      float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)      float GetBY(float* v){return TrkParams::GetBY(v);} ///< By (kGauss)
420      float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)      float GetBZ(float* v){return TrkParams::GetBZ(v);} ///< Bz (kGauss)
421      Float_t GetZTrk(Int_t);      Float_t GetZTrk(Int_t);
422      Float_t GetXTrkLeft(){return XMAGNLOW;};      Float_t GetXTrkLeft(){return XMAGNLOW;}
423      Float_t GetXTrkRight(){return XMAGNHIGH;};      Float_t GetXTrkRight(){return XMAGNHIGH;}
424      Float_t GetYTrkLeft(){return YMAGNLOW;};      Float_t GetYTrkLeft(){return YMAGNLOW;}
425      Float_t GetYTrkRight(){return YMAGNHIGH;};      Float_t GetYTrkRight(){return YMAGNHIGH;}
426            
427      Bool_t IsMaskedVK(int,int);      Bool_t IsMaskedVK(int,int);
428      Bool_t GetVKMask(int,int);      Bool_t GetVKMask(int,int);
# Line 409  class TrkLevel2 : public TObject { Line 432  class TrkLevel2 : public TObject {
432      TrkSinglet   *GetSingletY(int);      TrkSinglet   *GetSingletY(int);
433            
434      TrkTrack     *GetStoredTrack(int i);      TrkTrack     *GetStoredTrack(int i);
435      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
436    
437      TRefArray *GetTracks_NFitSorted();      TRefArray *GetTracks_NFitSorted();
438      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();};      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();}
439            
440      Int_t     GetNTracks();      Int_t     GetNTracks();
441      TrkTrack* GetTrack(int i);      TrkTrack* GetTrack(int i);
442      TrkTrack* GetTrackImage(int i);      TrkTrack* GetTrackImage(int i);
443            
444      TrkLevel2*    GetTrkLevel2(){return this;}      TrkLevel2*    GetTrkLevel2(){return this;}
445      TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array      TClonesArray* GetTrackArray(){return Track;}///< returns pointer to the track array
446            
447      void   StatusDump(int view);      void   StatusDump(int view);
448      Bool_t StatusCheck(int view, int flagmask);      Bool_t StatusCheck(int view, int flagmask);
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.23