AliPhysics  b555aef (b555aef)
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 
69  void LoadBackgroundScalingHistogram(const char* path = "alien:///alice/cern.ch/user/j/jmulliga/scaleFactorEMCalLHC15o.root", const char* name1 = "hEtaPhiSFCorrection", const char* name2 = "hEtaPhiJetPtCorrection");
70 
71  // Setters
73  void SetMaxPt(Double_t d) { fMaxPt = d; }
78  void SetDoGeometricalMatching(Bool_t b, Double_t r, Double_t trackThresh, Double_t clusThresh)
88  void SetNEtaBins(Int_t n) { fNEtaBins = n; }
89  void SetNPhiBins(Int_t n) { fNPhiBins = n; }
90 
91  protected:
92  void ExecOnce() ;
94  Bool_t Run() ;
96  void RunChanged(Int_t run) ;
97 
98  // Analysis and plotting functions
99  void GenerateHistoBins() ;
106  void FindDijet(AliJetContainer* jetCont, Int_t leadingHadronCutBin);
107  void ComputeBackground() ;
108  void DoMomentumBalance(TString histname) ;
109  void DoGeometricalMatching() ;
110  void DoTriggerSimulation() ;
111  void FindMatchingDijet(AliJetContainer* jetCont) ;
114  void FillMomentumBalanceHistograms(TString histname, Double_t deltaPhi, Double_t trackPt, Double_t balancePt);
116  void FillTriggerSimHistograms() ;
117 
118  // Utility functions
122  Double_t GetDeltaR(AliTLorentzVector* part, Double_t etaRef, Double_t phiRef);
124 
125  // Analysis parameters
139 
140  // Analysis configuration and plotting options
149 
150  // Plotting parameters
156 
157  // Event selection
159  AliEventCuts fEventCuts;
162 
163  // Trigger parameters
168 
169  // Embedding parameters
171 
172  // Hist manager
174 
175  private:
178 
182 };
183 #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
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
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.