/[PAMELA software]/DarthVader/ToFLevel2/inc/ToFLevel2.h
ViewVC logotype

Diff of /DarthVader/ToFLevel2/inc/ToFLevel2.h

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by mocchiut, Fri Jun 30 09:22:02 2006 UTC revision 1.20 by mocchiut, Mon Mar 3 09:51:01 2008 UTC
# Line 1  Line 1 
1    /**
2     * \file ToFLevel2.h
3     * \author Gianfranca DeRosa / Wolfgang Menn
4     */
5    
6  #ifndef ToFLevel2_h  #ifndef ToFLevel2_h
7  #define ToFLevel2_h  #define ToFLevel2_h
8  //  //
9  #include <TObject.h>  #include <TObject.h>
10    #include <TArrayI.h>
11    #include <TArrayF.h>
12  #include <TClonesArray.h>  #include <TClonesArray.h>
13    
14    #include <math.h> // EMILIANO
15    
16    #include <ToFStruct.h>
17    
18    
19  //  //
20  // class which contains track related variables  // class which contains track related variables
21  //  //
# Line 13  Line 26 
26  #define ZTOF31 -23.49  #define ZTOF31 -23.49
27  #define ZTOF32 -24.34  #define ZTOF32 -24.34
28    
29    
30    /**
31     * \brief  Class which contains the PMT data
32     *
33     * If there is a valid ADC or a TDC value (value<4095) for a PMT, both ADC and TDC data
34     * are stored in the PMT class.
35     * Look in the ToFLevel2Ex.cxx example in the repository how to read the PMT class.
36     */
37    class ToFPMT : public TObject {
38    
39     private:
40    
41     public:
42        Int_t pmt_id;     ///<the identification number of the PMT from 0 to 47
43        Float_t adc;      ///<raw ADC value for this PMT
44        Float_t tdc;      ///<raw TDC value for this PMT
45        Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT
46        //
47        ToFPMT();
48        ToFPMT(const ToFPMT&);
49        //
50        ToFPMT* GetToFPMT(){return this;};
51        void Clear(Option_t *t="");
52    
53    
54    
55        ClassDef(ToFPMT,2);
56    };
57    
58    
59    /**
60     * \brief Class which contains the tracker related variables
61     *
62     * We can use the ToF standalone to find hitted paddles, calculate beta, etc..
63     * These results are then stored with the "trkseqno" = -1.
64     * If we use the track from the tracker, then the penetration points in the
65     * scintillators are calculated, which defines the hitted paddles. For these paddles
66     * we calculate then all the output.
67     * Note: The artificial ADC values are stored as dEdx in the output, the dEdx will be
68     * by definition = 1. However, the artificial TDC values are just used internally
69     * and not stored in the output. But one can see in both cases which PMT has artificial
70     * values using "adcflag" and "tdcflag".
71     * Look in the ToFLevel2Ex.cxx example in the repository how to read the tracker related
72     * variables.
73     */
74  class ToFTrkVar : public TObject {  class ToFTrkVar : public TObject {
75      
76   private:   private:
77    
78   public:   public:
   Int_t trkseqno; // tracker entry coming from tracker, 100 if image track is used, -100 if the track is not consistent with MyDetector2 one  
79    //    //
80    Float_t beta_a[13];    Int_t trkseqno; ///< tracker sequ. number: -1=ToF standalone, 0=first Tracker track, ...
81    Float_t adc_c[12][4];    //
82        Int_t npmttdc;  ///<number of the TDC measurements used to evaluate beta
83      TArrayI pmttdc;  ///<contains the ID (0..47) for the PMT used to evaluate beta
84      TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value
85    
86    /**
87     * \brief beta, 12 measurements for the 12  combinations, beta[13] is modified weighted mean
88     *
89     * The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for
90     * S2-S3 and S1-S2.
91     * The calculation of beta[13] is now modified:
92     * We check the individual weights for artificial TDC values, then calculate
93     * am mean beta for the first time. In a second step we loop again through
94     * the single measurements, checking for the residual from the mean
95     * The cut on the residual reject measurements > "x"-sigma. A chi2 value is
96     * calculated, furthermore a "quality" value by adding the weights which
97     * are finally used. If all measurements are taken, "quality" will be = 505.
98     * A chi2 cut around 3-4 and a quality-cut > 400 is needed for clean beta
99     * The Level2 beta[12] which is derived in the fortran routines uses: 10.,200.,20.
100     * This is not a very high quality measurement. One can re-calculate a new beta[13]
101     * using the L2-method "CalcBeta"
102     */
103      Float_t beta[13];
104      //
105      Int_t npmtadc;  ///<number of the ADC measurements used for dEdx evaluation
106      TArrayI pmtadc; ///<contains the ID (0..47) for the PMT used to evaluate dEdx
107      TArrayI adcflag; ///<flag for artificial ADCs, "0" if normal ADC value
108      TArrayF dedx;    ///<energy loss for this PMT in mip
109      //
110      Float_t xtofpos[3];  ///<x-measurement using the TDC values and the calibration from S12, S21, S32
111      Float_t ytofpos[3];  ///<x-measurement using the TDC values and the calibration from S11, S22, S31
112      //
113      Float_t xtr_tof[6];  ///<x-measurement in the ToF layers from tracker
114      Float_t ytr_tof[6];  ///<x-measurement in the ToF layers from tracker
115    //    //
116    ToFTrkVar();    ToFTrkVar();
117    ToFTrkVar(const ToFTrkVar&);    ToFTrkVar(const ToFTrkVar&);
118      
119    ToFTrkVar* GetToFTrkVar(){return this;};    ToFTrkVar* GetToFTrkVar(){return this;};
120        void Clear(Option_t *t="");
121    
122    ClassDef(ToFTrkVar,1);    ClassDef(ToFTrkVar,1);
123    //    //
124  };  };
125    
126    /**
127     * \brief  Class to describe ToF LEVEL2 data
128     *
129     */
130    
131  class ToFLevel2 : public TObject {  class ToFLevel2 : public TObject {
132   private:   private:
133      
134   public:   public:
   Float_t xtofpos[3];  
   Float_t ytofpos[3];  
   Int_t  tof_i_flag[6];  
   Int_t tof_j_flag[6];  
   Float_t betatof_a[13];  
   Float_t adctof_c[12][4];  
   Float_t tdc_c[12][4];  
   TClonesArray *ToFTrk; // track related variable class  
135    //    //
136      TClonesArray *PMT; ///<class needed to store PMT hit informations
137      TClonesArray *ToFTrk; ///<track related variable class
138      Int_t tof_j_flag[6];  ///<number of hitted paddle(s) for each ToF layer: flag = flag + 2**(paddlenumber-1)
139      //
140      Int_t unpackError;///< zero if no error presente
141      Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise
142      //
143      Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07
144    
145      Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  wm feb 08
146    
147      //
148    //  Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t chi2cut, Float_t qualitycut);   // wm feb 08
149    //    //
150    // methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables    // methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables
151    //    //
152    Int_t ntrk(){return ToFTrk->GetEntries();};    Int_t ntrk(){return ToFTrk->GetEntries();};
153        Int_t npmt(){return PMT->GetEntries();};
154      
155      //
156      void GetLevel2Struct(cToFLevel2 *) const;
157      //
158      ToFTrkVar *GetToFTrkVar(Int_t notrack);
159      ToFPMT *GetToFPMT(Int_t nohit);
160      Int_t GetPMTid(Int_t gg, Int_t hh);
161      TString GetPMTName(Int_t ind);
162      
163      Int_t GetPlaneIndex(Int_t pmt_id);
164      void GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]);
165      void GetPMTIndex(Int_t pmt_id, Int_t &gg, Int_t &hh);
166    
167      // gf Apr 07
168      void GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
169      TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);
170      Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07
171      void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07
172      void GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright); // gf Apr 07
173      void GetPaddleGeometry(Int_t plane, Int_t paddle, Float_t &xleft, Float_t &xright, Float_t &yleft, Float_t &yright); // gf Apr 07
174      Int_t GetPaddleid(Int_t plane, Int_t paddle);
175      void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle);
176      Int_t GetNPaddle(Int_t plane);
177    //    //
178    ToFTrkVar *GetToFTrkVar(Int_t notrack);      //
179      
180    //    //
181    // constructor    // constructor
182    //    //
183    ToFLevel2();    ToFLevel2();
184      ~ToFLevel2(){Delete();}; //ELENA
185      void Delete(Option_t *t=""); //ELENA
186      void Set();//ELENA
187    //    //
188    //    //
189    ToFLevel2*   GetToFLevel2(){return this;};    ToFLevel2*   GetToFLevel2(){return this;};
190      Float_t      GetZTOF(Int_t plane_id){  
191    /**
192     * Method to get the z-position of the 6 TOF layers from the plane ID
193     * @param plane_id Plane ID (11 12 21 22 31 32)
194     */
195      Float_t      GetZTOF(Int_t plane_id){
196        switch(plane_id){        switch(plane_id){
197        case 11: return ZTOF11;        case 11: return ZTOF11;
198        case 12: return ZTOF12;        case 12: return ZTOF12;
# Line 71  class ToFLevel2 : public TObject { Line 202  class ToFLevel2 : public TObject {
202        case 32: return ZTOF32;        case 32: return ZTOF32;
203        default: return 0.;        default: return 0.;
204        };        };
205      };    };
206        
207      //      //
208      // Paddles position      // Paddles position
209      //      //
# Line 84  class ToFLevel2 : public TObject { Line 215  class ToFLevel2 : public TObject {
215        S31 3 paddles  15.0 x 6.0 cm        S31 3 paddles  15.0 x 6.0 cm
216        S32 3 paddles  18.0 x 5.0 cm        S32 3 paddles  18.0 x 5.0 cm
217      */      */
218        
219      Int_t  GetToFPlaneID(Int_t ip);      Int_t  GetToFPlaneID(Int_t ip);
220      Int_t  GetToFPlaneIndex(Int_t plane_id);      Int_t  GetToFPlaneIndex(Int_t plane_id);
221      Bool_t HitPaddle(Int_t ,Int_t);      Bool_t HitPaddle(Int_t ,Int_t);
222      Int_t  GetNHitPaddles(Int_t plane);      Int_t  GetNHitPaddles(Int_t plane);
223      void Clear();      void Clear(Option_t *t="");
224      //      //
225      ClassDef(ToFLevel2,1);      ClassDef(ToFLevel2,3);
226  };  };
227    
228  #endif  #endif
229    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.23