@@ -690,54 +690,54 @@ struct Pi0EtaToGammaGammaMC {
690690 }
691691 if constexpr (requires { g2.omegaMBWeight (); }) {
692692 wpair *= g2.omegaMBWeight ();
693- float alphaMeson = std::fabs (g1.e () - g2.e ()) / (g1.e () + g2.e ());
694- float alphaCut = 999 .f ;
695- switch (static_cast <AlphaMesonCutOption>(cfgAlphaMesonCut.value )) {
696- case AlphaMesonCutOption::Off:
697- break ;
698- case AlphaMesonCutOption::SpecificValue:
699- alphaCut = cfgAlphaMeson;
700- break ;
701- case AlphaMesonCutOption::PTDependent: {
702- alphaCut = cfgAlphaMesonA * std::tanh (cfgAlphaMesonB * v12.pt ());
703- break ;
693+ float alphaMeson = std::fabs (g1.e () - g2.e ()) / (g1.e () + g2.e ());
694+ float alphaCut = 999 .f ;
695+ switch (static_cast <AlphaMesonCutOption>(cfgAlphaMesonCut.value )) {
696+ case AlphaMesonCutOption::Off:
697+ break ;
698+ case AlphaMesonCutOption::SpecificValue:
699+ alphaCut = cfgAlphaMeson;
700+ break ;
701+ case AlphaMesonCutOption::PTDependent: {
702+ alphaCut = cfgAlphaMesonA * std::tanh (cfgAlphaMesonB * v12.pt ());
703+ break ;
704+ }
705+ default :
706+ LOGF (error, " Invalid option for alpha meson cut. No alpha cut will be applied." );
704707 }
705- default :
706- LOGF (error, " Invalid option for alpha meson cut. No alpha cut will be applied." );
707- }
708- if (alphaMeson > alphaCut) {
709- continue ;
710- }
711-
712- if (pairtype == o2::aod::pwgem::photonmeson::photonpair::PairType::kEMCEMC ) {
713- float openingAngle = std::acos (v1.Vect ().Dot (v2.Vect ()) / (v1.P () * v2.P ()));
714- if (openingAngle < emccuts.minOpenAngle ) {
708+ if (alphaMeson > alphaCut) {
715709 continue ;
716710 }
717- }
718711
719- if (g1mc.globalIndex () == g2mc.globalIndex ()) {
720- if (o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (g1mc, mcparticles) == 111 )
721- fRegistry .fill (HIST (" Pair/Pi0/hs_FromSameGamma" ), v12.M (), v12.Pt (), wpair);
722- else if (o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (g1mc, mcparticles) == 221 )
723- fRegistry .fill (HIST (" Pair/Eta/hs_FromSameGamma" ), v12.M (), v12.Pt (), wpair);
724- continue ;
725- }
712+ if (pairtype == o2::aod::pwgem::photonmeson::photonpair::PairType::kEMCEMC ) {
713+ float openingAngle = std::acos (v1.Vect ().Dot (v2.Vect ()) / (v1.P () * v2.P ()));
714+ if (openingAngle < emccuts.minOpenAngle ) {
715+ continue ;
716+ }
717+ }
726718
727- if (pi0id > 0 ) {
728- auto pi0mc = mcparticles.iteratorAt (pi0id);
729- if (cfgRequireTrueAssociation && (pi0mc.emmceventId () != collision.emmceventId ())) {
719+ if (g1mc.globalIndex () == g2mc.globalIndex ()) {
720+ if (o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (g1mc, mcparticles) == 111 )
721+ fRegistry .fill (HIST (" Pair/Pi0/hs_FromSameGamma" ), v12.M (), v12.Pt (), wpair);
722+ else if (o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (g1mc, mcparticles) == 221 )
723+ fRegistry .fill (HIST (" Pair/Eta/hs_FromSameGamma" ), v12.M (), v12.Pt (), wpair);
730724 continue ;
731725 }
732- o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo (&fRegistry , v12, pi0mc, mcparticles, mccollisions, f1fd_k0s_to_pi0, wpair);
733- } else if (etaid > 0 ) {
734- auto etamc = mcparticles.iteratorAt (etaid);
735- if (cfgRequireTrueAssociation && (etamc.emmceventId () != collision.emmceventId ())) {
736- continue ;
726+
727+ if (pi0id > 0 ) {
728+ auto pi0mc = mcparticles.iteratorAt (pi0id);
729+ if (cfgRequireTrueAssociation && (pi0mc.emmceventId () != collision.emmceventId ())) {
730+ continue ;
731+ }
732+ o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo (&fRegistry , v12, pi0mc, mcparticles, mccollisions, f1fd_k0s_to_pi0, wpair);
733+ } else if (etaid > 0 ) {
734+ auto etamc = mcparticles.iteratorAt (etaid);
735+ if (cfgRequireTrueAssociation && (etamc.emmceventId () != collision.emmceventId ())) {
736+ continue ;
737+ }
738+ o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo (&fRegistry , v12, etamc, mcparticles, mccollisions, f1fd_k0s_to_pi0, wpair);
737739 }
738- o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo (&fRegistry , v12, etamc, mcparticles, mccollisions, f1fd_k0s_to_pi0, wpair);
739- }
740- } // end of pairing loop
740+ } // end of pairing loop
741741 } else if constexpr (pairtype == o2::aod::pwgem::photonmeson::photonpair::PairType::kPCMDalitzEE ) {
742742
743743 auto photons1_per_collision = photons1.sliceByCached (TDetectorTag1::perCollision (), collision.globalIndex (), cache);
0 commit comments