HWWAnalyzer Class Reference

#include <HWWAnalyzer.h>

Inheritance diagram for HWWAnalyzer:

Public Member Functions
	HWWAnalyzer (const edm::ParameterSet &)
virtual	~HWWAnalyzer ()
Public Member Functions inherited from thread_unsafe::DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
void	FillHistograms ()

Private Attributes
MonitorElement *	cutflowHist [4]
bool	doTest
EGammaMvaEleEstimator *	egammaMvaEleEstimator
ElectronMaker	electronMaker
EventMaker	eventMaker
EventMonitor	eventMonitor
GSFTrackMaker	gsfTrackMaker
HypDilepMaker	hypDilepMaker
MuonMaker	muonMaker
MuonMVAEstimator *	muonMVAEstimator
MVAJetIdMaker	mvaJetIdMaker
PFCandidateMaker	pfCandidateMaker
PFElectronMaker	pfElectronMaker
PFElToElAssMaker	pfElToElAssMaker
PFJetMaker	pfJetMaker
PFMETMaker	pfMETMaker
RecoConversionMaker	recoConversionMaker
RhoMaker	rhoMaker
TrackMaker	trackMaker
TrkMETMaker	trkMETMaker
VertexMaker	vertexMaker

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 41 of file HWWAnalyzer.h.

Constructor & Destructor Documentation

HWWAnalyzer::HWWAnalyzer

(

const edm::ParameterSet &

iConfig

)

Definition at line 6 of file HWWAnalyzer.cc.

References doTest, egammaMvaEleEstimator, contentValuesFiles::fullPath, edm::ParameterSet::getParameter(), EGammaMvaEleEstimator::initialize(), MuonMVAEstimator::initialize(), MuonMVAEstimator::kIsoRings, EGammaMvaEleEstimator::kTrig, muonMVAEstimator, and AlCaHLTBitMon_QueryRunRegistry::string.

            :vertexMaker        (iConfig, consumesCollector()),
             trackMaker         (iConfig, consumesCollector()),
             electronMaker      (iConfig, consumesCollector()),
             muonMaker          (iConfig, consumesCollector()),
             pfJetMaker         (iConfig, consumesCollector()),
             pfCandidateMaker   (iConfig, consumesCollector()),
             pfElectronMaker    (iConfig, consumesCollector()),
             gsfTrackMaker      (iConfig, consumesCollector()),
             recoConversionMaker(iConfig, consumesCollector()),
             rhoMaker           (iConfig, consumesCollector()),
             pfMETMaker         (iConfig, consumesCollector()),
             mvaJetIdMaker      (iConfig, consumesCollector())
{

  doTest = iConfig.getParameter<bool>("doTest");
  if(doTest) edm::LogInfo("OutputInfo") << "running with doTest==True";

  egammaMvaEleEstimator = 0;
  muonMVAEstimator = 0;


  // --------------- EGamma Id MVA  --------------------------
  std::string egammaweights1 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights1").fullPath();
  std::string egammaweights2 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights2").fullPath();
  std::string egammaweights3 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights3").fullPath();
  std::string egammaweights4 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights4").fullPath();
  std::string egammaweights5 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights5").fullPath();
  std::string egammaweights6 = iConfig.getParameter<edm::FileInPath> ("InputEGammaWeights6").fullPath();
  std::vector<std::string> egammaweights = {
    egammaweights1,
    egammaweights2,
    egammaweights3,
    egammaweights4,
    egammaweights5,
    egammaweights6
  };
  egammaMvaEleEstimator = new EGammaMvaEleEstimator();
  egammaMvaEleEstimator->initialize("BDT", EGammaMvaEleEstimator::kTrig, true, egammaweights );

  // --------------- Muon RingIso MVA  --------------------------
  std::string muonisoweights1 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights1").fullPath();
  std::string muonisoweights2 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights2").fullPath();
  std::string muonisoweights3 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights3").fullPath();
  std::string muonisoweights4 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights4").fullPath();
  std::string muonisoweights5 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights5").fullPath();
  std::string muonisoweights6 = iConfig.getParameter<edm::FileInPath> ("InputMuonIsoWeights6").fullPath();
  std::vector<std::string> muonisoweights = {
    muonisoweights1,
    muonisoweights2,
    muonisoweights3,
    muonisoweights4,
    muonisoweights5,
    muonisoweights6
  };
  muonMVAEstimator = new MuonMVAEstimator();
  muonMVAEstimator->initialize( "MuonIso_BDTG_IsoRings", MuonMVAEstimator::kIsoRings, true, muonisoweights );

}

HWWAnalyzer::~HWWAnalyzer ( )

virtual

Definition at line 67 of file HWWAnalyzer.cc.

                         { 
}

Member Function Documentation

void HWWAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 71 of file HWWAnalyzer.cc.

References HWWFunctions::bestHypothesis(), EventMonitor::hypo_monitor::count(), HWWFunctions::doCutFlow(), doTest, HWWFunctions::EE, egammaMvaEleEstimator, electronMaker, HWWFunctions::EM, eventMaker, eventMonitor, FillHistograms(), gsfTrackMaker, hypDilepMaker, i, HWWFunctions::MM, EventMonitor::monitor, muonMaker, muonMVAEstimator, mvaJetIdMaker, HWWFunctions::passBaseline(), HWWFunctions::passFirstCuts(), pfCandidateMaker, pfElectronMaker, pfElToElAssMaker, pfJetMaker, pfMETMaker, recoConversionMaker, rhoMaker, EventMaker::SetVars(), TrkMETMaker::SetVars(), PFElToElAssMaker::SetVars(), HypDilepMaker::SetVars(), RhoMaker::SetVars(), TrackMaker::SetVars(), GSFTrackMaker::SetVars(), RecoConversionMaker::SetVars(), ElectronMaker::SetVars(), PFMETMaker::SetVars(), VertexMaker::SetVars(), PFJetMaker::SetVars(), PFElectronMaker::SetVars(), MVAJetIdMaker::SetVars(), PFCandidateMaker::SetVars(), MuonMaker::SetVars(), trackMaker, trkMETMaker, and vertexMaker.

                                                                            {

  using namespace HWWFunctions;

  HWW hww;

  //count total events
  eventMonitor.monitor.count(MM, "total events", 1.0);
  eventMonitor.monitor.count(EE, "total events", 1.0);
  eventMonitor.monitor.count(EM, "total events", 1.0);
  eventMonitor.monitor.count(ME, "total events", 1.0);

  //if doTest flag is true, all we do is access all the collections 
  //without having to make it through the cutflow.
  if(doTest){

    eventMaker          .SetVars(hww, iEvent, iSetup);
    vertexMaker         .SetVars(hww, iEvent, iSetup);
    trackMaker          .SetVars(hww, iEvent, iSetup);
    electronMaker       .SetVars(hww, iEvent, iSetup);
    muonMaker           .SetVars(hww, iEvent, iSetup);
    pfJetMaker          .SetVars(hww, iEvent, iSetup);
    hypDilepMaker       .SetVars(hww, iEvent, iSetup);
    pfCandidateMaker    .SetVars(hww, iEvent, iSetup);
    pfElectronMaker     .SetVars(hww, iEvent, iSetup);
    pfElToElAssMaker    .SetVars(hww, iEvent, iSetup);
    gsfTrackMaker       .SetVars(hww, iEvent, iSetup);
    recoConversionMaker .SetVars(hww, iEvent, iSetup);
    rhoMaker            .SetVars(hww, iEvent, iSetup);
    pfMETMaker          .SetVars(hww, iEvent, iSetup);
    trkMETMaker         .SetVars(hww, iEvent, iSetup);
    mvaJetIdMaker       .SetVars(hww, iEvent, iSetup);

    return;

  }

  //get variables
  eventMaker    .SetVars(hww, iEvent, iSetup);
  vertexMaker   .SetVars(hww, iEvent, iSetup);
  trackMaker    .SetVars(hww, iEvent, iSetup);

  if(hww.trks_trk_p4().size() < 2) return;

  electronMaker .SetVars(hww, iEvent, iSetup);
  muonMaker     .SetVars(hww, iEvent, iSetup);
  pfJetMaker    .SetVars(hww, iEvent, iSetup);
  hypDilepMaker .SetVars(hww, iEvent, iSetup);

  //check some basic event requirements
  std::vector<int> goodHyps;
  for(unsigned int i=0; i < hww.hyp_p4().size(); i++){
    if(!passFirstCuts(hww, i)) continue;
    goodHyps.push_back(i);
  }
  
  //no need to continue if event didn't pass basic requirements
  if(goodHyps.size() > 0){

    //get variables
    pfCandidateMaker    .SetVars(hww, iEvent, iSetup);
    pfElectronMaker     .SetVars(hww, iEvent, iSetup);
    pfElToElAssMaker    .SetVars(hww, iEvent, iSetup);
    gsfTrackMaker       .SetVars(hww, iEvent, iSetup);
    recoConversionMaker .SetVars(hww, iEvent, iSetup);
    rhoMaker            .SetVars(hww, iEvent, iSetup);

    //to hold indices of candidate lepton pairs
    std::vector<int> candidates;

    //get lepton pairs that pass baseline selection
    for(unsigned int i=0; i < goodHyps.size(); i++){
      if(!passBaseline(hww, goodHyps.at(i), egammaMvaEleEstimator, muonMVAEstimator)) continue;
      candidates.push_back(i);      
    }

    if(candidates.size()>0){

      //get variables
      pfMETMaker            .SetVars(hww, iEvent, iSetup);
      trkMETMaker           .SetVars(hww, iEvent, iSetup);
      mvaJetIdMaker         .SetVars(hww, iEvent, iSetup);

      //find best lepton pair
      int bestHyp = bestHypothesis(hww, candidates);

      //perform remaining selections
      doCutFlow(hww, bestHyp, eventMonitor.monitor, egammaMvaEleEstimator, muonMVAEstimator);

    }
  }  

  FillHistograms();

}//end analyze

void HWWAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements thread_unsafe::DQMEDAnalyzer.

Definition at line 168 of file HWWAnalyzer.cc.

References DQMStore::IBooker::book1D(), cutflowHist, and DQMStore::IBooker::setCurrentFolder().

                                                                                               {

  ibooker.setCurrentFolder("PhysicsHWW");

  cutflowHist[0] = ibooker.book1D("cutflow_mm", "HWW cutflow mm", 21, 0, 21);   
  cutflowHist[1] = ibooker.book1D("cutflow_ee", "HWW cutflow ee", 21, 0, 21);   
  cutflowHist[2] = ibooker.book1D("cutflow_em", "HWW cutflow em", 21, 0, 21);   
  cutflowHist[3] = ibooker.book1D("cutflow_me", "HWW cutflow me", 21, 0, 21);   
  
}

void HWWAnalyzer::FillHistograms ( )

private

Definition at line 179 of file HWWAnalyzer.cc.

References EventMonitor::hypo_monitor::counters, cutflowHist, eventMonitor, i, j, EventMonitor::monitor, MonitorElement::setBinContent(), and MonitorElement::setBinLabel().

Referenced by analyze().

                                {

  for (unsigned int i=0; i<4; i++){
    for (unsigned int j=0; j<eventMonitor.monitor.counters.size(); j++){
      cutflowHist[i]->setBinContent(j+1, eventMonitor.monitor.counters[j].nevt[i] - 1); //the "- 1" accounts for the initial count in the EventMonitor constructor
      cutflowHist[i]->setBinLabel(j+1, eventMonitor.monitor.counters[j].name.c_str(), 1);
    }
  }

}

Member Data Documentation

MonitorElement* HWWAnalyzer::cutflowHist[4]

private

Definition at line 79 of file HWWAnalyzer.h.

Referenced by bookHistograms(), and FillHistograms().

bool HWWAnalyzer::doTest

private

Definition at line 55 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

EGammaMvaEleEstimator* HWWAnalyzer::egammaMvaEleEstimator

private

Definition at line 57 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

ElectronMaker HWWAnalyzer::electronMaker

private

Definition at line 65 of file HWWAnalyzer.h.

Referenced by analyze().

EventMaker HWWAnalyzer::eventMaker

private

Definition at line 63 of file HWWAnalyzer.h.

Referenced by analyze().

EventMonitor HWWAnalyzer::eventMonitor

private

Definition at line 60 of file HWWAnalyzer.h.

Referenced by analyze(), and FillHistograms().

GSFTrackMaker HWWAnalyzer::gsfTrackMaker

private

Definition at line 72 of file HWWAnalyzer.h.

Referenced by analyze().

HypDilepMaker HWWAnalyzer::hypDilepMaker

private

Definition at line 68 of file HWWAnalyzer.h.

Referenced by analyze().

MuonMaker HWWAnalyzer::muonMaker

private

Definition at line 66 of file HWWAnalyzer.h.

Referenced by analyze().

MuonMVAEstimator* HWWAnalyzer::muonMVAEstimator

private

Definition at line 58 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

MVAJetIdMaker HWWAnalyzer::mvaJetIdMaker

private

Definition at line 77 of file HWWAnalyzer.h.

Referenced by analyze().

PFCandidateMaker HWWAnalyzer::pfCandidateMaker

private

Definition at line 69 of file HWWAnalyzer.h.

Referenced by analyze().

PFElectronMaker HWWAnalyzer::pfElectronMaker

private

Definition at line 70 of file HWWAnalyzer.h.

Referenced by analyze().

PFElToElAssMaker HWWAnalyzer::pfElToElAssMaker

private

Definition at line 71 of file HWWAnalyzer.h.

Referenced by analyze().

PFJetMaker HWWAnalyzer::pfJetMaker

private

Definition at line 67 of file HWWAnalyzer.h.

Referenced by analyze().

PFMETMaker HWWAnalyzer::pfMETMaker

private

Definition at line 75 of file HWWAnalyzer.h.

Referenced by analyze().

RecoConversionMaker HWWAnalyzer::recoConversionMaker

private

Definition at line 73 of file HWWAnalyzer.h.

Referenced by analyze().

RhoMaker HWWAnalyzer::rhoMaker

private

Definition at line 74 of file HWWAnalyzer.h.

Referenced by analyze().

TrackMaker HWWAnalyzer::trackMaker

private

Definition at line 64 of file HWWAnalyzer.h.

Referenced by analyze().

TrkMETMaker HWWAnalyzer::trkMETMaker

private

Definition at line 76 of file HWWAnalyzer.h.

Referenced by analyze().

VertexMaker HWWAnalyzer::vertexMaker

private

Definition at line 62 of file HWWAnalyzer.h.

Referenced by analyze().