PamelaLevel2/inc/PamLevel2.h



/**
 * \file PamLevel2.h
 * \author Elena Vannuccini
 */
#ifndef pamlevel2_h
#define pamlevel2_h

#include <TRefArray.h>

#include <TrkLevel1.h>
#include <TrkLevel2.h>
#include <TrkHough.h>
#include <CaloLevel1.h>
#include <CaloLevel2.h>
#include <ToFLevel2.h>
#include <TrigLevel2.h>
#include <S4Level2.h>
#include <NDLevel2.h>
#include <AcLevel2.h>
#include <OrbitalInfo.h>
#include <OrbitalInfo.h>
#include <RunInfo.h>
#include <GLTables.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 <iostream>
using namespace std;

/**
 * \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).
 *
 * This class inherites from TrkTrack, CaloTrkVar ecc... 
 * 
 */
class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {

 private:

    TrkTrack*   trk_track;
    CaloTrkVar* calo_track;
    ToFTrkVar*  tof_track;

 public:

    PamTrack();

    PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );

    ClassDef(PamTrack,1);
};

/**
 * \brief Class to describe PAMELA LEVEL2 data.
 *
 * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods. 
 * This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.
 */
                
class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {

 private:
 public:

    TrkLevel1*   trk1_obj;
    TrkLevel2*   trk2_obj;
    TrkHough*    trkh_obj;

    CaloLevel1*  calo1_obj;
    CaloLevel2*  calo2_obj;

    ToFLevel2*   tof_obj;
    TrigLevel2*  trig_obj;
    S4Level2*    s4_obj;
    NDLevel2*    nd_obj;
    AcLevel2*    ac_obj;
    OrbitalInfo* orb_obj;
        
    GL_RUN*      run_obj;
    
/*     TrkLevel1*   trk1_clone; */
/*     TrkLevel2*   trk2_clone; */
/*     TrkHough*    trkh_clone; */

/*     CaloLevel1*  calo1_clone; */
/*     CaloLevel2*  calo2_clone; */

/*     ToFLevel2*   tof_clone; */
/*     TrigLevel2*  trig_clone; */
/*     S4Level2*    s4_clone; */
/*     NDLevel2*    nd_clone; */
/*     AcLevel2*    ac_clone; */
/*     OrbitalInfo* orb_clone; */

/*     TTree*       T_clone; //tracker */
/*     TTree*       C_clone; //calo */
/*     TTree*       O_clone; //tof */
/*     TTree*       R_clone; //trigger */
/*     TTree*       S_clone; //s4 */
/*     TTree*       A_clone; //anti */
/*     TTree*       N_clone; //nd */
/*     TTree*       B_clone; //orbital */

    TTree*       tree_clone[8];

    TRefArray*   sorted_tracks; 
        
//    TChain*      Tout;
        

    Bool_t CAL2; //level2 branch
    Bool_t CAL1; //level1 branch

    Bool_t TRK2; //level2 branch
    Bool_t TRK1; //level1 branch
    Bool_t TRKh; //hough branch

    Bool_t TRG;
    Bool_t TOF;
    Bool_t S4;
    Bool_t ND;
    Bool_t AC;
    Bool_t ORB;
         
    PamLevel2();
    ~PamLevel2();
        
    void Clear();

/*     TrkLevel1*   GetTrkLevel1()  { return trk1_obj;  }; */
/*     TrkLevel2*   GetTrkLevel2()  { return trk2_obj;  }; */
/*     TrkHough*    GetTrkHough()   { return trkh_obj;  }; */
/*     CaloLevel1*  GetCaloLevel1() { return calo1_obj; }; */
/*     CaloLevel2*  GetCaloLevel2() { return calo2_obj; }; */
/*     ToFLevel2*   GetToFLevel2()  { return tof_obj;  }; */
/*     TrigLevel2*  GetTrigLevel2() { return trig_obj; }; */
/*     S4Level2*    GetS4Level2()   { return s4_obj;   }; */
/*     NDLevel2*    GetNDLevel2()   { return nd_obj;   }; */
/*     AcLevel2*    GetAcLevel2()   { return ac_obj;   }; */
/*     OrbitalInfo* GetOrbitalInfo(){ return orb_obj;  }; */
/*     TrkLevel1*   GetTrkLevel1()  { return ; }; */
/*     TrkLevel2*   GetTrkLevel2()  { return TrkLevel2::GetTrkLevel2(); }; */
/*     TrkHough*    GetTrkHough()   { return TrkHough::GetTrkHough(); }; */
/*     CaloLevel1*  GetCaloLevel1() { return CaloLevel1::GetCaloLevel1(); }; */
/*     CaloLevel2*  GetCaloLevel2() { return CaloLevel2::GetCaloLevel2(); }; */
/*     ToFLevel2*   GetToFLevel2()  { return ToFLevel2::GetToFLevel2(); }; */
/*     TrigLevel2*  GetTrigLevel2() { return TrigLevel2::GetTrigLevel2(); }; */
/*     S4Level2*    GetS4Level2()   { return S4Level2::GetS4Level2(); }; */
/*     NDLevel2*    GetNDLevel2()   { return NDLevel2::GetNDLevel2(); }; */
/*     AcLevel2*    GetAcLevel2()   { return AcLevel2::GetAcLevel2(); }; */
/*     OrbitalInfo* GetOrbitalInfo(){ return OrbitalInfo::GetOrbitalInfo(); }; */
    GL_RUN*      GetRunInfo()    { return run_obj; };
        
//  =========================================
//  methods needed to read/write Level2 trees
//  =========================================
    void    SetWhichTrees(TString);
    TList*  GetListOfLevel2Files(TString, TString);
    Bool_t  CheckLevel2File(TString);
    TTree*  GetPamTree(TFile*,TString);
    TTree*  GetPamTree(TFile* f){return GetPamTree(f,"+ALL");};
    TChain* GetPamTree(TList*,TString);
    TTree*  GetRunTree(TFile*);
    TChain* GetRunTree(TList*);

//    TTree*  GetNewPamTree();
//    void    FillNewPamTree(TTree*);
//    void    FillNewPamTree();
    void    CreateCloneTrees(TChain *fChain);
    void    FillCloneTrees();
    TTree  *GetCloneTree(TString);
    void    WriteCloneTrees();

    void    SetBranchAddress(TTree *);
    void    SetBranchAddress(TChain *);

//  =====================================
//  methods needed to set brach addresses
//  =====================================
/*     void*       GetPointerToTrk(); */
/*     void*       GetPointerToTrk(int s); */
/*     void*       GetPointerToCalo(); */
/*     void*       GetPointerToCalo(int s); */
/*     void*       GetPointerToToF(); */
/*     void*       GetPointerToTrig(); */
/*     void*       GetPointerToS4(); */
/*     void*       GetPointerToND(); */
/*     void*       GetPointerToAc(); */
/*     void*       GetPointerToOrb(); */
/*     void*       GetPointerToRun(); */
    void*       GetPointerTo(const char*);

//  ==============================================
//  methods to retrieve subdetectors stored tracks
//  ==============================================
    TrkTrack*   GetTrkStoredTrack(int seqno){return TrkLevel2::GetStoredTrack(seqno); };
    CaloTrkVar* GetCaloStoredTrack(int);
    ToFTrkVar*  GetToFStoredTrack(int);
    //
    PamTrack*   GetPamTrackAlong(TrkTrack*);
//      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
//      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
        
//  =================================
//  methods to retrieve pamela tracks
//  =================================
    PamTrack*     GetStoredTrack(Int_t);
    //
    void          SortTracks(TString);          //assign value to variable sorted_tracks
    TRefArray*    GetTracks();                  //return sorted_tracks
    PamTrack*     GetTrack(int);
    PamTrack*     GetTrackImage(int);
        
//  =================================
//  other methods
//  =================================
    Int_t ntrk(){ return TrkLevel2::ntrk(); };


    ClassDef(PamLevel2,1);
};

//===============================================================================================
//
//
//
//
//===============================================================================================


#endif

1	pam-fi	1.1
2
3			/**
4			* \file PamLevel2.h
5			* \author Elena Vannuccini
6			*/
7			#ifndef pamlevel2_h
8			#define pamlevel2_h
9
10	pam-fi	1.2	#include <TRefArray.h>
11	pam-fi	1.1
12	pam-fi	1.4	#include <TrkLevel1.h>
13	pam-fi	1.1	#include <TrkLevel2.h>
14	pam-fi	1.6	#include <TrkHough.h>
15			#include <CaloLevel1.h>
16	pam-fi	1.1	#include <CaloLevel2.h>
17			#include <ToFLevel2.h>
18			#include <TrigLevel2.h>
19			#include <S4Level2.h>
20			#include <NDLevel2.h>
21			#include <AcLevel2.h>
22	pam-fi	1.2	#include <OrbitalInfo.h>
23	pam-fi	1.4	#include <OrbitalInfo.h>
24			#include <RunInfo.h>
25			#include <GLTables.h>
26
27	pam-fi	1.1
28			#include <TTree.h>
29	pam-fi	1.7	#include <TFriendElement.h>
30	pam-fi	1.2	#include <TChain.h>
31	pam-fi	1.1	#include <TFile.h>
32	pam-fi	1.2	#include <TList.h>
33			#include <TKey.h>
34			#include <TSystemFile.h>
35			#include <TSystemDirectory.h>
36	pam-fi	1.1
37			#include <iostream>
38			using namespace std;
39
40			/**
41			* \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).
42			*
43			* This class inherites from TrkTrack, CaloTrkVar ecc...
44			*
45			*/
46			class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
47
48			private:
49
50			TrkTrack* trk_track;
51			CaloTrkVar* calo_track;
52			ToFTrkVar* tof_track;
53
54			public:
55
56			PamTrack();
57
58			PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );
59
60			ClassDef(PamTrack,1);
61			};
62
63			/**
64			* \brief Class to describe PAMELA LEVEL2 data.
65			*
66			* This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.
67			* This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.
68			*/
69
70	pam-fi	1.7	class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {
71	pam-fi	1.1
72			private:
73	pam-fi	1.7	public:
74	pam-fi	1.1
75	pam-fi	1.6	TrkLevel1* trk1_obj;
76			TrkLevel2* trk2_obj;
77			TrkHough* trkh_obj;
78
79			CaloLevel1* calo1_obj;
80			CaloLevel2* calo2_obj;
81
82	pam-fi	1.2	ToFLevel2* tof_obj;
83			TrigLevel2* trig_obj;
84			S4Level2* s4_obj;
85			NDLevel2* nd_obj;
86			AcLevel2* ac_obj;
87	pam-fi	1.6	OrbitalInfo* orb_obj;
88	pam-fi	1.4
89	pam-fi	1.6	GL_RUN* run_obj;
90	pam-fi	1.1
91	pam-fi	1.7	/* TrkLevel1* trk1_clone; */
92			/* TrkLevel2* trk2_clone; */
93			/* TrkHough* trkh_clone; */
94
95			/* CaloLevel1* calo1_clone; */
96			/* CaloLevel2* calo2_clone; */
97
98			/* ToFLevel2* tof_clone; */
99			/* TrigLevel2* trig_clone; */
100			/* S4Level2* s4_clone; */
101			/* NDLevel2* nd_clone; */
102			/* AcLevel2* ac_clone; */
103			/* OrbitalInfo* orb_clone; */
104
105			/* TTree* T_clone; //tracker */
106			/* TTree* C_clone; //calo */
107			/* TTree* O_clone; //tof */
108			/* TTree* R_clone; //trigger */
109			/* TTree* S_clone; //s4 */
110			/* TTree* A_clone; //anti */
111			/* TTree* N_clone; //nd */
112			/* TTree* B_clone; //orbital */
113
114			TTree* tree_clone[8];
115
116	pam-fi	1.6	TRefArray* sorted_tracks;
117	pam-fi	1.4
118	pam-fi	1.6	// TChain* Tout;
119	pam-fi	1.2
120	pam-fi	1.1
121	pam-fi	1.6	Bool_t CAL2; //level2 branch
122			Bool_t CAL1; //level1 branch
123
124			Bool_t TRK2; //level2 branch
125			Bool_t TRK1; //level1 branch
126			Bool_t TRKh; //hough branch
127
128			Bool_t TRG;
129			Bool_t TOF;
130			Bool_t S4;
131			Bool_t ND;
132			Bool_t AC;
133			Bool_t ORB;
134
135	pam-fi	1.1	PamLevel2();
136	pam-fi	1.6	~PamLevel2();
137	pam-fi	1.2
138	pam-fi	1.6	void Clear();
139	pam-fi	1.2
140	pam-fi	1.7	/* TrkLevel1* GetTrkLevel1() { return trk1_obj; }; */
141			/* TrkLevel2* GetTrkLevel2() { return trk2_obj; }; */
142			/* TrkHough* GetTrkHough() { return trkh_obj; }; */
143			/* CaloLevel1* GetCaloLevel1() { return calo1_obj; }; */
144			/* CaloLevel2* GetCaloLevel2() { return calo2_obj; }; */
145			/* ToFLevel2* GetToFLevel2() { return tof_obj; }; */
146			/* TrigLevel2* GetTrigLevel2() { return trig_obj; }; */
147			/* S4Level2* GetS4Level2() { return s4_obj; }; */
148			/* NDLevel2* GetNDLevel2() { return nd_obj; }; */
149			/* AcLevel2* GetAcLevel2() { return ac_obj; }; */
150			/* OrbitalInfo* GetOrbitalInfo(){ return orb_obj; }; */
151			/* TrkLevel1* GetTrkLevel1() { return ; }; */
152			/* TrkLevel2* GetTrkLevel2() { return TrkLevel2::GetTrkLevel2(); }; */
153			/* TrkHough* GetTrkHough() { return TrkHough::GetTrkHough(); }; */
154			/* CaloLevel1* GetCaloLevel1() { return CaloLevel1::GetCaloLevel1(); }; */
155			/* CaloLevel2* GetCaloLevel2() { return CaloLevel2::GetCaloLevel2(); }; */
156			/* ToFLevel2* GetToFLevel2() { return ToFLevel2::GetToFLevel2(); }; */
157			/* TrigLevel2* GetTrigLevel2() { return TrigLevel2::GetTrigLevel2(); }; */
158			/* S4Level2* GetS4Level2() { return S4Level2::GetS4Level2(); }; */
159			/* NDLevel2* GetNDLevel2() { return NDLevel2::GetNDLevel2(); }; */
160			/* AcLevel2* GetAcLevel2() { return AcLevel2::GetAcLevel2(); }; */
161			/* OrbitalInfo* GetOrbitalInfo(){ return OrbitalInfo::GetOrbitalInfo(); }; */
162			GL_RUN* GetRunInfo() { return run_obj; };
163
164			// =========================================
165			// methods needed to read/write Level2 trees
166			// =========================================
167	pam-fi	1.6	void SetWhichTrees(TString);
168			TList* GetListOfLevel2Files(TString, TString);
169			Bool_t CheckLevel2File(TString);
170			TTree* GetPamTree(TFile*,TString);
171			TTree* GetPamTree(TFile* f){return GetPamTree(f,"+ALL");};
172			TChain* GetPamTree(TList*,TString);
173			TTree* GetRunTree(TFile*);
174			TChain* GetRunTree(TList*);
175	pam-fi	1.7
176			// TTree* GetNewPamTree();
177			// void FillNewPamTree(TTree*);
178			// void FillNewPamTree();
179			void CreateCloneTrees(TChain *fChain);
180			void FillCloneTrees();
181			TTree *GetCloneTree(TString);
182			void WriteCloneTrees();
183
184	pam-fi	1.6	void SetBranchAddress(TTree *);
185			void SetBranchAddress(TChain *);
186
187	pam-fi	1.2	// =====================================
188			// methods needed to set brach addresses
189			// =====================================
190	pam-fi	1.7	/* void* GetPointerToTrk(); */
191			/* void* GetPointerToTrk(int s); */
192			/* void* GetPointerToCalo(); */
193			/* void* GetPointerToCalo(int s); */
194			/* void* GetPointerToToF(); */
195			/* void* GetPointerToTrig(); */
196			/* void* GetPointerToS4(); */
197			/* void* GetPointerToND(); */
198			/* void* GetPointerToAc(); */
199			/* void* GetPointerToOrb(); */
200			/* void* GetPointerToRun(); */
201	pam-fi	1.6	void* GetPointerTo(const char*);
202	pam-fi	1.3
203			// ==============================================
204			// methods to retrieve subdetectors stored tracks
205			// ==============================================
206	pam-fi	1.6	TrkTrack* GetTrkStoredTrack(int seqno){return TrkLevel2::GetStoredTrack(seqno); };
207			CaloTrkVar* GetCaloStoredTrack(int);
208			ToFTrkVar* GetToFStoredTrack(int);
209			//
210			PamTrack* GetPamTrackAlong(TrkTrack*);
211	pam-fi	1.3	// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve
212			// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve
213
214	pam-fi	1.1	// =================================
215			// methods to retrieve pamela tracks
216			// =================================
217	pam-fi	1.3	PamTrack* GetStoredTrack(Int_t);
218	pam-fi	1.6	//
219			void SortTracks(TString); //assign value to variable sorted_tracks
220			TRefArray* GetTracks(); //return sorted_tracks
221			PamTrack* GetTrack(int);
222	pam-fi	1.1	PamTrack* GetTrackImage(int);
223	pam-fi	1.2
224	pam-fi	1.6	// =================================
225			// other methods
226			// =================================
227			Int_t ntrk(){ return TrkLevel2::ntrk(); };
228
229
230	pam-fi	1.1	ClassDef(PamLevel2,1);
231			};
232
233			//===============================================================================================
234			//
235			//
236			//
237			//
238			//===============================================================================================
239
240
241			#endif
242