CaloNuclei/inc/CaloNuclei.h



/**
 * \file CaloNuclei.h
 * \author Emiliano Mocchiutti
 */
#ifndef calonuclei_h
#define calonuclei_h

#include <PamLevel2.h>

#include <TTree.h>
#include <TFriendElement.h>
#include <TChain.h>
#include <TFile.h>
#include <TList.h>
#include <TKey.h>
#include <TSystemFile.h>
#include <TSystemDirectory.h>
#include <TSQLServer.h>

#include <iostream>

using namespace std;

/**
 *
 * Class to store and calculate variables useful for nuclei analysis
 */
class CaloNuclei : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    Bool_t usetrack;
    const char*  trkAlg;
    
    //
    // needed to avoid reprocessing the same event over and over to obtain the variables
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    Int_t tr;
    Int_t sntr;
    //
    Bool_t usepl18x;
    //
    Int_t interplane; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Int_t N; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
    Int_t R; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
    Int_t UN; ///< Number of dE/dx measurements really used to calculate qpremeanN
    Float_t preq; ///< Energy release (MIP) up to the interaction plane (included)
    Float_t postq; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
    Float_t stdedx1; ///< Energy release (MIP) on the first Silicon detector (Y EVEN) around the strip with maximum energy release (no track information).
    Float_t dedx1; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
    Float_t dedx3; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
    Float_t qpremean; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
    Float_t qpremeanN; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
    Float_t qNmin1; ///< Truncated mean (MIP) along the track using N-1 measurements before of the interaction plane
    Float_t maxrel; ///<Energy maximum release on first Calorimeter plane (dedx of strip with maximum release)
    Float_t ethr; ///< Threshold (MIP) needed to find the interaction plane
    Bool_t multhit; ///< True if the interaction plane has been determined by multiple hit counting
    Bool_t gap; ///< True if determining the interaction plane a big (>5 planes) gap has been found between a point and another along the track
    Float_t charge_siegen1; 
    Float_t ZCalo_dedx_b;  //Z from Calo using dedx in first Calorimeter plane vs. beta
    Float_t ZCalo_maxrel_b; //Z from Calo using maximum release in first Calorimeter plane vs. beta 
    Float_t ZCalo_dedx_defl; //Z from Calo using dedx in first Calorimeter plane vs. rigidity
    Float_t ZCalo_Nmin1_defl;  //Z from Calo using truncated mean on N-1 Calorimeter planes (plane N+1 is the interaction plane) vs. rigidity
   //
    Float_t qNmin1_w;
    Int_t S2;
 public:
    //
    //
    //char* version();

    Int_t Get_interplane(){ Process(); return interplane;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Int_t Get_N(){ return N;}; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
    Int_t Get_UsedN(){ return UN;}; ///< Number of dE/dx measurements really used to calculate qpremeanN
    Int_t Get_R(){ return R;}; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3 
    Float_t Get_preq(){ Process(); return preq;}; ///< Energy release (MIP) up to the interaction plane (included)
    Float_t Get_postq(){ Process(); return postq;}; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
    Float_t Get_StdEdx1(){ Process(); return stdedx1;}; ///< Energy release (MIP) on the first Silicon detector (Y EVEN) around the strip with maximum energy release (no track information, 3 strips in total).
    Float_t Get_dEdx1(){ Process(); return dedx1;}; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
    Float_t Get_dEdx3(){ Process(); return dedx3;}; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
    Float_t Get_qpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
    Float_t Get_qpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
    Float_t Get_qNmin1(){ Process(); return qNmin1;}; ///< Truncated mean (MIP) along the track using N-1 measurements before of the interaction plane
    Float_t Get_maxrel(){ Process(); return maxrel;}; ///<Energy maximum release on first Calorimeter plane (dedx of strip with maximum release)
    Float_t Get_ethr(){ Process(); return ethr;}; ///< Threshold (MIP) needed to find the interaction plane
    Bool_t IsMulthit(){ Process(); return multhit;}; ///< True if the interaction plane has been determined by multiple hit counting
    //
    Float_t Get_charge_siegen1(){ Process(); return charge_siegen1;}; ///< charge Siegen method stdedx1 vs. beta
    Float_t Get_ZCalo_dedx_b(){ Process(); return ZCalo_dedx_b;};  //Z from Calo using dedx (or StdEdx) in first Calorimeter plane vs. beta
    Float_t Get_ZCalo_maxrel_b(){ Process(); return ZCalo_maxrel_b;}; //Z from Calo using maximum release in first Calorimeter plane vs. beta 
    Float_t Get_ZCalo_dedx_defl(){ Process(); return ZCalo_dedx_defl;}; //Z from Calo using dedx in first Calorimeter plane vs. rigidity
    Float_t Get_ZCalo_Nmin1_defl(){ Process(); return ZCalo_Nmin1_defl;};  //Z from Calo using truncated mean on N-1 Calorimeter planes (plane N+1 is the interaction plane) vs. rigidity
    
    //
    void Set_N(Int_t n){ N=n;};
    void Set_R(Int_t r){ R=r;};
    //
//    CaloNuclei();
    CaloNuclei(PamLevel2 *L2,const char* alg);
    ~CaloNuclei(){ Delete(); };
    //
    void SetDebug(Bool_t d){ debug=d; };
    void UseTrack(Bool_t d){ usetrack=d; };


    void UsePlane18X(Bool_t use){usepl18x = use;};
    //
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    void Delete(Option_t *option){Delete();};
    //
    void Process(); ///< Process data for track number 0
    void Process(Int_t ntr); ///< Process data for track number ntr
    void Print();
    void Print(Option_t *option){Print();};
    //
    ClassDef(CaloNuclei,3);
};

#endif

1	mocchiut	1.1
2
3			/**
4			* \file CaloNuclei.h
5			* \author Emiliano Mocchiutti
6			*/
7			#ifndef calonuclei_h
8			#define calonuclei_h
9
10			#include <PamLevel2.h>
11
12			#include <TTree.h>
13			#include <TFriendElement.h>
14			#include <TChain.h>
15			#include <TFile.h>
16			#include <TList.h>
17			#include <TKey.h>
18			#include <TSystemFile.h>
19			#include <TSystemDirectory.h>
20			#include <TSQLServer.h>
21
22			#include <iostream>
23
24			using namespace std;
25
26			/**
27			*
28			* Class to store and calculate variables useful for nuclei analysis
29			*/
30			class CaloNuclei : public TObject {
31
32			private:
33			//
34			PamLevel2 *L2;
35	mocchiut	1.4	Bool_t debug;
36	mocchiut	1.5	Bool_t usetrack;
37	pam-fi	1.10	const char* trkAlg;
38
39	mocchiut	1.1	//
40			// needed to avoid reprocessing the same event over and over to obtain the variables
41			//
42			UInt_t OBT;
43			UInt_t PKT;
44			UInt_t atime;
45	mocchiut	1.3	Int_t tr;
46	mocchiut	1.5	Int_t sntr;
47	mocchiut	1.1	//
48	mocchiut	1.8	Bool_t usepl18x;
49			//
50	mocchiut	1.1	Int_t interplane; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
51			Int_t N; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
52	mocchiut	1.2	Int_t R; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
53	mocchiut	1.6	Int_t UN; ///< Number of dE/dx measurements really used to calculate qpremeanN
54	mocchiut	1.1	Float_t preq; ///< Energy release (MIP) up to the interaction plane (included)
55			Float_t postq; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
56	mocchiut	1.5	Float_t stdedx1; ///< Energy release (MIP) on the first Silicon detector (Y EVEN) around the strip with maximum energy release (no track information).
57	mocchiut	1.1	Float_t dedx1; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
58			Float_t dedx3; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
59			Float_t qpremean; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
60			Float_t qpremeanN; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
61	malvezzi	1.7	Float_t qNmin1; ///< Truncated mean (MIP) along the track using N-1 measurements before of the interaction plane
62			Float_t maxrel; ///<Energy maximum release on first Calorimeter plane (dedx of strip with maximum release)
63	mocchiut	1.1	Float_t ethr; ///< Threshold (MIP) needed to find the interaction plane
64			Bool_t multhit; ///< True if the interaction plane has been determined by multiple hit counting
65			Bool_t gap; ///< True if determining the interaction plane a big (>5 planes) gap has been found between a point and another along the track
66	malvezzi	1.7	Float_t charge_siegen1;
67			Float_t ZCalo_dedx_b; //Z from Calo using dedx in first Calorimeter plane vs. beta
68			Float_t ZCalo_maxrel_b; //Z from Calo using maximum release in first Calorimeter plane vs. beta
69			Float_t ZCalo_dedx_defl; //Z from Calo using dedx in first Calorimeter plane vs. rigidity
70			Float_t ZCalo_Nmin1_defl; //Z from Calo using truncated mean on N-1 Calorimeter planes (plane N+1 is the interaction plane) vs. rigidity
71			//
72			Float_t qNmin1_w;
73			Int_t S2;
74	mocchiut	1.1	public:
75			//
76			//
77	malvezzi	1.7	//char* version();
78
79	mocchiut	1.3	Int_t Get_interplane(){ Process(); return interplane;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
80	mocchiut	1.6	Int_t Get_N(){ return N;}; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
81			Int_t Get_UsedN(){ return UN;}; ///< Number of dE/dx measurements really used to calculate qpremeanN
82	mocchiut	1.3	Int_t Get_R(){ return R;}; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
83			Float_t Get_preq(){ Process(); return preq;}; ///< Energy release (MIP) up to the interaction plane (included)
84			Float_t Get_postq(){ Process(); return postq;}; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
85	mocchiut	1.5	Float_t Get_StdEdx1(){ Process(); return stdedx1;}; ///< Energy release (MIP) on the first Silicon detector (Y EVEN) around the strip with maximum energy release (no track information, 3 strips in total).
86	mocchiut	1.3	Float_t Get_dEdx1(){ Process(); return dedx1;}; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
87			Float_t Get_dEdx3(){ Process(); return dedx3;}; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
88			Float_t Get_qpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
89			Float_t Get_qpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
90	malvezzi	1.7	Float_t Get_qNmin1(){ Process(); return qNmin1;}; ///< Truncated mean (MIP) along the track using N-1 measurements before of the interaction plane
91			Float_t Get_maxrel(){ Process(); return maxrel;}; ///<Energy maximum release on first Calorimeter plane (dedx of strip with maximum release)
92	mocchiut	1.3	Float_t Get_ethr(){ Process(); return ethr;}; ///< Threshold (MIP) needed to find the interaction plane
93	mocchiut	1.1	Bool_t IsMulthit(){ Process(); return multhit;}; ///< True if the interaction plane has been determined by multiple hit counting
94			//
95	malvezzi	1.7	Float_t Get_charge_siegen1(){ Process(); return charge_siegen1;}; ///< charge Siegen method stdedx1 vs. beta
96			Float_t Get_ZCalo_dedx_b(){ Process(); return ZCalo_dedx_b;}; //Z from Calo using dedx (or StdEdx) in first Calorimeter plane vs. beta
97			Float_t Get_ZCalo_maxrel_b(){ Process(); return ZCalo_maxrel_b;}; //Z from Calo using maximum release in first Calorimeter plane vs. beta
98			Float_t Get_ZCalo_dedx_defl(){ Process(); return ZCalo_dedx_defl;}; //Z from Calo using dedx in first Calorimeter plane vs. rigidity
99			Float_t Get_ZCalo_Nmin1_defl(){ Process(); return ZCalo_Nmin1_defl;}; //Z from Calo using truncated mean on N-1 Calorimeter planes (plane N+1 is the interaction plane) vs. rigidity
100
101			//
102	mocchiut	1.3	void Set_N(Int_t n){ N=n;};
103			void Set_R(Int_t r){ R=r;};
104	mocchiut	1.1	//
105	pam-fi	1.10	// CaloNuclei();
106			CaloNuclei(PamLevel2 L2,const char alg);
107	mocchiut	1.1	~CaloNuclei(){ Delete(); };
108			//
109	mocchiut	1.4	void SetDebug(Bool_t d){ debug=d; };
110	mocchiut	1.5	void UseTrack(Bool_t d){ usetrack=d; };
111	mocchiut	1.8
112
113			void UsePlane18X(Bool_t use){usepl18x = use;};
114			//
115	mocchiut	1.4	//
116	mocchiut	1.1	void Clear();
117			void Clear(Option_t *option){Clear();};
118			void Delete();
119	mocchiut	1.6	void Delete(Option_t *option){Delete();};
120	mocchiut	1.1	//
121	mocchiut	1.3	void Process(); ///< Process data for track number 0
122	mocchiut	1.5	void Process(Int_t ntr); ///< Process data for track number ntr
123	mocchiut	1.1	void Print();
124	mocchiut	1.6	void Print(Option_t *option){Print();};
125	mocchiut	1.1	//
126	mocchiut	1.9	ClassDef(CaloNuclei,3);
127	mocchiut	1.1	};
128
129			#endif
130