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

Comparing UserCode/MitPhysics/Utils/src/JetIDMVA.cc (file contents):
Revision 1.6 by pharris, Mon Apr 23 15:38:00 2012 UTC vs.
Revision 1.9 by ceballos, Sat May 12 06:17:43 2012 UTC

# Line 15 | Line 15 | using namespace mithep;
15  
16   //--------------------------------------------------------------------------------------------------
17   JetIDMVA::JetIDMVA() :
18 <  fJetPtMin(0)  , //We need to lower this
18 >  fJetPtMin(0.)  , //We need to lower this
19    fDZCut   (0.2),
20    fLowPtMethodName ("JetIDMVALowPt" ),
21    fHighPtMethodName("JetIDMVAHighPt"),
# Line 154 | Line 154 | Bool_t JetIDMVA::pass(const PFJet *iJet,
154                        const PileupEnergyDensityCol *iPileupEnergyDensity) {
155    
156    if(!JetTools::passPFLooseId(iJet))                 return false;
157 <  if(iJet->Pt() < fJetPtMin)                         return false;
158 <  if(fabs(iJet->Eta()) > 4.99)                       return true; //==> Castor
159 <  //if(iJet->Pt() > 50)                                return true; //==> we can raise this
157 >  if(fabs(iJet->Eta()) > 4.99)                       return false;
158    
159    double lMVA = MVAValue   (iJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity);
160    double lPt  = correctedPt(iJet,                  iJetCorrector,iPileupEnergyDensity);
161 <
161 >  if(lPt < fJetPtMin)                         return false;
162 >  
163    int lPtId = 0;
164 <  if(lPt > 10 && iJet->Pt() < 20) lPtId = 1;
165 <  if(lPt > 20 && iJet->Pt() < 30) lPtId = 2;
164 >  if(lPt > 10 && lPt < 20) lPtId = 1;
165 >  if(lPt > 20 && lPt < 30) lPtId = 2;
166    if(lPt > 30                   ) lPtId = 3;
167    
168    int lEtaId = 0;
# Line 184 | Line 183 | Bool_t JetIDMVA::pass(const PFJet *iJet,
183   Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex,const VertexCol *iVertices) {
184    if(!JetTools::passPFLooseId(iJet))                 return false;
185    if(iJet->Pt()        < fJetPtMin) return false;
186 <  //if(iJet->Pt()        > 50)        return true; //==> we can raise this
188 <  if(fabs(iJet->Eta()) > 4.99)     return true; //==> Castor
186 >  if(fabs(iJet->Eta()) > 4.99)      return false;
187    double lMVA = MVAValue(iJet,iVertex,iVertices);
188    
189    int lPtId = 0;

Diff Legend

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