AliPhysics  d497547 (d497547)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliEmcalJetUtilityGenSubtractor.cxx
Go to the documentation of this file.
2 #include "AliEmcalJet.h"
3 #include "AliRhoParameter.h"
4 #include "AliEmcalJetTask.h"
5 
7 
8 //______________________________________________________________________________
11  fDoGenericSubtractionJetMass(kFALSE),
12  fDoGenericSubtractionGR(kFALSE),
13  fDoGenericSubtractionExtraJetShapes(kFALSE),
14  fDoGenericSubtractionNsubjettiness(kFALSE),
15  fUseExternalBkg(kFALSE),
16  fRhoName(""),
17  fRhomName(""),
18  fRho(0),
19  fRhom(0),
20  fRMax(0.4),
21  fDRStep(0.04),
22  fPtMinGR(40.),
23  fRhoParam(0),
24  fRhomParam(0)
25 {
26  // Dummy constructor.
27 
28 }
29 
30 //______________________________________________________________________________
32  AliEmcalJetUtility(name),
33  fDoGenericSubtractionJetMass(kFALSE),
34  fDoGenericSubtractionGR(kFALSE),
35  fDoGenericSubtractionExtraJetShapes(kFALSE),
36  fDoGenericSubtractionNsubjettiness(kFALSE),
37  fUseExternalBkg(kFALSE),
38  fRhoName(""),
39  fRhomName(""),
40  fRho(0),
41  fRhom(0),
42  fRMax(0.4),
43  fDRStep(0.04),
44  fPtMinGR(40.),
45  fRhoParam(0),
46  fRhomParam(0)
47 {
48  // Default constructor.
49 }
50 
51 //______________________________________________________________________________
53  AliEmcalJetUtility(other),
54  fDoGenericSubtractionJetMass(other.fDoGenericSubtractionJetMass),
55  fDoGenericSubtractionGR(other.fDoGenericSubtractionGR),
56  fDoGenericSubtractionExtraJetShapes(other.fDoGenericSubtractionExtraJetShapes),
57  fDoGenericSubtractionNsubjettiness(other.fDoGenericSubtractionNsubjettiness),
58  fUseExternalBkg(other.fUseExternalBkg),
59  fRhoName(other.fRhoName),
60  fRhomName(other.fRhomName),
61  fRho(other.fRho),
62  fRhom(other.fRhom),
63  fRMax(other.fRMax),
64  fDRStep(other.fDRStep),
65  fPtMinGR(other.fPtMinGR),
66  fRhoParam(other.fRhoParam),
67  fRhomParam(other.fRhomParam)
68 {
69  // Copy constructor.
70 }
71 
72 //______________________________________________________________________________
74 {
75  // Assignment.
76 
77  if (&other == this) return *this;
84  fRhoName = other.fRhoName;
85  fRhomName = other.fRhomName;
86  fRho = other.fRho;
87  fRhom = other.fRhom;
88  fRMax = other.fRMax;
89  fDRStep = other.fDRStep;
90  fPtMinGR = other.fPtMinGR;
91  fRhoParam = other.fRhoParam;
92  fRhomParam = other.fRhomParam;
93  return *this;
94 }
95 
96 //______________________________________________________________________________
98 {
99  // Initialize the utility.
100 
101  if (!fRhoName.IsNull() && !fRhoParam) { // get rho from the event
102  if(!fJetTask) return;
103  fRhoParam = dynamic_cast<AliRhoParameter*>(fJetTask->InputEvent()->FindListObject(fRhoName));
104  if (!fRhoParam) {
105  AliError(Form("%s: Could not retrieve rho %s!", GetName(), fRhoName.Data()));
106  return;
107  }
108  }
109 
110  if (!fRhomName.IsNull() && !fRhomParam) { // get rhom from the event
111  fRhomParam = dynamic_cast<AliRhoParameter*>(fJetTask->InputEvent()->FindListObject(fRhomName));
112  if (!fRhomParam) {
113  AliError(Form("%s: Could not retrieve rho_m %s!", GetName(), fRhomName.Data()));
114  return;
115  }
116  }
117 
118  fInit = kTRUE;
119 }
120 
121 
122 //______________________________________________________________________________
124 {
125  // Prepare the utility.
126 
127 }
128 
129 //______________________________________________________________________________
131 {
132  // Prepare the utility.
133 
134  if (!fInit) return;
135 
136  if (fRhoParam) fRho = fRhoParam->GetVal();
137  if (fRhomParam) fRhom = fRhomParam->GetVal();
138 
139  //run generic subtractor
143  }
144 
153  }
154 
170  }
171 }
172 
173 //______________________________________________________________________________
175 {
176  // Proceess each jet.
177 
178  if (!fInit) return;
179 
180 #ifdef FASTJET_VERSION
181 
183  std::vector<fastjet::contrib::GenericSubtractorInfo> jetMassInfo = fjw.GetGenSubtractorInfoJetMass();
184  Int_t n = (Int_t)jetMassInfo.size();
185  if(n > ij && n > 0) {
186  jet->GetShapeProperties()->SetFirstDerivative(jetMassInfo[ij].first_derivative());
187  jet->GetShapeProperties()->SetSecondDerivative(jetMassInfo[ij].second_derivative());
188  jet->GetShapeProperties()->SetFirstOrderSubtracted(jetMassInfo[ij].first_order_subtracted());
189  jet->GetShapeProperties()->SetSecondOrderSubtracted(jetMassInfo[ij].second_order_subtracted());
190  }
191  }
192 
193  //here do generic subtraction for angular structure function
194  Double_t ptcorr = jet->Pt()-fjw.GetJetArea(ij)*fRho;
195  if (fDoGenericSubtractionGR && ptcorr>fPtMinGR) {
197  fRMax = fJetTask->GetRadius()+0.2;
199  fjw.DoGenericSubtractionGR(ij);
200  std::vector<double> num = fjw.GetGRNumerator();
201  std::vector<double> den = fjw.GetGRDenominator();
202  std::vector<double> nums = fjw.GetGRNumeratorSub();
203  std::vector<double> dens = fjw.GetGRDenominatorSub();
204  //pass this to AliEmcalJet
205  jet->GetShapeProperties()->SetGRNumSize(num.size());
206  jet->GetShapeProperties()->SetGRDenSize(den.size());
207  jet->GetShapeProperties()->SetGRNumSubSize(nums.size());
208  jet->GetShapeProperties()->SetGRDenSubSize(dens.size());
209  Int_t nsize = (Int_t)num.size();
210  for (Int_t g = 0; g < nsize; ++g) {
211  jet->GetShapeProperties()->AddGRNumAt(num[g],g);
212  jet->GetShapeProperties()->AddGRNumSubAt(nums[g],g);
213  }
214  Int_t dsize = (Int_t)den.size();
215  for (Int_t g = 0; g < dsize; ++g) {
216  jet->GetShapeProperties()->AddGRDenAt(den[g], g);
217  jet->GetShapeProperties()->AddGRDenSubAt(dens[g], g);
218  }
219  }
220 
222  std::vector<fastjet::contrib::GenericSubtractorInfo> jetAngularityInfo = fjw.GetGenSubtractorInfoJetAngularity();
223  Int_t na = (Int_t)jetAngularityInfo.size();
224  if(na > ij && na > 0) {
225  jet->GetShapeProperties()->SetFirstDerivativeAngularity(jetAngularityInfo[ij].first_derivative());
226  jet->GetShapeProperties()->SetSecondDerivativeAngularity(jetAngularityInfo[ij].second_derivative());
227  jet->GetShapeProperties()->SetFirstOrderSubtractedAngularity(jetAngularityInfo[ij].first_order_subtracted());
228  jet->GetShapeProperties()->SetSecondOrderSubtractedAngularity(jetAngularityInfo[ij].second_order_subtracted());
229  }
230 
231  std::vector<fastjet::contrib::GenericSubtractorInfo> jetpTDInfo = fjw.GetGenSubtractorInfoJetpTD();
232  Int_t np = (Int_t)jetpTDInfo.size();
233  if(np > ij && np > 0) {
234  jet->GetShapeProperties()->SetFirstDerivativepTD(jetpTDInfo[ij].first_derivative());
235  jet->GetShapeProperties()->SetSecondDerivativepTD(jetpTDInfo[ij].second_derivative());
236  jet->GetShapeProperties()->SetFirstOrderSubtractedpTD(jetpTDInfo[ij].first_order_subtracted());
237  jet->GetShapeProperties()->SetSecondOrderSubtractedpTD(jetpTDInfo[ij].second_order_subtracted());
238  }
239 
240  std::vector<fastjet::contrib::GenericSubtractorInfo> jetCircularityInfo = fjw.GetGenSubtractorInfoJetCircularity();
241  Int_t nc = (Int_t)jetCircularityInfo.size();
242  if(nc > ij && nc > 0) {
243  jet->GetShapeProperties()->SetFirstDerivativeCircularity(jetCircularityInfo[ij].first_derivative());
244  jet->GetShapeProperties()->SetSecondDerivativeCircularity(jetCircularityInfo[ij].second_derivative());
245  jet->GetShapeProperties()->SetFirstOrderSubtractedCircularity(jetCircularityInfo[ij].first_order_subtracted());
246  jet->GetShapeProperties()->SetSecondOrderSubtractedCircularity(jetCircularityInfo[ij].second_order_subtracted());
247  }
248 
249  std::vector<fastjet::contrib::GenericSubtractorInfo> jetSigma2Info = fjw.GetGenSubtractorInfoJetSigma2();
250  Int_t ns = (Int_t)jetSigma2Info.size();
251  if (ns > ij && ns > 0) {
252  jet->GetShapeProperties()->SetFirstDerivativeSigma2(jetSigma2Info[ij].first_derivative());
253  jet->GetShapeProperties()->SetSecondDerivativeSigma2(jetSigma2Info[ij].second_derivative());
254  jet->GetShapeProperties()->SetFirstOrderSubtractedSigma2(jetSigma2Info[ij].first_order_subtracted());
255  jet->GetShapeProperties()->SetSecondOrderSubtractedSigma2(jetSigma2Info[ij].second_order_subtracted());
256  }
257 
258 
259  std::vector<fastjet::contrib::GenericSubtractorInfo> jetConstituentInfo = fjw.GetGenSubtractorInfoJetConstituent();
260  Int_t nco = (Int_t)jetConstituentInfo.size();
261  if(nco > ij && nco > 0) {
262  jet->GetShapeProperties()->SetFirstDerivativeConstituent(jetConstituentInfo[ij].first_derivative());
263  jet->GetShapeProperties()->SetSecondDerivativeConstituent(jetConstituentInfo[ij].second_derivative());
264  jet->GetShapeProperties()->SetFirstOrderSubtractedConstituent(jetConstituentInfo[ij].first_order_subtracted());
265  jet->GetShapeProperties()->SetSecondOrderSubtractedConstituent(jetConstituentInfo[ij].second_order_subtracted());
266  }
267 
268  std::vector<fastjet::contrib::GenericSubtractorInfo> jetLeSubInfo = fjw.GetGenSubtractorInfoJetLeSub();
269  Int_t nlsub = (Int_t)jetLeSubInfo.size();
270  if(nlsub > ij && nlsub > 0) {
271  jet->GetShapeProperties()->SetFirstDerivativeLeSub(jetLeSubInfo[ij].first_derivative());
272  jet->GetShapeProperties()->SetSecondDerivativeLeSub(jetLeSubInfo[ij].second_derivative());
273  jet->GetShapeProperties()->SetFirstOrderSubtractedLeSub(jetLeSubInfo[ij].first_order_subtracted());
274  jet->GetShapeProperties()->SetSecondOrderSubtractedLeSub(jetLeSubInfo[ij].second_order_subtracted());
275  }
276  }
277 
279  std::vector<fastjet::contrib::GenericSubtractorInfo> jet1subjettinessktInfo = fjw.GetGenSubtractorInfoJet1subjettiness_kt();
280  Int_t n1subjettiness_kt = (Int_t)jet1subjettinessktInfo.size();
281  if(n1subjettiness_kt > ij && n1subjettiness_kt > 0) {
282  jet->GetShapeProperties()->SetFirstDerivative1subjettiness_kt(jet1subjettinessktInfo[ij].first_derivative());
283  jet->GetShapeProperties()->SetSecondDerivative1subjettiness_kt(jet1subjettinessktInfo[ij].second_derivative());
284  jet->GetShapeProperties()->SetFirstOrderSubtracted1subjettiness_kt(jet1subjettinessktInfo[ij].first_order_subtracted());
285  jet->GetShapeProperties()->SetSecondOrderSubtracted1subjettiness_kt(jet1subjettinessktInfo[ij].second_order_subtracted());
286  }
287 
288  std::vector<fastjet::contrib::GenericSubtractorInfo> jet2subjettinessktInfo = fjw.GetGenSubtractorInfoJet2subjettiness_kt();
289  Int_t n2subjettiness_kt = (Int_t)jet2subjettinessktInfo.size();
290  if(n2subjettiness_kt > ij && n2subjettiness_kt > 0) {
291  jet->GetShapeProperties()->SetFirstDerivative2subjettiness_kt(jet2subjettinessktInfo[ij].first_derivative());
292  jet->GetShapeProperties()->SetSecondDerivative2subjettiness_kt(jet2subjettinessktInfo[ij].second_derivative());
293  jet->GetShapeProperties()->SetFirstOrderSubtracted2subjettiness_kt(jet2subjettinessktInfo[ij].first_order_subtracted());
294  jet->GetShapeProperties()->SetSecondOrderSubtracted2subjettiness_kt(jet2subjettinessktInfo[ij].second_order_subtracted());
295  }
296 
297  std::vector<fastjet::contrib::GenericSubtractorInfo> jet3subjettinessktInfo = fjw.GetGenSubtractorInfoJet3subjettiness_kt();
298  Int_t n3subjettiness_kt = (Int_t)jet3subjettinessktInfo.size();
299  if(n3subjettiness_kt > ij && n3subjettiness_kt > 0) {
300  jet->GetShapeProperties()->SetFirstDerivative3subjettiness_kt(jet3subjettinessktInfo[ij].first_derivative());
301  jet->GetShapeProperties()->SetSecondDerivative3subjettiness_kt(jet3subjettinessktInfo[ij].second_derivative());
302  jet->GetShapeProperties()->SetFirstOrderSubtracted3subjettiness_kt(jet3subjettinessktInfo[ij].first_order_subtracted());
303  jet->GetShapeProperties()->SetSecondOrderSubtracted3subjettiness_kt(jet3subjettinessktInfo[ij].second_order_subtracted());
304  }
305 
306  std::vector<fastjet::contrib::GenericSubtractorInfo> jetOpeningAnglektInfo = fjw.GetGenSubtractorInfoJetOpeningAngle_kt();
307  Int_t nOpeningAngle_kt = (Int_t)jetOpeningAnglektInfo.size();
308  if(nOpeningAngle_kt > ij && nOpeningAngle_kt > 0) {
309  jet->GetShapeProperties()->SetFirstDerivativeOpeningAngle_kt(jetOpeningAnglektInfo[ij].first_derivative());
310  jet->GetShapeProperties()->SetSecondDerivativeOpeningAngle_kt(jetOpeningAnglektInfo[ij].second_derivative());
311  jet->GetShapeProperties()->SetFirstOrderSubtractedOpeningAngle_kt(jetOpeningAnglektInfo[ij].first_order_subtracted());
312  jet->GetShapeProperties()->SetSecondOrderSubtractedOpeningAngle_kt(jetOpeningAnglektInfo[ij].second_order_subtracted());
313  }
314  std::vector<fastjet::contrib::GenericSubtractorInfo> jet1subjettinesscaInfo = fjw.GetGenSubtractorInfoJet1subjettiness_ca();
315  Int_t n1subjettiness_ca = (Int_t)jet1subjettinesscaInfo.size();
316  if(n1subjettiness_ca > ij && n1subjettiness_ca > 0) {
317  jet->GetShapeProperties()->SetFirstDerivative1subjettiness_ca(jet1subjettinesscaInfo[ij].first_derivative());
318  jet->GetShapeProperties()->SetSecondDerivative1subjettiness_ca(jet1subjettinesscaInfo[ij].second_derivative());
319  jet->GetShapeProperties()->SetFirstOrderSubtracted1subjettiness_ca(jet1subjettinesscaInfo[ij].first_order_subtracted());
320  jet->GetShapeProperties()->SetSecondOrderSubtracted1subjettiness_ca(jet1subjettinesscaInfo[ij].second_order_subtracted());
321  }
322 
323  std::vector<fastjet::contrib::GenericSubtractorInfo> jet2subjettinesscaInfo = fjw.GetGenSubtractorInfoJet2subjettiness_ca();
324  Int_t n2subjettiness_ca = (Int_t)jet2subjettinesscaInfo.size();
325  if(n2subjettiness_ca > ij && n2subjettiness_ca > 0) {
326  jet->GetShapeProperties()->SetFirstDerivative2subjettiness_ca(jet2subjettinesscaInfo[ij].first_derivative());
327  jet->GetShapeProperties()->SetSecondDerivative2subjettiness_ca(jet2subjettinesscaInfo[ij].second_derivative());
328  jet->GetShapeProperties()->SetFirstOrderSubtracted2subjettiness_ca(jet2subjettinesscaInfo[ij].first_order_subtracted());
329  jet->GetShapeProperties()->SetSecondOrderSubtracted2subjettiness_ca(jet2subjettinesscaInfo[ij].second_order_subtracted());
330  }
331 
332  std::vector<fastjet::contrib::GenericSubtractorInfo> jetOpeningAnglecaInfo = fjw.GetGenSubtractorInfoJetOpeningAngle_ca();
333  Int_t nOpeningAngle_ca = (Int_t)jetOpeningAnglecaInfo.size();
334  if(nOpeningAngle_ca > ij && nOpeningAngle_ca > 0) {
335  jet->GetShapeProperties()->SetFirstDerivativeOpeningAngle_ca(jetOpeningAnglecaInfo[ij].first_derivative());
336  jet->GetShapeProperties()->SetSecondDerivativeOpeningAngle_ca(jetOpeningAnglecaInfo[ij].second_derivative());
337  jet->GetShapeProperties()->SetFirstOrderSubtractedOpeningAngle_ca(jetOpeningAnglecaInfo[ij].first_order_subtracted());
338  jet->GetShapeProperties()->SetSecondOrderSubtractedOpeningAngle_ca(jetOpeningAnglecaInfo[ij].second_order_subtracted());
339  }
340  std::vector<fastjet::contrib::GenericSubtractorInfo> jet1subjettinessakt02Info = fjw.GetGenSubtractorInfoJet1subjettiness_akt02();
341  Int_t n1subjettiness_akt02 = (Int_t)jet1subjettinessakt02Info.size();
342  if(n1subjettiness_akt02 > ij && n1subjettiness_akt02 > 0) {
343  jet->GetShapeProperties()->SetFirstDerivative1subjettiness_akt02(jet1subjettinessakt02Info[ij].first_derivative());
344  jet->GetShapeProperties()->SetSecondDerivative1subjettiness_akt02(jet1subjettinessakt02Info[ij].second_derivative());
345  jet->GetShapeProperties()->SetFirstOrderSubtracted1subjettiness_akt02(jet1subjettinessakt02Info[ij].first_order_subtracted());
346  jet->GetShapeProperties()->SetSecondOrderSubtracted1subjettiness_akt02(jet1subjettinessakt02Info[ij].second_order_subtracted());
347  }
348 
349  std::vector<fastjet::contrib::GenericSubtractorInfo> jet2subjettinessakt02Info = fjw.GetGenSubtractorInfoJet2subjettiness_akt02();
350  Int_t n2subjettiness_akt02 = (Int_t)jet2subjettinessakt02Info.size();
351  if(n2subjettiness_akt02 > ij && n2subjettiness_akt02 > 0) {
352  jet->GetShapeProperties()->SetFirstDerivative2subjettiness_akt02(jet2subjettinessakt02Info[ij].first_derivative());
353  jet->GetShapeProperties()->SetSecondDerivative2subjettiness_akt02(jet2subjettinessakt02Info[ij].second_derivative());
354  jet->GetShapeProperties()->SetFirstOrderSubtracted2subjettiness_akt02(jet2subjettinessakt02Info[ij].first_order_subtracted());
355  jet->GetShapeProperties()->SetSecondOrderSubtracted2subjettiness_akt02(jet2subjettinessakt02Info[ij].second_order_subtracted());
356  }
357 
358  std::vector<fastjet::contrib::GenericSubtractorInfo> jetOpeningAngleakt02Info = fjw.GetGenSubtractorInfoJetOpeningAngle_akt02();
359  Int_t nOpeningAngle_akt02 = (Int_t)jetOpeningAngleakt02Info.size();
360  if(nOpeningAngle_akt02 > ij && nOpeningAngle_akt02 > 0) {
361  jet->GetShapeProperties()->SetFirstDerivativeOpeningAngle_akt02(jetOpeningAngleakt02Info[ij].first_derivative());
362  jet->GetShapeProperties()->SetSecondDerivativeOpeningAngle_akt02(jetOpeningAngleakt02Info[ij].second_derivative());
363  jet->GetShapeProperties()->SetFirstOrderSubtractedOpeningAngle_akt02(jetOpeningAngleakt02Info[ij].first_order_subtracted());
364  jet->GetShapeProperties()->SetSecondOrderSubtractedOpeningAngle_akt02(jetOpeningAngleakt02Info[ij].second_order_subtracted());
365  }
366  std::vector<fastjet::contrib::GenericSubtractorInfo> jet1subjettinesscasdInfo = fjw.GetGenSubtractorInfoJet1subjettiness_casd();
367  Int_t n1subjettiness_casd = (Int_t)jet1subjettinesscasdInfo.size();
368  if(n1subjettiness_casd > ij && n1subjettiness_casd > 0) {
369  jet->GetShapeProperties()->SetFirstDerivative1subjettiness_casd(jet1subjettinesscasdInfo[ij].first_derivative());
370  jet->GetShapeProperties()->SetSecondDerivative1subjettiness_casd(jet1subjettinesscasdInfo[ij].second_derivative());
371  jet->GetShapeProperties()->SetFirstOrderSubtracted1subjettiness_casd(jet1subjettinesscasdInfo[ij].first_order_subtracted());
372  jet->GetShapeProperties()->SetSecondOrderSubtracted1subjettiness_casd(jet1subjettinesscasdInfo[ij].second_order_subtracted());
373  }
374 
375  std::vector<fastjet::contrib::GenericSubtractorInfo> jet2subjettinesscasdInfo = fjw.GetGenSubtractorInfoJet2subjettiness_casd();
376  Int_t n2subjettiness_casd = (Int_t)jet2subjettinesscasdInfo.size();
377  if(n2subjettiness_casd > ij && n2subjettiness_casd > 0) {
378  jet->GetShapeProperties()->SetFirstDerivative2subjettiness_casd(jet2subjettinesscasdInfo[ij].first_derivative());
379  jet->GetShapeProperties()->SetSecondDerivative2subjettiness_casd(jet2subjettinesscasdInfo[ij].second_derivative());
380  jet->GetShapeProperties()->SetFirstOrderSubtracted2subjettiness_casd(jet2subjettinesscasdInfo[ij].first_order_subtracted());
381  jet->GetShapeProperties()->SetSecondOrderSubtracted2subjettiness_casd(jet2subjettinesscasdInfo[ij].second_order_subtracted());
382  }
383 
384  std::vector<fastjet::contrib::GenericSubtractorInfo> jetOpeningAnglecasdInfo = fjw.GetGenSubtractorInfoJetOpeningAngle_casd();
385  Int_t nOpeningAngle_casd = (Int_t)jetOpeningAnglecasdInfo.size();
386  if(nOpeningAngle_casd > ij && nOpeningAngle_casd > 0) {
387  jet->GetShapeProperties()->SetFirstDerivativeOpeningAngle_casd(jetOpeningAnglecasdInfo[ij].first_derivative());
388  jet->GetShapeProperties()->SetSecondDerivativeOpeningAngle_casd(jetOpeningAnglecasdInfo[ij].second_derivative());
389  jet->GetShapeProperties()->SetFirstOrderSubtractedOpeningAngle_casd(jetOpeningAnglecasdInfo[ij].first_order_subtracted());
390  jet->GetShapeProperties()->SetSecondOrderSubtractedOpeningAngle_casd(jetOpeningAnglecasdInfo[ij].second_order_subtracted());
391  }
392  }
393 
394 #endif
395 }
396 
397 //______________________________________________________________________________
399 {
400  // Run termination of the utility (after each event).
401 }
void SetSecondOrderSubtractedOpeningAngle_akt02(Double_t d)
void SetSecondDerivative1subjettiness_ca(Double_t d)
void SetRMaxAndStep(Double_t rmax, Double_t dr)
Definition: AliFJWrapper.h:140
void SetUseExternalBkg(Bool_t b, Double_t rho, Double_t rhom)
Definition: AliFJWrapper.h:139
void SetSecondOrderSubtracted2subjettiness_casd(Double_t d)
void AddGRNumAt(Float_t num, Int_t idx)
void SetSecondOrderSubtracted1subjettiness_casd(Double_t d)
double Double_t
Definition: External.C:58
void SetFirstDerivativeOpeningAngle_ca(Double_t d)
void SetFirstOrderSubtracted2subjettiness_akt02(Double_t d)
void SetFirstDerivative1subjettiness_kt(Double_t d)
AliEmcalJetUtility & operator=(const AliEmcalJetUtility &jet)
void SetFirstOrderSubtractedAngularity(Double_t d)
virtual Int_t DoGenericSubtractionJetOpeningAngle_kt()
virtual Int_t DoGenericSubtractionJetpTD()
void SetFirstOrderSubtracted1subjettiness_kt(Double_t d)
void SetSecondOrderSubtractedOpeningAngle_kt(Double_t d)
void SetSecondDerivativeOpeningAngle_kt(Double_t d)
void SetFirstDerivativeOpeningAngle_akt02(Double_t d)
void SetSecondOrderSubtracted3subjettiness_kt(Double_t d)
void SetSecondOrderSubtracted1subjettiness_akt02(Double_t d)
virtual Int_t DoGenericSubtractionJetOpeningAngle_ca()
virtual Int_t DoGenericSubtractionJetAngularity()
virtual Int_t DoGenericSubtractionGR(Int_t ijet)
virtual Int_t DoGenericSubtractionJetMass()
Double_t GetRadius()
void SetFirstDerivative2subjettiness_akt02(Double_t d)
void SetSecondOrderSubtracted2subjettiness_akt02(Double_t d)
void SetFirstDerivative2subjettiness_ca(Double_t d)
virtual Int_t DoGenericSubtractionJet2subjettiness_casd()
void SetSecondDerivativeOpeningAngle_akt02(Double_t d)
void SetSecondOrderSubtractedOpeningAngle_ca(Double_t d)
void SetSecondOrderSubtractedCircularity(Double_t d)
void AddGRDenSubAt(Float_t den, Int_t idx)
virtual Int_t DoGenericSubtractionJetOpeningAngle_akt02()
void AddGRDenAt(Float_t den, Int_t idx)
void SetFirstOrderSubtractedOpeningAngle_akt02(Double_t d)
virtual Int_t DoGenericSubtractionJetOpeningAngle_casd()
int Int_t
Definition: External.C:63
void SetFirstDerivativeOpeningAngle_kt(Double_t d)
void ProcessJet(AliEmcalJet *jet, Int_t ij, AliFJWrapper &fjw)
void SetFirstOrderSubtracted3subjettiness_kt(Double_t d)
Double_t GetJetArea(UInt_t idx) const
void SetSecondOrderSubtracted2subjettiness_kt(Double_t d)
void SetFirstDerivative1subjettiness_akt02(Double_t d)
void AddGRNumSubAt(Float_t num, Int_t idx)
virtual Int_t DoGenericSubtractionJetCircularity()
void SetFirstOrderSubtracted1subjettiness_akt02(Double_t d)
void SetFirstDerivative2subjettiness_casd(Double_t d)
void SetSecondDerivativeOpeningAngle_casd(Double_t d)
void SetSecondDerivative1subjettiness_kt(Double_t d)
void SetFirstOrderSubtractedOpeningAngle_ca(Double_t d)
AliEmcalJetUtilityGenSubtractor & operator=(const AliEmcalJetUtilityGenSubtractor &jet)
void SetFirstOrderSubtracted2subjettiness_kt(Double_t d)
virtual std::vector< double > GetGRNumeratorSub() const
Definition: AliFJWrapper.h:85
void SetFirstOrderSubtracted1subjettiness_ca(Double_t d)
virtual Int_t DoGenericSubtractionJet1subjettiness_kt()
virtual Int_t DoGenericSubtractionJet1subjettiness_casd()
Double_t Pt() const
Definition: AliEmcalJet.h:102
void SetSecondDerivative2subjettiness_kt(Double_t d)
virtual Int_t DoGenericSubtractionJet2subjettiness_ca()
AliEmcalJetTask * fJetTask
virtual Int_t DoGenericSubtractionJet1subjettiness_akt02()
void SetFirstOrderSubtractedCircularity(Double_t d)
virtual Int_t DoGenericSubtractionJetSigma2()
void SetFirstOrderSubtracted1subjettiness_casd(Double_t d)
virtual std::vector< double > GetGRDenominatorSub() const
Definition: AliFJWrapper.h:86
virtual Int_t DoGenericSubtractionJet2subjettiness_akt02()
void SetFirstOrderSubtracted2subjettiness_casd(Double_t d)
void SetSecondDerivative1subjettiness_akt02(Double_t d)
void SetFirstOrderSubtractedOpeningAngle_casd(Double_t d)
void SetSecondOrderSubtracted1subjettiness_ca(Double_t d)
void SetSecondOrderSubtractedConstituent(Double_t d)
void SetSecondDerivative2subjettiness_ca(Double_t d)
virtual std::vector< double > GetGRDenominator() const
Definition: AliFJWrapper.h:84
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:44
void SetSecondDerivative2subjettiness_akt02(Double_t d)
void SetFirstDerivative1subjettiness_ca(Double_t d)
void SetSecondDerivative2subjettiness_casd(Double_t d)
void SetSecondDerivativeOpeningAngle_ca(Double_t d)
void SetFirstOrderSubtractedOpeningAngle_kt(Double_t d)
virtual Int_t DoGenericSubtractionJet3subjettiness_kt()
void SetFirstOrderSubtracted2subjettiness_ca(Double_t d)
void SetFirstDerivative3subjettiness_kt(Double_t d)
AliEmcalJetShapeProperties * GetShapeProperties() const
Definition: AliEmcalJet.h:256
virtual Int_t DoGenericSubtractionJet1subjettiness_ca()
void SetSecondOrderSubtractedOpeningAngle_casd(Double_t d)
virtual std::vector< double > GetGRNumerator() const
Definition: AliFJWrapper.h:83
virtual Int_t DoGenericSubtractionJetLeSub()
void SetFirstDerivative2subjettiness_kt(Double_t d)
virtual Int_t DoGenericSubtractionJetConstituent()
void SetFirstOrderSubtractedConstituent(Double_t d)
void SetSecondOrderSubtracted2subjettiness_ca(Double_t d)
void SetFirstDerivative1subjettiness_casd(Double_t d)
void SetSecondDerivative3subjettiness_kt(Double_t d)
virtual Int_t DoGenericSubtractionJet2subjettiness_kt()
void SetSecondOrderSubtracted1subjettiness_kt(Double_t d)
void SetFirstDerivativeOpeningAngle_casd(Double_t d)
void SetSecondDerivative1subjettiness_casd(Double_t d)
void SetSecondOrderSubtractedAngularity(Double_t d)