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

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

Parent Directory | Revision Log | View Patch Patch

-revision 1.30 by pam-fi,
Wed Aug 22 07:03:44 2007 UTC
+revision 1.35 by pam-fi,
Wed Mar  5 17:00:18 2008 UTC
 Line 30
  #define XMAGNLOW -8.07
  #define YMAGNHIGH 6.57
  #define YMAGNLOW -6.57
+ // tof planes
+ #define ZS11  53.74
+ #define ZS12  53.04
+ #define ZS21  23.94
+ #define ZS22  23.44
+ #define ZS31 -23.49
+ #define ZS32 -24.34
  // (mechanical) x/y-coordinates of magnet cavity
- #define XTRKL -8.1
+ /* #define XTRKL -8.1 */
- #define XTRKR  8.1
+ /* #define XTRKR  8.1 */
- #define YTRKL -6.6
+ /* #define YTRKL -6.6 */
- #define YTRKR  6.6
+ /* #define YTRKR  6.6 */
  /**
   * \brief Class to describe, by points, a particle trajectory in the apparatus.
-Line 49 
 class Trajectory : public TObject{
+Line 57 
 class Trajectory : public TObject{
   public:
      int npoint; ///< number of evaluated points along the trajectory
-     float* x;   ///< x coordinates
+     float* x;   //[npoint]
-     float* y;   ///< y coordinates
+     float* y;   //[npoint]
-     float* z;   ///< z coordinates
+     float* z;   //[npoint]
-     float* thx; ///< x projected angle
+     float* thx; //[npoint]
-     float* thy; ///< y projected angle
+     float* thy; //[npoint]
-     float* tl;  ///< track length
+     float* tl;  //[npoint]
      Trajectory();
      Trajectory(int n);
-Line 67 
 class Trajectory : public TObject{
+Line 75 
 class Trajectory : public TObject{
      float GetLength(){float l=0; for(int i=0; i<npoint;i++)l=l+tl[i]; return l;};
      float GetLength(int,int);
-     ClassDef(Trajectory,2);
+     ClassDef(Trajectory,3);
  };
  /**
-Line 80 
 class Trajectory : public TObject{
+Line 88 
 class Trajectory : public TObject{
   * Cluster flags: xgood[6], ygood[6]
   *
   * xgood/ygood = +/- 0lsccccccc
-  *                |   |||------- ID (1-7483647) of the included cluster
+  * ccccccc ID (1-7483647) of the included cluster
-  *                |   ||-------- sensor number (1,2   - increasing y)
+  * s       sensor number (1,2   - increasing y)
-  *                |   |--------- ladder number (1,2,3 - increasing x)
+  * l       ladder number (1,2,3 - increasing x)
-  *                |------------- does-not/does include bad strips
+  * +/-     does-not/does include bad strips
+  *
   */
  // ==================================================================
  class TrkTrack : public TObject {
  private:
+ public:
      int   seqno;           ///<stored track sequential number
      int   image;           ///<sequential number of track-image
- public:
      float al[5];           ///<TRACK STATE VECTOR
      float coval[5][5];     ///<covariance matrix
-     int   xgood[6];        ///<cluster flag 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)
-     int   ygood[6];        ///<cluster flag 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)
      float xm[6];           ///<measured x coordinates
      float ym[6];           ///<measured y coordinates
      float zm[6];           ///<measured z coordinates
-Line 107 
 public:
+Line 116 
 public:
      float tailx[6];        ///<spatial resolution tail on X view
      float taily[6];        ///<spatial resolution tail on y view
      float chi2;            ///<chi2
-     int   nstep;           ///<n. step
+     int   nstep;           ///<n.step
      float xv[6];           ///<calculated x coordinates
      float yv[6];           ///<calculated y coordinates
      float zv[6];           ///<calculated z coordinates
-Line 115 
 public:
+Line 124 
 public:
      float ayv[6];          ///<calculated angles (deg) on y view
      float dedx_x[6];       ///<dE/dx in MIP (<0 if saturated)
      float dedx_y[6];       ///<dE/dx in MIP (<0 if saturated)
+     int   multmaxx[6];     ///<cluster multiplicity and strip of maximum on x view
+     int   multmaxy[6];     ///<cluster multiplicity and strip of maximum on y view
+     float seedx[6];        ///< seed of the cluster x
+     float seedy[6];        ///< seed of the cluster y
+     float xpu[6];          ///< x coordinate in pitch units
+     float ypu[6];          ///< y coordinate in pitch units
+     float xGF[14];         ///<calculated x coordinates on GF reference planes
+     float yGF[14];         ///<calculated y coordinates on GF reference planes
      TrkTrack();
      TrkTrack(const TrkTrack&);
-Line 131 
 public:
+Line 149 
 public:
      Int_t  GetSeqNo(){return seqno;}        ///< Returns the track sequential number
      Int_t  GetImageSeqNo(){return image;}   ///< Returns the track image sequential number
      Bool_t HasImage(){return !(image==-1);} ///< Returns true if the track has an image
-     int DoTrack(Trajectory* t);                         ///< Evaluates the trajectory in the apparatus.
+     int DoTrack(Trajectory* t);             ///< Evaluates the trajectory in the apparatus.
-     int DoTrack2(Trajectory* t);                        ///< Evaluates the trajectory in the apparatus.
+     int DoTrack2(Trajectory* t);            ///< Evaluates the trajectory in the apparatus.
-     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.
      Int_t GetNX(){Int_t n=0; for(Int_t i=0; i<6; i++)n+=(Int_t)XGood(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;};
      Int_t GetNtot(){return GetNX()+GetNY();};
-Line 167 
 public:
+Line 185 
 public:
      void SetTrackingMode(int trackmode);
      void SetPrecisionFactor(double fact);
      void SetStepMin(int istepmin);
-     Bool_t IsInsideCavity();
+     void SetDeltaB(int id, double db);
+     Bool_t IsInsideCavity(float);
+     Bool_t IsInsideCavity(){ return IsInsideCavity(0.); };
+     Bool_t IsInsideAcceptance();
      Bool_t EvaluateClusterPositions();
-Line 191 
 public:
+Line 213 
 public:
      Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); };
      Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); };
+     Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); };
+     Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); };
+     Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); };
+     Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); };
+     Float_t GetClusterX_Seed(int ip){ return seedx[ip]; };
+     Float_t GetClusterY_Seed(int ip){ return seedy[ip]; };
+ /*     Float_t GetClusterX_CoordinatePU(int ip); */
+ /*     Float_t GetClusterY_CoordinatePU(int ip); */
+     Float_t GetYav();
+     Float_t GetXav();
+     Float_t GetZav();
+     Int_t GetNColumns();
+     Float_t GetDEDX_max(int ip, int iv);
+     Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); };
+     Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); };
+     Float_t GetDEDX_min(int ip, int iv);
+     Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); };
+     Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); };
+     Float_t GetResidual_max(int ip, int iv);
+     Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); };
+     Float_t GetResidual_max(){ return GetResidual_max(-1,-1); };
+     Int_t GetClusterX_Multiplicity_max();
+     Int_t GetClusterX_Multiplicity_min();
+     Int_t GetClusterY_Multiplicity_max();
+     Int_t GetClusterY_Multiplicity_min();
+     Float_t GetClusterX_Seed_min();
+     Float_t GetClusterY_Seed_min();
      TrkTrack* GetTrkTrack(){return this;};
      friend class TrkLevel2;
-     ClassDef(TrkTrack,3);
+     ClassDef(TrkTrack,5);
  };
  /**
-Line 213 
 public:
+Line 269 
 public:
      int plane;       ///<plane
      float coord[2];  ///<coordinate (on sensor 1 and 2)
      float sgnl;      ///<cluster signal in MIP (<0 if saturated)
+     int multmax;     ///<cluster multiplicity and strip of maximum
      TrkSinglet();
      TrkSinglet(const TrkSinglet&);
-Line 224 
 public:
+Line 281 
 public:
      void Delete(){Clear();};
      Float_t GetSignal(){return fabs(sgnl);}
      Bool_t IsSaturated(){return (sgnl<0); };
+     Bool_t IsBad()                 { return multmax<=0; };
+     Int_t GetCluster_Multiplicity(){ return (Int_t)(abs(multmax)/10000); };
+     Int_t GetCluster_MaxStrip()    { return (Int_t)(abs(multmax)%10000); };
      friend class TrkLevel2;
-     ClassDef(TrkSinglet,3);
+     ClassDef(TrkSinglet,4);
  };
-Line 245 
 public:
+Line 307 
 public:
   * The event status indicates the processing status of data from each DSP, according to the following
   * notation:
   *
-  *     xxxx xxxx xxxx xxxx xxxx xxxx
+  * LSB --> 0 missing packet
-  *     |||| |||| |||| |||| |||| ||||_ 0 missing packet
+  *         1 CRC error
-  *     |||| |||| |||| |||| |||| |||__ 1 CRC error
+  *         2 on-line software alarm (latch-up, timeout ecc...)
-  *     |||| |||| |||| |||| |||| ||___ 2 on-line software alarm (latch-up, timeout ecc...)
+  *         3 jump in the trigger counter
-  *     |||| |||| |||| |||| |||| |____ 3 jump in the trigger counter
+  *         4 decode error
-  *     |||| |||| |||| |||| ||||______ 4 decode error
+  *         5 n.clusters > maximum number (level1 processing)
-  *     |||| |||| |||| |||| |||_______ 5 n.clusters > maximum number (level1 processing)
+  *         6
-  *     |||| |||| |||| |||| ||________ 6
+  *         7
-  *     |||| |||| |||| |||| |_________ 7
+  *         8 n.clusters > maximum value (level2 processing)
-  *     |||| |||| |||| ||||___________ 8 n.clusters > maximum value (level2 processing)
+  *         9 n.couples per plane > maximum values (vector dimention)
-  *     |||| |||| |||| |||____________ 9 n.couples per plane > maximum values (vector dimention)
+  *         10 n.doublets > maximum values
-  *     |||| |||| |||| ||_____________ 10 n.doublets > maximum values
+  *         11 n.triplets > maximum values
-  *     |||| |||| |||| |______________ 11 n.triplets > maximum values
+  *         12 n.yz-clouds > maximum values
-  *     |||| |||| ||||________________ 12 n.yz-clouds > maximum values
+  *         13 n.xz-clouds > maximum values
-  *     |||| |||| |||_________________ 13 n.xz-clouds > maximum values
+  *         14 n.candidate-tracks > maximum values
-  *     |||| |||| ||__________________ 14 n.candidate-tracks > maximum values
+  *         15 n.couples per plane > maximum values (for Hough transform)
-  *     |||| |||| |___________________ 15 n.couples per plane > maximum values (for Hough transform)
+  * MSB --> 16
-  *     |||| ||||_____________________ 16
   *
   *
   * For all data processed before June 2007 the event status was coded according to
-Line 326 
 class TrkLevel2 : public TObject {
+Line 387 
 class TrkLevel2 : public TObject {
      float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss)
      float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss)
      Float_t GetZTrk(Int_t);
-     Float_t GetXTrkLeft(){return XTRKL;};
+     Float_t GetXTrkLeft(){return XMAGNLOW;};
-     Float_t GetXTrkRight(){return XTRKR;};
+     Float_t GetXTrkRight(){return XMAGNHIGH;};
-     Float_t GetYTrkLeft(){return YTRKL;};
+     Float_t GetYTrkLeft(){return YMAGNLOW;};
-     Float_t GetYTrkRight(){return YTRKR;};
+     Float_t GetYTrkRight(){return YMAGNHIGH;};
      Bool_t IsMaskedVK(int,int);
      Bool_t GetVKMask(int,int);

 Legend:



Removed from v.1.30
 


changed lines


 
Added in v.1.35
 Legend:



Removed from v.1.30
 


changed lines


 
Added in v.1.35
-Removed from v.1.30
+Added in v.1.35

	ViewVC Help
Powered by ViewVC 1.1.23