288 |
|
|
289 |
|
for(std::vector<susy::Photon>::iterator it = photonVector.begin(); |
290 |
|
it != photonVector.end(); ++it ) { |
291 |
< |
if( !it->isEB() && skim ) |
291 |
> |
if( !(it->isEE() || it->isEB()) it->isEBEtaGap() && it->isEBPhiGap() && it->isEERingGap() && it->isEEDeeGap() && it->isEBEEGap() && skim ) |
292 |
|
continue; |
293 |
|
tree::Photon thisphoton; |
294 |
|
thisphoton.pt = getPtFromMatchedJet( *it, jetVector, loggingVerbosity ); |
305 |
|
&& thisphoton.chargedIso<2.6 |
306 |
|
&& thisphoton.neutralIso<3.5+0.04*thisphoton.pt |
307 |
|
&& thisphoton.photonIso<1.3+0.005*thisphoton.pt; |
308 |
< |
/*bool loose_photon_endcap = thisphoton.pt > 20 |
308 |
> |
bool loose_photon_endcap = thisphoton.pt > 20 |
309 |
|
&& it->isEE() |
310 |
|
&& it->passelectronveto |
311 |
|
&& it->hadTowOverEm<0.05 |
312 |
|
&& it->sigmaIetaIeta<0.034 |
313 |
|
&& thisphoton.chargedIso<2.3 |
314 |
|
&& thisphoton.neutralIso<2.9+0.04*thisphoton.pt; |
315 |
< |
*/ |
316 |
< |
if(!(loose_photon_barrel || thisphoton.pt > 75 ) && skim ) |
315 |
> |
|
316 |
> |
if(!(loose_photon_endcap || loose_photon_barrel || thisphoton.pt > 75 ) && skim ) |
317 |
|
continue; |
318 |
|
thisphoton.eta = it->momentum.Eta(); |
319 |
|
thisphoton.phi = it->momentum.Phi(); |