#include <RazorVarAnalyzer.h>

Inheritance diagram for RazorVarAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	RazorVarAnalyzer (const edm::ParameterSet &)

	~RazorVarAnalyzer () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::InputTag	m_eleCollectionTag

edm::EDGetTokenT< reco::ElectronCollection >	m_eleCollectionTagToken_

edm::InputTag	m_jetCollectionTag

edm::EDGetTokenT< reco::CaloJetCollection >	m_jetCollectionTagToken_

edm::InputTag	m_muCollectionTag

edm::EDGetTokenT< std::vector< reco::RecoChargedCandidate > >	m_muCollectionTagToken_

edm::InputTag	m_razorVarCollectionTag

edm::EDGetTokenT< std::vector< double > >	m_razorVarCollectionTagToken_

MonitorElement *	m_rsqMREle

MonitorElement *	m_rsqMREleEle

MonitorElement *	m_rsqMREleMJ

MonitorElement *	m_rsqMREleMu

MonitorElement *	m_rsqMRFullyInc

MonitorElement *	m_rsqMRHad

MonitorElement *	m_rsqMRHadMJ

MonitorElement *	m_rsqMRInc10J

MonitorElement *	m_rsqMRInc12J

MonitorElement *	m_rsqMRInc14J

MonitorElement *	m_rsqMRInc4J

MonitorElement *	m_rsqMRInc6J

MonitorElement *	m_rsqMRInc8J

MonitorElement *	m_rsqMRMu

MonitorElement *	m_rsqMRMuMJ

MonitorElement *	m_rsqMRMuMu

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from ScoutingAnalyzerBase
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

MonitorElement *	bookH1 (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")

MonitorElement *	bookH1andDivide (DQMStore::IBooker &, const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH1BinArray (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, float xbinsize, const std::string &titleX="", const std::string &titleY="Events", Option_t option="E1 P")

MonitorElement *	bookH1withSumw2 (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")

MonitorElement *	bookH1withSumw2BinArray (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, float xbinsize, const std::string &titleX="", const std::string &titleY="Events", Option_t option="E1 P")

MonitorElement *	bookH2 (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")

MonitorElement *	bookH2andDivide (DQMStore::IBooker &, const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH2withSumw2 (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")

MonitorElement *	bookP1 (DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="E1 P")

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

std::string	newName (const std::string &name)

void	prepareBooking (DQMStore::IBooker &)

MonitorElement *	profileX (DQMStore::IBooker &, MonitorElement *me2d, const std::string &title="", const std::string &titleX="", const std::string &titleY="", Double_t minimum=-1111, Double_t maximum=-1111)

MonitorElement *	profileY (DQMStore::IBooker &, MonitorElement *me2d, const std::string &title="", const std::string &titleX="", const std::string &titleY="", Double_t minimum=-1111, Double_t maximum=-1111)

	ScoutingAnalyzerBase (const edm::ParameterSet &conf)

	~ScoutingAnalyzerBase () 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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 10 of file RazorVarAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file RazorVarAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_eleCollectionTagToken_, m_jetCollectionTagToken_, m_muCollectionTagToken_, and m_razorVarCollectionTagToken_.

                                                                 :
   ScoutingAnalyzerBase(conf),
   m_eleCollectionTag(conf.getUntrackedParameter<edm::InputTag>("eleCollectionName",edm::InputTag("hltPixelMatchElectronsActivity"))),
   m_jetCollectionTag(conf.getUntrackedParameter<edm::InputTag>("jetCollectionName",edm::InputTag("hltCaloJetIDPassed"))),
   m_muCollectionTag(conf.getUntrackedParameter<edm::InputTag>("muCollectionName",edm::InputTag("hltL3MuonCandidates"))),
   m_razorVarCollectionTag(conf.getUntrackedParameter<edm::InputTag>("razorVarCollectionName")) {
     //set Token(-s)
     m_jetCollectionTagToken_ = consumes<reco::CaloJetCollection>(conf.getUntrackedParameter<edm::InputTag>("jetCollectionName",edm::InputTag("hltCaloJetIDPassed")));
     m_muCollectionTagToken_ = consumes<std::vector<reco::RecoChargedCandidate> >(conf.getUntrackedParameter<edm::InputTag>("muCollectionName",edm::InputTag("hltL3MuonCandidates")));
     m_eleCollectionTagToken_ = consumes<reco::ElectronCollection>(conf.getUntrackedParameter<edm::InputTag>("eleCollectionName",edm::InputTag("hltPixelMatchElectronsActivity")));
     m_razorVarCollectionTagToken_ = consumes<std::vector<double> >(conf.getUntrackedParameter<edm::InputTag>("razorVarCollectionName"));
 }

RazorVarAnalyzer::~RazorVarAnalyzer ( )

override

Definition at line 28 of file RazorVarAnalyzer.cc.

28 {}

Member Function Documentation

void RazorVarAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	c
	)

override

Definition at line 169 of file RazorVarAnalyzer.cc.

References MonitorElement::Fill(), edm::Event::getByToken(), edm::HandleBase::isValid(), m_eleCollectionTagToken_, m_jetCollectionTagToken_, m_muCollectionTagToken_, m_razorVarCollectionTagToken_, m_rsqMREle, m_rsqMREleEle, m_rsqMREleMJ, m_rsqMREleMu, m_rsqMRFullyInc, m_rsqMRHad, m_rsqMRHadMJ, m_rsqMRInc10J, m_rsqMRInc12J, m_rsqMRInc14J, m_rsqMRInc4J, m_rsqMRInc6J, m_rsqMRInc8J, m_rsqMRMu, m_rsqMRMuMJ, m_rsqMRMuMu, and dttmaxenums::R.

                                                                                  {
   //count the number of jets with a minimal selection
   edm::Handle<reco::CaloJetCollection> calojets_handle;
   iEvent.getByToken(m_jetCollectionTagToken_, calojets_handle);
 
   unsigned int njets = 0;
   for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it){
     if(it->pt() >= 30. && fabs(it->eta()) <= 3.0){
       njets++;
     }
   }
 
   //count the number of muons
   edm::Handle<std::vector<reco::RecoChargedCandidate> > muon_handle;
   iEvent.getByToken(m_muCollectionTagToken_, muon_handle);
 
   unsigned int nmu_loose = 0;
   unsigned int nmu_tight = 0;
   if(muon_handle.isValid()){
     for(std::vector<reco::RecoChargedCandidate>::const_iterator it = muon_handle->begin(); it != muon_handle->end(); ++it){
       if(it->pt() >= 15 && fabs(it->eta()) <= 2.1) nmu_tight++;
       if(it->pt() >= 10 && fabs(it->eta()) <= 2.4) nmu_loose++;
     }
   }
 
   //count the number of electrons
   edm::Handle<reco::ElectronCollection> ele_handle;
   iEvent.getByToken(m_eleCollectionTagToken_, ele_handle);
 
   unsigned int nele_loose = 0;
   unsigned int nele_tight = 0;
   if(ele_handle.isValid()){
     for(reco::ElectronCollection::const_iterator it = ele_handle->begin(); it != ele_handle->end(); ++it){
       if(it->pt() >= 20 && fabs(it->eta()) <= 2.5) nele_tight++;
       if(it->pt() >= 10 && fabs(it->eta()) <= 2.5) nele_loose++;
     }
   }
 
   //now get the box number: {'MuEle':0,'MuMu':1,'EleEle':2,'Mu':3,'Ele':4,'Had':5}
   unsigned int box_num = 5;
   if(nmu_tight > 0 && nele_tight > 0){
     box_num = 0;
   }else if(nmu_tight > 0 && nmu_loose > 1){
     box_num = 1;
   }else if(nele_tight > 0 && nele_loose > 1){
     box_num = 2;
   }else if(nmu_tight > 0){
     box_num = 3;
   }else if(nele_tight > 0){
     box_num = 4;
   }
 
   edm::Handle<std::vector<double> > razorvar_handle;
   iEvent.getByToken(m_razorVarCollectionTagToken_, razorvar_handle);
   if(razorvar_handle->size() > 1){
     const double MR = razorvar_handle->at(0);
     const double R = razorvar_handle->at(1);
     m_rsqMRFullyInc->Fill(MR,R*R);
     if(njets >= 4) m_rsqMRInc4J->Fill(MR,R*R);
     if(njets >= 6) m_rsqMRInc6J->Fill(MR,R*R);
     if(njets >= 8) m_rsqMRInc8J->Fill(MR,R*R);
     if(njets >= 10) m_rsqMRInc10J->Fill(MR,R*R);
     if(njets >= 12) m_rsqMRInc12J->Fill(MR,R*R);
     if(njets >= 14) m_rsqMRInc14J->Fill(MR,R*R);
 
     //now fill the boxes
     if(box_num == 0) m_rsqMREleMu->Fill(MR,R*R);
     if(box_num == 1) m_rsqMRMuMu->Fill(MR,R*R);
     if(box_num == 2) m_rsqMREleEle->Fill(MR,R*R);
     if(box_num == 3) m_rsqMRMu->Fill(MR,R*R);
     if(box_num == 4) m_rsqMREle->Fill(MR,R*R);
     if(box_num == 5) m_rsqMRHad->Fill(MR,R*R);
 
     //finally the multijet boxes - think ttbar
     //muon boxes: muons are not in jets
     if( box_num == 3 && njets >= 4) m_rsqMRMuMJ->Fill(MR,R*R);
     //ele boxes: electrons are in jets
     else if( box_num == 4 && njets >= 5) m_rsqMREleMJ->Fill(MR,R*R);
     //fill the Had box
     else if( box_num == 5 && njets >= 6) m_rsqMRHadMJ->Fill(MR,R*R);
   }
 } 

void RazorVarAnalyzer::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 31 of file RazorVarAnalyzer.cc.

References ScoutingAnalyzerBase::bookH2withSumw2(), m_rsqMREle, m_rsqMREleEle, m_rsqMREleMJ, m_rsqMREleMu, m_rsqMRFullyInc, m_rsqMRHad, m_rsqMRHadMJ, m_rsqMRInc10J, m_rsqMRInc12J, m_rsqMRInc14J, m_rsqMRInc4J, m_rsqMRInc6J, m_rsqMRInc8J, m_rsqMRMu, m_rsqMRMuMJ, m_rsqMRMuMu, and ScoutingAnalyzerBase::prepareBooking().

                                                                                                       {
   ScoutingAnalyzerBase::prepareBooking(iBooker);
   //the full inclusive histograms
   m_rsqMRFullyInc = bookH2withSumw2(
       iBooker,
       "rsqMRFullyInc",
       "M_{R} vs R^{2} (All Events)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc4J = bookH2withSumw2(
       iBooker,
       "rsqMRInc4J",
       "M_{R} vs R^{2} (>= 4j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc6J = bookH2withSumw2(
       iBooker,
       "rsqMRInc6J",
       "M_{R} vs R^{2} (>= 6j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc8J = bookH2withSumw2(
       iBooker,
       "rsqMRInc8J",
       "M_{R} vs R^{2} (>= 8j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc10J = bookH2withSumw2(
       iBooker,
       "rsqMRInc10J",
       "M_{R} vs R^{2} (>= 10j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc12J = bookH2withSumw2(
       iBooker,
       "rsqMRInc12J",
       "M_{R} vs R^{2} (>= 12j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRInc14J = bookH2withSumw2(
       iBooker,
       "rsqMRInc14J",
       "M_{R} vs R^{2} (>= 14j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
 
   //the by box histograms
   m_rsqMREleMu = bookH2withSumw2(
       iBooker,
       "rsqMREleMu",
       "M_{R} vs R^{2} (EleMu box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRMuMu = bookH2withSumw2(
       iBooker,
       "rsqMRMuMu",
       "M_{R} vs R^{2} (MuMu box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMREleEle = bookH2withSumw2(
       iBooker,
       "rsqMREleEle",
       "M_{R} vs R^{2} (EleEle box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRMu = bookH2withSumw2(
       iBooker,
       "rsqMRMu",
       "M_{R} vs R^{2} (Mu box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMREle = bookH2withSumw2(
       iBooker,
       "rsqMREle",
       "M_{R} vs R^{2} (Ele box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRHad = bookH2withSumw2(
       iBooker,
       "rsqMRHad",
       "M_{R} vs R^{2} (Had box)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
 
   //the by box histograms
   m_rsqMRMuMJ = bookH2withSumw2(
       iBooker,
       "rsqMRMuMJ",
       "M_{R} vs R^{2} (Mu box >= 4j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMREleMJ = bookH2withSumw2(
       iBooker,
       "rsqMREleMJ",
       "M_{R} vs R^{2} (Ele box >= 5j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
   m_rsqMRHadMJ = bookH2withSumw2(
       iBooker,
       "rsqMRHadMJ",
       "M_{R} vs R^{2} (Had box >= 6j)",
       400,0.,4000.,
       50,0.,1.,
       "M_{R} [GeV]",
       "R^{2}");
 }