CaloPreSampler/inc/CaloPreSampler.h

/**
 * \file CaloPreSampler.h
 * \author Emiliano Mocchiutti (2007/07/18)
 */
#ifndef calopresampler_h
#define calopresampler_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 presampler analysis
 */
class CaloPreSampler : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    //
    // needed to avoid reprocessing the same event over and over to obtain the variables;
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    Bool_t ssel;
    //
    Bool_t usepl18x;
    //
    Bool_t simulation; ///< True when using simulated data, false by default;
    //
    CaloLevel0 *event; ///< CaloLevel0 event needed to process data from level0/1 to level2
    CaloStrip *cstrip; ///< CaloStrip object needed to easily recover parameter files
    CaloLevel2 *pcalo; ///< CaloLevel2 object which contains variables calculated with a reduced calorimeter.
    //   CaloLevel1 *c1; ///< CaloLevel1 object which contains variables calculated with a reduced calorimeter.
    //
    Int_t N; ///< Number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used;
    Int_t NC; ///< Number of W planes to be used as calorimeter.
    Bool_t sel; ///< Selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true;
    Bool_t cont; ///< Contamination mode: planes from N to 22 are used;
    Bool_t emulate18; ///< If true emulates the non-working plane 18;
    Bool_t nox; ///< If true do not use X-view
    Bool_t noy; ///< If true do not use Y-view
    Bool_t forcecalo;
    Int_t mask[2][22];
    //
    Float_t rigdefault; ///< Rigidity to be used if no track is available or in standalone mode
    Bool_t withtrk; ///< Use or not the tracker
    Int_t forcefitmode;
    //

 public:
    //
    //
    // CaloLevel1* GetCaloLevel1(){ Process(); return c1;}; ///< Retrieve CaloLevel1 pointer which contains variables calculated with a reduced calorimeter.
    CaloLevel2* GetCaloLevel2(){ Process(); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
    CaloLevel2* GetLevel2Pointer(){ printf(" OBSOLETE, use GetCaloLevel2Pointer()!\n"); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
    CaloLevel2* GetCaloLevel2Pointer(){ return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
    //
    CaloTrkVar* GetCaloTrack(Int_t t);
    //
    void SetNoWpreSampler(Int_t n); ///< Set the number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used. Default: N = 4;
    void SetNoWcalo(Int_t n); ///< Set the number of W planes to be used as calorimeter. Default: NC = 22-N (N presampler);
    void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo);
    void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT);
    void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used.
    void EmulatePlane18(Bool_t emu){emulate18 = emu;} ///< Emulates the non-working plane 18; Default = true;
    void Simulation(Bool_t sim){simulation = sim;}; ///< Set simulated data flag. Default FALSE;
    //
    CaloPreSampler();
    CaloPreSampler(PamLevel2 *L2);
    ~CaloPreSampler(){ Delete(); };
    //
    Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler. 
    Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter.
    void SetDebug(Bool_t d){ debug=d; };
    void UseTracker(Bool_t ch){ withtrk = ch; };
    void SetDefaultRig(Float_t rig){ rigdefault = rig; };
    void SetForceFitMode(Int_t mode){ forcefitmode = mode; };
    void MaskX(Bool_t mskx){ nox = mskx; };
    void MaskY(Bool_t msky){ noy = msky; };
    void Mask(Int_t view, Int_t plane){ mask[view][plane] = 1;};
    void ForceCaloFit(){forcecalo=true;};
    void ForceNextProc(){atime=0;};
    //
    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 
    void Print();
    void Print(Option_t *option){Print();};
    //
    CaloTrkVar* AddCaloTrkVar(float *al,int trktag);//ELENA

    ClassDef(CaloPreSampler,3);
};

#endif

1	mocchiut	1.1	/**
2			* \file CaloPreSampler.h
3			* \author Emiliano Mocchiutti (2007/07/18)
4			*/
5			#ifndef calopresampler_h
6			#define calopresampler_h
7
8			#include <PamLevel2.h>
9
10			#include <TTree.h>
11			#include <TFriendElement.h>
12			#include <TChain.h>
13			#include <TFile.h>
14			#include <TList.h>
15			#include <TKey.h>
16			#include <TSystemFile.h>
17			#include <TSystemDirectory.h>
18			#include <TSQLServer.h>
19
20			#include <iostream>
21
22			using namespace std;
23
24			/**
25			*
26			* Class to store and calculate variables useful for presampler analysis
27			*/
28			class CaloPreSampler : public TObject {
29
30			private:
31			//
32			PamLevel2 *L2;
33			Bool_t debug;
34			//
35			// needed to avoid reprocessing the same event over and over to obtain the variables;
36			//
37			UInt_t OBT;
38			UInt_t PKT;
39			UInt_t atime;
40			Bool_t ssel;
41			//
42	mocchiut	1.15	Bool_t usepl18x;
43			//
44	mocchiut	1.1	Bool_t simulation; ///< True when using simulated data, false by default;
45			//
46			CaloLevel0 *event; ///< CaloLevel0 event needed to process data from level0/1 to level2
47			CaloStrip *cstrip; ///< CaloStrip object needed to easily recover parameter files
48			CaloLevel2 *pcalo; ///< CaloLevel2 object which contains variables calculated with a reduced calorimeter.
49	mocchiut	1.2	// CaloLevel1 *c1; ///< CaloLevel1 object which contains variables calculated with a reduced calorimeter.
50	mocchiut	1.1	//
51			Int_t N; ///< Number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used;
52	mocchiut	1.5	Int_t NC; ///< Number of W planes to be used as calorimeter.
53	mocchiut	1.1	Bool_t sel; ///< Selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true;
54			Bool_t cont; ///< Contamination mode: planes from N to 22 are used;
55			Bool_t emulate18; ///< If true emulates the non-working plane 18;
56	mocchiut	1.4	Bool_t nox; ///< If true do not use X-view
57			Bool_t noy; ///< If true do not use Y-view
58	mocchiut	1.6	Bool_t forcecalo;
59	mocchiut	1.7	Int_t mask[2][22];
60	mocchiut	1.1	//
61	mocchiut	1.3	Float_t rigdefault; ///< Rigidity to be used if no track is available or in standalone mode
62			Bool_t withtrk; ///< Use or not the tracker
63	mocchiut	1.10	Int_t forcefitmode;
64	mocchiut	1.3	//
65	mocchiut	1.1
66			public:
67			//
68			//
69	mocchiut	1.2	// CaloLevel1* GetCaloLevel1(){ Process(); return c1;}; ///< Retrieve CaloLevel1 pointer which contains variables calculated with a reduced calorimeter.
70	mocchiut	1.1	CaloLevel2* GetCaloLevel2(){ Process(); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
71	mocchiut	1.12	CaloLevel2* GetLevel2Pointer(){ printf(" OBSOLETE, use GetCaloLevel2Pointer()!\n"); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
72			CaloLevel2* GetCaloLevel2Pointer(){ return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
73	mocchiut	1.1	//
74	mocchiut	1.13	CaloTrkVar* GetCaloTrack(Int_t t);
75			//
76	mocchiut	1.5	void SetNoWpreSampler(Int_t n); ///< Set the number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used. Default: N = 4;
77			void SetNoWcalo(Int_t n); ///< Set the number of W planes to be used as calorimeter. Default: NC = 22-N (N presampler);
78	mocchiut	1.8	void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo);
79	mocchiut	1.1	void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT);
80			void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used.
81			void EmulatePlane18(Bool_t emu){emulate18 = emu;} ///< Emulates the non-working plane 18; Default = true;
82			void Simulation(Bool_t sim){simulation = sim;}; ///< Set simulated data flag. Default FALSE;
83			//
84			CaloPreSampler();
85			CaloPreSampler(PamLevel2 *L2);
86			~CaloPreSampler(){ Delete(); };
87			//
88	mocchiut	1.5	Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler.
89			Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter.
90	mocchiut	1.1	void SetDebug(Bool_t d){ debug=d; };
91	mocchiut	1.3	void UseTracker(Bool_t ch){ withtrk = ch; };
92			void SetDefaultRig(Float_t rig){ rigdefault = rig; };
93	mocchiut	1.10	void SetForceFitMode(Int_t mode){ forcefitmode = mode; };
94	mocchiut	1.4	void MaskX(Bool_t mskx){ nox = mskx; };
95			void MaskY(Bool_t msky){ noy = msky; };
96	mocchiut	1.7	void Mask(Int_t view, Int_t plane){ mask[view][plane] = 1;};
97	mocchiut	1.6	void ForceCaloFit(){forcecalo=true;};
98	mocchiut	1.14	void ForceNextProc(){atime=0;};
99	mocchiut	1.1	//
100	mocchiut	1.15	void UsePlane18X(Bool_t use){usepl18x = use;};
101			//
102	mocchiut	1.1	void Clear();
103			void Clear(Option_t *option){Clear();};
104			void Delete();
105	mocchiut	1.4	void Delete(Option_t *option){Delete();};
106	mocchiut	1.1	//
107			void Process(); ///< Process data
108			void Print();
109	mocchiut	1.4	void Print(Option_t *option){Print();};
110	mocchiut	1.1	//
111	pam-fi	1.9	CaloTrkVar* AddCaloTrkVar(float *al,int trktag);//ELENA
112
113	mocchiut	1.13	ClassDef(CaloPreSampler,3);
114	mocchiut	1.1	};
115
116			#endif
117