ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/Selection/src/ReferenceSelection.cc
(Generate patch)

Comparing UserCode/MitHzz4l/Selection/src/ReferenceSelection.cc (file contents):
Revision 1.12 by anlevin, Sun Jul 1 13:03:43 2012 UTC vs.
Revision 1.13 by dkralph, Fri Jul 6 15:27:57 2012 UTC

# Line 233 | Line 233 | EventData apply_HZZ4L_reference_selectio
233    //----------------------------------------------------
234    for(int i=0; i<muonArr->GetEntries(); i++)
235      {
236 <      const mithep::Muon *mu = (mithep::Muon*)((*muonArr)[i]);
237 <
236 >      const mithep::Muon *mu = (mithep::Muon*)((*muonArr)[i]);      
237 >      
238        SelectionStatus musel;
239        if(ctrl.debug) cout << "musel.status  before anything: " << musel.getStatus() << endl;
240  
# Line 460 | Line 460 | EventData apply_HZZ4L_reference_selectio
460        } else {
461          passingLeptons.push_back(lepvec[i]);
462        }
463 – #ifndef SYNC
463        if( !already_pushed && fabs(lepvec[i].ip3dSig)>4 )  
464 <        pt_of_leptons_to_erase.push_back(lepvec[i].vec.Pt());
466 < #endif      
464 >        pt_of_leptons_to_erase.push_back(lepvec[i].vec.Pt());
465      }
466 <    for( int i=0; i<pt_of_leptons_to_erase.size(); i++ ) {
466 >   for( int i=0; i<pt_of_leptons_to_erase.size(); i++ ) {
467        for( vector<SimpleLepton>::iterator it=lepvec.begin();
468             it != lepvec.end(); it++ ) {
469          SimpleLepton flep = *it;
# Line 496 | Line 494 | EventData apply_HZZ4L_reference_selectio
494          vector<SimpleLepton> lepvec_i = lepvec;
495          vector<SimpleLepton> lepvec_j = lepvec;
496  
497 <        if( ctrl.doFSR ) {
498 <          cout << endl;
499 <          cout << "----------------> FSR ("<<i<<","<<j<<") <----------------------" << endl;
497 >        if( ctrl.doFSR ) {
498 >            cout << endl;
499 >            cout << "----------------> FSR ("<<i<<","<<j<<") <----------------------" << endl;
500            photonsToVeto.clear();
501            float old_pt_i = lepvec[i].vec.Pt();
502            float old_pt_j = lepvec[j].vec.Pt();
# Line 551 | Line 549 | EventData apply_HZZ4L_reference_selectio
549              if( recover_typeII_Photon( ctrl, newmu, j, lepvec_j, pfCandidates ) ) {
550                if(ctrl.debug) cout << "FSR TYPEII :: oldpt: " << old_pt_j
551                                    << "\tnewpt: " << lepvec_j[j].vec.Pt()
554 –
552                                    << "\tindex: " << j
553                                    << endl;
554                //lepvec[j].fsrRecoveryAttempted=true;
# Line 836 | Line 833 | EventData apply_HZZ4L_reference_selectio
833        (ZCandidatesLeptons[best_Z2_index].first.vec) +
834        (ZCandidatesLeptons[best_Z2_index].second.vec);
835      
836 <    if( zzvec.M() > 70. ) {
836 >    if( zzvec.M() > 70. ) {
837 >      if(ctrl.debug) cout << "passes mzz > 70, mzz: " << zzvec.M() << endl;
838        ret.status.selectionBits.flip(PASS_m4l_GT_70);
839        cutvec[PASS_m4l_GT_70] +=1;
840        zzcutvec[ZZtype][PASS_m4l_GT_70] += 1;
# Line 846 | Line 844 | EventData apply_HZZ4L_reference_selectio
844      }
845      
846      if( zzvec.M() > 100. &&  (ret.Z2leptons[0].vec+ret.Z2leptons[1].vec).M() > 12) {
847 +      if(ctrl.debug) cout << "passes mzz > 100, mzz: " << zzvec.M() << "\t mZ2 > 12" << endl;
848        ret.status.selectionBits.flip(PASS_m4l_GT_100);
849        cutvec[PASS_m4l_GT_100] +=1;
850        zzcutvec[ZZtype][PASS_m4l_GT_100] += 1;

Diff Legend

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