PamelaLevel2/inc/PamLevel2.h



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

#include <TRefArray.h>

#include <TrkStruct.h>
#include <TrkParams.h>
#include <TrkLevel0.h>
#include <TrkLevel1.h>
#include <TrkLevel2.h>
#include <TrkHough.h>
#include <CaloLevel0.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 <RunInfo.h>
#include <GLTables.h>


#include <TTree.h>
#include <TFriendElement.h>
#include <TChain.h>
#include <TChainElement.h>
#include <TFile.h>
#include <TList.h>
#include <TKey.h>
#include <TSystemFile.h>
#include <TSystemDirectory.h>
#include <TSQLServer.h>

#include <stdlib.h>
#include <iostream>
#include <iomanip>
using namespace std;

#define NCLONES 10

using namespace std;


/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////
/**
 * \brief Class to describe GP output (to be moved to PamelaDigitizer)
 * 
 */

class GPamela : public TObject {

 private:

 public:

    // Declaration of leave types
    Int_t           Irun;
    Int_t           Ievnt;
    Int_t           Ipa;
    Float_t         X0;
    Float_t         Y0;
    Float_t         Z0;
    Float_t         Theta;
    Float_t         Phi;
    Float_t         P0;
    Int_t           Nthtof;
    UChar_t         Ipltof[100];   //[Nthtof]
    UChar_t         Ipaddle[100];   //[Nthtof]
    UChar_t         Ipartof[100];   //[Nthtof]
    Float_t         Xintof[100];   //[Nthtof]
    Float_t         Yintof[100];   //[Nthtof]
    Float_t         Zintof[100];   //[Nthtof]
    Float_t         Xouttof[100];   //[Nthtof]
    Float_t         Youttof[100];   //[Nthtof]
    Float_t         Zouttof[100];   //[Nthtof]
    Float_t         Ereltof[100];   //[Nthtof]
    Float_t         Timetof[100];   //[Nthtof]
    Float_t         Pathtof[100];   //[Nthtof]
    Float_t         P0tof[100];   //[Nthtof]
    Int_t           Nthcat;
    UChar_t         Iparcat[50];   //[Nthcat]
    UChar_t         Icat[50];   //[Nthcat]
    Float_t         Xincat[50];   //[Nthcat]
    Float_t         Yincat[50];   //[Nthcat]
    Float_t         Zincat[50];   //[Nthcat]
    Float_t         Xoutcat[50];   //[Nthcat]
    Float_t         Youtcat[50];   //[Nthcat]
    Float_t         Zoutcat[50];   //[Nthcat]
    Float_t         Erelcat[50];   //[Nthcat]
    Float_t         Timecat[50];   //[Nthcat]
    Float_t         Pathcat[50];   //[Nthcat]
    Float_t         P0cat[50];   //[Nthcat]
    Int_t           Nthcas;
    UChar_t         Iparcas[50];   //[Nthcas]
    UChar_t         Icas[50];   //[Nthcas]
    Float_t         Xincas[50];   //[Nthcas]
    Float_t         Yincas[50];   //[Nthcas]
    Float_t         Zincas[50];   //[Nthcas]
    Float_t         Xoutcas[50];   //[Nthcas]
    Float_t         Youtcas[50];   //[Nthcas]
    Float_t         Zoutcas[50];   //[Nthcas]
    Float_t         Erelcas[50];   //[Nthcas]
    Float_t         Timecas[50];   //[Nthcas]
    Float_t         Pathcas[50];   //[Nthcas]
    Float_t         P0cas[50];   //[Nthcas]
    Int_t           Nthspe;
    UChar_t         Iparspe[100];   //[Nthspe]
    UChar_t         Itrpb[100];   //[Nthspe]
    UChar_t         Itrsl[100];   //[Nthspe]
    UChar_t         Itspa[100];   //[Nthspe]
    Float_t         Xinspe[100];   //[Nthspe]
    Float_t         Yinspe[100];   //[Nthspe]
    Float_t         Zinspe[100];   //[Nthspe]
    Float_t         Xoutspe[100];   //[Nthspe]
    Float_t         Youtspe[100];   //[Nthspe]
    Float_t         Zoutspe[100];   //[Nthspe]
    Float_t         Xavspe[100];   //[Nthspe]
    Float_t         Yavspe[100];   //[Nthspe]
    Float_t         Zavspe[100];   //[Nthspe]
    Float_t         Erelspe[100];   //[Nthspe]
    Float_t         Pathspe[100];   //[Nthspe]
    Float_t         P0spe[100];   //[Nthspe]
    UChar_t         Nxmult[100];   //[Nthspe]
    UChar_t         Nymult[100];   //[Nthspe]
    Int_t           Nstrpx;
    UChar_t         Npstripx[1000];   //[Nstrpx]
    UChar_t         Ntstripx[1000];   //[Nstrpx]
    UShort_t        Istripx[1000];   //[Nstrpx]
    //UInt_t        Istripx[1000];   //[Nstrpx]
    Float_t         Qstripx[1000];   //[Nstrpx]
    Float_t         Xstripx[1000];   //[Nstrpx]
    Int_t           Nstrpy;
    UChar_t         Npstripy[1000];   //[Nstrpy]
    UChar_t         Ntstripy[1000];   //[Nstrpy]
    UShort_t        Istripy[1000];   //[Nstrpy]
    //UInt_t        Istripy[1000];   //[Nstrpy]
    Float_t         Qstripy[1000];   //[Nstrpy]
    Float_t         Ystripy[1000];   //[Nstrpy]
    Int_t           Nthcali;
    UChar_t         Icaplane[4224];   //[Nthcali]
    UChar_t         Icastrip[4224];   //[Nthcali]
    Int_t           Icamod[4224];   //[Nthcali]
    Float_t         Enestrip[4224];   //[Nthcali]
    Int_t           Nthcal;
    UChar_t         Icapl[1000];   //[Nthcal]
    UChar_t         Icasi[1000];   //[Nthcal]
    UChar_t         Icast[1000];   //[Nthcal]
    Float_t         Xincal[1000];   //[Nthcal]
    Float_t         Yincal[1000];   //[Nthcal]
    Float_t         Zincal[1000];   //[Nthcal]
    Float_t         Erelcal[1000];   //[Nthcal]
    Int_t           Nthnd;
    UChar_t         Itubend[200];   //[Nthnd]
    UChar_t         Iparnd[200];   //[Nthnd]
    Float_t         Xinnd[200];   //[Nthnd]
    Float_t         Yinnd[200];   //[Nthnd]
    Float_t         Zinnd[200];   //[Nthnd]
    Float_t         Xoutnd[200];   //[Nthnd]
    Float_t         Youtnd[200];   //[Nthnd]
    Float_t         Zoutnd[200];   //[Nthnd]
    Float_t         Erelnd[200];   //[Nthnd]
    Float_t         Timend[200];   //[Nthnd]
    Float_t         Pathnd[200];   //[Nthnd]
    Float_t         P0nd[200];   //[Nthnd]
    Int_t           Nthcard;
    UChar_t         Iparcard[100];   //[Nthcard]
    UChar_t         Icard[100];   //[Nthcard]
    Float_t         Xincard[100];   //[Nthcard]
    Float_t         Yincard[100];   //[Nthcard]
    Float_t         Zincard[100];   //[Nthcard]
    Float_t         Xoutcard[100];   //[Nthcard]
    Float_t         Youtcard[100];   //[Nthcard]
    Float_t         Zoutcard[100];   //[Nthcard]
    Float_t         Erelcard[100];   //[Nthcard]
    Float_t         Timecard[100];   //[Nthcard]
    Float_t         Pathcard[100];   //[Nthcard]
    Float_t         P0card[100];   //[Nthcard]

    GPamela(){ Clear(); };
//    ~GPamela(){ Delete(); };
    ~GPamela(){ Clear(); };

    void Clear();
    void Delete();
    void SetBranchAddress( TChain* );
    void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};

    GPamela*    GetGPamela(){return this;};

    ClassDef(GPamela,1);

};
/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////


/**
 * \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 {
class PamTrack : public TObject {

 private:

    TrkTrack*   trk_track;
    CaloTrkVar* calo_track;
    ToFTrkVar*  tof_track;
    Int_t pscore;
    Int_t iscore;
    bool candeleteobj;

//    TClonesArray *Cluster;        ///< clusters (if stored)

 public:

    PamTrack();
    PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );
    PamTrack(const PamTrack&);
    ~PamTrack(){ Delete(); };

    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();

    Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
    Int_t GetScore(){return (pscore-iscore);};
    Int_t GetPScore(){return (pscore);};
    Int_t GetIScore(){return (iscore);};
    void SetPScore(Int_t p){pscore=p;};
    void SetIScore(Int_t i){iscore=i;};

    TrkTrack*   GetTrkTrack(){return trk_track;};
    CaloTrkVar* GetCaloTrack(){return calo_track;};
    ToFTrkVar*  GetToFTrack(){return tof_track;};

    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 {
class PamLevel2 : public TObject {

 private:

    // ------------------------------
    // variables to retrieve run info
    // ------------------------------
    Long64_t  irun;          //current run entry, relative to the TChain
    Long64_t  irunt;         //current run entry, relative to the TTree in the TChain
    Long64_t  runfirstentry; //first entry of run, relative to the TChain
    Long64_t  runlastentry;  //last  entry of run, relative to the TChain
    Long64_t  irunentry;     //position of current entry, relative to the run
    TChain*   run_tree;
    TTree*    run_tree_clone;
    TChain*   sel_tree;
    TTree*    sel_tree_clone;
    int *irunoffset;

    // -------------------------------------------
    // variables to store total live and dead time
    // -------------------------------------------
    ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events   

    // ---------------------------------
    // variables to retrieve level0 info
    // ---------------------------------
    TFile*      l0_file;
    TTree*      l0_tree;
    Int_t       iroot;      
    Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time
    Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10 
    // --------------------------------------------
    // variables to store Level0 info (if required)
    // --------------------------------------------
    EventHeader* h0_obj;
    TrkLevel0*   trk0_obj;
    CaloLevel0*  calo0_obj;

    // ---------------------------
    // variables to connect the DB
    // ---------------------------
    TSQLServer* dbc;
    TString host;
    TString user;
    TString psw;

    // ---------------------------------------------
    // variables to retrieve and clone level1-2 info
    // ---------------------------------------------
    TChain* pam_tree;
    TTree*  pam_tree_clone[NCLONES];

    // --------------------------------------------
    // variables to store Level1-2 info 
    // --------------------------------------------
    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;
    SoftInfo*    soft_obj; // Emiliano

    // --------------------------------------------
    // variables to handle GP
    // --------------------------------------------
    GPamela*     gp_obj;

    // -----------------------------------
    // variables to sort the track
    // -----------------------------------
    TClonesArray *tsorted;
    TClonesArray *timage;
    PamTrack    *track; //store the last required 
    TString howtosort;
    Float_t sortthr;

    // --------------------------------------
    // flags to load level0-1-2 trees/braches
    // --------------------------------------
    Bool_t CAL2; //level2 branch
    Bool_t CAL1; //level1 branch
    Bool_t CAL0; //yoda

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

    Bool_t TOF;
    Bool_t TOF0;

    Bool_t TRG;
    Bool_t S4;
    Bool_t ND;
    Bool_t AC;
    Bool_t ORB;

    Bool_t GP; //GP simulation tree

    // -----------
    // other flags
    // -----------
    Bool_t RUN;   // RunInfo tree
    Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
    Int_t  SELLI; // tree of selected events (there is a selection list)

    void Initialize();


 public:
         
    PamLevel2(TString ddir,TString list,TString detlist);
    PamLevel2(TString ddir,TString list);
    PamLevel2();
    ~PamLevel2(){Delete();};
        
    void Clear();
    void Reset();
    void Delete();

         
//  =============================================
//  methods needed to read/write Level0-1-2 trees
//  =============================================
    void    SetWhichTrees(TString);
    void    GetWhichTrees(TFile*);

    TTree*  GetPamTree(TFile*,TString);
    TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};    

    TList*  GetListOfLevel2Files(TString, TString);
    Bool_t  CheckLevel2File(TString);
    TChain* GetPamTree(TList*,TString);
    TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };
    TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
    TChain* AddFriend(TString);

    TTree*  GetRunTree(TFile*);
    TChain* GetRunTree(TList*);
    TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
    Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
    Bool_t  UpdateRunInfo(TTree *,Long64_t);
    Bool_t  UpdateRunInfo(TChain *,Long64_t);

    void    CreateCloneTrees0(TChain *fChain, TFile* file);
    void    CreateCloneTrees(TFile* file);
    void    CreateCloneTrees(){CreateCloneTrees(NULL);};
    void    FillCloneTrees();
    TTree  *GetCloneTree(TString);
    void    WriteCloneTrees();

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

//    Int_t   GetEntry(Int_t);
    Int_t   GetEntry(Long64_t);
    Int_t   GetYodaEntry(); 

    void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
    
    Long64_t GetEntries(){return pam_tree->GetEntries();};
    TChain* GetPamTree(){return pam_tree;};
    TChain* GetRunTree(){return run_tree;};
    TTree*  GetYodaTree();

    TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things

//  =============
//  DB connection
//  =============
    void SetDBHost(TString str){host=str;};
    void SetDBUser(TString str){user=str;};
    void SetDBPsw(TString str){psw=str;};
    TSQLServer* GetSQLServer(){return dbc;};
    Bool_t SetDBConnection();

//  ========================================================
//  methods needed to retrieve subdetector general variables
//  ========================================================
    void*       GetPointerTo(const char*);
    EventHeader  *GetEventHeader()  { return h0_obj; };

/*     TrkLevel0    *GetTrkLevel0()   { return trk0_obj; }; */
/*     CaloLevel0   *GetCaloLevel0()  { return calo0_obj; }; */
    TrkLevel0    *GetTrkLevel0();
    CaloLevel0   *GetCaloLevel0();


    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;  };
    GL_RUN       *GetRunInfo()    { return run_obj;  };
    SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
    GPamela      *GetGPamela()    { return gp_obj; };

//  ==============================================
//  methods to retrieve subdetectors stored tracks
//  ==============================================
    TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->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
    void          SortTracks();          //assign value to variable sorted_tracks
//    TRefArray*    GetTracks();                  //return sorted_tracks
    TClonesArray* GetTracks();                  //return sorted_tracks
    PamTrack*     GetTrack(int);
    PamTrack*     GetTrackImage(int);
    void          SetSortingMethod(TString how){ howtosort = how; return;};
    TString       GetSortingMethod(){return howtosort;};
    void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
    Float_t       GetNucleiSortingThreshold(){ return sortthr;};
        
//  =================================
//  other methods
//  =================================
//    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };

    // -----
    // flags methods
    // -----
    Bool_t IsCAL2(){return CAL2;}; //level2 branch
    Bool_t IsCAL1(){return CAL1;}; //level1 branch
    Bool_t IsCAL0(){return CAL0;}; //yoda

    Bool_t IsTRK2(){return TRK2;}; //level2 branch
    Bool_t IsTRK1(){return TRK1;}; //level1 branch
    Bool_t IsTRKh(){return TRKh;}; //hough branch
    Bool_t IsTRK0(){return TRK0;}; //yoda

    Bool_t IsTOF(){return TOF;};
    Bool_t IsTOF0(){return TOF0;};

    Bool_t IsTRG(){return TRG;};
    Bool_t IsS4(){return S4;};
    Bool_t IsND(){return ND;};
    Bool_t IsAC(){return AC;};
    Bool_t IsORB(){return ORB;};

    Bool_t IsRUN(){return RUN;};

    Bool_t IsGood();

    Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)

    Bool_t IsGP(){return ISGP;};

    void   SetGP(Bool_t tt){ISGP=tt;};


    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.13	#include <TrkStruct.h>
13			#include <TrkParams.h>
14			#include <TrkLevel0.h>
15	pam-fi	1.4	#include <TrkLevel1.h>
16	pam-fi	1.1	#include <TrkLevel2.h>
17	pam-fi	1.6	#include <TrkHough.h>
18	mocchiut	1.16	#include <CaloLevel0.h>
19	pam-fi	1.6	#include <CaloLevel1.h>
20	pam-fi	1.1	#include <CaloLevel2.h>
21			#include <ToFLevel2.h>
22			#include <TrigLevel2.h>
23			#include <S4Level2.h>
24			#include <NDLevel2.h>
25			#include <AcLevel2.h>
26	pam-fi	1.2	#include <OrbitalInfo.h>
27	pam-fi	1.4	#include <RunInfo.h>
28			#include <GLTables.h>
29
30	pam-fi	1.1
31			#include <TTree.h>
32	pam-fi	1.7	#include <TFriendElement.h>
33	pam-fi	1.2	#include <TChain.h>
34	pam-fi	1.26	#include <TChainElement.h>
35	pam-fi	1.1	#include <TFile.h>
36	pam-fi	1.2	#include <TList.h>
37			#include <TKey.h>
38			#include <TSystemFile.h>
39			#include <TSystemDirectory.h>
40	pam-fi	1.13	#include <TSQLServer.h>
41	pam-fi	1.1
42	pam-fi	1.26	#include <stdlib.h>
43	pam-fi	1.1	#include <iostream>
44	pam-fi	1.26	#include <iomanip>
45			using namespace std;
46	pam-fi	1.13
47			#define NCLONES 10
48
49	pam-fi	1.1	using namespace std;
50
51	pam-fi	1.30
52			/////////////////////////////////////////////////////////////////////////////
53			/////////////////////////////////////////////////////////////////////////////
54			/////////////////////////////////////////////////////////////////////////////
55			/////////////////////////////////////////////////////////////////////////////
56			/**
57			* \brief Class to describe GP output (to be moved to PamelaDigitizer)
58			*
59			*/
60
61			class GPamela : public TObject {
62
63			private:
64
65			public:
66
67			// Declaration of leave types
68			Int_t Irun;
69			Int_t Ievnt;
70			Int_t Ipa;
71			Float_t X0;
72			Float_t Y0;
73			Float_t Z0;
74			Float_t Theta;
75			Float_t Phi;
76			Float_t P0;
77			Int_t Nthtof;
78			UChar_t Ipltof[100]; //[Nthtof]
79			UChar_t Ipaddle[100]; //[Nthtof]
80			UChar_t Ipartof[100]; //[Nthtof]
81			Float_t Xintof[100]; //[Nthtof]
82			Float_t Yintof[100]; //[Nthtof]
83			Float_t Zintof[100]; //[Nthtof]
84			Float_t Xouttof[100]; //[Nthtof]
85			Float_t Youttof[100]; //[Nthtof]
86			Float_t Zouttof[100]; //[Nthtof]
87			Float_t Ereltof[100]; //[Nthtof]
88			Float_t Timetof[100]; //[Nthtof]
89			Float_t Pathtof[100]; //[Nthtof]
90			Float_t P0tof[100]; //[Nthtof]
91			Int_t Nthcat;
92			UChar_t Iparcat[50]; //[Nthcat]
93			UChar_t Icat[50]; //[Nthcat]
94			Float_t Xincat[50]; //[Nthcat]
95			Float_t Yincat[50]; //[Nthcat]
96			Float_t Zincat[50]; //[Nthcat]
97			Float_t Xoutcat[50]; //[Nthcat]
98			Float_t Youtcat[50]; //[Nthcat]
99			Float_t Zoutcat[50]; //[Nthcat]
100			Float_t Erelcat[50]; //[Nthcat]
101			Float_t Timecat[50]; //[Nthcat]
102			Float_t Pathcat[50]; //[Nthcat]
103			Float_t P0cat[50]; //[Nthcat]
104			Int_t Nthcas;
105			UChar_t Iparcas[50]; //[Nthcas]
106			UChar_t Icas[50]; //[Nthcas]
107			Float_t Xincas[50]; //[Nthcas]
108			Float_t Yincas[50]; //[Nthcas]
109			Float_t Zincas[50]; //[Nthcas]
110			Float_t Xoutcas[50]; //[Nthcas]
111			Float_t Youtcas[50]; //[Nthcas]
112			Float_t Zoutcas[50]; //[Nthcas]
113			Float_t Erelcas[50]; //[Nthcas]
114			Float_t Timecas[50]; //[Nthcas]
115			Float_t Pathcas[50]; //[Nthcas]
116			Float_t P0cas[50]; //[Nthcas]
117			Int_t Nthspe;
118			UChar_t Iparspe[100]; //[Nthspe]
119			UChar_t Itrpb[100]; //[Nthspe]
120			UChar_t Itrsl[100]; //[Nthspe]
121			UChar_t Itspa[100]; //[Nthspe]
122			Float_t Xinspe[100]; //[Nthspe]
123			Float_t Yinspe[100]; //[Nthspe]
124			Float_t Zinspe[100]; //[Nthspe]
125			Float_t Xoutspe[100]; //[Nthspe]
126			Float_t Youtspe[100]; //[Nthspe]
127			Float_t Zoutspe[100]; //[Nthspe]
128			Float_t Xavspe[100]; //[Nthspe]
129			Float_t Yavspe[100]; //[Nthspe]
130			Float_t Zavspe[100]; //[Nthspe]
131			Float_t Erelspe[100]; //[Nthspe]
132			Float_t Pathspe[100]; //[Nthspe]
133			Float_t P0spe[100]; //[Nthspe]
134			UChar_t Nxmult[100]; //[Nthspe]
135			UChar_t Nymult[100]; //[Nthspe]
136			Int_t Nstrpx;
137			UChar_t Npstripx[1000]; //[Nstrpx]
138			UChar_t Ntstripx[1000]; //[Nstrpx]
139			UShort_t Istripx[1000]; //[Nstrpx]
140			//UInt_t Istripx[1000]; //[Nstrpx]
141			Float_t Qstripx[1000]; //[Nstrpx]
142			Float_t Xstripx[1000]; //[Nstrpx]
143			Int_t Nstrpy;
144			UChar_t Npstripy[1000]; //[Nstrpy]
145			UChar_t Ntstripy[1000]; //[Nstrpy]
146			UShort_t Istripy[1000]; //[Nstrpy]
147			//UInt_t Istripy[1000]; //[Nstrpy]
148			Float_t Qstripy[1000]; //[Nstrpy]
149			Float_t Ystripy[1000]; //[Nstrpy]
150			Int_t Nthcali;
151			UChar_t Icaplane[4224]; //[Nthcali]
152			UChar_t Icastrip[4224]; //[Nthcali]
153			Int_t Icamod[4224]; //[Nthcali]
154			Float_t Enestrip[4224]; //[Nthcali]
155			Int_t Nthcal;
156			UChar_t Icapl[1000]; //[Nthcal]
157			UChar_t Icasi[1000]; //[Nthcal]
158			UChar_t Icast[1000]; //[Nthcal]
159			Float_t Xincal[1000]; //[Nthcal]
160			Float_t Yincal[1000]; //[Nthcal]
161			Float_t Zincal[1000]; //[Nthcal]
162			Float_t Erelcal[1000]; //[Nthcal]
163			Int_t Nthnd;
164			UChar_t Itubend[200]; //[Nthnd]
165			UChar_t Iparnd[200]; //[Nthnd]
166			Float_t Xinnd[200]; //[Nthnd]
167			Float_t Yinnd[200]; //[Nthnd]
168			Float_t Zinnd[200]; //[Nthnd]
169			Float_t Xoutnd[200]; //[Nthnd]
170			Float_t Youtnd[200]; //[Nthnd]
171			Float_t Zoutnd[200]; //[Nthnd]
172			Float_t Erelnd[200]; //[Nthnd]
173			Float_t Timend[200]; //[Nthnd]
174			Float_t Pathnd[200]; //[Nthnd]
175			Float_t P0nd[200]; //[Nthnd]
176			Int_t Nthcard;
177			UChar_t Iparcard[100]; //[Nthcard]
178			UChar_t Icard[100]; //[Nthcard]
179			Float_t Xincard[100]; //[Nthcard]
180			Float_t Yincard[100]; //[Nthcard]
181			Float_t Zincard[100]; //[Nthcard]
182			Float_t Xoutcard[100]; //[Nthcard]
183			Float_t Youtcard[100]; //[Nthcard]
184			Float_t Zoutcard[100]; //[Nthcard]
185			Float_t Erelcard[100]; //[Nthcard]
186			Float_t Timecard[100]; //[Nthcard]
187			Float_t Pathcard[100]; //[Nthcard]
188			Float_t P0card[100]; //[Nthcard]
189
190			GPamela(){ Clear(); };
191			// ~GPamela(){ Delete(); };
192			~GPamela(){ Clear(); };
193
194			void Clear();
195			void Delete();
196			void SetBranchAddress( TChain* );
197			void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};
198
199			GPamela* GetGPamela(){return this;};
200
201			ClassDef(GPamela,1);
202
203			};
204			/////////////////////////////////////////////////////////////////////////////
205			/////////////////////////////////////////////////////////////////////////////
206			/////////////////////////////////////////////////////////////////////////////
207			/////////////////////////////////////////////////////////////////////////////
208
209
210
211	pam-fi	1.1	/**
212			* \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).
213			*
214			* This class inherites from TrkTrack, CaloTrkVar ecc...
215			*
216			*/
217	pam-fi	1.10	//class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
218			class PamTrack : public TObject {
219	pam-fi	1.1
220			private:
221
222			TrkTrack* trk_track;
223			CaloTrkVar* calo_track;
224			ToFTrkVar* tof_track;
225	mocchiut	1.27	Int_t pscore;
226			Int_t iscore;
227	pam-fi	1.10	bool candeleteobj;
228
229			// TClonesArray *Cluster; ///< clusters (if stored)
230
231	pam-fi	1.1	public:
232
233			PamTrack();
234	pam-fi	1.10	PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );
235			PamTrack(const PamTrack&);
236			~PamTrack(){ Delete(); };
237	pam-fi	1.1
238	pam-fi	1.10	void Clear();
239			void Clear(Option_t *option){Clear();};
240			void Delete();
241	pam-fi	1.1
242	mocchiut	1.27	Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
243			Int_t GetScore(){return (pscore-iscore);};
244			Int_t GetPScore(){return (pscore);};
245			Int_t GetIScore(){return (iscore);};
246			void SetPScore(Int_t p){pscore=p;};
247			void SetIScore(Int_t i){iscore=i;};
248
249	pam-fi	1.10	TrkTrack* GetTrkTrack(){return trk_track;};
250	pam-fi	1.8	CaloTrkVar* GetCaloTrack(){return calo_track;};
251	pam-fi	1.10	ToFTrkVar* GetToFTrack(){return tof_track;};
252	pam-fi	1.8
253	pam-fi	1.1	ClassDef(PamTrack,1);
254			};
255
256			/**
257			* \brief Class to describe PAMELA LEVEL2 data.
258			*
259			* This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.
260			* This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.
261			*/
262
263	pam-fi	1.10	//class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {
264			class PamLevel2 : public TObject {
265	pam-fi	1.1
266			private:
267
268	pam-fi	1.13	// ------------------------------
269			// variables to retrieve run info
270			// ------------------------------
271	pam-fi	1.17	Long64_t irun; //current run entry, relative to the TChain
272	pam-fi	1.31	Long64_t irunt; //current run entry, relative to the TTree in the TChain
273	pam-fi	1.17	Long64_t runfirstentry; //first entry of run, relative to the TChain
274			Long64_t runlastentry; //last entry of run, relative to the TChain
275			Long64_t irunentry; //position of current entry, relative to the run
276	pam-fi	1.13	TChain* run_tree;
277			TTree* run_tree_clone;
278			TChain* sel_tree;
279			TTree* sel_tree_clone;
280	pam-fi	1.32	int *irunoffset;
281	pam-fi	1.26
282			// -------------------------------------------
283			// variables to store total live and dead time
284			// -------------------------------------------
285			ULong64_t totdltime[3];// 0=live 1=dead 2=n.events
286
287	pam-fi	1.13	// ---------------------------------
288			// variables to retrieve level0 info
289			// ---------------------------------
290			TFile* l0_file;
291			TTree* l0_tree;
292			Int_t iroot;
293	mocchiut	1.16	Int_t prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time
294	mocchiut	1.29	Int_t maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
295	pam-fi	1.13	// --------------------------------------------
296			// variables to store Level0 info (if required)
297			// --------------------------------------------
298			EventHeader* h0_obj;
299			TrkLevel0* trk0_obj;
300	pam-fi	1.17	CaloLevel0* calo0_obj;
301	pam-fi	1.13
302			// ---------------------------
303			// variables to connect the DB
304			// ---------------------------
305			TSQLServer* dbc;
306			TString host;
307			TString user;
308			TString psw;
309
310			// ---------------------------------------------
311			// variables to retrieve and clone level1-2 info
312			// ---------------------------------------------
313			TChain* pam_tree;
314			TTree* pam_tree_clone[NCLONES];
315	mocchiut	1.18
316	pam-fi	1.13	// --------------------------------------------
317	pam-fi	1.30	// variables to store Level1-2 info
318	pam-fi	1.13	// --------------------------------------------
319	pam-fi	1.6	TrkLevel1* trk1_obj;
320			TrkLevel2* trk2_obj;
321			TrkHough* trkh_obj;
322			CaloLevel1* calo1_obj;
323			CaloLevel2* calo2_obj;
324	pam-fi	1.2	ToFLevel2* tof_obj;
325			TrigLevel2* trig_obj;
326			S4Level2* s4_obj;
327			NDLevel2* nd_obj;
328			AcLevel2* ac_obj;
329	pam-fi	1.6	OrbitalInfo* orb_obj;
330			GL_RUN* run_obj;
331	mocchiut	1.11	SoftInfo* soft_obj; // Emiliano
332	pam-fi	1.10
333	pam-fi	1.30	// --------------------------------------------
334			// variables to handle GP
335			// --------------------------------------------
336			GPamela* gp_obj;
337
338	pam-fi	1.13	// -----------------------------------
339			// variables to sort the track
340			// -----------------------------------
341	pam-fi	1.10	TClonesArray *tsorted;
342			TClonesArray *timage;
343			PamTrack *track; //store the last required
344	mocchiut	1.20	TString howtosort;
345			Float_t sortthr;
346	pam-fi	1.7
347	pam-fi	1.30	// --------------------------------------
348			// flags to load level0-1-2 trees/braches
349			// --------------------------------------
350	pam-fi	1.6	Bool_t CAL2; //level2 branch
351			Bool_t CAL1; //level1 branch
352	pam-fi	1.13	Bool_t CAL0; //yoda
353	pam-fi	1.6
354			Bool_t TRK2; //level2 branch
355			Bool_t TRK1; //level1 branch
356			Bool_t TRKh; //hough branch
357	pam-fi	1.13	Bool_t TRK0; //yoda
358
359			Bool_t TOF;
360			Bool_t TOF0;
361	pam-fi	1.6
362			Bool_t TRG;
363			Bool_t S4;
364			Bool_t ND;
365			Bool_t AC;
366			Bool_t ORB;
367	pam-fi	1.9
368	pam-fi	1.30	Bool_t GP; //GP simulation tree
369	pam-fi	1.25
370	pam-fi	1.30	// -----------
371			// other flags
372			// -----------
373			Bool_t RUN; // RunInfo tree
374			Bool_t ISGP; // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
375			Int_t SELLI; // tree of selected events (there is a selection list)
376	pam-fi	1.13
377	mocchiut	1.18	void Initialize();
378
379	pam-fi	1.13
380			public:
381	pam-fi	1.6
382	mocchiut	1.18	PamLevel2(TString ddir,TString list,TString detlist);
383			PamLevel2(TString ddir,TString list);
384	pam-fi	1.1	PamLevel2();
385	pam-fi	1.10	~PamLevel2(){Delete();};
386	pam-fi	1.2
387	pam-fi	1.6	void Clear();
388	mocchiut	1.18	void Reset();
389	pam-fi	1.10	void Delete();
390	pam-fi	1.2
391	mocchiut	1.12
392	pam-fi	1.13	// =============================================
393			// methods needed to read/write Level0-1-2 trees
394			// =============================================
395	pam-fi	1.6	void SetWhichTrees(TString);
396	pam-fi	1.8	void GetWhichTrees(TFile*);
397
398			TTree* GetPamTree(TFile*,TString);
399	pam-fi	1.26	TTree* GetPamTree(TFile* f){return GetPamTree(f,"");};
400	pam-fi	1.8
401	pam-fi	1.6	TList* GetListOfLevel2Files(TString, TString);
402			Bool_t CheckLevel2File(TString);
403			TChain* GetPamTree(TList*,TString);
404	pam-fi	1.8	TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };
405			TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
406	pam-fi	1.26	TChain* AddFriend(TString);
407	pam-fi	1.8
408	pam-fi	1.6	TTree* GetRunTree(TFile*);
409			TChain* GetRunTree(TList*);
410	mocchiut	1.12	TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
411	pam-fi	1.28	Bool_t UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
412	pam-fi	1.17	Bool_t UpdateRunInfo(TTree *,Long64_t);
413			Bool_t UpdateRunInfo(TChain *,Long64_t);
414	pam-fi	1.7
415	pam-fi	1.8	void CreateCloneTrees0(TChain fChain, TFile file);
416			void CreateCloneTrees(TFile* file);
417	pam-fi	1.23	void CreateCloneTrees(){CreateCloneTrees(NULL);};
418	pam-fi	1.7	void FillCloneTrees();
419			TTree *GetCloneTree(TString);
420			void WriteCloneTrees();
421
422	pam-fi	1.6	void SetBranchAddress(TTree *);
423			void SetBranchAddress(TChain *);
424
425	pam-fi	1.17	// Int_t GetEntry(Int_t);
426			Int_t GetEntry(Long64_t);
427	pam-fi	1.13	Int_t GetYodaEntry();
428	mocchiut	1.21
429	mocchiut	1.29	void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
430
431	mocchiut	1.21	Long64_t GetEntries(){return pam_tree->GetEntries();};
432	pam-fi	1.13	TChain* GetPamTree(){return pam_tree;};
433			TChain* GetRunTree(){return run_tree;};
434			TTree* GetYodaTree();
435
436	mocchiut	1.16	TFile* GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
437
438	pam-fi	1.13	// =============
439			// DB connection
440			// =============
441			void SetDBHost(TString str){host=str;};
442			void SetDBUser(TString str){user=str;};
443			void SetDBPsw(TString str){psw=str;};
444	pam-fi	1.24	TSQLServer* GetSQLServer(){return dbc;};
445			Bool_t SetDBConnection();
446	pam-fi	1.13
447	pam-fi	1.10	// ========================================================
448			// methods needed to retrieve subdetector general variables
449			// ========================================================
450	pam-fi	1.6	void* GetPointerTo(const char*);
451	pam-fi	1.15	EventHeader *GetEventHeader() { return h0_obj; };
452	pam-fi	1.24
453			/* TrkLevel0 GetTrkLevel0() { return trk0_obj; }; /
454			/* CaloLevel0 GetCaloLevel0() { return calo0_obj; }; /
455			TrkLevel0 *GetTrkLevel0();
456			CaloLevel0 *GetCaloLevel0();
457
458	pam-fi	1.13
459			TrkLevel1 *GetTrkLevel1() { return trk1_obj; };
460			TrkLevel2 *GetTrkLevel2() { return trk2_obj; };
461			TrkHough *GetTrkHough() { return trkh_obj; };
462	pam-fi	1.10	CaloLevel1 *GetCaloLevel1() { return calo1_obj;};
463			CaloLevel2 *GetCaloLevel2() { return calo2_obj;};
464	pam-fi	1.13	ToFLevel2 *GetToFLevel2() { return tof_obj; };
465			TrigLevel2 *GetTrigLevel2() { return trig_obj; };
466			S4Level2 *GetS4Level2() { return s4_obj; };
467			NDLevel2 *GetNDLevel2() { return nd_obj; };
468			AcLevel2 *GetAcLevel2() { return ac_obj; };
469			OrbitalInfo *GetOrbitalInfo(){ return orb_obj; };
470			GL_RUN *GetRunInfo() { return run_obj; };
471	mocchiut	1.11	SoftInfo *GetSoftInfo() { return soft_obj; }; // Emiliano
472	pam-fi	1.30	GPamela *GetGPamela() { return gp_obj; };
473	pam-fi	1.3
474			// ==============================================
475			// methods to retrieve subdetectors stored tracks
476			// ==============================================
477	pam-fi	1.10	TrkTrack* GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
478	pam-fi	1.6	CaloTrkVar* GetCaloStoredTrack(int);
479			ToFTrkVar* GetToFStoredTrack(int);
480			//
481			PamTrack* GetPamTrackAlong(TrkTrack*);
482	pam-fi	1.3	// PamTrack* GetPamTrackAlong(CaloTrkVar*); //se serve
483			// PamTrack* GetPamTrackAlong(ToFTrkVar*); //se serve
484
485	pam-fi	1.1	// =================================
486			// methods to retrieve pamela tracks
487			// =================================
488	pam-fi	1.3	PamTrack* GetStoredTrack(Int_t);
489	pam-fi	1.6	//
490			void SortTracks(TString); //assign value to variable sorted_tracks
491	mocchiut	1.20	void SortTracks(); //assign value to variable sorted_tracks
492	pam-fi	1.10	// TRefArray* GetTracks(); //return sorted_tracks
493			TClonesArray* GetTracks(); //return sorted_tracks
494	pam-fi	1.6	PamTrack* GetTrack(int);
495	pam-fi	1.1	PamTrack* GetTrackImage(int);
496	mocchiut	1.20	void SetSortingMethod(TString how){ howtosort = how; return;};
497			TString GetSortingMethod(){return howtosort;};
498			void SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
499			Float_t GetNucleiSortingThreshold(){ return sortthr;};
500	pam-fi	1.2
501	pam-fi	1.6	// =================================
502			// other methods
503			// =================================
504	pam-fi	1.10	// Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
505	pam-fi	1.6
506	pamela	1.14	// -----
507			// flags methods
508			// -----
509			Bool_t IsCAL2(){return CAL2;}; //level2 branch
510			Bool_t IsCAL1(){return CAL1;}; //level1 branch
511			Bool_t IsCAL0(){return CAL0;}; //yoda
512
513			Bool_t IsTRK2(){return TRK2;}; //level2 branch
514			Bool_t IsTRK1(){return TRK1;}; //level1 branch
515			Bool_t IsTRKh(){return TRKh;}; //hough branch
516			Bool_t IsTRK0(){return TRK0;}; //yoda
517
518			Bool_t IsTOF(){return TOF;};
519			Bool_t IsTOF0(){return TOF0;};
520
521			Bool_t IsTRG(){return TRG;};
522			Bool_t IsS4(){return S4;};
523			Bool_t IsND(){return ND;};
524			Bool_t IsAC(){return AC;};
525			Bool_t IsORB(){return ORB;};
526
527			Bool_t IsRUN(){return RUN;};
528
529	mocchiut	1.22	Bool_t IsGood();
530
531	pamela	1.14	Int_t IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
532	pam-fi	1.6
533	pam-fi	1.25	Bool_t IsGP(){return ISGP;};
534
535			void SetGP(Bool_t tt){ISGP=tt;};
536
537
538	pam-fi	1.1	ClassDef(PamLevel2,1);
539			};
540
541			//===============================================================================================
542			//
543			//
544			//
545			//
546			//===============================================================================================
547
548
549			#endif
550