diff --git a/PWGEM/PhotonMeson/TableProducer/associateMCinfoPhoton.cxx b/PWGEM/PhotonMeson/TableProducer/associateMCinfoPhoton.cxx index daa5e7011b8..f2f2a3eed2d 100644 --- a/PWGEM/PhotonMeson/TableProducer/associateMCinfoPhoton.cxx +++ b/PWGEM/PhotonMeson/TableProducer/associateMCinfoPhoton.cxx @@ -66,7 +66,7 @@ struct AssociateMCInfoPhoton { Configurable max_rxy_gen{"max_rxy_gen", 100, "max rxy to store generated information"}; Configurable requireGammaGammaDecay{"requireGammaGammaDecay", false, "require gamma gamma decay for generated pi0 and eta meson"}; Configurable cfg_max_eta_photon{"cfg_max_eta_gamma", 0.8, "max eta for photons at PV"}; - + HistogramRegistry registry{"EMMCEvent"}; void init(o2::framework::InitContext&) @@ -172,7 +172,7 @@ struct AssociateMCInfoPhoton { for (const auto& mcParticle : groupedMcParticles) { // store necessary information for denominator of efficiency if ((mcParticle.isPhysicalPrimary() || mcParticle.producedByGenerator()) && std::fabs(mcParticle.y()) < 0.9f && mcParticle.pt() < 20.f) { auto binNumber = hBinFinder->FindBin(mcParticle.pt(), std::fabs(mcParticle.y())); // caution: pack - size_t lNDaughters = 0; + size_t lNDaughters = 0; size_t mesonAccepted = 0; if (mcParticle.has_daughters()) { auto lDaughters = mcParticle.daughters_as(); @@ -180,11 +180,11 @@ struct AssociateMCInfoPhoton { auto lDaughter0 = lDaughters.begin(); if (lNDaughters == 2) { auto lDaughter1 = lDaughters.iteratorAt(1); - if ( std::fabs(lDaughter0.eta())< cfg_max_eta_photon && std::fabs(lDaughter1.eta()) < cfg_max_eta_photon){ + if (std::fabs(lDaughter0.eta()) < cfg_max_eta_photon && std::fabs(lDaughter1.eta()) < cfg_max_eta_photon) { mesonAccepted = 1; } - } - } + } + } switch (std::abs(mcParticle.pdgCode())) { case PDG_t::kGamma: registry.fill(HIST("Generated/h2PtY_Gamma"), mcParticle.pt(), std::fabs(mcParticle.y())); @@ -195,7 +195,7 @@ struct AssociateMCInfoPhoton { continue; registry.fill(HIST("Generated/h2PtY_Pi0"), mcParticle.pt(), std::fabs(mcParticle.y())); genPi0[binNumber]++; - if (mesonAccepted == 1){ + if (mesonAccepted == 1) { registry.fill(HIST("Generated/h2PtY_Accepted_Pi0"), mcParticle.pt(), std::fabs(mcParticle.y())); } break; @@ -204,7 +204,7 @@ struct AssociateMCInfoPhoton { continue; registry.fill(HIST("Generated/h2PtY_Eta"), mcParticle.pt(), std::fabs(mcParticle.y())); genEta[binNumber]++; - if (mesonAccepted == 1){ + if (mesonAccepted == 1) { registry.fill(HIST("Generated/h2PtY_Accepted_Eta"), mcParticle.pt(), std::fabs(mcParticle.y())); } break;