Skip to content

Commit 2376515

Browse files
blacwovieblacwovie
authored
Update HadNucleiFemto.cxx
1 parent 2008c1f commit 2376515

File tree

1 file changed

+23
-23
lines changed

1 file changed

+23
-23
lines changed

‎PWGCF/Femto/FemtoNuclei/TableProducer/HadNucleiFemto.cxx‎

Lines changed: 23 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -92,7 +92,7 @@ enum Selections {
9292
kAll
9393
};
9494

95-
float MassHad = 0;
95+
float MassHad = 0;
9696

9797
} // namespace
9898

@@ -280,7 +280,7 @@ struct HadNucleiFemto {
280280
{"hNuPhi", "phi distribution; phi(Nu)", {HistType::kTH1F, {{600, -4.0f, 4.0f}}}},
281281
{"hHadPhi", "phi distribution; phi(had)", {HistType::kTH1F, {{600, -4.0f, 4.0f}}}},
282282
{"h2dEdxNucandidates", "dEdx distribution; #it{p} (GeV/#it{c}); dE/dx (a.u.)", {HistType::kTH2F, {{200, -5.0f, 5.0f}, {100, 0.0f, 2000.0f}}}},
283-
{"h2dEdxHadcandidates", "dEdx distribution; #it{p} (GeV/#it{c}); dE/dx (a.u.)", {HistType::kTH2F, {{200, -5.0f, 5.0f}, {100, 0.0f, 2000.0f}}}},
283+
{"h2dEdxHadcandidates", "dEdx distribution; #it{p} (GeV/#it{c}); dE/dx (a.u.)", {HistType::kTH2F, {{200, -5.0f, 5.0f}, {100, 0.0f, 2000.0f}}}},
284284
{"h2dEdx", "dEdx distribution; #it{p} (GeV/#it{c}); dE/dx (a.u.)", {HistType::kTH2F, {{200, -5.0f, 5.0f}, {100, 0.0f, 2000.0f}}}},
285285
{"h2NsigmaNuTPC", "NsigmaNu TPC distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{TPC}(Nu)", {HistType::kTH2F, {{100, -2.0f, 2.0f}, {200, -5.0f, 5.0f}}}},
286286
{"h2NsigmaNuComb", "NsigmaNu TPCTOF comb distribution; #it{p}_{T} (GeV/#it{c}); n#sigma_{comb}(Nu)", {HistType::kTH2F, {{100, -2.0f, 2.0f}, {100, 0.0f, 5.0f}}}},
@@ -507,7 +507,7 @@ struct HadNucleiFemto {
507507
if (std::abs(candidate.pt()) < settingCutHadptMin || std::abs(candidate.pt()) > settingCutHadptMax)
508508
return false;
509509

510-
// reject protons and pions
510+
// reject protons and pions
511511
if (std::abs(candidate.tpcNSigmaPr()) < settingCutNsigTPCPrMin || std::abs(candidate.tpcNSigmaPi()) < settingCutNsigTPCPiMin)
512512
return false;
513513
mQaRegistry.fill(HIST("h2NsigmaHadPrTPC"), candidate.tpcNSigmaPr());
@@ -532,14 +532,14 @@ struct HadNucleiFemto {
532532
}
533533
mQaRegistry.fill(HIST("h2NsigmaHadTPC"), candidate.sign() * candidate.pt(), tpcNSigmaKa);
534534
mQaRegistry.fill(HIST("h2NsigmaHadTOF"), candidate.sign() * candidate.pt(), tofNSigmaKa);
535-
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
535+
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
536536
return true;
537537
} else if (candidate.tpcInnerParam() < settingCutPinMinTOFHad) {
538538
if (std::abs(tpcNSigmaKa) > settingCutNsigmaTPCHad) {
539539
return false;
540540
}
541541
mQaRegistry.fill(HIST("h2NsigmaHadTPC"), candidate.sign() * candidate.pt(), tpcNSigmaKa);
542-
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
542+
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
543543
return true;
544544
}
545545
return false;
@@ -570,9 +570,9 @@ struct HadNucleiFemto {
570570
auto combNsigma = std::sqrt(tofNSigmaPi * tofNSigmaPi + tpcNSigmaPi * tpcNSigmaPi);
571571

572572
mQaRegistry.fill(HIST("h2NsigmaHadTOF_preselection"), candidate.pt(), tofNSigmaPi);
573-
//if (combNsigma > settingCutNsigmaTOFTPCHad) {
574-
// return false;
575-
//}
573+
// if (combNsigma > settingCutNsigmaTOFTPCHad) {
574+
// return false;
575+
// }
576576
if (std::abs(tofNSigmaPi) > settingCutNsigmaTOFHad) {
577577
return false;
578578
}
@@ -582,14 +582,14 @@ struct HadNucleiFemto {
582582
mQaRegistry.fill(HIST("h2NsigmaHadTPC"), candidate.sign() * candidate.pt(), tpcNSigmaPi);
583583
mQaRegistry.fill(HIST("h2NsigmaHadTOF"), candidate.sign() * candidate.pt(), tofNSigmaPi);
584584
mQaRegistry.fill(HIST("h2NsigmaHadComb"), candidate.sign() * candidate.pt(), combNsigma);
585-
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
585+
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
586586
return true;
587587
} else if (candidate.tpcInnerParam() < settingCutPinMinTOFHad) {
588588
if (std::abs(tpcNSigmaPi) > settingCutNsigmaTPCHad) {
589589
return false;
590590
}
591591
mQaRegistry.fill(HIST("h2NsigmaHadTPC"), candidate.sign() * candidate.pt(), tpcNSigmaPi);
592-
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
592+
mQaRegistry.fill(HIST("h2dEdxHadcandidates"), candidate.sign() * candidate.tpcInnerParam(), candidate.tpcSignal());
593593
return true;
594594
}
595595
return false;
@@ -598,17 +598,17 @@ struct HadNucleiFemto {
598598
template <typename Ttrack>
599599
bool selectionPIDHadron(const Ttrack& candidate)
600600
{
601-
bool PID = false;
602-
if (settingHadPDGCode == PDG_t::kPiPlus) {
603-
PID = selectionPIDPion(candidate);
604-
MassHad = o2::constants::physics::MassPiPlus;
605-
} else if (settingHadPDGCode == PDG_t::kKPlus) {
606-
PID = selectionPIDKaon(candidate);
607-
MassHad = o2::constants::physics::MassKPlus;
608-
} else {
609-
LOG(info) << "invalid PDG code";
610-
}
611-
return PID;
601+
bool PID = false;
602+
if (settingHadPDGCode == PDG_t::kPiPlus) {
603+
PID = selectionPIDPion(candidate);
604+
MassHad = o2::constants::physics::MassPiPlus;
605+
} else if (settingHadPDGCode == PDG_t::kKPlus) {
606+
PID = selectionPIDKaon(candidate);
607+
MassHad = o2::constants::physics::MassKPlus;
608+
} else {
609+
LOG(info) << "invalid PDG code";
610+
}
611+
return PID;
612612
}
613613

614614
template <typename Ttrack>
@@ -760,7 +760,7 @@ struct HadNucleiFemto {
760760
hadNucand.momNu = std::array{trackDe.px(), trackDe.py(), trackDe.pz()};
761761
hadNucand.momHad = std::array{trackHad.px(), trackHad.py(), trackHad.pz()};
762762
float invMass = 0;
763-
invMass = RecoDecay::m(std::array<std::array<float, 3>, 2>{hadNucand.momNu,hadNucand.momHad}, std::array<float, 2>{static_cast<float>(o2::constants::physics::MassDeuteron),MassHad});
763+
invMass = RecoDecay::m(std::array<std::array<float, 3>, 2>{hadNucand.momNu, hadNucand.momHad}, std::array<float, 2>{static_cast<float>(o2::constants::physics::MassDeuteron), MassHad});
764764
if (settingCutInvMass > 0 && invMass > settingCutInvMass) {
765765
mQaRegistry.fill(HIST("hSkipReasons"), 3);
766766
return false;
@@ -857,7 +857,7 @@ struct HadNucleiFemto {
857857
hadHypercand.momHad = std::array{trackHad.px(), trackHad.py(), trackHad.pz()};
858858

859859
float invMass = 0;
860-
invMass = RecoDecay::m(std::array<std::array<float, 3>, 2>{hadHypercand.momNu,hadHypercand.momHad}, std::array<float, 2>{static_cast<float>(o2::constants::physics::MassHelium3),MassHad});
860+
invMass = RecoDecay::m(std::array<std::array<float, 3>, 2>{hadHypercand.momNu, hadHypercand.momHad}, std::array<float, 2>{static_cast<float>(o2::constants::physics::MassHelium3), MassHad});
861861
if (settingCutInvMass > 0 && invMass > settingCutInvMass) {
862862
return false;
863863
}

0 commit comments

Comments
 (0)