AliPhysics  8d00e07 (8d00e07)
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 SetMaxPt(Double_t d) { fMaxPt = d; }
124  void SetPlotDCal(Bool_t b) { fPlotDCal = b; }
125 
126  protected:
127  void ExecOnce() ;
129  Bool_t Run() ;
131  void RunChanged(Int_t run) ;
132 
133  // Analysis and plotting functions
134  void GenerateHistoBins() ;
135  void AllocateJetHistograms() ;
141  void FillJetHistograms() ;
142  void FillClusterHistograms() ;
144  void FillParticleCompositionClusterHistograms(const AliMCEvent* mcevent);
145  void FillParticleCompositionJetHistograms(const AliMCEvent* mcevent);
146  void SetParticleTypeLabels(TAxis* axis) ;
147  void ComputeBackground() ;
148  void DoTriggerSimulation() ;
149  void FillTriggerSimHistograms() ;
150  void FillMatchedJetHistograms() ;
151 
152  // Utility functions
153  Double_t GetJetPt(const AliEmcalJet* jet, Double_t rho);
154  Double_t GetDeltaR(const AliTLorentzVector* part, Double_t etaRef, Double_t phiRef);
155  Double_t GetJetType(const AliEmcalJet* jet);
156  ContributorType GetContributorType(const AliVCluster* clus, const AliMCEvent* mcevent, Int_t label);
157  Bool_t IsHadron(const ContributorType contributor);
158 
159  // Jet matching functions (for embedding)
160  const AliEmcalJet* GetMatchedPartLevelJet(const AliJetContainer * jets, const AliEmcalJet * jet, const std::string & histName);
161 
162  // Analysis parameters
171 
172  // Plotting parameters
184 
185  // Track matching parameters (for cluster histogram plots)
188 
189  // Trigger parameters
194 
195  // Embedding parameters
201 
202  // Event selection
204  AliEventCuts fEventCuts;
207 
208  // MC options
210 
211  // Hist manager
213 
214  private:
217 
221 };
222 #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 &)
TString part
use mixed event to constrain combinatorial background
Definition: InvMassFit.C:52
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 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.
float Float_t
Definition: External.C:68
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 fPlotJetHistograms
Set whether to enable inclusive jet histograms.
Int_t fNEoverPBins
! number of variable E/p bins
AliEventCuts fEventCuts
event selection utility
Bool_t fPlotMatchedJetHistograms
Set whether to plot matched jet histograms (must run ResponseMaker first)
Bool_t fkEMCEJE
! flag telling whether the event is "triggered" or not in "simulation"
const AliEmcalJet * GetMatchedPartLevelJet(const AliJetContainer *jets, const AliEmcalJet *jet, const std::string &histName)
Photon+Hadron are the only contributors, with Photon leading.
Bool_t fUseResponseMaker
Flag to use Response Maker rather than JetTagger.
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) ...
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.
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
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 fMaxMatchedJetDistance
Maximum distance between two matched jets.