/[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.47 by pam-ts, Wed Oct 15 08:45:51 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.23 //-22.22  //Aprile 2014... trovata differenza con mech_pos.dat
21  #define ZTRK5 -13.31  #define ZTRK5 -13.32 //-13.31  // ...speriamo bene... no comment
22  #define ZTRK4 -4.41  #define ZTRK4 -4.42//-4.41
23  #define ZTRK3 4.49  #define ZTRK3  4.48//4.49
24  #define ZTRK2 13.39  #define ZTRK2 13.38//13.39
25  #define ZTRK1 22.29  #define ZTRK1 22.28//22.29
26  // magnet cavity dimensions  // magnet cavity dimensions
27  #define ZMAGNHIGH 21.83  #define ZMAGNHIGH 21.83
28  #define ZMAGNLOW -21.83  #define ZMAGNLOW -21.83
# 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        Int_t GetNhit() ;
176      Float_t GetRigidity();      Float_t GetRigidity();
177      Float_t GetDeflection();      Float_t GetDeflection();
178      Bool_t IsSaturated(int,int);      Bool_t IsSaturated(int,int);
# Line 164  public: Line 182  public:
182      Float_t GetDEDX();      Float_t GetDEDX();
183      Float_t GetDEDX(int ip);      Float_t GetDEDX(int ip);
184      Float_t GetDEDX(int ip,int iv);      Float_t GetDEDX(int ip,int iv);
185        Int_t GetLeverArmXY();
186      Int_t GetLeverArmX();      Int_t GetLeverArmX();
187      Int_t GetLeverArmY();      Int_t GetLeverArmY();
188      Float_t GetChi2X();      Float_t GetChi2X();
189      Float_t GetChi2Y();      Float_t GetChi2Y();
190      Float_t GetLnLX();      Float_t GetLnLX();
191      Float_t GetLnLY();      Float_t GetLnLY();
192        Float_t GetMDR(){return (coval[4][4]>0 ? 1./sqrt(coval[4][4]) : 0.); };
193    
194      Float_t GetEffectiveAngle(int ip, int iv);      Float_t GetEffectiveAngle(int ip, int iv);
195            
# Line 180  public: Line 200  public:
200      void SetGood(int *xg, int *yg);      void SetGood(int *xg, int *yg);
201      void LoadField(TString s);      void LoadField(TString s);
202      void Fit(double pfixed, int& fail, int iprint, int froml1);      void Fit(double pfixed, int& fail, int iprint, int froml1);
203      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); }
204      void FitReset();      void FitReset();
205      void SetTrackingMode(int trackmode);      void SetTrackingMode(int trackmode);
206      void SetPrecisionFactor(double fact);      void SetPrecisionFactor(double fact);
# Line 188  public: Line 208  public:
208      void SetDeltaB(int id, double db);      void SetDeltaB(int id, double db);
209    
210      Bool_t IsInsideCavity(float);      Bool_t IsInsideCavity(float);
211      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); };      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); }
212      Bool_t IsInsideAcceptance();      Bool_t IsInsideAcceptance(float);
213        Bool_t IsInsideAcceptance(){ return IsInsideAcceptance(0.); }
214        Bool_t IsInsideGFSurface(const char*,float);
215        Bool_t IsInsideGFSurface(const char* surf){ return IsInsideGFSurface(surf,0.); }
216    
217      Bool_t EvaluateClusterPositions();      Bool_t EvaluateClusterPositions();
218    
219      void FillMiniStruct(cMini2track&);      void FillMiniStruct(cMini2track&);
220      void SetFromMiniStruct(cMini2track*);      void SetFromMiniStruct(cMini2track*);
221        void FillMiniStruct(){ extern cMini2track track_; FillMiniStruct(track_); };
222        void SetFromMiniStruct(){extern cMini2track track_; SetFromMiniStruct(&track_);};
223            
224      Int_t GetClusterX_ID(int ip);      Int_t GetClusterX_ID(int ip);
225      Int_t GetClusterY_ID(int ip);      Int_t GetClusterY_ID(int ip);
226      Int_t GetLadder(int ip);      Int_t GetLadder(int ip);
227      Int_t GetSensor(int ip);      Int_t GetSensor(int ip);
228      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; };      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; }
229      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; };      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; }
230      void ResetXGoo(int ip){ xgood[ip]=0; };      void ResetXGood(int ip){ xgood[ip]=0; }
231      void ResetYGood(int ip){ ygood[ip]=0; };      void ResetYGood(int ip){ ygood[ip]=0; }
232  /*     void SetXGood(int ip, int clid, int is); */  /*     void SetXGood(int ip, int clid, int is); */
233  /*     void SetYGood(int ip, int clid, int is); */  /*     void SetYGood(int ip, int clid, int is); */
234      void SetXGood(int ip, int clid, int il, int is, bool bad);      void SetXGood(int ip, int clid, int il, int is, bool bad);
235      void SetYGood(int ip, int clid, int il, int is, bool bad);      void SetYGood(int ip, int clid, int il, int is, bool bad);
236      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); }
237      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); }
238    
239    
240      Bool_t BadClusterX(int ip){ return IsBad(ip,0); };      Bool_t BadClusterX(int ip){ return IsBad(ip,0); }
241      Bool_t BadClusterY(int ip){ return IsBad(ip,1); };      Bool_t BadClusterY(int ip){ return IsBad(ip,1); }
242    
243      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }
244      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }
245    
246      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); }
247      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); }
248      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); }
249      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); }
250      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; }
251      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; }
252  /*     Float_t GetClusterX_oordinatePU(int ip); */  /*     Float_t GetClusterX_oordinatePU(int ip); */
253  /*     Float_t GetClusterY_CoordinatePU(int ip); */  /*     Float_t GetClusterY_CoordinatePU(int ip); */
254            
# Line 234  public: Line 259  public:
259      Int_t GetNColumns();      Int_t GetNColumns();
260    
261      Float_t GetDEDX_max(int ip, int iv);      Float_t GetDEDX_max(int ip, int iv);
262      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); }
263      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); }
264      Float_t GetDEDX_min(int ip, int iv);      Float_t GetDEDX_min(int ip, int iv);
265      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); }
266      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); }
267    
268      Float_t GetResidual_max(int ip, int iv);      Float_t GetResidual_max(int ip, int iv);
269      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); }
270      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); }
271        Float_t GetResidual_av(int ip, int iv);
272        Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); }
273        Float_t GetResidual_av(){ return GetResidual_av(-1,-1); }
274    
275      Int_t GetClusterX_Multiplicity_max();      Int_t GetClusterX_Multiplicity_max();
276      Int_t GetClusterX_Multiplicity_min();      Int_t GetClusterX_Multiplicity_min();
# Line 252  public: Line 280  public:
280      Float_t GetClusterX_Seed_min();      Float_t GetClusterX_Seed_min();
281      Float_t GetClusterY_Seed_min();      Float_t GetClusterY_Seed_min();
282    
283      TrkTrack* GetTrkTrack(){return this;};      TrkTrack* GetTrkTrack(){return this;}
284    
285      friend class TrkLevel2;      friend class TrkLevel2;
286    
287      ClassDef(TrkTrack,5);      ClassDef(TrkTrack,6);
288    
289  };  };
290  /**  /**
# Line 278  public: Line 306  public:
306    
307      TrkSinglet();      TrkSinglet();
308      TrkSinglet(const TrkSinglet&);      TrkSinglet(const TrkSinglet&);
309      ~TrkSinglet(){Delete();};      ~TrkSinglet(){Delete();}
310    
311      void Dump();      void Dump();
312      void Clear();      void Clear();
313      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
314      void Delete(){Clear();};      void Delete(){Clear();};
315      Float_t GetSignal(){return fabs(sgnl);}      Float_t GetSignal(){return fabs(sgnl);}
316      Bool_t IsSaturated(){return (sgnl<0); };      Bool_t IsSaturated(){return (sgnl<0); }
317    
318      Bool_t IsBad()                 { return multmax<=0; };      Bool_t IsBad()                 { return multmax<=0; }
319      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); }
320      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); };      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); }
321    
322    
323      friend class TrkLevel2;      friend class TrkLevel2;
# Line 370  class TrkLevel2 : public TObject { Line 398  class TrkLevel2 : public TObject {
398    
399      TrkLevel2();      TrkLevel2();
400  //    TrkLevel2(cTrkLevel2 *);  //    TrkLevel2(cTrkLevel2 *);
401      ~TrkLevel2(){Delete();};      ~TrkLevel2(){Delete();}
402                    
403      void Clear();      void Clear();
404      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
405      void Delete();      void Delete();
406      void Set();      void Set();
407        void SetTrackArray(TClonesArray *track);///<set pointer to the track array
408        TClonesArray*  GetTrackArray(){return Track;}///< returns pointer to the track array
409        TClonesArray** GetPointerToTrackArray(){return &Track;}///< returns pointer to pointer to the track array
410        int UnpackError(){ for(int i=0; i<12; i++)if(!StatusCheck(i,0x12))return 1; return 0;}
411            
412      int ntrk() {return Track->GetEntries();}    ///< number of stored track      int ntrk() {return Track->GetEntries();}    ///< number of stored track
413      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 415  class TrkLevel2 : public TObject {
415    
416      void Dump();      void Dump();
417      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);
418      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          };      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          }
419      void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); };          void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); }    
420      void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     };          void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     }    
421      void GetLevel2Struct(cTrkLevel2 *) const;      void GetLevel2Struct(cTrkLevel2 *) const;
422      void LoadField(TString);      void LoadField(TString);
423      float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss)      float GetBX(float* v){return TrkParams::GetBX(v);} ///< Bx (kGauss)
424      float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)      float GetBY(float* v){return TrkParams::GetBY(v);} ///< By (kGauss)
425      float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)      float GetBZ(float* v){return TrkParams::GetBZ(v);} ///< Bz (kGauss)
426      Float_t GetZTrk(Int_t);      Float_t GetZTrk(Int_t);
427      Float_t GetXTrkLeft(){return XMAGNLOW;};      Float_t GetXTrkLeft(){return XMAGNLOW;}
428      Float_t GetXTrkRight(){return XMAGNHIGH;};      Float_t GetXTrkRight(){return XMAGNHIGH;}
429      Float_t GetYTrkLeft(){return YMAGNLOW;};      Float_t GetYTrkLeft(){return YMAGNLOW;}
430      Float_t GetYTrkRight(){return YMAGNHIGH;};      Float_t GetYTrkRight(){return YMAGNHIGH;}
431            
432      Bool_t IsMaskedVK(int,int);      Bool_t IsMaskedVK(int,int);
433      Bool_t GetVKMask(int,int);      Bool_t GetVKMask(int,int);
# Line 405  class TrkLevel2 : public TObject { Line 437  class TrkLevel2 : public TObject {
437      TrkSinglet   *GetSingletY(int);      TrkSinglet   *GetSingletY(int);
438            
439      TrkTrack     *GetStoredTrack(int i);      TrkTrack     *GetStoredTrack(int i);
440      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
441    
442      TRefArray *GetTracks_NFitSorted();      TRefArray *GetTracks_NFitSorted();
443      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();};      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();}
444            
445      Int_t     GetNTracks();      Int_t     GetNTracks();
446      TrkTrack* GetTrack(int i);      TrkTrack* GetTrack(int i);
447      TrkTrack* GetTrackImage(int i);      TrkTrack* GetTrackImage(int i);
448            
449      TrkLevel2*    GetTrkLevel2(){return this;}      TrkLevel2*     GetTrkLevel2(){return this;}
     TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array  
       
450      void   StatusDump(int view);      void   StatusDump(int view);
451      Bool_t StatusCheck(int view, int flagmask);      Bool_t StatusCheck(int view, int flagmask);
452    
453      ClassDef(TrkLevel2,3);      ClassDef(TrkLevel2,4);
454    
455  };  };
456    
Legend:
Removed from v.1.36  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.23