Skip to content

Commit 0773808

Browse files
lhusovaLucia Anna Tarasovicova
lhusova
and
Lucia Anna Tarasovicova
authored
event sel for MC (#6253)
Co-authored-by: Lucia Anna Tarasovicova <lucia.anna.husova@cern.ch>
1 parent c5f712d commit 0773808

1 file changed

Lines changed: 65 additions & 44 deletions

File tree

‎PWGLF/Tasks/Strangeness/derivedcascadeanalysis.cxx‎

Lines changed: 65 additions & 44 deletions
Original file line numberDiff line numberDiff line change
@@ -359,65 +359,86 @@ struct derivedCascadeAnalysis {
359359
histos.fill(HIST("hEventPVcontributorsVsCentralityBefCuts"), coll.centFT0C(), coll.multNTracksPVeta1());
360360
histos.fill(HIST("hEventGlobalTracksVsCentralityBefCuts"), coll.centFT0C(), coll.multNTracksGlobal());
361361
}
362+
if (isMC) {
363+
if (!coll.selection_bit(aod::evsel::kIsTriggerTVX)) {
364+
return false;
365+
}
366+
histos.fill(HIST("hEventSelection"), 1.5 /* MB trigger*/);
362367

363-
if (!sel) {
364-
return false;
365-
}
366-
histos.fill(HIST("hEventSelection"), 1.5 /* collisions after sel*/);
367-
if (TMath::Abs(coll.posZ()) > zVertexCut) {
368-
return false;
369-
}
368+
if (TMath::Abs(coll.posZ()) > zVertexCut) {
369+
return false;
370+
}
370371

371-
histos.fill(HIST("hEventSelection"), 2.5 /* collisions after sel pvz sel*/);
372+
histos.fill(HIST("hEventSelection"), 2.5 /* collisions after sel pvz sel*/);
372373

373-
if (coll.centFT0C() > centMax || coll.centFT0C() < centMin) {
374-
return false;
375-
}
374+
if (doTFeventCut && !coll.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
375+
return false;
376+
}
377+
histos.fill(HIST("hEventSelection"), 3.5 /* Not at TF border */);
378+
if (doITSFrameBorderCut && !coll.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
379+
return false;
380+
}
381+
histos.fill(HIST("hEventSelection"), 4.5 /* Not at ITS ROF border */);
382+
} else {
383+
if (!sel) {
384+
return false;
385+
}
386+
histos.fill(HIST("hEventSelection"), 1.5 /* collisions after sel*/);
387+
if (TMath::Abs(coll.posZ()) > zVertexCut) {
388+
return false;
389+
}
376390

377-
histos.fill(HIST("hEventSelection"), 3.5 /* collisions after centrality sel*/);
391+
histos.fill(HIST("hEventSelection"), 2.5 /* collisions after sel pvz sel*/);
378392

379-
if (doSameBunchPileUpEventCut && !coll.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
380-
return false;
381-
}
382-
histos.fill(HIST("hEventSelection"), 4.5 /* Not same Bunch pile up */);
393+
if (coll.centFT0C() > centMax || coll.centFT0C() < centMin) {
394+
return false;
395+
}
383396

384-
if (doGoodPVFT0EventCut && !coll.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
385-
return false;
386-
}
387-
histos.fill(HIST("hEventSelection"), 5.5 /* No large vertexZ difference from tracks and FT0*/);
397+
histos.fill(HIST("hEventSelection"), 3.5 /* collisions after centrality sel*/);
388398

389-
if (doITSTPCvertexEventCut && !coll.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) {
390-
return false;
391-
}
392-
histos.fill(HIST("hEventSelection"), 6.5 /* At least one ITS-TPC track in the event*/);
399+
if (doSameBunchPileUpEventCut && !coll.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
400+
return false;
401+
}
402+
histos.fill(HIST("hEventSelection"), 4.5 /* Not same Bunch pile up */);
393403

394-
if (doVertexTOFmatch && !coll.selection_bit(o2::aod::evsel::kIsVertexTOFmatched)) {
395-
return false;
396-
}
397-
histos.fill(HIST("hEventSelection"), 7.5 /* At least one of vertex contributors is matched to TOF*/);
404+
if (doGoodPVFT0EventCut && !coll.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
405+
return false;
406+
}
407+
histos.fill(HIST("hEventSelection"), 5.5 /* No large vertexZ difference from tracks and FT0*/);
398408

399-
if (doVertexTRDmatch && !coll.selection_bit(o2::aod::evsel::kIsVertexTRDmatched)) {
400-
return false;
401-
}
402-
histos.fill(HIST("hEventSelection"), 8.5 /* At least one of vertex contributors is matched to TRD*/);
409+
if (doITSTPCvertexEventCut && !coll.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) {
410+
return false;
411+
}
412+
histos.fill(HIST("hEventSelection"), 6.5 /* At least one ITS-TPC track in the event*/);
403413

404-
if (doITSFrameBorderCut && !coll.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
405-
return false;
406-
}
407-
histos.fill(HIST("hEventSelection"), 9.5 /* Not at ITS ROF border */);
414+
if (doVertexTOFmatch && !coll.selection_bit(o2::aod::evsel::kIsVertexTOFmatched)) {
415+
return false;
416+
}
417+
histos.fill(HIST("hEventSelection"), 7.5 /* At least one of vertex contributors is matched to TOF*/);
408418

409-
if (doTFeventCut && !coll.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
410-
return false;
411-
}
412-
histos.fill(HIST("hEventSelection"), 10.5 /* Not at TF border */);
419+
if (doVertexTRDmatch && !coll.selection_bit(o2::aod::evsel::kIsVertexTRDmatched)) {
420+
return false;
421+
}
422+
histos.fill(HIST("hEventSelection"), 8.5 /* At least one of vertex contributors is matched to TRD*/);
413423

414-
if (doMultiplicityCorrCut) {
415-
if (coll.multNTracksGlobal() < (1343.3 * TMath::Exp(-0.0443259 * coll.centFT0C()) - 50) || coll.multNTracksGlobal() > (2098.9 * TMath::Exp(-0.0332444 * coll.centFT0C())))
424+
if (doITSFrameBorderCut && !coll.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
416425
return false;
417-
if (coll.multNTracksPVeta1() < (3703 * TMath::Exp(-0.0455483 * coll.centFT0C()) - 150) || coll.multNTracksPVeta1() > (4937.33 * TMath::Exp(-0.0372668 * coll.centFT0C()) + 20))
426+
}
427+
histos.fill(HIST("hEventSelection"), 9.5 /* Not at ITS ROF border */);
428+
429+
if (doTFeventCut && !coll.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
418430
return false;
431+
}
432+
histos.fill(HIST("hEventSelection"), 10.5 /* Not at TF border */);
433+
434+
if (doMultiplicityCorrCut) {
435+
if (coll.multNTracksGlobal() < (1343.3 * TMath::Exp(-0.0443259 * coll.centFT0C()) - 50) || coll.multNTracksGlobal() > (2098.9 * TMath::Exp(-0.0332444 * coll.centFT0C())))
436+
return false;
437+
if (coll.multNTracksPVeta1() < (3703 * TMath::Exp(-0.0455483 * coll.centFT0C()) - 150) || coll.multNTracksPVeta1() > (4937.33 * TMath::Exp(-0.0372668 * coll.centFT0C()) + 20))
438+
return false;
439+
}
440+
histos.fill(HIST("hEventSelection"), 11.5 /* Remove outlyers */);
419441
}
420-
histos.fill(HIST("hEventSelection"), 11.5 /* Remove outlyers */);
421442

422443
histos.fill(HIST("hEventCentrality"), coll.centFT0C());
423444
histos.fill(HIST("hEventVertexZ"), coll.posZ());

0 commit comments

Comments
 (0)