TrackerLevel2/inc/TrkLevel1.h

/**
 * \file TrkLevel1.h
 * \author Elena Vannuccini
 */
#ifndef trklevel1_h
#define trklevel1_h

#include <TObject.h>
#include <TClonesArray.h>
#include <TRef.h>

#include <TrkStruct.h>

#define INC_CUT 4.
/**
 * \brief Class to describe tracker clusters. 
 *
 * A cluster is defined as a group of adjacent strips whose signals exceed noise by INC_CUT sigma, and at least one strip has signal/noise > SEED_CUT. Tipical values are SEED_CUT = 7 and INC_CUT = 4.
 * 
 */
// ==================================================================
class TrkCluster : public TObject {

private:


public:
        
    int view; ///< view
    int maxs; ///< strip number (1-3072) of cluster center
    int indmax; 
    
    Int_t         CLlength; ///< number of stored strip info (signal+sigma+adc+bad)
    Float_t      *clsignal; //[CLlength] 
    Float_t      *clsigma;  //[CLlength]
    Int_t        *cladc;    //[CLlength]
    Bool_t       *clbad;    //[CLlength]
    
    TrkCluster();
    ~TrkCluster(){Clear();};
    TrkCluster(const TrkCluster&);

    void Dump();
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete(){Clear();};
        
    Float_t GetSignal(Int_t, Float_t);    ///< cluster signal
    Float_t GetSignal(Float_t cut)  { return GetSignal( (Int_t)0, cut); };
    Float_t GetSignal(Int_t nstrip) { return GetSignal( nstrip, (Float_t)(-1000) ); };
    Float_t GetSignal()             { return GetSignal( (Int_t)0,(Float_t)INC_CUT); };
    
    Float_t GetSignalToNoise(Int_t,Float_t);  ///< cluster signal/noise
    Float_t GetSignalToNoise(Float_t cut)  { return GetSignalToNoise( (Int_t)0, cut); };
    Float_t GetSignalToNoise(Int_t nstrip) { return GetSignalToNoise( nstrip, (Float_t)(-1000) ); };
    Float_t GetSignalToNoise()             { return GetSignalToNoise( (Int_t)0, (Float_t)INC_CUT); };


    Int_t   GetMultiplicity(Float_t);         ///< cluster multiplicity 
    Int_t   GetMultiplicity()  { return GetMultiplicity(INC_CUT);  };
    
    cTrkLevel1* GetLevel1Struct();
    
    Float_t GetCOG(Int_t);
    Float_t GetCOG(Float_t angle);
    Float_t GetCOG(){ return GetCOG(0); };
    Float_t GetETA(Int_t,float);
    Float_t GetETA(float angle){ return GetETA(0,angle); };
    
    Bool_t IsBad(Int_t);                ///< bad-cluster flag
    Bool_t IsSaturated(Int_t);          ///< saturated-cluster flag

    Int_t GetPlane() { return (Int_t)((view+1)/2);};         ///< plane number (1-6)
    Int_t GetLadder(){ return 1+(Int_t)((maxs-1)/1024);  };  ///< ladder number (1-3)
    Bool_t IsY(){ return (Bool_t)view%2; };
    Bool_t IsX(){ return !((Bool_t)view%2); };
        
    TrkCluster* GetTrkCluster(){ return this; };

    friend class TrkLevel1;

    ClassDef(TrkCluster,2);

};

/**
 * \brief Class to describe tracker LEVEL1 data.
 *
 * A Level1 tracker event is defined as a collection of clusters ( TrkCluster objects ).
 * The result of common-noise computation for each viking is also stored ( cnev[24][12] ).
 * A general flag summarize the event status (missing sections, crc failures, decoding errors ans so on...).
 */
class TrkLevel1 : public TObject {

 private:
    
 public:

    Int_t         good[12];       ///< event status
    Float_t       cn[24][12];     ///< CN
//    Float_t       cnrms[24][12];  ///< CN RMS
    Int_t         cnn[24][12];    ///< number of strips for CN computation

        
    TClonesArray *Cluster;        ///< clusters
        
    TrkLevel1();
    ~TrkLevel1(){Delete();};
    
    int nclstr() {return Cluster->GetEntries();}    ///< number of stored clusters
    
    void Dump();
    void SetFromLevel1Struct(cTrkLevel1 *, Bool_t);
//      void GetLevel1Struct(cTrkLevel1 *) const;
    cTrkLevel1* GetLevel1Struct();
    void Clear();
    void Delete();
    
    TrkCluster*   GetCluster(int);
    
    TrkLevel1*    GetTrkLevel1(){return this;}
    TClonesArray* GetClusters(){return Cluster;};   ///< returns pointer to the cluster array
    
    int           LoadPfaParam(TString);
    
    ClassDef(TrkLevel1,2);

};


#endif
1	/**
2	* \file TrkLevel1.h
3	* \author Elena Vannuccini
4	*/
5	#ifndef trklevel1_h
6	#define trklevel1_h
7
8	#include <TObject.h>
9	#include <TClonesArray.h>
10	#include <TRef.h>
11
12	#include <TrkStruct.h>
13
14	#define INC_CUT 4.
15	/**
16	* \brief Class to describe tracker clusters.
17	*
18	* A cluster is defined as a group of adjacent strips whose signals exceed noise by INC_CUT sigma, and at least one strip has signal/noise > SEED_CUT. Tipical values are SEED_CUT = 7 and INC_CUT = 4.
19	*
20	*/
21	// ==================================================================
22	class TrkCluster : public TObject {
23
24	private:
25
26
27	public:
28
29	int view; ///< view
30	int maxs; ///< strip number (1-3072) of cluster center
31	int indmax;
32
33	Int_t CLlength; ///< number of stored strip info (signal+sigma+adc+bad)
34	Float_t *clsignal; //[CLlength]
35	Float_t *clsigma; //[CLlength]
36	Int_t *cladc; //[CLlength]
37	Bool_t *clbad; //[CLlength]
38
39	TrkCluster();
40	~TrkCluster(){Clear();};
41	TrkCluster(const TrkCluster&);
42
43	void Dump();
44	void Clear();
45	void Clear(Option_t *option){Clear();};
46	void Delete(){Clear();};
47
48	Float_t GetSignal(Int_t, Float_t); ///< cluster signal
49	Float_t GetSignal(Float_t cut) { return GetSignal( (Int_t)0, cut); };
50	Float_t GetSignal(Int_t nstrip) { return GetSignal( nstrip, (Float_t)(-1000) ); };
51	Float_t GetSignal() { return GetSignal( (Int_t)0,(Float_t)INC_CUT); };
52
53	Float_t GetSignalToNoise(Int_t,Float_t); ///< cluster signal/noise
54	Float_t GetSignalToNoise(Float_t cut) { return GetSignalToNoise( (Int_t)0, cut); };
55	Float_t GetSignalToNoise(Int_t nstrip) { return GetSignalToNoise( nstrip, (Float_t)(-1000) ); };
56	Float_t GetSignalToNoise() { return GetSignalToNoise( (Int_t)0, (Float_t)INC_CUT); };
57
58
59	Int_t GetMultiplicity(Float_t); ///< cluster multiplicity
60	Int_t GetMultiplicity() { return GetMultiplicity(INC_CUT); };
61
62	cTrkLevel1* GetLevel1Struct();
63
64	Float_t GetCOG(Int_t);
65	Float_t GetCOG(Float_t angle);
66	Float_t GetCOG(){ return GetCOG(0); };
67	Float_t GetETA(Int_t,float);
68	Float_t GetETA(float angle){ return GetETA(0,angle); };
69
70	Bool_t IsBad(Int_t); ///< bad-cluster flag
71	Bool_t IsSaturated(Int_t); ///< saturated-cluster flag
72
73	Int_t GetPlane() { return (Int_t)((view+1)/2);}; ///< plane number (1-6)
74	Int_t GetLadder(){ return 1+(Int_t)((maxs-1)/1024); }; ///< ladder number (1-3)
75	Bool_t IsY(){ return (Bool_t)view%2; };
76	Bool_t IsX(){ return !((Bool_t)view%2); };
77
78	TrkCluster* GetTrkCluster(){ return this; };
79
80	friend class TrkLevel1;
81
82	ClassDef(TrkCluster,2);
83
84	};
85
86	/**
87	* \brief Class to describe tracker LEVEL1 data.
88	*
89	* A Level1 tracker event is defined as a collection of clusters ( TrkCluster objects ).
90	* The result of common-noise computation for each viking is also stored ( cnev[24][12] ).
91	* A general flag summarize the event status (missing sections, crc failures, decoding errors ans so on...).
92	*/
93	class TrkLevel1 : public TObject {
94
95	private:
96
97	public:
98
99	Int_t good[12]; ///< event status
100	Float_t cn[24][12]; ///< CN
101	// Float_t cnrms[24][12]; ///< CN RMS
102	Int_t cnn[24][12]; ///< number of strips for CN computation
103
104
105	TClonesArray *Cluster; ///< clusters
106
107	TrkLevel1();
108	~TrkLevel1(){Delete();};
109
110	int nclstr() {return Cluster->GetEntries();} ///< number of stored clusters
111
112	void Dump();
113	void SetFromLevel1Struct(cTrkLevel1 *, Bool_t);
114	// void GetLevel1Struct(cTrkLevel1 *) const;
115	cTrkLevel1* GetLevel1Struct();
116	void Clear();
117	void Delete();
118
119	TrkCluster* GetCluster(int);
120
121	TrkLevel1* GetTrkLevel1(){return this;}
122	TClonesArray* GetClusters(){return Cluster;}; ///< returns pointer to the cluster array
123
124	int LoadPfaParam(TString);
125
126	ClassDef(TrkLevel1,2);
127
128	};
129
130
131	#endif