@@ -66,7 +66,7 @@ struct AssociateMCInfoPhoton {
6666 Configurable<float > max_rxy_gen{" max_rxy_gen" , 100 , " max rxy to store generated information" };
6767 Configurable<bool > requireGammaGammaDecay{" requireGammaGammaDecay" , false , " require gamma gamma decay for generated pi0 and eta meson" };
6868 Configurable<float > cfg_max_eta_photon{" cfg_max_eta_gamma" , 0.8 , " max eta for photons at PV" };
69-
69+
7070 HistogramRegistry registry{" EMMCEvent" };
7171
7272 void init (o2::framework::InitContext&)
@@ -172,19 +172,19 @@ struct AssociateMCInfoPhoton {
172172 for (const auto & mcParticle : groupedMcParticles) { // store necessary information for denominator of efficiency
173173 if ((mcParticle.isPhysicalPrimary () || mcParticle.producedByGenerator ()) && std::fabs (mcParticle.y ()) < 0 .9f && mcParticle.pt () < 20 .f ) {
174174 auto binNumber = hBinFinder->FindBin (mcParticle.pt (), std::fabs (mcParticle.y ())); // caution: pack
175- size_t lNDaughters = 0 ;
175+ size_t lNDaughters = 0 ;
176176 size_t mesonAccepted = 0 ;
177177 if (mcParticle.has_daughters ()) {
178178 auto lDaughters = mcParticle.daughters_as <aod::McParticles>();
179179 lNDaughters = lDaughters.size ();
180180 auto lDaughter0 = lDaughters.begin ();
181181 if (lNDaughters == 2 ) {
182182 auto lDaughter1 = lDaughters.iteratorAt (1 );
183- if ( std::fabs (lDaughter0.eta ())< cfg_max_eta_photon && std::fabs (lDaughter1.eta ()) < cfg_max_eta_photon){
183+ if (std::fabs (lDaughter0.eta ()) < cfg_max_eta_photon && std::fabs (lDaughter1.eta ()) < cfg_max_eta_photon) {
184184 mesonAccepted = 1 ;
185185 }
186- }
187- }
186+ }
187+ }
188188 switch (std::abs (mcParticle.pdgCode ())) {
189189 case PDG_t::kGamma :
190190 registry.fill (HIST (" Generated/h2PtY_Gamma" ), mcParticle.pt (), std::fabs (mcParticle.y ()));
@@ -195,7 +195,7 @@ struct AssociateMCInfoPhoton {
195195 continue ;
196196 registry.fill (HIST (" Generated/h2PtY_Pi0" ), mcParticle.pt (), std::fabs (mcParticle.y ()));
197197 genPi0[binNumber]++;
198- if (mesonAccepted == 1 ){
198+ if (mesonAccepted == 1 ) {
199199 registry.fill (HIST (" Generated/h2PtY_Accepted_Pi0" ), mcParticle.pt (), std::fabs (mcParticle.y ()));
200200 }
201201 break ;
@@ -204,7 +204,7 @@ struct AssociateMCInfoPhoton {
204204 continue ;
205205 registry.fill (HIST (" Generated/h2PtY_Eta" ), mcParticle.pt (), std::fabs (mcParticle.y ()));
206206 genEta[binNumber]++;
207- if (mesonAccepted == 1 ){
207+ if (mesonAccepted == 1 ) {
208208 registry.fill (HIST (" Generated/h2PtY_Accepted_Eta" ), mcParticle.pt (), std::fabs (mcParticle.y ()));
209209 }
210210 break ;
0 commit comments