ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/src/JetIDMod.cc
(Generate patch)

Comparing UserCode/MitPhysics/Mods/src/JetIDMod.cc (file contents):
Revision 1.2 by ceballos, Tue Nov 11 21:22:54 2008 UTC vs.
Revision 1.3 by sixie, Fri Nov 21 11:05:39 2008 UTC

# Line 18 | Line 18 | ClassImp(mithep::JetIDMod)
18    fGoodJetsName(Names::gkGoodJetsName),  
19    fJetIDType("HWWJets"),
20    fJets(0),
21 <  fNEventsProcessed(0)
21 >  fNEventsProcessed(0),
22 >  fUseJetCorrection(false),
23 >  fJetEtCut(15.0)
24   {
25    // Constructor.
26   }
# Line 52 | Line 54 | void JetIDMod::Process()
54      Jet *jet = fJets->At(i);
55      
56      const int nCuts = 3;
55 –    double cutValue[nCuts] = {15., 5.0, 0.2};
57      bool passCut[nCuts] = {false, false, false};
58      
59 <    if(jet->Et() > cutValue[0])           passCut[0] = true;
60 <    if(fabs(jet->Eta()) < cutValue[1])    passCut[1] = true;        
61 <    if(jet->Alpha() > cutValue[2] ||
62 <       jet->Et() > 20.)
63 <      passCut[2] = true;
59 >    if(fUseJetCorrection == false) {
60 >      if(jet->Et() > fJetEtCut)             passCut[0] = true;
61 >      if(fabs(jet->Eta()) < 5.0)            passCut[1] = true;        
62 >      if(jet->Alpha() > 0.2 ||
63 >         jet->Et() > 20.)
64 >        passCut[2] = true;
65 >    } else {
66 >      if(jet->Et()*
67 >         jet->L2RelativeCorrectionScale()*
68 >         jet->L3AbsoluteCorrectionScale() > fJetEtCut)            passCut[0] = true;
69 >      if(fabs(jet->Eta()) < 5.0)                                  passCut[1] = true;  
70 >                                                                  passCut[2] = true;
71 >    }
72      
73      // Final decision
74      bool passAllCuts = true;

Diff Legend

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