AliPhysics  v5-07-15-01 (b3d7633)
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros
AliRDHFCutsXictoeleXifromAODtracks.cxx
Go to the documentation of this file.
1 /**************************************************************************
2  * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
3  * *
4  * Author: The ALICE Off-line Project. *
5  * Contributors are mentioned in the code where appropriate. *
6  * *
7  * Permission to use, copy, modify and distribute this software and its *
8  * documentation strictly for non-commercial purposes is hereby granted *
9  * without fee, provided that the above copyright notice appears in all *
10  * copies and that both the copyright notice and this permission notice *
11  * appear in the supporting documentation. The authors make no claims *
12  * about the suitability of this software for any purpose. It is *
13  * provided "as is" without express or implied warranty. *
14  **************************************************************************/
15 
16 /* $Id$ */
17 
19 //
20 // Class for cuts on AOD reconstructed Xic->e+Xi
21 //
22 // Modified by Y.S Watanabe - wyosuke@cns.s.u-tokyo.ac.jp
23 //
25 
26 #include <Riostream.h>
27 
28 #include <TDatabasePDG.h>
29 #include <TMath.h>
30 
31 #include "AliAnalysisManager.h"
32 #include "AliInputEventHandler.h"
33 #include "AliPIDResponse.h"
34 #include "AliRDHFCutsXictoeleXifromAODtracks.h"
35 #include "AliAODRecoCascadeHF.h"
36 #include "AliAODTrack.h"
37 #include "AliESDtrack.h"
38 #include "AliESDVertex.h"
39 #include "AliAODVertex.h"
40 #include "AliAODv0.h"
41 #include "AliAODcascade.h"
42 #include "AliESDv0.h"
43 
44 using std::cout;
45 using std::endl;
46 
48 ClassImp(AliRDHFCutsXictoeleXifromAODtracks);
50 
51 
52 //--------------------------------------------------------------------------
53 AliRDHFCutsXictoeleXifromAODtracks::AliRDHFCutsXictoeleXifromAODtracks(const char* name) :
54 AliRDHFCuts(name),
55  fPIDStrategy(kNSigmaCuts),
56  fCombinedPIDThreshold(0.),
57  fUseCascadePID(kFALSE),
58  fPidObjCascPi(0),
59  fPidObjCascPr(0),
60  fProdTrackTPCNclsPIDMin(0),
61  fProdTrackTPCNclsRatioMin(0.0),
62  fProdUseAODFilterBit(kTRUE),
63  fProdMassTolLambda(0.010),
64  fProdMassTolXiRough(0.020),
65  fProdMassTolXi(0.008),
66  fProdMassRejOmega(0.008),
67  fProdRfidMinV0(0.6),
68  fProdRfidMaxV0(100.0),
69  fProdRfidMinXi(0.6),
70  fProdRfidMaxXi(100.0),
71  fProdCascProperDecayLengthMax(100.0),
72  fProdDcaXiDaughtersMax(100.),
73  fProdDcaV0DaughtersMax(100.),
74  fProdDcaBachToPrimVertexMin(0.),
75  fProdDcaV0ToPrimVertexMin(0.),
76  fProdDcaV0PrToPrimVertexMin(0.),
77  fProdDcaV0PiToPrimVertexMin(0.),
78  fProdXiCosineOfPoiningAngleMin(-1.),
79  fProdV0CosineOfPoiningAngleXiMin(-1.),
80  fProdCascNTPCClustersMin(0.0),
81  fProdCascEtaMin(-9999.),
82  fProdCascEtaMax(9999.),
83  fProdCascRapMin(-9999.),
84  fProdCascRapMax(9999.),
85  fProdRoughMassTol(0.25),
86  fProdRoughPtMin(0.0),
87  fExcludePionTPC(kFALSE),
88  fExcludeProtonTPC(kFALSE),
89  fExcludeKaonTPC(kFALSE),
90  fExcludenSigmaPionTPC(3.),
91  fExcludenSigmaProtonTPC(3.),
92  fExcludenSigmaKaonTPC(3.),
93  fSigmaElectronTPCMin(-9999.),
94  fSigmaElectronTPCPtDepPar0(-9999.),
95  fSigmaElectronTPCPtDepPar1(-9999.),
96  fSigmaElectronTPCPtDepPar2(0.),
97  fSigmaElectronTPCMax(9999.),
98  fSigmaElectronTOFMin(-9999.),
99  fSigmaElectronTOFMax(9999.)
100 {
101  //
102  // Default Constructor
103  //
104 
105  const Int_t nvars=2;
106  SetNVars(nvars);
107  TString varNames[nvars]={
108  "InvMass [GeV/c2]", //0
109  "cos(Opening angle) [cos(rad)]" //1
110  };
111 
112  Bool_t isUpperCut[nvars]={
113  kTRUE, //0
114  kFALSE //1
115  };
116  SetVarNames(nvars,varNames,isUpperCut);
117  Bool_t forOpt[nvars]={
118  kTRUE, //0
119  kTRUE, //1
120  };
121  SetVarsForOpt(nvars,forOpt);
122 
123  Float_t limits[2]={0,999999999.};
124  SetPtBins(2,limits);
125 }
126 //--------------------------------------------------------------------------
127 AliRDHFCutsXictoeleXifromAODtracks::AliRDHFCutsXictoeleXifromAODtracks(const AliRDHFCutsXictoeleXifromAODtracks &source) :
128  AliRDHFCuts(source),
129  fPIDStrategy(source.fPIDStrategy),
130  fCombinedPIDThreshold(source.fCombinedPIDThreshold),
131  fUseCascadePID(source.fUseCascadePID),
132  fPidObjCascPi(source.fPidObjCascPi),
133  fPidObjCascPr(source.fPidObjCascPr),
134  fProdTrackTPCNclsPIDMin(source.fProdTrackTPCNclsPIDMin),
135  fProdTrackTPCNclsRatioMin(source.fProdTrackTPCNclsRatioMin),
136  fProdUseAODFilterBit(source.fProdUseAODFilterBit),
137  fProdMassTolLambda(source.fProdMassTolLambda),
138  fProdMassTolXiRough(source.fProdMassTolXiRough),
139  fProdMassTolXi(source.fProdMassTolXi),
140  fProdMassRejOmega(source.fProdMassRejOmega),
141  fProdRfidMinV0(source.fProdRfidMinV0),
142  fProdRfidMaxV0(source.fProdRfidMaxV0),
143  fProdRfidMinXi(source.fProdRfidMinXi),
144  fProdRfidMaxXi(source.fProdRfidMaxXi),
145  fProdCascProperDecayLengthMax(source.fProdCascProperDecayLengthMax),
146  fProdDcaXiDaughtersMax(source.fProdDcaXiDaughtersMax),
147  fProdDcaV0DaughtersMax(source.fProdDcaV0DaughtersMax),
148  fProdDcaBachToPrimVertexMin(source.fProdDcaBachToPrimVertexMin),
149  fProdDcaV0ToPrimVertexMin(source.fProdDcaV0ToPrimVertexMin),
150  fProdDcaV0PrToPrimVertexMin(source.fProdDcaV0PrToPrimVertexMin),
151  fProdDcaV0PiToPrimVertexMin(source.fProdDcaV0PiToPrimVertexMin),
152  fProdXiCosineOfPoiningAngleMin(source.fProdXiCosineOfPoiningAngleMin),
153  fProdV0CosineOfPoiningAngleXiMin(source.fProdV0CosineOfPoiningAngleXiMin),
154  fProdCascNTPCClustersMin(source.fProdCascNTPCClustersMin),
155  fProdCascEtaMin(source.fProdCascEtaMin),
156  fProdCascEtaMax(source.fProdCascEtaMax),
157  fProdCascRapMin(source.fProdCascRapMin),
158  fProdCascRapMax(source.fProdCascRapMax),
159  fProdRoughMassTol(source.fProdRoughMassTol),
160  fProdRoughPtMin(source.fProdRoughPtMin),
161  fExcludePionTPC(source.fExcludePionTPC),
162  fExcludeProtonTPC(source.fExcludeProtonTPC),
163  fExcludeKaonTPC(source.fExcludeKaonTPC),
164  fExcludenSigmaPionTPC(source.fExcludenSigmaPionTPC),
165  fExcludenSigmaProtonTPC(source.fExcludenSigmaProtonTPC),
166  fExcludenSigmaKaonTPC(source.fExcludenSigmaKaonTPC),
167  fSigmaElectronTPCMin(source.fSigmaElectronTPCMin),
168  fSigmaElectronTPCPtDepPar0(source.fSigmaElectronTPCPtDepPar0),
169  fSigmaElectronTPCPtDepPar1(source.fSigmaElectronTPCPtDepPar1),
170  fSigmaElectronTPCPtDepPar2(source.fSigmaElectronTPCPtDepPar2),
171  fSigmaElectronTPCMax(source.fSigmaElectronTPCMax),
172  fSigmaElectronTOFMin(source.fSigmaElectronTOFMin),
173  fSigmaElectronTOFMax(source.fSigmaElectronTOFMax)
174 {
175  //
176  // Copy constructor
177  //
178 }
179 //--------------------------------------------------------------------------
180 AliRDHFCutsXictoeleXifromAODtracks &AliRDHFCutsXictoeleXifromAODtracks::operator=(const AliRDHFCutsXictoeleXifromAODtracks &source)
181 {
182  //
183  // assignment operator
184  //
185 
186  if (this != &source) {
187  AliRDHFCuts::operator=(source);
188  }
189 
190  fPIDStrategy = source.fPIDStrategy;
191  fCombinedPIDThreshold = source.fCombinedPIDThreshold;
192  fUseCascadePID = source.fUseCascadePID;
193  fPidObjCascPi = source.fPidObjCascPi;
194  fPidObjCascPr = source.fPidObjCascPr;
195  fProdTrackTPCNclsPIDMin = source.fProdTrackTPCNclsPIDMin;
196  fProdTrackTPCNclsRatioMin = source.fProdTrackTPCNclsRatioMin;
197  fProdUseAODFilterBit = source.fProdUseAODFilterBit;
198  fProdMassTolLambda = source.fProdMassTolLambda;
199  fProdMassTolXiRough = source.fProdMassTolXiRough;
200  fProdMassTolXi = source.fProdMassTolXi;
201  fProdMassRejOmega = source.fProdMassRejOmega;
202  fProdRfidMinV0 = source.fProdRfidMinV0;
203  fProdRfidMaxV0 = source.fProdRfidMaxV0;
204  fProdRfidMinXi = source.fProdRfidMinXi;
205  fProdRfidMaxXi = source.fProdRfidMaxXi;
206  fProdCascProperDecayLengthMax = source.fProdCascProperDecayLengthMax;
207  fProdDcaXiDaughtersMax = source.fProdDcaXiDaughtersMax;
208  fProdDcaV0DaughtersMax = source.fProdDcaV0DaughtersMax;
209  fProdDcaBachToPrimVertexMin = source.fProdDcaBachToPrimVertexMin;
210  fProdDcaV0ToPrimVertexMin = source.fProdDcaV0ToPrimVertexMin;
211  fProdDcaV0PrToPrimVertexMin = source.fProdDcaV0PrToPrimVertexMin;
212  fProdDcaV0PiToPrimVertexMin = source.fProdDcaV0PiToPrimVertexMin;
213  fProdXiCosineOfPoiningAngleMin = source.fProdXiCosineOfPoiningAngleMin;
214  fProdV0CosineOfPoiningAngleXiMin = source.fProdV0CosineOfPoiningAngleXiMin;
215  fProdCascNTPCClustersMin = source.fProdCascNTPCClustersMin;
216  fProdCascEtaMin = source.fProdCascEtaMin;
217  fProdCascEtaMax = source.fProdCascEtaMax;
218  fProdCascRapMin = source.fProdCascRapMin;
219  fProdCascRapMax = source.fProdCascRapMax;
220  fProdRoughMassTol = source.fProdRoughMassTol;
221  fProdRoughPtMin = source.fProdRoughPtMin;
222  fExcludePionTPC = source.fExcludePionTPC;
223  fExcludeProtonTPC = source.fExcludeProtonTPC;
224  fExcludeKaonTPC = source.fExcludeKaonTPC;
225  fExcludenSigmaPionTPC = source.fExcludenSigmaPionTPC;
226  fExcludenSigmaProtonTPC = source.fExcludenSigmaProtonTPC;
227  fExcludenSigmaKaonTPC = source.fExcludenSigmaKaonTPC;
228  fSigmaElectronTPCMin = source.fSigmaElectronTPCMin;
229  fSigmaElectronTPCPtDepPar0 = source.fSigmaElectronTPCPtDepPar0;
230  fSigmaElectronTPCPtDepPar1 = source.fSigmaElectronTPCPtDepPar1;
231  fSigmaElectronTPCPtDepPar2 = source.fSigmaElectronTPCPtDepPar2;
232  fSigmaElectronTPCMax = source.fSigmaElectronTPCMax;
233  fSigmaElectronTOFMin = source.fSigmaElectronTOFMin;
234  fSigmaElectronTOFMax = source.fSigmaElectronTOFMax;
235 
236 
237  return *this;
238 }
239 
240 //---------------------------------------------------------------------------
241 AliRDHFCutsXictoeleXifromAODtracks::~AliRDHFCutsXictoeleXifromAODtracks() {
242  //
243  // Default Destructor
244  //
245 
246 }
247 
248 //---------------------------------------------------------------------------
249 void AliRDHFCutsXictoeleXifromAODtracks::GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters) {
250  //
251  // Fills in vars the values of the variables
252  //
253 
254  if (pdgdaughters[0]==-9999) return; // dummy
255 
256  AliAODRecoCascadeHF* dd=(AliAODRecoCascadeHF*)d;
257  if(!dd){
258  AliError("No AliAODRecoCascadeHF object found\n");
259  return;
260  }
261 
262  if (nvars!=fnVarsForOpt) {
263  AliError("Cut object seems to have the wrong number of variables\n");
264  return;
265  }
266 
267  //Double_t ptD=d->Pt();
268  //Int_t ptbin=PtBin(ptD);
269  Int_t iter=-1;
270 
271  if(fVarsForOpt[0]){
272  iter++;
273  UInt_t pdgdg[2]={11,3312};
274  vars[iter]= dd->InvMass(2,pdgdg);
275  }
276  if(fVarsForOpt[1]){
277  iter++;
278  Double_t xipx = dd->PxProng(1);
279  Double_t xipy = dd->PyProng(1);
280  Double_t xipz = dd->PzProng(1);
281  Double_t epx = dd->PxProng(0);
282  Double_t epy = dd->PyProng(0);
283  Double_t epz = dd->PzProng(0);
284  vars[iter]= (xipx*epx+xipy*epy+xipz*epz)/sqrt(xipx*xipx+xipy*xipy+xipz*xipz)/sqrt(epx*epx+epy*epy+epz*epz);
285  }
286 
287  return;
288 }
289 
290 //---------------------------------------------------------------------------
291 Int_t AliRDHFCutsXictoeleXifromAODtracks::IsSelected(TObject* obj, Int_t selectionLevel) {
292  //
293  // Apply selection
294  //
295 
296  if (!fCutsRD) {
297  AliFatal("Cut matrix not inizialized. Exit...");
298  return 0;
299  }
300 
301  AliAODRecoCascadeHF* d=(AliAODRecoCascadeHF*)obj;
302  if(!d){
303  AliDebug(2,"AliAODRecoCascadeHF null");
304  return 0;
305  }
306 
307  Double_t ptD=d->Pt();
308  if(ptD<fMinPtCand) return 0;
309  if(ptD>fMaxPtCand) return 0;
310 
311  if (selectionLevel==AliRDHFCuts::kAll ||
312  selectionLevel==AliRDHFCuts::kTracks) {
313  //Performed in production stage
314  }
315 
316  Int_t returnvalueCuts=1;
317  // selection on candidate
318  if (selectionLevel==AliRDHFCuts::kAll ||
319  selectionLevel==AliRDHFCuts::kCandidate) {
320 
321  Double_t pt=d->Pt();
322  Int_t ptbin=PtBin(pt);
323  if (ptbin==-1) {
324  return 0;
325  }
326  Bool_t okcand=kTRUE;
327 
328  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
329  UInt_t pdgdg[2]={11,3312};
330  Double_t InvMassEleXi = d->InvMass(2,pdgdg);
331  Double_t xipx = d->PxProng(1);
332  Double_t xipy = d->PyProng(1);
333  Double_t xipz = d->PzProng(1);
334  Double_t epx = d->PxProng(0);
335  Double_t epy = d->PyProng(0);
336  Double_t epz = d->PzProng(0);
337  Double_t cosoa = (xipx*epx+xipy*epy+xipz*epz)/sqrt(xipx*xipx+xipy*xipy+xipz*xipz)/sqrt(epx*epx+epy*epy+epz*epz);
338 
339  if(InvMassEleXi > fCutsRD[GetGlobalIndex(0,ptbin)])
340  {
341  okcand = kFALSE;
342  }
343  if(cosoa < fCutsRD[GetGlobalIndex(1,ptbin)])
344  {
345  okcand = kFALSE;
346  }
347 
348  if(!okcand) return 0;
349  returnvalueCuts = 1;
350  }
351 
352  Int_t returnvaluePID=1;
353  if(selectionLevel==AliRDHFCuts::kAll ||
354  selectionLevel==AliRDHFCuts::kCandidate||
355  selectionLevel==AliRDHFCuts::kPID) {
356  //not used. applied at single level
357  }
358 
359  Int_t returnvalue = 0;
360  if(returnvalueCuts==1 && returnvaluePID==1) returnvalue=1;
361 
362  return returnvalue;
363 }
364 
365 //---------------------------------------------------------------------------
366 Int_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedPID(AliAODRecoDecayHF* obj)
367 {
368  //
369  // IsSelectedPID
370  //
371 
372  if(!fUsePID || !obj) return 1;
373 
374  AliAODRecoCascadeHF* dd=(AliAODRecoCascadeHF*)obj;
375  AliAODTrack *part = dd->GetBachelor();
376 
377  Int_t returnvalue=1;
378 
379  if(fPidHF->GetPidResponse()==0x0){
380  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
381  AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
382  AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
383  fPidHF->SetPidResponse(pidResp);
384  }
385 
386  Int_t isElectron=fPidHF->MakeRawPid(part,0);
387  if(isElectron<1) returnvalue = 0;
388 
389  return returnvalue;
390 }
391 
392 //---------------------------------------------------------------------------
393 Int_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedCombinedPID(AliAODRecoDecayHF* obj) {
394  //
395  // IsSelectedCombinedPID
396  //
397 
398  if(!fUsePID || !obj) {return 1;}
399 
400  AliAODRecoCascadeHF* dd=(AliAODRecoCascadeHF*)obj;
401  AliAODTrack *part = dd->GetBachelor();
402  if(!part) return 0;
403 
404  Int_t returnvalue=1;
405 
406  return returnvalue;
407 }
408 
409 //________________________________________________________________________
410 Bool_t AliRDHFCutsXictoeleXifromAODtracks::SingleTrkCuts(AliAODTrack *trk, AliAODVertex *primVert)
411 {
412  //
413  // Single Track Cut to be applied before object creation
414  //
415 
416  if(trk->GetStatus()&AliESDtrack::kITSpureSA) return kFALSE;
417  if(!(trk->GetStatus()&AliESDtrack::kITSin)) return kFALSE;
418  if(fProdUseAODFilterBit && !trk->TestFilterMask(BIT(4))) return kFALSE;
419 
420  Double_t pos[3]; primVert->GetXYZ(pos);
421  Double_t cov[6]; primVert->GetCovarianceMatrix(cov);
422  const AliESDVertex vESD(pos,cov,100.,100);
423  if(fTrackCuts&&!IsDaughterSelected(trk,&vESD,fTrackCuts)) return kFALSE;
424 
425  if(trk->GetTPCsignalN()<fProdTrackTPCNclsPIDMin) return kFALSE;
426  if(trk->GetTPCNclsF()>0){
427  Float_t tpcratio = (Float_t)trk->GetTPCncls()/(Float_t)trk->GetTPCNclsF();
428  if(tpcratio<fProdTrackTPCNclsRatioMin) return kFALSE;
429  }
430 
431  if(fUsePID)
432  {
433  if(fPidHF->GetPidResponse()==0x0){
434  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
435  AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
436  AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
437  fPidHF->SetPidResponse(pidResp);
438  }
439 
440  switch(fPIDStrategy){
441  case kNSigmaCuts:
442  return IsSelectedeID(trk);
443  break;
444  case kNSigmaCustomizedCuts:
445  return IsSelectedCustomizedeID(trk);
446  break;
447  case kNSigmaCustomizedPtDepCuts:
448  return IsSelectedCustomizedPtDepeID(trk);
449  break;
450  case kCombinedCuts:
451  return IsSelectedCombinedeID(trk);
452  break;
453  }
454 
455  }
456 
457  return kTRUE;
458 }
459 //________________________________________________________________________
460 Bool_t AliRDHFCutsXictoeleXifromAODtracks::SingleTrkCutsNoPID(AliAODTrack *trk, AliAODVertex *primVert)
461 {
462  //
463  // Single Track Cut to be applied before object creation
464  //
465 
466  if(trk->GetStatus()&AliESDtrack::kITSpureSA) return kFALSE;
467  if(!(trk->GetStatus()&AliESDtrack::kITSin)) return kFALSE;
468  if(fProdUseAODFilterBit && !trk->TestFilterMask(BIT(4))) return kFALSE;
469 
470  Double_t pos[3]; primVert->GetXYZ(pos);
471  Double_t cov[6]; primVert->GetCovarianceMatrix(cov);
472  const AliESDVertex vESD(pos,cov,100.,100);
473  if(fTrackCuts&&!IsDaughterSelected(trk,&vESD,fTrackCuts)) return kFALSE;
474 
475  if(trk->GetTPCsignalN()<fProdTrackTPCNclsPIDMin) return kFALSE;
476  if(trk->GetTPCNclsF()>0){
477  Float_t tpcratio = (Float_t)trk->GetTPCncls()/(Float_t)trk->GetTPCNclsF();
478  if(tpcratio<fProdTrackTPCNclsRatioMin) return kFALSE;
479  }
480 
481  return kTRUE;
482 }
483 //________________________________________________________________________
484 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedeID(AliAODTrack *trk)
485 {
486  //
487  // electron ID first shot
488  //
489 
490  Double_t nSigmaPion = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kPion);
491  if(fExcludePionTPC){
492  if(TMath::Abs(nSigmaPion)<fExcludenSigmaPionTPC){
493  return kFALSE;
494  }
495  }
496 
497  Double_t nSigmaProton = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kProton);
498  if(fExcludeProtonTPC){
499  if(TMath::Abs(nSigmaProton)<fExcludenSigmaProtonTPC){
500  return kFALSE;
501  }
502  }
503 
504  Double_t nSigmaKaon = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kKaon);
505  if(fExcludeKaonTPC){
506  if(TMath::Abs(nSigmaKaon)<fExcludenSigmaKaonTPC){
507  return kFALSE;
508  }
509  }
510 
511  Int_t isElectron=fPidHF->MakeRawPid(trk,0);
512  if(isElectron<1) return kFALSE;
513 
514  return kTRUE;
515 }
516 //________________________________________________________________________
517 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedCustomizedeID(AliAODTrack *trk)
518 {
519  //
520  // electron ID first shot
521  //
522 
523  Double_t nSigmaPion = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kPion);
524  if(fExcludePionTPC){
525  if(TMath::Abs(nSigmaPion)<fExcludenSigmaPionTPC){
526  return kFALSE;
527  }
528  }
529 
530  Double_t nSigmaProton = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kProton);
531  if(fExcludeProtonTPC){
532  if(TMath::Abs(nSigmaProton)<fExcludenSigmaProtonTPC){
533  return kFALSE;
534  }
535  }
536 
537  Double_t nSigmaKaon = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kKaon);
538  if(fExcludeKaonTPC){
539  if(TMath::Abs(nSigmaKaon)<fExcludenSigmaKaonTPC){
540  return kFALSE;
541  }
542  }
543 
544  Double_t nSigmaTPCele = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kElectron);
545  Double_t nSigmaTOFele = fPidHF->GetPidResponse()->NumberOfSigmasTOF(trk,AliPID::kElectron);
546 
547  if(nSigmaTPCele<fSigmaElectronTPCMin) return kFALSE;
548  if(nSigmaTPCele>fSigmaElectronTPCMax) return kFALSE;
549  if(fabs(fSigmaElectronTOFMin)<10.&&fabs(fSigmaElectronTOFMax)<10.){
550  if(nSigmaTOFele<fSigmaElectronTOFMin) return kFALSE;
551  if(nSigmaTOFele>fSigmaElectronTOFMax) return kFALSE;
552  }
553 
554  return kTRUE;
555 }
556 //________________________________________________________________________
557 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedCustomizedPtDepeID(AliAODTrack *trk)
558 {
559  //
560  // electron ID pt dependent
561  //
562 
563  Double_t nSigmaTPCele = fPidHF->GetPidResponse()->NumberOfSigmasTPC(trk,AliPID::kElectron);
564  Double_t nSigmaTOFele = fPidHF->GetPidResponse()->NumberOfSigmasTOF(trk,AliPID::kElectron);
565 
566  if(nSigmaTOFele<fSigmaElectronTOFMin) return kFALSE;
567  if(nSigmaTOFele>fSigmaElectronTOFMax) return kFALSE;
568 
569  Double_t pte = trk->Pt();
570  Double_t nsigmamin = fSigmaElectronTPCPtDepPar0+fSigmaElectronTPCPtDepPar1*pte+fSigmaElectronTPCPtDepPar2*pte*pte;
571  if(pte>5.) nsigmamin = fSigmaElectronTPCPtDepPar0+fSigmaElectronTPCPtDepPar1*5.+fSigmaElectronTPCPtDepPar2*25.;
572  if(nSigmaTPCele<nsigmamin) return kFALSE;
573  if(nSigmaTPCele>fSigmaElectronTPCMax) return kFALSE;
574 
575  return kTRUE;
576 }
577 //________________________________________________________________________
578 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSelectedCombinedeID(AliAODTrack *trk)
579 {
580  //
581  // electron ID Basyian not implemented
582  //
583 
584  return kTRUE;
585 }
586 //________________________________________________________________________
587 Bool_t AliRDHFCutsXictoeleXifromAODtracks::SingleCascadeCuts(AliAODcascade *casc,Double_t *primvert)
588 {
589  //
590  // Single Cascade Cut
591  //
592 
593  if(!casc) return kFALSE;
594  if(!(casc->GetSecondaryVtx())) return kFALSE;
595  if(!(casc->GetDecayVertexXi())) return kFALSE;
596 
597  AliAODTrack *ptrack = (AliAODTrack*) (casc->GetDaughter(0));
598  AliAODTrack *ntrack = (AliAODTrack*) (casc->GetDaughter(1));
599  AliAODTrack *btrack = (AliAODTrack*) (casc->GetDecayVertexXi()->GetDaughter(0));
600 
601  if(!ptrack||!ntrack||!btrack) return kFALSE;
602 
603  if(ptrack->Charge()<0 && ntrack->Charge()>0){
604  ptrack = (AliAODTrack*) (casc->GetDaughter(1));
605  ntrack = (AliAODTrack*) (casc->GetDaughter(0));
606  }
607 
608  if(ptrack->GetTPCClusterInfo(2,1)<fProdCascNTPCClustersMin) return kFALSE;
609  if(ntrack->GetTPCClusterInfo(2,1)<fProdCascNTPCClustersMin) return kFALSE;
610  if(btrack->GetTPCClusterInfo(2,1)<fProdCascNTPCClustersMin) return kFALSE;
611 
612 
613  Double_t mLPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
614  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
615  Double_t momegaPDG = TDatabasePDG::Instance()->GetParticle(3334)->Mass();
616 
617  Double_t massLambda = casc->MassLambda();
618  Double_t massAntiLambda = casc->MassAntiLambda();
619  if(TMath::Abs(massLambda-mLPDG)>fProdMassTolLambda && TMath::Abs(massAntiLambda-mLPDG)>fProdMassTolLambda)
620  return kFALSE;
621 
622  Bool_t isparticle = kTRUE;
623  //if(TMath::Abs(massAntiLambda-mLPDG)<fProdMassTolLambda) isparticle = kFALSE;
624  if(casc->ChargeXi()>0) isparticle = kFALSE;
625 
626  Double_t massXi = casc->MassXi();
627  if(TMath::Abs(massXi-mxiPDG)>fProdMassTolXiRough)
628  return kFALSE;
629 
630  Double_t massOmega = casc->MassOmega();
631  if(TMath::Abs(massOmega-momegaPDG)<fProdMassRejOmega)
632  return kFALSE;
633 
634 
635  Double_t lPosXi[3];
636  lPosXi[0] = casc->DecayVertexXiX();
637  lPosXi[1] = casc->DecayVertexXiY();
638  lPosXi[2] = casc->DecayVertexXiZ();
639  Double_t decayvertXi = TMath::Sqrt(lPosXi[0]*lPosXi[0]+lPosXi[1]*lPosXi[1]);
640  Double_t lPosV0[3];
641  lPosV0[0] = casc->DecayVertexV0X();
642  lPosV0[1] = casc->DecayVertexV0Y();
643  lPosV0[2] = casc->DecayVertexV0Z();
644  Double_t decayvertV0 = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
645 
646  if(decayvertV0<fProdRfidMinV0 || decayvertV0>fProdRfidMaxV0) return kFALSE;
647  if(decayvertXi<fProdRfidMinXi || decayvertXi>fProdRfidMaxXi) return kFALSE;
648 
649  Double_t ptotxi = TMath::Sqrt(pow(casc->MomXiX(),2)+pow(casc->MomXiY(),2)+pow(casc->MomXiZ(),2));
650  Double_t properdl = casc->DecayLengthXi(primvert[0],primvert[1],primvert[2])*mxiPDG/ptotxi;
651  if(properdl>fProdCascProperDecayLengthMax) return kFALSE;
652 
653  Double_t lDcaXiDaughters = casc->DcaXiDaughters();
654  Double_t lDcaV0Daughters = casc->DcaV0Daughters();
655  if(lDcaXiDaughters > fProdDcaXiDaughtersMax) return kFALSE;
656  if(lDcaV0Daughters > fProdDcaV0DaughtersMax) return kFALSE;
657 
658  Double_t lDcaBachToPrimVertex = casc->DcaBachToPrimVertex();
659  Double_t lDcaV0ToPrimVertex = casc->DcaV0ToPrimVertex();
660  Double_t lDcaPosToPrimVertex = casc->DcaPosToPrimVertex();
661  Double_t lDcaNegToPrimVertex = casc->DcaNegToPrimVertex();
662  if(lDcaBachToPrimVertex < fProdDcaBachToPrimVertexMin) return kFALSE;
663  if(lDcaV0ToPrimVertex < fProdDcaV0ToPrimVertexMin) return kFALSE;
664  if(isparticle){
665  if(lDcaPosToPrimVertex < fProdDcaV0PrToPrimVertexMin) return kFALSE;
666  if(lDcaNegToPrimVertex < fProdDcaV0PiToPrimVertexMin) return kFALSE;
667  }else{
668  if(lDcaPosToPrimVertex < fProdDcaV0PiToPrimVertexMin) return kFALSE;
669  if(lDcaNegToPrimVertex < fProdDcaV0PrToPrimVertexMin) return kFALSE;
670  }
671 
672  Double_t lXiCosineOfPointingAngle = casc->CosPointingAngleXi(primvert[0],primvert[1],primvert[2]);
673  Double_t lV0CosineOfPointingAngleXi = casc->CosPointingAngle(lPosXi);
674 
675  if(lXiCosineOfPointingAngle < fProdXiCosineOfPoiningAngleMin) return kFALSE;
676  if(lV0CosineOfPointingAngleXi < fProdV0CosineOfPoiningAngleXiMin) return kFALSE;
677 
678  if(fUseCascadePID)
679  {
680  if(fPidObjCascPi->GetPidResponse()==0x0){
681  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
682  AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
683  AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
684  fPidObjCascPi->SetPidResponse(pidResp);
685  }
686  if(fPidObjCascPr->GetPidResponse()==0x0){
687  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
688  AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
689  AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
690  fPidObjCascPr->SetPidResponse(pidResp);
691  }
692  if(isparticle){
693  Int_t isProton=fPidObjCascPr->MakeRawPid(ptrack,4);
694  Int_t isPion1 =fPidObjCascPi->MakeRawPid(btrack,2);
695  Int_t isPion2 =fPidObjCascPi->MakeRawPid(ntrack,2);
696  if(isProton<1) return kFALSE;
697  if(isPion1<1) return kFALSE;
698  if(isPion2<1) return kFALSE;
699  }else{
700  Int_t isProton=fPidObjCascPr->MakeRawPid(ntrack,4);
701  Int_t isPion1 =fPidObjCascPi->MakeRawPid(btrack,2);
702  Int_t isPion2 =fPidObjCascPi->MakeRawPid(ptrack,2);
703  if(isProton<1) return kFALSE;
704  if(isPion1<1) return kFALSE;
705  if(isPion2<1) return kFALSE;
706  }
707  }
708 
709  Double_t RapXi = casc->RapXi();
710  if(RapXi<fProdCascRapMin || RapXi>fProdCascRapMax) return kFALSE;
711 
712  Double_t EtaXi = 0.5*TMath::Log((ptotxi+casc->MomXiZ())/(ptotxi-casc->MomXiZ()));
713  if(EtaXi<fProdCascEtaMin || EtaXi>fProdCascEtaMax) return kFALSE;
714 
715  return kTRUE;
716 }
717 
718 //________________________________________________________________________
719 Bool_t AliRDHFCutsXictoeleXifromAODtracks::SelectWithRoughCuts(AliAODcascade *casc, AliAODTrack *part)
720 {
721  //
722  // Mass and pT Cut to be applied before object creation
723  // Not used now
724  //
725  if(!casc) return kFALSE;
726  if(!part) return kFALSE;
727 
728  return kTRUE;
729 }
730 
731 //________________________________________________________________________
732 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsPeakRegion(AliAODcascade *casc)
733 {
734  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
735  Double_t massXi = casc->MassXi();
736  if(TMath::Abs(massXi-mxiPDG)<fProdMassTolXi)
737  return kTRUE;
738  return kFALSE;
739 }
740 
741 //________________________________________________________________________
742 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsPeakRegion(TLorentzVector *casc)
743 {
744  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
745  Double_t massXi = casc->M();
746  if(TMath::Abs(massXi-mxiPDG)<fProdMassTolXi)
747  return kTRUE;
748  return kFALSE;
749 }
750 
751 //________________________________________________________________________
752 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSideBand(AliAODcascade *casc)
753 {
754  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
755  Double_t massXi = casc->MassXi();
756  Bool_t issideband = kFALSE;
757  if((massXi-mxiPDG)>fProdMassTolXiRough-fProdMassTolXi) issideband = kTRUE;
758  if((massXi-mxiPDG)<-fProdMassTolXiRough+fProdMassTolXi) issideband = kTRUE;
759  return issideband;
760 }
761 
762 //________________________________________________________________________
763 Bool_t AliRDHFCutsXictoeleXifromAODtracks::IsSideBand(TLorentzVector *casc)
764 {
765  Double_t mxiPDG = TDatabasePDG::Instance()->GetParticle(3312)->Mass();
766  Double_t massXi = casc->M();
767  Bool_t issideband = kFALSE;
768  if((massXi-mxiPDG)>fProdMassTolXiRough-fProdMassTolXi) issideband = kTRUE;
769  if((massXi-mxiPDG)<-fProdMassTolXiRough+fProdMassTolXi) issideband = kTRUE;
770  return issideband;
771 }
772 
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTPCMax
Double_t fSigmaElectronTPCMax
nSigma electron lower limit (par2)
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:189
AliRDHFCutsXictoeleXifromAODtracks::AliRDHFCutsXictoeleXifromAODtracks
AliRDHFCutsXictoeleXifromAODtracks(const char *name="CutsXictoeleXi")
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:53
ClassImp
ClassImp(AliAnalysisTaskTriggerRates) AliAnalysisTaskTriggerRates
Definition: AliAnalysisTaskTriggerRates.cxx:34
AliRDHFCutsXictoeleXifromAODtracks::fExcludePionTPC
Bool_t fExcludePionTPC
pT cut for Lc used before object creation
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:179
AliRDHFCutsXictoeleXifromAODtracks::fProdRfidMaxV0
Double_t fProdRfidMaxV0
Minimum Decay vertex of V0.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:159
AliRDHFCutsXictoeleXifromAODtracks::SingleCascadeCuts
Bool_t SingleCascadeCuts(AliAODcascade *casc, Double_t *vert)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:587
AliRDHFCutsXictoeleXifromAODtracks::operator=
AliRDHFCutsXictoeleXifromAODtracks & operator=(const AliRDHFCutsXictoeleXifromAODtracks &source)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:180
AliRDHFCutsXictoeleXifromAODtracks::fProdDcaV0PrToPrimVertexMin
Double_t fProdDcaV0PrToPrimVertexMin
Min Dca between v0 and PV.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:167
AliRDHFCutsXictoeleXifromAODtracks::fExcludeKaonTPC
Bool_t fExcludeKaonTPC
Flag wheter to exlude proton band.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:181
AliRDHFCutsXictoeleXifromAODtracks::fProdV0CosineOfPoiningAngleXiMin
Double_t fProdV0CosineOfPoiningAngleXiMin
Min Xi cos pointing angle to PV.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:170
AliRDHFCutsXictoeleXifromAODtracks::fExcludenSigmaPionTPC
Double_t fExcludenSigmaPionTPC
Flag wheter to exlude proton band.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:182
AliRDHFCutsXictoeleXifromAODtracks::GetCutVarsForOpt
virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d, Float_t *vars, Int_t nvars, Int_t *pdgdaughters)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:249
AliRDHFCutsXictoeleXifromAODtracks::fProdCascEtaMin
Double_t fProdCascEtaMin
Minimum number of TPC clusters.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:172
AliAODRecoCascadeHF.h
AliRDHFCutsXictoeleXifromAODtracks::fProdMassTolXi
Double_t fProdMassTolXi
Tolerance of Xi mass from PDG value (including sideband)
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:156
AliRDHFCutsXictoeleXifromAODtracks::fProdDcaBachToPrimVertexMin
Double_t fProdDcaBachToPrimVertexMin
Max Dca between V0 daughters.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:165
AliRDHFCutsXictoeleXifromAODtracks::fProdMassRejOmega
Double_t fProdMassRejOmega
Tolerance of Xi mass from PDG value.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:157
AliRDHFCutsXictoeleXifromAODtracks::fProdTrackTPCNclsRatioMin
Double_t fProdTrackTPCNclsRatioMin
Min. Number of TPC PID cluster.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:152
AliRDHFCutsXictoeleXifromAODtracks::fProdDcaV0DaughtersMax
Double_t fProdDcaV0DaughtersMax
Max Dca between Xi daughters.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:164
AliRDHFCutsXictoeleXifromAODtracks::IsSelectedPID
virtual Int_t IsSelectedPID(AliAODRecoDecayHF *obj)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:366
AliRDHFCutsXictoeleXifromAODtracks::fExcludenSigmaKaonTPC
Double_t fExcludenSigmaKaonTPC
nSigma to exclude for proton band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:184
AliRDHFCutsXictoeleXifromAODtracks::fProdRfidMaxXi
Double_t fProdRfidMaxXi
Minimum Decay vertex of Xi.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:161
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTPCMin
Double_t fSigmaElectronTPCMin
nSigma to exclude for Kaon band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:185
AliRDHFCutsXictoeleXifromAODtracks::fProdDcaV0ToPrimVertexMin
Double_t fProdDcaV0ToPrimVertexMin
Min Dca between Bachelor and PV.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:166
AliRDHFCutsXictoeleXifromAODtracks::fExcludeProtonTPC
Bool_t fExcludeProtonTPC
Flag wheter to exlude pion band.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:180
AliRDHFCuts::SetNVars
void SetNVars(Int_t nVars)
Definition: AliRDHFCuts.h:361
AliRDHFCutsXictoeleXifromAODtracks::fProdRoughMassTol
Double_t fProdRoughMassTol
Maximum rapidity of cascade.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:176
AliRDHFCutsXictoeleXifromAODtracks::fPidObjCascPr
AliAODPidHF * fPidObjCascPr
PID object for cascade-pion.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:148
AliRDHFCuts::operator=
AliRDHFCuts & operator=(const AliRDHFCuts &source)
Definition: AliRDHFCuts.cxx:199
AliRDHFCuts::fUsePID
Bool_t fUsePID
Definition: AliRDHFCuts.h:388
AliRDHFCutsXictoeleXifromAODtracks::fUseCascadePID
Bool_t fUseCascadePID
Threshold used in IsSelectedCombinedPID.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:146
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTOFMin
Double_t fSigmaElectronTOFMin
nSigma to exclude for Kaon band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:190
AliRDHFCutsXictoeleXifromAODtracks::fProdDcaV0PiToPrimVertexMin
Double_t fProdDcaV0PiToPrimVertexMin
Min Dca between v0-proton and PV.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:168
AliRDHFCutsXictoeleXifromAODtracks::fProdMassTolLambda
Double_t fProdMassTolLambda
Flag for AOD filter Bit used before object creation.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:154
AliRDHFCutsXictoeleXifromAODtracks::SingleTrkCuts
Bool_t SingleTrkCuts(AliAODTrack *trk, AliAODVertex *primvert)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:410
AliRDHFCuts::fnVarsForOpt
Int_t fnVarsForOpt
Definition: AliRDHFCuts.h:383
AliRDHFCutsXictoeleXifromAODtracks::SelectWithRoughCuts
Bool_t SelectWithRoughCuts(AliAODcascade *casc, AliAODTrack *trk1)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:719
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTPCPtDepPar1
Double_t fSigmaElectronTPCPtDepPar1
nSigma electron lower limit (par0)
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:187
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTPCPtDepPar0
Double_t fSigmaElectronTPCPtDepPar0
nSigma to exclude for Kaon band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:186
AliRDHFCuts::fMaxPtCand
Double_t fMaxPtCand
minimum pt of the candidate
Definition: AliRDHFCuts.h:410
AliRDHFCuts::fTrackCuts
AliESDtrackCuts * fTrackCuts
quality cuts on the daughter tracks
Definition: AliRDHFCuts.h:376
AliRDHFCutsXictoeleXifromAODtracks::IsSelected
virtual Int_t IsSelected(TObject *obj, Int_t selectionLevel)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:291
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTPCPtDepPar2
Double_t fSigmaElectronTPCPtDepPar2
nSigma electron lower limit (par1)
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:188
AliRDHFCuts::IsDaughterSelected
Bool_t IsDaughterSelected(AliAODTrack *track, const AliESDVertex *primary, AliESDtrackCuts *cuts) const
Definition: AliRDHFCuts.cxx:698
AliAODRecoCascadeHF::GetBachelor
AliAODTrack * GetBachelor() const
Definition: AliAODRecoCascadeHF.h:52
AliRDHFCutsXictoeleXifromAODtracks::fSigmaElectronTOFMax
Double_t fSigmaElectronTOFMax
nSigma to exclude for Kaon band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:191
AliRDHFCutsXictoeleXifromAODtracks::fProdTrackTPCNclsPIDMin
Int_t fProdTrackTPCNclsPIDMin
PID object for cascade-proton.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:151
AliAODPidHF::MakeRawPid
Int_t MakeRawPid(AliAODTrack *track, Int_t specie)
Definition: AliAODPidHF.cxx:735
AliAODPidHF::GetPidResponse
AliPIDResponse * GetPidResponse() const
Definition: AliAODPidHF.h:160
AliRDHFCuts::fCutsRD
Float_t * fCutsRD
fnVars*fnPtBins
Definition: AliRDHFCuts.h:386
AliRDHFCutsXictoeleXifromAODtracks::fProdMassTolXiRough
Double_t fProdMassTolXiRough
Tolerance of Lambda mass from PDG value.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:155
AliRDHFCutsXictoeleXifromAODtracks::fProdXiCosineOfPoiningAngleMin
Double_t fProdXiCosineOfPoiningAngleMin
Min Dca between v0-pion and PV.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:169
AliRDHFCutsXictoeleXifromAODtracks::fProdCascProperDecayLengthMax
Double_t fProdCascProperDecayLengthMax
Max Decay vertex of Xi.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:162
AliRDHFCuts::SetVarsForOpt
void SetVarsForOpt(Int_t nVars, Bool_t *forOpt)
Definition: AliRDHFCuts.cxx:912
AliRDHFCuts::SetVarNames
void SetVarNames(Int_t nVars, TString *varNames, Bool_t *isUpperCut)
Definition: AliRDHFCuts.cxx:889
AliRDHFCuts::fVarsForOpt
Bool_t * fVarsForOpt
number of cut vars to be optimized for candidates
Definition: AliRDHFCuts.h:384
AliRDHFCutsXictoeleXifromAODtracks::fProdCascNTPCClustersMin
Double_t fProdCascNTPCClustersMin
Min V0 cos pointing angle to Xi vertex.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:171
AliRDHFCutsXictoeleXifromAODtracks.h
AliRDHFCutsXictoeleXifromAODtracks::fProdCascRapMax
Double_t fProdCascRapMax
Minimum rapidity of cascade.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:175
AliRDHFCutsXictoeleXifromAODtracks::fProdUseAODFilterBit
Bool_t fProdUseAODFilterBit
Min. Number of TPC PID cluster.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:153
AliRDHFCutsXictoeleXifromAODtracks::fPidObjCascPi
AliAODPidHF * fPidObjCascPi
Use PID for cascade or not.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:147
AliRDHFCutsXictoeleXifromAODtracks::fProdRfidMinV0
Double_t fProdRfidMinV0
Rejection range of Omega mass from PDG value.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:158
AliRDHFCuts::SetPtBins
void SetPtBins(Int_t nPtBinLimits, Float_t *ptBinLimits)
Definition: AliRDHFCuts.cxx:866
AliRDHFCutsXictoeleXifromAODtracks::fExcludenSigmaProtonTPC
Double_t fExcludenSigmaProtonTPC
nSigma to exclude for pion band
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:183
AliRDHFCutsXictoeleXifromAODtracks::SingleTrkCutsNoPID
Bool_t SingleTrkCutsNoPID(AliAODTrack *trk, AliAODVertex *primvert)
Definition: AliRDHFCutsXictoeleXifromAODtracks.cxx:460
AliRDHFCutsXictoeleXifromAODtracks::fProdRoughPtMin
Double_t fProdRoughPtMin
Mass cut for Lc used before object creation.
Definition: AliRDHFCutsXictoeleXifromAODtracks.h:177
AliRDHFCuts::fPidHF
AliAODPidHF * fPidHF
enable AOD049 centrality cleanup
Definition: AliRDHFCuts.h:390
AliRDHFCuts::PtBin
Int_t PtBin(Double_t pt) const
Definition: AliRDHFCuts.cxx:1139
AliRDHFCuts::GetGlobalIndex
Int_t GetGlobalIndex(Int_t iVar, Int_t iPtBin) const
Definition: AliRDHFCuts.cxx:1120
AliRDHFCuts::fMinPtCand
Double_t fMinPtCand
outcome of PID selection
Definition: AliRDHFCuts.h:409
AliAODPidHF::SetPidResponse
void SetPidResponse(AliPIDResponse *pidResp)
Definition: AliAODPidHF.h:111