QcdUeDQM Class Reference

#include <QcdUeDQM.h>

Inheritance diagram for QcdUeDQM:

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	QcdUeDQM (const edm::ParameterSet &parameters)
	~QcdUeDQM () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
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)

Private Types
typedef std::vector< const reco::Track * >	container

Private Member Functions
void	book1D (DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, bool sumw2=true, bool sbox=true)
void	bookProfile (DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, double y1, double y2, bool sumw2=true, bool sbox=true)
void	fill1D (std::vector< TH1F * > &hs, double val, double w=1.)
void	fill1D (std::vector< MonitorElement * > &mes, double val, double w=1.)
void	fill2D (std::vector< TH2F * > &hs, double valx, double valy, double w=1.)
void	fill2D (std::vector< MonitorElement * > &mes, double valx, double valy, double w=1.)
void	fill3D (std::vector< TH3F * > &hs, int gbin, double w=1.)
void	fillChargedJetSpectra (const edm::Handle< reco::TrackJetCollection > trackJets)
void	fillHltBits (const edm::Event &iEvent, const edm::EventSetup &iSetup)
void	fillProfile (std::vector< TProfile * > &hs, double valx, double valy, double w=1.)
void	fillProfile (std::vector< MonitorElement * > &mes, double valx, double valy, double w=1.)
void	fillpTMaxRelated (const std::vector< const reco::Track * > &track)
void	fillUE_with_ChargedJets (const std::vector< const reco::Track * > &track, const edm::Handle< reco::TrackJetCollection > &trackJets)
void	fillUE_with_MaxpTtrack (const std::vector< const reco::Track * > &track)
bool	fillVtxPlots (const reco::BeamSpot *bs, const edm::Handle< reco::VertexCollection > vtxColl)
template<typename TYPE >
void	getProduct (const std::string name, edm::Handle< TYPE > &prod, const edm::Event &event) const
template<typename TYPE >
bool	getProductSafe (const std::string name, edm::Handle< TYPE > &prod, const edm::Event &event) const
void	setLabel1D (std::vector< MonitorElement * > &mes)
bool	trackSelection (const reco::Track &trk, const reco::BeamSpot *bs, const reco::Vertex &vtx, int sizevtx)

Private Attributes
std::vector< reco::TrackBase::TrackAlgorithm >	algorithm_
bool	allowTriplets_
edm::EDGetTokenT< reco::BeamSpot >	bsLabel_
double	bsPos_
bool	bsuse_
edm::EDGetTokenT< reco::CaloJetCollection >	caloJetLabel_
edm::EDGetTokenT< reco::TrackJetCollection >	chargedJetLabel_
double	diffvtxbs_
MonitorElement *	h2TrigCorr_
std::vector< MonitorElement * >	hBeamSpot_x_
std::vector< MonitorElement * >	hBeamSpot_y_
std::vector< MonitorElement * >	hBeamSpot_z_
std::vector< MonitorElement * >	hChargedJetConstituent_
std::vector< MonitorElement * >	hChargedJetMulti_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Away500_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Away900_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Toward500_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Toward900_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Transverse500_
std::vector< MonitorElement * >	hdNdEtadPhi_pTMax_Transverse900_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Away500_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Away900_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Toward500_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Toward900_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Transverse500_
std::vector< MonitorElement * >	hdNdEtadPhi_trackJet_Transverse900_
std::vector< MonitorElement * >	hdPhi_chargedJet_tracks_
std::vector< MonitorElement * >	hdPhi_maxpTTrack_tracks_
std::vector< MonitorElement * >	hGoodTrkEta500_
std::vector< MonitorElement * >	hGoodTrkEta900_
std::vector< MonitorElement * >	hGoodTrkPhi500_
std::vector< MonitorElement * >	hGoodTrkPhi900_
std::vector< MonitorElement * >	hGoodTrkPt500_
std::vector< MonitorElement * >	hGoodTrkPt900_
std::vector< MonitorElement * >	hLeadingChargedJet_etaSpectrum_
std::vector< MonitorElement * >	hLeadingChargedJet_phiSpectrum_
std::vector< MonitorElement * >	hLeadingChargedJet_pTSpectrum_
std::vector< MonitorElement * >	hLeadingTrack_etaSpectrum_
std::vector< MonitorElement * >	hLeadingTrack_phiSpectrum_
std::vector< MonitorElement * >	hLeadingTrack_pTSpectrum_
HLTConfigProvider	hltConfig
std::vector< std::string >	hltProcNames_
std::string	hltResName_
std::vector< int >	hltTrgBits_
std::vector< bool >	hltTrgDeci_
std::vector< std::string >	hltTrgNames_
std::vector< std::string >	hltTrgUsedNames_
std::string	hltUsedResName_
std::vector< MonitorElement * >	hNevts_
std::vector< MonitorElement * >	hNgoodTrk_
std::vector< MonitorElement * >	hNtrackerLayer_
std::vector< MonitorElement * >	hNtrackerPixelLayer_
std::vector< MonitorElement * >	hNtrackerStripPixelLayer_
std::vector< MonitorElement * >	hNvertices_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Away500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Away900_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Toward500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Toward900_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Transverse500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_pTMax_Transverse900_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Away500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Away900_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Toward500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Toward900_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Transverse500_
std::vector< MonitorElement * >	hpTSumdEtadPhi_trackJet_Transverse900_
std::vector< MonitorElement * >	hRatioDxySigmaDxyBS_
std::vector< MonitorElement * >	hRatioDxySigmaDxyPV_
std::vector< MonitorElement * >	hRatioDzSigmaDzBS_
std::vector< MonitorElement * >	hRatioDzSigmaDzPV_
std::vector< MonitorElement * >	hRatioPtErrorPt_
std::vector< MonitorElement * >	hTrkChi2_
std::vector< MonitorElement * >	hTrkEta_
std::vector< MonitorElement * >	hTrkNdof_
std::vector< MonitorElement * >	hTrkPhi_
std::vector< MonitorElement * >	hTrkPt_
std::vector< MonitorElement * >	hVertex_ndof_
std::vector< MonitorElement * >	hVertex_rho_
std::vector< MonitorElement * >	hVertex_x_
std::vector< MonitorElement * >	hVertex_y_
std::vector< MonitorElement * >	hVertex_z_
std::vector< MonitorElement * >	hVertex_z_bs_
bool	isHltConfigSuccessful_
double	lip_
double	maxChi2_
double	maxRapidity_
int	min3DHit_
int	minHit_
double	minRapidity_
double	ptErr_pt_
double	ptMin_
double	pxlLayerMinCut_
std::vector< reco::TrackBase::TrackQuality >	quality_
MonitorElement *	repSumMap_
MonitorElement *	repSummary_
bool	requirePIX1_
container	selected_
const TrackerGeometry *	tgeo_
double	tip_
edm::EDGetTokenT< reco::TrackCollection >	trackLabel_
int	verbose_
reco::Vertex	vtx1
edm::EDGetTokenT< reco::VertexCollection >	vtxLabel_
double	vtxntk_

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
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from 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)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 53 of file QcdUeDQM.h.

Member Typedef Documentation

typedef std::vector<const reco::Track *> QcdUeDQM::container

private

Definition at line 398 of file QcdUeDQM.h.

Constructor & Destructor Documentation

QcdUeDQM::QcdUeDQM

(

const edm::ParameterSet &

parameters

)

Definition at line 52 of file QcdUeDQM.cc.

References reco::TrackBase::algoByName(), qcdUeDQM_cfi::algorithm, algorithm_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltProcNames_, hltTrgNames_, isHltConfigSuccessful_, dqmiolumiharvest::j, qcdUeDQM_cfi::quality, quality_, and reco::TrackBase::qualityByName().

    : hltResName_(parameters.getUntrackedParameter<string>("hltTrgResults")),
      verbose_(parameters.getUntrackedParameter<int>("verbose", 3)),
      tgeo_(nullptr),
      repSumMap_(nullptr),
      repSummary_(nullptr),
      h2TrigCorr_(nullptr),
      ptMin_(parameters.getParameter<double>("ptMin")),
      minRapidity_(parameters.getParameter<double>("minRapidity")),
      maxRapidity_(parameters.getParameter<double>("maxRapidity")),
      tip_(parameters.getParameter<double>("tip")),
      lip_(parameters.getParameter<double>("lip")),
      diffvtxbs_(parameters.getParameter<double>("diffvtxbs")),
      ptErr_pt_(parameters.getParameter<double>("ptErr_pt")),
      vtxntk_(parameters.getParameter<double>("vtxntk")),
      minHit_(parameters.getParameter<int>("minHit")),
      pxlLayerMinCut_(parameters.getParameter<double>("pxlLayerMinCut")),
      requirePIX1_(parameters.getParameter<bool>("requirePIX1")),
      min3DHit_(parameters.getParameter<int>("min3DHit")),
      maxChi2_(parameters.getParameter<double>("maxChi2")),
      bsuse_(parameters.getParameter<bool>("bsuse")),
      allowTriplets_(parameters.getParameter<bool>("allowTriplets")),
      bsPos_(parameters.getParameter<double>("bsPos")),
      caloJetLabel_(consumes<reco::CaloJetCollection>(parameters.getUntrackedParameter<edm::InputTag>("caloJetTag"))),
      chargedJetLabel_(
          consumes<reco::TrackJetCollection>(parameters.getUntrackedParameter<edm::InputTag>("chargedJetTag"))),
      trackLabel_(consumes<reco::TrackCollection>(parameters.getUntrackedParameter<edm::InputTag>("trackTag"))),
      vtxLabel_(consumes<reco::VertexCollection>(parameters.getUntrackedParameter<edm::InputTag>("vtxTag"))),
      bsLabel_(consumes<reco::BeamSpot>(parameters.getParameter<edm::InputTag>("beamSpotTag"))) {
  // Constructor.
  std::vector<std::string> quality = parameters.getParameter<std::vector<std::string> >("quality");
  for (unsigned int j = 0; j < quality.size(); j++)
    quality_.push_back(reco::TrackBase::qualityByName(quality[j]));
  std::vector<std::string> algorithm = parameters.getParameter<std::vector<std::string> >("algorithm");
  for (unsigned int j = 0; j < algorithm.size(); j++)
    algorithm_.push_back(reco::TrackBase::algoByName(algorithm[j]));

  if (parameters.exists("hltTrgNames"))
    hltTrgNames_ = parameters.getUntrackedParameter<vector<string> >("hltTrgNames");

  if (parameters.exists("hltProcNames"))
    hltProcNames_ = parameters.getUntrackedParameter<vector<string> >("hltProcNames");
  else {
    //     hltProcNames_.push_back("FU");
    hltProcNames_.push_back("HLT");
  }

  isHltConfigSuccessful_ = false;  // init
}

QcdUeDQM::~QcdUeDQM ( )

override

Definition at line 102 of file QcdUeDQM.cc.

                    {
  // Destructor.
}

Member Function Documentation

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

override

Definition at line 636 of file QcdUeDQM.cc.

References pwdgSkimBPark_cfi::beamSpot, bsLabel_, caloJetLabel_, chargedJetLabel_, fill1D(), fillChargedJetSpectra(), fillHltBits(), fillpTMaxRelated(), fillUE_with_ChargedJets(), fillUE_with_MaxpTtrack(), fillVtxPlots(), edm::Event::getByToken(), hNevts_, isHltConfigSuccessful_, edm::Handle< T >::product(), selected_, trackLabel_, PDWG_EXOHSCP_cff::tracks, trackSelection(), vtx1, and vtxLabel_.

                                                                    {
  if (!isHltConfigSuccessful_)
    return;

  // Analyze the given event.

  edm::Handle<reco::BeamSpot> beamSpot;
  bool ValidBS_ = iEvent.getByToken(bsLabel_, beamSpot);
  if (!ValidBS_)
    return;

  edm::Handle<reco::TrackCollection> tracks;
  bool ValidTrack_ = iEvent.getByToken(trackLabel_, tracks);
  if (!ValidTrack_)
    return;

  edm::Handle<reco::TrackJetCollection> trkJets;
  bool ValidTrackJet_ = iEvent.getByToken(chargedJetLabel_, trkJets);
  if (!ValidTrackJet_)
    return;

  edm::Handle<reco::CaloJetCollection> calJets;
  bool ValidCaloJet_ = iEvent.getByToken(caloJetLabel_, calJets);
  if (!ValidCaloJet_)
    return;

  edm::Handle<reco::VertexCollection> vertexColl;
  bool ValidVtxColl_ = iEvent.getByToken(vtxLabel_, vertexColl);
  if (!ValidVtxColl_)
    return;

  reco::TrackCollection tracks_sort = *tracks;
  std::sort(tracks_sort.begin(), tracks_sort.end(), PtSorter());

  // get tracker geometry
  /*  ESHandle<TrackerGeometry> trackerHandle;
   iSetup.get<TrackerDigiGeometryRecord>().get(trackerHandle);
   tgeo_ = trackerHandle.product();
   if (!tgeo_)return;
 */
  selected_.clear();
  fillHltBits(iEvent, iSetup);
  // select good tracks
  if (fillVtxPlots(beamSpot.product(), vertexColl)) {
    fill1D(hNevts_, 1);
    for (reco::TrackCollection::const_iterator Trk = tracks_sort.begin(); Trk != tracks_sort.end(); ++Trk) {
      if (trackSelection(*Trk, beamSpot.product(), vtx1, vertexColl->size()))
        selected_.push_back(&*Trk);
    }

    fillpTMaxRelated(selected_);
    fillChargedJetSpectra(trkJets);
    //    fillCaloJetSpectra(calJets);
    fillUE_with_MaxpTtrack(selected_);
    if (!trkJets->empty())
      fillUE_with_ChargedJets(selected_, trkJets);
    // if(calJets->size()>0)fillUE_with_CaloJets(selected_,calJets);
  }
}

void QcdUeDQM::book1D ( DQMStore::IBooker &  , std::vector< MonitorElement * > &  mes, const std::string &  name, const std::string &  title, int  nx, double  x1, double  x2, bool  sumw2 = true, bool  sbox = true  )

private

Definition at line 696 of file QcdUeDQM.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), MillePedeFileConverter_cfg::e, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, dqm::impl::MonitorElement::getTH1(), hltTrgUsedNames_, mps_fire::i, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms().

                                 {
  // Book 1D histos.
  for (size_t i = 0; i < hltTrgUsedNames_.size(); ++i) {
    std::string folderName = "Physics/QcdUe/" + hltTrgUsedNames_.at(i);
    iBooker.setCurrentFolder(folderName);
    MonitorElement *e = iBooker.book1D(Form("%s_%s", name.c_str(), hltTrgUsedNames_.at(i).c_str()),
                                       Form("%s: %s", hltTrgUsedNames_.at(i).c_str(), title.c_str()),
                                       nx,
                                       x1,
                                       x2);
    TH1 *h1 = e->getTH1();
    if (sumw2) {
      if (0 == h1->GetSumw2N()) {  // protect against re-summing (would cause
                                   // exception)
        h1->Sumw2();
      }
    }
    h1->SetStats(sbox);
    mes.push_back(e);
  }
}

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 164 of file QcdUeDQM.cc.

References book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), bookProfile(), dqm::impl::MonitorElement::getTH1(), h, h2TrigCorr_, hBeamSpot_x_, hBeamSpot_y_, hBeamSpot_z_, hChargedJetMulti_, hdNdEtadPhi_pTMax_Away500_, hdNdEtadPhi_pTMax_Away900_, hdNdEtadPhi_pTMax_Toward500_, hdNdEtadPhi_pTMax_Toward900_, hdNdEtadPhi_pTMax_Transverse500_, hdNdEtadPhi_pTMax_Transverse900_, hdNdEtadPhi_trackJet_Away500_, hdNdEtadPhi_trackJet_Away900_, hdNdEtadPhi_trackJet_Toward500_, hdNdEtadPhi_trackJet_Toward900_, hdNdEtadPhi_trackJet_Transverse500_, hdNdEtadPhi_trackJet_Transverse900_, hdPhi_chargedJet_tracks_, hdPhi_maxpTTrack_tracks_, hGoodTrkEta500_, hGoodTrkEta900_, hGoodTrkPhi500_, hGoodTrkPhi900_, hGoodTrkPt500_, hGoodTrkPt900_, hLeadingChargedJet_etaSpectrum_, hLeadingChargedJet_phiSpectrum_, hLeadingChargedJet_pTSpectrum_, hLeadingTrack_etaSpectrum_, hLeadingTrack_phiSpectrum_, hLeadingTrack_pTSpectrum_, hltTrgUsedNames_, hNevts_, hNgoodTrk_, hNtrackerLayer_, hNtrackerPixelLayer_, hNtrackerStripPixelLayer_, hNvertices_, hpTSumdEtadPhi_pTMax_Away500_, hpTSumdEtadPhi_pTMax_Away900_, hpTSumdEtadPhi_pTMax_Toward500_, hpTSumdEtadPhi_pTMax_Toward900_, hpTSumdEtadPhi_pTMax_Transverse500_, hpTSumdEtadPhi_pTMax_Transverse900_, hpTSumdEtadPhi_trackJet_Away500_, hpTSumdEtadPhi_trackJet_Away900_, hpTSumdEtadPhi_trackJet_Toward500_, hpTSumdEtadPhi_trackJet_Toward900_, hpTSumdEtadPhi_trackJet_Transverse500_, hpTSumdEtadPhi_trackJet_Transverse900_, hRatioDxySigmaDxyBS_, hRatioDxySigmaDxyPV_, hRatioDzSigmaDzBS_, hRatioDzSigmaDzPV_, hRatioPtErrorPt_, hTrkChi2_, hTrkEta_, hTrkNdof_, hTrkPhi_, hTrkPt_, hVertex_ndof_, hVertex_rho_, hVertex_x_, hVertex_y_, hVertex_z_, hVertex_z_bs_, mps_fire::i, dqm::impl::MonitorElement::setBinLabel(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

                                                                                             {
  // Book histograms if needed.
  iBooker.setCurrentFolder("Physics/QcdUe");

  if (true) {
    const int Nx = hltTrgUsedNames_.size();
    const double x1 = -0.5;
    const double x2 = Nx - 0.5;
    h2TrigCorr_ = iBooker.book2D("h2TriCorr", "Trigger bit x vs y;y&&!x;x&&y", Nx, x1, x2, Nx, x1, x2);
    for (size_t i = 1; i <= hltTrgUsedNames_.size(); ++i) {
      h2TrigCorr_->setBinLabel(i, hltTrgUsedNames_.at(i - 1), 1);
      h2TrigCorr_->setBinLabel(i, hltTrgUsedNames_.at(i - 1), 2);
    }
    TH1 *h = h2TrigCorr_->getTH1();
    if (h)
      h->SetStats(false);
  }
  book1D(iBooker, hNevts_, "hNevts", "number of events", 2, 0, 2);
  book1D(iBooker, hNtrackerLayer_, "hNtrackerLayer", "number of tracker layers;multiplicity", 20, -0.5, 19.5);
  book1D(iBooker, hNtrackerPixelLayer_, "hNtrackerPixelLayer", "number of pixel layers;multiplicity", 10, -0.5, 9.5);
  book1D(iBooker,
         hNtrackerStripPixelLayer_,
         "hNtrackerStripPixelLayer",
         "number of strip + pixel layers;multiplicity",
         30,
         -0.5,
         39.5);
  book1D(iBooker, hRatioPtErrorPt_, "hRatioPtErrorPt", "ratio of pT error and track pT", 25, 0., 5.);
  book1D(iBooker, hTrkPt_, "hTrkPt", "pT of all tracks", 50, 0., 50.);
  book1D(iBooker, hTrkEta_, "hTrkEta", "eta of all tracks", 40, -4., 4.);
  book1D(iBooker, hTrkPhi_, "hTrkPhi", "phi of all tracks", 40, -4., 4.);
  book1D(iBooker,
         hRatioDxySigmaDxyBS_,
         "hRatioDxySigmaDxyBS",
         "ratio of transverse impact parameter and its significance wrt beam spot",
         60,
         -10.,
         10);
  book1D(iBooker,
         hRatioDxySigmaDxyPV_,
         "hRatioDxySigmaDxyPV",
         "ratio of transverse impact parameter and its significance wrt PV",
         60,
         -10.,
         10);
  book1D(iBooker,
         hRatioDzSigmaDzBS_,
         "hRatioDzSigmaDzBS",
         "ratio of longitudinal impact parameter and its significance wrt beam "
         "spot",
         80,
         -20.,
         20);
  book1D(iBooker,
         hRatioDzSigmaDzPV_,
         "hRatioDzSigmaDzPV",
         "ratio of longitudinal impact parameter and its significance wrt PV",
         80,
         -20.,
         20);
  book1D(iBooker, hTrkChi2_, "hTrkChi2", "track chi2", 30, 0., 30);
  book1D(iBooker, hTrkNdof_, "hTrkNdof", "track NDOF", 100, 0, 100);

  book1D(iBooker, hNgoodTrk_, "hNgoodTrk", "number of good tracks", 50, -0.5, 49.5);

  book1D(iBooker, hGoodTrkPt500_, "hGoodTrkPt500", "pT of all good tracks with pT > 500 MeV", 50, 0., 50.);
  book1D(iBooker, hGoodTrkEta500_, "hGoodTrkEta500", "eta of all good tracks pT > 500 MeV", 40, -4., 4.);
  book1D(iBooker, hGoodTrkPhi500_, "hGoodTrkPhi500", "phi of all good tracks pT > 500 MeV", 40, -4., 4.);

  book1D(iBooker, hGoodTrkPt900_, "hGoodTrkPt900", "pT of all good tracks with pT > 900 MeV", 50, 0., 50.);
  book1D(iBooker, hGoodTrkEta900_, "hGoodTrkEta900", "eta of all good tracks pT > 900 MeV", 40, -4., 4.);
  book1D(iBooker, hGoodTrkPhi900_, "hGoodTrkPhi900", "phi of all good tracks pT > 900 MeV", 40, -4., 4.);

  book1D(iBooker, hNvertices_, "hNvertices", "number of vertices", 5, -0.5, 4.5);
  book1D(iBooker, hVertex_z_, "hVertex_z", "z position of vertex; z[cm]", 200, -50, 50);
  book1D(iBooker, hVertex_y_, "hVertex_y", "y position of vertex; y[cm]", 100, -5, 5);
  book1D(iBooker, hVertex_x_, "hVertex_x", "x position of vertex; x[cm]", 100, -5, 5);
  book1D(iBooker, hVertex_ndof_, "hVertex_ndof", "ndof of vertex", 100, 0, 100);
  book1D(iBooker, hVertex_rho_, "hVertex_rho", "rho of vertex", 100, 0, 5);
  book1D(iBooker, hVertex_z_bs_, "hVertex_z_bs", "z position of vertex from beamspot; z[cm]", 200, -50, 50);

  book1D(iBooker, hBeamSpot_z_, "hBeamSpot_z", "z position of beamspot; z[cm]", 100, -20, 20);
  book1D(iBooker, hBeamSpot_y_, "hBeamSpot_y", "y position of beamspot; y[cm]", 50, -10, 10);
  book1D(iBooker, hBeamSpot_x_, "hBeamSpot_x", "x position of beamspot; x[cm]", 50, -10, 10);

  if (true) {
    const int Nx = 25;
    const double x1 = 0.0;
    const double x2 = 50.0;
    book1D(iBooker,
           hLeadingTrack_pTSpectrum_,
           "hLeadingTrack_pTSpectrum",
           "pT spectrum of leading track;pT(GeV/c)",
           Nx,
           x1,
           x2);
    book1D(iBooker,
           hLeadingChargedJet_pTSpectrum_,
           "hLeadingChargedJet_pTSpectrum",
           "pT spectrum of leading track jet;pT(GeV/c)",
           Nx,
           x1,
           x2);
  }

  if (true) {
    const int Nx = 24;
    const double x1 = -4.;
    const double x2 = 4.;
    book1D(iBooker,
           hLeadingTrack_phiSpectrum_,
           "hLeadingTrack_phiSpectrum",
           "#phi spectrum of leading track;#phi",
           Nx,
           x1,
           x2);
    book1D(iBooker,
           hLeadingChargedJet_phiSpectrum_,
           "hLeadingChargedJet_phiSpectrum",
           "#phi spectrum of leading track jet;#phi",
           Nx,
           x1,
           x2);
  }

  if (true) {
    const int Nx = 24;
    const double x1 = -4.;
    const double x2 = 4.;
    book1D(iBooker,
           hLeadingTrack_etaSpectrum_,
           "hLeadingTrack_etaSpectrum",
           "#eta spectrum of leading track;#eta",
           Nx,
           x1,
           x2);
    book1D(iBooker,
           hLeadingChargedJet_etaSpectrum_,
           "hLeadingChargedJet_etaSpectrum",
           "#eta spectrum of leading track jet;#eta",
           Nx,
           x1,
           x2);
  }

  if (true) {
    const int Nx = 75;
    const double x1 = 0.0;
    const double x2 = 75.0;
    const double y1 = 0.;
    const double y2 = 10.;
    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Toward500_,
                "hdNdEtadPhi_pTMax_Toward500",
                "Average number of tracks (pT > 500 MeV) in toward region vs "
                "leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Transverse500_,
                "hdNdEtadPhi_pTMax_Transverse500",
                "Average number of tracks (pT > 500 MeV) in transverse region "
                "vs leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Away500_,
                "hdNdEtadPhi_pTMax_Away500",
                "Average number of tracks (pT > 500 MeV) in away region vs "
                "leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Toward500_,
                "hdNdEtadPhi_trackJet_Toward500",
                "Average number of tracks (pT > 500 MeV) in toward region vs "
                "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Transverse500_,
                "hdNdEtadPhi_trackJet_Transverse500",
                "Average number of tracks (pT > 500 MeV) in transverse region "
                "vs leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Away500_,
                "hdNdEtadPhi_trackJet_Away500",
                "Average number of tracks (pT > 500 MeV) in away region vs "
                "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);

    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Toward500_,
                "hpTSumdEtadPhi_pTMax_Toward500",
                "Average number of tracks (pT > 500 MeV) in toward region vs "
                "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Transverse500_,
                "hpTSumdEtadPhi_pTMax_Transverse500",
                "Average number of tracks (pT > 500 MeV) in transverse region "
                "vs leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Away500_,
                "hpTSumdEtadPhi_pTMax_Away500",
                "Average number of tracks (pT > 500 MeV) in away region vs "
                "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Toward500_,
                "hpTSumdEtadPhi_trackJet_Toward500",
                "Average number of tracks (pT > 500 MeV) in toward region vs "
                "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Transverse500_,
                "hpTSumdEtadPhi_trackJet_Transverse500",
                "Average number of tracks (pT > 500 MeV) in transverse region "
                "vs leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Away500_,
                "hpTSumdEtadPhi_trackJet_Away500",
                "Average number of tracks (pT > 500 MeV) in away region vs "
                "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);

    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Toward900_,
                "hdNdEtadPhi_pTMax_Toward900",
                "Average number of tracks (pT > 900 MeV) in toward region vs "
                "leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Transverse900_,
                "hdNdEtadPhi_pTMax_Transverse900",
                "Average number of tracks (pT > 900 MeV) in transverse region "
                "vs leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_pTMax_Away900_,
                "hdNdEtadPhi_pTMax_Away900",
                "Average number of tracks (pT > 900 MeV) in away region vs "
                "leading track pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Toward900_,
                "hdNdEtadPhi_trackJet_Toward900",
                "Average number of tracks (pT > 900 MeV) in toward region vs "
                "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Transverse900_,
                "hdNdEtadPhi_trackJet_Transverse900",
                "Average number of tracks (pT > 900 MeV) in transverse region "
                "vs leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hdNdEtadPhi_trackJet_Away900_,
                "hdNdEtadPhi_trackJet_Away900",
                "Average number of tracks (pT > 900 MeV) in away region vs "
                "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);

    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Toward900_,
                "hpTSumdEtadPhi_pTMax_Toward900",
                "Average number of tracks (pT > 900 MeV) in toward region vs "
                "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Transverse900_,
                "hpTSumdEtadPhi_pTMax_Transverse900",
                "Average number of tracks (pT > 900 MeV) in transverse region "
                "vs leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_pTMax_Away900_,
                "hpTSumdEtadPhi_pTMax_Away900",
                "Average number of tracks (pT > 900 MeV) in away region vs "
                "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Toward900_,
                "hpTSumdEtadPhi_trackJet_Toward900",
                "Average number of tracks (pT > 900 MeV) in toward region vs "
                "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Transverse900_,
                "hpTSumdEtadPhi_trackJet_Transverse900",
                "Average number of tracks (pT > 900 MeV) in transverse region "
                "vs leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
    bookProfile(iBooker,
                hpTSumdEtadPhi_trackJet_Away900_,
                "hpTSumdEtadPhi_trackJet_Away900",
                "Average number of tracks (pT > 900 MeV) in away region vs "
                "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
                Nx,
                x1,
                x2,
                y1,
                y2,
                false,
                false);
  }

  if (true) {
    const int Nx = 20;
    const double x1 = 0.0;
    const double x2 = 20.0;

    book1D(iBooker, hChargedJetMulti_, "hChargedJetMulti", "Charged jet multiplicity;multiplicities", Nx, x1, x2);
  }

  if (true) {
    const int Nx = 60;
    const double x1 = -180.0;
    const double x2 = 180.0;

    book1D(iBooker,
           hdPhi_maxpTTrack_tracks_,
           "hdPhi_maxpTTrack_tracks",
           "delta phi between leading tracks and other "
           "tracks;#Delta#phi(leading track-track)",
           Nx,
           x1,
           x2);
    book1D(iBooker,
           hdPhi_chargedJet_tracks_,
           "hdPhi_chargedJet_tracks",
           "delta phi between leading charged jet  and "
           "tracks;#Delta#phi(leading charged jet-track)",
           Nx,
           x1,
           x2);
  }
}

void QcdUeDQM::bookProfile ( DQMStore::IBooker &  , std::vector< MonitorElement * > &  mes, const std::string &  name, const std::string &  title, int  nx, double  x1, double  x2, double  y1, double  y2, bool  sumw2 = true, bool  sbox = true  )

private

Definition at line 726 of file QcdUeDQM.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), MillePedeFileConverter_cfg::e, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, hltTrgUsedNames_, mps_fire::i, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms().

                                      {
  // Book Profile histos.

  for (size_t i = 0; i < hltTrgUsedNames_.size(); ++i) {
    std::string folderName = "Physics/QcdUe/" + hltTrgUsedNames_.at(i);
    iBooker.setCurrentFolder(folderName);
    MonitorElement *e = iBooker.bookProfile(Form("%s_%s", name.c_str(), hltTrgUsedNames_.at(i).c_str()),
                                            Form("%s: %s", hltTrgUsedNames_.at(i).c_str(), title.c_str()),
                                            nx,
                                            x1,
                                            x2,
                                            y1,
                                            y2,
                                            " ");
    mes.push_back(e);
  }
}

void QcdUeDQM::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  iSetup  )

override

Definition at line 106 of file QcdUeDQM.cc.

References gather_cfg::cout, CP, newFWLiteAna::found, hltConfig, hltProcNames_, hltResName_, hltTrgBits_, hltTrgDeci_, hltTrgNames_, hltTrgUsedNames_, hltUsedResName_, mps_fire::i, HLTConfigProvider::init(), isHltConfigSuccessful_, dqmiolumiharvest::j, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().

                                                                   {
  // indicating change of HLT cfg at run boundries
  // for HLTConfigProvider::init()
  bool isHltCfgChange = false;
  isHltConfigSuccessful_ = false;  // init

  string teststr;
  for (size_t i = 0; i < hltProcNames_.size(); ++i) {
    if (i > 0)
      teststr += ", ";
    teststr += hltProcNames_.at(i);
    if (hltConfig.init(run, iSetup, hltProcNames_.at(i), isHltCfgChange)) {
      isHltConfigSuccessful_ = true;
      hltUsedResName_ = hltResName_;
      if (hltResName_.find(':') == string::npos)
        hltUsedResName_ += "::";
      else
        hltUsedResName_ += ":";
      hltUsedResName_ += hltProcNames_.at(i);
      break;
    }
  }

  if (!isHltConfigSuccessful_)
    return;

  // setup "Any" bit
  hltTrgBits_.clear();
  hltTrgBits_.push_back(-1);
  hltTrgDeci_.clear();
  hltTrgDeci_.push_back(true);
  hltTrgUsedNames_.clear();
  hltTrgUsedNames_.push_back("Any");

  // figure out relation of trigger name to trigger bit and store used trigger
  // names/bits
  for (size_t i = 0; i < hltTrgNames_.size(); ++i) {
    const string &n1(hltTrgNames_.at(i));
    // unsigned int hlt_prescale = hltConfig.prescaleValue(iSetup, n1);
    // cout<<"trigger=="<<n1<<"presc=="<<hlt_prescale<<endl;
    bool found = false;
    for (size_t j = 0; j < hltConfig.size(); ++j) {
      const string &n2(hltConfig.triggerName(j));
      if (n2 == n1) {
        hltTrgBits_.push_back(j);
        hltTrgUsedNames_.push_back(n1);
        hltTrgDeci_.push_back(false);
        found = true;
        break;
      }
    }
    if (!found) {
      CP(2) cout << "Could not find trigger bit" << endl;
    }
  }
  isHltConfigSuccessful_ = true;
}

void QcdUeDQM::fill1D ( std::vector< TH1F * > &  hs, double  val, double  w = 1.  )

private

Definition at line 754 of file QcdUeDQM.cc.

References hltTrgDeci_, mps_fire::i, heppy_batch::val, and w.

Referenced by analyze(), fillChargedJetSpectra(), fillpTMaxRelated(), fillUE_with_ChargedJets(), fillUE_with_MaxpTtrack(), fillVtxPlots(), and trackSelection().

                                                                 {
  // Loop over histograms and fill if trigger has fired.

  for (size_t i = 0; i < hs.size(); ++i) {
    if (!hltTrgDeci_.at(i))
      continue;
    hs.at(i)->Fill(val, w);
  }
}

void QcdUeDQM::fill1D ( std::vector< MonitorElement * > &  mes, double  val, double  w = 1.  )

private

void QcdUeDQM::fill2D ( std::vector< TH2F * > &  hs, double  valx, double  valy, double  w = 1.  )

private

Definition at line 788 of file QcdUeDQM.cc.

References hltTrgDeci_, mps_fire::i, and w.

                                                                               {
  // Loop over histograms and fill if trigger has fired.

  for (size_t i = 0; i < hs.size(); ++i) {
    if (!hltTrgDeci_.at(i))
      continue;
    hs.at(i)->Fill(valx, valy, w);
  }
}

void QcdUeDQM::fill2D ( std::vector< MonitorElement * > &  mes, double  valx, double  valy, double  w = 1.  )

private

void QcdUeDQM::fill3D ( std::vector< TH3F * > &  hs, int  gbin, double  w = 1.  )

private

void QcdUeDQM::fillChargedJetSpectra ( const edm::Handle< reco::TrackJetCollection >  trackJets )

private

Definition at line 976 of file QcdUeDQM.cc.

References f, fill1D(), hChargedJetMulti_, hLeadingChargedJet_etaSpectrum_, hLeadingChargedJet_phiSpectrum_, and hLeadingChargedJet_pTSpectrum_.

Referenced by analyze().

                                                                                      {
  fill1D(hChargedJetMulti_, trackJets->size());
  for (reco::TrackJetCollection::const_iterator f = trackJets->begin(); f != trackJets->end(); f++) {
    if (f != trackJets->begin())
      continue;
    fill1D(hLeadingChargedJet_pTSpectrum_, f->pt());
    fill1D(hLeadingChargedJet_etaSpectrum_, f->eta());
    fill1D(hLeadingChargedJet_phiSpectrum_, f->phi());
  }
}

void QcdUeDQM::fillHltBits ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 1251 of file QcdUeDQM.cc.

References edm::HLTGlobalStatus::accept(), dqm::impl::MonitorElement::Fill(), getProduct(), h2TrigCorr_, hltTrgBits_, hltTrgDeci_, hltUsedResName_, mps_fire::i, dqmiolumiharvest::j, and findQualityFiles::size.

Referenced by analyze().

                                                                        {
  // Fill HLT trigger bits.

  Handle<TriggerResults> triggerResultsHLT;
  getProduct(hltUsedResName_, triggerResultsHLT, iEvent);

  /*  const unsigned int ntrigs(triggerResultsHLT.product()->size());
    if( ntrigs != 0 ) {
    const edm::TriggerNames & triggerNames =
  iEvent.triggerNames(*triggerResultsHLT);
    for (unsigned int j=0; j!=ntrigs; j++) {
  //    if(triggerResultsHLT->accept(j)){
  //    unsigned int hlt_prescale = hltConfig_.prescaleValue(iEvent, iSetup,
  triggerName);
      cout<<"trigger fired"<<triggerNames.triggerName(j)<<endl;
      for(unsigned int itrigName = 0; itrigName < hltTrgNames_.size();
  itrigName++ ) {
       unsigned int hlt_prescale = hltConfig.prescaleValue(iEvent, iSetup,
  hltTrgNames_[itrigName]);

       if(triggerNames.triggerIndex(hltTrgNames_[itrigName]) >= (unsigned
  int)ntrigs ) continue;
  //    if(
  triggerResultsHLT->accept(triggerNames.triggerIndex(hltTrgNames_[itrigName]))
  )cout<<hltTrgNames_[itrigName]<<endl;

   }
      }
     }
   // }
  */
  for (size_t i = 0; i < hltTrgBits_.size(); ++i) {
    if (hltTrgBits_.at(i) < 0)
      continue;  // ignore unknown trigger
    size_t tbit = hltTrgBits_.at(i);
    if (tbit < triggerResultsHLT->size()) {
      hltTrgDeci_[i] = triggerResultsHLT->accept(tbit);
    }
  }
  // fill correlation histogram
  for (size_t i = 0; i < hltTrgBits_.size(); ++i) {
    if (hltTrgDeci_.at(i))
      h2TrigCorr_->Fill(i, i);
    for (size_t j = i + 1; j < hltTrgBits_.size(); ++j) {
      if (hltTrgDeci_.at(i) && hltTrgDeci_.at(j))
        h2TrigCorr_->Fill(i, j);
      if (hltTrgDeci_.at(i) && !hltTrgDeci_.at(j))
        h2TrigCorr_->Fill(j, i);
    }
  }
}

void QcdUeDQM::fillProfile ( std::vector< TProfile * > &  hs, double  valx, double  valy, double  w = 1.  )

private

Definition at line 809 of file QcdUeDQM.cc.

References hltTrgDeci_, mps_fire::i, w, and y.

Referenced by fillUE_with_ChargedJets(), and fillUE_with_MaxpTtrack().

                                                                                        {
  // Loop over histograms and fill if trigger has fired.

  for (size_t i = 0; i < hs.size(); ++i) {
    if (!hltTrgDeci_.at(i))
      continue;
    hs.at(i)->Fill(valx, valy, w);
  }
}

void QcdUeDQM::fillProfile ( std::vector< MonitorElement * > &  mes, double  valx, double  valy, double  w = 1.  )

private

void QcdUeDQM::fillpTMaxRelated ( const std::vector< const reco::Track * > &  track )

private

Definition at line 957 of file QcdUeDQM.cc.

References PVValHelper::eta, fill1D(), hGoodTrkEta500_, hGoodTrkEta900_, hGoodTrkPhi500_, hGoodTrkPhi900_, hGoodTrkPt500_, hGoodTrkPt900_, hLeadingTrack_etaSpectrum_, hLeadingTrack_phiSpectrum_, hLeadingTrack_pTSpectrum_, hNgoodTrk_, mps_fire::i, phi, and DiDispStaMuonMonitor_cfi::pt.

Referenced by analyze().

                                                                         {
  fill1D(hNgoodTrk_, track.size());
  if (!track.empty()) {
    fill1D(hLeadingTrack_pTSpectrum_, track[0]->pt());
    fill1D(hLeadingTrack_phiSpectrum_, track[0]->phi());
    fill1D(hLeadingTrack_etaSpectrum_, track[0]->eta());
  }
  for (size_t i = 0; i < track.size(); i++) {
    fill1D(hGoodTrkPt500_, track[i]->pt());
    fill1D(hGoodTrkEta500_, track[i]->eta());
    fill1D(hGoodTrkPhi500_, track[i]->phi());
    if (track[i]->pt() > 0.9) {
      fill1D(hGoodTrkPt900_, track[i]->pt());
      fill1D(hGoodTrkEta900_, track[i]->eta());
      fill1D(hGoodTrkPhi900_, track[i]->phi());
    }
  }
}

void QcdUeDQM::fillUE_with_ChargedJets ( const std::vector< const reco::Track * > &  track, const edm::Handle< reco::TrackJetCollection > &  trackJets  )

private

Definition at line 1082 of file QcdUeDQM.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, fill1D(), fillProfile(), hdNdEtadPhi_trackJet_Away500_, hdNdEtadPhi_trackJet_Away900_, hdNdEtadPhi_trackJet_Toward500_, hdNdEtadPhi_trackJet_Toward900_, hdNdEtadPhi_trackJet_Transverse500_, hdNdEtadPhi_trackJet_Transverse900_, hdPhi_chargedJet_tracks_, hpTSumdEtadPhi_trackJet_Away500_, hpTSumdEtadPhi_trackJet_Away900_, hpTSumdEtadPhi_trackJet_Toward500_, hpTSumdEtadPhi_trackJet_Toward900_, hpTSumdEtadPhi_trackJet_Transverse500_, hpTSumdEtadPhi_trackJet_Transverse900_, mps_fire::i, reco::btau::jetPhi, phi, PI, DiDispStaMuonMonitor_cfi::pt, and w.

Referenced by analyze().

                                                                                           {
  double w = 0.119;
  double nTrk500_TransReg = 0;
  double nTrk500_AwayReg = 0;
  double nTrk500_TowardReg = 0;

  double pTSum500_TransReg = 0;
  double pTSum500_AwayReg = 0;
  double pTSum500_TowardReg = 0;

  double nTrk900_TransReg = 0;
  double nTrk900_AwayReg = 0;
  double nTrk900_TowardReg = 0;

  double pTSum900_TransReg = 0;
  double pTSum900_AwayReg = 0;
  double pTSum900_TowardReg = 0;

  if (!(trackJets->empty()) && (trackJets->begin())->pt() > 1.) {
    double jetPhi = (trackJets->begin())->phi();
    for (size_t i = 0; i < track.size(); i++) {
      double dphi = (180. / PI) * (deltaPhi(jetPhi, track[i]->phi()));
      fill1D(hdPhi_chargedJet_tracks_, dphi);
      if (fabs(dphi) > 60. && fabs(dphi) < 120.) {
        pTSum500_TransReg = pTSum500_TransReg + track[i]->pt();
        nTrk500_TransReg++;
        if (track[i]->pt() > 0.9) {
          pTSum900_TransReg = pTSum900_TransReg + track[i]->pt();
          nTrk900_TransReg++;
        }
      }

      if (fabs(dphi) > 120. && fabs(dphi) < 180.) {
        pTSum500_AwayReg = pTSum500_AwayReg + track[i]->pt();
        nTrk500_AwayReg++;
        if (track[i]->pt() > 0.9) {
          pTSum900_AwayReg = pTSum900_AwayReg + track[i]->pt();
          nTrk900_AwayReg++;
        }
      }
      if (fabs(dphi) < 60.) {
        pTSum500_TowardReg = pTSum500_TowardReg + track[i]->pt();
        nTrk500_TowardReg++;
        if (track[i]->pt() > 0.9) {
          pTSum900_TowardReg = pTSum900_TowardReg + track[i]->pt();
          nTrk900_TowardReg++;
        }
      }
    }  // tracks loop

  }  // leading track jet

  fillProfile(hdNdEtadPhi_trackJet_Toward500_, (trackJets->begin())->pt(), nTrk500_TowardReg, w);
  fillProfile(hdNdEtadPhi_trackJet_Transverse500_, (trackJets->begin())->pt(), nTrk500_TransReg, w);
  fillProfile(hdNdEtadPhi_trackJet_Away500_, (trackJets->begin())->pt(), nTrk500_AwayReg, w);

  fillProfile(hpTSumdEtadPhi_trackJet_Toward500_, (trackJets->begin())->pt(), pTSum500_TowardReg, w);
  fillProfile(hpTSumdEtadPhi_trackJet_Transverse500_, (trackJets->begin())->pt(), pTSum500_TransReg, w);
  fillProfile(hpTSumdEtadPhi_trackJet_Away500_, (trackJets->begin())->pt(), pTSum500_AwayReg, w);

  fillProfile(hdNdEtadPhi_trackJet_Toward900_, (trackJets->begin())->pt(), nTrk900_TowardReg, w);
  fillProfile(hdNdEtadPhi_trackJet_Transverse900_, (trackJets->begin())->pt(), nTrk900_TransReg, w);
  fillProfile(hdNdEtadPhi_trackJet_Away900_, (trackJets->begin())->pt(), nTrk900_AwayReg, w);

  fillProfile(hpTSumdEtadPhi_trackJet_Toward900_, (trackJets->begin())->pt(), pTSum900_TowardReg, w);
  fillProfile(hpTSumdEtadPhi_trackJet_Transverse900_, (trackJets->begin())->pt(), pTSum900_TransReg, w);
  fillProfile(hpTSumdEtadPhi_trackJet_Away900_, (trackJets->begin())->pt(), pTSum900_AwayReg, w);
}

void QcdUeDQM::fillUE_with_MaxpTtrack ( const std::vector< const reco::Track * > &  track )

private

Definition at line 1012 of file QcdUeDQM.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, fill1D(), fillProfile(), hdNdEtadPhi_pTMax_Away500_, hdNdEtadPhi_pTMax_Away900_, hdNdEtadPhi_pTMax_Toward500_, hdNdEtadPhi_pTMax_Toward900_, hdNdEtadPhi_pTMax_Transverse500_, hdNdEtadPhi_pTMax_Transverse900_, hdPhi_maxpTTrack_tracks_, hpTSumdEtadPhi_pTMax_Away500_, hpTSumdEtadPhi_pTMax_Away900_, hpTSumdEtadPhi_pTMax_Toward500_, hpTSumdEtadPhi_pTMax_Toward900_, hpTSumdEtadPhi_pTMax_Transverse500_, hpTSumdEtadPhi_pTMax_Transverse900_, mps_fire::i, phi, PI, DiDispStaMuonMonitor_cfi::pt, and w.

Referenced by analyze().

                                                                               {
  double w = 0.119;
  // double w = 1.;
  double nTrk500_TransReg = 0;
  double nTrk500_AwayReg = 0;
  double nTrk500_TowardReg = 0;

  double pTSum500_TransReg = 0;
  double pTSum500_AwayReg = 0;
  double pTSum500_TowardReg = 0;

  double nTrk900_TransReg = 0;
  double nTrk900_AwayReg = 0;
  double nTrk900_TowardReg = 0;

  double pTSum900_TransReg = 0;
  double pTSum900_AwayReg = 0;
  double pTSum900_TowardReg = 0;
  if (!track.empty()) {
    if (track[0]->pt() > 1.) {
      for (size_t i = 1; i < track.size(); i++) {
        double dphi = (180. / PI) * (deltaPhi(track[0]->phi(), track[i]->phi()));
        fill1D(hdPhi_maxpTTrack_tracks_, dphi);
        if (fabs(dphi) > 60. && fabs(dphi) < 120.) {
          pTSum500_TransReg = pTSum500_TransReg + track[i]->pt();
          nTrk500_TransReg++;
          if (track[i]->pt() > 0.9) {
            pTSum900_TransReg = pTSum900_TransReg + track[i]->pt();
            nTrk900_TransReg++;
          }
        }

        if (fabs(dphi) > 120. && fabs(dphi) < 180.) {
          pTSum500_AwayReg = pTSum500_AwayReg + track[i]->pt();
          nTrk500_AwayReg++;
          if (track[i]->pt() > 0.9) {
            pTSum900_AwayReg = pTSum900_AwayReg + track[i]->pt();
            nTrk900_AwayReg++;
          }
        }

        if (fabs(dphi) < 60.) {
          pTSum500_TowardReg = pTSum500_TowardReg + track[i]->pt();
          nTrk500_TowardReg++;
          if (track[i]->pt() > 0.9) {
            pTSum900_TowardReg = pTSum900_TowardReg + track[i]->pt();
            nTrk900_TowardReg++;
          }
        }
      }  // track loop
    }    // leading track
         // non empty collection
    fillProfile(hdNdEtadPhi_pTMax_Toward500_, track[0]->pt(), nTrk500_TowardReg, w);
    fillProfile(hdNdEtadPhi_pTMax_Transverse500_, track[0]->pt(), nTrk500_TransReg, w);
    fillProfile(hdNdEtadPhi_pTMax_Away500_, track[0]->pt(), nTrk500_AwayReg, w);

    fillProfile(hpTSumdEtadPhi_pTMax_Toward500_, track[0]->pt(), pTSum500_TowardReg, w);
    fillProfile(hpTSumdEtadPhi_pTMax_Transverse500_, track[0]->pt(), pTSum500_TransReg, w);
    fillProfile(hpTSumdEtadPhi_pTMax_Away500_, track[0]->pt(), pTSum500_AwayReg, w);

    fillProfile(hdNdEtadPhi_pTMax_Toward900_, track[0]->pt(), nTrk900_TowardReg, w);
    fillProfile(hdNdEtadPhi_pTMax_Transverse900_, track[0]->pt(), nTrk900_TransReg, w);
    fillProfile(hdNdEtadPhi_pTMax_Away900_, track[0]->pt(), nTrk900_AwayReg, w);

    fillProfile(hpTSumdEtadPhi_pTMax_Toward900_, track[0]->pt(), pTSum900_TowardReg, w);
    fillProfile(hpTSumdEtadPhi_pTMax_Transverse900_, track[0]->pt(), pTSum900_TransReg, w);
    fillProfile(hpTSumdEtadPhi_pTMax_Away900_, track[0]->pt(), pTSum900_AwayReg, w);
  }
}

bool QcdUeDQM::fillVtxPlots ( const reco::BeamSpot *  bs, const edm::Handle< reco::VertexCollection >  vtxColl  )

private

Definition at line 934 of file QcdUeDQM.cc.

References diffvtxbs_, fill1D(), hNvertices_, hVertex_ndof_, hVertex_rho_, hVertex_x_, hVertex_y_, hVertex_z_, hVertex_z_bs_, edm::Handle< T >::product(), vtx1, and reco::BeamSpot::z0().

Referenced by analyze().

                                                                                                 {
  const reco::VertexCollection theVertices = *(vtxColl.product());
  bool goodVtx = false;
  fill1D(hNvertices_, theVertices.size());
  for (reco::VertexCollection::const_iterator vertexIt = theVertices.begin(); vertexIt != theVertices.end();
       ++vertexIt) {
    fill1D(hVertex_z_, vertexIt->z());
    fill1D(hVertex_y_, vertexIt->y());
    fill1D(hVertex_x_, vertexIt->x());
    fill1D(hVertex_ndof_, vertexIt->ndof());
    fill1D(hVertex_rho_, vertexIt->position().rho());
    fill1D(hVertex_z_bs_, (vertexIt->z() - bs->z0()));

    if (fabs(vertexIt->z() - bs->z0()) < diffvtxbs_ && vertexIt->ndof() >= 4 && vertexIt->position().rho() <= 2.0) {
      goodVtx = true;
      vtx1 = (*vertexIt);

      break;
    }
  }  // Loop over vertcies
  return goodVtx;
}

template<typename TYPE >

void QcdUeDQM::getProduct ( const std::string  name, edm::Handle< TYPE > &  prod, const edm::Event &  event  ) const

inlineprivate

Definition at line 405 of file QcdUeDQM.h.

References edm::errors::Configuration, and edm::HandleBase::isValid().

Referenced by fillHltBits().

                                                                                                         {
  // Try to access data collection from EDM file. We check if we really get just
  // one
  // product with the given name. If not we throw an exception.

  event.getByLabel(edm::InputTag(name), prod);
  if (!prod.isValid())
    throw edm::Exception(edm::errors::Configuration, "QcdUeDQM::GetProduct()\n")
        << "Collection with label " << name << " is not valid" << std::endl;
}

template<typename TYPE >

bool QcdUeDQM::getProductSafe ( const std::string  name, edm::Handle< TYPE > &  prod, const edm::Event &  event  ) const

inlineprivate

Definition at line 418 of file QcdUeDQM.h.

References edm::HandleBase::isValid().

                                                                                                             {
  // Try to safely access data collection from EDM file. We check if we really
  // get just one
  // product with the given name. If not, we return false.

  if (name.empty())
    return false;

  try {
    event.getByLabel(edm::InputTag(name), prod);
    if (!prod.isValid())
      return false;
  } catch (...) {
    return false;
  }
  return true;
}

void QcdUeDQM::setLabel1D ( std::vector< MonitorElement * > &  mes )

private

Definition at line 776 of file QcdUeDQM.cc.

References TkAlMuonSelectors_cfi::cut, hltTrgDeci_, mps_fire::i, and dqmiolumiharvest::j.

                                                          {
  // Loop over histograms and fill if trigger has fired.
  string cut[5] = {"Nevt", "vtx!=bmspt", "Zvtx<10cm", "pT>1GeV", "trackFromVtx"};
  for (size_t i = 0; i < mes.size(); ++i) {
    if (!hltTrgDeci_.at(i))
      continue;
    for (size_t j = 1; j < 6; j++)
      mes.at(i)->setBinLabel(j, cut[j - 1], 1);
  }
}

bool QcdUeDQM::trackSelection ( const reco::Track &  trk, const reco::BeamSpot *  bs, const reco::Vertex &  vtx, int  sizevtx  )

private

Definition at line 832 of file QcdUeDQM.cc.

References reco::TrackBase::algo(), algorithm_, allowTriplets_, bsuse_, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::TrackBase::eta(), fill1D(), spr::find(), reco::HitPattern::getHitPattern(), reco::HitPattern::getLayer(), hBeamSpot_x_, hBeamSpot_y_, hBeamSpot_z_, reco::TrackBase::hitPattern(), hNtrackerLayer_, hNtrackerPixelLayer_, hNtrackerStripPixelLayer_, hRatioDxySigmaDxyBS_, hRatioDxySigmaDxyPV_, hRatioDzSigmaDzBS_, hRatioDzSigmaDzPV_, hRatioPtErrorPt_, hTrkChi2_, hTrkEta_, hTrkNdof_, hTrkPhi_, hTrkPt_, mps_fire::i, lip_, maxChi2_, maxRapidity_, min3DHit_, minHit_, minRapidity_, reco::TrackBase::ndof(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfAllHits(), reco::HitPattern::numberOfValidStripLayersWithMonoAndStereo(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::phi(), reco::HitPattern::pixelHitFilter(), reco::HitPattern::pixelLayersWithMeasurement(), reco::BeamSpot::position(), reco::Vertex::position(), reco::TrackBase::pt(), ptErr_pt_, reco::TrackBase::ptError(), ptMin_, pxlLayerMinCut_, reco::TrackBase::quality(), quality_, requirePIX1_, tip_, reco::HitPattern::TRACK_HITS, reco::HitPattern::trackerLayersWithMeasurement(), reco::HitPattern::validHitFilter(), reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by analyze().

                                                                                                              {
  bool goodTrk = false;

  if (sizevtx != 1)
    return false;  // selection events with only a vertex

  // Fill basic information of all the tracks
  fill1D(hNtrackerLayer_, trk.hitPattern().trackerLayersWithMeasurement());
  fill1D(hNtrackerPixelLayer_, trk.hitPattern().pixelLayersWithMeasurement());

  fill1D(
      hNtrackerStripPixelLayer_,
      (trk.hitPattern().pixelLayersWithMeasurement() + trk.hitPattern().numberOfValidStripLayersWithMonoAndStereo()));

  fill1D(hRatioPtErrorPt_, (trk.ptError() / trk.pt()));
  fill1D(hTrkPt_, trk.pt());
  fill1D(hTrkEta_, trk.eta());
  fill1D(hTrkPhi_, trk.phi());
  fill1D(hRatioDxySigmaDxyBS_, (trk.dxy(bs->position()) / trk.dxyError()));
  fill1D(hRatioDxySigmaDxyPV_, (trk.dxy(vtx.position()) / trk.dxyError()));
  fill1D(hRatioDzSigmaDzBS_, (trk.dz(bs->position()) / trk.dzError()));
  fill1D(hRatioDzSigmaDzPV_, (trk.dz(vtx.position()) / trk.dzError()));
  fill1D(hTrkChi2_, trk.normalizedChi2());
  fill1D(hTrkNdof_, trk.ndof());

  fill1D(hBeamSpot_x_, bs->x0());
  fill1D(hBeamSpot_y_, bs->y0());
  fill1D(hBeamSpot_z_, bs->z0());

  // number of layers
  bool layerMinCutbool = (trk.hitPattern().trackerLayersWithMeasurement() >= minHit_ ||
                          (trk.hitPattern().trackerLayersWithMeasurement() == 3 &&
                           trk.hitPattern().pixelLayersWithMeasurement() == 3 && allowTriplets_));

  // number of pixel layers
  bool pxlLayerMinCutbool = (trk.hitPattern().pixelLayersWithMeasurement() >= pxlLayerMinCut_);

  // cut on the hits in pixel layers
  bool hasPIX1 = false;
  if (requirePIX1_) {
    const reco::HitPattern &p = trk.hitPattern();
    for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
      uint32_t hit = p.getHitPattern(reco::HitPattern::TRACK_HITS, i);
      if (reco::HitPattern::validHitFilter(hit) && reco::HitPattern::pixelHitFilter(hit) &&
          reco::HitPattern::getLayer(hit) == 1) {
        hasPIX1 = true;
        break;
      }
    }
  } else {
    hasPIX1 = true;
  }

  // cut on the pT error
  bool ptErrorbool =
      (trk.ptError() / trk.pt() < ptErr_pt_ || (trk.hitPattern().trackerLayersWithMeasurement() == 3 &&
                                                trk.hitPattern().pixelLayersWithMeasurement() == 3 && allowTriplets_));

  // quality cut
  bool quality_ok = true;
  if (!quality_.empty()) {
    quality_ok = false;
    for (unsigned int i = 0; i < quality_.size(); ++i) {
      if (trk.quality(quality_[i])) {
        quality_ok = true;
        break;
      }
    }
  }
  //-----
  bool algo_ok = true;
  if (!algorithm_.empty()) {
    if (std::find(algorithm_.begin(), algorithm_.end(), trk.algo()) == algorithm_.end())
      algo_ok = false;
  }

  if (bsuse_ == 1) {
    if (hasPIX1 && pxlLayerMinCutbool && layerMinCutbool &&
        (trk.hitPattern().pixelLayersWithMeasurement() +
         trk.hitPattern().numberOfValidStripLayersWithMonoAndStereo()) >= min3DHit_ &&
        ptErrorbool && fabs(trk.pt()) >= ptMin_ && trk.eta() >= minRapidity_ && trk.eta() <= maxRapidity_ &&
        fabs(trk.dxy(bs->position()) / trk.dxyError()) < tip_ && fabs(trk.dz(bs->position()) / trk.dzError()) < lip_ &&
        trk.normalizedChi2() <= maxChi2_ && quality_ok && algo_ok) {
      goodTrk = true;
    }
  }

  if (bsuse_ == 0) {
    if (hasPIX1 && pxlLayerMinCutbool && layerMinCutbool &&
        (trk.hitPattern().pixelLayersWithMeasurement() +
         trk.hitPattern().numberOfValidStripLayersWithMonoAndStereo()) >= min3DHit_ &&
        ptErrorbool && fabs(trk.pt()) >= ptMin_ && trk.eta() >= minRapidity_ && trk.eta() <= maxRapidity_ &&
        fabs(trk.dxy(vtx.position()) / trk.dxyError()) < tip_ && fabs(trk.dz(vtx.position()) / trk.dzError()) < lip_ &&
        trk.normalizedChi2() <= maxChi2_ && quality_ok && algo_ok) {
      goodTrk = true;
    }
  }

  return goodTrk;
}

Member Data Documentation

std::vector<reco::TrackBase::TrackAlgorithm> QcdUeDQM::algorithm_

private

Definition at line 397 of file QcdUeDQM.h.

Referenced by QcdUeDQM(), and trackSelection().

bool QcdUeDQM::allowTriplets_

private

Definition at line 389 of file QcdUeDQM.h.

Referenced by trackSelection().

edm::EDGetTokenT<reco::BeamSpot> QcdUeDQM::bsLabel_

private

Definition at line 395 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::bsPos_

private

Definition at line 390 of file QcdUeDQM.h.

bool QcdUeDQM::bsuse_

private

Definition at line 388 of file QcdUeDQM.h.

Referenced by trackSelection().

edm::EDGetTokenT<reco::CaloJetCollection> QcdUeDQM::caloJetLabel_

private

Definition at line 391 of file QcdUeDQM.h.

Referenced by analyze().

edm::EDGetTokenT<reco::TrackJetCollection> QcdUeDQM::chargedJetLabel_

private

Definition at line 392 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::diffvtxbs_

private

Definition at line 380 of file QcdUeDQM.h.

Referenced by fillVtxPlots().

MonitorElement* QcdUeDQM::h2TrigCorr_

private

Definition at line 121 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillHltBits().

std::vector<MonitorElement *> QcdUeDQM::hBeamSpot_x_

private

Definition at line 155 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hBeamSpot_y_

private

Definition at line 154 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hBeamSpot_z_

private

Definition at line 153 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hChargedJetConstituent_

private

Definition at line 165 of file QcdUeDQM.h.

std::vector<MonitorElement *> QcdUeDQM::hChargedJetMulti_

private

Definition at line 164 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Away500_

private

Definition at line 214 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Away900_

private

Definition at line 301 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Toward500_

private

Definition at line 201 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Toward900_

private

Definition at line 288 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Transverse500_

private

Definition at line 207 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Transverse900_

private

Definition at line 294 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Away500_

private

Definition at line 235 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Away900_

private

Definition at line 322 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Toward500_

private

Definition at line 225 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Toward900_

private

Definition at line 312 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Transverse500_

private

Definition at line 232 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Transverse900_

private

Definition at line 319 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdPhi_chargedJet_tracks_

private

Definition at line 188 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hdPhi_maxpTTrack_tracks_

private

Definition at line 183 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkEta500_

private

Definition at line 133 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkEta900_

private

Definition at line 136 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkPhi500_

private

Definition at line 134 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkPhi900_

private

Definition at line 137 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkPt500_

private

Definition at line 132 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hGoodTrkPt900_

private

Definition at line 135 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_etaSpectrum_

private

Definition at line 192 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_phiSpectrum_

private

Definition at line 196 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_pTSpectrum_

private

Definition at line 168 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

std::vector<MonitorElement *> QcdUeDQM::hLeadingTrack_etaSpectrum_

private

Definition at line 159 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hLeadingTrack_phiSpectrum_

private

Definition at line 161 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hLeadingTrack_pTSpectrum_

private

Definition at line 157 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

HLTConfigProvider QcdUeDQM::hltConfig

private

Definition at line 107 of file QcdUeDQM.h.

Referenced by dqmBeginRun().

std::vector<std::string> QcdUeDQM::hltProcNames_

private

Definition at line 110 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and QcdUeDQM().

std::string QcdUeDQM::hltResName_

private

Definition at line 109 of file QcdUeDQM.h.

Referenced by dqmBeginRun().

std::vector<int> QcdUeDQM::hltTrgBits_

private

Definition at line 113 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and fillHltBits().

std::vector<bool> QcdUeDQM::hltTrgDeci_

private

Definition at line 114 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), fill1D(), fill2D(), fillHltBits(), fillProfile(), and setLabel1D().

std::vector<std::string> QcdUeDQM::hltTrgNames_

private

Definition at line 111 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and QcdUeDQM().

std::vector<std::string> QcdUeDQM::hltTrgUsedNames_

private

Definition at line 115 of file QcdUeDQM.h.

Referenced by book1D(), bookHistograms(), bookProfile(), and dqmBeginRun().

std::string QcdUeDQM::hltUsedResName_

private

Definition at line 116 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and fillHltBits().

std::vector<MonitorElement *> QcdUeDQM::hNevts_

private

Definition at line 123 of file QcdUeDQM.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> QcdUeDQM::hNgoodTrk_

private

Definition at line 131 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

std::vector<MonitorElement *> QcdUeDQM::hNtrackerLayer_

private

Definition at line 124 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hNtrackerPixelLayer_

private

Definition at line 125 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hNtrackerStripPixelLayer_

private

Definition at line 126 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hNvertices_

private

Definition at line 145 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Away500_

private

Definition at line 252 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Away900_

private

Definition at line 339 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Toward500_

private

Definition at line 242 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Toward900_

private

Definition at line 329 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Transverse500_

private

Definition at line 249 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Transverse900_

private

Definition at line 336 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Away500_

private

Definition at line 279 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Away900_

private

Definition at line 366 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Toward500_

private

Definition at line 265 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Toward900_

private

Definition at line 352 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Transverse500_

private

Definition at line 276 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Transverse900_

private

Definition at line 363 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

std::vector<MonitorElement *> QcdUeDQM::hRatioDxySigmaDxyBS_

private

Definition at line 138 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hRatioDxySigmaDxyPV_

private

Definition at line 139 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hRatioDzSigmaDzBS_

private

Definition at line 140 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hRatioDzSigmaDzPV_

private

Definition at line 141 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hRatioPtErrorPt_

private

Definition at line 127 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hTrkChi2_

private

Definition at line 142 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hTrkEta_

private

Definition at line 129 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hTrkNdof_

private

Definition at line 143 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hTrkPhi_

private

Definition at line 130 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hTrkPt_

private

Definition at line 128 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

std::vector<MonitorElement *> QcdUeDQM::hVertex_ndof_

private

Definition at line 149 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hVertex_rho_

private

Definition at line 150 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hVertex_x_

private

Definition at line 148 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hVertex_y_

private

Definition at line 147 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hVertex_z_

private

Definition at line 146 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

std::vector<MonitorElement *> QcdUeDQM::hVertex_z_bs_

private

Definition at line 151 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

bool QcdUeDQM::isHltConfigSuccessful_

private

Definition at line 62 of file QcdUeDQM.h.

Referenced by analyze(), dqmBeginRun(), and QcdUeDQM().

double QcdUeDQM::lip_

private

Definition at line 379 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::maxChi2_

private

Definition at line 387 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::maxRapidity_

private

Definition at line 377 of file QcdUeDQM.h.

Referenced by trackSelection().

int QcdUeDQM::min3DHit_

private

Definition at line 386 of file QcdUeDQM.h.

Referenced by trackSelection().

int QcdUeDQM::minHit_

private

Definition at line 383 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::minRapidity_

private

Definition at line 376 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::ptErr_pt_

private

Definition at line 381 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::ptMin_

private

Definition at line 375 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::pxlLayerMinCut_

private

Definition at line 384 of file QcdUeDQM.h.

Referenced by trackSelection().

std::vector<reco::TrackBase::TrackQuality> QcdUeDQM::quality_

private

Definition at line 396 of file QcdUeDQM.h.

Referenced by QcdUeDQM(), and trackSelection().

MonitorElement* QcdUeDQM::repSumMap_

private

Definition at line 119 of file QcdUeDQM.h.

MonitorElement* QcdUeDQM::repSummary_

private

Definition at line 120 of file QcdUeDQM.h.

bool QcdUeDQM::requirePIX1_

private

Definition at line 385 of file QcdUeDQM.h.

Referenced by trackSelection().

container QcdUeDQM::selected_

private

Definition at line 399 of file QcdUeDQM.h.

Referenced by analyze().

const TrackerGeometry* QcdUeDQM::tgeo_

private

Definition at line 118 of file QcdUeDQM.h.

double QcdUeDQM::tip_

private

Definition at line 378 of file QcdUeDQM.h.

Referenced by trackSelection().

edm::EDGetTokenT<reco::TrackCollection> QcdUeDQM::trackLabel_

private

Definition at line 393 of file QcdUeDQM.h.

Referenced by analyze().

int QcdUeDQM::verbose_

private

Definition at line 117 of file QcdUeDQM.h.

reco::Vertex QcdUeDQM::vtx1

private

Definition at line 400 of file QcdUeDQM.h.

Referenced by analyze(), and fillVtxPlots().

edm::EDGetTokenT<reco::VertexCollection> QcdUeDQM::vtxLabel_

private

Definition at line 394 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::vtxntk_

private

Definition at line 382 of file QcdUeDQM.h.