AliPhysics  5364b50 (5364b50)
AliAnalysisTaskEmcalDijetImbalance.h
Go to the documentation of this file.
1 #ifndef AliAnalysisTaskEmcalDijetImbalance_H
2 #define AliAnalysisTaskEmcalDijetImbalance_H
3 
14 /* Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
15  * See cxx source for full Copyright notice */
16 
17 class AliPHOSGeometry;
18 
19 #include "THistManager.h"
20 
21 #include "AliEventCuts.h"
23 #include "AliEmcalEmbeddingQA.h"
24 
26  public:
27 
28  // Struct to store all relevant info of a di-jet pair
29  struct Dijet_t {
31  assJet(0), assJetPt(0), assJetPhi(0), assJetEta(0), isAccepted(kFALSE), deltaPhi(0), deltaEta(0), AJ(0), xJ(0), kTy(0) {}
32 
34 
39 
44 
46 
52 
53  void clear() {leadingHadronCutType=0; trigJet=0; trigJetPt=0; trigJetPhi=0; trigJetEta=0;
54  assJet=0; assJetPt=0; assJetPhi=0; assJetEta=0; isAccepted=kFALSE; deltaPhi=0; deltaEta=0; AJ=0; xJ=0; kTy=0;}
55  };
56 
57  enum ClusterType {
58  kNA = -1,
59  kEMCal = 0,
60  kDCal = 1,
61  kPHOS = 2
62  };
63 
65  AliAnalysisTaskEmcalDijetImbalance(const char *name) ;
67 
68  static AliAnalysisTaskEmcalDijetImbalance* AddTaskEmcalDijetImbalance(const char *ntracks = "usedefault",
69  const char *nclusters = "usedefault",
70  const Double_t deltaPhiMin = 2*TMath::Pi()/3,
71  const Bool_t doGeomMatching = kFALSE,
72  const Double_t minTrPtHardCore = 3.0,
73  const Double_t minClPtHardCore = 3.0,
74  const Double_t jetR = 0.2,
75  const Bool_t includePHOS = kTRUE,
76  const Double_t minTrPt = 0.15,
77  const Double_t minClPt = 0.30,
78  const char *suffix = "");
79 
81  void LoadBackgroundScalingHistogram(const char* path = "alien:///alice/cern.ch/user/j/jmulliga/scaleFactorEMCalLHC15o.root", const char* name1 = "hEtaPhiSFCorrection", const char* name2 = "hEtaPhiJetPtCorrection");
82 
83  // Setters
85  void SetMaxPt(Double_t d) { fMaxPt = d; }
90  void SetDoGeometricalMatching(Bool_t b, Double_t r, Double_t trackThresh, Double_t clusThresh)
100  void SetNEtaBins(Int_t n) { fNEtaBins = n; }
101  void SetNPhiBins(Int_t n) { fNPhiBins = n; }
102 
103  protected:
104  void ExecOnce() ;
106  Bool_t Run() ;
108  void RunChanged(Int_t run) ;
109 
110  // Analysis and plotting functions
111  void GenerateHistoBins() ;
118  void FindDijet(AliJetContainer* jetCont, Int_t leadingHadronCutBin);
119  void ComputeBackground() ;
120  void DoMomentumBalance(TString histname) ;
121  void DoGeometricalMatching() ;
122  void DoTriggerSimulation() ;
123  void FindMatchingDijet(AliJetContainer* jetCont) ;
126  void FillMomentumBalanceHistograms(TString histname, Double_t deltaPhi, Double_t trackPt, Double_t balancePt);
128  void FillTriggerSimHistograms() ;
129 
130  // Utility functions
136 
137  // Analysis parameters
151 
152  // Analysis configuration and plotting options
161 
162  // Plotting parameters
168 
169  // Event selection
171  AliEventCuts fEventCuts;
174 
175  // Trigger parameters
180 
181  // Embedding parameters
183 
184  // Hist manager
186 
187  private:
190 
194 };
195 #endif
AliEmcalJet * GetLeadingJet(AliJetContainer *jetCont)
Int_t fNPhiBins
Number of phi bins in DCal region (for background/correction)
Int_t fNPtHistBins
! number of variable pt bins
double Double_t
Definition: External.C:58
Bool_t fUseManualEventCuts
Flag to use manual event cuts.
Bool_t fUseAliEventCuts
Flag to use AliEventCuts (otherwise AliAnalysisTaskEmcal will be used)
Int_t fNEtaBins
Number of eta bins in DCal region (for background/correction)
TList * fEventCutList
! Output list for event cut histograms
TH2D * fGapJetScalingWeights
Histogram storing eta-phi weights scaling jets near the gap region.
Double_t fDijetLeadingHadronPt
leading hadron pT threshold for leading jet in dijet
void SetDoGeometricalMatching(Bool_t b, Double_t r, Double_t trackThresh, Double_t clusThresh)
Bool_t fkEMCEJE
! flag telling whether the event is "triggered" or not in "simulation"
Double_t fMedianEMCal
! median patch energy in EMCal, per event
Double_t fMedianDCal
! median patch energy in DCal, per event
Double_t fMinAssJetPt
subleading jet min pT in a dijet pair, for it to be accepted
TString part
use mixed event to constrain combinatorial background
Definition: InvMassFit.C:52
QA Class for EMCal Embedding Framework.
Bool_t fDoMomentumBalance
Set whether to enable momentum balance study.
void FindDijet(AliJetContainer *jetCont, Int_t leadingHadronCutBin)
Bool_t fPlotDijetImbalanceHistograms
Set whether to enable dijet imbalance histograms.
AliEventCuts fEventCuts
event selection utility
Dijet_t fDijet
! dijet candidate (per event)
int Int_t
Definition: External.C:63
Double_t GetDeltaR(AliEmcalJet *jet1, AliEmcalJet *jet2)
float Float_t
Definition: External.C:68
Double_t fClusterConstituentThreshold
constituent threshold for matching study
Bool_t fDoTriggerSimulation
Set whether to perform a simple trigger simulation.
Definition: External.C:228
Double_t fDeltaPhiMin
minimum delta phi between di-jets
Double_t fMBUpscaleFactor
! inverse of downscale factor, for MB trigger
Bool_t fComputeBackground
Set whether to enable study of background.
void LoadBackgroundScalingHistogram(const char *path="alien:///alice/cern.ch/user/j/jmulliga/scaleFactorEMCalLHC15o.root", const char *name1="hEtaPhiSFCorrection", const char *name2="hEtaPhiJetPtCorrection")
Double_t GetJetPt(AliJetContainer *jetCont, AliEmcalJet *jet)
void FillMomentumBalanceHistograms(TString histname, Double_t deltaPhi, Double_t trackPt, Double_t balancePt)
Bool_t fComputeMBDownscaling
Set whether to compute and plot MB downscaling factors.
AliAnalysisTaskEmcalDijetImbalance & operator=(const AliAnalysisTaskEmcalDijetImbalance &)
Bool_t fLoadBackgroundScalingWeights
Flag to load eta-phi weights for full-jet background scale factors.
Base task in the EMCAL jet framework.
Bool_t fDoGeometricalMatching
Set whether to enable constituent study with geometrical matching.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
Container class for histograms.
Definition: THistManager.h:99
Dijet_t fMatchingDijet
! low-threshold matching dijet, for matching study
TH2D * fBackgroundScalingWeights
Histogram storing eta-phi weights for full-jet background scale factors.
AliEmcalEmbeddingQA fEmbeddingQA
! QA hists for embedding (will only be added if embedding)
bool Bool_t
Definition: External.C:53
Double_t fTrackConstituentThreshold
constituent threshold for matching study
static AliAnalysisTaskEmcalDijetImbalance * AddTaskEmcalDijetImbalance(const char *ntracks="usedefault", const char *nclusters="usedefault", const Double_t deltaPhiMin=2 *TMath::Pi()/3, const Bool_t doGeomMatching=kFALSE, const Double_t minTrPtHardCore=3.0, const Double_t minClPtHardCore=3.0, const Double_t jetR=0.2, const Bool_t includePHOS=kTRUE, const Double_t minTrPt=0.15, const Double_t minClPt=0.30, const char *suffix="")
Double_t fMinTrigJetPt
leading jet min pT in a dijet pair
Container for jet within the EMCAL jet framework.
Bool_t fPlotDijetCandHistograms
Set whether to enable dijet pair histograms.