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
2
3	/**
4	* \file PamLevel2.h
5	* \author Elena Vannuccini
6	*/
7	#ifndef pamlevel2_h
8	#define pamlevel2_h
9
10	#include <TRefArray.h>
11
12	#include <TrkLevel1.h>
13	#include <TrkLevel2.h>
14	#include <TrkHough.h>
15	#include <CaloLevel1.h>
16	#include <CaloLevel2.h>
17	#include <ToFLevel2.h>
18	#include <TrigLevel2.h>
19	#include <S4Level2.h>
20	#include <NDLevel2.h>
21	#include <AcLevel2.h>
22	#include <OrbitalInfo.h>
23	#include <OrbitalInfo.h>
24	#include <RunInfo.h>
25	#include <GLTables.h>
26
27
28	#include <TTree.h>
29	#include <TFriendElement.h>
30	#include <TChain.h>
31	#include <TFile.h>
32	#include <TList.h>
33	#include <TKey.h>
34	#include <TSystemFile.h>
35	#include <TSystemDirectory.h>
36
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	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
72	private:
73	public:
74
75	TrkLevel1* trk1_obj;
76	TrkLevel2* trk2_obj;
77	TrkHough* trkh_obj;
78
79	CaloLevel1* calo1_obj;
80	CaloLevel2* calo2_obj;
81
82	ToFLevel2* tof_obj;
83	TrigLevel2* trig_obj;
84	S4Level2* s4_obj;
85	NDLevel2* nd_obj;
86	AcLevel2* ac_obj;
87	OrbitalInfo* orb_obj;
88
89	GL_RUN* run_obj;
90
91	/* 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	TRefArray* sorted_tracks;
117
118	// TChain* Tout;
119
120
121	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	PamLevel2();
136	~PamLevel2();
137
138	void Clear();
139
140	/* 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	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
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	void SetBranchAddress(TTree *);
185	void SetBranchAddress(TChain *);
186
187	// =====================================
188	// methods needed to set brach addresses
189	// =====================================
190	/* 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	void* GetPointerTo(const char*);
202
203	// ==============================================
204	// methods to retrieve subdetectors stored tracks
205	// ==============================================
206	TrkTrack* GetTrkStoredTrack(int seqno){return TrkLevel2::GetStoredTrack(seqno); };
207	CaloTrkVar* GetCaloStoredTrack(int);
208	ToFTrkVar* GetToFStoredTrack(int);
209	//
210	PamTrack* GetPamTrackAlong(TrkTrack*);
211	// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve
212	// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve
213
214	// =================================
215	// methods to retrieve pamela tracks
216	// =================================
217	PamTrack* GetStoredTrack(Int_t);
218	//
219	void SortTracks(TString); //assign value to variable sorted_tracks
220	TRefArray* GetTracks(); //return sorted_tracks
221	PamTrack* GetTrack(int);
222	PamTrack* GetTrackImage(int);
223
224	// =================================
225	// other methods
226	// =================================
227	Int_t ntrk(){ return TrkLevel2::ntrk(); };
228
229
230	ClassDef(PamLevel2,1);
231	};
232
233	//===============================================================================================
234	//
235	//
236	//
237	//
238	//===============================================================================================
239
240
241	#endif
242