/[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.16 by mocchiut, Mon Apr 30 15:46:28 2007 UTC revision 1.22 by pamela, Fri Jun 20 12:22:11 2008 UTC
# Line 11  Line 11 
11  #include <TArrayF.h>  #include <TArrayF.h>
12  #include <TClonesArray.h>  #include <TClonesArray.h>
13    
14    #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    //
29    // Declaration of the core fortran routines
30    //
31    #define tofl2com tofl2com_
32    extern "C" int tofl2com();
33    #define toftrk toftrk_
34    extern "C" int toftrk();
35    #define rdtofcal rdtofcal_
36    extern "C" int rdtofcal(char [], int *);
37    
38  //  //
39  // class which contains track related variables  // class which contains track related variables
# Line 39  class ToFPMT : public TObject { Line 60  class ToFPMT : public TObject {
60   public:   public:
61      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
62      Float_t adc;      ///<raw ADC value for this PMT      Float_t adc;      ///<raw ADC value for this PMT
63        Float_t tdc;      ///<raw TDC value for this PMT
64      Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT      Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT
65      //      //
66      ToFPMT();      ToFPMT();
67      ToFPMT(const ToFPMT&);      ToFPMT(const ToFPMT&);
68      //      //
69      ToFPMT* GetToFPMT(){return this;};      ToFPMT* GetToFPMT(){return this;};
70      void Clear();      void Clear(Option_t *t="");
71    
72    
73    
74      ClassDef(ToFPMT,1);      ClassDef(ToFPMT,2);
75  };  };
76    
77    
# Line 81  class ToFTrkVar : public TObject { Line 103  class ToFTrkVar : public TObject {
103    TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value    TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value
104    
105  /**  /**
106   * \brief beta, 12 measurements for the 12  combinations, beta[13] is weighted mean   * \brief beta, 12 measurements for the 12  combinations, beta[13] is modified weighted mean
107   *   *
108   * The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for   * The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for
109   * S2-S3 and S1-S2.   * S2-S3 and S1-S2.
110   * In the moment all measurements are taken and the weighted mean is calculated.   * The calculation of beta[13] is now modified:
111   * Note that the weights are just simple overall results for S1-S3, S2-S3, and S1-S2.   * We check the individual weights for artificial TDC values, then calculate
112   * Artificial measurments are not treated correct, (since there is only one real   * am mean beta for the first time. In a second step we loop again through
113   * measurment the weight should be different then for two meassurments).   * the single measurements, checking for the residual from the mean
114   * The beta calculation will be improved in the next release.       * The cut on the residual reject measurements > "x"-sigma. A chi2 value is
115     * calculated, furthermore a "quality" value by adding the weights which
116     * are finally used. If all measurements are taken, "quality" will be = 505.
117     * A chi2 cut around 3-4 and a quality-cut > 400 is needed for clean beta
118     * The Level2 beta[12] which is derived in the fortran routines uses: 10.,200.,20.
119     * This is not a very high quality measurement. One can re-calculate a new beta[13]
120     * using the L2-method "CalcBeta"
121   */   */
122    Float_t beta[13];    Float_t beta[13];
123    //    //
# Line 108  class ToFTrkVar : public TObject { Line 136  class ToFTrkVar : public TObject {
136    ToFTrkVar(const ToFTrkVar&);    ToFTrkVar(const ToFTrkVar&);
137    
138    ToFTrkVar* GetToFTrkVar(){return this;};    ToFTrkVar* GetToFTrkVar(){return this;};
139    void Clear();    void Clear(Option_t *t="");
140    
141    ClassDef(ToFTrkVar,1);    ClassDef(ToFTrkVar,1);
142    //    //
# Line 127  class ToFLevel2 : public TObject { Line 155  class ToFLevel2 : public TObject {
155    TClonesArray *PMT; ///<class needed to store PMT hit informations    TClonesArray *PMT; ///<class needed to store PMT hit informations
156    TClonesArray *ToFTrk; ///<track related variable class    TClonesArray *ToFTrk; ///<track related variable class
157    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)
158      //
159    Int_t unpackError;    Int_t unpackError;///< zero if no error presente
160      Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise
161    //    //
162    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
163    
164      Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  wm feb 08
165    
166      //
167    //  Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t chi2cut, Float_t qualitycut);   // wm feb 08
168    //    //
169    // 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
170    //    //
# Line 153  class ToFLevel2 : public TObject { Line 187  class ToFLevel2 : public TObject {
187    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
188    TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);    TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);
189    Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07    Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07
190      Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin); // wm jun 2008
191    void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07    void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07
192    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
193    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
# Line 160  class ToFLevel2 : public TObject { Line 195  class ToFLevel2 : public TObject {
195    void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle);    void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle);
196    Int_t GetNPaddle(Int_t plane);    Int_t GetNPaddle(Int_t plane);
197    //    //
198      //
199      //
200      Int_t Process(TrkLevel2 *trk, TrigLevel2 *trg, GL_RUN *run, OrbitalInfo *orb, Bool_t force); // Emiliano
201        
202    //    //
203    // constructor    // constructor
204    //    //
205    ToFLevel2();    ToFLevel2();
206    ~ToFLevel2(){Delete();}; //ELENA    ~ToFLevel2(){Delete();}; //ELENA
207    void Delete(); //ELENA    void Delete(Option_t *t=""); //ELENA
208    void Set();//ELENA    void Set();//ELENA
209    //    //
210    //    //
# Line 204  class ToFLevel2 : public TObject { Line 242  class ToFLevel2 : public TObject {
242      Int_t  GetToFPlaneIndex(Int_t plane_id);      Int_t  GetToFPlaneIndex(Int_t plane_id);
243      Bool_t HitPaddle(Int_t ,Int_t);      Bool_t HitPaddle(Int_t ,Int_t);
244      Int_t  GetNHitPaddles(Int_t plane);      Int_t  GetNHitPaddles(Int_t plane);
245      void Clear();      void Clear(Option_t *t="");
246      //      //
247      ClassDef(ToFLevel2,2);      ClassDef(ToFLevel2,4);
248  };  };
249    
250  #endif  #endif
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.23