ToFLevel2/inc/ToFLevel2.h

/**
 * \file ToFLevel2.h
 * \author Gianfranca DeRosa / Wolfgang Menn / Rita Carbone with E. M. supervision
 */

#ifndef ToFLevel2_h
#define ToFLevel2_h
//
#include <TObject.h>
#include <TArrayI.h>
#include <TArrayF.h>
#include <TClonesArray.h>

#include <math.h> // EMILIANO
#include <iostream> // from ToFLevel2.cpp
#include <fstream> // Emiliano
#include <sstream> // Emiliano 
#include <string> // Emiliano


#include <ToFStruct.h>

#include <TrkLevel2.h> // Emiliano
#include <TrigLevel2.h> // Emiliano 
#include <GLTables.h> // Emiliano
#include <OrbitalInfo.h> // Emiliano
#include <ToFCore.h> // Emiliano
#include <physics/tof/TofEvent.h>

//
// Declaration of the core fortran routines
//
#define tofl2com tofl2com_
extern "C" int tofl2com();
#define toftrk toftrk_
extern "C" int toftrk();
#define rdtofcal rdtofcal_
//extern "C" int rdtofcal(char [], int *);
extern "C" int rdtofcal(const char *, int *);

//
// class which contains track related variables
//
#define ZTOF11 53.74
#define ZTOF12 53.04
#define ZTOF21 23.94
#define ZTOF22 23.44
#define ZTOF31 -23.49
#define ZTOF32 -24.34


class ToFGeom : public TObject {

private:
  TArrayI  ePlane, eXY;

public:
  ToFGeom() {
    int plane[24] = { 
      0, 0, 0, 0, 0, 0, 0, 0,
      1, 1, 1, 1, 1, 1, 
      2, 2,             
      3, 3,             
      4, 4, 4,
      5, 5, 5
    };
    int plXY[6]= {  2,  1,  1,  2,  2,  1 }; // X==1, Y==2 */
    ePlane.Set(24,plane);
    eXY.Set(6,plXY);
  }

  int GetPad(   int idpmt)  { return (int)((idpmt+0.5)/2.); }
  int GetPlane( int idpmt)  { return ePlane[ GetPad(idpmt) ];  }
  int GetXY(    int idpmt)  { return eXY[ GetPlane(idpmt) ]; }

   ClassDef(ToFGeom,1);

};


/**
 * \brief  Class which contains the PMT data
 *
 * If there is a valid ADC or a TDC value (value<4095) for a PMT, both ADC and TDC data
 * are stored in the PMT class.
 * Look in the ToFLevel2Ex.cxx example in the repository how to read the PMT class. 
 */
class ToFPMT : public TObject {

 private:

 public:
    Int_t pmt_id;     ///<the identification number of the PMT from 0 to 47
    Float_t adc;      ///<raw ADC value for this PMT
    Float_t tdc;      ///<raw TDC value for this PMT
    Float_t tdc_tw;   ///<time-walk corrected TDC value for this PMT
    Float_t l0flag_adc; ///< warning flags from unpacking
    Float_t l0flag_tdc; ///< warning flags from unpacking
    //
    ToFPMT();
    ToFPMT(const ToFPMT&);
    //
    ToFPMT* GetToFPMT(){return this;};
    void Clear(Option_t *t="");

    ClassDef(ToFPMT,3);
};

/**
 * \brief  Class used to calibrate adc to dEdx for each PMT
 *
 * Class used to calibrate adc to dEdx for each PMT
 */
class ToFdEdx : public TObject {

 private:
    //
    ToFGeom  eGeom;  // ToF geometry
    //
    Float_t  adc_he;
    TArrayF *eDEDXpmt;  // 0-47  pmt dEdx
    // parameters:
    TArrayF PMTsat;  // 0-47  saturation parameters
    Float_t adc[48];
    //

    TArrayF parAtt[48];      // 48 x 6
    TArrayF parPos[48];      // 48 x 4
    TArrayF parDesatBB[48];  // 48 x 3
    TArrayF parBBneg[48];    // 48 x 3
    TArrayF  parBBpos;    // 48 x 1
 
    double f_adcPC( float x );
    double f_BB( TArrayF &p, float x );
    double f_BB5B( float x );
    double f_att( TArrayF &p, float x ) ;
    double f_att5B( float x );
    double f_desatBB( TArrayF &p, float x );  
    double f_desatBB5B( float x );
    double f_pos( TArrayF &p, float x );
    double f_pos5B( float x );
    float Get_adc_he( int id, float pl_x[6], float pl_y[6]);

    Bool_t conn[12]; 

    UInt_t ts[12];
    UInt_t te[12];    
   
    
 public:
    ToFdEdx();   // class constructor
    ~ToFdEdx();  // class distructor
    //
    void Clear(Option_t *option="");
    void Delete(Option_t *option="");

    void Init(pamela::tof::TofEvent *tofl0 );  // init parameters
    void Init(Int_t i, Int_t j, Float_t adce);
    void Define_PMTsat();

    void ReadParAtt(            const char *fname  );
    void ReadParPos(            const char *fname  );
    void ReadParBBneg(          const char *fname  );
    void ReadParBBpos(          const char *fname  );
    void ReadParDesatBB(        const char *fname  );

    void CheckConnectors(UInt_t atime, GL_PARAM *glparam, TSQLServer *dbc);

    void Process( UInt_t atime, Float_t betamean, Float_t *xtr_tof, Float_t *ytr_tof, Int_t exitat=-1); //
    void Print(Option_t *option="");

    //    Float_t GetdEdx_pmt(Int_t ipmt) { return (Float_t)eDEDXpmt[ipmt]; }  // 0-47 dEdx for each PMT for tracked events
    Float_t GetdEdx_pmt(Int_t ipmt) { return eDEDXpmt->At(ipmt); }  // 0-47 dEdx for each PMT for tracked events
    //
    ToFdEdx* GetToFdEdx(){return this;};
    ClassDef(ToFdEdx,3);
};


/**
 * \brief Class which contains the tracker related variables
 *
 * We can use the ToF standalone to find hitted paddles, calculate beta, etc..
 * These results are then stored with the "trkseqno" = -1.
 * If we use the track from the tracker, then the penetration points in the
 * scintillators are calculated, which defines the hitted paddles. For these paddles
 * we calculate then all the output.
 * Note: The artificial ADC values are stored as dEdx in the output, the dEdx will be
 * by definition = 1. However, the artificial TDC values are just used internally
 * and not stored in the output. But one can see in both cases which PMT has artificial
 * values using "adcflag" and "tdcflag".
 * Look in the ToFLevel2Ex.cxx example in the repository how to read the tracker related
 * variables. 
 */
class ToFTrkVar : public TObject {

 private:

 public:
  //
  Int_t trkseqno; ///< tracker sequ. number: -1=ToF standalone, 0=first Tracker track, ...
  //
  Int_t npmttdc;  ///<number of the TDC measurements used to evaluate beta
  TArrayI pmttdc;  ///<contains the ID (0..47) for the PMT used to evaluate beta
  TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value

/**
 * \brief beta, 12 measurements for the 12  combinations, beta[13] is modified weighted mean
 *
 * The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for
 * S2-S3 and S1-S2. 
 * The calculation of beta[13] is now modified:
 * We check the individual weights for artificial TDC values, then calculate
 * am mean beta for the first time. In a second step we loop again through
 * the single measurements, checking for the residual from the mean
 * The cut on the residual reject measurements > "x"-sigma. A chi2 value is
 * calculated, furthermore a "quality" value by adding the weights which
 * are finally used. If all measurements are taken, "quality" will be = 505.
 * A chi2 cut around 3-4 and a quality-cut > 400 is needed for clean beta
 * The Level2 beta[12] which is derived in the fortran routines uses: 10.,200.,20.
 * This is not a very high quality measurement. One can re-calculate a new beta[13]
 * using the L2-method "CalcBeta"
 */
  Float_t beta[13]; 
  //
  Int_t npmtadc;  ///<number of the ADC measurements used for dEdx evaluation
  TArrayI pmtadc; ///<contains the ID (0..47) for the PMT used to evaluate dEdx
  TArrayI adcflag; ///<flag for artificial ADCs, "0" if normal ADC value
  TArrayF dedx;    ///<energy loss for this PMT in mip
  //
  Float_t xtofpos[3];  ///<x-measurement using the TDC values and the calibration from S12, S21, S32
  Float_t ytofpos[3];  ///<x-measurement using the TDC values and the calibration from S11, S22, S31
  //
  Float_t xtr_tof[6];  ///<x-measurement in the ToF layers from tracker
  Float_t ytr_tof[6];  ///<x-measurement in the ToF layers from tracker
  //
  ToFTrkVar();
  ToFTrkVar(const ToFTrkVar&);

  ToFTrkVar* GetToFTrkVar(){return this;};
  void Clear(Option_t *t="");

  Float_t CalcBeta(Float_t resmax, Float_t qualitycut, Float_t chi2cut);  //  Elena 2015

  ClassDef(ToFTrkVar,2);
  //
};

/**
 * \brief  Class to describe ToF LEVEL2 data
 *
 */

class ToFLevel2 : public TObject {
 private:

 public:
  //
  TClonesArray *PMT; ///<class needed to store PMT hit informations
  TClonesArray *ToFTrk; ///<track related variable class
  Int_t tof_j_flag[6];  ///<number of hitted paddle(s) for each ToF layer: flag = flag + 2**(paddlenumber-1)
  //
  Int_t unpackError;///< zero if no errors
  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
  Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise
  //
  Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07
  Float_t GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl); // gf Apr 07

  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 chi2cut, Float_t qualitycut);   // wm feb 08
  //
  // methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables
  //
  Int_t ntrk(){return ToFTrk->GetEntries();};
  Int_t npmt(){return PMT->GetEntries();};

  //
  void GetLevel2Struct(cToFLevel2 *) const;
  //
  ToFTrkVar *GetToFTrkVar(Int_t notrack);
  ToFTrkVar* GetToFStoredTrack(Int_t seqno);///< returns pointer to the track set related to the seqno number
  ToFPMT *GetToFPMT(Int_t nohit);
  static Int_t GetPMTid(Int_t gg, Int_t hh);
  static TString GetPMTName(Int_t ind);
  
  static Int_t GetPlaneIndex(Int_t pmt_id);
  void GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]);
  static void GetPMTIndex(Int_t pmt_id, Int_t &gg, Int_t &hh);

  // gf Apr 07
  void GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
  void GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
  static TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);
  Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07 //EMXX
  Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin); // wm jun 2008
  static void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07
  void GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright); // 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
  Int_t GetPaddleid(Int_t plane, Int_t paddle);
  void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle);
  Int_t GetNPaddle(Int_t plane);
  //
  //
  //
  Int_t Process(TrkLevel2 *trk, TrigLevel2 *trg, GL_RUN *run, OrbitalInfo *orb, Bool_t force); // Emiliano

  //
  //
  bool bit(int decimal, char pos);
  bool checkPMT(TString givenpmt);
  bool checkPMTpatternPMThit(TrigLevel2 *trg, int &pmtpattern, int &pmtnosignal);
  bool checkPMTpmttrig(TrigLevel2 *trg);
  void printPMT();

  //
  // constructor
  //
  ToFLevel2();
  ~ToFLevel2(){Delete();}; //ELENA
  void Delete(Option_t *t=""); //ELENA
  void Set();//ELENA
  //
  //
  ToFLevel2*   GetToFLevel2(){return this;};

/**
 * Method to get the z-position of the 6 TOF layers from the plane ID
 * @param plane_id Plane ID (11 12 21 22 31 32)
 */
  Float_t      GetZTOF(Int_t plane_id){
      switch(plane_id){
      case 11: return ZTOF11;
      case 12: return ZTOF12;
      case 21: return ZTOF21;
      case 22: return ZTOF22;
      case 31: return ZTOF31;
      case 32: return ZTOF32;
      default: return 0.;
      };
  };

    //
    // Paddles position
    //
    /*
      S11 8 paddles  33.0 x 5.1 cm
      S12 6 paddles  40.8 x 5.5 cm
      S21 2 paddles  18.0 x 7.5 cm
      S22 2 paddles  15.0 x 9.0 cm
      S31 3 paddles  15.0 x 6.0 cm
      S32 3 paddles  18.0 x 5.0 cm
    */

    Int_t  GetToFPlaneID(Int_t ip);
    Int_t  GetToFPlaneIndex(Int_t plane_id);
    Bool_t HitPaddle(Int_t ,Int_t); // EMXX 
    Int_t  GetNHitPaddles(Int_t plane);
    Int_t  GetTrueNHitPaddles(Int_t plane);
    void Clear(Option_t *t="");

    //
    TClonesArray *GetTrackArray(){return ToFTrk;} ///< returns a pointer to the track related variables array
    TClonesArray** GetPointerToTrackArray(){return &ToFTrk;}///< returns pointer to pointer to the track array
    void SetTrackArray(TClonesArray *track);///<set pointer to the track array


    //
    ClassDef(ToFLevel2,7);
};

#endif

1	pam-de	1.14	/**
2			* \file ToFLevel2.h
3	mocchiut	1.24	* \author Gianfranca DeRosa / Wolfgang Menn / Rita Carbone with E. M. supervision
4	pam-de	1.14	*/
5
6	mocchiut	1.1	#ifndef ToFLevel2_h
7			#define ToFLevel2_h
8			//
9			#include <TObject.h>
10	mocchiut	1.4	#include <TArrayI.h>
11			#include <TArrayF.h>
12	mocchiut	1.1	#include <TClonesArray.h>
13	pam-fi	1.6
14	mocchiut	1.18	#include <math.h> // EMILIANO
15	mocchiut	1.21	#include <iostream> // from ToFLevel2.cpp
16			#include <fstream> // Emiliano
17			#include <sstream> // Emiliano
18			#include <string> // Emiliano
19
20	mocchiut	1.18
21	pam-fi	1.6	#include <ToFStruct.h>
22	pam-de	1.14
23	mocchiut	1.21	#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	carbone	1.23	#include <physics/tof/TofEvent.h>
29
30	mocchiut	1.21	//
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	mocchiut	1.28	//extern "C" int rdtofcal(char [], int *);
39			extern "C" int rdtofcal(const char , int );
40	pam-de	1.14
41	mocchiut	1.1	//
42			// class which contains track related variables
43			//
44			#define ZTOF11 53.74
45			#define ZTOF12 53.04
46			#define ZTOF21 23.94
47			#define ZTOF22 23.44
48			#define ZTOF31 -23.49
49			#define ZTOF32 -24.34
50
51	mocchiut	1.4
52	carbone	1.23	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	pam-de	1.14	/**
82			* \brief Class which contains the PMT data
83			*
84	pam-de	1.15	* If there is a valid ADC or a TDC value (value<4095) for a PMT, both ADC and TDC data
85			* are stored in the PMT class.
86			* Look in the ToFLevel2Ex.cxx example in the repository how to read the PMT class.
87	pam-de	1.14	*/
88	mocchiut	1.4	class ToFPMT : public TObject {
89	pam-de	1.14
90	mocchiut	1.4	private:
91
92			public:
93	pam-de	1.14	Int_t pmt_id; ///<the identification number of the PMT from 0 to 47
94	mocchiut	1.25	Float_t adc; ///<raw ADC value for this PMT
95	mocchiut	1.17	Float_t tdc; ///<raw TDC value for this PMT
96	pam-de	1.15	Float_t tdc_tw; ///<time-walk corrected TDC value for this PMT
97	mocchiut	1.30	Float_t l0flag_adc; ///< warning flags from unpacking
98			Float_t l0flag_tdc; ///< warning flags from unpacking
99	mocchiut	1.4	//
100			ToFPMT();
101			ToFPMT(const ToFPMT&);
102			//
103			ToFPMT* GetToFPMT(){return this;};
104	mocchiut	1.19	void Clear(Option_t *t="");
105	mocchiut	1.4
106	mocchiut	1.30	ClassDef(ToFPMT,3);
107	mocchiut	1.4	};
108
109	carbone	1.23	/**
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	mocchiut	1.28	TArrayF *eDEDXpmt; // 0-47 pmt dEdx
122	carbone	1.23	// parameters:
123			TArrayF PMTsat; // 0-47 saturation parameters
124	mocchiut	1.24	Float_t adc[48];
125			//
126	carbone	1.23
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	mocchiut	1.29	~ToFdEdx(); // class distructor
153	carbone	1.23	//
154			void Clear(Option_t *option="");
155	mocchiut	1.29	void Delete(Option_t *option="");
156	carbone	1.23
157	mocchiut	1.24	void Init(pamela::tof::TofEvent *tofl0 ); // init parameters
158			void Init(Int_t i, Int_t j, Float_t adce);
159	carbone	1.23	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	mocchiut	1.27	void Process( UInt_t atime, Float_t betamean, Float_t xtr_tof, Float_t ytr_tof, Int_t exitat=-1); //
170	carbone	1.23	void Print(Option_t *option="");
171
172	mocchiut	1.28	// 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	carbone	1.23	//
175			ToFdEdx* GetToFdEdx(){return this;};
176	mocchiut	1.27	ClassDef(ToFdEdx,3);
177	carbone	1.23	};
178
179	mocchiut	1.4
180	pam-de	1.14	/**
181	pam-de	1.15	* \brief Class which contains the tracker related variables
182	pam-de	1.14	*
183	pam-de	1.15	* We can use the ToF standalone to find hitted paddles, calculate beta, etc..
184			* These results are then stored with the "trkseqno" = -1.
185			* If we use the track from the tracker, then the penetration points in the
186			* scintillators are calculated, which defines the hitted paddles. For these paddles
187			* we calculate then all the output.
188			* Note: The artificial ADC values are stored as dEdx in the output, the dEdx will be
189			* by definition = 1. However, the artificial TDC values are just used internally
190			* and not stored in the output. But one can see in both cases which PMT has artificial
191			* values using "adcflag" and "tdcflag".
192			* Look in the ToFLevel2Ex.cxx example in the repository how to read the tracker related
193			* variables.
194	pam-de	1.14	*/
195	mocchiut	1.1	class ToFTrkVar : public TObject {
196	pam-de	1.14
197	mocchiut	1.1	private:
198
199			public:
200	mocchiut	1.13	//
201	pam-de	1.14	Int_t trkseqno; ///< tracker sequ. number: -1=ToF standalone, 0=first Tracker track, ...
202	mocchiut	1.1	//
203	pam-de	1.15	Int_t npmttdc; ///<number of the TDC measurements used to evaluate beta
204			TArrayI pmttdc; ///<contains the ID (0..47) for the PMT used to evaluate beta
205			TArrayI tdcflag; ///<flag for artificial TDC, "0" if normal TDC value
206
207			/**
208	mocchiut	1.20	* \brief beta, 12 measurements for the 12 combinations, beta[13] is modified weighted mean
209	pam-de	1.15	*
210			* The 12 measurements are S11-S31, S11-S32, S12-S31, S12-S32, and then analogue for
211			* S2-S3 and S1-S2.
212	mocchiut	1.20	* The calculation of beta[13] is now modified:
213			* We check the individual weights for artificial TDC values, then calculate
214			* am mean beta for the first time. In a second step we loop again through
215			* the single measurements, checking for the residual from the mean
216			* The cut on the residual reject measurements > "x"-sigma. A chi2 value is
217			* calculated, furthermore a "quality" value by adding the weights which
218			* are finally used. If all measurements are taken, "quality" will be = 505.
219			* A chi2 cut around 3-4 and a quality-cut > 400 is needed for clean beta
220			* The Level2 beta[12] which is derived in the fortran routines uses: 10.,200.,20.
221			* This is not a very high quality measurement. One can re-calculate a new beta[13]
222			* using the L2-method "CalcBeta"
223	pam-de	1.15	*/
224			Float_t beta[13];
225			//
226			Int_t npmtadc; ///<number of the ADC measurements used for dEdx evaluation
227			TArrayI pmtadc; ///<contains the ID (0..47) for the PMT used to evaluate dEdx
228			TArrayI adcflag; ///<flag for artificial ADCs, "0" if normal ADC value
229			TArrayF dedx; ///<energy loss for this PMT in mip
230	mocchiut	1.4	//
231	pam-de	1.15	Float_t xtofpos[3]; ///<x-measurement using the TDC values and the calibration from S12, S21, S32
232			Float_t ytofpos[3]; ///<x-measurement using the TDC values and the calibration from S11, S22, S31
233	mocchiut	1.16	//
234			Float_t xtr_tof[6]; ///<x-measurement in the ToF layers from tracker
235			Float_t ytr_tof[6]; ///<x-measurement in the ToF layers from tracker
236	mocchiut	1.1	//
237			ToFTrkVar();
238			ToFTrkVar(const ToFTrkVar&);
239	pam-de	1.14
240	mocchiut	1.1	ToFTrkVar* GetToFTrkVar(){return this;};
241	mocchiut	1.19	void Clear(Option_t *t="");
242	mocchiut	1.4
243	pam-fi	1.35	Float_t CalcBeta(Float_t resmax, Float_t qualitycut, Float_t chi2cut); // Elena 2015
244
245	mocchiut	1.36	ClassDef(ToFTrkVar,2);
246	mocchiut	1.1	//
247			};
248
249	pam-de	1.14	/**
250			* \brief Class to describe ToF LEVEL2 data
251			*
252			*/
253
254	mocchiut	1.1	class ToFLevel2 : public TObject {
255			private:
256	pam-de	1.14
257	mocchiut	1.1	public:
258	mocchiut	1.4	//
259	pam-de	1.14	TClonesArray *PMT; ///<class needed to store PMT hit informations
260			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)
262	mocchiut	1.17	//
263	mocchiut	1.30	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	mocchiut	1.17	Int_t default_calib; ///< one if the default calibration has been used to process the data, zero otherwise
266	mocchiut	1.1	//
267	mocchiut	1.16	Float_t GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl); // gf Apr 07
268	mocchiut	1.32	Float_t GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl); // gf Apr 07
269	mocchiut	1.20
270			Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut); // wm feb 08
271
272			//
273			// Float_t CalcBeta(Int_t notrack, Float_t resmax, Float_t chi2cut, Float_t qualitycut); // wm feb 08
274	mocchiut	1.1	//
275			// methods to make life simplier during the analysis, returns a pointer to the ToFTrkVar class containing track related variables
276			//
277			Int_t ntrk(){return ToFTrk->GetEntries();};
278	mocchiut	1.4	Int_t npmt(){return PMT->GetEntries();};
279	pam-de	1.14
280	mocchiut	1.1	//
281	pam-fi	1.6	void GetLevel2Struct(cToFLevel2 *) const;
282			//
283	pam-de	1.14	ToFTrkVar *GetToFTrkVar(Int_t notrack);
284	pam-fi	1.34	ToFTrkVar* GetToFStoredTrack(Int_t seqno);///< returns pointer to the track set related to the seqno number
285	pam-de	1.14	ToFPMT *GetToFPMT(Int_t nohit);
286	pam-fi	1.35	static Int_t GetPMTid(Int_t gg, Int_t hh);
287			static TString GetPMTName(Int_t ind);
288	mocchiut	1.16
289	pam-fi	1.35	static Int_t GetPlaneIndex(Int_t pmt_id);
290	mocchiut	1.4	void GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]);
291	pam-fi	1.35	static void GetPMTIndex(Int_t pmt_id, Int_t &gg, Int_t &hh);
292	mocchiut	1.16
293			// gf Apr 07
294			void GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
295	mocchiut	1.32	void GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning); // gf Apr 07
296	pam-fi	1.35	static TString GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt);
297	mocchiut	1.31	Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane); // gf Apr 07 //EMXX
298	pamela	1.22	Int_t GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin); // wm jun 2008
299	pam-fi	1.35	static void GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle); // gf Apr 07
300	mocchiut	1.16	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
302			Int_t GetPaddleid(Int_t plane, Int_t paddle);
303			void GetPaddlePlane(Int_t padid, Int_t &plane, Int_t &paddle);
304			Int_t GetNPaddle(Int_t plane);
305			//
306	mocchiut	1.20	//
307	mocchiut	1.21	//
308			Int_t Process(TrkLevel2 trk, TrigLevel2 trg, GL_RUN run, OrbitalInfo orb, Bool_t force); // Emiliano
309	mocchiut	1.30
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	mocchiut	1.1	//
319			// constructor
320			//
321			ToFLevel2();
322	mocchiut	1.13	~ToFLevel2(){Delete();}; //ELENA
323	mocchiut	1.19	void Delete(Option_t *t=""); //ELENA
324	mocchiut	1.13	void Set();//ELENA
325	mocchiut	1.1	//
326			//
327			ToFLevel2* GetToFLevel2(){return this;};
328	pam-de	1.14
329			/**
330	pam-de	1.15	* Method to get the z-position of the 6 TOF layers from the plane ID
331	pam-de	1.14	* @param plane_id Plane ID (11 12 21 22 31 32)
332			*/
333	mocchiut	1.4	Float_t GetZTOF(Int_t plane_id){
334	mocchiut	1.1	switch(plane_id){
335			case 11: return ZTOF11;
336			case 12: return ZTOF12;
337			case 21: return ZTOF21;
338			case 22: return ZTOF22;
339			case 31: return ZTOF31;
340			case 32: return ZTOF32;
341			default: return 0.;
342			};
343	mocchiut	1.4	};
344	pam-de	1.14
345	mocchiut	1.1	//
346			// Paddles position
347			//
348			/*
349			S11 8 paddles 33.0 x 5.1 cm
350			S12 6 paddles 40.8 x 5.5 cm
351			S21 2 paddles 18.0 x 7.5 cm
352			S22 2 paddles 15.0 x 9.0 cm
353			S31 3 paddles 15.0 x 6.0 cm
354			S32 3 paddles 18.0 x 5.0 cm
355			*/
356	pam-de	1.14
357	mocchiut	1.1	Int_t GetToFPlaneID(Int_t ip);
358			Int_t GetToFPlaneIndex(Int_t plane_id);
359	mocchiut	1.31	Bool_t HitPaddle(Int_t ,Int_t); // EMXX
360	mocchiut	1.1	Int_t GetNHitPaddles(Int_t plane);
361	mocchiut	1.30	Int_t GetTrueNHitPaddles(Int_t plane);
362	mocchiut	1.19	void Clear(Option_t *t="");
363	pam-ts	1.33
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	pam-fi	1.34	void SetTrackArray(TClonesArray *track);///<set pointer to the track array
368	pam-ts	1.33
369
370	mocchiut	1.1	//
371	mocchiut	1.36	ClassDef(ToFLevel2,7);
372	mocchiut	1.1	};
373
374			#endif
375	pam-de	1.14