/[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.20 by mocchiut, Mon Mar 3 09:51:01 2008 UTC revision 1.35 by pam-fi, Tue Mar 24 11:09:54 2015 UTC
# Line 1  Line 1 
1  /**  /**
2   * \file ToFLevel2.h   * \file ToFLevel2.h
3   * \author Gianfranca DeRosa / Wolfgang Menn   * \author Gianfranca DeRosa / Wolfgang Menn / Rita Carbone with E. M. supervision
4   */   */
5    
6  #ifndef ToFLevel2_h  #ifndef ToFLevel2_h
# Line 12  Line 12 
12  #include <TClonesArray.h>  #include <TClonesArray.h>
13    
14  #include <math.h> // EMILIANO  #include <math.h> // EMILIANO
15    #include <iostream> // from ToFLevel2.cpp
16    #include <fstream> // Emiliano
17    #include <sstream> // Emiliano
18    #include <string> // Emiliano
19    
20    
21  #include <ToFStruct.h>  #include <ToFStruct.h>
22    
23    #include <TrkLevel2.h> // Emiliano
24    #include <TrigLevel2.h> // Emiliano
25    #include <GLTables.h> // Emiliano
26    #include <OrbitalInfo.h> // Emiliano
27    #include <ToFCore.h> // Emiliano
28    #include <physics/tof/TofEvent.h>
29    
30    //
31    // Declaration of the core fortran routines
32    //
33    #define tofl2com tofl2com_
34    extern "C" int tofl2com();
35    #define toftrk toftrk_
36    extern "C" int toftrk();
37    #define rdtofcal rdtofcal_
38    //extern "C" int rdtofcal(char [], int *);
39    extern "C" int rdtofcal(const char *, int *);
40    
41  //  //
42  // class which contains track related variables  // class which contains track related variables
# Line 27  Line 49 
49  #define ZTOF32 -24.34  #define ZTOF32 -24.34
50    
51    
52    class ToFGeom : public TObject {
53    
54    private:
55      TArrayI  ePlane, eXY;
56    
57    public:
58      ToFGeom() {
59        int plane[24] = {
60          0, 0, 0, 0, 0, 0, 0, 0,
61          1, 1, 1, 1, 1, 1,
62          2, 2,            
63          3, 3,            
64          4, 4, 4,
65          5, 5, 5
66        };
67        int plXY[6]= {  2,  1,  1,  2,  2,  1 }; // X==1, Y==2 */
68        ePlane.Set(24,plane);
69        eXY.Set(6,plXY);
70      }
71    
72      int GetPad(   int idpmt)  { return (int)((idpmt+0.5)/2.); }
73      int GetPlane( int idpmt)  { return ePlane[ GetPad(idpmt) ];  }
74      int GetXY(    int idpmt)  { return eXY[ GetPlane(idpmt) ]; }
75    
76       ClassDef(ToFGeom,1);
77    
78    };
79    
80    
81  /**  /**
82   * \brief  Class which contains the PMT data   * \brief  Class which contains the PMT data
83   *   *
# Line 40  class ToFPMT : public TObject { Line 91  class ToFPMT : public TObject {
91    
92   public:   public:
93      Int_t pmt_id;     ///<the identification number of the PMT from 0 to 47      Int_t pmt_id;     ///<the identification number of the PMT from 0 to 47
94      Float_t adc;      ///<raw ADC value for this PMT      Float_t adc;      ///<raw ADC value for this PMT
95      Float_t tdc;      ///<raw TDC value for this PMT      Float_t tdc;      ///<raw TDC value for this PMT
96      Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT      Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT
97        Float_t l0flag_adc; ///< warning flags from unpacking
98        Float_t l0flag_tdc; ///< warning flags from unpacking
99      //      //
100      ToFPMT();      ToFPMT();
101      ToFPMT(const ToFPMT&);      ToFPMT(const ToFPMT&);
# Line 50  class ToFPMT : public TObject { Line 103  class ToFPMT : public TObject {
103      ToFPMT* GetToFPMT(){return this;};      ToFPMT* GetToFPMT(){return this;};
104      void Clear(Option_t *t="");      void Clear(Option_t *t="");
105    
106        ClassDef(ToFPMT,3);
107    };
108    
109    /**
110     * \brief  Class used to calibrate adc to dEdx for each PMT
111     *
112     * Class used to calibrate adc to dEdx for each PMT
113     */
114    class ToFdEdx : public TObject {
115    
116     private:
117        //
118        ToFGeom  eGeom;  // ToF geometry
119        //
120        Float_t  adc_he;
121        TArrayF *eDEDXpmt;  // 0-47  pmt dEdx
122        // parameters:
123        TArrayF PMTsat;  // 0-47  saturation parameters
124        Float_t adc[48];
125        //
126    
127        TArrayF parAtt[48];      // 48 x 6
128        TArrayF parPos[48];      // 48 x 4
129        TArrayF parDesatBB[48];  // 48 x 3
130        TArrayF parBBneg[48];    // 48 x 3
131        TArrayF  parBBpos;    // 48 x 1
132    
133        double f_adcPC( float x );
134        double f_BB( TArrayF &p, float x );
135        double f_BB5B( float x );
136        double f_att( TArrayF &p, float x ) ;
137        double f_att5B( float x );
138        double f_desatBB( TArrayF &p, float x );  
139        double f_desatBB5B( float x );
140        double f_pos( TArrayF &p, float x );
141        double f_pos5B( float x );
142        float Get_adc_he( int id, float pl_x[6], float pl_y[6]);
143    
144        Bool_t conn[12];
145    
146        UInt_t ts[12];
147        UInt_t te[12];    
148      
149        
150     public:
151        ToFdEdx();   // class constructor
152        ~ToFdEdx();  // class distructor
153        //
154        void Clear(Option_t *option="");
155        void Delete(Option_t *option="");
156    
157        void Init(pamela::tof::TofEvent *tofl0 );  // init parameters
158        void Init(Int_t i, Int_t j, Float_t adce);
159        void Define_PMTsat();
160    
161        void ReadParAtt(            const char *fname  );
162        void ReadParPos(            const char *fname  );
163        void ReadParBBneg(          const char *fname  );
164        void ReadParBBpos(          const char *fname  );
165        void ReadParDesatBB(        const char *fname  );
166    
167        void CheckConnectors(UInt_t atime, GL_PARAM *glparam, TSQLServer *dbc);
168    
169        void Process( UInt_t atime, Float_t betamean, Float_t *xtr_tof, Float_t *ytr_tof, Int_t exitat=-1); //
170        void Print(Option_t *option="");
171    
172      ClassDef(ToFPMT,2);      //    Float_t GetdEdx_pmt(Int_t ipmt) { return (Float_t)eDEDXpmt[ipmt]; }  // 0-47 dEdx for each PMT for tracked events
173        Float_t GetdEdx_pmt(Int_t ipmt) { return eDEDXpmt->At(ipmt); }  // 0-47 dEdx for each PMT for tracked events
174        //
175        ToFdEdx* GetToFdEdx(){return this;};
176        ClassDef(ToFdEdx,3);
177  };  };
178    
179    
# Line 119  class ToFTrkVar : public TObject { Line 240  class ToFTrkVar : public TObject {
240    ToFTrkVar* GetToFTrkVar(){return this;};    ToFTrkVar* GetToFTrkVar(){return this;};
241    void Clear(Option_t *t="");    void Clear(Option_t *t="");
242    
243      Float_t CalcBeta(Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  Elena 2015
244    
245    ClassDef(ToFTrkVar,1);    ClassDef(ToFTrkVar,1);
246    //    //
247  };  };
# Line 137  class ToFLevel2 : public TObject { Line 260  class ToFLevel2 : public TObject {
260    TClonesArray *ToFTrk; ///<track related variable class    TClonesArray *ToFTrk; ///<track related variable class
261    Int_t tof_j_flag[6];  ///<number of hitted paddle(s) for each ToF layer: flag = flag + 2**(paddlenumber-1)    Int_t tof_j_flag[6];  ///<number of hitted paddle(s) for each ToF layer: flag = flag + 2**(paddlenumber-1)
262    //    //
263    Int_t unpackError;///< zero if no error presente    Int_t unpackError;///< zero if no errors
264      Int_t unpackWarning;///< zero if no warnings | bitwise: ... 3 2 1 0 --> bit 0 = discharging flag on at least one PMT (adc OR tdc), bit 1 = there is at least one PMT off (no tdc nor adc) which is included in the patterntrig,  bit 2 = PMTs hit are not enough to match trigger configuration, bit 3 = pattern trig does not match trigger configuration
265    Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise    Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise
266    //    //
267    Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07    Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07
268      Float_t GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl); // gf Apr 07
269    
270    Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  wm feb 08    Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  wm feb 08
271    
# Line 156  class ToFLevel2 : public TObject { Line 281  class ToFLevel2 : public TObject {
281    void GetLevel2Struct(cToFLevel2 *) const;    void GetLevel2Struct(cToFLevel2 *) const;
282    //    //
283    ToFTrkVar *GetToFTrkVar(Int_t notrack);    ToFTrkVar *GetToFTrkVar(Int_t notrack);
284      ToFTrkVar* GetToFStoredTrack(Int_t seqno);///< returns pointer to the track set related to the seqno number
285    ToFPMT *GetToFPMT(Int_t nohit);    ToFPMT *GetToFPMT(Int_t nohit);
286    Int_t GetPMTid(Int_t gg, Int_t hh);    static Int_t GetPMTid(Int_t gg, Int_t hh);
287    TString GetPMTName(Int_t ind);    static TString GetPMTName(Int_t ind);
288        
289    Int_t GetPlaneIndex(Int_t pmt_id);    static Int_t GetPlaneIndex(Int_t pmt_id);
290    void GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]);    void GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]);
291    void GetPMTIndex(Int_t pmt_id, Int_t &gg, Int_t &hh);    static void GetPMTIndex(Int_t pmt_id, Int_t &gg, Int_t &hh);
292    
293    // gf Apr 07    // gf Apr 07
294    void GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07    void GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
295    TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);    void GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
296    Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07    static TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);
297    void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07    Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07 //EMXX
298      Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin); // wm jun 2008
299      static void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07
300    void GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright); // gf Apr 07    void GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright); // gf Apr 07
301    void GetPaddleGeometry(Int_t plane, Int_t paddle, Float_t &xleft, Float_t &xright, Float_t &yleft, Float_t &yright); // gf Apr 07    void GetPaddleGeometry(Int_t plane, Int_t paddle, Float_t &xleft, Float_t &xright, Float_t &yleft, Float_t &yright); // gf Apr 07
302    Int_t GetPaddleid(Int_t plane, Int_t paddle);    Int_t GetPaddleid(Int_t plane, Int_t paddle);
# Line 176  class ToFLevel2 : public TObject { Line 304  class ToFLevel2 : public TObject {
304    Int_t GetNPaddle(Int_t plane);    Int_t GetNPaddle(Int_t plane);
305    //    //
306    //    //
307        //
308      Int_t Process(TrkLevel2 *trk, TrigLevel2 *trg, GL_RUN *run, OrbitalInfo *orb, Bool_t force); // Emiliano
309    
310      //
311      //
312      bool bit(int decimal, char pos);
313      bool checkPMT(TString givenpmt);
314      bool checkPMTpatternPMThit(TrigLevel2 *trg, int &pmtpattern, int &pmtnosignal);
315      bool checkPMTpmttrig(TrigLevel2 *trg);
316      void printPMT();
317    
318    //    //
319    // constructor    // constructor
320    //    //
# Line 218  class ToFLevel2 : public TObject { Line 356  class ToFLevel2 : public TObject {
356    
357      Int_t  GetToFPlaneID(Int_t ip);      Int_t  GetToFPlaneID(Int_t ip);
358      Int_t  GetToFPlaneIndex(Int_t plane_id);      Int_t  GetToFPlaneIndex(Int_t plane_id);
359      Bool_t HitPaddle(Int_t ,Int_t);      Bool_t HitPaddle(Int_t ,Int_t); // EMXX
360      Int_t  GetNHitPaddles(Int_t plane);      Int_t  GetNHitPaddles(Int_t plane);
361        Int_t  GetTrueNHitPaddles(Int_t plane);
362      void Clear(Option_t *t="");      void Clear(Option_t *t="");
363    
364        //
365        TClonesArray *GetTrackArray(){return ToFTrk;} ///< returns a pointer to the track related variables array
366        TClonesArray** GetPointerToTrackArray(){return &ToFTrk;}///< returns pointer to pointer to the track array
367        void SetTrackArray(TClonesArray *track);///<set pointer to the track array
368    
369    
370      //      //
371      ClassDef(ToFLevel2,3);      ClassDef(ToFLevel2,6);
372  };  };
373    
374  #endif  #endif
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.23