Inheritance diagram for RazorMonitor:

Public Types
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Types inherited from TriggerDQMBase
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

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

	~RazorMonitor () override throw ()

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

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

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

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

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const double ymin, const double ymax, const bool bookDen=true)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const uint nbinsY, const double ymin, const double ymax, const bool bookDen=true)

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

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

	TriggerDQMBase ()=default

virtual	~TriggerDQMBase ()=default

Static Public Member Functions
static double	CalcMR (const math::XYZTLorentzVector &ja, const math::XYZTLorentzVector &jb)

static double	CalcR (double MR, const math::XYZTLorentzVector &ja, const math::XYZTLorentzVector &jb, const edm::Handle< std::vector< reco::PFMET > > &met)

static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
std::unique_ptr < GenericTriggerEventFlag >	den_genTriggerEventFlag_

std::vector< double >	dphiR_binning_

ObjME	dPhiR_ME_

const std::string	folderName_

bool	hltPathsAreValid_

StringCutObjectSelector < reco::PFJet, true >	jetSelection_

edm::EDGetTokenT < reco::PFJetCollection >	jetToken_

StringCutObjectSelector < reco::MET, true >	metSelection_

edm::EDGetTokenT < reco::PFMETCollection >	metToken_

std::vector< double >	mr_binning_

ObjME	MR_ME_

float	mrCut_

ObjME	MRVsRsq_ME_

unsigned int	njets_

std::unique_ptr < GenericTriggerEventFlag >	num_genTriggerEventFlag_

const bool	requireValidHLTPaths_

std::vector< double >	rsq_binning_

ObjME	Rsq_ME_

float	rsqCut_

edm::EDGetTokenT< std::vector < math::XYZTLorentzVector > >	theHemispheres_

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 31 of file RazorMonitor.cc.

Member Typedef Documentation

typedef dqm::reco::DQMStore RazorMonitor::DQMStore

Definition at line 34 of file RazorMonitor.cc.

typedef dqm::reco::MonitorElement RazorMonitor::MonitorElement

Definition at line 33 of file RazorMonitor.cc.

Constructor & Destructor Documentation

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

Definition at line 89 of file RazorMonitor.cc.

     : folderName_(iConfig.getParameter<std::string>("FolderName")),
       requireValidHLTPaths_(iConfig.getParameter<bool>("requireValidHLTPaths")),
       hltPathsAreValid_(false),
       metToken_(consumes<reco::PFMETCollection>(iConfig.getParameter<edm::InputTag>("met"))),
       jetToken_(mayConsume<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("jets"))),
       theHemispheres_(
           consumes<std::vector<math::XYZTLorentzVector> >(iConfig.getParameter<edm::InputTag>("hemispheres"))),
       rsq_binning_(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("rsqBins")),
       mr_binning_(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("mrBins")),
       dphiR_binning_(
           iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("dphiRBins")),
       num_genTriggerEventFlag_(new GenericTriggerEventFlag(
           iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
       den_genTriggerEventFlag_(new GenericTriggerEventFlag(
           iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)),
       metSelection_(iConfig.getParameter<std::string>("metSelection")),
       jetSelection_(iConfig.getParameter<std::string>("jetSelection")),
       njets_(iConfig.getParameter<unsigned int>("njets")),
       rsqCut_(iConfig.getParameter<double>("rsqCut")),
       mrCut_(iConfig.getParameter<double>("mrCut")) {}

RazorMonitor::~RazorMonitor	(	)
throw	(
	)

override

Definition at line 111 of file RazorMonitor.cc.

References den_genTriggerEventFlag_, and num_genTriggerEventFlag_.

                                     {
   if (num_genTriggerEventFlag_) {
     num_genTriggerEventFlag_.reset();
   }
   if (den_genTriggerEventFlag_) {
     den_genTriggerEventFlag_.reset();
   }
 }

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 168 of file RazorMonitor.cc.

References funct::abs(), CalcMR(), CalcR(), srCondWrite_cfg::deltaPhi, den_genTriggerEventFlag_, TriggerDQMBase::ObjME::denominator, dPhiR_ME_, relativeConstraints::empty, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), razorHemispheres_cff::hemispheres, hltPathsAreValid_, edm::HandleBase::isValid(), dqmiolumiharvest::j, fwrapper::jets, jetSelection_, jetToken_, metSelection_, metToken_, MR_ME_, mrCut_, MRVsRsq_ME_, njets_, num_genTriggerEventFlag_, TriggerDQMBase::ObjME::numerator, dttmaxenums::R, requireValidHLTPaths_, Rsq_ME_, rsqCut_, and theHemispheres_.

                                                                               {
   // if valid HLT paths are required,
   // analyze event only if all paths are valid
   if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
     return;
   }
 
   // Filter out events if Trigger Filtering is requested
   if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
     return;
 
   //met collection
   edm::Handle<reco::PFMETCollection> metHandle;
   iEvent.getByToken(metToken_, metHandle);
   reco::PFMET pfmet = metHandle->front();
   if (!metSelection_(pfmet))
     return;
 
   //jet collection, track # of jets for two working points
   edm::Handle<reco::PFJetCollection> jetHandle;
   iEvent.getByToken(jetToken_, jetHandle);
   std::vector<reco::PFJet> jets;
   if (jetHandle->size() < njets_)
     return;
   for (auto const& j : *jetHandle) {
     if (jetSelection_(j))
       jets.push_back(j);
   }
   if (jets.size() < njets_)
     return;
 
   //razor hemisphere clustering from previous step
   edm::Handle<vector<math::XYZTLorentzVector> > hemispheres;
   iEvent.getByToken(theHemispheres_, hemispheres);
 
   if (not hemispheres.isValid()) {
     return;
   }
 
   if (hemispheres
           ->empty()) {  // the Hemisphere Maker will produce an empty collection of hemispheres if # of jets is too high
     edm::LogError("DQM_HLT_Razor") << "Cannot calculate M_R and R^2 because there are too many jets! (trigger passed "
                                       "automatically without forming the hemispheres)"
                                    << endl;
     return;
   }
 
   // should always have 2 hemispheres -- no muons included (c. 2017), if not return invalid hemisphere collection
   // retaining check for hemisphere size 5 or 10 which correspond to the one or two muon case for possible future use
   if (!hemispheres->empty() && hemispheres->size() != 2 && hemispheres->size() != 5 && hemispheres->size() != 10) {
     edm::LogError("DQM_HLT_Razor") << "Invalid hemisphere collection!  hemispheres->size() = " << hemispheres->size()
                                    << endl;
     return;
   }
 
   //calculate razor variables, with hemispheres pT-ordered
   double MR = 0, R = 0;
   if (hemispheres->at(1).Pt() > hemispheres->at(0).Pt()) {
     MR = CalcMR(hemispheres->at(1), hemispheres->at(0));
     R = CalcR(MR, hemispheres->at(1), hemispheres->at(0), metHandle);
   } else {
     MR = CalcMR(hemispheres->at(0), hemispheres->at(1));
     R = CalcR(MR, hemispheres->at(0), hemispheres->at(1), metHandle);
   }
 
   double Rsq = R * R;
   double dPhiR = abs(deltaPhi(hemispheres->at(0).Phi(), hemispheres->at(1).Phi()));
 
   //apply offline selection cuts
   if (Rsq < rsqCut_ && MR < mrCut_)
     return;
 
   // applying selection for denominator
   if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
     return;
 
   // filling histograms (denominator)
   if (Rsq >= rsqCut_) {
     MR_ME_.denominator->Fill(MR);
   }
 
   if (MR >= mrCut_) {
     Rsq_ME_.denominator->Fill(Rsq);
   }
 
   dPhiR_ME_.denominator->Fill(dPhiR);
 
   MRVsRsq_ME_.denominator->Fill(MR, Rsq);
 
   // applying selection for numerator
   if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup))
     return;
 
   // filling histograms (numerator)
   if (Rsq >= rsqCut_) {
     MR_ME_.numerator->Fill(MR);
   }
 
   if (MR >= mrCut_) {
     Rsq_ME_.numerator->Fill(Rsq);
   }
 
   dPhiR_ME_.numerator->Fill(dPhiR);
 
   MRVsRsq_ME_.numerator->Fill(MR, Rsq);
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 120 of file RazorMonitor.cc.

References TriggerDQMBase::bookME(), den_genTriggerEventFlag_, dphiR_binning_, dPhiR_ME_, folderName_, hltPathsAreValid_, mr_binning_, MR_ME_, MRVsRsq_ME_, num_genTriggerEventFlag_, requireValidHLTPaths_, rsq_binning_, Rsq_ME_, dqm::implementation::NavigatorBase::setCurrentFolder(), TriggerDQMBase::setMETitle(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                            {
   // 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);
 
   // check if every HLT path specified in numerator and denominator has a valid match in the HLT Menu
   hltPathsAreValid_ = (num_genTriggerEventFlag_ && den_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() &&
                        den_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->allHLTPathsAreValid() &&
                        den_genTriggerEventFlag_->allHLTPathsAreValid());
 
   // if valid HLT paths are required,
   // create DQM outputs only if all paths are valid
   if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
     return;
   }
 
   std::string histname, histtitle;
 
   std::string currentFolder = folderName_;
   ibooker.setCurrentFolder(currentFolder);
 
   // 1D hist, MR
   histname = "MR";
   histtitle = "PF MR";
   bookME(ibooker, MR_ME_, histname, histtitle, mr_binning_);
   setMETitle(MR_ME_, "PF M_{R} [GeV]", "events / [GeV]");
 
   // 1D hist, Rsq
   histname = "Rsq";
   histtitle = "PF Rsq";
   bookME(ibooker, Rsq_ME_, histname, histtitle, rsq_binning_);
   setMETitle(Rsq_ME_, "PF R^{2}", "events");
 
   // 1D hist, dPhiR
   histname = "dPhiR";
   histtitle = "dPhiR";
   bookME(ibooker, dPhiR_ME_, histname, histtitle, dphiR_binning_);
   setMETitle(dPhiR_ME_, "dPhi_{R}", "events");
 
   // 2D hist, MR & Rsq
   histname = "MRVsRsq";
   histtitle = "PF MR vs PF Rsq";
   bookME(ibooker, MRVsRsq_ME_, histname, histtitle, mr_binning_, rsq_binning_);
   setMETitle(MRVsRsq_ME_, "M_{R} [GeV]", "R^{2}");
 }

double RazorMonitor::CalcMR	(	const math::XYZTLorentzVector &	ja,
		const math::XYZTLorentzVector &	jb
	)

static

Definition at line 314 of file RazorMonitor.cc.

References A, B, and mathSSE::sqrt().

Referenced by analyze().

                                                                                             {
   if (ja.Pt() <= 0.1)
     return -1;
 
   double A = ja.P();
   double B = jb.P();
   double az = ja.Pz();
   double bz = jb.Pz();
   TVector3 jaT, jbT;
   jaT.SetXYZ(ja.Px(), ja.Py(), 0.0);
   jbT.SetXYZ(jb.Px(), jb.Py(), 0.0);
   double ATBT = (jaT + jbT).Mag2();
 
   double MR = sqrt((A + B) * (A + B) - (az + bz) * (az + bz) -
                    (jbT.Dot(jbT) - jaT.Dot(jaT)) * (jbT.Dot(jbT) - jaT.Dot(jaT)) / (jaT + jbT).Mag2());
   double mybeta = (jbT.Dot(jbT) - jaT.Dot(jaT)) / sqrt(ATBT * ((A + B) * (A + B) - (az + bz) * (az + bz)));
 
   double mygamma = 1. / sqrt(1. - mybeta * mybeta);
 
   //use gamma times MRstar
   return MR * mygamma;
 }

double RazorMonitor::CalcR	(	double	MR,
		const math::XYZTLorentzVector &	ja,
		const math::XYZTLorentzVector &	jb,
		const edm::Handle< std::vector< reco::PFMET > > &	met
	)

static

Definition at line 337 of file RazorMonitor.cc.

References objects.METAnalyzer::met, and mathSSE::sqrt().

Referenced by analyze().

                                                                              {
   //now we can calculate MTR
   const math::XYZVector met = (inputMet->front()).momentum();
 
   double MTR = sqrt(0.5 * (met.R() * (ja.Pt() + jb.Pt()) - met.Dot(ja.Vect() + jb.Vect())));
 
   //filter events
   return float(MTR) / float(MR);  //R
 }

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

static

Definition at line 275 of file RazorMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, GenericTriggerEventFlag::fillPSetDescription(), HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                               {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("FolderName", "HLT/SUSY/Razor");
   desc.add<bool>("requireValidHLTPaths", true);
 
   desc.add<edm::InputTag>("met", edm::InputTag("pfMet"));
   desc.add<edm::InputTag>("jets", edm::InputTag("ak4PFJetsCHS"));
   desc.add<edm::InputTag>("hemispheres", edm::InputTag("hemispheresDQM"))
       ->setComment("hemisphere jets used to compute razor variables");
   desc.add<std::string>("metSelection", "pt > 0");
 
   // from 2016 offline selection
   desc.add<std::string>("jetSelection", "pt > 80");
   desc.add<unsigned int>("njets", 2);
   desc.add<double>("mrCut", 300);
   desc.add<double>("rsqCut", 0.15);
 
   edm::ParameterSetDescription genericTriggerEventPSet;
   GenericTriggerEventFlag::fillPSetDescription(genericTriggerEventPSet);
   desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
   desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
 
   //binning from 2016 offline selection
   edm::ParameterSetDescription histoPSet;
   std::vector<double> mrbins = {0., 100., 200., 300., 400., 500., 575., 650., 750., 900., 1200., 1600., 2500., 4000.};
   histoPSet.add<std::vector<double> >("mrBins", mrbins);
 
   std::vector<double> rsqbins = {0., 0.05, 0.1, 0.15, 0.2, 0.25, 0.30, 0.41, 0.52, 0.64, 0.8, 1.5};
   histoPSet.add<std::vector<double> >("rsqBins", rsqbins);
 
   std::vector<double> dphirbins = {0., 0.5, 1.0, 1.5, 2.0, 2.5, 2.8, 3.0, 3.2};
   histoPSet.add<std::vector<double> >("dphiRBins", dphirbins);
 
   desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
 
   descriptions.add("razorMonitoring", desc);
 }