#include <DiJetMonitor.h>

Inheritance diagram for DiJetMonitor:

Public Member Functions
	DiJetMonitor (const edm::ParameterSet &)

	~DiJetMonitor () override throw ()

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

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)

Public Member Functions inherited from TriggerDQMBase
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbins, double xmin, double xmax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbinsX, double xmin, double xmax, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbinsX, double xmin, double xmax, unsigned nbinsY, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX, const std::vector< double > &binningY)

void	setMETitle (ObjME &me, const std::string &titleX, const std::string &titleY)

	TriggerDQMBase ()=default

virtual	~TriggerDQMBase ()=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from TriggerDQMBase
static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)

static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)

static MEbinning	getHistoLSPSet (const edm::ParameterSet &pset)

static MEbinning	getHistoPSet (const edm::ParameterSet &pset)

Protected Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override

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

bool	dijet_selection (double eta_1, double phi_1, double eta_2, double phi_2, double pt_1, double pt_2, int &tag_id, int &probe_id, int Event)

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
MEbinning	asy_binning

std::unique_ptr< GenericTriggerEventFlag >	den_genTriggerEventFlag_

MEbinning	dijet_eta_binning

const double	DiJet_MAX_ETA = 5

const double	DiJet_MAX_PHI = 3.2

unsigned int	DiJet_N_ETA = 20

unsigned int	DiJet_N_PHI = 32

MEbinning	dijet_phi_binning

MEbinning	dijetpt_binning_

MEbinning	dijetptThr_binning_

edm::EDGetTokenT< reco::PFJetCollection >	dijetSrc_

edm::EDGetTokenT< reco::GsfElectronCollection >	eleToken_

std::string	folderName_

std::string	histoSuffix_

ObjME	jetAsyEtaME_

ObjME	jetEtaPhiME_

ObjME	jetetaPrbME_

ObjME	jetetaTagME_

ObjME	jetphiPrbME_

ObjME	jetpt1ME_

ObjME	jetpt2ME_

ObjME	jetptAsyME_

ObjME	jetptAvgaME_

ObjME	jetptAvgaThrME_

ObjME	jetptAvgbME_

ObjME	jetptPrbME_

ObjME	jetptTagME_

const double	MAX_asy = 1

edm::EDGetTokenT< reco::PFMETCollection >	metToken_

const double	MIN_asy = -1

edm::EDGetTokenT< reco::MuonCollection >	muoToken_

unsigned int	N_asy = 50

int	nmuons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

double	ptcut_

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 49 of file DiJetMonitor.h.

Constructor & Destructor Documentation

DiJetMonitor::DiJetMonitor ( const edm::ParameterSet & iConfig )

Definition at line 14 of file DiJetMonitor.cc.

References dijetpt_binning_, dijetptThr_binning_, dijetSrc_, folderName_, TriggerDQMBase::getHistoPSet(), edm::ParameterSet::getParameter(), ptcut_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                           :
 num_genTriggerEventFlag_ ( new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"),consumesCollector(), *this))
 ,den_genTriggerEventFlag_ ( new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"),consumesCollector(), *this))
 {
   folderName_              = iConfig.getParameter<std::string>("FolderName"); 
   dijetSrc_                = mayConsume<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("dijetSrc"));//jet 
   
   dijetpt_binning_         = getHistoPSet   (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>   ("dijetPSet")    );
   dijetptThr_binning_      = getHistoPSet   (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>   ("dijetPtThrPSet")    );
   
   ptcut_      = iConfig.getParameter<double>("ptcut" ); 
 
 }

DiJetMonitor::~DiJetMonitor	(	)
throw	(
	)

inlineoverride

Definition at line 53 of file DiJetMonitor.h.

References analyze(), bookHistograms(), dijet_selection(), fillDescriptions(), and iEvent.

53 {};

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 101 of file DiJetMonitor.cc.

References den_genTriggerEventFlag_, TriggerDQMBase::ObjME::denominator, dijet_selection(), dijetSrc_, edm::EventID::event(), MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), jetAsyEtaME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetphiPrbME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, num_genTriggerEventFlag_, TriggerDQMBase::ObjME::numerator, and ptcut_.

Referenced by ~DiJetMonitor().

                                                                                {
   int Event = -999;
   Event = iEvent.id().event();
   // Filter out events if Trigger Filtering is requested
   if (den_genTriggerEventFlag_->on() && ! den_genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
  
   std::vector<double> v_jetpt;
   std::vector<double> v_jeteta;
   std::vector<double> v_jetphi;
     
   v_jetpt.clear();
   v_jeteta.clear();
   v_jetphi.clear();
 
   edm::Handle<reco::PFJetCollection > offjets;
   iEvent.getByToken( dijetSrc_, offjets );
   if (!offjets.isValid()){
       edm::LogWarning("DiJetMonitor") << "DiJet handle not valid \n";
       return;
   }
   for ( reco::PFJetCollection::const_iterator ibegin = offjets->begin(), iend = offjets->end(), ijet = ibegin; ijet != iend; ++ijet ) {
     if (ijet->pt()< ptcut_) {continue;}
     v_jetpt.push_back(ijet->pt()); 
     v_jeteta.push_back(ijet->eta()); 
     v_jetphi.push_back(ijet->phi()); 
   }
   if (v_jetpt.size() < 2) {return;}
   double pt_1 = v_jetpt[0];
   double eta_1 = v_jeteta[0];
   double phi_1 = v_jetphi[0];
   double pt_2 = v_jetpt[1];
   double eta_2 = v_jeteta[1];
   double phi_2 = v_jetphi[1];
   double pt_avg_b = (pt_1 + pt_2)*0.5;
   int tag_id = -999, probe_id = -999;
 
   jetpt1ME_.denominator    -> Fill(pt_1);
   jetpt2ME_.denominator    -> Fill(pt_2);
   jetptAvgbME_.denominator -> Fill(pt_avg_b);
 
   if( dijet_selection(eta_1, phi_1, eta_2, phi_2, pt_1, pt_2,tag_id, probe_id,Event) ){ 
   
     if(tag_id == 0 && probe_id == 1) {
       double pt_asy = (pt_2 - pt_1)/(pt_1 + pt_2);
       double pt_avg = (pt_1 + pt_2)*0.5;
       jetptAvgaME_.denominator -> Fill(pt_avg);
       jetptAvgaThrME_.denominator -> Fill(pt_avg);
       jetptTagME_.denominator  -> Fill(pt_1);
       jetptPrbME_.denominator  -> Fill(pt_2);
       jetetaPrbME_.denominator -> Fill(eta_2);
       jetetaTagME_.denominator -> Fill(eta_1);
       jetptAsyME_.denominator  ->Fill(pt_asy);
       jetphiPrbME_.denominator -> Fill(phi_2);
       jetAsyEtaME_.denominator -> Fill(pt_asy,eta_2);
       jetEtaPhiME_.denominator -> Fill(eta_2,phi_2);
     }
     if(tag_id == 1 && probe_id == 0) {
       double pt_asy = (pt_1 - pt_2)/(pt_2 + pt_1);
       double pt_avg = (pt_2 + pt_1)*0.5;
       jetptAvgaME_.denominator -> Fill(pt_avg);
       jetptAvgaThrME_.denominator -> Fill(pt_avg);
       jetptTagME_.denominator  -> Fill(pt_2);
       jetptPrbME_.denominator  -> Fill(pt_1);
       jetetaPrbME_.denominator -> Fill(eta_1);
       jetetaTagME_.denominator -> Fill(eta_2);
       jetptAsyME_.denominator  ->Fill(pt_asy);
       jetphiPrbME_.denominator -> Fill(phi_1);
       jetAsyEtaME_.denominator -> Fill(pt_asy,eta_1);
       jetEtaPhiME_.denominator -> Fill(eta_1,phi_1);
     }
 
   if (num_genTriggerEventFlag_->on() && ! num_genTriggerEventFlag_->accept( iEvent, iSetup) ) return; 
 
      jetpt1ME_.numerator   -> Fill(pt_1);
      jetpt2ME_.numerator   -> Fill(pt_2);
      jetptAvgbME_.numerator -> Fill(pt_avg_b);
 
     if(tag_id == 0 && probe_id == 1) {
       double pt_asy = (pt_2 - pt_1)/(pt_1 + pt_2);
       double pt_avg = (pt_1 + pt_2)*0.5;
       jetptAvgaME_.numerator -> Fill(pt_avg);
       jetptAvgaThrME_.numerator -> Fill(pt_avg);
       jetptTagME_.numerator  -> Fill(pt_1);
       jetptPrbME_.numerator  -> Fill(pt_2);
       jetetaPrbME_.numerator -> Fill(eta_2);
       jetetaTagME_.numerator -> Fill(eta_1);
       jetptAsyME_.numerator  ->Fill(pt_asy);
       jetphiPrbME_.numerator -> Fill(phi_2);
       jetAsyEtaME_.numerator -> Fill(pt_asy,eta_2);
       jetEtaPhiME_.numerator -> Fill(eta_2,phi_2);
     }
     if(tag_id == 1 && probe_id == 0) {
       double pt_asy = (pt_1 - pt_2)/(pt_2 + pt_1);
       double pt_avg = (pt_2 + pt_1)*0.5;
       jetptAvgaME_.numerator -> Fill(pt_avg);
       jetptAvgaThrME_.numerator -> Fill(pt_avg);
       jetptTagME_.numerator  -> Fill(pt_2);
       jetptPrbME_.numerator  -> Fill(pt_1);
       jetetaPrbME_.numerator -> Fill(eta_1);
       jetetaTagME_.numerator -> Fill(eta_2);
       jetptAsyME_.numerator  ->Fill(pt_asy);
       jetphiPrbME_.numerator -> Fill(phi_1);
       jetAsyEtaME_.numerator -> Fill(pt_asy,eta_1);
       jetEtaPhiME_.numerator -> Fill(eta_1,phi_1);
     }
   }
 }

void DiJetMonitor::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 28 of file DiJetMonitor.cc.

References asy_binning, TriggerDQMBase::bookME(), den_genTriggerEventFlag_, dijet_eta_binning, dijet_phi_binning, dijetpt_binning_, dijetptThr_binning_, folderName_, jetAsyEtaME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetphiPrbME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, TriggerDQMBase::MEbinning::nbins, num_genTriggerEventFlag_, DQMStore::IBooker::setCurrentFolder(), TriggerDQMBase::setMETitle(), AlCaHLTBitMon_QueryRunRegistry::string, TriggerDQMBase::MEbinning::xmax, and TriggerDQMBase::MEbinning::xmin.

Referenced by ~DiJetMonitor().

 {  
   
   std::string histname, histtitle;
   std::string currentFolder = folderName_ ;
   ibooker.setCurrentFolder(currentFolder);
 
   histname = "jetpt1"; histtitle = "leading Jet Pt";
   bookME(ibooker,jetpt1ME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetpt1ME_,"Pt_1 [GeV]","events");
 
   histname = "jetpt2"; histtitle = "second leading Jet Pt";
   bookME(ibooker,jetpt2ME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetpt2ME_,"Pt_2 [GeV]","events");
 
   histname = "jetptAvgB"; histtitle = "Pt average before offline selection";
   bookME(ibooker,jetptAvgbME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptAvgbME_,"(pt_1 + pt_2)*0.5 [GeV]","events");
 
   histname = "jetptAvgA"; histtitle = "Pt average after offline selection";
   bookME(ibooker,jetptAvgaME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptAvgaME_,"(pt_1 + pt_2)*0.5 [GeV]","events");
 
   histname = "jetptAvgAThr"; histtitle = "Pt average after offline selection";
   bookME(ibooker,jetptAvgaThrME_,histname,histtitle,dijetptThr_binning_.nbins,dijetptThr_binning_.xmin, dijetptThr_binning_.xmax);
   setMETitle(jetptAvgaThrME_,"(pt_1 + pt_2)*0.5 [GeV]","events");
 
   histname = "jetptTag"; histtitle = "Tag Jet Pt";
   bookME(ibooker,jetptTagME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptTagME_,"Pt_tag [GeV]","events ");
 
   histname = "jetptPrb"; histtitle = "Probe Jet Pt";
   bookME(ibooker,jetptPrbME_,histname,histtitle,dijetpt_binning_.nbins,dijetpt_binning_.xmin,dijetpt_binning_.xmax);
   setMETitle(jetptPrbME_,"Pt_prb [GeV]","events");
 
   histname = "jetptAsym"; histtitle = "Jet Pt Asymetry";
   bookME(ibooker,jetptAsyME_,histname,histtitle,asy_binning.nbins,asy_binning.xmin, asy_binning.xmax);
   setMETitle(jetptAsyME_,"(pt_prb - pt_tag)/(pt_prb + pt_tag)","events");
 
   histname = "jetetaPrb"; histtitle = "Probe Jet eta";
   bookME(ibooker,jetetaPrbME_,histname,histtitle,dijet_eta_binning.nbins,dijet_eta_binning.xmin, dijet_eta_binning.xmax);
   setMETitle(jetetaPrbME_,"Eta_probe #eta","events");
 
   histname = "jetetaTag"; histtitle = "Tag Jet eta";
   bookME(ibooker,jetetaTagME_,histname,histtitle,dijet_eta_binning.nbins,dijet_eta_binning.xmin, dijet_eta_binning.xmax);
   setMETitle(jetetaTagME_,"Eta_tag #eta","events");
 
   histname = "ptAsymVSetaPrb"; histtitle = "Pt_Asym vs eta_prb";
   bookME(ibooker,jetAsyEtaME_,histname,histtitle,asy_binning.nbins, asy_binning.xmin, asy_binning.xmax, dijet_eta_binning.nbins, dijet_eta_binning.xmin,dijet_eta_binning.xmax);
   setMETitle(jetAsyEtaME_,"(pt_prb - pt_tag)/(pt_prb + pt_tag)","Eta_probe #eta");
 
   histname = "etaPrbVSphiPrb"; histtitle = "eta_prb vs phi_prb";
   bookME(ibooker,jetEtaPhiME_,histname,histtitle,dijet_eta_binning.nbins, dijet_eta_binning.xmin, dijet_eta_binning.xmax, dijet_phi_binning.nbins, dijet_phi_binning.xmin,dijet_phi_binning.xmax);
   setMETitle(jetEtaPhiME_,"Eta_probe #eta","Phi_probe #phi");
 
   histname = "jetphiPrb"; histtitle = "Probe Jet phi";
   bookME(ibooker,jetphiPrbME_,histname,histtitle,dijet_phi_binning.nbins,dijet_phi_binning.xmin, dijet_phi_binning.xmax);
   setMETitle(jetphiPrbME_,"Phi_probe #phi","events");
 
   // Initialize the GenericTriggerEventFlag
   // Initialize the GenericTriggerEventFlag
   if ( num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() ) num_genTriggerEventFlag_->initRun( iRun, iSetup );
   if ( den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on() ) den_genTriggerEventFlag_->initRun( iRun, iSetup );
 
 }

bool DiJetMonitor::dijet_selection	(	double	eta_1,
		double	phi_1,
		double	eta_2,
		double	phi_2,
		double	pt_1,
		double	pt_2,
		int &	tag_id,
		int &	probe_id,
		int	Event
	)

protected

Definition at line 257 of file DiJetMonitor.cc.

References funct::abs(), and DEFINE_FWK_MODULE.

Referenced by analyze(), and ~DiJetMonitor().

                                                                                                                                                          {
   double etacut = 1.7;
   double phicut = 2.7;
     
   bool passeta = (std::abs(eta_1)< etacut || std::abs(eta_2) < etacut );//check that one of the jets in the barrel
  
   float delta_phi_1_2= (phi_1 - phi_2);
   bool other_cuts = (std::abs(delta_phi_1_2) >= phicut); //check that jets are back to back
 
   if(std::abs(eta_1)<etacut && std::abs(eta_2)>etacut) {
       tag_id = 0; 
       probe_id = 1;
   }
   else 
   if(std::abs(eta_2)<etacut && std::abs(eta_1)>etacut) {
       tag_id = 1; 
       probe_id = 0;
   }
   else 
   if(std::abs(eta_2)<etacut && std::abs(eta_1)<etacut){
       int numb = Event % 2;
       if(numb==0){
         tag_id = 0; 
         probe_id = 1;
       }
       if(numb==1){
         tag_id = 1; 
         probe_id = 0;
       }
   }
   if (passeta && other_cuts)
     return true;
   else
     return false;
 }

void DiJetMonitor::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 209 of file DiJetMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), create_public_pileup_plots::bins, TriggerDQMBase::fillHistoLSPSetDescription(), TriggerDQMBase::fillHistoPSetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~DiJetMonitor().

 {
   edm::ParameterSetDescription desc;
   desc.add<std::string>  ( "FolderName", "HLT/JME/Jets/AK4/PF" );
 
   desc.add<edm::InputTag>( "met",      edm::InputTag("pfMet") );
   desc.add<edm::InputTag>( "dijetSrc",     edm::InputTag("ak4PFJets") );
   desc.add<edm::InputTag>( "electrons",edm::InputTag("gedGsfElectrons") );
   desc.add<edm::InputTag>( "muons",    edm::InputTag("muons") );
   desc.add<int>("njets",      0);
   desc.add<int>("nelectrons", 0);
   desc.add<double>("ptcut",   20);
 
   edm::ParameterSetDescription genericTriggerEventPSet;
   genericTriggerEventPSet.add<bool>("andOr");
   genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi") );
   genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions",{});
   genericTriggerEventPSet.add<bool>("andOrDcs", false);
   genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
   genericTriggerEventPSet.add<std::string>("dbLabel","");
   genericTriggerEventPSet.add<bool>("andOrHlt", true);
   genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT") );
   genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths",{});
   genericTriggerEventPSet.add<bool>("errorReplyHlt",false);
   genericTriggerEventPSet.add<unsigned int>("verbosityLevel",1);
 
   desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
   desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
 
   edm::ParameterSetDescription histoPSet;
   edm::ParameterSetDescription dijetPSet;
   edm::ParameterSetDescription dijetPtThrPSet;
   fillHistoPSetDescription(dijetPSet);
   fillHistoPSetDescription(dijetPtThrPSet);
   histoPSet.add<edm::ParameterSetDescription>("dijetPSet", dijetPSet);  
   histoPSet.add<edm::ParameterSetDescription>("dijetPtThrPSet", dijetPtThrPSet);  
   std::vector<double> bins = {0.,20.,40.,60.,80.,90.,100.,110.,120.,130.,140.,150.,160.,170.,180.,190.,200.,220.,240.,260.,280.,300.,350.,400.,450.,1000.}; // DiJet pT Binning
   histoPSet.add<std::vector<double> >("jetptBinning", bins);
 
   edm::ParameterSetDescription lsPSet;
   fillHistoLSPSetDescription(lsPSet);
   histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
   desc.add<edm::ParameterSetDescription>("histoPSet",histoPSet);
 
   descriptions.add("dijetMonitoring", desc);
 }

Member Data Documentation

MEbinning DiJetMonitor::asy_binning

private

Initial value:

{
    N_asy, MIN_asy, MAX_asy
  }

Definition at line 116 of file DiJetMonitor.h.

Referenced by bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::den_genTriggerEventFlag_

private

Definition at line 92 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

MEbinning DiJetMonitor::dijet_eta_binning

private

Initial value:

{
    DiJet_N_ETA, -DiJet_MAX_ETA, DiJet_MAX_ETA
  }

Definition at line 108 of file DiJetMonitor.h.

Referenced by bookHistograms().

const double DiJetMonitor::DiJet_MAX_ETA = 5

private

Definition at line 105 of file DiJetMonitor.h.

const double DiJetMonitor::DiJet_MAX_PHI = 3.2

private

Definition at line 98 of file DiJetMonitor.h.

unsigned int DiJetMonitor::DiJet_N_ETA = 20

private

Definition at line 107 of file DiJetMonitor.h.

unsigned int DiJetMonitor::DiJet_N_PHI = 32

private

Definition at line 100 of file DiJetMonitor.h.

MEbinning DiJetMonitor::dijet_phi_binning

private

Initial value:

{
    DiJet_N_PHI, -DiJet_MAX_PHI, DiJet_MAX_PHI
  }

Definition at line 101 of file DiJetMonitor.h.

Referenced by bookHistograms().

MEbinning DiJetMonitor::dijetpt_binning_

private

Definition at line 73 of file DiJetMonitor.h.

Referenced by bookHistograms(), and DiJetMonitor().

MEbinning DiJetMonitor::dijetptThr_binning_

private

Definition at line 74 of file DiJetMonitor.h.

Referenced by bookHistograms(), and DiJetMonitor().

edm::EDGetTokenT<reco::PFJetCollection> DiJetMonitor::dijetSrc_

private

Definition at line 70 of file DiJetMonitor.h.

Referenced by analyze(), and DiJetMonitor().

edm::EDGetTokenT<reco::GsfElectronCollection> DiJetMonitor::eleToken_

private

Definition at line 68 of file DiJetMonitor.h.

std::string DiJetMonitor::folderName_

private

Definition at line 64 of file DiJetMonitor.h.

Referenced by bookHistograms(), and DiJetMonitor().

std::string DiJetMonitor::histoSuffix_

private

Definition at line 65 of file DiJetMonitor.h.

ObjME DiJetMonitor::jetAsyEtaME_

private

Definition at line 88 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetEtaPhiME_

private

Definition at line 89 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetetaPrbME_

private

Definition at line 85 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetetaTagME_

private

Definition at line 86 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetphiPrbME_

private

Definition at line 87 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetpt1ME_

private

Definition at line 77 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetpt2ME_

private

Definition at line 78 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptAsyME_

private

Definition at line 84 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptAvgaME_

private

Definition at line 79 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptAvgaThrME_

private

Definition at line 80 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptAvgbME_

private

Definition at line 81 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptPrbME_

private

Definition at line 83 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

ObjME DiJetMonitor::jetptTagME_

private

Definition at line 82 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

const double DiJetMonitor::MAX_asy = 1

private

Definition at line 112 of file DiJetMonitor.h.

edm::EDGetTokenT<reco::PFMETCollection> DiJetMonitor::metToken_

private

Definition at line 67 of file DiJetMonitor.h.

const double DiJetMonitor::MIN_asy = -1

private

Definition at line 113 of file DiJetMonitor.h.

edm::EDGetTokenT<reco::MuonCollection> DiJetMonitor::muoToken_

private

Definition at line 69 of file DiJetMonitor.h.

unsigned int DiJetMonitor::N_asy = 50

private

Definition at line 115 of file DiJetMonitor.h.

int DiJetMonitor::nmuons_

private

Definition at line 94 of file DiJetMonitor.h.

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::num_genTriggerEventFlag_

private

Definition at line 91 of file DiJetMonitor.h.

Referenced by analyze(), and bookHistograms().

double DiJetMonitor::ptcut_

private

Definition at line 95 of file DiJetMonitor.h.

Referenced by analyze(), and DiJetMonitor().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation