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;
    //
    // 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 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
    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 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 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
    //

 public:
    //
    //
    Int_t GetInterplane(){ Process(); return interplane;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Int_t GetN(){ return N;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Int_t GetR(){ return R;}; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3 
    Float_t GetPreq(){ Process(); return preq;}; ///< Energy release (MIP) up to the interaction plane (included)
    Float_t GetPostq(){ Process(); return postq;}; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
    Float_t GetDedx1(){ Process(); return dedx1;}; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
    Float_t GetDedx3(){ Process(); return dedx3;}; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
    Float_t GetQpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
    Float_t GetQpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
    Float_t GetEthr(){ 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
    //
    void SetN(Int_t n){ N=n;};
    void SetR(Int_t r){ R=r;};
    //
    CaloNuclei();
    CaloNuclei(PamLevel2 *L2);
    ~CaloNuclei(){ Delete(); };
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    //
    void Process();
    void Print();
    //
    ClassDef(CaloNuclei,1);
};

#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			//
36			// needed to avoid reprocessing the same event over and over to obtain the variables
37			//
38			UInt_t OBT;
39			UInt_t PKT;
40			UInt_t atime;
41			//
42			Int_t interplane; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
43			Int_t N; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
44	mocchiut	1.2	Int_t R; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
45	mocchiut	1.1	Float_t preq; ///< Energy release (MIP) up to the interaction plane (included)
46			Float_t postq; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
47			Float_t dedx1; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
48			Float_t dedx3; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
49			Float_t qpremean; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
50			Float_t qpremeanN; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
51			Float_t ethr; ///< Threshold (MIP) needed to find the interaction plane
52			Bool_t multhit; ///< True if the interaction plane has been determined by multiple hit counting
53			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
54			//
55
56			public:
57			//
58			//
59			Int_t GetInterplane(){ Process(); return interplane;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
60			Int_t GetN(){ return N;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
61	mocchiut	1.2	Int_t GetR(){ return R;}; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
62	mocchiut	1.1	Float_t GetPreq(){ Process(); return preq;}; ///< Energy release (MIP) up to the interaction plane (included)
63			Float_t GetPostq(){ Process(); return postq;}; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
64			Float_t GetDedx1(){ Process(); return dedx1;}; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
65			Float_t GetDedx3(){ Process(); return dedx3;}; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
66	mocchiut	1.2	Float_t GetQpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
67			Float_t GetQpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
68	mocchiut	1.1	Float_t GetEthr(){ Process(); return ethr;}; ///< Threshold (MIP) needed to find the interaction plane
69			Bool_t IsMulthit(){ Process(); return multhit;}; ///< True if the interaction plane has been determined by multiple hit counting
70			//
71			void SetN(Int_t n){ N=n;};
72	mocchiut	1.2	void SetR(Int_t r){ R=r;};
73	mocchiut	1.1	//
74			CaloNuclei();
75			CaloNuclei(PamLevel2 *L2);
76			~CaloNuclei(){ Delete(); };
77			//
78			void Clear();
79			void Clear(Option_t *option){Clear();};
80			void Delete();
81			//
82			void Process();
83			void Print();
84			//
85			ClassDef(CaloNuclei,1);
86			};
87
88			#endif
89