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 |
|
|
291 |
// ------------------------------------------- |
// ------------------------------------------- |
292 |
// variables to store total live and dead time |
// variables to store total live and dead time |
293 |
// ------------------------------------------- |
// ------------------------------------------- |
294 |
ULong64_t totdltime[3];// 0=live 1=dead 2=n.events |
ULong64_t totdltime[3];// 0=live 1=dead 2=n.events |
295 |
|
|
296 |
// --------------------------------- |
// --------------------------------- |
297 |
// variables to retrieve level0 info |
// variables to retrieve level0 info |
298 |
// --------------------------------- |
// --------------------------------- |
299 |
TFile* l0_file; |
TFile* l0_file; |
300 |
TTree* l0_tree; |
TTree* l0_tree; |
301 |
Int_t iroot; |
Int_t iroot; |
302 |
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 |
303 |
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 |
304 |
// -------------------------------------------- |
// -------------------------------------------- |
305 |
// variables to store Level0 info (if required) |
// variables to store Level0 info (if required) |
306 |
// -------------------------------------------- |
// -------------------------------------------- |
323 |
TTree* pam_tree_clone[NCLONES]; |
TTree* pam_tree_clone[NCLONES]; |
324 |
|
|
325 |
// -------------------------------------------- |
// -------------------------------------------- |
326 |
// variables to store Level1-2 info |
// variables to store Level1-2 info |
327 |
// -------------------------------------------- |
// -------------------------------------------- |
328 |
TrkLevel1* trk1_obj; |
TrkLevel1* trk1_obj; |
329 |
TrkLevel2* trk2_obj; |
TrkLevel2* trk2_obj; |
349 |
// ----------------------------------- |
// ----------------------------------- |
350 |
TClonesArray *tsorted; |
TClonesArray *tsorted; |
351 |
TClonesArray *timage; |
TClonesArray *timage; |
352 |
PamTrack *track; //store the last required |
PamTrack *track; //store the last required |
353 |
TString howtosort; |
TString howtosort; |
354 |
Float_t sortthr; |
Float_t sortthr; |
355 |
|
|
356 |
|
Bool_t issorted; // true if the current event has already been sorted with the current sort criteria |
357 |
|
Long64_t lastsorted; // the last event whose tracks has been sorted |
358 |
// -------------------------------------- |
// -------------------------------------- |
359 |
// flags to load level0-1-2 trees/braches |
// flags to load level0-1-2 trees/braches |
360 |
// -------------------------------------- |
// -------------------------------------- |
391 |
|
|
392 |
|
|
393 |
public: |
public: |
394 |
|
|
395 |
PamLevel2(TString ddir,TList *list,TString detlist); |
PamLevel2(TString ddir,TList *list,TString detlist); |
396 |
PamLevel2(TString ddir,TString list,TString detlist); |
PamLevel2(TString ddir,TString list,TString detlist); |
397 |
PamLevel2(TString ddir,TString list); |
PamLevel2(TString ddir,TString list); |
398 |
PamLevel2(); |
PamLevel2(); |
399 |
~PamLevel2(){Delete();}; |
~PamLevel2(){Delete();}; |
400 |
|
|
401 |
void Clear(); |
void Clear(); |
402 |
void Reset(); |
void Reset(); |
403 |
void Delete(); |
void Delete(); |
404 |
|
|
405 |
|
|
406 |
// ============================================= |
// ============================================= |
407 |
// methods needed to read/write Level0-1-2 trees |
// methods needed to read/write Level0-1-2 trees |
408 |
// ============================================= |
// ============================================= |
410 |
void GetWhichTrees(TFile*); |
void GetWhichTrees(TFile*); |
411 |
|
|
412 |
TTree* GetPamTree(TFile*,TString); |
TTree* GetPamTree(TFile*,TString); |
413 |
TTree* GetPamTree(TFile* f){return GetPamTree(f,"");}; |
TTree* GetPamTree(TFile* f){return GetPamTree(f,"");}; |
414 |
|
|
415 |
TList* GetListOfLevel2Files(TString, TString); |
TList* GetListOfLevel2Files(TString, TString); |
416 |
Bool_t CheckLevel2File(TString); |
Bool_t CheckLevel2File(TString); |
438 |
|
|
439 |
// Int_t GetEntry(Int_t); |
// Int_t GetEntry(Int_t); |
440 |
Int_t GetEntry(Long64_t); |
Int_t GetEntry(Long64_t); |
441 |
Int_t GetYodaEntry(); |
Long64_t GetReadEntry(); |
442 |
|
Int_t GetYodaEntry(); |
443 |
|
|
444 |
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 |
445 |
|
|
446 |
Long64_t GetEntries(){return pam_tree->GetEntries();}; |
Long64_t GetEntries(){return pam_tree->GetEntries();}; |
447 |
TChain* GetPamTree(){return pam_tree;}; |
TChain* GetPamTree(){return pam_tree;}; |
448 |
TChain* GetRunTree(){return run_tree;}; |
TChain* GetRunTree(){return run_tree;}; |
497 |
PamTrack* GetPamTrackAlong(TrkTrack*); |
PamTrack* GetPamTrackAlong(TrkTrack*); |
498 |
// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve |
// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve |
499 |
// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve |
// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve |
500 |
|
|
501 |
// ================================= |
// ================================= |
502 |
// methods to retrieve pamela tracks |
// methods to retrieve pamela tracks |
503 |
// ================================= |
// ================================= |
509 |
TClonesArray* GetTracks(); //return sorted_tracks |
TClonesArray* GetTracks(); //return sorted_tracks |
510 |
PamTrack* GetTrack(int); |
PamTrack* GetTrack(int); |
511 |
PamTrack* GetTrackImage(int); |
PamTrack* GetTrackImage(int); |
512 |
void SetSortingMethod(TString how){ howtosort = how; return;}; |
void SetSortingMethod(TString how); |
513 |
TString GetSortingMethod(){return howtosort;}; |
TString GetSortingMethod(){return howtosort;}; |
514 |
void SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;}; |
void SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;}; |
515 |
Float_t GetNucleiSortingThreshold(){ return sortthr;}; |
Float_t GetNucleiSortingThreshold(){ return sortthr;}; |
516 |
|
|
517 |
// ================================= |
// ================================= |
518 |
// other methods |
// other methods |
519 |
// ================================= |
// ================================= |
549 |
Bool_t IsGood(Bool_t strict=true); |
Bool_t IsGood(Bool_t strict=true); |
550 |
|
|
551 |
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) |
552 |
void SetSELLI(Int_t selli){SELLI=selli;}; |
void SetSELLI(Int_t selli){SELLI=selli;}; |
553 |
|
|
554 |
Bool_t IsGP(){return ISGP;}; |
Bool_t IsGP(){return ISGP;}; |
555 |
|
|
557 |
|
|
558 |
void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); }; |
void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); }; |
559 |
|
|
560 |
|
|
561 |
ClassDef(PamLevel2,1); |
ClassDef(PamLevel2,1); |
562 |
}; |
}; |
563 |
|
|