BPHMonitor Class Reference

#include <BPHMonitor.h>

Inheritance diagram for BPHMonitor:

Public Member Functions
	BPHMonitor (const edm::ParameterSet &)
void	case11_selection (const float &dimuonCL, const float &jpsi_cos, const GlobalPoint &displacementFromBeamspotJpsi, const GlobalError &jerr, const edm::Handle< reco::TrackCollection > &trHandle, const std::string &hltpath, const edm::Handle< trigger::TriggerEvent > &handleTriggerEvent, const reco::Muon &m, const reco::Muon &m1, const edm::ESHandle< MagneticField > &bFieldHandle, const reco::BeamSpot &vertexBeamSpot, MonitorElement *phi1, MonitorElement *eta1, MonitorElement *pT1, MonitorElement *phi2, MonitorElement *eta2, MonitorElement *pT2)
	~BPHMonitor () override
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
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () 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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)
static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)
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
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, int &nbins, double &xmin, double &xmax)
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, std::vector< double > binningX)
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, int &nbinsX, double &xmin, double &xmax, double &ymin, double &ymax)
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, int &nbinsX, double &xmin, double &xmax, int &nbinsY, double &ymin, double &ymax)
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, std::vector< double > binningX, std::vector< double > binningY)
void	bookME (DQMStore::IBooker &, METME &me, std::string &histname, std::string &histtitle, MEbinning &binning)
template<typename T >
bool	matchToTrigger (const std::string &theTriggerName, T t)
double	Prescale (const std::string num, const std::string den, edm::Event const &iEvent, edm::EventSetup const &iSetup, HLTPrescaleProvider *hltPrescale_)
void	setMETitle (METME &me, std::string titleX, std::string titleY)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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
std::string	getTriggerName (std::string partialName)

Static Private Member Functions
static MEbinning	getHistoPSet (edm::ParameterSet pset)

Private Attributes
METME	BMass_
MEbinning	Bmass_binning_
edm::InputTag	bsInputTag_
edm::EDGetTokenT< reco::BeamSpot >	bsToken_
MEbinning	cos_binning_
MEbinning	d0_binning_
MEbinning	dca_binning_
GenericTriggerEventFlag *	den_genTriggerEventFlag_
METME	DiMuDCA_
METME	DiMudR_
METME	DiMuDS_
METME	DiMuEta_
METME	DiMuMass_
METME	DiMuPhi_
METME	DiMuProb_
METME	DiMuPt_
METME	DiMuPVcos_
int	displaced_
StringCutObjectSelector< reco::Candidate::LorentzVector, true >	DMSelection_ref
MEbinning	dMu_pt_binning_
MEbinning	dR_binning_
MEbinning	ds_binning_
int	enum_
MEbinning	eta_binning_
std::string	folderName_
edm::Handle< trigger::TriggerEvent >	handleTriggerEvent
std::string	histoSuffix_
HLTConfigProvider	hltConfig_
edm::EDGetTokenT< trigger::TriggerEvent >	hltInputTag_
edm::InputTag	hltInputTag_1
std::vector< std::string >	hltpaths_den
std::vector< std::string >	hltpaths_num
HLTPrescaleProvider *	hltPrescale_
edm::Handle< edm::TriggerResults >	HLTR
edm::EDGetTokenT< edm::TriggerResults >	hltTrigResTag_
int	Jpsi_
double	kaon_mass
int	L3_
MEbinning	mass_binning_
double	max_dR
double	maxmass_
double	maxmassJpsi
double	maxmassJpsiTk
double	maxmassTkTk
double	maxmassUpsilon
double	min_dR
double	mincos
double	minDS
double	minmass_
double	minmassJpsi
double	minmassJpsiTk
double	minmassTkTk
double	minmassUpsilon
double	minprob
METME	mu1d0_
METME	mu1Eta_
METME	mu1Phi_
METME	mu1Pt_
METME	mu1z0_
METME	mu2d0_
METME	mu2Eta_
METME	mu2Phi_
METME	mu2Pt_
METME	mu2z0_
METME	mu3d0_
METME	mu3Eta_
METME	mu3Phi_
METME	mu3Pt_
METME	mu3z0_
double	mu_mass
METME	mud0_
METME	muEta_
edm::InputTag	muoInputTag_
StringCutObjectSelector< reco::Muon, true >	muoSelection_
StringCutObjectSelector< reco::Muon, true >	muoSelection_probe
StringCutObjectSelector< reco::Muon, true >	muoSelection_ref
StringCutObjectSelector< reco::Muon, true >	muoSelection_tag
edm::EDGetTokenT< reco::MuonCollection >	muoToken_
METME	muPhi_
METME	muPt_
METME	muz0_
int	nmuons_
GenericTriggerEventFlag *	num_genTriggerEventFlag_
METME	phEta_
MEbinning	phi_binning_
edm::InputTag	phInputTag_
METME	phPhi_
METME	phPt_
edm::EDGetTokenT< reco::PhotonCollection >	phToken_
MEbinning	prob_binning_
MEbinning	pt_binning_
int	ptCut_
int	seagull_
bool	tnp_
edm::InputTag	trInputTag_
int	trOrMu_
StringCutObjectSelector< reco::Track, true >	trSelection_
StringCutObjectSelector< reco::Track, true >	trSelection_ref
edm::EDGetTokenT< reco::TrackCollection >	trToken_
int	Upsilon_
edm::InputTag	vtxInputTag_
edm::EDGetTokenT< reco::VertexCollection >	vtxToken_
std::vector< bool >	warningPrinted4token_
MEbinning	z0_binning_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
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 76 of file BPHMonitor.h.

Constructor & Destructor Documentation

BPHMonitor::BPHMonitor

(

const edm::ParameterSet &

iConfig

)

Definition at line 11 of file BPHMonitor.cc.

References BMass_, METME::denominator, DiMuDCA_, DiMudR_, DiMuDS_, DiMuEta_, DiMuMass_, DiMuPhi_, DiMuProb_, DiMuPt_, DiMuPVcos_, mu1Eta_, mu1Phi_, mu1Pt_, mu2Eta_, mu2Phi_, mu2Pt_, mu3Eta_, mu3Phi_, mu3Pt_, mud0_, muEta_, muPhi_, muPt_, muz0_, METME::numerator, phEta_, phPhi_, phPt_, and warningPrinted4token_.

    : folderName_(iConfig.getParameter<std::string>("FolderName")),
      muoInputTag_(iConfig.getParameter<edm::InputTag>("muons")),
      bsInputTag_(iConfig.getParameter<edm::InputTag>("beamSpot")),
      trInputTag_(iConfig.getParameter<edm::InputTag>("tracks")),
      phInputTag_(iConfig.getParameter<edm::InputTag>("photons")),
      vtxInputTag_(iConfig.getParameter<edm::InputTag>("offlinePVs")),
      muoToken_(mayConsume<reco::MuonCollection>(muoInputTag_)),
      bsToken_(mayConsume<reco::BeamSpot>(bsInputTag_)),
      trToken_(mayConsume<reco::TrackCollection>(trInputTag_)),
      phToken_(mayConsume<reco::PhotonCollection>(phInputTag_)),
      vtxToken_(mayConsume<reco::VertexCollection>(vtxInputTag_)),
      phi_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("phiPSet"))),
      pt_binning_(
          getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("ptPSet"))),
      dMu_pt_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dMu_ptPSet"))),
      eta_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("etaPSet"))),
      d0_binning_(
          getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("d0PSet"))),
      z0_binning_(
          getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("z0PSet"))),
      dR_binning_(
          getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dRPSet"))),
      mass_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("massPSet"))),
      Bmass_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("BmassPSet"))),
      dca_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dcaPSet"))),
      ds_binning_(
          getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dsPSet"))),
      cos_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("cosPSet"))),
      prob_binning_(getHistoPSet(
          iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("probPSet"))),
      num_genTriggerEventFlag_(new GenericTriggerEventFlag(
          iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
      den_genTriggerEventFlag_(new GenericTriggerEventFlag(
          iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)),
      hltPrescale_(new HLTPrescaleProvider(iConfig, consumesCollector(), *this)),
      muoSelection_(iConfig.getParameter<std::string>("muoSelection")),
      muoSelection_ref(iConfig.getParameter<std::string>("muoSelection_ref")),
      muoSelection_tag(iConfig.getParameter<std::string>("muoSelection_tag")),
      muoSelection_probe(iConfig.getParameter<std::string>("muoSelection_probe")),
      nmuons_(iConfig.getParameter<int>("nmuons")),
      tnp_(iConfig.getParameter<bool>("tnp")),
      L3_(iConfig.getParameter<int>("L3")),
      ptCut_(iConfig.getParameter<int>("ptCut")),
      displaced_(iConfig.getParameter<int>("displaced")),
      trOrMu_(iConfig.getParameter<int>("trOrMu")),
      Jpsi_(iConfig.getParameter<int>("Jpsi")),
      Upsilon_(iConfig.getParameter<int>("Upsilon"))  // if ==1 path with Upsilon constraint
      ,
      enum_(iConfig.getParameter<int>("enum")),
      seagull_(iConfig.getParameter<int>("seagull")),
      maxmass_(iConfig.getParameter<double>("maxmass")),
      minmass_(iConfig.getParameter<double>("minmass")),
      maxmassJpsi(iConfig.getParameter<double>("maxmassJpsi")),
      minmassJpsi(iConfig.getParameter<double>("minmassJpsi")),
      maxmassUpsilon(iConfig.getParameter<double>("maxmassUpsilon")),
      minmassUpsilon(iConfig.getParameter<double>("minmassUpsilon")),
      maxmassTkTk(iConfig.getParameter<double>("maxmassTkTk")),
      minmassTkTk(iConfig.getParameter<double>("minmassTkTk")),
      maxmassJpsiTk(iConfig.getParameter<double>("maxmassJpsiTk")),
      minmassJpsiTk(iConfig.getParameter<double>("minmassJpsiTk")),
      kaon_mass(iConfig.getParameter<double>("kaon_mass")),
      mu_mass(iConfig.getParameter<double>("mu_mass")),
      min_dR(iConfig.getParameter<double>("min_dR")),
      max_dR(iConfig.getParameter<double>("max_dR")),
      minprob(iConfig.getParameter<double>("minprob")),
      mincos(iConfig.getParameter<double>("mincos")),
      minDS(iConfig.getParameter<double>("minDS")),
      hltInputTag_1(iConfig.getParameter<edm::InputTag>("hltTriggerSummaryAOD")),
      hltInputTag_(mayConsume<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("hltTriggerSummaryAOD"))),
      hltpaths_num(iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet")
                       .getParameter<std::vector<std::string>>("hltPaths")),
      hltpaths_den(iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet")
                       .getParameter<std::vector<std::string>>("hltPaths")),
      trSelection_(iConfig.getParameter<std::string>("muoSelection")),
      trSelection_ref(iConfig.getParameter<std::string>("trSelection_ref")),
      DMSelection_ref(iConfig.getParameter<std::string>("DMSelection_ref")) {
  muPhi_.numerator = nullptr;
  muPhi_.denominator = nullptr;
  muEta_.numerator = nullptr;
  muEta_.denominator = nullptr;
  muPt_.numerator = nullptr;
  muPt_.denominator = nullptr;
  mud0_.numerator = nullptr;
  mud0_.denominator = nullptr;
  muz0_.numerator = nullptr;
  muz0_.denominator = nullptr;

  mu1Phi_.numerator = nullptr;
  mu1Phi_.denominator = nullptr;
  mu1Eta_.numerator = nullptr;
  mu1Eta_.denominator = nullptr;
  mu1Pt_.numerator = nullptr;
  mu1Pt_.denominator = nullptr;

  mu2Phi_.numerator = nullptr;
  mu2Phi_.denominator = nullptr;
  mu2Eta_.numerator = nullptr;
  mu2Eta_.denominator = nullptr;
  mu2Pt_.numerator = nullptr;
  mu2Pt_.denominator = nullptr;

  mu3Phi_.numerator = nullptr;
  mu3Phi_.denominator = nullptr;
  mu3Eta_.numerator = nullptr;
  mu3Eta_.denominator = nullptr;
  mu3Pt_.numerator = nullptr;
  mu3Pt_.denominator = nullptr;

  phPhi_.numerator = nullptr;
  phPhi_.denominator = nullptr;
  phEta_.numerator = nullptr;
  phEta_.denominator = nullptr;
  phPt_.numerator = nullptr;
  phPt_.denominator = nullptr;

  DiMuPhi_.numerator = nullptr;
  DiMuPhi_.denominator = nullptr;
  DiMuEta_.numerator = nullptr;
  DiMuEta_.denominator = nullptr;
  DiMuPt_.numerator = nullptr;
  DiMuPt_.denominator = nullptr;
  DiMuPVcos_.numerator = nullptr;
  DiMuPVcos_.denominator = nullptr;
  DiMuProb_.numerator = nullptr;
  DiMuProb_.denominator = nullptr;
  DiMuDS_.numerator = nullptr;
  DiMuDS_.denominator = nullptr;
  DiMuDCA_.numerator = nullptr;
  DiMuDCA_.denominator = nullptr;
  DiMuMass_.numerator = nullptr;
  DiMuMass_.denominator = nullptr;
  BMass_.numerator = nullptr;
  BMass_.denominator = nullptr;
  DiMudR_.numerator = nullptr;
  DiMudR_.denominator = nullptr;

  // this vector has to be alligned to the the number of Tokens accessed by this module
  warningPrinted4token_.push_back(false);  // MuonCollection
  warningPrinted4token_.push_back(false);  // BeamSpot
  warningPrinted4token_.push_back(false);  // TrackCollection
  warningPrinted4token_.push_back(false);  // PhotonCollection
  warningPrinted4token_.push_back(false);  // VertexCollection
}

BPHMonitor::~BPHMonitor ( )

override

Definition at line 163 of file BPHMonitor.cc.

References den_genTriggerEventFlag_, hltPrescale_, and num_genTriggerEventFlag_.

                        {
  if (num_genTriggerEventFlag_)
    delete num_genTriggerEventFlag_;
  if (den_genTriggerEventFlag_)
    delete den_genTriggerEventFlag_;
  delete hltPrescale_;
}

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 436 of file BPHMonitor.cc.

References GenericTriggerEventFlag::accept(), pwdgSkimBPark_cfi::beamSpot, BMass_, bsInputTag_, bsToken_, ClosestApproachInRPhi::calculate(), reco::Vertex::chi2(), SiPixelRawToDigiRegional_cfi::deltaPhi, reco::deltaR(), den_genTriggerEventFlag_, METME::denominator, DiMuDCA_, DiMudR_, DiMuDS_, DiMuEta_, DiMuMass_, DiMuPhi_, DiMuProb_, DiMuPt_, DiMuPVcos_, displaced_, ClosestApproachInRPhi::distance(), DMSelection_ref, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), StorageManager_cfg::e1, SiPixelPhase1Clusters_cfi::e3, vertexPlots::e4, enum_, runTheMatrix::err, dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), getTriggerName(), handleTriggerEvent, reco::TrackBase::highPurity, hltInputTag_, hltpaths_den, hltpaths_num, hltPrescale_, mps_fire::i, createfilelist::int, edm::HandleBase::isValid(), TrajectoryStateClosestToPoint::isValid(), TransientVertex::isValid(), Jpsi_, kaon_mass, edm::InputTag::label(), visualization-live-secondInstance_cfg::m, matchToTrigger(), maxmass_, maxmassJpsi, maxmassJpsiTk, maxmassTkTk, maxmassUpsilon, min_dR, mincos, minDS, minmass_, minmassJpsi, minmassJpsiTk, minmassTkTk, minmassUpsilon, minprob, reco::TrackBase::momentum(), mu1Eta_, mu1Phi_, mu1Pt_, mu2Eta_, mu2Phi_, mu2Pt_, mu3Eta_, mu3Phi_, mu3Pt_, mu_mass, muEta_, muoInputTag_, muoSelection_, muoSelection_ref, muoToken_, muPhi_, muPt_, reco::Vertex::ndof(), num_genTriggerEventFlag_, METME::numerator, GenericTriggerEventFlag::on(), AlCaHLTBitMon_ParallelJobs::p, p1, p2, p3, p4, PV3DBase< T, PVType, FrameType >::perp(), phEta_, VtxSmearedParameters_cfi::Phi, phInputTag_, phPhi_, phPt_, phToken_, TransientVertex::position(), TransientVertex::positionError(), Prescale(), DiDispStaMuonMonitor_cfi::pt, ptCut_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), reco::TrackBase::quality(), GlobalErrorBase< T, ErrorWeightType >::rerr(), seagull_, singleTopDQM_miniAOD_cfi::secondaryVertex, trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, OrderedSet::t, RandomServiceHelper::t1, TrajectoryStateClosestToPoint::theState(), tnp_, trInputTag_, trSelection_ref, trToken_, Upsilon_, bphysicsOniaDQM_cfi::vertex, KalmanVertexFitter::vertex(), warningPrinted4token_, PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x0(), PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y0(), PV3DBase< T, PVType, FrameType >::z(), and reco::BeamSpot::z0().

                                                                            {
  edm::Handle<reco::BeamSpot> beamSpot;
  iEvent.getByToken(bsToken_, beamSpot);
  if (!beamSpot.isValid()) {
    if (!warningPrinted4token_[0]) {
      warningPrinted4token_[0] = true;
      edm::LogWarning("BPHMonitor") << "skipping events because the collection " << bsInputTag_.label().c_str()
                                    << " is not available";
    }
    return;
  }

  edm::Handle<reco::MuonCollection> muoHandle;
  iEvent.getByToken(muoToken_, muoHandle);
  if (!muoHandle.isValid()) {
    if (!warningPrinted4token_[1]) {
      warningPrinted4token_[1] = true;
      edm::LogWarning("BPHMonitor") << "skipping events because the collection " << muoInputTag_.label().c_str()
                                    << " is not available";
    }
    return;
  }

  edm::Handle<reco::TrackCollection> trHandle;
  iEvent.getByToken(trToken_, trHandle);
  if (!trHandle.isValid()) {
    if (!warningPrinted4token_[2]) {
      warningPrinted4token_[2] = true;
      edm::LogWarning("BPHMonitor") << "skipping events because the collection " << trInputTag_.label().c_str()
                                    << " is not available";
    }
    return;
  }

  edm::Handle<reco::PhotonCollection> phHandle;
  iEvent.getByToken(phToken_, phHandle);

  edm::Handle<edm::TriggerResults> handleTriggerTrigRes;

  edm::ESHandle<MagneticField> bFieldHandle;

  const std::string& hltpath = getTriggerName(hltpaths_den[0]);
  const std::string& hltpath1 = getTriggerName(hltpaths_num[0]);

  double PrescaleWeight = 1.0;
  if (den_genTriggerEventFlag_->on() && den_genTriggerEventFlag_->accept(iEvent, iSetup) &&
      num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup))
    PrescaleWeight = Prescale(hltpath1, hltpath, iEvent, iSetup, hltPrescale_);

  if (tnp_ > 0) {  //TnP method

    if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
      return;
    iEvent.getByToken(hltInputTag_, handleTriggerEvent);
    if (handleTriggerEvent->sizeFilters() == 0)
      return;

    std::vector<reco::Muon> tagMuons;
    for (auto const& m : *muoHandle) {  // applying tag selection
      if (!matchToTrigger(hltpath, m))
        continue;
      if (muoSelection_ref(m))
        tagMuons.push_back(m);
    }

    for (int i = 0; i < int(tagMuons.size()); i++) {
      for (auto const& m : *muoHandle) {
        if (!matchToTrigger(hltpath, m))
          continue;
        if ((tagMuons[i].pt() == m.pt()))
          continue;  //not the same
        if ((tagMuons[i].p4() + m.p4()).M() > minmass_ &&
            (tagMuons[i].p4() + m.p4()).M() < maxmass_) {  //near to J/psi mass
          muPhi_.denominator->Fill(m.phi());
          muEta_.denominator->Fill(m.eta());
          muPt_.denominator->Fill(m.pt());
          if (muoSelection_(m) && num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
            muPhi_.numerator->Fill(m.phi(), PrescaleWeight);
            muEta_.numerator->Fill(m.eta(), PrescaleWeight);
            muPt_.numerator->Fill(m.pt(), PrescaleWeight);
          }
        }
      }
    }

  } else {  // reference method

    if (den_genTriggerEventFlag_->on() && (!den_genTriggerEventFlag_->accept(iEvent, iSetup)))
      return;

    iEvent.getByToken(hltInputTag_, handleTriggerEvent);
    if (handleTriggerEvent->sizeFilters() == 0)
      return;

    for (auto const& m : *muoHandle) {
      if (!muoSelection_ref(m))
        continue;
      if (!matchToTrigger(hltpath, m))
        continue;

      for (auto const& m1 : *muoHandle) {
        if (!(m1.pt() > m.pt()))
          continue;
        if (ptCut_) {
          if (!muoSelection_(m1))
            continue;
        } else if (!muoSelection_ref(m1))
          continue;
        if (!matchToTrigger(hltpath, m1))
          continue;

        if (enum_ != 10) {
          if (!DMSelection_ref(m1.p4() + m.p4()))
            continue;
          if (m.charge() * m1.charge() > 0)
            continue;
        }

        // dimuon vertex reconstruction
        iSetup.get<IdealMagneticFieldRecord>().get(bFieldHandle);
        const reco::BeamSpot& vertexBeamSpot = *beamSpot;
        std::vector<reco::TransientTrack> j_tks;
        reco::TransientTrack mu1TT(m.track(), &(*bFieldHandle));
        reco::TransientTrack mu2TT(m1.track(), &(*bFieldHandle));
        j_tks.push_back(mu1TT);
        j_tks.push_back(mu2TT);
        KalmanVertexFitter jkvf;
        TransientVertex jtv = jkvf.vertex(j_tks);
        if (!jtv.isValid())
          continue;
        reco::Vertex jpsivertex = jtv;
        float dimuonCL = 0;
        if ((jpsivertex.chi2() >= 0) && (jpsivertex.ndof() > 0))
          dimuonCL = TMath::Prob(jpsivertex.chi2(), jpsivertex.ndof());
        math::XYZVector jpperp(m.px() + m1.px(), m.py() + m1.py(), 0.);
        GlobalPoint jVertex = jtv.position();
        GlobalError jerr = jtv.positionError();
        GlobalPoint displacementFromBeamspotJpsi(
            -1 * ((vertexBeamSpot.x0() - jVertex.x()) + (jVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dxdz()),
            -1 * ((vertexBeamSpot.y0() - jVertex.y()) + (jVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dydz()),
            0);
        reco::Vertex::Point vperpj(displacementFromBeamspotJpsi.x(), displacementFromBeamspotJpsi.y(), 0.);
        float jpsi_cos = vperpj.Dot(jpperp) / (vperpj.R() * jpperp.R());
        TrajectoryStateClosestToPoint mu1TS = mu1TT.impactPointTSCP();
        TrajectoryStateClosestToPoint mu2TS = mu2TT.impactPointTSCP();
        ClosestApproachInRPhi cApp;
        if (mu1TS.isValid() && mu2TS.isValid()) {
          if (!cApp.calculate(mu1TS.theState(), mu2TS.theState()))
            continue;
        } else
          continue;
        double DiMuMass = (m1.p4() + m.p4()).M();

        switch (
            enum_) {  // enum_ = 1...9, represents different sets of variables for different paths, we want to have different hists for different paths

          case 1:

            tnp_ = true;  // already filled hists for tnp method
            [[fallthrough]];
          case 2:

            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;
            if (dimuonCL < minprob)
              continue;

            mu1Phi_.denominator->Fill(m.phi());
            mu1Eta_.denominator->Fill(m.eta());
            mu1Pt_.denominator->Fill(m.pt());
            mu2Phi_.denominator->Fill(m1.phi());
            mu2Eta_.denominator->Fill(m1.eta());
            mu2Pt_.denominator->Fill(m1.pt());
            DiMuPt_.denominator->Fill((m1.p4() + m.p4()).Pt());
            DiMuEta_.denominator->Fill((m1.p4() + m.p4()).Eta());
            DiMuPhi_.denominator->Fill((m1.p4() + m.p4()).Phi());

            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
              if (!matchToTrigger(hltpath1, m1))
                continue;
              if (!matchToTrigger(hltpath1, m))
                continue;
              mu1Phi_.numerator->Fill(m.phi(), PrescaleWeight);
              mu1Eta_.numerator->Fill(m.eta(), PrescaleWeight);
              mu1Pt_.numerator->Fill(m.pt(), PrescaleWeight);
              mu2Phi_.numerator->Fill(m1.phi(), PrescaleWeight);
              mu2Eta_.numerator->Fill(m1.eta(), PrescaleWeight);
              mu2Pt_.numerator->Fill(m1.pt(), PrescaleWeight);
              DiMuPt_.numerator->Fill((m1.p4() + m.p4()).Pt(), PrescaleWeight);
              DiMuEta_.numerator->Fill((m1.p4() + m.p4()).Eta(), PrescaleWeight);
              DiMuPhi_.numerator->Fill((m1.p4() + m.p4()).Phi(), PrescaleWeight);
            }

            break;

          case 3:

            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;
            if (dimuonCL < minprob)
              continue;

            mu1Phi_.denominator->Fill(m.phi());
            mu1Eta_.denominator->Fill(m.eta());
            mu1Pt_.denominator->Fill(m.pt());
            mu2Phi_.denominator->Fill(m1.phi());
            mu2Eta_.denominator->Fill(m1.eta());
            mu2Pt_.denominator->Fill(m1.pt());

            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
              if (!matchToTrigger(hltpath1, m1))
                continue;
              if (!matchToTrigger(hltpath1, m))
                continue;
              mu1Phi_.numerator->Fill(m.phi(), PrescaleWeight);
              mu1Eta_.numerator->Fill(m.eta(), PrescaleWeight);
              mu1Pt_.numerator->Fill(m.pt(), PrescaleWeight);
              mu2Phi_.numerator->Fill(m1.phi(), PrescaleWeight);
              mu2Eta_.numerator->Fill(m1.eta(), PrescaleWeight);
              mu2Pt_.numerator->Fill(m1.pt(), PrescaleWeight);
            }

            break;

          case 4:

            if (dimuonCL < minprob)
              continue;

            // fill mass plots without selecting mass region
            DiMuMass_.denominator->Fill(DiMuMass);
            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup) &&
                !(seagull_ && m.charge() * deltaPhi(m.phi(), m1.phi()) > 0) && matchToTrigger(hltpath1, m1) &&
                matchToTrigger(hltpath1, m))
              DiMuMass_.numerator->Fill(DiMuMass, PrescaleWeight);

            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;

            mu1Phi_.denominator->Fill(m.phi());
            mu1Eta_.denominator->Fill(m.eta());
            mu1Pt_.denominator->Fill(m.pt());
            mu2Phi_.denominator->Fill(m1.phi());
            mu2Eta_.denominator->Fill(m1.eta());
            mu2Pt_.denominator->Fill(m1.pt());
            DiMuPt_.denominator->Fill((m1.p4() + m.p4()).Pt());
            DiMuEta_.denominator->Fill((m1.p4() + m.p4()).Eta());
            DiMuPhi_.denominator->Fill((m1.p4() + m.p4()).Phi());
            DiMudR_.denominator->Fill(reco::deltaR(m, m1));

            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
              if (seagull_ && m.charge() * deltaPhi(m.phi(), m1.phi()) > 0)
                continue;
              if (!matchToTrigger(hltpath1, m1))
                continue;
              if (!matchToTrigger(hltpath1, m))
                continue;
              mu1Phi_.numerator->Fill(m.phi(), PrescaleWeight);
              mu1Eta_.numerator->Fill(m.eta(), PrescaleWeight);
              mu1Pt_.numerator->Fill(m.pt(), PrescaleWeight);
              mu2Phi_.numerator->Fill(m1.phi(), PrescaleWeight);
              mu2Eta_.numerator->Fill(m1.eta(), PrescaleWeight);
              mu2Pt_.numerator->Fill(m1.pt(), PrescaleWeight);
              DiMuPt_.numerator->Fill((m1.p4() + m.p4()).Pt(), PrescaleWeight);
              DiMuEta_.numerator->Fill((m1.p4() + m.p4()).Eta(), PrescaleWeight);
              DiMuPhi_.numerator->Fill((m1.p4() + m.p4()).Phi(), PrescaleWeight);
              DiMudR_.numerator->Fill(reco::deltaR(m, m1), PrescaleWeight);
            }

            break;

          case 5:

            if (dimuonCL < minprob)
              continue;
            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;

            mu1Phi_.denominator->Fill(m.phi());
            mu1Eta_.denominator->Fill(m.eta());
            mu1Pt_.denominator->Fill(m.pt());
            mu2Phi_.denominator->Fill(m1.phi());
            mu2Eta_.denominator->Fill(m1.eta());
            mu2Pt_.denominator->Fill(m1.pt());
            DiMuPt_.denominator->Fill((m1.p4() + m.p4()).Pt());
            DiMuEta_.denominator->Fill((m1.p4() + m.p4()).Eta());
            DiMuPhi_.denominator->Fill((m1.p4() + m.p4()).Phi());
            DiMudR_.denominator->Fill(reco::deltaR(m, m1));

            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
              if (seagull_ && m.charge() * deltaPhi(m.phi(), m1.phi()) > 0)
                continue;
              if (!matchToTrigger(hltpath1, m1))
                continue;
              if (!matchToTrigger(hltpath1, m))
                continue;
              mu1Phi_.numerator->Fill(m.phi(), PrescaleWeight);
              mu1Eta_.numerator->Fill(m.eta(), PrescaleWeight);
              mu1Pt_.numerator->Fill(m.pt(), PrescaleWeight);
              mu2Phi_.numerator->Fill(m1.phi(), PrescaleWeight);
              mu2Eta_.numerator->Fill(m1.eta(), PrescaleWeight);
              mu2Pt_.numerator->Fill(m1.pt(), PrescaleWeight);
              DiMuPt_.numerator->Fill((m1.p4() + m.p4()).Pt(), PrescaleWeight);
              DiMuEta_.numerator->Fill((m1.p4() + m.p4()).Eta(), PrescaleWeight);
              DiMuPhi_.numerator->Fill((m1.p4() + m.p4()).Phi(), PrescaleWeight);
              DiMudR_.numerator->Fill(reco::deltaR(m, m1), PrescaleWeight);
            }

            break;

          case 6:

            if (dimuonCL < minprob)
              continue;
            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;

            for (auto const& m2 : *muoHandle) {
              if (m2.pt() == m.pt())
                continue;  // remove duplicates but do not introduce ordering
              if (m2.pt() == m1.pt())
                continue;  // -> m2 will be the non-resonant and non-vertexing muon in the triplet
              if (!matchToTrigger(hltpath, m2))
                continue;

              mu1Phi_.denominator->Fill(m.phi());
              mu1Eta_.denominator->Fill(m.eta());
              mu1Pt_.denominator->Fill(m.pt());
              mu2Phi_.denominator->Fill(m1.phi());
              mu2Eta_.denominator->Fill(m1.eta());
              mu2Pt_.denominator->Fill(m1.pt());
              mu3Phi_.denominator->Fill(m2.phi());
              mu3Eta_.denominator->Fill(m2.eta());
              mu3Pt_.denominator->Fill(m2.pt());

              if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
                if (!matchToTrigger(hltpath1, m1))
                  continue;
                if (!matchToTrigger(hltpath1, m))
                  continue;
                if (!matchToTrigger(hltpath1, m2))
                  continue;
                mu1Phi_.numerator->Fill(m.phi(), PrescaleWeight);
                mu1Eta_.numerator->Fill(m.eta(), PrescaleWeight);
                mu1Pt_.numerator->Fill(m.pt(), PrescaleWeight);
                mu2Phi_.numerator->Fill(m1.phi(), PrescaleWeight);
                mu2Eta_.numerator->Fill(m1.eta(), PrescaleWeight);
                mu2Pt_.numerator->Fill(m1.pt(), PrescaleWeight);
                mu3Phi_.numerator->Fill(m2.phi(), PrescaleWeight);
                mu3Eta_.numerator->Fill(m2.eta(), PrescaleWeight);
                mu3Pt_.numerator->Fill(m2.pt(), PrescaleWeight);
              }
            }

            break;

          case 7:

            if (phHandle.isValid()) {
              if (!phHandle->empty())
                for (auto const& p : *phHandle) {
                  if (!matchToTrigger(hltpath, p))
                    continue;

                  phPhi_.denominator->Fill(p.phi());
                  phEta_.denominator->Fill(p.eta());
                  phPt_.denominator->Fill(p.pt());

                  if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
                    if (!matchToTrigger(hltpath1, p))
                      continue;
                    if (!matchToTrigger(hltpath1, m))
                      continue;
                    if (!matchToTrigger(hltpath1, m1))
                      continue;
                    phPhi_.numerator->Fill(p.phi(), PrescaleWeight);
                    phEta_.numerator->Fill(p.eta(), PrescaleWeight);
                    phPt_.numerator->Fill(p.pt(), PrescaleWeight);
                  }
                }
            } else {
              if (!warningPrinted4token_[3]) {
                warningPrinted4token_[3] = true;
                if (phInputTag_.label().empty())
                  edm::LogWarning("BPHMonitor") << "PhotonCollection not set";
                else
                  edm::LogWarning("BPHMonitor") << "skipping events because the collection "
                                                << phInputTag_.label().c_str() << " is not available";
              }
              // if Handle is not valid, because the InputTag has been mis-configured, then skip the event
              if (!phInputTag_.label().empty())
                return;
            }

            break;

          case 8:  //vtx monitoring, filling probability, DS, DCA, cos of pointing angle to the PV, eta, pT of dimuon

            if (displaced_)
              if ((displacementFromBeamspotJpsi.perp() / sqrt(jerr.rerr(displacementFromBeamspotJpsi))) < minDS)
                continue;
            if ((Jpsi_) && (!Upsilon_))
              if (DiMuMass > maxmassJpsi || DiMuMass < minmassJpsi)
                continue;
            if ((!Jpsi_) && (Upsilon_))
              if (DiMuMass > maxmassUpsilon || DiMuMass < minmassUpsilon)
                continue;

            // fill vtx-prob plots before selecting on this variable
            DiMuProb_.denominator->Fill(dimuonCL);
            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup) &&
                matchToTrigger(hltpath1, m1) && matchToTrigger(hltpath1, m))
              DiMuProb_.numerator->Fill(dimuonCL, PrescaleWeight);

            if (dimuonCL < minprob)
              continue;

            DiMuDS_.denominator->Fill(displacementFromBeamspotJpsi.perp() /
                                      sqrt(jerr.rerr(displacementFromBeamspotJpsi)));
            DiMuPVcos_.denominator->Fill(jpsi_cos);
            DiMuPt_.denominator->Fill((m1.p4() + m.p4()).Pt());
            DiMuEta_.denominator->Fill((m1.p4() + m.p4()).Eta());
            DiMuPhi_.denominator->Fill((m1.p4() + m.p4()).Phi());
            DiMuDCA_.denominator->Fill(cApp.distance());

            if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
              if (!matchToTrigger(hltpath1, m1))
                continue;
              if (!matchToTrigger(hltpath1, m))
                continue;
              DiMuDS_.numerator->Fill(
                  displacementFromBeamspotJpsi.perp() / sqrt(jerr.rerr(displacementFromBeamspotJpsi)), PrescaleWeight);
              DiMuPVcos_.numerator->Fill(jpsi_cos, PrescaleWeight);
              DiMuPt_.numerator->Fill((m1.p4() + m.p4()).Pt(), PrescaleWeight);
              DiMuEta_.numerator->Fill((m1.p4() + m.p4()).Eta(), PrescaleWeight);
              DiMuPhi_.numerator->Fill((m1.p4() + m.p4()).Phi(), PrescaleWeight);
              DiMuDCA_.numerator->Fill(cApp.distance(), PrescaleWeight);
            }

            break;

          case 9:

            if (dimuonCL < minprob)
              continue;
            if (fabs(jpsi_cos) < mincos)
              continue;
            if ((displacementFromBeamspotJpsi.perp() / sqrt(jerr.rerr(displacementFromBeamspotJpsi))) < minDS)
              continue;

            if (trHandle.isValid())
              for (auto const& t : *trHandle) {
                if (!trSelection_ref(t))
                  continue;
                const reco::Track& itrk1 = t;
                if (reco::deltaR(t, m1) <= min_dR)
                  continue;
                if (reco::deltaR(t, m) <= min_dR)
                  continue;
                if (!itrk1.quality(reco::TrackBase::highPurity))
                  continue;

                // reconstruct B+ hadron
                reco::Particle::LorentzVector pB, p1, p2, p3;
                double trackMass2 = kaon_mass * kaon_mass;
                double MuMass2 = mu_mass * mu_mass;
                double e1 = sqrt(m.momentum().Mag2() + MuMass2);
                double e2 = sqrt(m1.momentum().Mag2() + MuMass2);
                double e3 = sqrt(itrk1.momentum().Mag2() + trackMass2);
                p1 = reco::Particle::LorentzVector(m.px(), m.py(), m.pz(), e1);
                p2 = reco::Particle::LorentzVector(m1.px(), m1.py(), m1.pz(), e2);
                p3 = reco::Particle::LorentzVector(itrk1.px(), itrk1.py(), itrk1.pz(), e3);
                pB = p1 + p2 + p3;
                if (pB.mass() > maxmassJpsiTk || pB.mass() < minmassJpsiTk)
                  continue;
                reco::TransientTrack trTT(itrk1, &(*bFieldHandle));
                std::vector<reco::TransientTrack> t_tks;
                t_tks.push_back(mu1TT);
                t_tks.push_back(mu2TT);
                t_tks.push_back(trTT);
                KalmanVertexFitter kvf;
                TransientVertex tv = kvf.vertex(t_tks);
                reco::Vertex vertex = tv;
                if (!tv.isValid())
                  continue;
                float JpsiTkCL = 0;
                if ((vertex.chi2() >= 0.0) && (vertex.ndof() > 0))
                  JpsiTkCL = TMath::Prob(vertex.chi2(), vertex.ndof());
                math::XYZVector pperp(m.px() + m1.px() + itrk1.px(), m.py() + m1.py() + itrk1.py(), 0.);
                GlobalPoint secondaryVertex = tv.position();
                GlobalError err = tv.positionError();
                GlobalPoint displacementFromBeamspot(
                    -1 * ((vertexBeamSpot.x0() - secondaryVertex.x()) +
                          (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dxdz()),
                    -1 * ((vertexBeamSpot.y0() - secondaryVertex.y()) +
                          (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dydz()),
                    0);
                reco::Vertex::Point vperp(displacementFromBeamspot.x(), displacementFromBeamspot.y(), 0.);
                float jpsiKcos = vperp.Dot(pperp) / (vperp.R() * pperp.R());
                if (JpsiTkCL < minprob)
                  continue;
                if (fabs(jpsiKcos) < mincos)
                  continue;
                if ((displacementFromBeamspot.perp() / sqrt(err.rerr(displacementFromBeamspot))) < minDS)
                  continue;

                muPhi_.denominator->Fill(t.phi());
                muEta_.denominator->Fill(t.eta());
                muPt_.denominator->Fill(t.pt());
                BMass_.denominator->Fill(pB.mass());

                if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
                  if (!matchToTrigger(hltpath1, m1))
                    continue;
                  if (!matchToTrigger(hltpath1, m))
                    continue;
                  if (!matchToTrigger(hltpath1, t))
                    continue;
                  muPhi_.numerator->Fill(t.phi(), PrescaleWeight);
                  muEta_.numerator->Fill(t.eta(), PrescaleWeight);
                  muPt_.numerator->Fill(t.pt(), PrescaleWeight);
                  BMass_.numerator->Fill(pB.mass(), PrescaleWeight);
                }
              }

            break;

          case 10:

            if (trHandle.isValid())
              for (auto const& t : *trHandle) {
                if (!trSelection_ref(t))
                  continue;
                const reco::Track& itrk1 = t;
                if (reco::deltaR(t, m1) <= min_dR)
                  continue;
                if (reco::deltaR(t, m) <= min_dR)
                  continue;
                if (!itrk1.quality(reco::TrackBase::highPurity))
                  continue;

                // reconstruct Mu+TkMu structure
                reco::Particle::LorentzVector pB, p2, p3;
                double trackMass2 = kaon_mass * kaon_mass;
                double MuMass2 = mu_mass * mu_mass;
                double e2 = sqrt(m1.momentum().Mag2() + MuMass2);
                double e3 = sqrt(itrk1.momentum().Mag2() + trackMass2);
                p2 = reco::Particle::LorentzVector(m1.px(), m1.py(), m1.pz(), e2);
                p3 = reco::Particle::LorentzVector(itrk1.px(), itrk1.py(), itrk1.pz(), e3);
                pB = p2 + p3;
                if (pB.mass() > maxmassJpsiTk || pB.mass() < minmassJpsiTk)
                  continue;
                reco::TransientTrack trTT(itrk1, &(*bFieldHandle));
                std::vector<reco::TransientTrack> t_tks;
                t_tks.push_back(mu2TT);
                t_tks.push_back(trTT);
                KalmanVertexFitter kvf;
                TransientVertex tv = kvf.vertex(t_tks);
                reco::Vertex vertex = tv;
                if (!tv.isValid())
                  continue;
                float JpsiTkCL = 0;
                if ((vertex.chi2() >= 0.0) && (vertex.ndof() > 0))
                  JpsiTkCL = TMath::Prob(vertex.chi2(), vertex.ndof());
                math::XYZVector pperp(m1.px() + itrk1.px(), m1.py() + itrk1.py(), 0.);
                GlobalPoint secondaryVertex = tv.position();
                GlobalError err = tv.positionError();
                GlobalPoint displacementFromBeamspot(
                    -1 * ((vertexBeamSpot.x0() - secondaryVertex.x()) +
                          (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dxdz()),
                    -1 * ((vertexBeamSpot.y0() - secondaryVertex.y()) +
                          (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dydz()),
                    0);
                reco::Vertex::Point vperp(displacementFromBeamspot.x(), displacementFromBeamspot.y(), 0.);
                if (JpsiTkCL < minprob)
                  continue;

                muPhi_.denominator->Fill(m1.phi());
                muEta_.denominator->Fill(m1.eta());
                muPt_.denominator->Fill(m1.pt());

                if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
                  if (!matchToTrigger(hltpath1, m1))
                    continue;
                  if (!matchToTrigger(hltpath1, m))
                    continue;
                  if (!matchToTrigger(hltpath1, t))
                    continue;
                  muPhi_.numerator->Fill(m1.phi(), PrescaleWeight);
                  muEta_.numerator->Fill(m1.eta(), PrescaleWeight);
                  muPt_.numerator->Fill(m1.pt(), PrescaleWeight);
                }
              }

            break;

          case 11:

            if (dimuonCL < minprob)
              continue;
            if (fabs(jpsi_cos) < mincos)
              continue;
            if (displacementFromBeamspotJpsi.perp() / sqrt(jerr.rerr(displacementFromBeamspotJpsi)) < minDS)
              continue;

            if (trHandle.isValid())
              for (auto const& t : *trHandle) {
                if (!trSelection_ref(t))
                  continue;
                if ((reco::deltaR(t, m) <= min_dR))
                  continue;
                if ((reco::deltaR(t, m1) <= min_dR))
                  continue;

                for (auto const& t1 : *trHandle) {
                  if (&t - &(*trHandle)[0] >= &t1 - &(*trHandle)[0])
                    continue;  // not enough, need the following DeltaR checks
                  if (!trSelection_ref(t1))
                    continue;
                  if ((reco::deltaR(t1, m) <= min_dR))
                    continue;
                  if ((reco::deltaR(t1, m1) <= min_dR))
                    continue;
                  if ((reco::deltaR(t, t1) <= min_dR))
                    continue;
                  const reco::Track& itrk1 = t;
                  const reco::Track& itrk2 = t1;
                  if (!itrk1.quality(reco::TrackBase::highPurity))
                    continue;
                  if (!itrk2.quality(reco::TrackBase::highPurity))
                    continue;

                  // reconstruct Bs candidate
                  reco::Particle::LorentzVector pB, pTkTk, p1, p2, p3, p4;
                  double trackMass2 = kaon_mass * kaon_mass;
                  double MuMass2 = mu_mass * mu_mass;
                  double e1 = sqrt(m.momentum().Mag2() + MuMass2);
                  double e2 = sqrt(m1.momentum().Mag2() + MuMass2);
                  double e3 = sqrt(itrk1.momentum().Mag2() + trackMass2);
                  double e4 = sqrt(itrk2.momentum().Mag2() + trackMass2);
                  p1 = reco::Particle::LorentzVector(m.px(), m.py(), m.pz(), e1);
                  p2 = reco::Particle::LorentzVector(m1.px(), m1.py(), m1.pz(), e2);
                  p3 = reco::Particle::LorentzVector(itrk1.px(), itrk1.py(), itrk1.pz(), e3);
                  p4 = reco::Particle::LorentzVector(itrk2.px(), itrk2.py(), itrk2.pz(), e4);
                  pTkTk = p3 + p4;
                  if (pTkTk.mass() > maxmassTkTk || pTkTk.mass() < minmassTkTk)
                    continue;
                  pB = p1 + p2 + p3 + p4;
                  if (pB.mass() > maxmassJpsiTk || pB.mass() < minmassJpsiTk)
                    continue;
                  reco::TransientTrack mu1TT(m.track(), &(*bFieldHandle));
                  reco::TransientTrack mu2TT(m1.track(), &(*bFieldHandle));
                  reco::TransientTrack trTT(itrk1, &(*bFieldHandle));
                  reco::TransientTrack tr1TT(itrk2, &(*bFieldHandle));
                  std::vector<reco::TransientTrack> t_tks;
                  t_tks.push_back(mu1TT);
                  t_tks.push_back(mu2TT);
                  t_tks.push_back(trTT);
                  t_tks.push_back(tr1TT);
                  KalmanVertexFitter kvf;
                  TransientVertex tv = kvf.vertex(t_tks);  // this will compare the tracks
                  reco::Vertex vertex = tv;
                  if (!tv.isValid())
                    continue;
                  float JpsiTkCL = 0;
                  if ((vertex.chi2() >= 0.0) && (vertex.ndof() > 0))
                    JpsiTkCL = TMath::Prob(vertex.chi2(), vertex.ndof());
                  math::XYZVector pperp(
                      m.px() + m1.px() + itrk1.px() + itrk2.px(), m.py() + m1.py() + itrk1.py() + itrk2.py(), 0.);
                  GlobalPoint secondaryVertex = tv.position();
                  GlobalError err = tv.positionError();
                  GlobalPoint displacementFromBeamspot(
                      -1 * ((vertexBeamSpot.x0() - secondaryVertex.x()) +
                            (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dxdz()),
                      -1 * ((vertexBeamSpot.y0() - secondaryVertex.y()) +
                            (secondaryVertex.z() - vertexBeamSpot.z0()) * vertexBeamSpot.dydz()),
                      0);
                  reco::Vertex::Point vperp(displacementFromBeamspot.x(), displacementFromBeamspot.y(), 0.);
                  float jpsiKcos = vperp.Dot(pperp) / (vperp.R() * pperp.R());
                  if (JpsiTkCL < minprob)
                    continue;
                  if (fabs(jpsiKcos) < mincos)
                    continue;
                  if ((displacementFromBeamspot.perp() / sqrt(err.rerr(displacementFromBeamspot))) < minDS)
                    continue;

                  mu1Phi_.denominator->Fill(t.phi());
                  mu1Eta_.denominator->Fill(t.eta());
                  mu1Pt_.denominator->Fill(t.pt());
                  mu2Phi_.denominator->Fill(t1.phi());
                  mu2Eta_.denominator->Fill(t1.eta());
                  mu2Pt_.denominator->Fill(t1.pt());
                  BMass_.denominator->Fill(pB.mass());

                  if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(iEvent, iSetup)) {
                    if (!matchToTrigger(hltpath1, m))
                      continue;
                    if (!matchToTrigger(hltpath1, m1))
                      continue;
                    if (!matchToTrigger(hltpath1, t))
                      continue;
                    if (!matchToTrigger(hltpath1, t1))
                      continue;
                    mu1Phi_.numerator->Fill(t.phi(), PrescaleWeight);
                    mu1Eta_.numerator->Fill(t.eta(), PrescaleWeight);
                    mu1Pt_.numerator->Fill(t.pt(), PrescaleWeight);
                    mu2Phi_.numerator->Fill(t1.phi(), PrescaleWeight);
                    mu2Eta_.numerator->Fill(t1.eta(), PrescaleWeight);
                    mu2Pt_.numerator->Fill(t1.pt(), PrescaleWeight);
                    BMass_.numerator->Fill(pB.mass(), PrescaleWeight);
                  }

                }  // for (auto const & t1 : *trHandle)
              }    // for (auto const & t : *trHandle)

            break;
        }
      }
    }
  }
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 279 of file BPHMonitor.cc.

References BMass_, Bmass_binning_, bookME(), cos_binning_, dca_binning_, den_genTriggerEventFlag_, DiMuDCA_, DiMudR_, DiMuDS_, DiMuEta_, DiMuMass_, DiMuPhi_, DiMuProb_, DiMuPt_, DiMuPVcos_, dMu_pt_binning_, dR_binning_, ds_binning_, enum_, eta_binning_, folderName_, hltConfig_, HLTPrescaleProvider::hltConfigProvider(), hltPrescale_, HLTPrescaleProvider::init(), GenericTriggerEventFlag::initRun(), mass_binning_, mu1Eta_, mu1Phi_, mu1Pt_, mu2Eta_, mu2Phi_, mu2Pt_, mu3Eta_, mu3Phi_, mu3Pt_, muEta_, muPhi_, muPt_, num_genTriggerEventFlag_, GenericTriggerEventFlag::on(), phi_binning_, prob_binning_, pt_binning_, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), setMETitle(), AlCaHLTBitMon_QueryRunRegistry::string, tnp_, and trOrMu_.

                                                                                                         {
  std::string histname, histtitle, istnp, trMuPh;
  bool Ph_ = false;
  if (enum_ == 7)
    Ph_ = true;
  if (tnp_)
    istnp = "Tag_and_Probe/";
  else
    istnp = "";
  std::string currentFolder = folderName_ + istnp;
  ibooker.setCurrentFolder(currentFolder);
  if (trOrMu_)
    trMuPh = "tr";
  else if (Ph_)
    trMuPh = "ph";
  else
    trMuPh = "mu";

  if (enum_ == 7 || enum_ == 1 || enum_ == 9 || enum_ == 10) {
    histname = trMuPh + "Pt";
    histtitle = trMuPh + "_P_{t}";
    bookME(ibooker, muPt_, histname, histtitle, pt_binning_);
    setMETitle(muPt_, trMuPh + "_Pt[GeV]", "events / 1 GeV");

    histname = trMuPh + "Phi";
    histtitle = trMuPh + "Phi";
    bookME(ibooker, muPhi_, histname, histtitle, phi_binning_);
    setMETitle(muPhi_, trMuPh + "_#phi", "events / 0.1 rad");

    histname = trMuPh + "Eta";
    histtitle = trMuPh + "_Eta";
    bookME(ibooker, muEta_, histname, histtitle, eta_binning_);
    setMETitle(muEta_, trMuPh + "_#eta", "events / 0.2");

    if (enum_ == 9) {
      histname = "BMass";
      histtitle = "BMass";
      bookME(ibooker, BMass_, histname, histtitle, Bmass_binning_);
      setMETitle(BMass_, "B_#mass", "events /");
    }
  } else {
    if (enum_ != 8) {
      histname = trMuPh + "1Pt";
      histtitle = trMuPh + "1_P_{t}";
      bookME(ibooker, mu1Pt_, histname, histtitle, pt_binning_);
      setMETitle(mu1Pt_, trMuPh + "_Pt[GeV]", "events / 1 GeV");

      histname = trMuPh + "1Phi";
      histtitle = trMuPh + "1Phi";
      bookME(ibooker, mu1Phi_, histname, histtitle, phi_binning_);
      setMETitle(mu1Phi_, trMuPh + "_#phi", "events / 0.1 rad");

      histname = trMuPh + "1Eta";
      histtitle = trMuPh + "1_Eta";
      bookME(ibooker, mu1Eta_, histname, histtitle, eta_binning_);
      setMETitle(mu1Eta_, trMuPh + "_#eta", "events / 0.2");

      histname = trMuPh + "2Pt";
      histtitle = trMuPh + "2_P_{t}";
      bookME(ibooker, mu2Pt_, histname, histtitle, pt_binning_);
      setMETitle(mu2Pt_, trMuPh + "_Pt[GeV]", "events / 1 GeV");

      histname = trMuPh + "2Phi";
      histtitle = trMuPh + "2Phi";
      bookME(ibooker, mu2Phi_, histname, histtitle, phi_binning_);
      setMETitle(mu2Phi_, trMuPh + "_#phi", "events / 0.1 rad");

      histname = trMuPh + "2Eta";
      histtitle = trMuPh + "2_Eta";
      bookME(ibooker, mu2Eta_, histname, histtitle, eta_binning_);
      setMETitle(mu2Eta_, trMuPh + "_#eta", "events / 0.2");
      if (enum_ == 11) {
        histname = "BMass";
        histtitle = "BMass";
        bookME(ibooker, BMass_, histname, histtitle, Bmass_binning_);
        setMETitle(BMass_, "B_#mass", "events /");
      }
    }
    if (enum_ == 6) {
      histname = trMuPh + "3Eta";
      histtitle = trMuPh + "3Eta";
      bookME(ibooker, mu3Eta_, histname, histtitle, eta_binning_);
      setMETitle(mu3Eta_, trMuPh + "3#eta", "events / 0.2");

      histname = trMuPh + "3Pt";
      histtitle = trMuPh + "3_P_{t}";
      bookME(ibooker, mu3Pt_, histname, histtitle, pt_binning_);
      setMETitle(mu3Pt_, trMuPh + "3_Pt[GeV]", "events / 1 GeV");

      histname = trMuPh + "3Phi";
      histtitle = trMuPh + "3Phi";
      bookME(ibooker, mu3Phi_, histname, histtitle, phi_binning_);
      setMETitle(mu3Phi_, trMuPh + "3_#phi", "events / 0.1 rad");

    } else if (enum_ == 2 || enum_ == 4 || enum_ == 5 || enum_ == 8) {
      histname = "DiMuEta";
      histtitle = "DiMuEta";
      bookME(ibooker, DiMuEta_, histname, histtitle, eta_binning_);
      setMETitle(DiMuEta_, "DiMu#eta", "events / 0.2");

      histname = "DiMuPt";
      histtitle = "DiMu_P_{t}";
      bookME(ibooker, DiMuPt_, histname, histtitle, dMu_pt_binning_);
      setMETitle(DiMuPt_, "DiMu_Pt[GeV]", "events / 1 GeV");

      histname = "DiMuPhi";
      histtitle = "DiMuPhi";
      bookME(ibooker, DiMuPhi_, histname, histtitle, phi_binning_);
      setMETitle(DiMuPhi_, "DiMu_#phi", "events / 0.1 rad");

      if (enum_ == 4 || enum_ == 5) {
        histname = "DiMudR";
        histtitle = "DiMudR";
        bookME(ibooker, DiMudR_, histname, histtitle, dR_binning_);
        setMETitle(DiMudR_, "DiMu_#dR", "events /");

        if (enum_ == 4) {
          histname = "DiMuMass";
          histtitle = "DiMuMass";
          bookME(ibooker, DiMuMass_, histname, histtitle, mass_binning_);
          setMETitle(DiMuMass_, "DiMu_#mass", "events /");
        }
      } else if (enum_ == 8) {
        histname = "DiMuProb";
        histtitle = "DiMuProb";
        bookME(ibooker, DiMuProb_, histname, histtitle, prob_binning_);
        setMETitle(DiMuProb_, "DiMu_#prob", "events /");

        histname = "DiMuPVcos";
        histtitle = "DiMuPVcos";
        bookME(ibooker, DiMuPVcos_, histname, histtitle, cos_binning_);
        setMETitle(DiMuPVcos_, "DiMu_#cosPV", "events /");

        histname = "DiMuDS";
        histtitle = "DiMuDS";
        bookME(ibooker, DiMuDS_, histname, histtitle, ds_binning_);
        setMETitle(DiMuDS_, "DiMu_#ds", "events /");

        histname = "DiMuDCA";
        histtitle = "DiMuDCA";
        bookME(ibooker, DiMuDCA_, histname, histtitle, dca_binning_);
        setMETitle(DiMuDCA_, "DiMu_#dca", "events /");
      }
    }  // if (enum_ == 2 || enum_ == 4 || enum_ == 5 || enum_ == 8)
  }

  // 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 changed = true;

  hltPrescale_->init(iRun, iSetup, "HLT", changed);
  hltConfig_ = hltPrescale_->hltConfigProvider();
}

void BPHMonitor::bookME ( DQMStore::IBooker &  ibooker, METME &  me, std::string &  histname, std::string &  histtitle, int &  nbins, double &  xmin, double &  xmax  )

protected

Definition at line 201 of file BPHMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), METME::denominator, and METME::numerator.

Referenced by bookHistograms(), and bookME().

                                     {
  me.numerator = ibooker.book1D(histname + "_numerator", histtitle + " (numerator)", nbins, min, max);
  me.denominator = ibooker.book1D(histname + "_denominator", histtitle + " (denominator)", nbins, min, max);
}

void BPHMonitor::bookME ( DQMStore::IBooker &  ibooker, METME &  me, std::string &  histname, std::string &  histtitle, std::vector< double >  binningX  )

protected

Definition at line 211 of file BPHMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), METME::denominator, LaserClient_cfi::nbins, and METME::numerator.

                                                                                                               {
  int nbins = binning.size() - 1;
  std::vector<float> fbinning(binning.begin(), binning.end());
  float* arr = &fbinning[0];
  me.numerator = ibooker.book1D(histname + "_numerator", histtitle + " (numerator)", nbins, arr);
  me.denominator = ibooker.book1D(histname + "_denominator", histtitle + " (denominator)", nbins, arr);
}

void BPHMonitor::bookME ( DQMStore::IBooker &  ibooker, METME &  me, std::string &  histname, std::string &  histtitle, int &  nbinsX, double &  xmin, double &  xmax, double &  ymin, double &  ymax  )

protected

Definition at line 219 of file BPHMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), METME::denominator, and METME::numerator.

                                      {
  me.numerator =
      ibooker.bookProfile(histname + "_numerator", histtitle + " (numerator)", nbinsX, xmin, xmax, ymin, ymax);
  me.denominator =
      ibooker.bookProfile(histname + "_denominator", histtitle + " (denominator)", nbinsX, xmin, xmax, ymin, ymax);
}

void BPHMonitor::bookME ( DQMStore::IBooker &  ibooker, METME &  me, std::string &  histname, std::string &  histtitle, int &  nbinsX, double &  xmin, double &  xmax, int &  nbinsY, double &  ymin, double &  ymax  )

protected

Definition at line 233 of file BPHMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), METME::denominator, and METME::numerator.

                                      {
  me.numerator =
      ibooker.book2D(histname + "_numerator", histtitle + " (numerator)", nbinsX, xmin, xmax, nbinsY, ymin, ymax);
  me.denominator =
      ibooker.book2D(histname + "_denominator", histtitle + " (denominator)", nbinsX, xmin, xmax, nbinsY, ymin, ymax);
}

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

protected

Definition at line 248 of file BPHMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), METME::denominator, hlt_dqm_clientPB-live_cfg::nbinsX, hlt_dqm_clientPB-live_cfg::nbinsY, and METME::numerator.

                                                    {
  int nbinsX = binningX.size() - 1;
  std::vector<float> fbinningX(binningX.begin(), binningX.end());
  float* arrX = &fbinningX[0];
  int nbinsY = binningY.size() - 1;
  std::vector<float> fbinningY(binningY.begin(), binningY.end());
  float* arrY = &fbinningY[0];

  me.numerator = ibooker.book2D(histname + "_numerator", histtitle + " (numerator)", nbinsX, arrX, nbinsY, arrY);
  me.denominator = ibooker.book2D(histname + "_denominator", histtitle + " (denominator)", nbinsX, arrX, nbinsY, arrY);
}

void BPHMonitor::bookME ( DQMStore::IBooker &  ibooker, METME &  me, std::string &  histname, std::string &  histtitle, MEbinning &  binning  )

protected

Definition at line 265 of file BPHMonitor.cc.

References bookME().

                                              {
  // If the vector in the binning is filled use the bins defined there
  // otherwise use a linear binning between min and max
  if (binning.edges.empty()) {
    this->bookME(ibooker, me, histname, histtitle, binning.nbins, binning.xmin, binning.xmax);
  } else {
    this->bookME(ibooker, me, histname, histtitle, binning.edges);
  }
}

void BPHMonitor::case11_selection	(	const float &	dimuonCL,
		const float &	jpsi_cos,
		const GlobalPoint &	displacementFromBeamspotJpsi,
		const GlobalError &	jerr,
		const edm::Handle< reco::TrackCollection > &	trHandle,
		const std::string &	hltpath,
		const edm::Handle< trigger::TriggerEvent > &	handleTriggerEvent,
		const reco::Muon &	m,
		const reco::Muon &	m1,
		const edm::ESHandle< MagneticField > &	bFieldHandle,
		const reco::BeamSpot &	vertexBeamSpot,
		MonitorElement *	phi1,
		MonitorElement *	eta1,
		MonitorElement *	pT1,
		MonitorElement *	phi2,
		MonitorElement *	eta2,
		MonitorElement *	pT2
	)

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

static

Definition at line 1186 of file BPHMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), BPHMonitor_cfi::BmassPSet, BPHMonitor_cfi::cosPSet, BPHMonitor_cfi::d0PSet, BPHMonitor_cfi::dcaPSet, BPHMonitor_cfi::dMu_ptPSet, BPHMonitor_cff::dRPSet, BPHMonitor_cfi::dsPSet, V0Monitor_cfi::etaPSet, fillHistoPSetDescription(), SiStripSourceConfigTier0_cff::genericTriggerEventPSet, LumiMonitor_cff::histoPSet, HLT_2018_cff::InputTag, V0Monitor_cff::massPSet, BPHMonitor_cff::phiPSet, BPHMonitor_cfi::probPSet, V0Monitor_cfi::ptPSet, AlCaHLTBitMon_QueryRunRegistry::string, and BPHMonitor_cfi::z0PSet.

                                                                            {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("FolderName", "HLT/BPH/");
  desc.add<edm::InputTag>("tracks", edm::InputTag("generalTracks"));
  desc.add<edm::InputTag>("photons", edm::InputTag("photons"));
  desc.add<edm::InputTag>("offlinePVs", edm::InputTag("offlinePrimaryVertices"));
  desc.add<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot"));
  desc.add<edm::InputTag>("muons", edm::InputTag("muons"));
  desc.add<edm::InputTag>("hltTriggerSummaryAOD", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
  desc.add<std::string>("muoSelection", "");
  desc.add<std::string>("muoSelection_ref",
                        "isPFMuon & isGlobalMuon  & innerTrack.hitPattern.trackerLayersWithMeasurement>5 & "
                        "innerTrack.hitPattern.numberOfValidPixelHits> 0");
  desc.add<std::string>(
      "muoSelection_tag",
      "isGlobalMuon && isPFMuon && isTrackerMuon && abs(eta) < 2.4 && innerTrack.hitPattern.numberOfValidPixelHits > 0 "
      "&& innerTrack.hitPattern.trackerLayersWithMeasurement > 5 && globalTrack.hitPattern.numberOfValidMuonHits > 0 "
      "&& globalTrack.normalizedChi2 < 10");  // tight selection for tag muon
  desc.add<std::string>("muoSelection_probe",
                        "isPFMuon & isGlobalMuon  & innerTrack.hitPattern.trackerLayersWithMeasurement>5 & "
                        "innerTrack.hitPattern.numberOfValidPixelHits> 0");
  desc.add<std::string>("trSelection_ref", "");
  desc.add<std::string>("DMSelection_ref", "Pt>4 & abs(eta)");

  desc.add<int>("nmuons", 1);
  desc.add<bool>("tnp", false);
  desc.add<int>("L3", 0);
  desc.add<int>("ptCut", 0);
  desc.add<int>("displaced", 0);
  desc.add<int>("trOrMu", 0);  // if =0, track param monitoring
  desc.add<int>("Jpsi", 0);
  desc.add<int>("Upsilon", 0);
  desc.add<int>("enum", 1);  // 1...9, 9 sets of variables to be filled, depends on the hlt path
  desc.add<int>("seagull", 1);
  desc.add<double>("maxmass", 3.596);
  desc.add<double>("minmass", 2.596);
  desc.add<double>("maxmassJpsi", 3.2);
  desc.add<double>("minmassJpsi", 3.);
  desc.add<double>("maxmassUpsilon", 10.0);
  desc.add<double>("minmassUpsilon", 8.8);
  desc.add<double>("maxmassTkTk", 10);
  desc.add<double>("minmassTkTk", 0);
  desc.add<double>("maxmassJpsiTk", 5.46);
  desc.add<double>("minmassJpsiTk", 5.1);
  desc.add<double>("kaon_mass", 0.493677);
  desc.add<double>("mu_mass", 0.1056583745);
  desc.add<double>("min_dR", 0.001);
  desc.add<double>("max_dR", 1.4);
  desc.add<double>("minprob", 0.005);
  desc.add<double>("mincos", 0.95);
  desc.add<double>("minDS", 3.);

  edm::ParameterSetDescription genericTriggerEventPSet;
  genericTriggerEventPSet.add<bool>("andOr");
  genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
  genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
  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<bool>("andOrL1", true);
  genericTriggerEventPSet.add<std::vector<std::string>>("hltPaths", {});
  genericTriggerEventPSet.add<std::vector<std::string>>("l1Algorithms", {});
  genericTriggerEventPSet.add<std::string>("hltDBKey", "");
  genericTriggerEventPSet.add<bool>("errorReplyHlt", false);
  genericTriggerEventPSet.add<bool>("errorReplyL1", true);
  genericTriggerEventPSet.add<bool>("l1BeforeMask", true);
  genericTriggerEventPSet.add<unsigned int>("verbosityLevel", 0);
  desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
  desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);

  edm::ParameterSetDescription histoPSet;
  edm::ParameterSetDescription phiPSet;
  edm::ParameterSetDescription etaPSet;
  edm::ParameterSetDescription ptPSet;
  edm::ParameterSetDescription dMu_ptPSet;
  edm::ParameterSetDescription d0PSet;
  edm::ParameterSetDescription z0PSet;
  edm::ParameterSetDescription dRPSet;
  edm::ParameterSetDescription massPSet;
  edm::ParameterSetDescription BmassPSet;
  edm::ParameterSetDescription dcaPSet;
  edm::ParameterSetDescription dsPSet;
  edm::ParameterSetDescription cosPSet;
  edm::ParameterSetDescription probPSet;
  edm::ParameterSetDescription TCoPSet;
  edm::ParameterSetDescription PUPSet;
  fillHistoPSetDescription(phiPSet);
  fillHistoPSetDescription(ptPSet);
  fillHistoPSetDescription(dMu_ptPSet);
  fillHistoPSetDescription(etaPSet);
  fillHistoPSetDescription(z0PSet);
  fillHistoPSetDescription(d0PSet);
  fillHistoPSetDescription(dRPSet);
  fillHistoPSetDescription(massPSet);
  fillHistoPSetDescription(BmassPSet);
  fillHistoPSetDescription(dcaPSet);
  fillHistoPSetDescription(dsPSet);
  fillHistoPSetDescription(cosPSet);
  fillHistoPSetDescription(probPSet);
  histoPSet.add<edm::ParameterSetDescription>("d0PSet", d0PSet);
  histoPSet.add<edm::ParameterSetDescription>("etaPSet", etaPSet);
  histoPSet.add<edm::ParameterSetDescription>("phiPSet", phiPSet);
  histoPSet.add<edm::ParameterSetDescription>("ptPSet", ptPSet);
  histoPSet.add<edm::ParameterSetDescription>("dMu_ptPSet", dMu_ptPSet);
  histoPSet.add<edm::ParameterSetDescription>("z0PSet", z0PSet);
  histoPSet.add<edm::ParameterSetDescription>("dRPSet", dRPSet);
  histoPSet.add<edm::ParameterSetDescription>("massPSet", massPSet);
  histoPSet.add<edm::ParameterSetDescription>("BmassPSet", BmassPSet);
  histoPSet.add<edm::ParameterSetDescription>("dcaPSet", dcaPSet);
  histoPSet.add<edm::ParameterSetDescription>("dsPSet", dsPSet);
  histoPSet.add<edm::ParameterSetDescription>("cosPSet", cosPSet);
  histoPSet.add<edm::ParameterSetDescription>("probPSet", probPSet);
  desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
  descriptions.add("bphMonitoring", desc);
}

void BPHMonitor::fillHistoLSPSetDescription ( edm::ParameterSetDescription &  pset )

static

Definition at line 1184 of file BPHMonitor.cc.

References edm::ParameterSetDescription::add().

{ pset.add<int>("nbins", 2500); }

void BPHMonitor::fillHistoPSetDescription ( edm::ParameterSetDescription &  pset )

static

Definition at line 1178 of file BPHMonitor.cc.

References edm::ParameterSetDescription::addNode().

Referenced by fillDescriptions().

                                                                          {
  pset.addNode((edm::ParameterDescription<int>("nbins", true) and edm::ParameterDescription<double>("xmin", true) and
                edm::ParameterDescription<double>("xmax", true)) xor
               edm::ParameterDescription<std::vector<double>>("edges", true));
}

MEbinning BPHMonitor::getHistoPSet ( edm::ParameterSet  pset )

staticprivate

Definition at line 171 of file BPHMonitor.cc.

References edm::ParameterSet::existsAs(), and edm::ParameterSet::getParameter().

                                                       {
  // Due to the setup of the fillDescription only one of the
  // two cases is possible at this point.
  if (pset.existsAs<std::vector<double>>("edges")) {
    return MEbinning{pset.getParameter<std::vector<double>>("edges")};
  }

  return MEbinning{
      pset.getParameter<int32_t>("nbins"),
      pset.getParameter<double>("xmin"),
      pset.getParameter<double>("xmax"),
  };
}

std::string BPHMonitor::getTriggerName ( std::string  partialName )

private

Definition at line 1304 of file BPHMonitor.cc.

References hltConfig_, mps_fire::i, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerName().

Referenced by analyze(), and matchToTrigger().

                                                          {
  const std::string trigger_name_tmp = partialName.substr(0, partialName.find("v*"));
  const unsigned int Ntriggers(hltConfig_.size());
  std::string trigger_name = "";
  for (unsigned int i = 0; i < Ntriggers; i++) {
    trigger_name = hltConfig_.triggerName(i);
    if (trigger_name.find(trigger_name_tmp) != std::string::npos)
      break;
  }

  return trigger_name;
}

template<typename T >

bool BPHMonitor::matchToTrigger ( const std::string &  theTriggerName, T  t  )

protected

Definition at line 1318 of file BPHMonitor.cc.

References Abs(), reco::deltaR(), trigger::TriggerObject::eta(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), getTriggerName(), handleTriggerEvent, hltConfig_, hltInputTag_1, mps_fire::i, HLTConfigProvider::inited(), muonTagProbeFilters_cff::matched, HLTConfigProvider::moduleLabels(), trigger::TriggerObject::phi(), edm::InputTag::process(), trigger::TriggerObject::pt(), HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerIndex().

Referenced by analyze().

                                                                    {
  bool matched = false;
  //validity check
  if (!hltConfig_.inited())
    return false;

  //Find the precise trigger name
  std::string trigger_name = getTriggerName(theTriggerName);
  const unsigned int trigger_index = hltConfig_.triggerIndex(trigger_name);

  //loop over all the modules for this trigger
  //by default use the last one
  unsigned int Nmodules = hltConfig_.size(trigger_index);
  const vector<string>& moduleLabels(hltConfig_.moduleLabels(trigger_index));
  unsigned int fIdx = 0;
  for (unsigned int i = 0; i < Nmodules; i++) {
    const unsigned int tmp_fIdx =
        handleTriggerEvent->filterIndex(edm::InputTag(moduleLabels[i], "", hltInputTag_1.process()));
    if (tmp_fIdx < handleTriggerEvent->sizeFilters())  //index of not used filters are set to sizeFilters()
    {
      fIdx = tmp_fIdx;
    }  //good index
  }

  //loop over all the objects in the filter of choice
  const trigger::Keys& KEYS(handleTriggerEvent->filterKeys(fIdx));
  const trigger::size_type nK(KEYS.size());
  const trigger::TriggerObjectCollection& TOC(handleTriggerEvent->getObjects());
  for (trigger::size_type i = 0; i != nK; ++i) {
    const trigger::TriggerObject& TO(TOC[KEYS[i]]);
    //perform matching: deltaR and pt check
    if ((reco::deltaR(t.eta(), t.phi(), TO.eta(), TO.phi()) <= 0.2) && (TMath::Abs(t.pt() - TO.pt()) < 0.12)) {
      matched = true;
    }
  }

  return matched;
}

double BPHMonitor::Prescale ( const std::string  num, const std::string  den, edm::Event const &  iEvent, edm::EventSetup const &  iSetup, HLTPrescaleProvider *  hltPrescale_  )

protected

Definition at line 1357 of file BPHMonitor.cc.

References DEFINE_FWK_MODULE, dqmdumpme::first, RemoveAddSevLevel::flag, mps_fire::i, HLTPrescaleProvider::prescaleValuesInDetail(), edm::second(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

                                                               {
  int PrescaleHLT_num = 1;
  int PrescaleHLT_den = 1;
  double Prescale_num = 1;
  double L1P = 1, HLTP = 1;
  bool flag = true;
  std::vector<bool> theSame_den;
  std::vector<bool> theSame_num;
  //retrieving HLT prescale
  PrescaleHLT_den = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).second;
  PrescaleHLT_num = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).second;
  if (PrescaleHLT_den > 0 && PrescaleHLT_num > 0)
    HLTP = PrescaleHLT_num / std::__gcd(PrescaleHLT_num, PrescaleHLT_den);

  //retrieving L1 prescale
  //Checking if we have the same l1 seeds in den and num
  //taking into account that they can be written in different order in num and den
  //and some of them can be also switched off

  //check if for each den l1 there is the same l1 seed in num
  if (!(hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.empty()) {
    for (size_t iSeed = 0; iSeed < (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.size();
         ++iSeed) {
      std::string l1_den = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.at(iSeed).first;
      int l1_denp = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.at(iSeed).second;
      if (l1_denp < 1)
        continue;
      flag = false;
      for (size_t iSeed1 = 0; iSeed1 < (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.size();
           ++iSeed1) {
        std::string l1_num = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.at(iSeed1).first;
        int l1_nump = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.at(iSeed1).second;
        if (l1_num == l1_den && l1_nump >= 1)  //the same seed
        {
          flag = true;
          break;
        }
      }
      theSame_den.push_back(flag);
    }
  }
  //check if for each num l1 there is the same l1 seed in den
  if (!(hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.empty()) {
    for (size_t iSeed = 0; iSeed < (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.size();
         ++iSeed) {
      std::string l1_num = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.at(iSeed).first;
      int l1_nump = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.at(iSeed).second;
      if (l1_nump < 1)
        continue;
      flag = false;
      for (size_t iSeed1 = 0; iSeed1 < (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.size();
           ++iSeed1) {
        std::string l1_den = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.at(iSeed1).first;
        int l1_denp = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath)).first.at(iSeed1).second;
        if (l1_den == l1_num && l1_denp >= 1)  //the same seed
        {
          flag = true;
          break;
        }
      }
      theSame_num.push_back(flag);
    }
  }
  flag = true;

  if (theSame_num.size() == theSame_den.size()) {
    for (size_t i = 0; i < theSame_num.size(); ++i) {
      if ((!theSame_num.at(i)) || (!theSame_den.at(i))) {
        flag = false;
        break;
      }
    }
  }

  if (flag && (theSame_num.size() == theSame_den.size())) {
    L1P = 1;  //den and num have the same set of l1 seeds
  } else {
    if (!(hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.empty()) {
      Prescale_num = 1;
      for (size_t iSeed = 0; iSeed < (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.size();
           ++iSeed) {
        int l1 = (hltPrescale_->prescaleValuesInDetail(iEvent, iSetup, hltpath1)).first.at(iSeed).second;
        if (l1 < 1)
          continue;
        if (l1 == 1) {
          Prescale_num = 1;
          break;
        } else
          Prescale_num *= 1 - (1.0 / (l1));
      }
      if (Prescale_num != 1)
        Prescale_num = 1.0 / (1 - Prescale_num);
    }
    L1P = Prescale_num;
  }

  return L1P * HLTP;
}

void BPHMonitor::setMETitle ( METME &  me, std::string  titleX, std::string  titleY  )

protected

Definition at line 194 of file BPHMonitor.cc.

References METME::denominator, METME::numerator, and dqm::impl::MonitorElement::setAxisTitle().

Referenced by bookHistograms().

                                                                         {
  me.numerator->setAxisTitle(titleX, 1);
  me.numerator->setAxisTitle(titleY, 2);
  me.denominator->setAxisTitle(titleX, 1);
  me.denominator->setAxisTitle(titleY, 2);
}

Member Data Documentation

METME BPHMonitor::BMass_

private

Definition at line 219 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

MEbinning BPHMonitor::Bmass_binning_

private

Definition at line 180 of file BPHMonitor.h.

Referenced by bookHistograms().

edm::InputTag BPHMonitor::bsInputTag_

private

Definition at line 162 of file BPHMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::BeamSpot> BPHMonitor::bsToken_

private

Definition at line 168 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::cos_binning_

private

Definition at line 183 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::d0_binning_

private

Definition at line 176 of file BPHMonitor.h.

MEbinning BPHMonitor::dca_binning_

private

Definition at line 181 of file BPHMonitor.h.

Referenced by bookHistograms().

GenericTriggerEventFlag* BPHMonitor::den_genTriggerEventFlag_

private

Definition at line 223 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and ~BPHMonitor().

METME BPHMonitor::DiMuDCA_

private

Definition at line 217 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMudR_

private

Definition at line 220 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuDS_

private

Definition at line 216 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuEta_

private

Definition at line 212 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuMass_

private

Definition at line 218 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuPhi_

private

Definition at line 211 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuProb_

private

Definition at line 215 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuPt_

private

Definition at line 213 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::DiMuPVcos_

private

Definition at line 214 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

int BPHMonitor::displaced_

private

Definition at line 233 of file BPHMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Candidate::LorentzVector, true> BPHMonitor::DMSelection_ref

private

Definition at line 264 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::dMu_pt_binning_

private

Definition at line 174 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::dR_binning_

private

Definition at line 178 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::ds_binning_

private

Definition at line 182 of file BPHMonitor.h.

Referenced by bookHistograms().

int BPHMonitor::enum_

private

Definition at line 237 of file BPHMonitor.h.

Referenced by analyze(), and bookHistograms().

MEbinning BPHMonitor::eta_binning_

private

Definition at line 175 of file BPHMonitor.h.

Referenced by bookHistograms().

std::string BPHMonitor::folderName_

private

Definition at line 158 of file BPHMonitor.h.

Referenced by bookHistograms().

edm::Handle<trigger::TriggerEvent> BPHMonitor::handleTriggerEvent

private

Definition at line 266 of file BPHMonitor.h.

Referenced by analyze(), and matchToTrigger().

std::string BPHMonitor::histoSuffix_

private

Definition at line 159 of file BPHMonitor.h.

HLTConfigProvider BPHMonitor::hltConfig_

private

Definition at line 268 of file BPHMonitor.h.

Referenced by bookHistograms(), getTriggerName(), and matchToTrigger().

edm::EDGetTokenT<trigger::TriggerEvent> BPHMonitor::hltInputTag_

private

Definition at line 259 of file BPHMonitor.h.

Referenced by analyze().

edm::InputTag BPHMonitor::hltInputTag_1

private

Definition at line 258 of file BPHMonitor.h.

Referenced by matchToTrigger().

std::vector<std::string> BPHMonitor::hltpaths_den

private

Definition at line 261 of file BPHMonitor.h.

Referenced by analyze().

std::vector<std::string> BPHMonitor::hltpaths_num

private

Definition at line 260 of file BPHMonitor.h.

Referenced by analyze().

HLTPrescaleProvider* BPHMonitor::hltPrescale_

private

Definition at line 224 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and ~BPHMonitor().

edm::Handle<edm::TriggerResults> BPHMonitor::HLTR

private

Definition at line 269 of file BPHMonitor.h.

edm::EDGetTokenT<edm::TriggerResults> BPHMonitor::hltTrigResTag_

private

Definition at line 257 of file BPHMonitor.h.

int BPHMonitor::Jpsi_

private

Definition at line 235 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::kaon_mass

private

Definition at line 249 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::L3_

private

Definition at line 231 of file BPHMonitor.h.

MEbinning BPHMonitor::mass_binning_

private

Definition at line 179 of file BPHMonitor.h.

Referenced by bookHistograms().

double BPHMonitor::max_dR

private

Definition at line 252 of file BPHMonitor.h.

double BPHMonitor::maxmass_

private

Definition at line 239 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassJpsi

private

Definition at line 241 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassJpsiTk

private

Definition at line 247 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassTkTk

private

Definition at line 245 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassUpsilon

private

Definition at line 243 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::min_dR

private

Definition at line 251 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::mincos

private

Definition at line 255 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minDS

private

Definition at line 256 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmass_

private

Definition at line 240 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassJpsi

private

Definition at line 242 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassJpsiTk

private

Definition at line 248 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassTkTk

private

Definition at line 246 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassUpsilon

private

Definition at line 244 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minprob

private

Definition at line 254 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::mu1d0_

private

Definition at line 195 of file BPHMonitor.h.

METME BPHMonitor::mu1Eta_

private

Definition at line 193 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu1Phi_

private

Definition at line 192 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu1Pt_

private

Definition at line 194 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu1z0_

private

Definition at line 196 of file BPHMonitor.h.

METME BPHMonitor::mu2d0_

private

Definition at line 200 of file BPHMonitor.h.

METME BPHMonitor::mu2Eta_

private

Definition at line 198 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu2Phi_

private

Definition at line 197 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu2Pt_

private

Definition at line 199 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu2z0_

private

Definition at line 201 of file BPHMonitor.h.

METME BPHMonitor::mu3d0_

private

Definition at line 205 of file BPHMonitor.h.

METME BPHMonitor::mu3Eta_

private

Definition at line 203 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu3Phi_

private

Definition at line 202 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu3Pt_

private

Definition at line 204 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::mu3z0_

private

Definition at line 206 of file BPHMonitor.h.

double BPHMonitor::mu_mass

private

Definition at line 250 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::mud0_

private

Definition at line 189 of file BPHMonitor.h.

Referenced by BPHMonitor().

METME BPHMonitor::muEta_

private

Definition at line 187 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

edm::InputTag BPHMonitor::muoInputTag_

private

Definition at line 161 of file BPHMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Muon, true> BPHMonitor::muoSelection_

private

Definition at line 225 of file BPHMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Muon, true> BPHMonitor::muoSelection_probe

private

Definition at line 228 of file BPHMonitor.h.

StringCutObjectSelector<reco::Muon, true> BPHMonitor::muoSelection_ref

private

Definition at line 226 of file BPHMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Muon, true> BPHMonitor::muoSelection_tag

private

Definition at line 227 of file BPHMonitor.h.

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

private

Definition at line 167 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::muPhi_

private

Definition at line 186 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::muPt_

private

Definition at line 188 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and BPHMonitor().

METME BPHMonitor::muz0_

private

Definition at line 190 of file BPHMonitor.h.

Referenced by BPHMonitor().

int BPHMonitor::nmuons_

private

Definition at line 229 of file BPHMonitor.h.

GenericTriggerEventFlag* BPHMonitor::num_genTriggerEventFlag_

private

Definition at line 222 of file BPHMonitor.h.

Referenced by analyze(), bookHistograms(), and ~BPHMonitor().

METME BPHMonitor::phEta_

private

Definition at line 209 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

MEbinning BPHMonitor::phi_binning_

private

Definition at line 172 of file BPHMonitor.h.

Referenced by bookHistograms().

edm::InputTag BPHMonitor::phInputTag_

private

Definition at line 164 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::phPhi_

private

Definition at line 208 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

METME BPHMonitor::phPt_

private

Definition at line 210 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

edm::EDGetTokenT<reco::PhotonCollection> BPHMonitor::phToken_

private

Definition at line 170 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::prob_binning_

private

Definition at line 184 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::pt_binning_

private

Definition at line 173 of file BPHMonitor.h.

Referenced by bookHistograms().

int BPHMonitor::ptCut_

private

Definition at line 232 of file BPHMonitor.h.

Referenced by analyze(), and badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons().

int BPHMonitor::seagull_

private

Definition at line 238 of file BPHMonitor.h.

Referenced by analyze().

bool BPHMonitor::tnp_

private

Definition at line 230 of file BPHMonitor.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag BPHMonitor::trInputTag_

private

Definition at line 163 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::trOrMu_

private

Definition at line 234 of file BPHMonitor.h.

Referenced by bookHistograms().

StringCutObjectSelector<reco::Track, true> BPHMonitor::trSelection_

private

Definition at line 262 of file BPHMonitor.h.

StringCutObjectSelector<reco::Track, true> BPHMonitor::trSelection_ref

private

Definition at line 263 of file BPHMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::TrackCollection> BPHMonitor::trToken_

private

Definition at line 169 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::Upsilon_

private

Definition at line 236 of file BPHMonitor.h.

Referenced by analyze().

edm::InputTag BPHMonitor::vtxInputTag_

private

Definition at line 165 of file BPHMonitor.h.

edm::EDGetTokenT<reco::VertexCollection> BPHMonitor::vtxToken_

private

Definition at line 171 of file BPHMonitor.h.

std::vector<bool> BPHMonitor::warningPrinted4token_

private

Definition at line 272 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

MEbinning BPHMonitor::z0_binning_

private

Definition at line 177 of file BPHMonitor.h.