Skip to content

Commit 43d5edc

Browse files
committed
Bugfixes for prior commit
1 parent 69c306f commit 43d5edc

1 file changed

Lines changed: 10 additions & 5 deletions

File tree

PWGJE/Tasks/statPromptPhoton.cxx

Lines changed: 10 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -89,7 +89,12 @@ struct statPromptPhoton {
8989
Configurable<bool> cfgGenReqRec{"cfgGenReqRec", false, "Only consider generated events which are successfully reconstructed"};
9090
Configurable<bool> cfgReqRecPS_REC{"cfgReqRecPS_REC", false, "Only consider reconstructed photons within the EMCAl acceptence"};
9191
Configurable<bool> cfgReqRecPS_GEN{"cfgReqRecPS_GEN", false, "Only consider generated photons within the EMCAl acceptence"};
92-
Configurable<double> cfgMCptNbins{"cfgMCptNbins", 200, "number of ptbins in MC QA plots"};
92+
Configurable<float> cfgEMClowPSphi{"cfgEMClowPSphi", 1.42, "lower limit of the EMC acceptance if Rec PS is required"};
93+
Configurable<float> cfgEMChighPSphi{"cfgEMChighPSphi", 3.26, "higher limit of the EMC acceptance if Rec PS is required"};
94+
Configurable<float> cfgEMChighPSeta{"cfgEMChighPSeta", 0.62, "symmetric eta cut if Rec PS is required"};
95+
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"};
97+
Configurable<int> cfgMCptNbins{"cfgMCptNbins", 200, "number of ptbins in MC QA plots"};
9398
Configurable<double> cfgMCptbinLow{"cfgMCptbinLow", 5, "lower limit of ptbins in MC QA plots"};
9499
Configurable<double> cfgMCptbinHigh{"cfgMCptbinHigh", 200, "upper limit of ptbins in MC QA plots"};
95100
Configurable<bool> cfgRecHistograms{"cfgRecHistograms", false, "Enables Reconstructed histograms"};
@@ -1522,10 +1527,10 @@ struct statPromptPhoton {
15221527
if (std::fabs(mcParticle.getGenStatusCode()) >= 81 || !mcParticle.isPhysicalPrimary())
15231528
continue;
15241529
if(cfgReqRecPS_GEN){
1525-
if(std::fabs(mcParticle.eta())>0.62)
1530+
if(std::fabs(mcParticle.eta())>cfgEMChighPSeta)
15261531
continue;
15271532
bool insideCalPhi = false;
1528-
if((mcParticle.phi()> 1.42 && mcParticle.phi()<3.26) || (mcParticle.phi()> 4.56 && mcParticle.phi()<5.70))
1533+
if((mcParticle.phi()> cfgEMClowPSphi && mcParticle.phi()<cfgEMChighPSphi) || (mcParticle.phi()> cfgDClowPSphi && mcParticle.phi()<cfgDChighPSphi))
15291534
insideCalPhi=true;
15301535
if(!insideCalPhi)
15311536
continue;
@@ -1648,10 +1653,10 @@ struct statPromptPhoton {
16481653
if (std::fabs(clusterparticle.getGenStatusCode()) >= 81)
16491654
continue;
16501655
if(cfgReqRecPS_REC){
1651-
if(std::fabs(clusterparticle.eta())>0.62)
1656+
if(std::fabs(clusterparticle.eta())>cfgEMChighPSeta)
16521657
continue;
16531658
bool insideCalPhi = false;
1654-
if((clusterparticle.phi()> 1.42 && clusterparticle.phi()<3.26) || (clusterparticle.phi()> 4.56 && clusterparticle.phi()<5.70))
1659+
if((clusterparticle.phi()> cfgEMClowPSphi && clusterparticle.phi()<cfgEMChighPSphi) || (clusterparticle.phi()> cfgDClowPSphi && clusterparticle.phi()<cfgDChighPSphi))
16551660
insideCalPhi=true;
16561661
if(!insideCalPhi)
16571662
continue;

0 commit comments

Comments
 (0)