--- DarthVader/TrackerLevel2/inc/TrkLevel2.h 2007/03/29 17:06:39 1.23 +++ DarthVader/TrackerLevel2/inc/TrkLevel2.h 2007/08/22 07:03:44 1.30 @@ -11,22 +11,25 @@ #include #include -//#include #include #include // z-coordinate of track state-vector reference-plane #define ZINI 23.5 -// upper and lower (mechanical) z-coordinate of the tracker -//#define ZTRKUP 22.29 -//#define ZTRKDW -22.22 // (mechanical) z-coordinate of the tracker planes -#define ZTRK6 -22.23 -#define ZTRK5 -13.32 -#define ZTRK4 -4.42 -#define ZTRK3 4.48 -#define ZTRK2 13.38 -#define ZTRK1 22.28 +#define ZTRK6 -22.22 +#define ZTRK5 -13.31 +#define ZTRK4 -4.41 +#define ZTRK3 4.49 +#define ZTRK2 13.39 +#define ZTRK1 22.29 +// magnet cavity dimensions +#define ZMAGNHIGH 21.83 +#define ZMAGNLOW -21.83 +#define XMAGNHIGH 8.07 +#define XMAGNLOW -8.07 +#define YMAGNHIGH 6.57 +#define YMAGNLOW -6.57 // (mechanical) x/y-coordinates of magnet cavity #define XTRKL -8.1 #define XTRKR 8.1 @@ -73,6 +76,14 @@ * A track is defined by the measured coordinates associated to it, the * track status vector, plus other quantities. * A track may have an "image", due to the ambiguity in the y view. + * + * Cluster flags: xgood[6], ygood[6] + * + * xgood/ygood = +/- 0lsccccccc + * | |||------- ID (1-7483647) of the included cluster + * | ||-------- sensor number (1,2 - increasing y) + * | |--------- ladder number (1,2,3 - increasing x) + * |------------- does-not/does include bad strips */ // ================================================================== class TrkTrack : public TObject { @@ -81,23 +92,20 @@ int seqno; /// 0;}; - Bool_t YGood(int ip){ return ygood[ip]> 0;}; - + Float_t GetDEDX(int ip); + Float_t GetDEDX(int ip,int iv); + Int_t GetLeverArmX(); + Int_t GetLeverArmY(); + Float_t GetChi2X(); + Float_t GetChi2Y(); + Float_t GetLnLX(); + Float_t GetLnLY(); + Float_t GetEffectiveAngle(int ip, int iv); + void SetMeasure(double *xmeas, double *ymeas, double *zmeas); void SetResolution(double *rx, double *ry); + void SetTail(double *tx, double *ty, double factor); + void SetStudentParam(int flag); void SetGood(int *xg, int *yg); void LoadField(TString s); - void Fit(double pfixed, int& fail, int iprint); + void Fit(double pfixed, int& fail, int iprint, int froml1); + void Fit(double pfixed, int& fail, int iprint){ Fit(pfixed,fail,iprint,0); }; void FitReset(); void SetTrackingMode(int trackmode); void SetPrecisionFactor(double fact); void SetStepMin(int istepmin); + Bool_t IsInsideCavity(); + + Bool_t EvaluateClusterPositions(); void FillMiniStruct(cMini2track&); void SetFromMiniStruct(cMini2track*); - TrkCluster *GetClusterX(int ip); - TrkCluster *GetClusterY(int ip); + Int_t GetClusterX_ID(int ip); + Int_t GetClusterY_ID(int ip); + Int_t GetLadder(int ip); + Int_t GetSensor(int ip); + Bool_t XGood(int ip){ return GetClusterX_ID(ip)!=-1; }; + Bool_t YGood(int ip){ return GetClusterY_ID(ip)!=-1; }; + void ResetXGood(int ip){ xgood[ip]=0; }; + void ResetYGood(int ip){ ygood[ip]=0; }; + void SetXGood(int ip, int clid, int is); + void SetYGood(int ip, int clid, int is); + + Bool_t BadClusterX(int ip){ return IsBad(ip,0); }; + Bool_t BadClusterY(int ip){ return IsBad(ip,1); }; + + Bool_t SaturatedClusterX(int ip){ return IsSaturated(ip,0); }; + Bool_t SaturatedClusterY(int ip){ return IsSaturated(ip,1); }; - Int_t GetClusterX_ID(int ip){return xgood[ip];}; - Int_t GetClusterY_ID(int ip){return ygood[ip];}; - TrkTrack* GetTrkTrack(){return this;}; friend class TrkLevel2; - ClassDef(TrkTrack,2); + ClassDef(TrkTrack,3); }; /** @@ -179,11 +210,9 @@ public: - TRef cls; - int plane; /// maximum number (level1 processing) + * |||| |||| |||| |||| ||________ 6 + * |||| |||| |||| |||| |_________ 7 + * |||| |||| |||| ||||___________ 8 n.clusters > maximum value (level2 processing) + * |||| |||| |||| |||____________ 9 n.couples per plane > maximum values (vector dimention) + * |||| |||| |||| ||_____________ 10 n.doublets > maximum values + * |||| |||| |||| |______________ 11 n.triplets > maximum values + * |||| |||| ||||________________ 12 n.yz-clouds > maximum values + * |||| |||| |||_________________ 13 n.xz-clouds > maximum values + * |||| |||| ||__________________ 14 n.candidate-tracks > maximum values + * |||| |||| |___________________ 15 n.couples per plane > maximum values (for Hough transform) + * |||| ||||_____________________ 16 + * + * + * For all data processed before June 2007 the event status was coded according to + * a different rule: + * + * Status of level1 processing + * 0 -- OK + * 1 -- missing packet + * 2 -- 1 CRC error + * 3 -- 2 on-line software alarm (latch-up flags asserted or n.transmitted-words = 0) + * 4 -- 3 jump in the trigger counter + * 10 -- 4 decode error + * 11 -- 5 n.clusters > maximum number (for level1 processing) + * Status of level2 processing + * 21 -- 0 n.clusters > maximum value (for level2 processing) + * 22 -- 1 n.couples per plane > maximum values (vector dimention) + * 23 -- 2 n.doublets > maximum values + * 24 -- 3 n.triplets > maximum values + * 25 -- 4 n.yz-clouds > maximum values + * 26 -- 5 n.xz-clouds > maximum values + * 27 -- 6 n.candidate-tracks > maximum values + * 28 -- 7 n.couples per plane > maximum values (for Hough transform) + * + * */ class TrkLevel2 : public TObject { @@ -220,6 +295,8 @@ public: Int_t good[12]; ///< event status + UInt_t VKmask[12]; ///< Viking-chip mask + UInt_t VKflag[12]; ///< Viking-chip flag TClonesArray *Track; ///< fitted tracks TClonesArray *SingletX; ///< x singlets @@ -245,27 +322,28 @@ void SetFromLevel2Struct() { SetFromLevel2Struct(&level2event_); }; void GetLevel2Struct(cTrkLevel2 *) const; void LoadField(TString); - float GetBX(float*); - float GetBY(float*); - float GetBZ(float*); + float GetBX(float* v){return TrkParams::GetBX(v);};///< Bx (kGauss) + 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 GetXTrkRight(){return XTRKR;}; Float_t GetYTrkLeft(){return YTRKL;}; Float_t GetYTrkRight(){return YTRKR;}; + Bool_t IsMaskedVK(int,int); + Bool_t GetVKMask(int,int); + Bool_t GetVKFlag(int,int); + TrkSinglet *GetSingletX(int); TrkSinglet *GetSingletY(int); TrkTrack *GetStoredTrack(int i); Int_t GetSeqNo(Int_t i) {return (((TrkTrack *)Track->At(i))->seqno);}; ///< Returns track sequential number -// TClonesArray *GetTracks_Chi2Sorted(); -// TClonesArray *GetTracks_NFitSorted(); -// TClonesArray *GetTracks(); + TRefArray *GetTracks_NFitSorted(); TRefArray *GetTracks(){return this->GetTracks_NFitSorted();}; -// int GetNTracks(){return this->GetTracks()->GetEntries();} Int_t GetNTracks(); TrkTrack* GetTrack(int i); TrkTrack* GetTrackImage(int i); @@ -273,7 +351,10 @@ TrkLevel2* GetTrkLevel2(){return this;} TClonesArray* GetTrackArray(){return Track;};///< returns pointer to the track array - ClassDef(TrkLevel2,2); + void StatusDump(int view); + Bool_t StatusCheck(int view, int flagmask); + + ClassDef(TrkLevel2,3); };