Skip to content

Commit 543f551

Browse files
authored
[PWGHF] added Lambda0 processes to calculate data driven efficiency
1 parent f40ea3e commit 543f551

File tree

1 file changed

+34
-34
lines changed

1 file changed

+34
-34
lines changed

PWGHF/HFC/TableProducer/correlatorLcScHadrons.cxx

Lines changed: 34 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -472,7 +472,7 @@ struct HfCorrelatorLcScHadrons {
472472
}
473473

474474
template <typename TCollision, typename V0>
475-
bool SelectionV0(TCollision const& collision, V0 const& candidate)
475+
bool selectionV0(TCollision const& collision, V0 const& candidate)
476476
{
477477
if (candidate.v0radius() < cfgV0.cfgV0radiusMin){
478478
return false;
@@ -534,65 +534,65 @@ struct HfCorrelatorLcScHadrons {
534534
}
535535

536536
template <bool isMcRec = false, typename CollType, typename V0, typename TrackType>
537-
void FillV0Histograms(CollType const& collV0, V0 const& v0s, TrackType const&)
537+
void fillV0Histograms(CollType const& collV0, V0 const& v0s, TrackType const&)
538538
{
539-
for (auto& v0 : v0s) {
539+
for (const auto& v0 : v0s) {
540540

541-
if (!SelectionV0(collV0, v0)){
541+
if (!selectionV0(collV0, v0)){
542542
continue;
543543
}
544544

545-
auto posTrack_V0 = v0.template posTrack_as<TrackType>();
546-
auto negTrack_V0 = v0.template negTrack_as<TrackType>();
545+
auto posTrackV0 = v0.template posTrack_as<TrackType>();
546+
auto negTrackV0 = v0.template negTrack_as<TrackType>();
547547

548-
if (isSelectedV0Daughter(posTrack_V0, kProton) && isSelectedV0Daughter(negTrack_V0, kPiPlus)) {
548+
if (isSelectedV0Daughter(posTrackV0, kProton) && isSelectedV0Daughter(negTrackV0, kPiPlus)) {
549549
if (std::abs(massLambda - v0.mLambda()) < cfgV0.cfgHypMassWindow){
550-
registry.fill(HIST("hV0Lambda"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
551-
registry.fill(HIST("hV0LambdaRefl"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
550+
registry.fill(HIST("hV0Lambda"), v0.mLambda(), v0.pt(), posTrackV0.pt());
551+
registry.fill(HIST("hV0LambdaRefl"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
552552

553-
if(passPIDSelection(posTrack_V0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
554-
registry.fill(HIST("hV0LambdaPiKRej"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
555-
registry.fill(HIST("hV0LambdaReflPiKRej"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
553+
if(passPIDSelection(posTrackV0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
554+
registry.fill(HIST("hV0LambdaPiKRej"), v0.mLambda(), v0.pt(), posTrackV0.pt());
555+
registry.fill(HIST("hV0LambdaReflPiKRej"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
556556
}
557557
}
558558

559559
}
560-
if (isSelectedV0Daughter(negTrack_V0, kProton) && isSelectedV0Daughter(posTrack_V0, kPiPlus)) {
560+
if (isSelectedV0Daughter(negTrackV0, kProton) && isSelectedV0Daughter(posTrackV0, kPiPlus)) {
561561
if (std::abs(massLambda - v0.mAntiLambda()) > cfgV0.cfgHypMassWindow){
562-
registry.fill(HIST("hV0Lambda"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
563-
registry.fill(HIST("hV0LambdaRefl"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
562+
registry.fill(HIST("hV0Lambda"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
563+
registry.fill(HIST("hV0LambdaRefl"), v0.mLambda(), v0.pt(), posTrackV0.pt());
564564

565-
if(passPIDSelection(negTrack_V0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
566-
registry.fill(HIST("hV0LambdaPiKRej"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
567-
registry.fill(HIST("hV0LambdaReflPiKRej"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
565+
if(passPIDSelection(negTrackV0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
566+
registry.fill(HIST("hV0LambdaPiKRej"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
567+
registry.fill(HIST("hV0LambdaReflPiKRej"), v0.mLambda(), v0.pt(), posTrackV0.pt());
568568
}
569569
}
570570
}
571571
if constexpr (isMcRec){
572-
if (!v0.has_mcParticle() || !posTrack_V0.has_mcParticle() || !negTrack_V0.has_mcParticle()){
572+
if (!v0.has_mcParticle() || !posTrackV0.has_mcParticle() || !negTrackV0.has_mcParticle()){
573573
continue;
574574
}
575575
auto v0Mc = v0.mcParticle();
576-
auto posTrack = posTrack_V0.mcParticle();
577-
auto negTrack = negTrack_V0.mcParticle();
576+
auto posTrack = posTrackV0.mcParticle();
577+
auto negTrack = negTrackV0.mcParticle();
578578

579579
if(std::abs(v0Mc.pdgCode()) == kLambda0) {
580580
if (std::abs(posTrack.pdgCode()) == kProton){
581-
registry.fill(HIST("hV0LambdaMcRec"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
582-
registry.fill(HIST("hV0LambdaReflMcRec"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
581+
registry.fill(HIST("hV0LambdaMcRec"), v0.mLambda(), v0.pt(), posTrackV0.pt());
582+
registry.fill(HIST("hV0LambdaReflMcRec"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
583583

584-
if(passPIDSelection(posTrack_V0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
585-
registry.fill(HIST("hV0LambdaPiKRejMcRec"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
586-
registry.fill(HIST("hV0LambdaReflPiKRejMcRec"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
584+
if(passPIDSelection(posTrackV0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
585+
registry.fill(HIST("hV0LambdaPiKRejMcRec"), v0.mLambda(), v0.pt(), posTrackV0.pt());
586+
registry.fill(HIST("hV0LambdaReflPiKRejMcRec"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
587587
}
588588
}
589589
if (std::abs(negTrack.pdgCode()) == kProton){
590-
registry.fill(HIST("hV0LambdaMcRec"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
591-
registry.fill(HIST("hV0LambdaReflMcRec"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
592-
593-
if(passPIDSelection(negTrack_V0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
594-
registry.fill(HIST("hV0LambdaPiKRejMcRec"), v0.mAntiLambda(), v0.pt(), negTrack_V0.pt());
595-
registry.fill(HIST("hV0LambdaReflPiKRejMcRec"), v0.mLambda(), v0.pt(), posTrack_V0.pt());
590+
registry.fill(HIST("hV0LambdaMcRec"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
591+
registry.fill(HIST("hV0LambdaReflMcRec"), v0.mLambda(), v0.pt(), posTrackV0.pt());
592+
593+
if(passPIDSelection(negTrackV0, trkPIDspecies, pidTPCMax, pidTOFMax, tofPIDThreshold, forceTOF)){
594+
registry.fill(HIST("hV0LambdaPiKRejMcRec"), v0.mAntiLambda(), v0.pt(), negTrackV0.pt());
595+
registry.fill(HIST("hV0LambdaReflPiKRejMcRec"), v0.mLambda(), v0.pt(), posTrackV0.pt());
596596
}
597597
}
598598
}
@@ -1394,7 +1394,7 @@ struct HfCorrelatorLcScHadrons {
13941394
}
13951395
registry.fill(HIST("hEventLambdaV0"), 1.5);
13961396

1397-
FillV0Histograms<false>(collision, V0s, tracks);
1397+
fillV0Histograms<false>(collision, V0s, tracks);
13981398
}
13991399
PROCESS_SWITCH(HfCorrelatorLcScHadrons, processDataLambdaV0, "Data process for v0 lambda", false);
14001400

@@ -1407,7 +1407,7 @@ struct HfCorrelatorLcScHadrons {
14071407
}
14081408
registry.fill(HIST("hEventLambdaV0"), 1.5);
14091409

1410-
FillV0Histograms<true>(collision, V0s, tracks);
1410+
fillV0Histograms<true>(collision, V0s, tracks);
14111411
}
14121412
PROCESS_SWITCH(HfCorrelatorLcScHadrons, processMcLambdaV0, "Mc process for v0 lambda", false);
14131413

0 commit comments

Comments
 (0)