class TBrowser;
class AliEMCALJetMicroDst: public TNamed {
-
-
public:
- AliEMCALJetMicroDst(char *name="jetMicroDst",
- char *tit="jet Micro Dst for preparation of proposal");
+ AliEMCALJetMicroDst(const char *name="jetMicroDst",
+ const char *tit="jet Micro Dst for preparation of proposal");
+ AliEMCALJetMicroDst(const AliEMCALJetMicroDst &) : TNamed("", ""){
+ Fatal("cpy ctor", "not implemented") ; }
virtual ~AliEMCALJetMicroDst();
Bool_t Create(TFile *file);
Bool_t Create(const char *fname);
void FillJetsControl(); // 18-jan-2003
- Bool_t Open(const Int_t mode=1) {return Open(DefineName(mode));} // *MENU*
+ Bool_t Open(Int_t mode=1) {return Open(DefineName(mode));} // *MENU*
Bool_t Open(const char *fname); // *MENU*
- const Char_t* DefineName(const Int_t mode=1); // *MENU*
+ const Char_t* DefineName(Int_t mode=1); // *MENU*
Bool_t Initialize(TFile *file);
void Print(Option_t* option="") const; // *MENU*
Int_t GetEntry(Int_t entry);
void Test();
- Int_t GetNpart() {return fnpart;}
- Bool_t GetParton(Int_t i, Float_t& pt, Float_t& eta, Float_t& phi);
- Bool_t GetParton(Int_t i, TVector3& vec);
- Int_t GetNjet() {return fnjet;}
- Bool_t GetJet(Int_t i,Int_t mode, Float_t& pt,Float_t& eta,Float_t& phi);
- Bool_t GetJet(Int_t i,Int_t mode, TVector3& vec);
+ Int_t GetNpart() const {return fnpart;}
+ Bool_t GetParton(Int_t i, Float_t& pt, Float_t& eta, Float_t& phi) const ;
+ Bool_t GetParton(Int_t i, TVector3& vec) const ;
+ Int_t GetNjet() const {return fnjet;}
+ Bool_t GetJet(Int_t i,Int_t mode, Float_t& pt,Float_t& eta,Float_t& phi) const ;
+ Bool_t GetJet(Int_t i,Int_t mode, TVector3& vec) const ;
static void FillVector(Float_t pt, Float_t eta, Float_t phi, TVector3& vec);
- void GetEtaPhi(Int_t id, Double_t &eta, Double_t &phi);
- TVector3& GetCellVector(Int_t i);
- TVector3& GetGridVector(Int_t i);
+ void GetEtaPhi(Int_t id, Double_t &eta, Double_t &phi) const ;
+ TVector3& GetCellVector(Int_t i) const ;
+ TVector3& GetGridVector(Int_t i) const ;
// 13-apr-2003
- Double_t GetSumInCone(TVector3 &jet, Int_t nc, Float_t *et,Float_t *eta,Float_t *phi, Double_t cellEtCut, Double_t rJet);
- Double_t GetEmcalEtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5);
- Double_t GetTpcPtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5);
- Double_t GetSum(Int_t n, Float_t *ar, Double_t cut=0.0);
+ Double_t GetSumInCone(TVector3 &jet, Int_t nc, Float_t *et,Float_t *eta,Float_t *phi, Double_t cellEtCut, Double_t rJet) const ;
+ Double_t GetEmcalEtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5) ;
+ Double_t GetTpcPtInCone(TVector3 &jet, Double_t cellEtCut=0.0, Double_t rJet=0.5) ;
+ Double_t GetSum(Int_t n, Float_t *ar, Double_t cut=0.0) const ;
Double_t GetSumEmcal(Double_t cut=0.0) {return GetSum(fncell, fetcell, cut);}
Double_t GetSumTpc(Double_t cut=0.0) {return GetSum(fnchp, fppt, cut);}
void SetDebug(Int_t flag) {fDebug = flag;}
Float_t GetDebug() const {return fDebug;}
- TTree* GetTree() {return fTree;}
- TFile* GetFile() {return fFile;}
+ TTree* GetTree() const {return fTree;}
+ TFile* GetFile() const {return fFile;}
void Close();
- Bool_t IsPythiaDst();
+ Bool_t IsPythiaDst() const ;
virtual Bool_t IsFolder() const;
- virtual void Browse(TBrowser* b);
+ virtual void Browse(TBrowser* b) const ;
+
+ // service routine
+ static TList *MoveHistsToList(const char* name="ListOfHists", Bool_t putToBrowser=kTRUE);
+ static void FillH1(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t w=1.);
+ static void FillH2(TList *l=0, Int_t ind=0, Double_t x=-99999., Double_t y=-99999., Double_t w=1.);
+ static int SaveListOfHists(TList *list=0, const char* name="test", Bool_t kSingleKey=kFALSE,
+ const char* opt="RECREATE");
+ static void Sgpdge(Int_t pdgId, Int_t &gId);
- static TList *MoveHistsToList(char* name="List of Hist", Bool_t putToBrowser=kTRUE);
+ AliEMCALJetMicroDst & operator = (const AliEMCALJetMicroDst &) {
+ Fatal("operator =", "not implemented") ; return *this ; }
+
private:
Float_t fpphi[20000]; //[nchp]
Int_t fDebug; // debug flag
TFile* fFile; // filename
TTree* fTree; // Tree pointer
- TString fName; // name
TList* fListHist; //!
TString fFileName; // for convenience
Int_t fpid[20000]; //[nchp]
Float_t fppt[20000]; //[nchp]
Float_t fpeta[20000]; //[nchp]
- TH1F* hPtPart; //hist
- TH1F* hNJet; //hist
- TH1F* hPtJet; //hist
- TH2F* hEtaPhiPart; //hist
- TH2F* hEtaPhiJet; //hist
- TH1F* hNcell; //hist
- TH1F* hCellId; //hist
- TH1F* hCellEt; //hist
+ TH1F* fhPtPart; //hist
+ TH1F* fhNJet; //hist
+ TH1F* fhPtJet; //hist
+ TH2F* fhEtaPhiPart; //hist
+ TH2F* fhEtaPhiJet; //hist
+ TH1F* fhNcell; //hist
+ TH1F* fhCellId; //hist
+ TH1F* fhCellEt; //hist
- TH1F* hSumEt; //hist
- TH1F* hNgrid; //hist
- TH1F* hGridId; //hist
- TH1F* hGridEt; //hist
- TH1F* hSumEtGrForJF; //hist
+ TH1F* fhSumEt; //hist
+ TH1F* fhNgrid; //hist
+ TH1F* fhGridId; //hist
+ TH1F* fhGridEt; //hist
+ TH1F* fhSumEtGrForJF; //hist
- ClassDef(AliEMCALJetMicroDst,2) // Micro Dst for jet analysis
+ ClassDef(AliEMCALJetMicroDst,3) // Micro Dst for jet analysis
};
#endif // AliEMCALJETMICRODST_H
+
+typedef AliEMCALJetMicroDst sv; // for convinience
+
/*
What to do
1. Common info about event