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.6 by pharris, Mon Apr 23 15:38:00 2012 UTC

# Line 86 | Line 86 | void JetIDMVA::Initialize( JetIDMVA::Cut
86    std::string lCutType = "Tight";
87    if(fCutType == kMedium) lCutType = "Medium";
88    if(fCutType == kLoose ) lCutType = "Loose";
89 +  if(fCutType == kMET   ) lCutType = "MET";
90    std::vector<double> lPt010  = lConfig.getParameter<std::vector<double> >(("Pt010_" +lCutType).c_str());
91    std::vector<double> lPt1020 = lConfig.getParameter<std::vector<double> >(("Pt1020_"+lCutType).c_str());
92    std::vector<double> lPt2030 = lConfig.getParameter<std::vector<double> >(("Pt2030_"+lCutType).c_str());
# Line 94 | Line 95 | void JetIDMVA::Initialize( JetIDMVA::Cut
95    for(int i0 = 0; i0 < 4; i0++) fMVACut[1][i0] = lPt1020[i0];
96    for(int i0 = 0; i0 < 4; i0++) fMVACut[2][i0] = lPt2030[i0];
97    for(int i0 = 0; i0 < 4; i0++) fMVACut[3][i0] = lPt3050[i0];
98 +  //std::cout << " Working Points : << " << std::endl;
99 +  //for(int i0 = 0; i0 < 4; i0++) for(int i1 = 0; i1 < 4; i1++)
100 +  //  std::cout << " ==> " << i0 << " -- " << i1 << " -- " << fMVACut[i0][i1] << std::endl;
101   }
102  
103   //--------------------------------------------------------------------------------------------------
# Line 151 | Line 155 | Bool_t JetIDMVA::pass(const PFJet *iJet,
155    
156    if(!JetTools::passPFLooseId(iJet))                 return false;
157    if(iJet->Pt() < fJetPtMin)                         return false;
158 <  if(iJet->Pt() > 50)                                return true; //==> we can raise this
159 <  if(fabs(iJet->Eta()) > 4.99)                       return false;//Castor
160 <
158 >  if(fabs(iJet->Eta()) > 4.99)                       return true; //==> Castor
159 >  //if(iJet->Pt() > 50)                                return true; //==> we can raise this
160 >  
161    double lMVA = MVAValue   (iJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity);
162    double lPt  = correctedPt(iJet,                  iJetCorrector,iPileupEnergyDensity);
163  
# Line 179 | Line 183 | Bool_t JetIDMVA::pass(const PFJet *iJet,
183   //--------------------------------------------------------------------------------------------------
184   Bool_t JetIDMVA::pass(const PFJet *iJet,const Vertex *iVertex,const VertexCol *iVertices) {
185    if(!JetTools::passPFLooseId(iJet))                 return false;
186 <  if(iJet->Pt() < fJetPtMin)                         return false;
187 <  if(iJet->Pt() > 50)                                return true; //==> we can raise this
188 <  if(fabs(iJet->Eta()) > 4.99)                       return false;//Castor
185 <
186 >  if(iJet->Pt()        < fJetPtMin) return false;
187 >  //if(iJet->Pt()        > 50)        return true; //==> we can raise this
188 >  if(fabs(iJet->Eta()) > 4.99)     return true; //==> Castor
189    double lMVA = MVAValue(iJet,iVertex,iVertices);
190    
191    int lPtId = 0;
# Line 237 | Line 240 | Double_t JetIDMVA::MVAValue(const PFJet
240    double lMVA = 0;
241    if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
242    if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
240 –  
243    if (printDebug == kTRUE) {
244      std::cout << "Debug Jet MVA: "
245                << fNVtx      << " "
# Line 294 | Line 296 | Double_t JetIDMVA::MVAValue(const PFJet
296    double lMVA = 0;
297    if(fJPt1 < 10) lMVA = fReader->EvaluateMVA( fLowPtMethodName  );
298    if(fJPt1 > 10) lMVA = fReader->EvaluateMVA( fHighPtMethodName );
299 <
299 >  
300    if (printDebug == kTRUE) {
301      std::cout << "Debug Jet MVA: "
302                << fNVtx      << " "

Diff Legend

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