AliPhysics  9b6b435 (9b6b435)
AliAnalysisTaskEmcalQGTagging.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKEMCALQGTAGGING_H
2 #define ALIANALYSISTASKEMCALQGTAGGING_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;
17 #include "AliFJWrapper.h"
18 
19 
21  public:
22 
23  enum JetShapeType {
24  kMCTrue = 0, // generated jets only
25  kTrueDet =1, // detector and generated jets
26  kData = 2, // raw data
27  kDetEmb = 3, //detector embedded jets
32  };
33  enum JetShapeSub {
34  kNoSub = 0,
35  kConstSub = 1,
36  kDerivSub = 2
37  };
40  kRecoil = 1
41  };
42 
46  };
47 
49  AliAnalysisTaskEmcalQGTagging(const char *name);
51 
53  void Terminate(Option_t *option);
54 
55  //Setters
62  void SetRMatching(Float_t f) { fRMatching = f ;}
64  void SetPtTriggerSelections(Float_t minpT, Float_t maxpT) { fminpTTrig = minpT; fmaxpTTrig = maxpT; }
76  void SetMinCentrality(Float_t t) { fCentMin = t ; }
77  void SetMaxCentrality(Float_t t) { fCentMax = t ; }
79  void SetHolePos(Float_t poshole) { fHolePos = poshole; }
80  void SetHoleWidth(Float_t holewidth) { fHoleWidth = holewidth; }
82  protected:
84  Bool_t Run();
86 
87  Float_t GetJetMass(AliEmcalJet *jet,Int_t jetContNb);
88  Float_t Angularity(AliEmcalJet *jet, Int_t jetContNb);
89  Float_t GetJetAngularity(AliEmcalJet *jet, Int_t jetContNb);
90  Float_t Coronna(AliEmcalJet *jet, Int_t jetContNb);
91  Float_t GetJetCoronna(AliEmcalJet *jet, Int_t jetContNb);
92  Float_t PTD(AliEmcalJet *jet, Int_t jetContNb);
93  Float_t GetJetpTD(AliEmcalJet *jet, Int_t jetContNb);
94  Float_t Circularity(AliEmcalJet *jet, Int_t jetContNb);
95  Float_t GetJetCircularity(AliEmcalJet *jet, Int_t jetContNb);
96  Float_t LeSub(AliEmcalJet *jet, Int_t jetContNb);
97  Float_t GetJetLeSub(AliEmcalJet *jet, Int_t jetContNb);
99  Float_t GetSigma2(AliEmcalJet *jet, Int_t jetContNb);
100  Float_t Sigma2(AliEmcalJet *jet, Int_t jetContNb);
101 
102  Int_t SelectTrigger(Float_t minpT, Float_t maxpT);
104  void RecursiveParents(AliEmcalJet *fJet,AliJetContainer *fJetCont);
105  void RecursiveParentsMCAverage(AliEmcalJet *fJet,Int_t km, Double_t &aver1, Double_t &aver2);
106  void CheckSubjetResolution(AliEmcalJet *fJet,AliJetContainer *fJetCont, AliEmcalJet *fJetM,AliJetContainer *fJetContM);
107  Bool_t CheckClosePartner(Int_t index, AliEmcalJet *fJet,AliVParticle *fTrack, AliParticleContainer *fTrackCont);
108  Int_t fContainer; // jets to be analyzed 0 for Base, 1 for subtracted.
109  Float_t fMinFractionShared; // only fill histos for jets if shared fraction larger than X
110  JetShapeType fJetShapeType; // jet type to be used
111  JetShapeSub fJetShapeSub; // jet subtraction to be used
112  JetSelectionType fJetSelection; // Jet selection: inclusive/recoil jet
113  Float_t fShapesVar[8]; // jet shapes used for the tagging
116  Int_t fSelectedShapes; //chose set of shapes
117  Float_t fminpTTrig; //min - max pT for trigger particle in case of recoil jet
119  Float_t fangWindowRecoil; //angular window for btb recoil analysis
120  Int_t fSemigoodCorrect; //if==1 we run over semigood runs
121  Float_t fHolePos; //position in radians of the bad TPC sector
122  Float_t fHoleWidth; //width of the hole in radians
123  Bool_t fCentSelectOn; // switch on/off centrality selection
124  Float_t fCentMin; // min centrality value
125  Float_t fCentMax; // max centrality value
126  Bool_t fOneConstSelectOn; // switch on/off one constituent selection
127  Bool_t fTrackCheckPlots; //switch on qa plots
128  Bool_t fCheckResolution; //check subjet energy resolution
129  Float_t fSubjetCutoff; //angular cutoff for subjets at det/gen level
130  Float_t fMinPtConst; //constituent pt cutoff
131  Float_t fHardCutoff; //hard cutoff in the iterative declustering
132  Bool_t fDoTwoTrack; //switch to consider 2 track effects
133  Float_t fPhiCutValue; //cuts from HBT
134  Float_t fEtaCutValue; //cuts from HBT
135  Float_t fMagFieldPolarity; //polarity, to calculate phimin
137 
138 
146  TH1F *fPtJet;
147  TH2F *fhpTjetpT; //control p[lot fo the recoil analysis
148  TH1F *fhPt;
149  TH1F *fhPhi;
150  TH1F *fhTrackPhi; //control plot for tracks
151  THnSparse *fHLundIterative;// iterative declustering
152  THnSparse *fHCheckResolutionSubjets;// to evaluate energy resolution of subjets as function fo apperture angle
154 
155  TTree *fTreeObservableTagging; // Tree with tagging variables subtracted MC or true MC or raw
156 
157  private:
160 
161  ClassDef(AliAnalysisTaskEmcalQGTagging, 6)
162 };
163 #endif
164 
Int_t SelectTrigger(Float_t minpT, Float_t maxpT)
Float_t Coronna(AliEmcalJet *jet, Int_t jetContNb)
double Double_t
Definition: External.C:58
Definition: External.C:260
Definition: External.C:236
Bool_t RetrieveEventObjects()
Retrieve common objects from event.
Bool_t FillHistograms()
Function filling histograms.
Definition: External.C:244
Float_t GetJetMass(AliEmcalJet *jet, Int_t jetContNb)
void SetPtTriggerSelections(Float_t minpT, Float_t maxpT)
AliAnalysisTaskEmcalQGTagging & operator=(const AliAnalysisTaskEmcalQGTagging &)
TCanvas * c
Definition: TestFitELoss.C:172
Float_t GetJetAngularity(AliEmcalJet *jet, Int_t jetContNb)
Container for particles within the EMCAL framework.
Float_t GetJetCoronna(AliEmcalJet *jet, Int_t jetContNb)
int Int_t
Definition: External.C:63
Float_t GetJetCircularity(AliEmcalJet *jet, Int_t jetContNb)
Float_t PTD(AliEmcalJet *jet, Int_t jetContNb)
float Float_t
Definition: External.C:68
Float_t GetSigma2(AliEmcalJet *jet, Int_t jetContNb)
Float_t Angularity(AliEmcalJet *jet, Int_t jetContNb)
Float_t LeSub(AliEmcalJet *jet, Int_t jetContNb)
void RecursiveParentsMCAverage(AliEmcalJet *fJet, Int_t km, Double_t &aver1, Double_t &aver2)
Definition: External.C:220
Float_t Sigma2(AliEmcalJet *jet, Int_t jetContNb)
Double_t RelativePhi(Double_t mphi, Double_t vphi)
Float_t GetJetNumberOfConstituents(AliEmcalJet *jet, Int_t jetContNb)
void RecursiveParents(AliEmcalJet *fJet, AliJetContainer *fJetCont)
Bool_t CheckClosePartner(Int_t index, AliEmcalJet *fJet, AliVParticle *fTrack, AliParticleContainer *fTrackCont)
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
Float_t GetJetpTD(AliEmcalJet *jet, Int_t jetContNb)
const char Option_t
Definition: External.C:48
Float_t Circularity(AliEmcalJet *jet, Int_t jetContNb)
bool Bool_t
Definition: External.C:53
Float_t GetJetLeSub(AliEmcalJet *jet, Int_t jetContNb)
Container for jet within the EMCAL jet framework.
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: External.C:196
void CheckSubjetResolution(AliEmcalJet *fJet, AliJetContainer *fJetCont, AliEmcalJet *fJetM, AliJetContainer *fJetContM)