31 |
#define YMAGNHIGH 6.57 |
#define YMAGNHIGH 6.57 |
32 |
#define YMAGNLOW -6.57 |
#define YMAGNLOW -6.57 |
33 |
// (mechanical) x/y-coordinates of magnet cavity |
// (mechanical) x/y-coordinates of magnet cavity |
34 |
#define XTRKL -8.1 |
/* #define XTRKL -8.1 */ |
35 |
#define XTRKR 8.1 |
/* #define XTRKR 8.1 */ |
36 |
#define YTRKL -6.6 |
/* #define YTRKL -6.6 */ |
37 |
#define YTRKR 6.6 |
/* #define YTRKR 6.6 */ |
38 |
|
|
39 |
/** |
/** |
40 |
* \brief Class to describe, by points, a particle trajectory in the apparatus. |
* \brief Class to describe, by points, a particle trajectory in the apparatus. |
90 |
|
|
91 |
private: |
private: |
92 |
|
|
93 |
|
public: |
94 |
|
|
95 |
int seqno; ///<stored track sequential number |
int seqno; ///<stored track sequential number |
96 |
int image; ///<sequential number of track-image |
int image; ///<sequential number of track-image |
97 |
|
|
|
public: |
|
|
|
|
98 |
float al[5]; ///<TRACK STATE VECTOR |
float al[5]; ///<TRACK STATE VECTOR |
99 |
float coval[5][5]; ///<covariance matrix |
float coval[5][5]; ///<covariance matrix |
100 |
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) |
101 |
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) |
102 |
float xm[6]; ///<measured x coordinates |
float xm[6]; ///<measured x coordinates |
103 |
float ym[6]; ///<measured y coordinates |
float ym[6]; ///<measured y coordinates |
104 |
float zm[6]; ///<measured z coordinates |
float zm[6]; ///<measured z coordinates |
107 |
float tailx[6]; ///<spatial resolution tail on X view |
float tailx[6]; ///<spatial resolution tail on X view |
108 |
float taily[6]; ///<spatial resolution tail on y view |
float taily[6]; ///<spatial resolution tail on y view |
109 |
float chi2; ///<chi2 |
float chi2; ///<chi2 |
110 |
int nstep; ///<n. step |
int nstep; ///<n.step |
111 |
float xv[6]; ///<calculated x coordinates |
float xv[6]; ///<calculated x coordinates |
112 |
float yv[6]; ///<calculated y coordinates |
float yv[6]; ///<calculated y coordinates |
113 |
float zv[6]; ///<calculated z coordinates |
float zv[6]; ///<calculated z coordinates |
115 |
float ayv[6]; ///<calculated angles (deg) on y view |
float ayv[6]; ///<calculated angles (deg) on y view |
116 |
float dedx_x[6]; ///<dE/dx in MIP (<0 if saturated) |
float dedx_x[6]; ///<dE/dx in MIP (<0 if saturated) |
117 |
float dedx_y[6]; ///<dE/dx in MIP (<0 if saturated) |
float dedx_y[6]; ///<dE/dx in MIP (<0 if saturated) |
118 |
|
int multmaxx[6]; ///<cluster multiplicity and strip of maximum on x view |
119 |
|
int multmaxy[6]; ///<cluster multiplicity and strip of maximum on y view |
120 |
|
float seedx[6]; ///< seed of the cluster x |
121 |
|
float seedy[6]; ///< seed of the cluster y |
122 |
|
float xpu[6]; ///< x coordinate in pitch units |
123 |
|
float ypu[6]; ///< y coordinate in pitch units |
124 |
|
|
125 |
TrkTrack(); |
TrkTrack(); |
126 |
TrkTrack(const TrkTrack&); |
TrkTrack(const TrkTrack&); |
173 |
void SetTrackingMode(int trackmode); |
void SetTrackingMode(int trackmode); |
174 |
void SetPrecisionFactor(double fact); |
void SetPrecisionFactor(double fact); |
175 |
void SetStepMin(int istepmin); |
void SetStepMin(int istepmin); |
176 |
|
void SetDeltaB(int id, double db); |
177 |
|
|
178 |
Bool_t IsInsideCavity(); |
Bool_t IsInsideCavity(); |
179 |
|
|
180 |
Bool_t EvaluateClusterPositions(); |
Bool_t EvaluateClusterPositions(); |
199 |
Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }; |
Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }; |
200 |
Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }; |
Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }; |
201 |
|
|
202 |
|
Int_t GetClusterX_Multiplicity(int ip){ return (Int_t)(multmaxx[ip]/10000); }; |
203 |
|
Int_t GetClusterY_Multiplicity(int ip){ return (Int_t)(multmaxy[ip]/10000); }; |
204 |
|
Int_t GetClusterX_MaxStrip(int ip){ return (Int_t)(multmaxx[ip]%10000); }; |
205 |
|
Int_t GetClusterY_MaxStrip(int ip){ return (Int_t)(multmaxy[ip]%10000); }; |
206 |
|
Float_t GetClusterX_Seed(int ip){ return seedx[ip]; }; |
207 |
|
Float_t GetClusterY_Seed(int ip){ return seedy[ip]; }; |
208 |
|
/* Float_t GetClusterX_CoordinatePU(int ip); */ |
209 |
|
/* Float_t GetClusterY_CoordinatePU(int ip); */ |
210 |
|
|
211 |
|
Float_t GetYav(); |
212 |
|
Float_t GetXav(); |
213 |
|
Float_t GetZav(); |
214 |
|
|
215 |
|
Int_t GetNColumns(); |
216 |
|
|
217 |
|
Float_t GetDEDX_max(int ip, int iv); |
218 |
|
Float_t GetDEDX_max(int iv){ return GetDEDX_max(-1,iv); }; |
219 |
|
Float_t GetDEDX_max(){ return GetDEDX_max(-1,-1); }; |
220 |
|
Float_t GetDEDX_min(int ip, int iv); |
221 |
|
Float_t GetDEDX_min(int iv){ return GetDEDX_min(-1,iv); }; |
222 |
|
Float_t GetDEDX_min(){ return GetDEDX_min(-1,-1); }; |
223 |
|
|
224 |
|
Float_t GetResidual_max(int ip, int iv); |
225 |
|
Float_t GetResidual_max(int iv){ return GetResidual_max(-1,iv); }; |
226 |
|
Float_t GetResidual_max(){ return GetResidual_max(-1,-1); }; |
227 |
|
|
228 |
|
Int_t GetClusterX_Multiplicity_max(); |
229 |
|
Int_t GetClusterX_Multiplicity_min(); |
230 |
|
Int_t GetClusterY_Multiplicity_max(); |
231 |
|
Int_t GetClusterY_Multiplicity_min(); |
232 |
|
|
233 |
|
Float_t GetClusterX_Seed_min(); |
234 |
|
Float_t GetClusterY_Seed_min(); |
235 |
|
|
236 |
TrkTrack* GetTrkTrack(){return this;}; |
TrkTrack* GetTrkTrack(){return this;}; |
237 |
|
|
238 |
friend class TrkLevel2; |
friend class TrkLevel2; |
239 |
|
|
240 |
ClassDef(TrkTrack,3); |
ClassDef(TrkTrack,4); |
241 |
|
|
242 |
}; |
}; |
243 |
/** |
/** |
368 |
float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss) |
float GetBY(float* v){return TrkParams::GetBY(v);};///< By (kGauss) |
369 |
float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss) |
float GetBZ(float* v){return TrkParams::GetBZ(v);};///< Bz (kGauss) |
370 |
Float_t GetZTrk(Int_t); |
Float_t GetZTrk(Int_t); |
371 |
Float_t GetXTrkLeft(){return XTRKL;}; |
Float_t GetXTrkLeft(){return XMAGNLOW;}; |
372 |
Float_t GetXTrkRight(){return XTRKR;}; |
Float_t GetXTrkRight(){return XMAGNHIGH;}; |
373 |
Float_t GetYTrkLeft(){return YTRKL;}; |
Float_t GetYTrkLeft(){return YMAGNLOW;}; |
374 |
Float_t GetYTrkRight(){return YTRKR;}; |
Float_t GetYTrkRight(){return YMAGNHIGH;}; |
375 |
|
|
376 |
Bool_t IsMaskedVK(int,int); |
Bool_t IsMaskedVK(int,int); |
377 |
Bool_t GetVKMask(int,int); |
Bool_t GetVKMask(int,int); |