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 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 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			Bool_t simulation; ///< True when using simulated data, false by default;
43			//
44			CaloLevel0 *event; ///< CaloLevel0 event needed to process data from level0/1 to level2
45			CaloStrip *cstrip; ///< CaloStrip object needed to easily recover parameter files
46			CaloLevel2 *pcalo; ///< CaloLevel2 object which contains variables calculated with a reduced calorimeter.
47	mocchiut	1.2	// CaloLevel1 *c1; ///< CaloLevel1 object which contains variables calculated with a reduced calorimeter.
48	mocchiut	1.1	//
49			Int_t N; ///< Number of W planes to be used as presampler. NOTICE: none the silicon detectors attached to these planes will be used;
50	mocchiut	1.5	Int_t NC; ///< Number of W planes to be used as calorimeter.
51	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;
52			Bool_t cont; ///< Contamination mode: planes from N to 22 are used;
53			Bool_t emulate18; ///< If true emulates the non-working plane 18;
54	mocchiut	1.4	Bool_t nox; ///< If true do not use X-view
55			Bool_t noy; ///< If true do not use Y-view
56	mocchiut	1.6	Bool_t forcecalo;
57	mocchiut	1.7	Int_t mask[2][22];
58	mocchiut	1.1	//
59	mocchiut	1.3	Float_t rigdefault; ///< Rigidity to be used if no track is available or in standalone mode
60			Bool_t withtrk; ///< Use or not the tracker
61	mocchiut	1.10	Int_t forcefitmode;
62	mocchiut	1.3	//
63	mocchiut	1.1
64			public:
65			//
66			//
67	mocchiut	1.2	// CaloLevel1* GetCaloLevel1(){ Process(); return c1;}; ///< Retrieve CaloLevel1 pointer which contains variables calculated with a reduced calorimeter.
68	mocchiut	1.1	CaloLevel2* GetCaloLevel2(){ Process(); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
69	mocchiut	1.12	CaloLevel2* GetLevel2Pointer(){ printf(" OBSOLETE, use GetCaloLevel2Pointer()!\n"); return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
70			CaloLevel2* GetCaloLevel2Pointer(){ return pcalo;}; ///< Retrieve CaloLevel2 pointer which contains variables calculated with a reduced calorimeter.
71	mocchiut	1.1	//
72	mocchiut	1.13	CaloTrkVar* GetCaloTrack(Int_t t);
73			//
74	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;
75			void SetNoWcalo(Int_t n); ///< Set the number of W planes to be used as calorimeter. Default: NC = 22-N (N presampler);
76	mocchiut	1.8	void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo);
77	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);
78			void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used.
79			void EmulatePlane18(Bool_t emu){emulate18 = emu;} ///< Emulates the non-working plane 18; Default = true;
80			void Simulation(Bool_t sim){simulation = sim;}; ///< Set simulated data flag. Default FALSE;
81			//
82			CaloPreSampler();
83			CaloPreSampler(PamLevel2 *L2);
84			~CaloPreSampler(){ Delete(); };
85			//
86	mocchiut	1.5	Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler.
87			Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter.
88	mocchiut	1.1	void SetDebug(Bool_t d){ debug=d; };
89	mocchiut	1.3	void UseTracker(Bool_t ch){ withtrk = ch; };
90			void SetDefaultRig(Float_t rig){ rigdefault = rig; };
91	mocchiut	1.10	void SetForceFitMode(Int_t mode){ forcefitmode = mode; };
92	mocchiut	1.4	void MaskX(Bool_t mskx){ nox = mskx; };
93			void MaskY(Bool_t msky){ noy = msky; };
94	mocchiut	1.7	void Mask(Int_t view, Int_t plane){ mask[view][plane] = 1;};
95	mocchiut	1.6	void ForceCaloFit(){forcecalo=true;};
96	mocchiut	1.14	void ForceNextProc(){atime=0;};
97	mocchiut	1.1	//
98			void Clear();
99			void Clear(Option_t *option){Clear();};
100			void Delete();
101	mocchiut	1.4	void Delete(Option_t *option){Delete();};
102	mocchiut	1.1	//
103			void Process(); ///< Process data
104			void Print();
105	mocchiut	1.4	void Print(Option_t *option){Print();};
106	mocchiut	1.1	//
107	pam-fi	1.9	CaloTrkVar* AddCaloTrkVar(float *al,int trktag);//ELENA
108
109	mocchiut	1.13	ClassDef(CaloPreSampler,3);
110	mocchiut	1.1	};
111
112			#endif
113