193 |
|
|
194 |
|
selePtGammaPi0EndCap_region3_ = iConfig.getParameter<double> ("selePtGammaPi0EndCap_region3"); |
195 |
|
selePtPi0EndCap_region3_ = iConfig.getParameter<double> ("selePtPi0EndCap_region3"); |
196 |
< |
|
196 |
> |
selePtPi0MaxEndCap_region3_ = iConfig.getParameter<double> ("selePtPi0MaxEndCap_region3"); |
197 |
|
|
198 |
|
seleS4S9GammaEndCap_ = iConfig.getParameter<double> ("seleS4S9GammaEndCap"); |
199 |
|
seleMinvMaxPi0EndCap_ = iConfig.getParameter<double> ("seleMinvMaxPi0EndCap"); |
251 |
|
|
252 |
|
selePtGammaEtaEndCap_region3_ = iConfig.getParameter<double> ("selePtGammaEtaEndCap_region3"); |
253 |
|
selePtEtaEndCap_region3_ = iConfig.getParameter<double> ("selePtEtaEndCap_region3"); |
254 |
< |
|
254 |
> |
selePtEtaMaxEndCap_region3_ = iConfig.getParameter<double> ("selePtEtaMaxEndCap_region3"); |
255 |
|
|
256 |
+ |
|
257 |
|
seleS4S9GammaEtaEndCap_ = iConfig.getParameter<double> ("seleS4S9GammaEtaEndCap"); |
258 |
|
seleS9S25GammaEtaEndCap_ = iConfig.getParameter<double> ("seleS9S25GammaEtaEndCap"); |
259 |
|
|
881 |
|
|
882 |
|
}catch(std::exception& ex ){ |
883 |
|
nErrorPrinted++; |
884 |
< |
if(nErrorPrinted< maxErrorToPrint) cout<<"encapRecHitsEta NA.."<<endl; |
884 |
> |
if(nErrorPrinted< maxErrorToPrint) cout<<"encapRecHitsPi0 NA.."<<endl; |
885 |
|
} |
886 |
|
} |
887 |
|
|
2001 |
|
double ptpairCut_endcap_region1 = 0 ; |
2002 |
|
double ptpairCut_endcap_region2 = 0 ; |
2003 |
|
double ptpairCut_endcap_region3 = 0 ; |
2004 |
+ |
double ptpairMaxCut_endcap_region3 = 0 ; |
2005 |
+ |
|
2006 |
|
|
2007 |
|
double eta_endcap_region1 = 0 ; |
2008 |
|
double eta_endcap_region2 = 0 ; |
2019 |
|
ptpairCut_endcap_region1 = selePtPi0EndCap_region1_; |
2020 |
|
ptpairCut_endcap_region2 = selePtPi0EndCap_region2_; |
2021 |
|
ptpairCut_endcap_region3 = selePtPi0EndCap_region3_; |
2022 |
< |
|
2022 |
> |
ptpairMaxCut_endcap_region3 = selePtPi0MaxEndCap_region3_; |
2023 |
> |
|
2024 |
> |
|
2025 |
|
eta_endcap_region1 = region1_Pi0EndCap_; |
2026 |
|
eta_endcap_region2 = region2_Pi0EndCap_; |
2027 |
|
|
2034 |
|
ptpairCut_endcap_region1 = selePtEtaEndCap_region1_; |
2035 |
|
ptpairCut_endcap_region2 = selePtEtaEndCap_region2_; |
2036 |
|
ptpairCut_endcap_region3 = selePtEtaEndCap_region3_; |
2037 |
< |
|
2037 |
> |
ptpairMaxCut_endcap_region3 = selePtEtaMaxEndCap_region3_; |
2038 |
> |
|
2039 |
|
eta_endcap_region1 = region1_EtaEndCap_; |
2040 |
|
eta_endcap_region2 = region2_EtaEndCap_; |
2041 |
|
|
2172 |
|
if(ptmin < ptminCut_endcap_region2 || ptpair < ptpairCut_endcap_region2) continue; |
2173 |
|
}else{ |
2174 |
|
if(ptmin < ptminCut_endcap_region3 || ptpair < ptpairCut_endcap_region3) continue; |
2175 |
+ |
if( ptpair > ptpairMaxCut_endcap_region3 ) continue; |
2176 |
|
} |
2177 |
|
}else{ |
2178 |
|
|
2180 |
|
|
2181 |
|
} |
2182 |
|
|
2183 |
< |
|
2183 |
> |
|
2184 |
|
mpair = sqrt ( (eClus[i] + eClus[j])*(eClus[i] + eClus[j]) - (p0x+p1x)*(p0x+p1x) - (p0y+p1y)*(p0y+p1y) - (p0z+p1z)*(p0z+p1z) ); |
2185 |
|
|
2186 |
|
|