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.18 by kiesel, Tue Apr 23 11:32:05 2013 UTC vs.
Revision 1.20 by kiesel, Wed Apr 24 07:16:49 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() || it->isEE()) && skim )
291 >                        if( !it->isEB() && 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 <                        if(!(loose_photon_endcap || loose_photon_barrel || thisphoton.pt > 75 ) && skim )
315 >                        */
316 >                        if(!(loose_photon_barrel || thisphoton.pt > 75 ) && skim )
317                                  continue;
318                          thisphoton.eta = it->momentum.Eta();
319                          thisphoton.phi = it->momentum.Phi();
# Line 408 | Line 409 | void TreeWriter::Loop() {
409                          if( it->momentum.Pt() < 20  || it->momentum.Pt() > 1e6 )
410                                  continue; // spike rejection
411                          float iso = ( it->chargedHadronIso + max(it->neutralHadronIso+it->photonIso
412 <                                                                                                                - effectiveAreaElectron(it->momentum.Eta())*event->rho25, (Float_t)0. )
412 >                                                                                                                - effectiveAreaElectron(it->momentum.Eta())*event->rho25, (float)0. )
413                                                  ) / it->momentum.Pt();
414                          float d0 = d0correction( *it, *event );
415                          float dZ = std::abs( dZcorrection( *it, *event ) );

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines