ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/interface/MVAMet.h
(Generate patch)

Comparing UserCode/MitPhysics/Utils/interface/MVAMet.h (file contents):
Revision 1.1 by pharris, Wed Mar 21 18:56:25 2012 UTC vs.
Revision 1.4 by pharris, Sat Apr 7 09:36:32 2012 UTC

# Line 22 | Line 22
22   #include "MitCommon/MathTools/interface/MathUtils.h"
23  
24   #include "MitPhysics/Utils/interface/RecoilTools.h"
25 + #include "MitPhysics/Utils/interface/GBRForest.h"
26  
27   class TRandom3;
28  
# Line 34 | Line 35 | namespace mithep {
35        kBaseline = 0
36      };
37  
38 <    void    setVariables(TMVA::Reader *iReader,bool iScale);
39 <    void    Initialize( TString iU1MethodName="U1MVA",
40 <                        TString iPhiMethodName="PhiMVA",
41 <                        TString iJetMVAFile="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilPhiRegress_data_clean2_njet.weights.xml",
42 <                        TString iU1Weights="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilRegress_clean2_njet.weights.xml",
43 <                        TString iPhiWeights="$CMSSW_BASE/src/MitPhysics/data/mva_JetID.weights.xml",
44 <                        MVAType iType=kBaseline);
45 <
38 >    //void    setVariables(TMVA::Reader *iReader,bool iScale);
39 >    //void    Initialize( TString iU1MethodName="U1MVA",
40 >    //                    TString iPhiMethodName="PhiMVA",
41 >    //                  TString iJetMVAFile="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilPhiRegress_baseline.weights.xml",
42 >    //                  TString iU1Weights="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilRegress_baseline.weights.xml",
43 >    //                  TString iPhiWeights="$CMSSW_BASE/src/MitPhysics/data/mva_JetID.weights.xml",
44 >    //                  MVAType iType=kBaseline);
45 >    void    Initialize(
46 >                       TString iJetLowPtFile ="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilPhiRegress_baseline.weights.xml",
47 >                       TString iJetHighPtFile="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilPhiRegress_baseline.weights.xml",
48 >                       TString iJetCutFile   ="$CMSSW_BASE/src/MitPhysics/data/mva_RecoilPhiRegress_baseline.weights.xml",
49 >                       TString iU1Weights    ="$CMSSW_BASE/src/MitPhysics/data/gbrmet.root",
50 >                       TString iPhiWeights   ="$CMSSW_BASE/src/MitPhysics/data/gbrmetphi.root",
51 >                       MVAMet::MVAType  iType=kBaseline);
52 >        
53      Bool_t   IsInitialized() const { return fIsInitialized; }
54 +    Double_t evaluatePhi();
55 +    Double_t evaluateU1();
56      Double_t MVAValue(  bool iPhi,
57                          Float_t iPFSumEt,
58                          Float_t iU      ,
# Line 71 | Line 81 | namespace mithep {
81  
82      Met GetMet(         Bool_t iPhi,Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis,
83                          const PFMet            *iMet  ,
84 <                        const PFCandidateCol   *iCands,const Vertex *iVertex,
84 >                        const PFCandidateCol   *iCands,
85 >                        const Vertex *iVertex,const VertexCol *iVertices,
86                          const PFJetCol         *iJets ,
87                          FactorizedJetCorrector *iJetCorrector,
88                          const PileupEnergyDensityCol *iPileupEnergyDensity,
89                          int iNPV,
90                          Bool_t printDebug=false);
91  
92 +    Met GetMet(         Bool_t iPhi,Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis,
93 +                        const PFMet            *iMet  ,
94 +                        const PFCandidateCol   *iCands,
95 +                        const Vertex *iVertex,const VertexCol *iVertices,
96 +                        const PFJetCol         *iJets ,
97 +                        int iNPV,
98 +                        Bool_t printDebug=false);
99 +
100      Met GetMet(         Bool_t iPhi,
101                          Float_t iPt1,Float_t iPhi1,Float_t iEta1,
102                          Float_t iPt2,Float_t iPhi2,Float_t iEta2,
103                          const PFMet            *iMet  ,
104 <                        const PFCandidateCol   *iCands,const Vertex *iVertex,
104 >                        const PFCandidateCol   *iCands,
105 >                        const Vertex *iVertex,const VertexCol *iVertices,
106                          const PFJetCol         *iJets ,
107                          FactorizedJetCorrector *iJetCorrector,
108                          const PileupEnergyDensityCol *iPUEnergyDensity,
109                          int iNPV,
110                          Bool_t printDebug);
111 +
112 +    Met GetMet(         Bool_t iPhi,
113 +                        Float_t iPt1,Float_t iPhi1,Float_t iEta1,
114 +                        Float_t iPt2,Float_t iPhi2,Float_t iEta2,
115 +                        const PFMet            *iMet  ,
116 +                        const PFCandidateCol   *iCands,
117 +                        const Vertex *iVertex,const VertexCol *iVertices,
118 +                        const PFJetCol         *iJets ,
119 +                        int iNPV,
120 +                        Bool_t printDebug);
121 +    
122      RecoilTools *fRecoilTools;
123 <    protected:
123 >    
124 >  protected:
125      TString      fPhiMethodName;
126      TString      fU1MethodName;
127      Bool_t       fIsInitialized;
# Line 119 | Line 151 | namespace mithep {
151      Float_t fNAllJet;
152      Float_t fNPV    ;
153      Float_t fUPhiMVA;
154 <
155 <    TMVA::Reader* fPhiReader;
156 <    TMVA::Reader* fU1Reader;
154 >    
155 >    Float_t* fPhiVals;
156 >    Float_t* fU1Vals;
157 >    
158 >    
159 >    GBRForest *fPhiReader;
160 >    GBRForest *fU1Reader;
161 >    //TMVA::Reader* fPhiReader;
162 >    //TMVA::Reader* fU1Reader;
163      ClassDef(MVAMet,0)
164    };
165   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines