AliPhysics  d497547 (d497547)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskEmcalJetSubstructureTree.h
Go to the documentation of this file.
1 /************************************************************************************
2  * Copyright (C) 2017, Copyright Holders of the ALICE Collaboration *
3  * All rights reserved. *
4  * *
5  * Redistribution and use in source and binary forms, with or without *
6  * modification, are permitted provided that the following conditions are met: *
7  * * Redistributions of source code must retain the above copyright *
8  * notice, this list of conditions and the following disclaimer. *
9  * * Redistributions in binary form must reproduce the above copyright *
10  * notice, this list of conditions and the following disclaimer in the *
11  * documentation and/or other materials provided with the distribution. *
12  * * Neither the name of the <organization> nor the *
13  * names of its contributors may be used to endorse or promote products *
14  * derived from this software without specific prior written permission. *
15  * *
16  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND *
17  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED *
18  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE *
19  * DISCLAIMED. IN NO EVENT SHALL ALICE COLLABORATION BE LIABLE FOR ANY *
20  * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES *
21  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; *
22  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND *
23  * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT *
24  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
25  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
26  ************************************************************************************/
27 #ifndef ALIANALYSISTASKEMCALJETSUBSTRUCTURETREE_H
28 #define ALIANALYSISTASKEMCALJETSUBSTRUCTURETREE_H
29 
31 #include <exception>
32 #include <TString.h>
33 #include <fastjet/PseudoJet.hh>
34 #include <fastjet/JetDefinition.hh>
35 
36 class TTree;
38 class AliEmcalJet;
40 class AliTrackContainer;
41 
42 namespace EmcalTriggerJets {
43 
52 };
53 
62 };
63 
75 };
76 
85  fastjet::JetAlgorithm fRecluserAlgo;
86 };
87 
91 };
92 
96 };
97 
105 public:
106  class ReclusterizerException : public std::exception {
107  public:
108  ReclusterizerException() : std::exception() {}
109  virtual ~ReclusterizerException() throw() {}
110 
111  const char *what() const throw() { return "Error in reclusterizing in fastjet"; }
112  };
114  kCAAlgo = 0,
115  kKTAlgo = 1,
117  };
119  kTRadius = 0,
120  kTWeight = 1,
129  kTNEFRec = 10,
130  kTNEFSim = 11,
131  kTMassRec = 12,
132  kTMassSim = 13,
134  kTZgTrue = 15,
136  kTRgTrue = 17,
138  kTMgTrue = 19,
140  kTPtgTrue = 21,
148  kTPtDTrue = 29,
154  kTNVar = 35
155  };
157  AliAnalysisTaskEmcalJetSubstructureTree(const char *name);
159 
160  void SetTriggerBits(UInt_t triggersel) { fTriggerSelectionBits = triggersel; }
161  void SetTriggerString(TString triggerstring) { fTriggerSelectionString = triggerstring; }
162 
163  void SetNameRhoPtRec(const char *rhoname) { fRhoNamePtRec = rhoname; }
164  void SetNameRhoMassRec(const char *rhoname) { fRhoNameMassRec = rhoname; }
165  void SetNameRhoPtSim(const char *rhoname) { fRhoNamePtSim = rhoname; }
166  void SetNameRhoMassSim(const char *rhoname) { fRhoNameMassSim = rhoname; }
167 
168  void SetSoftdropDefiniion(Double_t zcut, Double_t betacut, Reclusterizer_t reclusterizer) {
169  fSDZCut = zcut;
170  fSDBetaCut = betacut;
171  fReclusterizer = reclusterizer;
172  }
173 
175 
176 protected:
177  virtual void UserCreateOutputObjects();
178  virtual bool Run();
179 
180  AliJetSubstructureData MakeJetSubstructure(const AliEmcalJet &jet, double jetradius, const AliParticleContainer *tracks, const AliClusterContainer *clusters, const AliJetSubstructureSettings &settings) const;
181 
182  AliSoftDropParameters MakeSoftDropParameters(const fastjet::PseudoJet &jet, const AliSoftdropDefinition &cut) const;
183 
184  AliNSubjettinessParameters MakeNsubjettinessParameters(const fastjet::PseudoJet &jet, const AliNSubjettinessDefinition &cut) const;
185 
186  Double_t MakeAngularity(const AliEmcalJet &jet, const AliParticleContainer *tracks, const AliClusterContainer *clusters) const;
187 
188  Double_t MakePtD(const AliEmcalJet &jet, const AliParticleContainer *const particles, const AliClusterContainer *const clusters) const;
189 
190  void FillTree(double r, double weight, const AliEmcalJet *datajet, const AliEmcalJet *mcjet, AliSoftDropParameters *dataSoftdrop, AliSoftDropParameters *mcsoftdrop, AliNSubjettinessParameters *dataSubjettiness, AliNSubjettinessParameters *mcSubjettiness, Double_t *angularity, Double_t *ptd, Double_t *rhoparameters);
191 
192 private:
195 
199 
202 
207 
211 };
212 
213 } /* namespace EmcalTriggerJets */
214 
215 #endif /* ALIANALYSISTASKEMCALJETSUBSTRUCTURETREE_H */
TString fRhoNameMassSim
Name of the rho parameter coming from the rho task for mc-level jet (mass)
double Double_t
Definition: External.C:58
Double_t MakePtD(const AliEmcalJet &jet, const AliParticleContainer *const particles, const AliClusterContainer *const clusters) const
void SetSoftdropDefiniion(Double_t zcut, Double_t betacut, Reclusterizer_t reclusterizer)
Container with name, TClonesArray and cuts for particles.
AliJetSubstructureData MakeJetSubstructure(const AliEmcalJet &jet, double jetradius, const AliParticleContainer *tracks, const AliClusterContainer *clusters, const AliJetSubstructureSettings &settings) const
TString fRhoNamePtSim
Name of the rho parameter coming from the rho task for mc-level jet (pt)
Structure for results from the soft drop algorithm.
virtual bool Run()
Run function. This is the core function of the analysis and contains the user code. Therefore users have to implement this function.
fastjet::JetAlgorithm fRecluserAlgo
Reclusterization algorithm.
Container for particles within the EMCAL framework.
int Int_t
Definition: External.C:63
unsigned int UInt_t
Definition: External.C:33
AliNSubjettinessParameters MakeNsubjettinessParameters(const fastjet::PseudoJet &jet, const AliNSubjettinessDefinition &cut) const
Definition for the algorithm obtaining the softdrop parameters.
Double_t MakeAngularity(const AliEmcalJet &jet, const AliParticleContainer *tracks, const AliClusterContainer *clusters) const
TString fRhoNamePtRec
Name of the rho parameter coming from the rho task for rec jet (pt)
Bool_t isMC
Base task in the EMCAL jet framework.
void FillTree(double r, double weight, const AliEmcalJet *datajet, const AliEmcalJet *mcjet, AliSoftDropParameters *dataSoftdrop, AliSoftDropParameters *mcsoftdrop, AliNSubjettinessParameters *dataSubjettiness, AliNSubjettinessParameters *mcSubjettiness, Double_t *angularity, Double_t *ptd, Double_t *rhoparameters)
static AliAnalysisTaskEmcalJetSubstructureTree * AddEmcalJetSubstructureTreeMaker(Bool_t isMC, Bool_t isData, Double_t jetradius, const char *name)
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:44
AliSoftDropParameters MakeSoftDropParameters(const fastjet::PseudoJet &jet, const AliSoftdropDefinition &cut) const
bool Bool_t
Definition: External.C:53
TString fRhoNameMassRec
Name of the rho parameter coming from the rho task for rec jet (mass)
Container structure for EMCAL clusters.