Skip to content

Commit 3d95ef4

Browse files
sawankumawatSawan Sawan
sawankumawat
and
Sawan Sawan
authored
[PWGLF] Minor corrections (#12697)
Co-authored-by: Sawan Sawan <sawan.sawan@cern.ch>
1 parent 136e75c commit 3d95ef4

File tree

1 file changed

+51
-47
lines changed

1 file changed

+51
-47
lines changed

‎PWGLF/Tasks/Resonances/higherMassResonances.cxx‎

Lines changed: 51 additions & 47 deletions
Original file line numberDiff line numberDiff line change
@@ -237,8 +237,8 @@ struct HigherMassResonances {
237237
if (config.qAevents) {
238238
rEventSelection.add("hVertexZRec", "hVertexZRec", {HistType::kTH1F, {vertexZAxis}});
239239
rEventSelection.add("hmultiplicity", "multiplicity percentile distribution", {HistType::kTH1F, {{150, 0.0f, 150.0f}}});
240-
hglue.add("htrackscheck_v0", "htrackscheck_v0", kTH1I, {{15, 0, 15}});
241-
hglue.add("htrackscheck_v0_daughters", "htrackscheck_v0_daughters", kTH1I, {{15, 0, 15}});
240+
rEventSelection.add("htrackscheck_v0", "htrackscheck_v0", kTH1I, {{15, 0, 15}});
241+
rEventSelection.add("htrackscheck_v0_daughters", "htrackscheck_v0_daughters", kTH1I, {{15, 0, 15}});
242242
hMChists.add("events_check", "No. of events in the generated MC", kTH1I, {{20, 0, 20}});
243243
hMChists.add("events_checkrec", "No. of events in the reconstructed MC", kTH1I, {{25, 0, 25}});
244244

@@ -282,7 +282,7 @@ struct HigherMassResonances {
282282
hv0DauLabel->GetXaxis()->SetBinLabel(7, "Eta");
283283
hv0DauLabel->GetXaxis()->SetBinLabel(8, "PID TPC");
284284

285-
std::shared_ptr<TH1> hv0labelmcrec = rEventSelection.get<TH1>(HIST("events_checkrec"));
285+
std::shared_ptr<TH1> hv0labelmcrec = hMChists.get<TH1>(HIST("events_checkrec"));
286286
hv0labelmcrec->GetXaxis()->SetBinLabel(1, "All Tracks");
287287
hv0labelmcrec->GetXaxis()->SetBinLabel(2, "V0Daughter Sel.");
288288
hv0labelmcrec->GetXaxis()->SetBinLabel(3, "V0 Sel.");
@@ -455,52 +455,52 @@ struct HigherMassResonances {
455455
if (config.correlation2Dhist)
456456
rKzeroShort.fill(HIST("mass_lambda_kshort_before"), candidate.mK0Short(), candidate.mLambda());
457457

458-
hglue.fill(HIST("htrackscheck_v0"), 0.5);
458+
rEventSelection.fill(HIST("htrackscheck_v0"), 0.5);
459459

460460
if (config.isApplyDCAv0topv && std::fabs(candidate.dcav0topv()) > config.cMaxV0DCA) {
461461
return false;
462462
}
463-
hglue.fill(HIST("htrackscheck_v0"), 1.5);
463+
rEventSelection.fill(HIST("htrackscheck_v0"), 1.5);
464464

465465
if (std::abs(candidate.rapidity(0)) >= config.confKsrapidity) {
466466
return false;
467467
}
468-
hglue.fill(HIST("htrackscheck_v0"), 2.5);
468+
rEventSelection.fill(HIST("htrackscheck_v0"), 2.5);
469469

470470
if (pT < config.confV0PtMin) {
471471
return false;
472472
}
473-
hglue.fill(HIST("htrackscheck_v0"), 3.5);
473+
rEventSelection.fill(HIST("htrackscheck_v0"), 3.5);
474474

475475
if (dcaDaughv0 > config.confV0DCADaughMax) {
476476
return false;
477477
}
478-
hglue.fill(HIST("htrackscheck_v0"), 4.5);
478+
rEventSelection.fill(HIST("htrackscheck_v0"), 4.5);
479479

480480
if (cpav0 < config.confV0CPAMin) {
481481
return false;
482482
}
483-
hglue.fill(HIST("htrackscheck_v0"), 5.5);
483+
rEventSelection.fill(HIST("htrackscheck_v0"), 5.5);
484484

485485
if (tranRad < config.confV0TranRadV0Min) {
486486
return false;
487487
}
488-
hglue.fill(HIST("htrackscheck_v0"), 6.5);
488+
rEventSelection.fill(HIST("htrackscheck_v0"), 6.5);
489489

490490
if (tranRad > config.confV0TranRadV0Max) {
491491
return false;
492492
}
493-
hglue.fill(HIST("htrackscheck_v0"), 7.5);
493+
rEventSelection.fill(HIST("htrackscheck_v0"), 7.5);
494494

495495
if (std::fabs(ctauK0s) > config.cMaxV0LifeTime) {
496496
return false;
497497
}
498-
hglue.fill(HIST("htrackscheck_v0"), 8.5);
498+
rEventSelection.fill(HIST("htrackscheck_v0"), 8.5);
499499

500500
if (config.isapplyCompetingcut && (std::abs(candidate.mLambda() - o2::constants::physics::MassLambda0) <= config.competingcascrejlambda || std::abs(candidate.mAntiLambda() - o2::constants::physics::MassLambda0) <= config.competingcascrejlambda)) {
501501
return false;
502502
}
503-
hglue.fill(HIST("htrackscheck_v0"), 9.5);
503+
rEventSelection.fill(HIST("htrackscheck_v0"), 9.5);
504504

505505
if (config.correlation2Dhist)
506506
rKzeroShort.fill(HIST("mass_lambda_kshort_after10"), candidate.mK0Short(), candidate.mLambda());
@@ -512,12 +512,12 @@ struct HigherMassResonances {
512512
if (config.isStandardV0 && candidate.v0Type() != 1) {
513513
return false; // Only standard V0s are selected
514514
}
515-
hglue.fill(HIST("htrackscheck_v0"), 10.5);
515+
rEventSelection.fill(HIST("htrackscheck_v0"), 10.5);
516516

517517
if (candidate.mK0Short() < lowmasscutks0 || candidate.mK0Short() > highmasscutks0) {
518518
return false;
519519
}
520-
hglue.fill(HIST("htrackscheck_v0"), 11.5);
520+
rEventSelection.fill(HIST("htrackscheck_v0"), 11.5);
521521

522522
return true;
523523
}
@@ -533,44 +533,44 @@ struct HigherMassResonances {
533533
const auto tpcNClsF = track.tpcNClsFound();
534534
const auto sign = track.sign();
535535

536-
hglue.fill(HIST("htrackscheck_v0_daughters"), 0.5);
536+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 0.5);
537537

538538
if (config.hasTPC && !track.hasTPC())
539539
return false;
540-
hglue.fill(HIST("htrackscheck_v0_daughters"), 1.5);
540+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 1.5);
541541

542542
if (track.tpcNClsCrossedRows() < config.tpcCrossedrows)
543543
return false;
544-
hglue.fill(HIST("htrackscheck_v0_daughters"), 2.5);
544+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 2.5);
545545

546546
if (track.tpcCrossedRowsOverFindableCls() < config.tpcCrossedrowsOverfcls)
547547
return false;
548-
hglue.fill(HIST("htrackscheck_v0_daughters"), 3.5);
548+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 3.5);
549549

550550
if (tpcNClsF < config.confDaughTPCnclsMin) {
551551
return false;
552552
}
553-
hglue.fill(HIST("htrackscheck_v0_daughters"), 4.5);
553+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 4.5);
554554

555555
if (charge < 0 && sign > 0) {
556556
return false;
557557
}
558-
hglue.fill(HIST("htrackscheck_v0_daughters"), 5.5);
558+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 5.5);
559559

560560
if (charge > 0 && sign < 0) {
561561
return false;
562562
}
563-
hglue.fill(HIST("htrackscheck_v0_daughters"), 6.5);
563+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 6.5);
564564

565565
if (std::abs(eta) > config.confDaughEta) {
566566
return false;
567567
}
568-
hglue.fill(HIST("htrackscheck_v0_daughters"), 7.5);
568+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 7.5);
569569

570570
if (std::abs(nsigmaV0DaughterTPC) > config.confDaughPIDCutTPC) {
571571
return false;
572572
}
573-
hglue.fill(HIST("htrackscheck_v0_daughters"), 8.5);
573+
rEventSelection.fill(HIST("htrackscheck_v0_daughters"), 8.5);
574574

575575
// if (std::abs())
576576

@@ -887,8 +887,12 @@ struct HigherMassResonances {
887887
}
888888
PROCESS_SWITCH(HigherMassResonances, processSE, "same event process", true);
889889

890-
using EventCandidatesDerivedData = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraStamps, aod::StraZDCSP>;
891-
using V0CandidatesDerivedData = soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras, aod::V0TOFPIDs, aod::V0TOFNSigmas>;
890+
// using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::MultZeqs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::CentFV0As, aod::Mults, aod::PVMults>>;
891+
// using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTOFFullPi>>;
892+
// using V0TrackCandidate = aod::V0Datas;
893+
894+
using EventCandidatesDerivedData = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraStamps>;
895+
using V0CandidatesDerivedData = soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras>;
892896
using DauTracks = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;
893897

894898
ConfigurableAxis axisVertex{"axisVertex", {20, -10, 10}, "vertex axis for ME mixing"};
@@ -1280,7 +1284,7 @@ struct HigherMassResonances {
12801284
}
12811285
PROCESS_SWITCH(HigherMassResonances, processRec, "Process Reconstructed", false);
12821286

1283-
void processSEderived(EventCandidatesDerivedData::iterator const& collision, TrackCandidates const& /*tracks*/, aod::V0Datas const& V0s)
1287+
void processSEderived(EventCandidatesDerivedData::iterator const& collision, V0CandidatesDerivedData const& V0s)
12841288
{
12851289
if (config.cSelectMultEstimator == kFT0M) {
12861290
multiplicity = collision.centFT0M();
@@ -1319,22 +1323,22 @@ struct HigherMassResonances {
13191323
continue;
13201324
}
13211325

1322-
auto postrack1 = v1.template posTrack_as<TrackCandidates>();
1323-
auto negtrack1 = v1.template negTrack_as<TrackCandidates>();
1324-
auto postrack2 = v2.template posTrack_as<TrackCandidates>();
1325-
auto negtrack2 = v2.template negTrack_as<TrackCandidates>();
1326+
// auto postrack1 = v1.template posTrack_as<DauTracks>();
1327+
// auto negtrack1 = v1.template negTrack_as<DauTracks>();
1328+
// auto postrack2 = v2.template posTrack_as<DauTracks>();
1329+
// auto negtrack2 = v2.template negTrack_as<DauTracks>();
13261330

1327-
double nTPCSigmaPos1{postrack1.tpcNSigmaPi()};
1328-
double nTPCSigmaNeg1{negtrack1.tpcNSigmaPi()};
1329-
double nTPCSigmaPos2{postrack2.tpcNSigmaPi()};
1330-
double nTPCSigmaNeg2{negtrack2.tpcNSigmaPi()};
1331+
// double nTPCSigmaPos1{postrack1.tpcNSigmaPi()};
1332+
// double nTPCSigmaNeg1{negtrack1.tpcNSigmaPi()};
1333+
// double nTPCSigmaPos2{postrack2.tpcNSigmaPi()};
1334+
// double nTPCSigmaNeg2{negtrack2.tpcNSigmaPi()};
13311335

1332-
if (!(isSelectedV0Daughter(negtrack1, -1, nTPCSigmaNeg1, v1) && isSelectedV0Daughter(postrack1, 1, nTPCSigmaPos1, v1))) {
1333-
continue;
1334-
}
1335-
if (!(isSelectedV0Daughter(postrack2, 1, nTPCSigmaPos2, v2) && isSelectedV0Daughter(negtrack2, -1, nTPCSigmaNeg2, v2))) {
1336-
continue;
1337-
}
1336+
// if (!(isSelectedV0Daughter(negtrack1, -1, nTPCSigmaNeg1, v1) && isSelectedV0Daughter(postrack1, 1, nTPCSigmaPos1, v1))) {
1337+
// continue;
1338+
// }
1339+
// if (!(isSelectedV0Daughter(postrack2, 1, nTPCSigmaPos2, v2) && isSelectedV0Daughter(negtrack2, -1, nTPCSigmaNeg2, v2))) {
1340+
// continue;
1341+
// }
13381342

13391343
if (std::find(v0indexes.begin(), v0indexes.end(), v1.globalIndex()) == v0indexes.end()) {
13401344
v0indexes.push_back(v1.globalIndex());
@@ -1356,12 +1360,12 @@ struct HigherMassResonances {
13561360
// rKzeroShort.fill(HIST("positive_phi"), postrack1.phi());
13571361
// }
13581362

1359-
if (postrack1.globalIndex() == postrack2.globalIndex()) {
1360-
continue;
1361-
}
1362-
if (negtrack1.globalIndex() == negtrack2.globalIndex()) {
1363-
continue;
1364-
}
1363+
// if (postrack1.globalIndex() == postrack2.globalIndex()) {
1364+
// continue;
1365+
// }
1366+
// if (negtrack1.globalIndex() == negtrack2.globalIndex()) {
1367+
// continue;
1368+
// }
13651369

13661370
if (!applyAngSep(v1, v2)) {
13671371
continue;

0 commit comments

Comments
 (0)