AliPhysics  4646b6b (4646b6b)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskSubJetFraction.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKSUBJETFRACTION_H
2 #define ALIANALYSISTASKSUBJETFRACTION_H
3 
4 class TH1;
5 class TH2;
6 class TH3;
7 class TH3F;
8 class TTree;
9 class THnSparse;
10 class TClonesArray;
11 class TArrayI;
12 class AliAnalysisManager;
13 class AliJetContainer;
14 class AliEmcalJetFinder;
15 class AliFJWrapper;
16 
17 #include "AliAnalysisTaskEmcalJet.h"
18 #include "AliFJWrapper.h"
19 #include "AliClusterContainer.h"
20 const Int_t nVar = 28;
21 class AliAnalysisTaskSubJetFraction : public AliAnalysisTaskEmcalJet {
22  public:
23 
24  enum JetShapeType {
25  kTrue = 0, // generated jets only
26  kTrueDet =1, // detector and generated jets
27  kData = 2, // raw data
28  kDetEmbPart = 3, //detector embedded jets
29  kSim = 4,
30  kGenOnTheFly=5 //fast truth level MC
31  };
32  enum JetShapeSub {
33  kNoSub = 0,
34  kConstSub = 1,
35  kDerivSub = 2
36  };
37  enum JetSelectionType {
38  kInclusive = 0,
39  kRecoil = 1
40  };
41  enum DerivSubtrOrder {
42  kSecondOrder = 0,
43  kFirstOrder = 1
44  };
45 
46  AliAnalysisTaskSubJetFraction();
47  AliAnalysisTaskSubJetFraction(const char *name);
48  virtual ~AliAnalysisTaskSubJetFraction();
49 
50  void UserCreateOutputObjects();
51  void Terminate(Option_t *option);
52 
53  //Setters
54  void SetJetContainer(Int_t c) { fContainer = c ;}
55  void SetMinFractionShared(Double_t f) { fMinFractionShared = f ;}
56  void SetJetShapeType(JetShapeType t) { fJetShapeType = t ;}
57  void SetJetShapeSub(JetShapeSub t) { fJetShapeSub = t ;}
58  void SetJetSelection(JetSelectionType t) { fJetSelection = t ;}
59  void SetJetPtThreshold(Float_t f) { fPtThreshold = f ;}
60  void SetJetRadius(Double_t JetRadius) {fJetRadius=JetRadius;}
61  void SetRMatching(Float_t f) { fRMatching = f ;}
62  void SetPtTriggerSelections(Float_t PtMin, Float_t PtMax) { fPtMinTriggerHadron = PtMin; fPtMaxTriggerHadron = PtMax;}
63  void SetAngularWindowRecoilJet (Float_t Window) {fRecoilAngularWindow = Window; }
64  Float_t GetMinPtTriggerSelection() {return fPtMinTriggerHadron;}
65  Float_t GetMaxPtTriggerSelection() {return fPtMaxTriggerHadron;}
66  void SetCentralitySelectionOn(Bool_t t) { fCentSelectOn = t;}
67  void SetMinCentrality(Float_t t) { fCentMin = t ;}
68  void SetMaxCentrality(Float_t t) { fCentMax = t ;}
69  void SetSemigoodCorrect(Int_t yesno) {fSemigoodCorrect=yesno;}
70  void SetHolePos(Float_t poshole) { fHolePos = poshole;}
71  void SetHoleWidth(Float_t holewidth) { fHoleWidth = holewidth;}
72  void SetSubJetAlgorithm(Int_t SubJetAlgorithm) {fSubJetAlgorithm=SubJetAlgorithm;}
73  void SetSubJetRadius(Float_t SubJetRadius) {fSubJetRadius=SubJetRadius;}
74  void SetSubJetMinPt(Float_t SubJetMinPt) {fSubJetMinPt=SubJetMinPt;}
75  void SetRMatched(Double_t RMatched) {fRMatched=RMatched;}
76  void SetSharedFractionPtMin(Double_t SharedFractionPtMin) {fSharedFractionPtMin=SharedFractionPtMin;}
77  void SetDerivativeSubtractionOrder(Int_t Order) {fDerivSubtrOrder = Order;}
78  void SetFullTree(Bool_t FullTree) {fFullTree = FullTree;}
79  void SetBetaSD(Double_t BetaSD) {fBeta_SD = BetaSD;}
80  void SetZCut(Double_t ZCut) {fZCut = ZCut;}
81  void SetReclusteringAlgorithm(Int_t ReclusteringAlgorithm) {fReclusteringAlgorithm = ReclusteringAlgorithm;}
82  void SetSoftDropOn(Int_t SoftDropOn) {fSoftDropOn = SoftDropOn;}
83  Int_t GetSoftDropOn() {return fSoftDropOn;}
84 
85  void SetNsubUnNormMeasure( Bool_t NsubMeasure) {fNsubMeasure= NsubMeasure;}
86 
87  protected:
88  Bool_t RetrieveEventObjects();
89  Bool_t Run();
90  Bool_t FillHistograms();
91 
92  Int_t SelectTriggerHadron(Float_t PtMin, Float_t PtMax);
93  Double_t RelativePhiEventPlane(Double_t EventPlane, Double_t Phi);
94  Double_t RelativePhi(Double_t Phi1, Double_t Phi2);
95  Double_t Angularity(AliEmcalJet *Jet, Int_t JetContNb);
96  Double_t PTD(AliEmcalJet *Jet, Int_t JetContNb);
97  AliEmcalJetFinder *Recluster(AliEmcalJet *Jet, Int_t JetContNb, Double_t SubJetRadius, Double_t SubJetMinPt, Int_t Algorithm, const char* Name);
98  Double_t SubJetOrdering(AliEmcalJet *Jet, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t Type, Bool_t Index);
99  Double_t SubJetFraction(AliEmcalJet *Jet, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t Type, Bool_t Add, Bool_t Loss);
100  Double_t NSubJettiness(AliEmcalJet *Jet, Int_t JetContNb, Double_t JetRadius, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t A, Int_t B);
101  Double_t FjNSubJettiness(AliEmcalJet *Jet, Int_t JetContNb, Int_t N, Int_t Algorithm, Double_t Beta, Int_t Option=0, Double_t Beta_SD=0, Double_t ZCut=0.1);
102  Double_t FjNSubJettinessFastJet(std::pair<fastjet::PseudoJet,fastjet::ClusterSequence *> Jet_ClusterSequence, Int_t JetContNb,Int_t N, Int_t Algorithm, Double_t Beta, Int_t Option=0, Double_t Beta_SD=0, Double_t ZCut=0.1);
103  std::pair<fastjet::PseudoJet,fastjet::ClusterSequence *> ModifyJet(AliEmcalJet* Jet, Int_t JetContNb, TString Modification);
104  std::vector<fastjet::PseudoJet> RandomiseTracks(AliEmcalJet *Jet,std::vector<fastjet::PseudoJet> fInputVectors);
105  std::vector<fastjet::PseudoJet> AddExtraProng(std::vector<fastjet::PseudoJet> fInputVectors, Double_t Distance, Double_t PtFrac);
106 
107  Int_t fContainer; // jets to be analyzed 0 for Base, 1 for subtracted.
108  Float_t fMinFractionShared; // only fill histos for jets if shared fraction larger than X
109  JetShapeType fJetShapeType; // jet type to be used
110  JetShapeSub fJetShapeSub; // jet subtraction to be used
111  JetSelectionType fJetSelection; // Jet selection: inclusive/recoil jet
112  Double_t fShapesVar[nVar]; // jet shapes used for the tagging
113  Float_t fPtThreshold;
114  Float_t fRMatching;
115 
116  Float_t fPtMinTriggerHadron; //min - max pT for trigger particle in case of recoil jet
117  Float_t fPtMaxTriggerHadron;
118  Float_t fRecoilAngularWindow; //angular window for btb recoil analysis
119  Int_t fSemigoodCorrect; //if==1 we run over semigood runs
120  Float_t fHolePos; //position in radians of the bad TPC sector
121  Float_t fHoleWidth; //width of the hole in radians
122  TRandom3 *fRandom;
123  Bool_t fCentSelectOn; // switch on/off centrality selection
124  Float_t fCentMin; // min centrality value
125  Float_t fCentMax; // max centrality value
126  Int_t fSubJetAlgorithm;
127  Float_t fSubJetRadius;
128  Float_t fSubJetMinPt;
129  Double_t fJetRadius;
130  Double_t fRMatched;
131  Double_t fSharedFractionPtMin;
132  Double_t Background_Median;
133  Double_t Background_Fluc;
134  Int_t fDerivSubtrOrder;
135  Bool_t fFullTree;
136  Double_t fBeta_SD;
137  Double_t fZCut;
138  Int_t fReclusteringAlgorithm;
139  Int_t fSoftDropOn;
140 
141  Bool_t fNsubMeasure;
142 
143  TH1F *fhPtTriggerHadron;
144  TH1F *fhJetPt;
145  TH1F *fhJetPt_1;
146  TH1F *fhJetPt_2;
147  TH1F *fhJetPhi;
148  TH1F *fhJetPhi_1;
149  TH1F *fhJetPhi_2;
150  TH1F *fhJetEta;
151  TH1F *fhJetEta_1;
152  TH1F *fhJetEta_2;
153  TH1F *fhJetMass;
154  TH1F *fhJetMass_1;
155  TH1F *fhJetMass_2;
156  TH1F *fhJetRadius;
157  TH1F *fhJetRadius_1;
158  TH1F *fhJetRadius_2;
159  TH1F *fhJetCounter;
160  TH1F *fhJetCounter_1;
161  TH1F *fhJetCounter_2;
162  TH1F *fhNumberOfJetTracks;
163  TH1F *fhNumberOfJetTracks_1;
164  TH1F *fhNumberOfJetTracks_2;
165  TH1F *fhSubJetPt;
166  TH1F *fhSubJetPt_1;
167  TH1F *fhSubJetPt_2;
168  TH1F *fhSubJetMass;
169  TH1F *fhSubJetMass_1;
170  TH1F *fhSubJetMass_2;
171  TH1F *fhSubJetRadius;
172  TH1F *fhSubJetRadius_1;
173  TH1F *fhSubJetRadius_2;
174  TH1F *fhSubJetCounter;
175  TH1F *fhSubJetCounter_1;
176  TH1F *fhSubJetCounter_2;
177  TH1F *fhNumberOfSubJetTracks;
178  TH1F *fhNumberOfSubJetTracks_1;
179  TH1F *fhNumberOfSubJetTracks_2;
180  TH1F *fhEventCounter;
181  TH1F *fhEventCounter_1;
182  TH1F *fhEventCounter_2;
183  TH1F *fhPhiTriggerHadronJet;
184  TH1F *fhPhiTriggerHadronEventPlane;
185  TH1F *fhPhiTriggerHadronEventPlaneTPC;
186  TH2F *fh2PtTriggerHadronJet;
187  TH2F *fh2PtRatio;
188  TH2D *fhSubJettiness1CheckRatio_FJ_AKT;
189  TH2D *fhSubJettiness1CheckRatio_FJ_KT;
190  TH2D *fhSubJettiness1CheckRatio_FJ_CA;
191  TH2D *fhSubJettiness1CheckRatio_FJ_WTA_KT;
192  TH2D *fhSubJettiness1CheckRatio_FJ_WTA_CA;
193  TH2D *fhSubJettiness1CheckRatio_FJ_OP_AKT;
194  TH2D *fhSubJettiness1CheckRatio_FJ_OP_KT;
195  TH2D *fhSubJettiness1CheckRatio_FJ_OP_CA;
196  TH2D *fhSubJettiness1CheckRatio_FJ_OP_WTA_KT;
197  TH2D *fhSubJettiness1CheckRatio_FJ_OP_WTA_CA;
198  TH2D *fhSubJettiness1CheckRatio_FJ_MIN;
199  TH2D *fhSubJettiness2CheckRatio_FJ_AKT;
200  TH2D *fhSubJettiness2CheckRatio_FJ_KT;
201  TH2D *fhSubJettiness2CheckRatio_FJ_CA;
202  TH2D *fhSubJettiness2CheckRatio_FJ_WTA_KT;
203  TH2D *fhSubJettiness2CheckRatio_FJ_WTA_CA;
204  TH2D *fhSubJettiness2CheckRatio_FJ_OP_AKT;
205  TH2D *fhSubJettiness2CheckRatio_FJ_OP_KT;
206  TH2D *fhSubJettiness2CheckRatio_FJ_OP_CA;
207  TH2D *fhSubJettiness2CheckRatio_FJ_OP_WTA_KT;
208  TH2D *fhSubJettiness2CheckRatio_FJ_OP_WTA_CA;
209  TH2D *fhSubJettiness2CheckRatio_FJ_MIN;
210  TH2D *fhSubJettiness2to1CheckRatio_FJ_AKT;
211  TH2D *fhSubJettiness2to1CheckRatio_FJ_KT;
212  TH2D *fhSubJettiness2to1CheckRatio_FJ_CA;
213  TH2D *fhSubJettiness2to1CheckRatio_FJ_WTA_KT;
214  TH2D *fhSubJettiness2to1CheckRatio_FJ_WTA_CA;
215  TH2D *fhSubJettiness2to1CheckRatio_FJ_OP_AKT;
216  TH2D *fhSubJettiness2to1CheckRatio_FJ_OP_KT;
217  TH2D *fhSubJettiness2to1CheckRatio_FJ_OP_CA;
218  TH2D *fhSubJettiness2to1CheckRatio_FJ_OP_WTA_KT;
219  TH2D *fhSubJettiness2to1CheckRatio_FJ_OP_WTA_CA;
220  TH2D *fhSubJettiness2to1CheckRatio_FJ_MIN;
221  TH1D *fhSubJettiness1_FJ_AKT;
222  TH1D *fhSubJettiness1_FJ_KT;
223  TH1D *fhSubJettiness1_FJ_CA;
224  TH1D *fhSubJettiness1_FJ_WTA_KT;
225  TH1D *fhSubJettiness1_FJ_WTA_CA;
226  TH1D *fhSubJettiness1_FJ_OP_AKT;
227  TH1D *fhSubJettiness1_FJ_OP_KT;
228  TH1D *fhSubJettiness1_FJ_OP_CA;
229  TH1D *fhSubJettiness1_FJ_OP_WTA_KT;
230  TH1D *fhSubJettiness1_FJ_OP_WTA_CA;
231  TH1D *fhSubJettiness1_FJ_MIN;
232  TH1D *fhSubJettiness2_FJ_AKT;
233  TH1D *fhSubJettiness2_FJ_KT;
234  TH1D *fhSubJettiness2_FJ_CA;
235  TH1D *fhSubJettiness2_FJ_WTA_KT;
236  TH1D *fhSubJettiness2_FJ_WTA_CA;
237  TH1D *fhSubJettiness2_FJ_OP_AKT;
238  TH1D *fhSubJettiness2_FJ_OP_KT;
239  TH1D *fhSubJettiness2_FJ_OP_CA;
240  TH1D *fhSubJettiness2_FJ_OP_WTA_KT;
241  TH1D *fhSubJettiness2_FJ_OP_WTA_CA;
242  TH1D *fhSubJettiness2_FJ_MIN;
243  TH1D *fhSubJettiness2to1_FJ_AKT;
244  TH1D *fhSubJettiness2to1_FJ_KT;
245  TH1D *fhSubJettiness2to1_FJ_CA;
246  TH1D *fhSubJettiness2to1_FJ_WTA_KT;
247  TH1D *fhSubJettiness2to1_FJ_WTA_CA;
248  TH1D *fhSubJettiness2to1_FJ_OP_AKT;
249  TH1D *fhSubJettiness2to1_FJ_OP_KT;
250  TH1D *fhSubJettiness2to1_FJ_OP_CA;
251  TH1D *fhSubJettiness2to1_FJ_OP_WTA_KT;
252  TH1D *fhSubJettiness2to1_FJ_OP_WTA_CA;
253  TH1D *fhSubJettiness2to1_FJ_MIN;
254  TTree *fTreeResponseMatrixAxis; //Tree with tagging variables subtracted MC or true MC or raw
255 
256  private:
257  AliAnalysisTaskSubJetFraction(const AliAnalysisTaskSubJetFraction&); // not implemented
258  AliAnalysisTaskSubJetFraction &operator=(const AliAnalysisTaskSubJetFraction&); // not implemented
259 
260  ClassDef(AliAnalysisTaskSubJetFraction, 6)
261 };
262 #endif
263 
AliAnalysisTaskSubJetFraction::SetNsubUnNormMeasure
void SetNsubUnNormMeasure(Bool_t NsubMeasure)
Definition: AliAnalysisTaskSubJetFraction.h:85
AliAnalysisTaskSubJetFraction::fCentSelectOn
Bool_t fCentSelectOn
Random number generator.
Definition: AliAnalysisTaskSubJetFraction.h:123
Double_t
double Double_t
Definition: External.C:58
TH3F
Definition: External.C:260
AliAnalysisTaskEmcalJet.h
TH2F
Definition: External.C:236
AliAnalysisTaskSubJetFraction::Run
Bool_t Run()
Run function. This is the core function of the analysis and contains the user code. Therefore users have to implement this function.
Definition: AliAnalysisTaskSubJetFraction.cxx:736
TH3
Definition: External.C:244
AliAnalysisTaskSubJetFraction::operator=
AliAnalysisTaskSubJetFraction & operator=(const AliAnalysisTaskSubJetFraction &)
AliAnalysisTaskSubJetFraction::SetPtTriggerSelections
void SetPtTriggerSelections(Float_t PtMin, Float_t PtMax)
Definition: AliAnalysisTaskSubJetFraction.h:62
AliAnalysisTaskSubJetFraction::ModifyJet
std::pair< fastjet::PseudoJet, fastjet::ClusterSequence * > ModifyJet(AliEmcalJet *Jet, Int_t JetContNb, TString Modification)
Definition: AliAnalysisTaskSubJetFraction.cxx:1800
AliAnalysisTaskSubJetFraction::PTD
Double_t PTD(AliEmcalJet *Jet, Int_t JetContNb)
Definition: AliAnalysisTaskSubJetFraction.cxx:1451
AliAnalysisTaskSubJetFraction::RetrieveEventObjects
Bool_t RetrieveEventObjects()
Retrieve common objects from event.
Definition: AliAnalysisTaskSubJetFraction.cxx:1743
AliAnalysisTaskSubJetFraction::AddExtraProng
std::vector< fastjet::PseudoJet > AddExtraProng(std::vector< fastjet::PseudoJet > fInputVectors, Double_t Distance, Double_t PtFrac)
Definition: AliAnalysisTaskSubJetFraction.cxx:1864
c
TCanvas * c
Definition: TestFitELoss.C:172
AliAnalysisTaskSubJetFraction::NSubJettiness
Double_t NSubJettiness(AliEmcalJet *Jet, Int_t JetContNb, Double_t JetRadius, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t A, Int_t B)
Definition: AliAnalysisTaskSubJetFraction.cxx:1584
AliAnalysisTaskSubJetFraction::SetReclusteringAlgorithm
void SetReclusteringAlgorithm(Int_t ReclusteringAlgorithm)
Definition: AliAnalysisTaskSubJetFraction.h:81
AliAnalysisTaskSubJetFraction::RelativePhi
Double_t RelativePhi(Double_t Phi1, Double_t Phi2)
Definition: AliAnalysisTaskSubJetFraction.cxx:1345
AliFJWrapper.h
AliAnalysisTaskSubJetFraction::FillHistograms
Bool_t FillHistograms()
Function filling histograms.
Definition: AliAnalysisTaskSubJetFraction.cxx:745
Int_t
int Int_t
Definition: External.C:63
Float_t
float Float_t
Definition: External.C:68
AliAnalysisTaskSubJetFraction::SelectTriggerHadron
Int_t SelectTriggerHadron(Float_t PtMin, Float_t PtMax)
Definition: AliAnalysisTaskSubJetFraction.cxx:1359
TH2D
Definition: External.C:228
TH1D
Definition: External.C:212
AliAnalysisTaskSubJetFraction::SetSubJetAlgorithm
void SetSubJetAlgorithm(Int_t SubJetAlgorithm)
Definition: AliAnalysisTaskSubJetFraction.h:72
Option
Definition: Option.C:68
AliAnalysisTaskSubJetFraction::SubJetFraction
Double_t SubJetFraction(AliEmcalJet *Jet, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t Type, Bool_t Add, Bool_t Loss)
Definition: AliAnalysisTaskSubJetFraction.cxx:1544
AliAnalysisTaskSubJetFraction::SetSharedFractionPtMin
void SetSharedFractionPtMin(Double_t SharedFractionPtMin)
Definition: AliAnalysisTaskSubJetFraction.h:76
AliAnalysisTaskSubJetFraction::SetSubJetRadius
void SetSubJetRadius(Float_t SubJetRadius)
Definition: AliAnalysisTaskSubJetFraction.h:73
AliAnalysisTaskSubJetFraction::Recluster
AliEmcalJetFinder * Recluster(AliEmcalJet *Jet, Int_t JetContNb, Double_t SubJetRadius, Double_t SubJetMinPt, Int_t Algorithm, const char *Name)
Definition: AliAnalysisTaskSubJetFraction.cxx:1472
AliAnalysisTaskSubJetFraction::FjNSubJettinessFastJet
Double_t FjNSubJettinessFastJet(std::pair< fastjet::PseudoJet, fastjet::ClusterSequence * > Jet_ClusterSequence, Int_t JetContNb, Int_t N, Int_t Algorithm, Double_t Beta, Int_t Option=0, Double_t Beta_SD=0, Double_t ZCut=0.1)
Definition: AliAnalysisTaskSubJetFraction.cxx:1889
AliAnalysisTaskSubJetFraction::FjNSubJettiness
Double_t FjNSubJettiness(AliEmcalJet *Jet, Int_t JetContNb, Int_t N, Int_t Algorithm, Double_t Beta, Int_t Option=0, Double_t Beta_SD=0, Double_t ZCut=0.1)
Definition: AliAnalysisTaskSubJetFraction.cxx:1626
AliAnalysisTaskSubJetFraction::RandomiseTracks
std::vector< fastjet::PseudoJet > RandomiseTracks(AliEmcalJet *Jet, std::vector< fastjet::PseudoJet > fInputVectors)
Definition: AliAnalysisTaskSubJetFraction.cxx:1840
TH2
Definition: External.C:220
nVar
const Int_t nVar
Definition: AliAnalysisTaskSubJetFraction.h:20
AliAnalysisTaskEmcalJet
Base task in the EMCAL jet framework.
Definition: AliAnalysisTaskEmcalJet.h:30
AliEmcalJet
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
AliAnalysisTaskSubJetFraction::Angularity
Double_t Angularity(AliEmcalJet *Jet, Int_t JetContNb)
Definition: AliAnalysisTaskSubJetFraction.cxx:1427
Option_t
const char Option_t
Definition: External.C:48
AliAnalysisTaskSubJetFraction::RelativePhiEventPlane
Double_t RelativePhiEventPlane(Double_t EventPlane, Double_t Phi)
Definition: AliAnalysisTaskSubJetFraction.cxx:1335
Bool_t
bool Bool_t
Definition: External.C:53
AliAnalysisTaskSubJetFraction::SubJetOrdering
Double_t SubJetOrdering(AliEmcalJet *Jet, AliEmcalJetFinder *Reclusterer, Int_t N, Int_t Type, Bool_t Index)
Definition: AliAnalysisTaskSubJetFraction.cxx:1497
AliJetContainer
Container for jet within the EMCAL jet framework.
Definition: AliJetContainer.h:37
TH1
Definition: External.C:196
AliClusterContainer.h