AliPhysics  d565ceb (d565ceb)
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; }
121 
122  protected:
123  void ExecOnce() ;
125  Bool_t Run() ;
127  void RunChanged(Int_t run) ;
128 
129  // Analysis and plotting functions
130  void GenerateHistoBins() ;
131  void AllocateJetHistograms() ;
137  void FillJetHistograms() ;
138  void FillClusterHistograms() ;
140  void FillParticleCompositionClusterHistograms(const AliMCEvent* mcevent);
141  void FillParticleCompositionJetHistograms(const AliMCEvent* mcevent);
142  void SetParticleTypeLabels(TAxis* axis) ;
143  void ComputeBackground() ;
144  void DoTriggerSimulation() ;
145  void FillTriggerSimHistograms() ;
146  void FillMatchedJetHistograms() ;
147 
148  // Utility functions
149  Double_t GetJetPt(const AliEmcalJet* jet, Double_t rho);
150  Double_t GetDeltaR(const AliTLorentzVector* part, Double_t etaRef, Double_t phiRef);
151  Double_t GetJetType(const AliEmcalJet* jet);
152  ContributorType GetContributorType(const AliVCluster* clus, const AliMCEvent* mcevent, Int_t label);
153  Bool_t IsHadron(const ContributorType contributor);
154 
155  // Analysis parameters
163 
164  // Plotting parameters
176 
177  // Track matching parameters (for cluster histogram plots)
180 
181  // Trigger parameters
186 
187  // Embedding parameters
189 
190  // Event selection
192  AliEventCuts fEventCuts;
195 
196  // MC options
198 
199  // Hist manager
201 
202  private:
205 
209 };
210 #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
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"
Photon+Hadron are the only contributors, with Photon leading.
Bool_t IsHadron(const ContributorType contributor)
Bool_t fComputeBackground
Set whether to enable study of background.
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.
Hadron+Photon are the only contributors, with Hadron leading.
void FillParticleCompositionClusterHistograms(const AliMCEvent *mcevent)