ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/kiesel/TreeWriter/treeWriter.cc
(Generate patch)

Comparing UserCode/kiesel/TreeWriter/treeWriter.cc (file contents):
Revision 1.20 by kiesel, Wed Apr 24 07:16:49 2013 UTC vs.
Revision 1.21 by kiesel, Wed Apr 24 10:07:51 2013 UTC

# Line 288 | Line 288 | void TreeWriter::Loop() {
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 );
# Line 305 | Line 305 | void TreeWriter::Loop() {
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();

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines