AliPhysics  64f4410 (64f4410)
AliAnalysisTaskEmcalJetPerformance.h
Go to the documentation of this file.
1 #ifndef AliAnalysisTaskEmcalJetPerformance_H
2 #define AliAnalysisTaskEmcalJetPerformance_H
3 
4 /**********************************************************************************
5  * Copyright (C) 2016, Copyright Holders of the ALICE Collaboration *
6  * All rights reserved. *
7  * *
8  * Redistribution and use in source and binary forms, with or without *
9  * modification, are permitted provided that the following conditions are met: *
10  * * Redistributions of source code must retain the above copyright *
11  * notice, this list of conditions and the following disclaimer. *
12  * * Redistributions in binary form must reproduce the above copyright *
13  * notice, this list of conditions and the following disclaimer in the *
14  * documentation and/or other materials provided with the distribution. *
15  * * Neither the name of the <organization> nor the *
16  * names of its contributors may be used to endorse or promote products *
17  * derived from this software without specific prior written permission. *
18  * *
19  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND *
20  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED *
21  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE *
22  * DISCLAIMED. IN NO EVENT SHALL ALICE COLLABORATION BE LIABLE FOR ANY *
23  * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES *
24  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; *
25  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND *
26  * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT *
27  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
28  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
29  * *********************************************************************************/
30 
41 #include "THistManager.h"
42 
43 #include "AliEventCuts.h"
45 #include "AliEmcalEmbeddingQA.h"
46 
48  public:
49 
51  AliAnalysisTaskEmcalJetPerformance(const char *name) ;
53 
55  const char *ntracks = "usedefault",
56  const char *nclusters = "usedefault",
57  const char *nGenLev = "mcparticles",
58  const Double_t minTrPt = 0.15, // Minimum track pT in standard track container
59  const Double_t minClPt = 0.30, // Minimum cluster E in standard cluster container
60  const char *suffix = "");
61 
62  // Truth-level particle types
64  kUndefined = -1,
65  kPhoton = 0,
67  kProton = 2,
70  kK0L = 5,
71  kNeutron = 6,
73  kElectron = 8,
74  kMuon = 9,
75  kOther = 10
76  };
77 
78  // Detector-level particle types, based on PhysicalPrimary contributors
79  enum ParticleType {
80  kNotDefined = -1,
87  kSingleK0L = 6,
93  kMergedPi0 = 12,
98  };
99 
100  enum ClusterType {
101  kNA = -1,
102  kEMCal = 0,
103  kDCal = 1,
104  };
105 
106  void UserCreateOutputObjects() ;
107 
108  // Setters
109  void SetMinPt(Double_t d) { fMinPt = d; }
110  void SetMaxPt(Double_t d) { fMaxPt = d; }
122  void SetPlotDCal(Bool_t b) { fPlotDCal = b; }
132 
133  protected:
134  void ExecOnce() ;
136  Bool_t Run() ;
138  void RunChanged(Int_t run) ;
139 
140  // Analysis and plotting functions
141  void GenerateHistoBins() ;
142  void AllocateJetHistograms() ;
148  void FillJetHistograms() ;
149  void FillClusterHistograms() ;
151  void FillParticleCompositionClusterHistograms(const AliMCEvent* mcevent);
152  void FillParticleCompositionJetHistograms(const AliMCEvent* mcevent);
153  void SetParticleTypeLabels(TAxis* axis) ;
154  void ComputeBackground() ;
155  void DoTriggerSimulation() ;
156  void FillTriggerSimHistograms() ;
157  void FillMatchedJetHistograms() ;
158  void ComputeJetMatches(AliJetContainer* jetCont1, AliJetContainer* jetCont2, Bool_t bUseJetCont2Acceptance);
160  const AliEmcalJet* GetMatchedPartLevelJet(const AliEmcalJet* jet, Double_t detJetPt);
161 
162  // Utility functions
163  Double_t GetJetPt(const AliEmcalJet* jet, Double_t rho);
164  Double_t GetDeltaR(const AliTLorentzVector* part, Double_t etaRef, Double_t phiRef);
165  Double_t GetJetType(const AliEmcalJet* jet);
166  ContributorType GetContributorType(const AliVCluster* clus, const AliMCEvent* mcevent, Int_t label);
167  Bool_t IsHadron(const ContributorType contributor);
168 
169  // Analysis parameters
179 
180  // Plotting parameters
193 
194  // Track matching parameters (for cluster histogram plots)
197 
198  // Trigger parameters
204 
205  // Embedding parameters
217 
218  // Event selection
220  AliEventCuts fEventCuts;
223 
224  // MC options
226 
227  // Hist manager
229 
230  private:
233 
237 };
238 #endif
static AliAnalysisTaskEmcalJetPerformance * AddTaskEmcalJetPerformance(const char *ntracks="usedefault", const char *nclusters="usedefault", const char *nGenLev="mcparticles", const Double_t minTrPt=0.15, const Double_t minClPt=0.30, const char *suffix="")
double Double_t
Definition: External.C:58
TList * fEventCutList
! Output list for event cut histograms
Double_t GetJetPt(const AliEmcalJet *jet, Double_t rho)
Bool_t fUseAliEventCuts
Flag to use AliEventCuts (otherwise AliAnalysisTaskEmcal will be used)
void FillParticleCompositionJetHistograms(const AliMCEvent *mcevent)
AliMCParticleContainer * fGeneratorLevel
! generator level container
AliAnalysisTaskEmcalJetPerformance & operator=(const AliAnalysisTaskEmcalJetPerformance &)
QA Class for EMCal Embedding Framework.
Bool_t fPlotClusterHistograms
Set whether to plot cluster histograms.
Double_t GetDeltaR(const AliTLorentzVector *part, Double_t etaRef, Double_t phiRef)
Bool_t fDoJetMatchingGeometrical
Do geometrical matching between det-level and truth-level jet container.
Bool_t fComputeMBDownscaling
Set whether to compute and plot MB downscaling factors.
AliEmcalEmbeddingQA fEmbeddingQA
! QA hists for embedding (will only be added if embedding)
int Int_t
Definition: External.C:63
AliJetContainer * fMCJetContainer
!Pointer to jet container of truth-level jets
Double_t fTrackMatchingDeltaEtaMax
Maximum delta-eta to consider a track to be matched to a cluster.
Double_t fJetMatchingR
Jet matching R threshold.
void ComputeJetMatches(AliJetContainer *jetCont1, AliJetContainer *jetCont2, Bool_t bUseJetCont2Acceptance)
Double_t fTrackMatchingDeltaPhiMax
Maximum delta-phi to consider a track to be matched to a cluster.
Photon+Photon (not from same pi0) are the only contributors.
Bool_t fUseManualEventCuts
Flag to use manual event cuts.
Photon (direct or decay) is the only contributor.
Two particles from merged pi0 are the only contributors.
Bool_t fRequireMatchedJetAccepted
Flag to require matched truth jet to be accepted (aside from geometrical acceptance) ...
Bool_t fPlotJetHistograms
Set whether to enable inclusive jet histograms.
Int_t fNEoverPBins
! number of variable E/p bins
AliEventCuts fEventCuts
event selection utility
Bool_t fDoJetMatchingMCFraction
Do MC-fraction based matching using PbPb det-level, pp det-level, and pp truth-level jet containers...
Bool_t fPlotMatchedJetHistograms
Set whether to plot matched jet histograms.
void SetJetClosestCandidate(AliEmcalJet *jet1, AliEmcalJet *jet2)
Bool_t fkEMCEJE
! flag telling whether the event is "triggered" or not in "simulation"
const AliEmcalJet * GetMatchedPartLevelJet(const AliEmcalJet *jet, Double_t detJetPt)
Photon+Hadron are the only contributors, with Photon leading.
Double_t fMCJetMinMatchingPt
Min jet pT for MC jets being matched, for when container criteria is not applied. ...
Bool_t IsHadron(const ContributorType contributor)
Bool_t fComputeBackground
Set whether to enable study of background.
Double_t fMinSharedMomentumFraction
Minimum shared momentum (pp det-level track pT in combined jet) / (pp det-level track pT) ...
AliJetContainer * fDetJetContainerPPIntermediate
!Pointer to jet container of intermediate pp det-level jets, if MC-fraction matching ...
ContributorType GetContributorType(const AliVCluster *clus, const AliMCEvent *mcevent, Int_t label)
Int_t fNPtHistBins
! number of variable pt bins
Bool_t fPlotParticleCompositionHistograms
Set whether to plot jet composition histograms.
Bool_t fDoClosureTest
Set whether to do thermal model closure test.
Bool_t fDoTriggerResponse
flag whether to compute max patch response, in case of MC
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
Container class for histograms.
Definition: THistManager.h:99
Int_t fNM02HistBins
! number of variable M02 bins
Double_t fMBUpscaleFactor
! inverse of downscale factor, for MB trigger
AliJetContainer * fDetJetContainer
!Pointer to jet container of det-level jets
Bool_t fDoTriggerSimulation
Set whether to perform a simple trigger simulation.
bool Bool_t
Definition: External.C:53
Double_t fMedianDCal
! median patch energy in DCal, per event
Double_t fMedianEMCal
! median patch energy in EMCal, per event
Container for MC-true particles within the EMCAL framework.
Container for jet within the EMCAL jet framework.
Hadron+Photon are the only contributors, with Hadron leading.
void FillParticleCompositionClusterHistograms(const AliMCEvent *mcevent)
Bool_t fPlotDCal
Set whether to enable several DCal-specific histograms.
Double_t fPlotJetMatchCandThresh
Threshold for jet R to count candidates, affects plotting only.