Skip to content

Commit 1162a09

Browse files
committed
candidateCreator2Prong.cxx: Simplify masses
1 parent bd6a438 commit 1162a09

File tree

1 file changed

+9
-21
lines changed

1 file changed

+9
-21
lines changed

PWGHF/TableProducer/candidateCreator2Prong.cxx

Lines changed: 9 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -115,22 +115,15 @@ struct HfCandidateCreator2Prong {
115115
o2::vertexing::DCAFitterN<2> df; // 2-prong vertex fitter
116116
Service<o2::ccdb::BasicCCDBManager> ccdb;
117117

118-
using TracksWCovExtraPidPiKa = soa::Join<aod::TracksWCovExtra, aod::TracksPidPi, aod::PidTpcTofFullPi, aod::TracksPidKa, aod::PidTpcTofFullKa>;
119-
120118
int runNumber{0};
121-
float toMicrometers = 10000.; // from cm to µm
122-
double massPi{0.};
123-
double massK{0.};
124-
double massE{0.};
125-
double massMu{0.};
126-
double massPiK{0.};
127-
double massKPi{0.};
128-
double massEE{0.};
129-
double massMuMu{0.};
130119
double bz{0.};
131120

121+
const float toMicrometers = 10000.; // from cm to µm
122+
132123
std::shared_ptr<TH1> hCandidates;
133124

125+
using TracksWCovExtraPidPiKa = soa::Join<aod::TracksWCovExtra, aod::TracksPidPi, aod::PidTpcTofFullPi, aod::TracksPidKa, aod::PidTpcTofFullKa>;
126+
134127
ConfigurableAxis axisMass{"axisMass", {500, 1.6, 2.1}, "axis for mass (GeV/c^2)"};
135128

136129
HistogramRegistry registry{"registry"};
@@ -192,11 +185,6 @@ struct HfCandidateCreator2Prong {
192185
// init HF event selection helper
193186
hfEvSel.init(registry);
194187

195-
massPi = MassPiPlus;
196-
massK = MassKPlus;
197-
massE = MassElectron;
198-
massMu = MassMuon;
199-
200188
if (std::accumulate(doprocessDF.begin(), doprocessDF.end(), 0) == 1) {
201189
registry.fill(HIST("hVertexerType"), aod::hf_cand::VertexerType::DCAFitter);
202190
// Configure DCAFitterN
@@ -366,11 +354,11 @@ struct HfCandidateCreator2Prong {
366354
// fill histograms
367355
if (fillHistograms) {
368356
// calculate invariant masses
369-
auto arrayMomenta = std::array{pvec0, pvec1};
370-
massPiK = RecoDecay::m(arrayMomenta, std::array{massPi, massK});
371-
massKPi = RecoDecay::m(arrayMomenta, std::array{massK, massPi});
372-
massEE = RecoDecay::m(arrayMomenta, std::array{massE, massE});
373-
massMuMu = RecoDecay::m(arrayMomenta, std::array{massMu, massMu});
357+
const auto arrayMomenta = std::array{pvec0, pvec1};
358+
const auto massPiK = RecoDecay::m(arrayMomenta, std::array{MassPiPlus, MassKPlus});
359+
const auto massKPi = RecoDecay::m(arrayMomenta, std::array{MassKPlus, MassPiPlus});
360+
const auto massEE = RecoDecay::m(arrayMomenta, std::array{MassElectron, MassElectron});
361+
const auto massMuMu = RecoDecay::m(arrayMomenta, std::array{MassMuon, MassMuon});
374362
registry.fill(HIST("hMass2"), massPiK);
375363
registry.fill(HIST("hMass2"), massKPi);
376364
registry.fill(HIST("hMassEE"), massEE);

0 commit comments

Comments
 (0)