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

Comparing UserCode/MitHzz4l/Selection/src/applySelection.cc (file contents):
Revision 1.32 by anlevin, Sat Jun 2 12:17:29 2012 UTC vs.
Revision 1.33 by khahn, Sun Jun 3 15:59:56 2012 UTC

# Line 59 | Line 59 | using namespace std;
59   #include "InfoStruct.h"
60   #include "GenInfoStruct.h"
61   #include "WeightStruct.h"
62 //#include "GlobalDataAndFuncs.h"
63 #include "RunLumiRangeMap.h"
64
62   #include "AngleTuple.h"
63   #include "FOTuple.h"
64  
65 + #include "RunLumiRangeMap.h"
66   #include "SampleWeight.h"
67   #include "EfficiencyWeightsInterface.h"
68  
# Line 219 | Line 217 | int main(int argc, char** argv)
217    nt.makeKinematicsBranch(kinematics);
218    nt.makeInfoBranch(evtinfo);
219  
220 +  
221    FOTuple foTree( nt.getFile(), (const char*)"FOtree");
222    foTree.makeFailedLeptonBranch(failingLeptons);
223 <  foTree.makeZLeptonBranch(passingLeptons);
223 >  foTree.makePassedLeptonBranch(passingLeptons);
224  
225    TH1F * h_wf11_hpt;
226    if(ctrl.mc) {
227      nt.makeWeightBranch(weights);
228 <    nt.makeGenInfoBranch(geninfo);
228 >    //    nt.makeGenInfoBranch(geninfo);
229      initEfficiencyWeights();
230      initPUWeights();
231  
# Line 245 | Line 244 | int main(int argc, char** argv)
244    }
245  
246    //  initMuonIDMVA();
247 <  initMuonIsoMVA();
249 <  initElectronIDMVA();
250 <  initElectronIsoMVA();
247 >  initElectronIDMVAV1();
248    initTrigger();
249  
250    // hack
# Line 392 | Line 389 | int main(int argc, char** argv)
389            fillTriggerNames(hltTable, hltTableStrings );
390          }
391          if( ctrl.debug ) cout << "file is : " << currentFile  << endl;
392 <        fillTriggerBits( hltTable, trigMask, triggerBits );
392 >        fillTriggerNames( hltTable, trigMask, triggerBits );
393          passes_HLT_MC = true; // passed to apply as extern global, just for sync
394          if( !passHLTMC(triggerBits, info->RunNum(), info->EvtNum(), k2012_MC ) ) {
395            passes_HLT_MC = false;
# Line 420 | Line 417 | int main(int argc, char** argv)
417          }
418          if( ctrl.debug ) cout << "file is : " << currentFile  << endl;
419          /*
420 <        fillTriggerBits( hltTable, trigMask, triggerBits );
420 >        fillTriggerNames( hltTable, trigMask, triggerBits );
421          if( !passHLT(triggerBits, info->RunNum(), info->EvtNum() ) ) {
422            if( ctrl.debug ) cout << "\tfails trigger ... " << endl;
423            continue;
# Line 438 | Line 435 | int main(int argc, char** argv)
435  
436          // reference
437          apply_HZZ4L_reference_selection(ctrl, info,
438 <                                         vtxArr,
439 <                                         pfArr,
438 >                                        vtxArr,
439 >                                        pfArr,
440   #ifdef HACKED_RHOS
441 <                                         rho,
441 >                                        rho,
442   #else
443 <                                         puDArr,
443 >                                        puDArr,
444   #endif
445 <                                         electronArr,
446 <                                         &electronReferencePreSelection,
447 <                                         &electronReferenceIDMVASelectionV1,
448 <                                         &electronReferenceIsoSelection,
449 <                                         muonArr,
450 <                                         &muonReferencePreSelection,
451 <                                         &muonIDPFSelection,
452 <                                         &muonReferenceIsoSelection);
453 <
457 <      /*
458 <        // evolved
459 <        apply_HZZ4L_reference2_selection(ctrl, info,
460 <                                         vtx,
461 <                                         pfArr,
462 <                                         //rho,
463 <                                         puDArr,
464 <                                         electronArr,
465 <                                         &electronReferencePreSelection,
466 <                                         //&electronPreSelection,
467 <                                         &electronReferenceIDMVASelection,
468 <                                         &electronIsoMVASelection,
469 <                                         muonArr,
470 <                                         &muonReferencePreSelection,
471 <                                         //&muonPreSelection,
472 <                                         &muonIDPFSelection,
473 <                                         &muonIsoMVASelection);
474 <      */
475 <
445 >                                        electronArr,
446 >                                        &electronReferencePreSelection,
447 >                                        &electronReferenceIDMVASelectionV1,
448 >                                        &electronReferenceIsoSelection,
449 >                                        muonArr,
450 >                                        &muonReferencePreSelection,
451 >                                        &muonIDPFSelection,
452 >                                        &muonReferenceIsoSelection);
453 >      
454        if( ctrl.debug ) cout << endl;
455  
456  
# Line 521 | Line 499 | int main(int argc, char** argv)
499          pass++;
500          //      if( pass > 3 ) break;
501  
502 <      } else if( ret4l.status.selectionBits.test(PASS_GOODZ1) ) { // save for fakes ...
502 >      } else if( ret4l.status.selectionBits.test(PASS_GOODZ1) &&
503 >                 ret4l.status.selectionBits.to_ulong() < PASS_4L) { // save for fakes ...
504          //      if(ctrl.debug) {
505            cout << "failingLeptons : " << failingLeptons.size() << endl;
506            for( int f=0; f<failingLeptons.size(); f++ ) {

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines