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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)
static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)

Protected Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
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)

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_

Detailed Description

Definition at line 76 of file BPHMonitor.h.

Constructor & Destructor Documentation

BPHMonitor::BPHMonitor

(

const edm::ParameterSet &

iConfig

)

Definition at line 12 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 149 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  )

overrideprotected

Definition at line 373 of file BPHMonitor.cc.

References GenericTriggerEventFlag::accept(), ecalDrivenElectronSeedsParameters_cff::beamSpot, BMass_, bsInputTag_, bsToken_, ClosestApproachInRPhi::calculate(), reco::Vertex::chi2(), hiPixelPairStep_cff::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(), SiPixelPhase1TrackClustersV_cfi::e3, vertexPlots::e4, enum_, reco::tau::disc::Eta(), 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(), funct::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_, colinearityKinematic::Phi, phInputTag_, phPhi_, phPt_, phToken_, TransientVertex::position(), TransientVertex::positionError(), Prescale(), EnergyCorrector::pt, reco::tau::disc::Pt(), ptCut_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), reco::TrackBase::quality(), GlobalErrorBase< T, ErrorWeightType >::rerr(), seagull_, trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, lumiQTWidget::t, TrajectoryStateClosestToPoint::theState(), tnp_, trInputTag_, trSelection_ref, trToken_, Upsilon_, 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

        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  )

overrideprotected

Definition at line 237 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_, 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 189 of file BPHMonitor.cc.

References 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 194 of file BPHMonitor.cc.

References DQMStore::IBooker::book1D(), METME::denominator, pileupCalc::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 202 of file BPHMonitor.cc.

References 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 207 of file BPHMonitor.cc.

References 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 212 of file BPHMonitor.cc.

References DQMStore::IBooker::book2D(), METME::denominator, 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 225 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 1015 of file BPHMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), funct::false, fillHistoPSetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  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 1010 of file BPHMonitor.cc.

References edm::ParameterSetDescription::add().

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

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

static

Definition at line 1002 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 156 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 1126 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 1140 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(), electrons_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 1182 of file BPHMonitor.cc.

References DEFINE_FWK_MODULE, plotBeamSpotDB::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 180 of file BPHMonitor.cc.

References METME::denominator, METME::numerator, and 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 174 of file BPHMonitor.h.

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

MEbinning BPHMonitor::Bmass_binning_

private

Definition at line 134 of file BPHMonitor.h.

Referenced by bookHistograms().

edm::InputTag BPHMonitor::bsInputTag_

private

Definition at line 116 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 122 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::cos_binning_

private

Definition at line 137 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::d0_binning_

private

Definition at line 130 of file BPHMonitor.h.

MEbinning BPHMonitor::dca_binning_

private

Definition at line 135 of file BPHMonitor.h.

Referenced by bookHistograms().

GenericTriggerEventFlag* BPHMonitor::den_genTriggerEventFlag_

private

Definition at line 179 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuDCA_

private

Definition at line 172 of file BPHMonitor.h.

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

METME BPHMonitor::DiMudR_

private

Definition at line 175 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuDS_

private

Definition at line 171 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuEta_

private

Definition at line 167 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuMass_

private

Definition at line 173 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuPhi_

private

Definition at line 166 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuProb_

private

Definition at line 170 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuPt_

private

Definition at line 168 of file BPHMonitor.h.

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

METME BPHMonitor::DiMuPVcos_

private

Definition at line 169 of file BPHMonitor.h.

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

int BPHMonitor::displaced_

private

Definition at line 189 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 220 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::dMu_pt_binning_

private

Definition at line 128 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::dR_binning_

private

Definition at line 132 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::ds_binning_

private

Definition at line 136 of file BPHMonitor.h.

Referenced by bookHistograms().

int BPHMonitor::enum_

private

Definition at line 193 of file BPHMonitor.h.

Referenced by analyze(), and bookHistograms().

MEbinning BPHMonitor::eta_binning_

private

Definition at line 129 of file BPHMonitor.h.

Referenced by bookHistograms().

std::string BPHMonitor::folderName_

private

Definition at line 112 of file BPHMonitor.h.

Referenced by bookHistograms().

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

private

Definition at line 222 of file BPHMonitor.h.

Referenced by analyze(), and matchToTrigger().

std::string BPHMonitor::histoSuffix_

private

Definition at line 113 of file BPHMonitor.h.

HLTConfigProvider BPHMonitor::hltConfig_

private

Definition at line 224 of file BPHMonitor.h.

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

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

private

Definition at line 215 of file BPHMonitor.h.

Referenced by analyze().

edm::InputTag BPHMonitor::hltInputTag_1

private

Definition at line 214 of file BPHMonitor.h.

Referenced by matchToTrigger().

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

private

Definition at line 217 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 216 of file BPHMonitor.h.

Referenced by analyze().

HLTPrescaleProvider* BPHMonitor::hltPrescale_

private

Definition at line 180 of file BPHMonitor.h.

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

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

private

Definition at line 225 of file BPHMonitor.h.

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

private

Definition at line 213 of file BPHMonitor.h.

int BPHMonitor::Jpsi_

private

Definition at line 191 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::kaon_mass

private

Definition at line 205 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::L3_

private

Definition at line 187 of file BPHMonitor.h.

MEbinning BPHMonitor::mass_binning_

private

Definition at line 133 of file BPHMonitor.h.

Referenced by bookHistograms().

double BPHMonitor::max_dR

private

Definition at line 208 of file BPHMonitor.h.

double BPHMonitor::maxmass_

private

Definition at line 195 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassJpsi

private

Definition at line 197 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassJpsiTk

private

Definition at line 203 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassTkTk

private

Definition at line 201 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::maxmassUpsilon

private

Definition at line 199 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::min_dR

private

Definition at line 207 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::mincos

private

Definition at line 211 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minDS

private

Definition at line 212 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmass_

private

Definition at line 196 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassJpsi

private

Definition at line 198 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassJpsiTk

private

Definition at line 204 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassTkTk

private

Definition at line 202 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minmassUpsilon

private

Definition at line 200 of file BPHMonitor.h.

Referenced by analyze().

double BPHMonitor::minprob

private

Definition at line 210 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::mu1d0_

private

Definition at line 149 of file BPHMonitor.h.

METME BPHMonitor::mu1Eta_

private

Definition at line 147 of file BPHMonitor.h.

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

METME BPHMonitor::mu1Phi_

private

Definition at line 146 of file BPHMonitor.h.

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

METME BPHMonitor::mu1Pt_

private

Definition at line 148 of file BPHMonitor.h.

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

METME BPHMonitor::mu1z0_

private

Definition at line 150 of file BPHMonitor.h.

METME BPHMonitor::mu2d0_

private

Definition at line 154 of file BPHMonitor.h.

METME BPHMonitor::mu2Eta_

private

Definition at line 152 of file BPHMonitor.h.

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

METME BPHMonitor::mu2Phi_

private

Definition at line 151 of file BPHMonitor.h.

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

METME BPHMonitor::mu2Pt_

private

Definition at line 153 of file BPHMonitor.h.

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

METME BPHMonitor::mu2z0_

private

Definition at line 155 of file BPHMonitor.h.

METME BPHMonitor::mu3d0_

private

Definition at line 159 of file BPHMonitor.h.

METME BPHMonitor::mu3Eta_

private

Definition at line 157 of file BPHMonitor.h.

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

METME BPHMonitor::mu3Phi_

private

Definition at line 156 of file BPHMonitor.h.

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

METME BPHMonitor::mu3Pt_

private

Definition at line 158 of file BPHMonitor.h.

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

METME BPHMonitor::mu3z0_

private

Definition at line 160 of file BPHMonitor.h.

double BPHMonitor::mu_mass

private

Definition at line 206 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::mud0_

private

Definition at line 143 of file BPHMonitor.h.

Referenced by BPHMonitor().

METME BPHMonitor::muEta_

private

Definition at line 141 of file BPHMonitor.h.

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

edm::InputTag BPHMonitor::muoInputTag_

private

Definition at line 115 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 181 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 184 of file BPHMonitor.h.

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

private

Definition at line 182 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 183 of file BPHMonitor.h.

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

private

Definition at line 121 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::muPhi_

private

Definition at line 140 of file BPHMonitor.h.

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

METME BPHMonitor::muPt_

private

Definition at line 142 of file BPHMonitor.h.

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

METME BPHMonitor::muz0_

private

Definition at line 144 of file BPHMonitor.h.

Referenced by BPHMonitor().

int BPHMonitor::nmuons_

private

Definition at line 185 of file BPHMonitor.h.

GenericTriggerEventFlag* BPHMonitor::num_genTriggerEventFlag_

private

Definition at line 178 of file BPHMonitor.h.

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

METME BPHMonitor::phEta_

private

Definition at line 164 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

MEbinning BPHMonitor::phi_binning_

private

Definition at line 126 of file BPHMonitor.h.

Referenced by bookHistograms().

edm::InputTag BPHMonitor::phInputTag_

private

Definition at line 118 of file BPHMonitor.h.

Referenced by analyze().

METME BPHMonitor::phPhi_

private

Definition at line 163 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

METME BPHMonitor::phPt_

private

Definition at line 165 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

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

private

Definition at line 124 of file BPHMonitor.h.

Referenced by analyze().

MEbinning BPHMonitor::prob_binning_

private

Definition at line 138 of file BPHMonitor.h.

Referenced by bookHistograms().

MEbinning BPHMonitor::pt_binning_

private

Definition at line 127 of file BPHMonitor.h.

Referenced by bookHistograms().

int BPHMonitor::ptCut_

private

Definition at line 188 of file BPHMonitor.h.

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

int BPHMonitor::seagull_

private

Definition at line 194 of file BPHMonitor.h.

Referenced by analyze().

bool BPHMonitor::tnp_

private

Definition at line 186 of file BPHMonitor.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag BPHMonitor::trInputTag_

private

Definition at line 117 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::trOrMu_

private

Definition at line 190 of file BPHMonitor.h.

Referenced by bookHistograms().

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

private

Definition at line 218 of file BPHMonitor.h.

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

private

Definition at line 219 of file BPHMonitor.h.

Referenced by analyze().

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

private

Definition at line 123 of file BPHMonitor.h.

Referenced by analyze().

int BPHMonitor::Upsilon_

private

Definition at line 192 of file BPHMonitor.h.

Referenced by analyze().

edm::InputTag BPHMonitor::vtxInputTag_

private

Definition at line 119 of file BPHMonitor.h.

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

private

Definition at line 125 of file BPHMonitor.h.

std::vector<bool> BPHMonitor::warningPrinted4token_

private

Definition at line 228 of file BPHMonitor.h.

Referenced by analyze(), and BPHMonitor().

MEbinning BPHMonitor::z0_binning_

private

Definition at line 131 of file BPHMonitor.h.