CMS 3D CMS Logo
List of all members | Classes | Public Types | Public Member Functions | Public Attributes | Private Attributes
PhotonMVAEstimatorRun2Spring15NonTrig Class Reference

#include <PhotonMVAEstimatorRun2Spring15NonTrig.h>

Inheritance diagram for PhotonMVAEstimatorRun2Spring15NonTrig:
AnyMVAEstimatorRun2Base

Classes

struct  AllVariables
 

Public Types

enum  mvaCategories { UNDEFINED = -1, CAT_EB = 0, CAT_EE = 1 }
 

Public Member Functions

void constrainMVAVariables (AllVariables &) const
 
std::unique_ptr< const GBRForest > createSingleReader (const int iCategory, const edm::FileInPath &weightFile)
 
std::vector< float > fillMVAVariables (const edm::Ptr< reco::Candidate > &particle, const edm::Event &iEvent) const override
 
int findCategory (const edm::Ptr< reco::Candidate > &particle) const override
 
const std::string & getName () const final
 
int getNCategories () const override
 
const std::string & getTag () const final
 
bool isEndcapCategory (int category) const
 
float mvaValue (const edm::Ptr< reco::Candidate > &particle, const edm::Event &) const override
 
 PhotonMVAEstimatorRun2Spring15NonTrig (const edm::ParameterSet &conf)
 
void setConsumes (edm::ConsumesCollector &&) const override
 
 ~PhotonMVAEstimatorRun2Spring15NonTrig () override
 
- Public Member Functions inherited from AnyMVAEstimatorRun2Base
 AnyMVAEstimatorRun2Base (const edm::ParameterSet &conf)
 
virtual void getEventContent (const edm::Event &iEvent) const final
 
template<typename... Args>
std::vector< float > packMVAVariables (const Args...args) const
 
virtual ~AnyMVAEstimatorRun2Base ()
 

Public Attributes

const int nCategories = 2
 
- Public Attributes inherited from AnyMVAEstimatorRun2Base
const edm::ParameterSet _conf
 

Private Attributes

AllVariables _allMVAVars
 
const edm::InputTag _esEffSigmaRRMapLabel
 
const edm::InputTag _full5x5E1x3MapLabel
 
const edm::InputTag _full5x5E2x2MapLabel
 
const edm::InputTag _full5x5E2x5MaxMapLabel
 
const edm::InputTag _full5x5E5x5MapLabel
 
const edm::InputTag _full5x5SigmaIEtaIEtaMapLabel
 
const edm::InputTag _full5x5SigmaIEtaIPhiMapLabel
 
std::vector< std::unique_ptr< const GBRForest > > _gbrForests
 
const std::string _MethodName
 
const std::string _name = "PhotonMVAEstimatorRun2Spring15NonTrig"
 
const edm::InputTag _phoChargedIsolationLabel
 
const edm::InputTag _phoPhotonIsolationLabel
 
const edm::InputTag _phoWorstChargedIsolationLabel
 
const edm::InputTag _rhoLabel
 
std::string _tag
 
const bool _useValueMaps
 

Detailed Description

Definition at line 21 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Member Enumeration Documentation

enum PhotonMVAEstimatorRun2Spring15NonTrig::mvaCategories
Enumerator
UNDEFINED 
CAT_EB 
CAT_EE 

Definition at line 28 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

28  {
29  UNDEFINED = -1,
30  CAT_EB = 0,
31  CAT_EE = 1
32  };

Constructor & Destructor Documentation

PhotonMVAEstimatorRun2Spring15NonTrig::PhotonMVAEstimatorRun2Spring15NonTrig ( const edm::ParameterSet &  conf)

Definition at line 7 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References _gbrForests, _tag, createSingleReader(), Exception, edm::ParameterSet::getParameter(), mps_fire::i, nCategories, AlCaHLTBitMon_QueryRunRegistry::string, electrons_cff::weightFile, mvaElectronID_PHYS14_PU20bx25_nonTrig_V1_cff::weightFileNames, and ~PhotonMVAEstimatorRun2Spring15NonTrig().

7  :
8  AnyMVAEstimatorRun2Base(conf),
9  _MethodName("BDTG method"),
10  _useValueMaps(conf.getParameter<bool>("useValueMaps")),
11  _full5x5SigmaIEtaIEtaMapLabel(conf.getParameter<edm::InputTag>("full5x5SigmaIEtaIEtaMap")),
12  _full5x5SigmaIEtaIPhiMapLabel(conf.getParameter<edm::InputTag>("full5x5SigmaIEtaIPhiMap")),
13  _full5x5E1x3MapLabel(conf.getParameter<edm::InputTag>("full5x5E1x3Map")),
14  _full5x5E2x2MapLabel(conf.getParameter<edm::InputTag>("full5x5E2x2Map")),
15  _full5x5E2x5MaxMapLabel(conf.getParameter<edm::InputTag>("full5x5E2x5MaxMap")),
16  _full5x5E5x5MapLabel(conf.getParameter<edm::InputTag>("full5x5E5x5Map")),
17  _esEffSigmaRRMapLabel(conf.getParameter<edm::InputTag>("esEffSigmaRRMap")),
18  _phoChargedIsolationLabel(conf.getParameter<edm::InputTag>("phoChargedIsolation")),
19  _phoPhotonIsolationLabel(conf.getParameter<edm::InputTag>("phoPhotonIsolation")),
20  _phoWorstChargedIsolationLabel(conf.getParameter<edm::InputTag>("phoWorstChargedIsolation")),
21  _rhoLabel(conf.getParameter<edm::InputTag>("rho"))
22 {
23 
24  //
25  // Construct the MVA estimators
26  //
27  _tag = conf.getParameter<std::string>("mvaTag");
28 
29  const std::vector <std::string> weightFileNames
30  = conf.getParameter<std::vector<std::string> >("weightFileNames");
31 
32  if( (int)(weightFileNames.size()) != nCategories )
33  throw cms::Exception("MVA config failure: ")
34  << "wrong number of weightfiles" << std::endl;
35 
36  _gbrForests.clear();
37  // The method name is just a key to retrieve this method later, it is not
38  // a control parameter for a reader (the full definition of the MVA type and
39  // everything else comes from the xml weight files).
40 
41  // Create a TMVA reader object for each category
42  for(int i=0; i<nCategories; i++){
43 
44  // Use unique_ptr so that all readers are properly cleaned up
45  // when the vector clear() is called in the destructor
46 
47  edm::FileInPath weightFile( weightFileNames[i] );
48  _gbrForests.push_back( createSingleReader(i, weightFile ) );
49 
50  }
51 
52 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AnyMVAEstimatorRun2Base::AnyMVAEstimatorRun2Base
AnyMVAEstimatorRun2Base(const edm::ParameterSet &conf)
Definition: AnyMVAEstimatorRun2Base.h:15
mps_fire.i
i
Definition: mps_fire.py:269
PhotonMVAEstimatorRun2Spring15NonTrig::createSingleReader
std::unique_ptr< const GBRForest > createSingleReader(const int iCategory, const edm::FileInPath &weightFile)
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.cc:142
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
mvaElectronID_PHYS14_PU20bx25_nonTrig_V1_cff.weightFileNames
weightFileNames
Definition: mvaElectronID_PHYS14_PU20bx25_nonTrig_V1_cff.py:90
electrons_cff.weightFile
weightFile
Definition: electrons_cff.py:140
PhotonMVAEstimatorRun2Spring15NonTrig::_gbrForests
std::vector< std::unique_ptr< const GBRForest > > _gbrForests
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.h:105
PhotonMVAEstimatorRun2Spring15NonTrig::~PhotonMVAEstimatorRun2Spring15NonTrig ( )
override

Definition at line 55 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References mvaValue().

Referenced by PhotonMVAEstimatorRun2Spring15NonTrig().

55  {
56 }

Member Function Documentation

void PhotonMVAEstimatorRun2Spring15NonTrig::constrainMVAVariables ( AllVariables &  ) const

Definition at line 366 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

Referenced by fillMVAVariables(), and getTag().

366  {
367 
368  // Check that variables do not have crazy values
369 
370  // This function is currently empty as this specific MVA was not
371  // developed with restricting variables to specific physical ranges.
372  return;
373 
374 }
std::unique_ptr< const GBRForest > PhotonMVAEstimatorRun2Spring15NonTrig::createSingleReader ( const int  iCategory,
const edm::FileInPath &  weightFile 
)

Definition at line 142 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References _allMVAVars, _MethodName, edm::FileInPath::fullPath(), isEndcapCategory(), PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varChIsoRaw, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE1x3overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE2x2overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE2x5overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varESEffSigmaRR, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varESEnOverRawE, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varEta, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPhi, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPhoIsoRaw, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPt, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varR9, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varRawE, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varRho, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCEta, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCEtaWidth, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCPhiWidth, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSieie, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSieip, and PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varWorstChRaw.

Referenced by isEndcapCategory(), and PhotonMVAEstimatorRun2Spring15NonTrig().

142  {
143 
144  //
145  // Create the reader
146  //
147  TMVA::Reader tmpTMVAReader( "!Color:Silent:Error" );
148 
149  //
150  // Configure all variables and spectators. Note: the order and names
151  // must match what is found in the xml weights file!
152  //
153 
154  tmpTMVAReader.AddVariable("recoPhi" , &_allMVAVars.varPhi);
155  tmpTMVAReader.AddVariable("r9" , &_allMVAVars.varR9);
156  tmpTMVAReader.AddVariable("sieieFull5x5", &_allMVAVars.varSieie);
157  tmpTMVAReader.AddVariable("sieipFull5x5", &_allMVAVars.varSieip);
158  tmpTMVAReader.AddVariable("e1x3Full5x5/e5x5Full5x5", &_allMVAVars.varE1x3overE5x5);
159  tmpTMVAReader.AddVariable("e2x2Full5x5/e5x5Full5x5", &_allMVAVars.varE2x2overE5x5);
160  tmpTMVAReader.AddVariable("e2x5Full5x5/e5x5Full5x5", &_allMVAVars.varE2x5overE5x5);
161  tmpTMVAReader.AddVariable("recoSCEta" , &_allMVAVars.varSCEta);
162  tmpTMVAReader.AddVariable("rawE" , &_allMVAVars.varRawE);
163  tmpTMVAReader.AddVariable("scEtaWidth", &_allMVAVars.varSCEtaWidth);
164  tmpTMVAReader.AddVariable("scPhiWidth", &_allMVAVars.varSCPhiWidth);
165 
166  // Endcap only variables
167  if( isEndcapCategory(iCategory) ){
168  tmpTMVAReader.AddVariable("esEn/rawE" , &_allMVAVars.varESEnOverRawE);
169  tmpTMVAReader.AddVariable("esRR" , &_allMVAVars.varESEffSigmaRR);
170  }
171 
172  // Pileup
173  tmpTMVAReader.AddVariable("rho" , &_allMVAVars.varRho);
174 
175  // Isolations
176  tmpTMVAReader.AddVariable("phoIsoRaw" , &_allMVAVars.varPhoIsoRaw);
177  tmpTMVAReader.AddVariable("chIsoRaw" , &_allMVAVars.varChIsoRaw);
178  tmpTMVAReader.AddVariable("chWorstRaw", &_allMVAVars.varWorstChRaw);
179 
180  // Spectators
181  tmpTMVAReader.AddSpectator("recoPt" , &_allMVAVars.varPt);
182  tmpTMVAReader.AddSpectator("recoEta", &_allMVAVars.varEta);
183 
184  //
185  // Book the method and set up the weights file
186  //
187  tmpTMVAReader.BookMVA(_MethodName , weightFile.fullPath());
188 
189  return std::make_unique<const GBRForest>(dynamic_cast<TMVA::MethodBDT*>( tmpTMVAReader.FindMVA(_MethodName) ) );
190 
191 }
edm::FileInPath::fullPath
std::string fullPath() const
Definition: FileInPath.cc:184
std::vector< float > PhotonMVAEstimatorRun2Spring15NonTrig::fillMVAVariables ( const edm::Ptr< reco::Candidate > &  particle,
const edm::Event &  iEvent 
) const
overridevirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 194 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References _esEffSigmaRRMapLabel, _full5x5E1x3MapLabel, _full5x5E2x2MapLabel, _full5x5E2x5MaxMapLabel, _full5x5E5x5MapLabel, _full5x5SigmaIEtaIEtaMapLabel, _full5x5SigmaIEtaIPhiMapLabel, _phoChargedIsolationLabel, _phoPhotonIsolationLabel, _phoWorstChargedIsolationLabel, _rhoLabel, _useValueMaps, constrainMVAVariables(), reco::Photon::ShowerShape::e1x3, heepElectronID_HEEPV50_CSA14_25ns_cff::e2x5, heepElectronID_HEEPV50_CSA14_25ns_cff::e5x5, mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::esEffSigmaRRMap, reco::LeafCandidate::eta(), Exception, findCategory(), reco::Photon::full5x5_showerShapeVariables(), mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::full5x5E1x3Map, mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::full5x5E2x2Map, mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::full5x5E2x5MaxMap, mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::full5x5E5x5Map, cutBasedElectronID_CSA14_50ns_V1_cff::full5x5SigmaIEtaIEtaMap, mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff::full5x5SigmaIEtaIPhiMap, edm::Event::getByLabel(), isEndcapCategory(), edm::Ptr< T >::isNull(), edm::HandleBase::isValid(), SiStripPI::max, eostools::move(), AnyMVAEstimatorRun2Base::packMVAVariables(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::Photon::r9(), rho, reco::Photon::superCluster(), PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varChIsoRaw, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE1x3overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE2x2overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varE2x5overE5x5, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varESEffSigmaRR, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varESEnOverRawE, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varEta, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPhi, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPhoIsoRaw, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varPt, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varR9, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varRawE, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varRho, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCEta, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCEtaWidth, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSCPhiWidth, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSieie, PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varSieip, and PhotonMVAEstimatorRun2Spring15NonTrig::AllVariables::varWorstChRaw.

Referenced by getTag(), and mvaValue().

194  {
195 
196  //
197  // Declare all value maps corresponding to the above tokens
198  //
199  edm::Handle<edm::ValueMap<float> > full5x5SigmaIEtaIEtaMap;
200  edm::Handle<edm::ValueMap<float> > full5x5SigmaIEtaIPhiMap;
201  edm::Handle<edm::ValueMap<float> > full5x5E1x3Map;
202  edm::Handle<edm::ValueMap<float> > full5x5E2x2Map;
203  edm::Handle<edm::ValueMap<float> > full5x5E2x5MaxMap;
204  edm::Handle<edm::ValueMap<float> > full5x5E5x5Map;
205  edm::Handle<edm::ValueMap<float> > esEffSigmaRRMap;
206  //
207  edm::Handle<edm::ValueMap<float> > phoChargedIsolationMap;
208  edm::Handle<edm::ValueMap<float> > phoPhotonIsolationMap;
209  edm::Handle<edm::ValueMap<float> > phoWorstChargedIsolationMap;
210 
211  // Rho will be pulled from the event content
212  edm::Handle<double> rho;
213 
214  // Get the full5x5 and ES maps
215  if( _useValueMaps ) {
216  iEvent.getByLabel(_full5x5SigmaIEtaIEtaMapLabel, full5x5SigmaIEtaIEtaMap);
217  iEvent.getByLabel(_full5x5SigmaIEtaIPhiMapLabel, full5x5SigmaIEtaIPhiMap);
218  iEvent.getByLabel(_full5x5E1x3MapLabel, full5x5E1x3Map);
219  iEvent.getByLabel(_full5x5E2x2MapLabel, full5x5E2x2Map);
220  iEvent.getByLabel(_full5x5E2x5MaxMapLabel, full5x5E2x5MaxMap);
221  iEvent.getByLabel(_full5x5E5x5MapLabel, full5x5E5x5Map);
222  iEvent.getByLabel(_esEffSigmaRRMapLabel, esEffSigmaRRMap);
223  }
224 
225  // Get the isolation maps
226  iEvent.getByLabel(_phoChargedIsolationLabel, phoChargedIsolationMap);
227  iEvent.getByLabel(_phoPhotonIsolationLabel, phoPhotonIsolationMap);
228  iEvent.getByLabel(_phoWorstChargedIsolationLabel, phoWorstChargedIsolationMap);
229 
230  // Get rho
231  iEvent.getByLabel(_rhoLabel,rho);
232 
233  // Make sure everything is retrieved successfully
234  if(! ( ( !_useValueMaps || ( full5x5SigmaIEtaIEtaMap.isValid()
235  && full5x5SigmaIEtaIPhiMap.isValid()
236  && full5x5E1x3Map.isValid()
237  && full5x5E2x2Map.isValid()
238  && full5x5E2x5MaxMap.isValid()
239  && full5x5E5x5Map.isValid()
240  && esEffSigmaRRMap.isValid() ) )
241  && phoChargedIsolationMap.isValid()
242  && phoPhotonIsolationMap.isValid()
243  && phoWorstChargedIsolationMap.isValid()
244  && rho.isValid() ) )
245  throw cms::Exception("MVA failure: ")
246  << "Failed to retrieve event content needed for this MVA"
247  << std::endl
248  << "Check python MVA configuration file and make sure all needed"
249  << std::endl
250  << "producers are running upstream" << std::endl;
251 
252  // Try to cast the particle into a reco particle.
253  // This should work for both reco and pat.
254  const edm::Ptr<reco::Photon> phoRecoPtr(particle);
255  if( phoRecoPtr.isNull() )
256  throw cms::Exception("MVA failure: ")
257  << " given particle is expected to be reco::Photon or pat::Photon," << std::endl
258  << " but appears to be neither" << std::endl;
259 
260  // Both pat and reco particles have exactly the same accessors.
261  auto superCluster = phoRecoPtr->superCluster();
262  // Full 5x5 cluster shapes. We could take some of this directly from
263  // the photon object, but some of these are not available.
264  float e1x3 = std::numeric_limits<float>::max();
265  float e2x2 = std::numeric_limits<float>::max();
266  float e2x5 = std::numeric_limits<float>::max();
267  float e5x5 = std::numeric_limits<float>::max();
268  float full5x5_sigmaIetaIeta = std::numeric_limits<float>::max();
269  float full5x5_sigmaIetaIphi = std::numeric_limits<float>::max();
270  float effSigmaRR = std::numeric_limits<float>::max();
271 
272  AllVariables allMVAVars;
273 
274  if( _useValueMaps ) { //(before 752)
275  // in principle, in the photon object as well
276  // not in the photon object
277  e1x3 = (*full5x5E1x3Map )[ phoRecoPtr ];
278  e2x2 = (*full5x5E2x2Map )[ phoRecoPtr ];
279  e2x5 = (*full5x5E2x5MaxMap)[ phoRecoPtr ];
280  e5x5 = (*full5x5E5x5Map )[ phoRecoPtr ];
281  full5x5_sigmaIetaIeta = (*full5x5SigmaIEtaIEtaMap)[ phoRecoPtr ];
282  full5x5_sigmaIetaIphi = (*full5x5SigmaIEtaIPhiMap)[ phoRecoPtr ];
283  effSigmaRR = (*esEffSigmaRRMap)[ phoRecoPtr ];
284  } else {
285  // from 753
286  const auto& full5x5_pss = phoRecoPtr->full5x5_showerShapeVariables();
287  e1x3 = full5x5_pss.e1x3;
288  e2x2 = full5x5_pss.e2x2;
289  e2x5 = full5x5_pss.e2x5Max;
290  e5x5 = full5x5_pss.e5x5;
291  full5x5_sigmaIetaIeta = full5x5_pss.sigmaIetaIeta;
292  full5x5_sigmaIetaIphi = full5x5_pss.sigmaIetaIphi;
293  effSigmaRR = full5x5_pss.effSigmaRR;
294  }
295 
296  allMVAVars.varPhi = phoRecoPtr->phi();
297  allMVAVars.varR9 = phoRecoPtr->r9() ;
298  allMVAVars.varSieie = full5x5_sigmaIetaIeta;
299  allMVAVars.varSieip = full5x5_sigmaIetaIphi;
300  allMVAVars.varE1x3overE5x5 = e1x3/e5x5;
301  allMVAVars.varE2x2overE5x5 = e2x2/e5x5;
302  allMVAVars.varE2x5overE5x5 = e2x5/e5x5;
303  allMVAVars.varSCEta = superCluster->eta();
304  allMVAVars.varRawE = superCluster->rawEnergy();
305  allMVAVars.varSCEtaWidth = superCluster->etaWidth();
306  allMVAVars.varSCPhiWidth = superCluster->phiWidth();
307  allMVAVars.varESEnOverRawE = superCluster->preshowerEnergy() / superCluster->rawEnergy();
308  allMVAVars.varESEffSigmaRR = effSigmaRR;
309  allMVAVars.varRho = *rho;
310  allMVAVars.varPhoIsoRaw = (*phoPhotonIsolationMap)[phoRecoPtr];
311  allMVAVars.varChIsoRaw = (*phoChargedIsolationMap)[phoRecoPtr];
312  allMVAVars.varWorstChRaw = (*phoWorstChargedIsolationMap)[phoRecoPtr];
313  // Declare spectator vars
314  allMVAVars.varPt = phoRecoPtr->pt();
315  allMVAVars.varEta = phoRecoPtr->eta();
316 
317  constrainMVAVariables(allMVAVars);
318 
319  std::vector<float> vars;
320  if( isEndcapCategory( findCategory( particle ) ) ) {
321  vars = std::move( packMVAVariables(allMVAVars.varPhi,
322  allMVAVars.varR9,
323  allMVAVars.varSieie,
324  allMVAVars.varSieip,
325  allMVAVars.varE1x3overE5x5,
326  allMVAVars.varE2x2overE5x5,
327  allMVAVars.varE2x5overE5x5,
328  allMVAVars.varSCEta,
329  allMVAVars.varRawE,
330  allMVAVars.varSCEtaWidth,
331  allMVAVars.varSCPhiWidth,
332  allMVAVars.varESEnOverRawE,
333  allMVAVars.varESEffSigmaRR,
334  allMVAVars.varRho,
335  allMVAVars.varPhoIsoRaw,
336  allMVAVars.varChIsoRaw,
337  allMVAVars.varWorstChRaw,
338  // Declare spectator vars
339  allMVAVars.varPt,
340  allMVAVars.varEta)
341  );
342  } else {
343  vars = std::move( packMVAVariables(allMVAVars.varPhi,
344  allMVAVars.varR9,
345  allMVAVars.varSieie,
346  allMVAVars.varSieip,
347  allMVAVars.varE1x3overE5x5,
348  allMVAVars.varE2x2overE5x5,
349  allMVAVars.varE2x5overE5x5,
350  allMVAVars.varSCEta,
351  allMVAVars.varRawE,
352  allMVAVars.varSCEtaWidth,
353  allMVAVars.varSCPhiWidth,
354  allMVAVars.varRho,
355  allMVAVars.varPhoIsoRaw,
356  allMVAVars.varChIsoRaw,
357  allMVAVars.varWorstChRaw,
358  // Declare spectator vars
359  allMVAVars.varPt,
360  allMVAVars.varEta)
361  );
362  }
363  return vars;
364 }
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.full5x5E2x5MaxMap
full5x5E2x5MaxMap
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:76
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.full5x5E2x2Map
full5x5E2x2Map
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:75
AnyMVAEstimatorRun2Base::packMVAVariables
std::vector< float > packMVAVariables(const Args...args) const
Definition: AnyMVAEstimatorRun2Base.h:43
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.full5x5SigmaIEtaIPhiMap
full5x5SigmaIEtaIPhiMap
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:73
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.full5x5E5x5Map
full5x5E5x5Map
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:77
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.full5x5E1x3Map
full5x5E1x3Map
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:74
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:475
heepElectronID_HEEPV50_CSA14_25ns_cff.e2x5
e2x5
Definition: heepElectronID_HEEPV50_CSA14_25ns_cff.py:54
cutBasedElectronID_CSA14_50ns_V1_cff.full5x5SigmaIEtaIEtaMap
full5x5SigmaIEtaIEtaMap
Definition: cutBasedElectronID_CSA14_50ns_V1_cff.py:48
heepElectronID_HEEPV50_CSA14_25ns_cff.e5x5
e5x5
Definition: heepElectronID_HEEPV50_CSA14_25ns_cff.py:55
mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.esEffSigmaRRMap
esEffSigmaRRMap
Definition: mvaPhotonID_PHYS14_PU20bx25_nonTrig_V1_cff.py:78
PhotonMVAEstimatorRun2Spring15NonTrig::findCategory
int findCategory(const edm::Ptr< reco::Candidate > &particle) const override
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.cc:101
eostools.move
def move(src, dest)
Definition: eostools.py:510
int PhotonMVAEstimatorRun2Spring15NonTrig::findCategory ( const edm::Ptr< reco::Candidate > &  particle) const
overridevirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 101 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References funct::abs(), CAT_EB, CAT_EE, PVValHelper::eta, isEndcapCategory(), edm::Ptr< T >::isNull(), reco::Photon::superCluster(), and UNDEFINED.

Referenced by fillMVAVariables(), getTag(), and mvaValue().

101  {
102 
103  // Try to cast the particle into a reco particle.
104  // This should work for both reco and pat.
105  const edm::Ptr<reco::Photon> phoRecoPtr = ( edm::Ptr<reco::Photon> )particle;
106  if( phoRecoPtr.isNull() )
107  throw cms::Exception("MVA failure: ")
108  << " given particle is expected to be reco::Photon or pat::Photon," << std::endl
109  << " but appears to be neither" << std::endl;
110 
111  float eta = phoRecoPtr->superCluster()->eta();
112 
113  //
114  // Determine the category
115  //
116  int iCategory = UNDEFINED;
117  const float ebeeSplit = 1.479; // division between barrel and endcap
118 
119  if ( std::abs(eta) < ebeeSplit)
120  iCategory = CAT_EB;
121 
122  if (std::abs(eta) >= ebeeSplit)
123  iCategory = CAT_EE;
124 
125  return iCategory;
126 }
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const override
Ref to SuperCluster.
edm::Ptr::isNull
bool isNull() const
Checks for null.
Definition: Ptr.h:164
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const std::string& PhotonMVAEstimatorRun2Spring15NonTrig::getName ( void  ) const
inlinefinalvirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 73 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

References _name.

Referenced by plotting.Plot::draw().

73 { return _name; }
int PhotonMVAEstimatorRun2Spring15NonTrig::getNCategories ( ) const
inlineoverridevirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 71 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

References python.rootplot.argparse::category, isEndcapCategory(), and nCategories.

71 { return nCategories; }
const std::string& PhotonMVAEstimatorRun2Spring15NonTrig::getTag ( ) const
inlinefinalvirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 74 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

References _tag, constrainMVAVariables(), fillMVAVariables(), findCategory(), iEvent, and setConsumes().

74 { return _tag; }
bool PhotonMVAEstimatorRun2Spring15NonTrig::isEndcapCategory ( int  category) const

Definition at line 129 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References CAT_EE, createSingleReader(), and GeomDetEnumerators::isEndcap().

Referenced by createSingleReader(), fillMVAVariables(), findCategory(), getNCategories(), and mvaValue().

129  {
130 
131  // For this specific MVA the function is trivial, but kept for possible
132  // future evolution to an MVA with more categories in eta
133  bool isEndcap = false;
134  if( category == CAT_EE )
135  isEndcap = true;
136 
137  return isEndcap;
138 }
GeomDetEnumerators::isEndcap
bool isEndcap(GeomDetEnumerators::SubDetector m)
Definition: GeomDetEnumerators.cc:40
float PhotonMVAEstimatorRun2Spring15NonTrig::mvaValue ( const edm::Ptr< reco::Candidate > &  particle,
const edm::Event &  iEvent 
) const
overridevirtual

Implements AnyMVAEstimatorRun2Base.

Definition at line 59 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References _gbrForests, debug, fillMVAVariables(), findCategory(), isEndcapCategory(), eostools::move(), and mps_fire::result.

Referenced by ~PhotonMVAEstimatorRun2Spring15NonTrig().

59  {
60 
61  const int iCategory = findCategory( particle );
62  const std::vector<float> vars = std::move( fillMVAVariables( particle, iEvent ) );
63 
64  const float result = _gbrForests.at(iCategory)->GetClassifier(vars.data());
65 
66  // DEBUG
67  const bool debug = false;
68  if( debug ){
69  printf("Printout of the photon variable inputs for MVA:\n");
70  printf(" varPhi_ %f\n", vars[0] );
71  printf(" varR9_ %f\n", vars[1] );
72  printf(" varSieie_ %f\n", vars[2] );
73  printf(" varSieip_ %f\n", vars[3] );
74  printf(" varE1x3overE5x5_ %f\n", vars[4] );
75  printf(" varE2x2overE5x5_ %f\n", vars[5] );
76  printf(" varE2x5overE5x5_ %f\n", vars[6] );
77  printf(" varSCEta_ %f\n", vars[7] );
78  printf(" varRawE_ %f\n", vars[8] );
79  printf(" varSCEtaWidth_ %f\n", vars[9] );
80  printf(" varSCPhiWidth_ %f\n", vars[10] );
81  printf(" varRho_ %f\n", vars[11] );
82  printf(" varPhoIsoRaw_ %f\n", vars[12] );
83  printf(" varChIsoRaw_ %f\n", vars[13] );
84  printf(" varWorstChRaw_ %f\n", vars[14] );
85  if( isEndcapCategory( iCategory ) ) {
86  printf(" varESEnOverRawE_ %f\n", vars[15] ); // for endcap MVA only
87  printf(" varESEffSigmaRR_ %f\n", vars[16] ); // for endcap MVA only
88  // The spectators
89  printf(" varPt_ %f\n", vars[17] );
90  printf(" varEta_ %f\n", vars[18] );
91  } else {
92  // The spectators
93  printf(" varPt_ %f\n", vars[15] );
94  printf(" varEta_ %f\n", vars[16] );
95  }
96  }
97 
98  return result;
99 }
mps_fire.result
result
Definition: mps_fire.py:226
debug
#define debug
Definition: HDRShower.cc:19
PhotonMVAEstimatorRun2Spring15NonTrig::findCategory
int findCategory(const edm::Ptr< reco::Candidate > &particle) const override
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.cc:101
PhotonMVAEstimatorRun2Spring15NonTrig::_gbrForests
std::vector< std::unique_ptr< const GBRForest > > _gbrForests
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.h:105
PhotonMVAEstimatorRun2Spring15NonTrig::fillMVAVariables
std::vector< float > fillMVAVariables(const edm::Ptr< reco::Candidate > &particle, const edm::Event &iEvent) const override
Definition: PhotonMVAEstimatorRun2Spring15NonTrig.cc:194
eostools.move
def move(src, dest)
Definition: eostools.py:510
void PhotonMVAEstimatorRun2Spring15NonTrig::setConsumes ( edm::ConsumesCollector &&  cc) const
overridevirtual

Reimplemented from AnyMVAEstimatorRun2Base.

Definition at line 376 of file PhotonMVAEstimatorRun2Spring15NonTrig.cc.

References _esEffSigmaRRMapLabel, _full5x5E1x3MapLabel, _full5x5E2x2MapLabel, _full5x5E2x5MaxMapLabel, _full5x5E5x5MapLabel, _full5x5SigmaIEtaIEtaMapLabel, _full5x5SigmaIEtaIPhiMapLabel, _phoChargedIsolationLabel, _phoPhotonIsolationLabel, _phoWorstChargedIsolationLabel, _rhoLabel, _useValueMaps, and DEFINE_EDM_PLUGIN.

Referenced by getTag().

376  {
377  if( _useValueMaps ) {
378  cc.consumes<edm::ValueMap<float> >(_full5x5SigmaIEtaIEtaMapLabel);
379  cc.consumes<edm::ValueMap<float> >(_full5x5SigmaIEtaIPhiMapLabel);
380  cc.consumes<edm::ValueMap<float> >(_full5x5E1x3MapLabel);
381  cc.consumes<edm::ValueMap<float> >(_full5x5E2x2MapLabel);
382  cc.consumes<edm::ValueMap<float> >(_full5x5E2x5MaxMapLabel);
383  cc.consumes<edm::ValueMap<float> >(_full5x5E5x5MapLabel);
384  cc.consumes<edm::ValueMap<float> >(_esEffSigmaRRMapLabel);
385  }
386  cc.consumes<edm::ValueMap<float> >(_phoChargedIsolationLabel);
387  cc.consumes<edm::ValueMap<float> >(_phoPhotonIsolationLabel);
388  cc.consumes<edm::ValueMap<float> >( _phoWorstChargedIsolationLabel);
389  cc.consumes<double>(_rhoLabel);
390 }
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:52

Member Data Documentation

AllVariables PhotonMVAEstimatorRun2Spring15NonTrig::_allMVAVars
private

Definition at line 109 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by createSingleReader().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_esEffSigmaRRMapLabel
private

Definition at line 125 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5E1x3MapLabel
private

Definition at line 121 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5E2x2MapLabel
private

Definition at line 122 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5E2x5MaxMapLabel
private

Definition at line 123 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5E5x5MapLabel
private

Definition at line 124 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5SigmaIEtaIEtaMapLabel
private

Definition at line 119 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_full5x5SigmaIEtaIPhiMapLabel
private

Definition at line 120 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

std::vector< std::unique_ptr<const GBRForest> > PhotonMVAEstimatorRun2Spring15NonTrig::_gbrForests
private

Definition at line 105 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by mvaValue(), and PhotonMVAEstimatorRun2Spring15NonTrig().

const std::string PhotonMVAEstimatorRun2Spring15NonTrig::_MethodName
private

Definition at line 108 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by createSingleReader().

const std::string PhotonMVAEstimatorRun2Spring15NonTrig::_name = "PhotonMVAEstimatorRun2Spring15NonTrig"
private

Definition at line 97 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by plotting.Subtract::__str__(), plotting.Transform::__str__(), plotting.FakeDuplicate::__str__(), plotting.CutEfficiency::__str__(), plotting.AggregateBins::__str__(), plotting.AggregateHistos::__str__(), plotting.ROC::__str__(), trackingPlots.TimePerEventPlot::__str__(), trackingPlots.TimePerTrackPlot::__str__(), trackingPlots.TimePerEventPlot::_create(), plotting.PlotGroup::_drawSeparate(), plotting.PlotGroup::_save(), plotting.Subtract::create(), plotting.Transform::create(), plotting.FakeDuplicate::create(), plotting.CutEfficiency::create(), plotting.AggregateBins::create(), plotting.AggregateHistos::create(), trackingPlots.TimePerTrackPlot::create(), plotting.Plot::create(), plotting.PlotGroup::draw(), getName(), plotting.Plot::getName(), plotting.PlotGroup::getName(), plotting.PlotterFolder::getName(), plotting.PlotterItem::getName(), plotting.PlotterFolder::getSelectionNameIterator(), validation.SimpleSample::name(), trackingPlots.Iteration::name(), plotting.PlotterItem::readDirs(), and plotting.PlotGroup::remove().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_phoChargedIsolationLabel
private

Definition at line 127 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_phoPhotonIsolationLabel
private

Definition at line 128 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_phoWorstChargedIsolationLabel
private

Definition at line 129 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const edm::InputTag PhotonMVAEstimatorRun2Spring15NonTrig::_rhoLabel
private

Definition at line 130 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

std::string PhotonMVAEstimatorRun2Spring15NonTrig::_tag
private

Definition at line 102 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by getTag(), and PhotonMVAEstimatorRun2Spring15NonTrig().

const bool PhotonMVAEstimatorRun2Spring15NonTrig::_useValueMaps
private

Definition at line 118 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by fillMVAVariables(), and setConsumes().

const int PhotonMVAEstimatorRun2Spring15NonTrig::nCategories = 2

Definition at line 27 of file PhotonMVAEstimatorRun2Spring15NonTrig.h.

Referenced by getNCategories(), and PhotonMVAEstimatorRun2Spring15NonTrig().

Generated for CMSSW Reference Manual by doxygen 1.8.11