/[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.41 by pam-fi, Tue Feb 3 13:57:14 2009 UTC revision 1.46 by pam-ts, Wed Jun 4 07:57:02 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 DoTrack(float* al, float zini);      int DoTrack(float* al, float zini);
75      int DoTrack(float* al){ return DoTrack(al,23.5); };      int DoTrack(float* al){ return DoTrack(al,23.5); }
76    
77      int DoTrack2(float* al, float zini);      int DoTrack2(float* al, float zini);
78      int DoTrack2(float* al){ return DoTrack2(al,23.5); };      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;};      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 109  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 142  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 156  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 GetNXY(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)YGood(i)*XGood(i); return n;};      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();};      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 187  public: Line 199  public:
199      void SetGood(int *xg, int *yg);      void SetGood(int *xg, int *yg);
200      void LoadField(TString s);      void LoadField(TString s);
201      void Fit(double pfixed, int& fail, int iprint, int froml1);      void Fit(double pfixed, int& fail, int iprint, int froml1);
202      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); }
203      void FitReset();      void FitReset();
204      void SetTrackingMode(int trackmode);      void SetTrackingMode(int trackmode);
205      void SetPrecisionFactor(double fact);      void SetPrecisionFactor(double fact);
# Line 195  public: Line 207  public:
207      void SetDeltaB(int id, double db);      void SetDeltaB(int id, double db);
208    
209      Bool_t IsInsideCavity(float);      Bool_t IsInsideCavity(float);
210      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); };      Bool_t IsInsideCavity(){ return IsInsideCavity(0.); }
211      Bool_t IsInsideAcceptance();      Bool_t IsInsideAcceptance(float);
212        Bool_t IsInsideAcceptance(){ return IsInsideAcceptance(0.); }
213        Bool_t IsInsideGFSurface(const char*,float);
214        Bool_t IsInsideGFSurface(const char* surf){ return IsInsideGFSurface(surf,0.); }
215    
216      Bool_t EvaluateClusterPositions();      Bool_t EvaluateClusterPositions();
217    
218      void FillMiniStruct(cMini2track&);      void FillMiniStruct(cMini2track&);
219      void SetFromMiniStruct(cMini2track*);      void SetFromMiniStruct(cMini2track*);
220        void FillMiniStruct(){ extern cMini2track track_; FillMiniStruct(track_); };
221        void SetFromMiniStruct(){extern cMini2track track_; SetFromMiniStruct(&track_);};
222            
223      Int_t GetClusterX_ID(int ip);      Int_t GetClusterX_ID(int ip);
224      Int_t GetClusterY_ID(int ip);      Int_t GetClusterY_ID(int ip);
225      Int_t GetLadder(int ip);      Int_t GetLadder(int ip);
226      Int_t GetSensor(int ip);      Int_t GetSensor(int ip);
227      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; };      Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; }
228      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; };      Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; }
229      void ResetXGood(int ip){ xgood[ip]=0; };      void ResetXGood(int ip){ xgood[ip]=0; }
230      void ResetYGood(int ip){ ygood[ip]=0; };      void ResetYGood(int ip){ ygood[ip]=0; }
231  /*     void SetXGood(int ip, int clid, int is); */  /*     void SetXGood(int ip, int clid, int is); */
232  /*     void SetYGood(int ip, int clid, int is); */  /*     void SetYGood(int ip, int clid, int is); */
233      void SetXGood(int ip, int clid, int il, int is, bool bad);      void SetXGood(int ip, int clid, int il, int is, bool bad);
234      void SetYGood(int ip, int clid, int il, int is, bool bad);      void SetYGood(int ip, int clid, int il, int is, bool bad);
235      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); }
236      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); }
237    
238    
239      Bool_t BadClusterX(int ip){ return IsBad(ip,0); };      Bool_t BadClusterX(int ip){ return IsBad(ip,0); }
240      Bool_t BadClusterY(int ip){ return IsBad(ip,1); };      Bool_t BadClusterY(int ip){ return IsBad(ip,1); }
241    
242      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }
243      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }
244    
245      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };      Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); }
246      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };      Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); }
247      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };      Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); }
248      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };      Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); }
249      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };      Float_t GetClusterX_Seed(int ip){ return seedx[ip]; }
250      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };      Float_t GetClusterY_Seed(int ip){ return seedy[ip]; }
251  /*     Float_t GetClusterX_oordinatePU(int ip); */  /*     Float_t GetClusterX_oordinatePU(int ip); */
252  /*     Float_t GetClusterY_CoordinatePU(int ip); */  /*     Float_t GetClusterY_CoordinatePU(int ip); */
253            
# Line 241  public: Line 258  public:
258      Int_t GetNColumns();      Int_t GetNColumns();
259    
260      Float_t GetDEDX_max(int ip, int iv);      Float_t GetDEDX_max(int ip, int iv);
261      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };      Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); }
262      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };      Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); }
263      Float_t GetDEDX_min(int ip, int iv);      Float_t GetDEDX_min(int ip, int iv);
264      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };      Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); }
265      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };      Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); }
266    
267      Float_t GetResidual_max(int ip, int iv);      Float_t GetResidual_max(int ip, int iv);
268      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };      Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); }
269      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };      Float_t GetResidual_max(){ return GetResidual_max(-1,-1); }
270      Float_t GetResidual_av(int ip, int iv);      Float_t GetResidual_av(int ip, int iv);
271      Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); };      Float_t GetResidual_av(int iv){ return GetResidual_av(-1,iv); }
272      Float_t GetResidual_av(){ return GetResidual_av(-1,-1); };      Float_t GetResidual_av(){ return GetResidual_av(-1,-1); }
273    
274      Int_t GetClusterX_Multiplicity_max();      Int_t GetClusterX_Multiplicity_max();
275      Int_t GetClusterX_Multiplicity_min();      Int_t GetClusterX_Multiplicity_min();
# Line 262  public: Line 279  public:
279      Float_t GetClusterX_Seed_min();      Float_t GetClusterX_Seed_min();
280      Float_t GetClusterY_Seed_min();      Float_t GetClusterY_Seed_min();
281    
282      TrkTrack* GetTrkTrack(){return this;};      TrkTrack* GetTrkTrack(){return this;}
283    
284      friend class TrkLevel2;      friend class TrkLevel2;
285    
# Line 288  public: Line 305  public:
305    
306      TrkSinglet();      TrkSinglet();
307      TrkSinglet(const TrkSinglet&);      TrkSinglet(const TrkSinglet&);
308      ~TrkSinglet(){Delete();};      ~TrkSinglet(){Delete();}
309    
310      void Dump();      void Dump();
311      void Clear();      void Clear();
312      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
313      void Delete(){Clear();};      void Delete(){Clear();};
314      Float_t GetSignal(){return fabs(sgnl);}      Float_t GetSignal(){return fabs(sgnl);}
315      Bool_t IsSaturated(){return (sgnl<0); };      Bool_t IsSaturated(){return (sgnl<0); }
316    
317      Bool_t IsBad()                 { return multmax<=0; };      Bool_t IsBad()                 { return multmax<=0; }
318      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };      Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); }
319      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); };      Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); }
320    
321    
322      friend class TrkLevel2;      friend class TrkLevel2;
# Line 380  class TrkLevel2 : public TObject { Line 397  class TrkLevel2 : public TObject {
397    
398      TrkLevel2();      TrkLevel2();
399  //    TrkLevel2(cTrkLevel2 *);  //    TrkLevel2(cTrkLevel2 *);
400      ~TrkLevel2(){Delete();};      ~TrkLevel2(){Delete();}
401                    
402      void Clear();      void Clear();
403      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();}
404      void Delete();      void Delete();
405      void Set();      void Set();
406      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;}
407            
408      int ntrk() {return Track->GetEntries();}    ///< number of stored track      int ntrk() {return Track->GetEntries();}    ///< number of stored track
409      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets      int nclsx(){return SingletX->GetEntries();} ///< number of x singlets
# Line 394  class TrkLevel2 : public TObject { Line 411  class TrkLevel2 : public TObject {
411    
412      void Dump();      void Dump();
413      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);      void SetFromLevel2Struct(cTrkLevel2 *, TrkLevel1 *);
414      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          };      void SetFromLevel2Struct(cTrkLevel2 *s2){ SetFromLevel2Struct(s2, NULL);          }
415      void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); };          void SetFromLevel2Struct(TrkLevel1 *l1) { SetFromLevel2Struct(&level2event_, l1); }    
416      void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     };          void SetFromLevel2Struct()              { SetFromLevel2Struct(&level2event_);     }    
417      void GetLevel2Struct(cTrkLevel2 *) const;      void GetLevel2Struct(cTrkLevel2 *) const;
418      void LoadField(TString);      void LoadField(TString);
419      float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss)      float GetBX(float* v){return TrkParams::GetBX(v);} ///< Bx (kGauss)
420      float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)      float GetBY(float* v){return TrkParams::GetBY(v);} ///< By (kGauss)
421      float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)      float GetBZ(float* v){return TrkParams::GetBZ(v);} ///< Bz (kGauss)
422      Float_t GetZTrk(Int_t);      Float_t GetZTrk(Int_t);
423      Float_t GetXTrkLeft(){return XMAGNLOW;};      Float_t GetXTrkLeft(){return XMAGNLOW;}
424      Float_t GetXTrkRight(){return XMAGNHIGH;};      Float_t GetXTrkRight(){return XMAGNHIGH;}
425      Float_t GetYTrkLeft(){return YMAGNLOW;};      Float_t GetYTrkLeft(){return YMAGNLOW;}
426      Float_t GetYTrkRight(){return YMAGNHIGH;};      Float_t GetYTrkRight(){return YMAGNHIGH;}
427            
428      Bool_t IsMaskedVK(int,int);      Bool_t IsMaskedVK(int,int);
429      Bool_t GetVKMask(int,int);      Bool_t GetVKMask(int,int);
# Line 416  class TrkLevel2 : public TObject { Line 433  class TrkLevel2 : public TObject {
433      TrkSinglet   *GetSingletY(int);      TrkSinglet   *GetSingletY(int);
434            
435      TrkTrack     *GetStoredTrack(int i);      TrkTrack     *GetStoredTrack(int i);
436      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
437    
438      TRefArray *GetTracks_NFitSorted();      TRefArray *GetTracks_NFitSorted();
439      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();};      TRefArray *GetTracks(){return this->GetTracks_NFitSorted();}
440            
441      Int_t     GetNTracks();      Int_t     GetNTracks();
442      TrkTrack* GetTrack(int i);      TrkTrack* GetTrack(int i);
443      TrkTrack* GetTrackImage(int i);      TrkTrack* GetTrackImage(int i);
444            
445      TrkLevel2*    GetTrkLevel2(){return this;}      TrkLevel2*     GetTrkLevel2(){return this;}
446      TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array      TClonesArray*  GetTrackArray(){return Track;}///< returns pointer to the track array
447            TClonesArray** GetPointerToTrackArray(){return &Track;}///< returns pointer to pointer to the track array
448      void   StatusDump(int view);      void   StatusDump(int view);
449      Bool_t StatusCheck(int view, int flagmask);      Bool_t StatusCheck(int view, int flagmask);
450    
Legend:
Removed from v.1.41  
changed lines
  Added in v.1.46
  ViewVC Help
Powered by ViewVC 1.1.23