| 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 |
|
|