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.4 by pharris, Sat Apr 7 09:36:33 2012 UTC vs.
Revision 1.5 by pharris, Fri Apr 13 14:24:00 2012 UTC

# Line 94 | Line 94 | void JetIDMVA::Initialize( JetIDMVA::Cut
94    for(int i0 = 0; i0 < 4; i0++) fMVACut[1][i0] = lPt1020[i0];
95    for(int i0 = 0; i0 < 4; i0++) fMVACut[2][i0] = lPt2030[i0];
96    for(int i0 = 0; i0 < 4; i0++) fMVACut[3][i0] = lPt3050[i0];
97 +  //std::cout << " Working Points : << " << std::endl;
98 +  //for(int i0 = 0; i0 < 4; i0++) for(int i1 = 0; i1 < 4; i1++)
99 +  //  std::cout << " ==> " << i0 << " -- " << i1 << " -- " << fMVACut[i0][i1] << std::endl;
100   }
101  
102   //--------------------------------------------------------------------------------------------------
# Line 151 | Line 154 | Bool_t JetIDMVA::pass(const PFJet *iJet,
154    
155    if(!JetTools::passPFLooseId(iJet))                 return false;
156    if(iJet->Pt() < fJetPtMin)                         return false;
157 <  if(iJet->Pt() > 50)                                return true; //==> we can raise this
158 <  if(fabs(iJet->Eta()) > 4.99)                       return false;//Castor
159 <
157 >  if(fabs(iJet->Eta()) > 4.99)                       return true; //==> Castor
158 >  //if(iJet->Pt() > 50)                                return true; //==> we can raise this
159 >  
160    double lMVA = MVAValue   (iJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity);
161    double lPt  = correctedPt(iJet,                  iJetCorrector,iPileupEnergyDensity);
162  
# Line 179 | Line 182 | Bool_t JetIDMVA::pass(const PFJet *iJet,
182   //--------------------------------------------------------------------------------------------------
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
187 <  if(fabs(iJet->Eta()) > 4.99)                       return false;//Castor
185 <
185 >  if(iJet->Pt()        < fJetPtMin) return false;
186 >  //if(iJet->Pt()        > 50)        return true; //==> we can raise this
187 >  if(fabs(iJet->Eta()) > 4.99)     return true; //==> Castor
188    double lMVA = MVAValue(iJet,iVertex,iVertices);
189    
190    int lPtId = 0;
# Line 237 | Line 239 | Double_t JetIDMVA::MVAValue(const PFJet
239    double lMVA = 0;
240    if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
241    if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
240 –  
242    if (printDebug == kTRUE) {
243      std::cout << "Debug Jet MVA: "
244                << fNVtx      << " "
# Line 294 | Line 295 | Double_t JetIDMVA::MVAValue(const PFJet
295    double lMVA = 0;
296    if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
297    if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
298 <
298 >  
299    if (printDebug == kTRUE) {
300      std::cout << "Debug Jet MVA: "
301                << fNVtx      << " "

Diff Legend

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