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.19 by anlevin, Wed Oct 17 01:27:23 2012 UTC vs.
Revision 1.20 by anlevin, Fri Oct 19 18:14:53 2012 UTC

# Line 125 | Line 125 | EventData apply_HZZ4L_reference_selectio
125      for(int i=0; i<muonArr->GetEntries(); i++) {
126        const mithep::Muon *const_mu = (mithep::Muon*)((*muonArr)[i]);
127        mithep::Muon *mu = const_cast<mithep::Muon*>(const_mu);
128 <      correct_muon_momentum(mu,info->RunNum());
128 >      correct_muon_momentum(ctrl, mu,info->RunNum());
129 >
130      }
131    }
132  
# Line 596 | Line 597 | EventData apply_HZZ4L_reference_selectio
597                                    << "\tnewpt: " << lepvec_j[j].vec.Pt()
598                                    << "\tindex: " << j
599                                    << endl;
600 < cout << "FSR TYPEI :: oldpt: " << old_pt_j
600 <                                  << "\tnewpt: " << lepvec_j[j].vec.Pt()
601 <                                  << "\tindex: " << j
602 <                                  << endl;
600 >
601                //lepvec[j].fsrRecoveryAttempted=true;
602              }
603              if( recover_typeII_Photon( ctrl, newmu, j, lepvec_j, pfCandidates ) ) {
# Line 607 | Line 605 | cout << "FSR TYPEI :: oldpt: " << old_pt
605                                    << "\tnewpt: " << lepvec_j[j].vec.Pt()
606                                    << "\tindex: " << j
607                                    << endl;
610 –              cout << "FSR TYPEII :: oldpt: " << old_pt_j
611 –                                  << "\tnewpt: " << lepvec_j[j].vec.Pt()
612 –                                  << "\tindex: " << j
613 –                                  << endl;
608                //lepvec[j].fsrRecoveryAttempted=true;
609              }
610            } else {
# Line 908 | Line 902 | cout << "FSR TYPEI :: oldpt: " << old_pt
902        return ret;
903      }
904      
905 <    if( zzvec.M() > 100. &&  (ret.Z2leptons[0].vec+ret.Z2leptons[1].vec).M() > 12) {
906 <      if(ctrl.debug) cout << "passes mzz > 100, mzz: " << zzvec.M() << "\t mZ2 > 12" << endl;
905 >    if( zzvec.M() > 70. &&  (ret.Z2leptons[0].vec+ret.Z2leptons[1].vec).M() > 12) {
906 >      if(ctrl.debug) cout << "passes mzz > 70, mzz: " << zzvec.M() << "\t mZ2 > 12" << endl;
907        ret.status.selectionBits.flip(PASS_m4l_GT_100);
908 <      cutvec[PASS_m4l_GT_100] +=1;
915 <      zzcutvec[ZZtype][PASS_m4l_GT_100] += 1;
916 <      increment("m4l>100,mZ2>12",theZ1type,theZ2type);
908 >      increment("m4l>70,mZ2>12",theZ1type,theZ2type);
909      } else {
910        ret.status.setStatus(SelectionStatus::FAIL);
911        return ret;
912      }
913  
914 +    if( zzvec.M() > 100. &&  (ret.Z2leptons[0].vec+ret.Z2leptons[1].vec).M() > 12) {
915 +      if(ctrl.debug) cout << "passes mzz > 100, mzz: " << zzvec.M() << "\t mZ2 > 12" << endl;
916 +      cutvec[PASS_m4l_GT_100] +=1;
917 +      zzcutvec[ZZtype][PASS_m4l_GT_100] += 1;
918 +      increment("m4l>100,mZ2>12",theZ1type,theZ2type);
919 +    }
920 +
921      //***************************************************************
922      // finish
923      //***************************************************************

Diff Legend

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