Skip to content

Commit ce79846

Browse files
author
Preet Pati
committed
Rectification in PID track selection
1 parent 50d9c10 commit ce79846

File tree

1 file changed

+30
-29
lines changed

1 file changed

+30
-29
lines changed

‎PWGCF/TwoParticleCorrelations/Tasks/pidDiHadron.cxx‎

Lines changed: 30 additions & 29 deletions
Original file line numberDiff line numberDiff line change
@@ -126,7 +126,7 @@ struct PidDiHadron {
126126

127127
// make the filters and cuts.
128128
Filter collisionFilter = (nabs(aod::collision::posZ) < cfgCutVtxZ);
129-
Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtMin) && (aod::track::pt < cfgCutPtMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
129+
Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtMin) && (aod::track::pt < cfgCutPtMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
130130
using FilteredCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSel, aod::CentFT0Cs, aod::CentFT0CVariant1s, aod::CentFT0Ms, aod::CentFV0As, aod::Mults>>;
131131
using FilteredTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr, aod::pidTOFbeta, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>>;
132132

@@ -193,19 +193,19 @@ struct PidDiHadron {
193193

194194
// Event Counter
195195
if (doprocessSame && cfgUseAdditionalEventCut) {
196-
registry.add("hEventCountSpecific", "Number of Event;; Count", {HistType::kTH1D, {{12, 0, 12}}});
197-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(1, "after sel8");
198-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(2, "kNoSameBunchPileup");
199-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(3, "kNoITSROFrameBorder");
200-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(4, "kNoTimeFrameBorder");
201-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(5, "kIsGoodZvtxFT0vsPV");
202-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(6, "kNoCollInTimeRangeStandard");
203-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(7, "kIsGoodITSLayersAll");
204-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(8, "kNoCollInRofStandard");
205-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(9, "kNoHighMultCollInPrevRof");
206-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(10, "occupancy");
207-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(11, "MultCorrelation");
208-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(12, "cfgEvSelV0AT0ACut");
196+
registry.add("hEventCount", "Number of Event;; Count", {HistType::kTH1D, {{12, 0, 12}}});
197+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(1, "after sel8");
198+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(2, "kNoSameBunchPileup");
199+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(3, "kNoITSROFrameBorder");
200+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(4, "kNoTimeFrameBorder");
201+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(5, "kIsGoodZvtxFT0vsPV");
202+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(6, "kNoCollInTimeRangeStandard");
203+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(7, "kIsGoodITSLayersAll");
204+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(8, "kNoCollInRofStandard");
205+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(9, "kNoHighMultCollInPrevRof");
206+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(10, "occupancy");
207+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(11, "MultCorrelation");
208+
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(12, "cfgEvSelV0AT0ACut");
209209
}
210210

211211
if (cfgUseAdditionalEventCut) {
@@ -311,9 +311,6 @@ struct PidDiHadron {
311311
template <typename TTrack>
312312
bool trackSelected(TTrack track)
313313
{
314-
if (cfgPIDParticle && getNsigmaPID(track) != cfgPIDParticle) {
315-
return false;
316-
}
317314
return ((track.tpcNClsFound() >= cfgCutTPCclu) && (track.tpcNClsCrossedRows() >= cfgCutTPCCrossedRows) && (track.itsNCls() >= cfgCutITSclu));
318315
}
319316

@@ -363,6 +360,8 @@ struct PidDiHadron {
363360
for (auto const& track1 : tracks) {
364361
if (!trackSelected(track1))
365362
continue;
363+
if (cfgPIDParticle && getNsigmaPID(track1) != cfgPIDParticle)
364+
continue; // if PID is selected, check if the track has the right PID
366365
if (!getEfficiencyCorrection(weff1, track1.eta(), track1.pt(), vtxz))
367366
continue;
368367
registry.fill(HIST("Phi"), RecoDecay::constrainAngle(track1.phi(), 0.0));
@@ -426,6 +425,8 @@ struct PidDiHadron {
426425

427426
if (!trackSelected(track1))
428427
continue;
428+
if (cfgPIDParticle && getNsigmaPID(track1) != cfgPIDParticle)
429+
continue; // if PID is selected, check if the track has the right PID
429430
if (!getEfficiencyCorrection(triggerWeight, track1.eta(), track1.pt(), posZ))
430431
continue;
431432
if (system == SameEvent) {
@@ -489,61 +490,61 @@ struct PidDiHadron {
489490
template <typename TCollision>
490491
bool eventSelected(TCollision collision, const int multTrk, const float centrality, const bool fillCounter)
491492
{
492-
registry.fill(HIST("hEventCountSpecific"), 0.5);
493+
registry.fill(HIST("hEventCount"), 0.5);
493494
if (cfgEvSelkNoSameBunchPileup && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
494495
// rejects collisions which are associated with the same "found-by-T0" bunch crossing
495496
// https://indico.cern.ch/event/1396220/#1-event-selection-with-its-rof
496497
return 0;
497498
}
498499
if (fillCounter && cfgEvSelkNoSameBunchPileup)
499-
registry.fill(HIST("hEventCountSpecific"), 1.5);
500+
registry.fill(HIST("hEventCount"), 1.5);
500501
if (cfgEvSelkNoITSROFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
501502
return 0;
502503
}
503504
if (fillCounter && cfgEvSelkNoITSROFrameBorder)
504-
registry.fill(HIST("hEventCountSpecific"), 2.5);
505+
registry.fill(HIST("hEventCount"), 2.5);
505506
if (cfgEvSelkNoTimeFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
506507
return 0;
507508
}
508509
if (fillCounter && cfgEvSelkNoTimeFrameBorder)
509-
registry.fill(HIST("hEventCountSpecific"), 3.5);
510+
registry.fill(HIST("hEventCount"), 3.5);
510511
if (cfgEvSelkIsGoodZvtxFT0vsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
511512
// removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference
512513
// use this cut at low multiplicities with caution
513514
return 0;
514515
}
515516
if (fillCounter && cfgEvSelkIsGoodZvtxFT0vsPV)
516-
registry.fill(HIST("hEventCountSpecific"), 4.5);
517+
registry.fill(HIST("hEventCount"), 4.5);
517518
if (cfgEvSelkNoCollInTimeRangeStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
518519
// no collisions in specified time range
519520
return 0;
520521
}
521522
if (fillCounter && cfgEvSelkNoCollInTimeRangeStandard)
522-
registry.fill(HIST("hEventCountSpecific"), 5.5);
523+
registry.fill(HIST("hEventCount"), 5.5);
523524
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
524525
// from Jan 9 2025 AOT meeting
525526
// cut time intervals with dead ITS staves
526527
return 0;
527528
}
528529
if (fillCounter && cfgEvSelkIsGoodITSLayersAll)
529-
registry.fill(HIST("hEventCountSpecific"), 6.5);
530+
registry.fill(HIST("hEventCount"), 6.5);
530531
if (cfgEvSelkNoCollInRofStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
531532
// no other collisions in this Readout Frame with per-collision multiplicity above threshold
532533
return 0;
533534
}
534535
if (fillCounter && cfgEvSelkNoCollInRofStandard)
535-
registry.fill(HIST("hEventCountSpecific"), 7.5);
536+
registry.fill(HIST("hEventCount"), 7.5);
536537
if (cfgEvSelkNoHighMultCollInPrevRof && !collision.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof)) {
537538
// veto an event if FT0C amplitude in previous ITS ROF is above threshold
538539
return 0;
539540
}
540541
if (fillCounter && cfgEvSelkNoHighMultCollInPrevRof)
541-
registry.fill(HIST("hEventCountSpecific"), 8.5);
542+
registry.fill(HIST("hEventCount"), 8.5);
542543
auto occupancy = collision.trackOccupancyInTimeRange();
543544
if (cfgEvSelOccupancy && (occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh))
544545
return 0;
545546
if (fillCounter && cfgEvSelOccupancy)
546-
registry.fill(HIST("hEventCountSpecific"), 9.5);
547+
registry.fill(HIST("hEventCount"), 9.5);
547548

548549
auto multNTracksPV = collision.multNTracksPV();
549550
if (cfgEvSelMultCorrelation) {
@@ -557,14 +558,14 @@ struct PidDiHadron {
557558
return 0;
558559
}
559560
if (fillCounter && cfgEvSelMultCorrelation)
560-
registry.fill(HIST("hEventCountSpecific"), 10.5);
561+
registry.fill(HIST("hEventCount"), 10.5);
561562

562563
// V0A T0A 5 sigma cut
563564
float sigma = 5.0;
564565
if (cfgEvSelV0AT0ACut && (std::fabs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > sigma * fT0AV0ASigma->Eval(collision.multFT0A())))
565566
return 0;
566567
if (fillCounter && cfgEvSelV0AT0ACut)
567-
registry.fill(HIST("hEventCountSpecific"), 11.5);
568+
registry.fill(HIST("hEventCount"), 11.5);
568569

569570
return 1;
570571
}

0 commit comments

Comments
 (0)