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 &ev, edm::EventSetup const &es) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) override
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
	~DQMEDAnalyzer () override=default
Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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 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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 54 of file QcdUeDQM.h.

Member Typedef Documentation

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

private

Definition at line 402 of file QcdUeDQM.h.

Constructor & Destructor Documentation

QcdUeDQM::QcdUeDQM

(

const edm::ParameterSet &

parameters

)

Definition at line 53 of file QcdUeDQM.cc.

References reco::TrackBase::algoByName(), electronCleaner_cfi::algorithm, algorithm_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltProcNames_, hltTrgNames_, isHltConfigSuccessful_, jets_cff::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 111 of file QcdUeDQM.cc.

                    {
  // Destructor.
}

Member Function Documentation

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

override

Definition at line 444 of file QcdUeDQM.cc.

References ecalDrivenElectronSeedsParameters_cff::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_, l1t::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 &  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 498 of file QcdUeDQM.cc.

References DQMStore::IBooker::book1D(), MillePedeFileConverter_cfg::e, MonitorElement::getTH1(), hltTrgUsedNames_, mps_fire::i, 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 171 of file QcdUeDQM.cc.

References book1D(), DQMStore::IBooker::book2D(), bookProfile(), 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, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), globals_cff::x1, and globals_cff::x2.

                                                     {
  // 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 &  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 522 of file QcdUeDQM.cc.

References DQMStore::IBooker::bookProfile(), MillePedeFileConverter_cfg::e, hltTrgUsedNames_, mps_fire::i, 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 115 of file QcdUeDQM.cc.

References gather_cfg::cout, CP, runEdmFileComparison::found, hltConfig, hltProcNames_, hltResName_, hltTrgBits_, hltTrgDeci_, hltTrgNames_, hltTrgUsedNames_, hltUsedResName_, mps_fire::i, HLTConfigProvider::init(), isHltConfigSuccessful_, 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 540 of file QcdUeDQM.cc.

References hltTrgDeci_, and mps_fire::i.

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

Definition at line 550 of file QcdUeDQM.cc.

References hltTrgDeci_, and mps_fire::i.

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

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

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

private

Definition at line 572 of file QcdUeDQM.cc.

References hltTrgDeci_, and mps_fire::i.

                                {
  // 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

Definition at line 583 of file QcdUeDQM.cc.

References hltTrgDeci_, and mps_fire::i.

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

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

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 773 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 1074 of file QcdUeDQM.cc.

References edm::HLTGlobalStatus::accept(), MonitorElement::Fill(), getProduct(), h2TrigCorr_, hltTrgBits_, hltTrgDeci_, hltUsedResName_, mps_fire::i, 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 593 of file QcdUeDQM.cc.

References hltTrgDeci_, and mps_fire::i.

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

Definition at line 604 of file QcdUeDQM.cc.

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

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

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

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

private

Definition at line 754 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 EnergyCorrector::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 892 of file QcdUeDQM.cc.

References hiPixelPairStep_cff::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, EnergyCorrector::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 810 of file QcdUeDQM.cc.

References hiPixelPairStep_cff::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, EnergyCorrector::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 729 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 409 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 424 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 561 of file QcdUeDQM.cc.

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

                                                          {
  // 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 616 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 401 of file QcdUeDQM.h.

Referenced by QcdUeDQM(), and trackSelection().

bool QcdUeDQM::allowTriplets_

private

Definition at line 393 of file QcdUeDQM.h.

Referenced by trackSelection().

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

private

Definition at line 399 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::bsPos_

private

Definition at line 394 of file QcdUeDQM.h.

bool QcdUeDQM::bsuse_

private

Definition at line 392 of file QcdUeDQM.h.

Referenced by trackSelection().

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

private

Definition at line 395 of file QcdUeDQM.h.

Referenced by analyze().

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

private

Definition at line 396 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::diffvtxbs_

private

Definition at line 384 of file QcdUeDQM.h.

Referenced by fillVtxPlots().

MonitorElement* QcdUeDQM::h2TrigCorr_

private

Definition at line 119 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillHltBits().

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

private

Definition at line 153 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 152 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 151 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 163 of file QcdUeDQM.h.

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

private

Definition at line 162 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

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

private

Definition at line 212 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 302 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 199 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 289 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 205 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 295 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 234 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 324 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 223 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 313 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 231 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 321 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 186 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 181 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 131 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 134 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 132 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 135 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 130 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 133 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 190 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

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

private

Definition at line 194 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

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

private

Definition at line 166 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillChargedJetSpectra().

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

private

Definition at line 157 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 159 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 155 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

HLTConfigProvider QcdUeDQM::hltConfig

private

Definition at line 105 of file QcdUeDQM.h.

Referenced by dqmBeginRun().

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

private

Definition at line 108 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and QcdUeDQM().

std::string QcdUeDQM::hltResName_

private

Definition at line 107 of file QcdUeDQM.h.

Referenced by dqmBeginRun().

std::vector<int> QcdUeDQM::hltTrgBits_

private

Definition at line 111 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and fillHltBits().

std::vector<bool> QcdUeDQM::hltTrgDeci_

private

Definition at line 112 of file QcdUeDQM.h.

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

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

private

Definition at line 109 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and QcdUeDQM().

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

private

Definition at line 113 of file QcdUeDQM.h.

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

std::string QcdUeDQM::hltUsedResName_

private

Definition at line 114 of file QcdUeDQM.h.

Referenced by dqmBeginRun(), and fillHltBits().

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

private

Definition at line 121 of file QcdUeDQM.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 129 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillpTMaxRelated().

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

private

Definition at line 122 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 123 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 124 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 143 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 342 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 241 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 331 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 339 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_MaxpTtrack().

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

private

Definition at line 280 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 370 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 355 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 277 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 367 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillUE_with_ChargedJets().

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

private

Definition at line 136 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 137 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 138 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 139 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 125 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 140 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 127 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 141 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 128 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 126 of file QcdUeDQM.h.

Referenced by bookHistograms(), and trackSelection().

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

private

Definition at line 147 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

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

private

Definition at line 148 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

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

private

Definition at line 146 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

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

private

Definition at line 145 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

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

private

Definition at line 144 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

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

private

Definition at line 149 of file QcdUeDQM.h.

Referenced by bookHistograms(), and fillVtxPlots().

bool QcdUeDQM::isHltConfigSuccessful_

private

Definition at line 64 of file QcdUeDQM.h.

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

double QcdUeDQM::lip_

private

Definition at line 383 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::maxChi2_

private

Definition at line 391 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::maxRapidity_

private

Definition at line 381 of file QcdUeDQM.h.

Referenced by trackSelection().

int QcdUeDQM::min3DHit_

private

Definition at line 390 of file QcdUeDQM.h.

Referenced by trackSelection().

int QcdUeDQM::minHit_

private

Definition at line 387 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::minRapidity_

private

Definition at line 380 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::ptErr_pt_

private

Definition at line 385 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::ptMin_

private

Definition at line 379 of file QcdUeDQM.h.

Referenced by trackSelection().

double QcdUeDQM::pxlLayerMinCut_

private

Definition at line 388 of file QcdUeDQM.h.

Referenced by trackSelection().

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

private

Definition at line 400 of file QcdUeDQM.h.

Referenced by QcdUeDQM(), and trackSelection().

MonitorElement* QcdUeDQM::repSumMap_

private

Definition at line 117 of file QcdUeDQM.h.

MonitorElement* QcdUeDQM::repSummary_

private

Definition at line 118 of file QcdUeDQM.h.

bool QcdUeDQM::requirePIX1_

private

Definition at line 389 of file QcdUeDQM.h.

Referenced by trackSelection().

container QcdUeDQM::selected_

private

Definition at line 403 of file QcdUeDQM.h.

Referenced by analyze().

const TrackerGeometry* QcdUeDQM::tgeo_

private

Definition at line 116 of file QcdUeDQM.h.

double QcdUeDQM::tip_

private

Definition at line 382 of file QcdUeDQM.h.

Referenced by trackSelection().

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

private

Definition at line 397 of file QcdUeDQM.h.

Referenced by analyze().

int QcdUeDQM::verbose_

private

Definition at line 115 of file QcdUeDQM.h.

reco::Vertex QcdUeDQM::vtx1

private

Definition at line 404 of file QcdUeDQM.h.

Referenced by analyze(), and fillVtxPlots().

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

private

Definition at line 398 of file QcdUeDQM.h.

Referenced by analyze().

double QcdUeDQM::vtxntk_

private

Definition at line 386 of file QcdUeDQM.h.