HWWAnalyzer Class Reference

#include <HWWAnalyzer.h>

Inheritance diagram for HWWAnalyzer:

Public Member Functions
	HWWAnalyzer (const edm::ParameterSet &)
virtual	~HWWAnalyzer ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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 Attributes
bool	doTest
EGammaMvaEleEstimator *	egammaMvaEleEstimator
ElectronMaker	electronMaker
EventMaker	eventMaker
std::unique_ptr< 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::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 42 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::stream::EDAnalyzerBase.

Definition at line 71 of file HWWAnalyzer.cc.

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

                                                                            {

  using namespace HWWFunctions;

  HWW hww;

  //count total events
  eventMonitor->count(MM, "total events", 1.0);
  eventMonitor->count(EE, "total events", 1.0);
  eventMonitor->count(EM, "total events", 1.0);
  eventMonitor->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, egammaMvaEleEstimator, muonMVAEstimator);

    }
  }  
}//end analyze

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 165 of file HWWAnalyzer.cc.

References eventMonitor.

                                                                                                {

  eventMonitor.reset(new EventMonitor(ibooker));
}

Member Data Documentation

bool HWWAnalyzer::doTest

private

Definition at line 56 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

EGammaMvaEleEstimator* HWWAnalyzer::egammaMvaEleEstimator

private

Definition at line 58 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

ElectronMaker HWWAnalyzer::electronMaker

private

Definition at line 64 of file HWWAnalyzer.h.

Referenced by analyze().

EventMaker HWWAnalyzer::eventMaker

private

Definition at line 62 of file HWWAnalyzer.h.

Referenced by analyze().

std::unique_ptr<EventMonitor> HWWAnalyzer::eventMonitor

private

Definition at line 78 of file HWWAnalyzer.h.

Referenced by analyze(), and bookHistograms().

GSFTrackMaker HWWAnalyzer::gsfTrackMaker

private

Definition at line 71 of file HWWAnalyzer.h.

Referenced by analyze().

HypDilepMaker HWWAnalyzer::hypDilepMaker

private

Definition at line 67 of file HWWAnalyzer.h.

Referenced by analyze().

MuonMaker HWWAnalyzer::muonMaker

private

Definition at line 65 of file HWWAnalyzer.h.

Referenced by analyze().

MuonMVAEstimator* HWWAnalyzer::muonMVAEstimator

private

Definition at line 59 of file HWWAnalyzer.h.

Referenced by analyze(), and HWWAnalyzer().

MVAJetIdMaker HWWAnalyzer::mvaJetIdMaker

private

Definition at line 76 of file HWWAnalyzer.h.

Referenced by analyze().

PFCandidateMaker HWWAnalyzer::pfCandidateMaker

private

Definition at line 68 of file HWWAnalyzer.h.

Referenced by analyze().

PFElectronMaker HWWAnalyzer::pfElectronMaker

private

Definition at line 69 of file HWWAnalyzer.h.

Referenced by analyze().

PFElToElAssMaker HWWAnalyzer::pfElToElAssMaker

private

Definition at line 70 of file HWWAnalyzer.h.

Referenced by analyze().

PFJetMaker HWWAnalyzer::pfJetMaker

private

Definition at line 66 of file HWWAnalyzer.h.

Referenced by analyze().

PFMETMaker HWWAnalyzer::pfMETMaker

private

Definition at line 74 of file HWWAnalyzer.h.

Referenced by analyze().

RecoConversionMaker HWWAnalyzer::recoConversionMaker

private

Definition at line 72 of file HWWAnalyzer.h.

Referenced by analyze().

RhoMaker HWWAnalyzer::rhoMaker

private

Definition at line 73 of file HWWAnalyzer.h.

Referenced by analyze().

TrackMaker HWWAnalyzer::trackMaker

private

Definition at line 63 of file HWWAnalyzer.h.

Referenced by analyze().

TrkMETMaker HWWAnalyzer::trkMETMaker

private

Definition at line 75 of file HWWAnalyzer.h.

Referenced by analyze().

VertexMaker HWWAnalyzer::vertexMaker

private

Definition at line 61 of file HWWAnalyzer.h.

Referenced by analyze().