24 |
#include <NDLevel2.h> |
#include <NDLevel2.h> |
25 |
#include <AcLevel2.h> |
#include <AcLevel2.h> |
26 |
#include <OrbitalInfo.h> |
#include <OrbitalInfo.h> |
27 |
#include <InclinationInfo.h> // EMILIANO, now this library is included in libDarthVader.so |
#include <InclinationInfo.h> // EMILIANO, now this library is included in libDarthVader.so |
28 |
// and the header installed in $PAM_INC. If you get the |
// and the header installed in $PAM_INC. If you get the |
29 |
// error "Cannot find InclinationInfo.h" update your DV installation |
// error "Cannot find InclinationInfo.h" update your DV installation |
30 |
#include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so |
#include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so |
31 |
// and the header installed in $PAM_INC. If you get the |
// and the header installed in $PAM_INC. If you get the |
32 |
// error "Cannot find InclinationInfo.h" update your DV installation |
// error "Cannot find InclinationInfo.h" update your DV installation |
33 |
#include <RunInfo.h> |
#include <RunInfo.h> |
34 |
#include <GLTables.h> |
#include <GLTables.h> |
61 |
///////////////////////////////////////////////////////////////////////////// |
///////////////////////////////////////////////////////////////////////////// |
62 |
/** |
/** |
63 |
* \brief Class to describe GP output (to be moved to PamelaDigitizer) |
* \brief Class to describe GP output (to be moved to PamelaDigitizer) |
64 |
* |
* |
65 |
*/ |
*/ |
66 |
|
|
67 |
class GPamela : public TObject { |
class GPamela : public TObject { |
217 |
/** |
/** |
218 |
* \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF). |
* \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF). |
219 |
* |
* |
220 |
* This class inherites from TrkTrack, CaloTrkVar ecc... |
* This class inherites from TrkTrack, CaloTrkVar ecc... |
221 |
* |
* |
222 |
*/ |
*/ |
223 |
//class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar { |
//class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar { |
224 |
class PamTrack : public TObject { |
class PamTrack : public TObject { |
264 |
/** |
/** |
265 |
* \brief Class to describe PAMELA LEVEL2 data. |
* \brief Class to describe PAMELA LEVEL2 data. |
266 |
* |
* |
267 |
* This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods. |
* This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods. |
268 |
* This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors. |
* This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors. |
269 |
*/ |
*/ |
270 |
|
|
271 |
//class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo { |
//class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo { |
272 |
class PamLevel2 : public TObject { |
class PamLevel2 : public TObject { |
273 |
|
|
286 |
TChain* sel_tree; |
TChain* sel_tree; |
287 |
TTree* sel_tree_clone; |
TTree* sel_tree_clone; |
288 |
int *irunoffset; |
int *irunoffset; |
289 |
|
GL_TIMESYNC *gltsync; // Emiliano |
290 |
|
Bool_t fUpdateRunInfo; // Emiliano |
291 |
|
|
292 |
// ------------------------------------------- |
// ------------------------------------------- |
293 |
// variables to store total live and dead time |
// variables to store total live and dead time |
294 |
// ------------------------------------------- |
// ------------------------------------------- |
295 |
ULong64_t totdltime[3];// 0=live 1=dead 2=n.events |
ULong64_t totdltime[3];// 0=live 1=dead 2=n.events |
296 |
|
|
297 |
// --------------------------------- |
// --------------------------------- |
298 |
// variables to retrieve level0 info |
// variables to retrieve level0 info |
299 |
// --------------------------------- |
// --------------------------------- |
300 |
TFile* l0_file; |
TFile* l0_file; |
301 |
TTree* l0_tree; |
TTree* l0_tree; |
302 |
Int_t iroot; |
Int_t iroot; |
303 |
Int_t prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time |
Int_t prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time |
304 |
Int_t maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10 |
Int_t maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10 |
305 |
// -------------------------------------------- |
// -------------------------------------------- |
306 |
// variables to store Level0 info (if required) |
// variables to store Level0 info (if required) |
307 |
// -------------------------------------------- |
// -------------------------------------------- |
324 |
TTree* pam_tree_clone[NCLONES]; |
TTree* pam_tree_clone[NCLONES]; |
325 |
|
|
326 |
// -------------------------------------------- |
// -------------------------------------------- |
327 |
// variables to store Level1-2 info |
// variables to store Level1-2 info |
328 |
// -------------------------------------------- |
// -------------------------------------------- |
329 |
TrkLevel1* trk1_obj; |
TrkLevel1* trk1_obj; |
330 |
TrkLevel2* trk2_obj; |
TrkLevel2* trk2_obj; |
350 |
// ----------------------------------- |
// ----------------------------------- |
351 |
TClonesArray *tsorted; |
TClonesArray *tsorted; |
352 |
TClonesArray *timage; |
TClonesArray *timage; |
353 |
PamTrack *track; //store the last required |
PamTrack *track; //store the last required |
354 |
TString howtosort; |
TString howtosort; |
355 |
Float_t sortthr; |
Float_t sortthr; |
356 |
|
|
357 |
|
Bool_t issorted; // true if the current event has already been sorted with the current sort criteria |
358 |
|
Long64_t lastsorted; // the last event whose tracks has been sorted |
359 |
// -------------------------------------- |
// -------------------------------------- |
360 |
// flags to load level0-1-2 trees/braches |
// flags to load level0-1-2 trees/braches |
361 |
// -------------------------------------- |
// -------------------------------------- |
392 |
|
|
393 |
|
|
394 |
public: |
public: |
395 |
|
|
396 |
PamLevel2(TString ddir,TList *list,TString detlist); |
PamLevel2(TString ddir,TList *list,TString detlist); |
397 |
PamLevel2(TString ddir,TString list,TString detlist); |
PamLevel2(TString ddir,TString list,TString detlist); |
398 |
PamLevel2(TString ddir,TString list); |
PamLevel2(TString ddir,TString list); |
399 |
PamLevel2(); |
PamLevel2(); |
400 |
~PamLevel2(){Delete();}; |
~PamLevel2(){Delete();}; |
401 |
|
|
402 |
void Clear(); |
void Clear(); |
403 |
void Reset(); |
void Reset(); |
404 |
void Delete(); |
void Delete(); |
405 |
|
|
406 |
|
|
407 |
// ============================================= |
// ============================================= |
408 |
// methods needed to read/write Level0-1-2 trees |
// methods needed to read/write Level0-1-2 trees |
409 |
// ============================================= |
// ============================================= |
411 |
void GetWhichTrees(TFile*); |
void GetWhichTrees(TFile*); |
412 |
|
|
413 |
TTree* GetPamTree(TFile*,TString); |
TTree* GetPamTree(TFile*,TString); |
414 |
TTree* GetPamTree(TFile* f){return GetPamTree(f,"");}; |
TTree* GetPamTree(TFile* f){return GetPamTree(f,"");}; |
415 |
|
|
416 |
TList* GetListOfLevel2Files(TString, TString); |
TList* GetListOfLevel2Files(TString, TString); |
417 |
Bool_t CheckLevel2File(TString); |
Bool_t CheckLevel2File(TString); |
427 |
Bool_t UpdateRunInfo(TTree *,Long64_t); |
Bool_t UpdateRunInfo(TTree *,Long64_t); |
428 |
Bool_t UpdateRunInfo(TChain *,Long64_t); |
Bool_t UpdateRunInfo(TChain *,Long64_t); |
429 |
|
|
430 |
|
void SkipRunInfoUpdate(){fUpdateRunInfo = false;}; // Emiliano |
431 |
|
void DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano |
432 |
|
void NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano |
433 |
|
|
434 |
void CreateCloneTrees0(TChain *fChain, TFile* file); |
void CreateCloneTrees0(TChain *fChain, TFile* file); |
435 |
void CreateCloneTrees(TFile* file); |
void CreateCloneTrees(TFile* file); |
436 |
void CreateCloneTrees(){CreateCloneTrees(NULL);}; |
void CreateCloneTrees(){CreateCloneTrees(NULL);}; |
443 |
|
|
444 |
// Int_t GetEntry(Int_t); |
// Int_t GetEntry(Int_t); |
445 |
Int_t GetEntry(Long64_t); |
Int_t GetEntry(Long64_t); |
446 |
Int_t GetYodaEntry(); |
Long64_t GetReadEntry(); |
447 |
|
Int_t GetYodaEntry(); |
448 |
|
|
449 |
void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable |
void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable |
450 |
|
|
451 |
Long64_t GetEntries(){return pam_tree->GetEntries();}; |
Long64_t GetEntries(){return pam_tree->GetEntries();}; |
452 |
TChain* GetPamTree(){return pam_tree;}; |
TChain* GetPamTree(){return pam_tree;}; |
453 |
TChain* GetRunTree(){return run_tree;}; |
TChain* GetRunTree(){return run_tree;}; |
502 |
PamTrack* GetPamTrackAlong(TrkTrack*); |
PamTrack* GetPamTrackAlong(TrkTrack*); |
503 |
// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve |
// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve |
504 |
// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve |
// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve |
505 |
|
|
506 |
// ================================= |
// ================================= |
507 |
// methods to retrieve pamela tracks |
// methods to retrieve pamela tracks |
508 |
// ================================= |
// ================================= |
514 |
TClonesArray* GetTracks(); //return sorted_tracks |
TClonesArray* GetTracks(); //return sorted_tracks |
515 |
PamTrack* GetTrack(int); |
PamTrack* GetTrack(int); |
516 |
PamTrack* GetTrackImage(int); |
PamTrack* GetTrackImage(int); |
517 |
void SetSortingMethod(TString how){ howtosort = how; return;}; |
void SetSortingMethod(TString how); |
518 |
TString GetSortingMethod(){return howtosort;}; |
TString GetSortingMethod(){return howtosort;}; |
519 |
void SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;}; |
void SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;}; |
520 |
Float_t GetNucleiSortingThreshold(){ return sortthr;}; |
Float_t GetNucleiSortingThreshold(){ return sortthr;}; |
521 |
|
|
522 |
// ================================= |
// ================================= |
523 |
// other methods |
// other methods |
524 |
// ================================= |
// ================================= |
554 |
Bool_t IsGood(Bool_t strict=true); |
Bool_t IsGood(Bool_t strict=true); |
555 |
|
|
556 |
Int_t IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list) |
Int_t IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list) |
557 |
|
void SetSELLI(Int_t selli){SELLI=selli;}; |
558 |
|
|
559 |
Bool_t IsGP(){return ISGP;}; |
Bool_t IsGP(){return ISGP;}; |
560 |
|
|
562 |
|
|
563 |
void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); }; |
void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); }; |
564 |
|
|
565 |
|
|
566 |
ClassDef(PamLevel2,1); |
ClassDef(PamLevel2,1); |
567 |
}; |
}; |
568 |
|
|