@@ -93,7 +93,7 @@ struct statPromptPhoton {
9393 Configurable<float > cfgEMChighPSphi{" cfgEMChighPSphi" , 3.26 , " higher limit of the EMC acceptance if Rec PS is required" };
9494 Configurable<float > cfgEMChighPSeta{" cfgEMChighPSeta" , 0.62 , " symmetric eta cut if Rec PS is required" };
9595 Configurable<float > cfgDClowPSphi{" cfgDClowPSphi" , 4.56 , " lower limit of the DCal acceptance if Rec PS is required" };
96- Configurable<float > cfgDChighPSphi{" cfgDChighPSphi" , 5.70 , " higher limit of the DCal acceptance if Rec PS is required" };
96+ Configurable<float > cfgDChighPSphi{" cfgDChighPSphi" , 5.70 , " higher limit of the DCal acceptance if Rec PS is required" };
9797 Configurable<int > cfgMCptNbins{" cfgMCptNbins" , 200 , " number of ptbins in MC QA plots" };
9898 Configurable<double > cfgMCptbinLow{" cfgMCptbinLow" , 5 , " lower limit of ptbins in MC QA plots" };
9999 Configurable<double > cfgMCptbinHigh{" cfgMCptbinHigh" , 200 , " upper limit of ptbins in MC QA plots" };
@@ -113,7 +113,6 @@ struct statPromptPhoton {
113113
114114 const AxisSpec MCptAxis = {cfgMCptNbins, cfgMCptbinLow, cfgMCptbinHigh};
115115
116-
117116 triggerMaskBits = jetderiveddatautilities::initialiseTriggerMaskBits (cfgTriggerMasks);
118117 if (cfgJETracks) {
119118 trackFilter = jetderiveddatautilities::initialiseTrackSelection (static_cast <std::string>(cfgTrackFilter));
@@ -1526,16 +1525,15 @@ struct statPromptPhoton {
15261525 continue ;
15271526 if (std::fabs (mcParticle.getGenStatusCode ()) >= 81 || !mcParticle.isPhysicalPrimary ())
15281527 continue ;
1529- if (cfgReqRecPS_GEN){
1530- if (std::fabs (mcParticle.eta ())>cfgEMChighPSeta)
1531- continue ;
1532- bool insideCalPhi = false ;
1533- if ((mcParticle.phi ()> cfgEMClowPSphi && mcParticle.phi ()<cfgEMChighPSphi) || (mcParticle.phi ()> cfgDClowPSphi && mcParticle.phi ()<cfgDChighPSphi))
1534- insideCalPhi=true ;
1535- if (!insideCalPhi)
1536- continue ;
1537- }
1538-
1528+ if (cfgReqRecPS_GEN) {
1529+ if (std::fabs (mcParticle.eta ()) > cfgEMChighPSeta)
1530+ continue ;
1531+ bool insideCalPhi = false ;
1532+ if ((mcParticle.phi () > cfgEMClowPSphi && mcParticle.phi () < cfgEMChighPSphi) || (mcParticle.phi () > cfgDClowPSphi && mcParticle.phi () < cfgDChighPSphi))
1533+ insideCalPhi = true ;
1534+ if (!insideCalPhi)
1535+ continue ;
1536+ }
15391537
15401538 // Chase this final-state photon upward
15411539 int chaseindex = -1 ;
@@ -1651,17 +1649,17 @@ struct statPromptPhoton {
16511649 if (clusterparticle.getGenStatusCode () < 0 )
16521650 continue ;
16531651 if (std::fabs (clusterparticle.getGenStatusCode ()) >= 81 )
1654- continue ;
1655- if (cfgReqRecPS_REC){
1656- if (std::fabs (clusterparticle.eta ())> cfgEMChighPSeta)
1657- continue ;
1658- bool insideCalPhi = false ;
1659- if ((clusterparticle.phi ()> cfgEMClowPSphi && clusterparticle.phi ()< cfgEMChighPSphi) || (clusterparticle.phi ()> cfgDClowPSphi && clusterparticle.phi ()< cfgDChighPSphi))
1660- insideCalPhi= true ;
1661- if (!insideCalPhi)
1662- continue ;
1663- }
1664-
1652+ continue ;
1653+ if (cfgReqRecPS_REC) {
1654+ if (std::fabs (clusterparticle.eta ()) > cfgEMChighPSeta)
1655+ continue ;
1656+ bool insideCalPhi = false ;
1657+ if ((clusterparticle.phi () > cfgEMClowPSphi && clusterparticle.phi () < cfgEMChighPSphi) || (clusterparticle.phi () > cfgDClowPSphi && clusterparticle.phi () < cfgDChighPSphi))
1658+ insideCalPhi = true ;
1659+ if (!insideCalPhi)
1660+ continue ;
1661+ }
1662+
16651663 int chaseindex = -1 ;
16661664 for (auto & mom : clusterparticle.mothers_as <aod::JMcParticles>()) {
16671665 chaseindex = mom.globalIndex ();
0 commit comments