@@ -251,9 +251,9 @@ struct Pi0EtaToGammaGammaMC {
251251 DefinePHOSCut ();
252252
253253 f1fd_k0s_to_pi0 = new TF1 (" f1fd_k0s_to_pi0" , TString (fd_k0s_to_pi0), 0 .f , 100 .f );
254-
254+
255255 fRegistry .add (" Event/hNrecPerMCCollision" , " Nrec per mc collision;N_{rec} collisions per MC collision" , o2::framework::HistType::kTH1F , {{21 , -0 .5f , 20 .5f }}, false );
256- if (cfgGGContaCheck){
256+ if (cfgGGContaCheck) {
257257 fRegistry .add (" Event/hNGGContamEta" , " Number of Eta from etaToGammaGamma; p_{T, #eta} (GeV/#it{c}); N" , o2::framework::HistType::kTH1F , {{40 , -0 .5f , 20 .5f }}, false );
258258 fRegistry .add (" Event/hNGGContamPion" , " Number of Pion from etaToGammaGamma; p_{T, #pi} (GeV/#it{c}); N" , o2::framework::HistType::kTH1F , {{40 , -0 .5f , 20 .5f }}, false );
259259 }
@@ -799,19 +799,19 @@ struct Pi0EtaToGammaGammaMC {
799799
800800 auto pos2mc = mcparticles.iteratorAt (pos2.emmcparticleId ());
801801 auto ele2mc = mcparticles.iteratorAt (ele2.emmcparticleId ());
802- if (cfgGGContaCheck){
803- photonid2 = o2::aod::pwgem::dilepton::utils::mcutil::FindCommonMotherFrom2Prongs (pos2mc, ele2mc, -11 , 11 , 22 , mcparticles); // check possible contamination
804- if (photonid2 > 0 ){
802+ if (cfgGGContaCheck) {
803+ photonid2 = o2::aod::pwgem::dilepton::utils::mcutil::FindCommonMotherFrom2Prongs (pos2mc, ele2mc, -11 , 11 , 22 , mcparticles); // check possible contamination
804+ if (photonid2 > 0 ) {
805805 auto photon2 = mcparticles.iteratorAt (photonid2);
806806 int photon2pdg = photon2.pdgCode ();
807807 int photon2mothid = photon2.mothersIds ()[0 ];
808808 auto photon2moth = mcparticles.iteratorAt (photon2mothid);
809- if (photon2pdg == 22 && (o2::aod::pwgem::photonmeson::utils::mcutil::isGammaGammaDecay (photon2moth, mcparticles))){
810- int mothID = o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (photon2, mcparticles);
811- if (mothID == 221 ) {
809+ if (photon2pdg == 22 && (o2::aod::pwgem::photonmeson::utils::mcutil::isGammaGammaDecay (photon2moth, mcparticles))) {
810+ int mothID = o2::aod::pwgem::dilepton::utils::mcutil::getMotherPDGCode (photon2, mcparticles);
811+ if (mothID == 221 ) {
812812 fRegistry .fill (HIST (" Event/hNGGContamEta" ), photon2moth.pt ());
813- }
814- if (mothID == 111 ) {
813+ }
814+ if (mothID == 111 ) {
815815 fRegistry .fill (HIST (" Event/hNGGContamPion" ), photon2moth.pt ());
816816 }
817817 }
@@ -844,8 +844,8 @@ struct Pi0EtaToGammaGammaMC {
844844 }
845845 o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo (&fRegistry , veeg, etamc, mcparticles, mccollisions, f1fd_k0s_to_pi0, weight);
846846 }
847- } // end of dielectron loop
848- } // end of pcm loop
847+ } // end of dielectron loop
848+ } // end of pcm loop
849849 } else { // PCM-EMC, PCM-PHOS.
850850 // TODO: implement proper functionality if we ever want to run this in Pb-Pb
851851 auto photons1_per_collision = photons1.sliceByCached (TDetectorTag1::perCollision (), collision.globalIndex (), cache);
@@ -892,8 +892,8 @@ struct Pi0EtaToGammaGammaMC {
892892 // o2::aod::pwgem::photonmeson::utils::nmhistogram::fillTruePairInfo(&fRegistry, v12, etamc, mcparticles, mccollisions, f1fd_k0s_to_pi0, weight);
893893 // }
894894 } // end of pairing loop
895- } // end of pairing in same event
896- } // end of collision loop
895+ } // end of pairing in same event
896+ } // end of collision loop
897897 }
898898
899899 template <int par_id, typename TBinnedData>
0 commit comments