AliPhysics  ba8894a (ba8894a)
AliAnalysisTaskJetCoreEmcal.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKJETCOREEMCAL_H
2 #define ALIANALYSISTASKJETCOREEMCAL_H
3 
4 //
5 // Implementation of jet core task in Emcal framework
6 // Basic functionality copied from AliAnalysisTaskEmcalJetSample
7 //
8 
9 class TH1F;
10 class TH1I;
11 class TH2F;
12 class TH3F;
13 class THnSparse;
14 class TRandom3;
15 
16 #include "AliEventCuts.h"
18 #include "THistManager.h"
19 
20 
22  public:
23 
25  AliAnalysisTaskJetCoreEmcal(const char *name) ;
27 
29  void Terminate(Option_t *option) ;
30 
31 
32  // Setters
33  virtual void SetJetShapeType(Int_t type){fJetShapeType=type;}
34  virtual void SetCentrality(Float_t cmin, Float_t cmax){fCentMin=cmin; fCentMax=cmax;}
35  virtual void SetTTLowRef(Float_t ttlow){fTTLowRef=ttlow;}
36  virtual void SetTTUpRef(Float_t ttup){fTTUpRef=ttup;}
37  virtual void SetTTLowSig(Float_t ttlows){fTTLowSig=ttlows;}
38  virtual void SetTTUpSig(Float_t ttups){fTTUpSig=ttups;}
39  virtual void SetNRPBins(Float_t nrpb){fNRPBins=nrpb;}
40  virtual void SetSignalFraction(Float_t sfrac){fFrac=sfrac;}
41  virtual void SetJetEtaMin(Float_t eta){fJetEtaMin=eta;}
42  virtual void SetJetEtaMax(Float_t eta){fJetEtaMax=eta;}
43  virtual void SetJetHadronDeltaPhi(Float_t delta){fJetHadronDeltaPhi=delta;}
45  virtual void SetMinEmbJetPt(Float_t min){fMinEmbJetPt=min;}
46  virtual void SetJetContName(TString cont){fJetContName=cont;}
47  virtual void SetJetContTrueName(TString cont){fJetContTrueName=cont;}
48  virtual void SetJetContPartName(TString cont){fJetContPartName=cont;}
54  virtual void SetPtHardBin(Int_t bin){fPtHardBin=bin;}
56 
57 // static AliAnalysisTaskJetCoreEmcal* AddTaskJetCoreEmcal(
58 // const char *ntracks = "usedefault",
59 // const char *nclusters = "usedefault",
60 // const char* ncells = "usedefault",
61 // const char *suffix = "");
62 
63  enum JetShapeType {
64  kMCTrue = 0, // generated jets only
65  kTrueDet =1, // detector and generated jets
66  kData = 2, // raw data
68 // kDetEmb = 3, //detector embedded jets
69 // kPythiaDef = 5,
70 // kDetEmbPartPythia=6,
71 // kGenOnTheFly = 7
72  };
73 
74  protected:
75  void ExecOnce() ;
77  Bool_t Run() ;
78 
79  void AllocateJetHistograms() ;
82  void AllocateCellHistograms() ;
84 
85  void DoJetLoop() ;
86  void DoTrackLoop() ;
87  void DoClusterLoop() ;
88  void DoCellLoop() ;
89  void DoJetCoreLoop() ;
90  void DoMatchingLoop()
91  ;
92 
93  Int_t SelectTrigger(TList *list,Double_t minT,Double_t maxT,Int_t &number);
96 
98 
99  // flags and selection
100  AliEventCuts fEventCuts;
125  //
126  TRandom3 *fRandom;
129  //histograms to fill
131  // recoil jet info contained in THnSparse
132  THnSparse *fHJetSpec;
133  // recoil histograms
143  // embedding histograms
144  // inclusive jets
150  TH1F *fhPtDet;
152  TH1F *fhResidual;
156  //recoil jets
171  // embedding trees
174 
175  private:
176  AliAnalysisTaskJetCoreEmcal(const AliAnalysisTaskJetCoreEmcal&) ; // not implemented
178 
180  ClassDef(AliAnalysisTaskJetCoreEmcal, 7);
182 };
183 #endif
virtual void SetFillTrackHistograms(Bool_t b)
Float_t fJetHadronDeltaPhi
max angle from pi (set <0 for no selection)
Float_t fTTLowRef
minimum reference trigger track pt
double Double_t
Definition: External.C:58
Definition: External.C:260
Float_t fFrac
fraction of events that are used to fill signal recoil jet population
Bool_t fFillRecoilTHnSparse
switch to fill recoil THnSparse for main analysis
Definition: External.C:236
Int_t SelectTrigger(TList *list, Double_t minT, Double_t maxT, Int_t &number)
virtual void SetTTLowRef(Float_t ttlow)
Float_t fTTUpRef
maximum reference trigger track pt
virtual void SetTTUpRef(Float_t ttup)
Bool_t fFillInclusiveTree
switch to fill embedding tree with inclusive jet info
virtual void SetFillRecoilTree(Bool_t b)
virtual void SetCentrality(Float_t cmin, Float_t cmax)
virtual void SetFillInclusiveTree(Bool_t b)
Float_t fMinEmbJetPt
min corrected jet pt to use in embedding
TString fJetContName
Base level jet container name.
Float_t fMinFractionSharedPt
min fraction of pt between hybrid / detector jets
TString fJetContTrueName
True pp (detector) level jet container name.
virtual void SetSignalFraction(Float_t sfrac)
virtual void SetJetContTrueName(TString cont)
int Int_t
Definition: External.C:63
Float_t fTTUpSig
maximum signal trigger track pt
Definition: External.C:204
float Float_t
Definition: External.C:68
Float_t fTTLowSig
minimum signal trigger track pt
virtual void SetJetContName(TString cont)
Definition: External.C:212
TString fJetContPartName
Particle(MC) level jet container name.
virtual void SetMinFractionSharedPt(Float_t min)
Float_t fCentMax
maximum centrality
virtual void SetJetEtaMin(Float_t eta)
Bool_t fFillJetHistograms
switch to fill jet histograms
virtual void SetJetHadronDeltaPhi(Float_t delta)
virtual void SetMinEmbJetPt(Float_t min)
Double_t RelativePhi(Double_t mphi, Double_t vphi)
virtual void SetTTLowSig(Float_t ttlows)
virtual void SetJetEtaMax(Float_t eta)
Int_t fPtHardBin
pt hard bin if running embedding
THistManager fHistManager
Histogram manager.
virtual void SetRejectionFactorInclusiveJets(Int_t f)
virtual void SetNRPBins(Float_t nrpb)
Base task in the EMCAL jet framework.
Container class for histograms.
Definition: THistManager.h:99
virtual void SetFillRecoilTHnSparse(Bool_t b)
const char Option_t
Definition: External.C:48
bool Bool_t
Definition: External.C:53
Bool_t fFillRecoilTree
switch to fill embedding tree with recoil jet info
virtual void SetJetShapeType(Int_t type)
virtual void SetFillJetHistograms(Bool_t b)
AliAnalysisTaskJetCoreEmcal & operator=(const AliAnalysisTaskJetCoreEmcal &)
Float_t fCentMin
minimum centrality
virtual void SetTTUpSig(Float_t ttups)
Bool_t fFillTrackHistograms
switch to fill track histograms
Int_t fRejectionFactorInclusiveJets
factor to reject inclusive jets, to reduce size of ttree
virtual void SetJetContPartName(TString cont)