TrackSplittingMonitor Class Reference

#include <DQM/TrackingMonitor/src/TrackSplittingMonitor.cc>

Inheritance diagram for TrackSplittingMonitor:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	TrackSplittingMonitor (const edm::ParameterSet &)
	~TrackSplittingMonitor () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	doProfileX (MonitorElement *th2m, MonitorElement *me)
void	doProfileX (TH2 *th2, MonitorElement *me)

Private Attributes
edm::ParameterSet	conf_
const CSCGeometry *	cscGeometry
const edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	cscGeomToken_
double	d0Cut_
MonitorElement *	dcurvAbsoluteResiduals_global_
MonitorElement *	dcurvAbsoluteResiduals_tracker_
MonitorElement *	dcurvNormalizedResiduals_global_
MonitorElement *	dcurvNormalizedResiduals_tracker_
MonitorElement *	ddxyAbsoluteResiduals_global_
MonitorElement *	ddxyAbsoluteResiduals_tracker_
MonitorElement *	ddxyNormalizedResiduals_global_
MonitorElement *	ddxyNormalizedResiduals_tracker_
MonitorElement *	ddzAbsoluteResiduals_global_
MonitorElement *	ddzAbsoluteResiduals_tracker_
MonitorElement *	ddzNormalizedResiduals_global_
MonitorElement *	ddzNormalizedResiduals_tracker_
MonitorElement *	dphiAbsoluteResiduals_global_
MonitorElement *	dphiAbsoluteResiduals_tracker_
MonitorElement *	dphiNormalizedResiduals_global_
MonitorElement *	dphiNormalizedResiduals_tracker_
MonitorElement *	dptAbsoluteResiduals_global_
MonitorElement *	dptAbsoluteResiduals_tracker_
MonitorElement *	dptNormalizedResiduals_global_
MonitorElement *	dptNormalizedResiduals_tracker_
DQMStore *	dqmStore_
const DTGeometry *	dtGeometry
const edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomToken_
MonitorElement *	dthetaAbsoluteResiduals_global_
MonitorElement *	dthetaAbsoluteResiduals_tracker_
MonitorElement *	dthetaNormalizedResiduals_global_
MonitorElement *	dthetaNormalizedResiduals_tracker_
double	dzCut_
std::string	histname
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	mfToken_
double	norchiCut_
int	pixelHitsPerLeg_
bool	plotMuons_
double	ptCut_
const RPCGeometry *	rpcGeometry
const edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeomToken_
edm::InputTag	splitMuons_
edm::EDGetTokenT< std::vector< reco::Muon > >	splitMuonsToken_
edm::InputTag	splitTracks_
edm::EDGetTokenT< std::vector< reco::Track > >	splitTracksToken_
const TrackerGeometry *	theGeometry
const MagneticField *	theMagField
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_
int	totalHitsPerLeg_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Monitoring source for general quantities related to tracks.

Definition at line 41 of file TrackSplittingMonitor.h.

Constructor & Destructor Documentation

TrackSplittingMonitor()

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

explicit

Definition at line 26 of file TrackSplittingMonitor.cc.

    : dqmStore_(edm::Service<DQMStore>().operator->()),
      conf_(iConfig),
      mfToken_(esConsumes()),
      tkGeomToken_(esConsumes()),
      dtGeomToken_(esConsumes()),
      cscGeomToken_(esConsumes()),
      rpcGeomToken_(esConsumes()) {
  splitTracks_ = conf_.getParameter<edm::InputTag>("splitTrackCollection");
  splitMuons_ = conf_.getParameter<edm::InputTag>("splitMuonCollection");
  splitTracksToken_ = consumes<std::vector<reco::Track> >(splitTracks_);
  splitMuonsToken_ = mayConsume<std::vector<reco::Muon> >(splitMuons_);
 
  plotMuons_ = conf_.getParameter<bool>("ifPlotMuons");
 
  // cuts
  pixelHitsPerLeg_ = conf_.getParameter<int>("pixelHitsPerLeg");
  totalHitsPerLeg_ = conf_.getParameter<int>("totalHitsPerLeg");
  d0Cut_ = conf_.getParameter<double>("d0Cut");
  dzCut_ = conf_.getParameter<double>("dzCut");
  ptCut_ = conf_.getParameter<double>("ptCut");
  norchiCut_ = conf_.getParameter<double>("norchiCut");
}

References conf_, d0Cut_, dzCut_, edm::ParameterSet::getParameter(), norchiCut_, pixelHitsPerLeg_, plotMuons_, ptCut_, splitMuons_, splitMuonsToken_, splitTracks_, splitTracksToken_, and totalHitsPerLeg_.

~TrackSplittingMonitor()

TrackSplittingMonitor::~TrackSplittingMonitor ( )

override

Definition at line 50 of file TrackSplittingMonitor.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 178 of file TrackSplittingMonitor.cc.

                                                                                       {
  theMagField = &iSetup.getData(mfToken_);
  theGeometry = &iSetup.getData(tkGeomToken_);
  dtGeometry = &iSetup.getData(dtGeomToken_);
  cscGeometry = &iSetup.getData(cscGeomToken_);
  rpcGeometry = &iSetup.getData(rpcGeomToken_);
 
  edm::Handle<std::vector<reco::Track> > splitTracks = iEvent.getHandle(splitTracksToken_);
  if (!splitTracks.isValid())
    return;
 
  edm::Handle<std::vector<reco::Muon> > splitMuons;
  if (plotMuons_) {
    splitMuons = iEvent.getHandle(splitMuonsToken_);
  }
 
  if (splitTracks->size() == 2) {
    // check that there are 2 tracks in split track collection
    edm::LogInfo("TrackSplittingMonitor") << "Split Track size: " << splitTracks->size();
 
    // split tracks calculations
    reco::Track track1 = splitTracks->at(0);
    reco::Track track2 = splitTracks->at(1);
 
    // -------------------------- basic selection ---------------------------
 
    // hit counting
    // looping through the hits for track 1
    double nRechits1 = 0;
    double nRechitinBPIX1 = 0;
    for (trackingRecHit_iterator iHit = track1.recHitsBegin(); iHit != track1.recHitsEnd(); ++iHit) {
      if ((*iHit)->isValid()) {
        nRechits1++;
        int type = (*iHit)->geographicalId().subdetId();
        if (type == int(PixelSubdetector::PixelBarrel)) {
          ++nRechitinBPIX1;
        }
      }
    }
    // looping through the hits for track 2
    double nRechits2 = 0;
    double nRechitinBPIX2 = 0;
    for (trackingRecHit_iterator iHit = track2.recHitsBegin(); iHit != track2.recHitsEnd(); ++iHit) {
      if ((*iHit)->isValid()) {
        nRechits2++;
        int type = (*iHit)->geographicalId().subdetId();
        if (type == int(PixelSubdetector::PixelBarrel)) {
          ++nRechitinBPIX2;
        }
      }
    }
 
    // DCA of each track
    double d01 = track1.d0();
    double dz1 = track1.dz();
    double d02 = track2.d0();
    double dz2 = track2.dz();
 
    // pT of each track
    double pt1 = track1.pt();
    double pt2 = track2.pt();
 
    // chi2 of each track
    double norchi1 = track1.normalizedChi2();
    double norchi2 = track2.normalizedChi2();
 
    // basic selection
    // pixel hits and total hits
    if ((nRechitinBPIX1 >= pixelHitsPerLeg_) && (nRechitinBPIX1 >= pixelHitsPerLeg_) &&
        (nRechits1 >= totalHitsPerLeg_) && (nRechits2 >= totalHitsPerLeg_)) {
      // dca cut
      if (((fabs(d01) < d0Cut_)) && (fabs(d02) < d0Cut_) && (fabs(dz2) < dzCut_) && (fabs(dz2) < dzCut_)) {
        // pt cut
        if ((pt1 + pt2) / 2 < ptCut_) {
          // chi2 cut
          if ((norchi1 < norchiCut_) && (norchi2 < norchiCut_)) {
            // passed all cuts...
            edm::LogInfo("TrackSplittingMonitor") << " Setected after all cuts ?";
 
            double ddxyVal = d01 - d02;
            double ddzVal = dz1 - dz2;
            double dphiVal = track1.phi() - track2.phi();
            double dthetaVal = track1.theta() - track2.theta();
            double dptVal = pt1 - pt2;
            double dcurvVal = (1 / pt1) - (1 / pt2);
 
            double d01ErrVal = track1.d0Error();
            double d02ErrVal = track2.d0Error();
            double dz1ErrVal = track1.dzError();
            double dz2ErrVal = track2.dzError();
            double phi1ErrVal = track1.phiError();
            double phi2ErrVal = track2.phiError();
            double theta1ErrVal = track1.thetaError();
            double theta2ErrVal = track2.thetaError();
            double pt1ErrVal = track1.ptError();
            double pt2ErrVal = track2.ptError();
 
            ddxyAbsoluteResiduals_tracker_->Fill(10000.0 * ddxyVal / sqrt(2.0));
            ddxyAbsoluteResiduals_tracker_->Fill(10000.0 * ddzVal / sqrt(2.0));
            ddxyAbsoluteResiduals_tracker_->Fill(1000.0 * dphiVal / sqrt(2.0));
            ddxyAbsoluteResiduals_tracker_->Fill(1000.0 * dthetaVal / sqrt(2.0));
            ddxyAbsoluteResiduals_tracker_->Fill(dptVal / sqrt(2.0));
            ddxyAbsoluteResiduals_tracker_->Fill(dcurvVal / sqrt(2.0));
 
            ddxyNormalizedResiduals_tracker_->Fill(ddxyVal / sqrt(d01ErrVal * d01ErrVal + d02ErrVal * d02ErrVal));
            ddxyNormalizedResiduals_tracker_->Fill(ddzVal / sqrt(dz1ErrVal * dz1ErrVal + dz2ErrVal * dz2ErrVal));
            ddxyNormalizedResiduals_tracker_->Fill(dphiVal / sqrt(phi1ErrVal * phi1ErrVal + phi2ErrVal * phi2ErrVal));
            ddxyNormalizedResiduals_tracker_->Fill(dthetaVal /
                                                   sqrt(theta1ErrVal * theta1ErrVal + theta2ErrVal * theta2ErrVal));
            ddxyNormalizedResiduals_tracker_->Fill(dptVal / sqrt(pt1ErrVal * pt1ErrVal + pt2ErrVal * pt2ErrVal));
            ddxyNormalizedResiduals_tracker_->Fill(
                dcurvVal / sqrt(pow(pt1ErrVal, 2) / pow(pt1, 4) + pow(pt2ErrVal, 2) / pow(pt2, 4)));
 
            // if do the same for split muons
            if (plotMuons_ && splitMuons.isValid()) {
              int gmCtr = 0;
              bool topGlobalMuonFlag = false;
              bool bottomGlobalMuonFlag = false;
              int topGlobalMuon = -1;
              int bottomGlobalMuon = -1;
              double topGlobalMuonNorchi2 = 1e10;
              double bottomGlobalMuonNorchi2 = 1e10;
 
              // check if usable split global muons
              for (std::vector<reco::Muon>::const_iterator gmI = splitMuons->begin(); gmI != splitMuons->end(); gmI++) {
                if (gmI->isTrackerMuon() && gmI->isStandAloneMuon() && gmI->isGlobalMuon()) {
                  reco::TrackRef trackerTrackRef1(splitTracks, 0);
                  reco::TrackRef trackerTrackRef2(splitTracks, 1);
 
                  if (gmI->innerTrack() == trackerTrackRef1) {
                    if (gmI->globalTrack()->normalizedChi2() < topGlobalMuonNorchi2) {
                      topGlobalMuonFlag = true;
                      topGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
                      topGlobalMuon = gmCtr;
                    }
                  }
                  if (gmI->innerTrack() == trackerTrackRef2) {
                    if (gmI->globalTrack()->normalizedChi2() < bottomGlobalMuonNorchi2) {
                      bottomGlobalMuonFlag = true;
                      bottomGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
                      bottomGlobalMuon = gmCtr;
                    }
                  }
                }
                gmCtr++;
              }
 
              if (bottomGlobalMuonFlag && topGlobalMuonFlag) {
                reco::Muon muonTop = splitMuons->at(topGlobalMuon);
                reco::Muon muonBottom = splitMuons->at(bottomGlobalMuon);
 
                reco::TrackRef glb1 = muonTop.globalTrack();
                reco::TrackRef glb2 = muonBottom.globalTrack();
 
                double ddxyValGlb = glb1->d0() - glb2->d0();
                double ddzValGlb = glb1->dz() - glb2->dz();
                double dphiValGlb = glb1->phi() - glb2->phi();
                double dthetaValGlb = glb1->theta() - glb2->theta();
                double dptValGlb = glb1->pt() - glb2->pt();
                double dcurvValGlb = (1 / glb1->pt()) - (1 / glb2->pt());
 
                double d01ErrValGlb = glb1->d0Error();
                double d02ErrValGlb = glb2->d0Error();
                double dz1ErrValGlb = glb1->dzError();
                double dz2ErrValGlb = glb2->dzError();
                double phi1ErrValGlb = glb1->phiError();
                double phi2ErrValGlb = glb2->phiError();
                double theta1ErrValGlb = glb1->thetaError();
                double theta2ErrValGlb = glb2->thetaError();
                double pt1ErrValGlb = glb1->ptError();
                double pt2ErrValGlb = glb2->ptError();
 
                ddxyAbsoluteResiduals_global_->Fill(10000.0 * ddxyValGlb / sqrt(2.0));
                ddxyAbsoluteResiduals_global_->Fill(10000.0 * ddzValGlb / sqrt(2.0));
                ddxyAbsoluteResiduals_global_->Fill(1000.0 * dphiValGlb / sqrt(2.0));
                ddxyAbsoluteResiduals_global_->Fill(1000.0 * dthetaValGlb / sqrt(2.0));
                ddxyAbsoluteResiduals_global_->Fill(dptValGlb / sqrt(2.0));
                ddxyAbsoluteResiduals_global_->Fill(dcurvValGlb / sqrt(2.0));
 
                ddxyNormalizedResiduals_global_->Fill(ddxyValGlb /
                                                      sqrt(d01ErrValGlb * d01ErrValGlb + d02ErrValGlb * d02ErrValGlb));
                ddxyNormalizedResiduals_global_->Fill(ddzValGlb /
                                                      sqrt(dz1ErrValGlb * dz1ErrValGlb + dz2ErrValGlb * dz2ErrValGlb));
                ddxyNormalizedResiduals_global_->Fill(
                    dphiValGlb / sqrt(phi1ErrValGlb * phi1ErrValGlb + phi2ErrValGlb * phi2ErrValGlb));
                ddxyNormalizedResiduals_global_->Fill(
                    dthetaValGlb / sqrt(theta1ErrValGlb * theta1ErrValGlb + theta2ErrValGlb * theta2ErrValGlb));
                ddxyNormalizedResiduals_global_->Fill(dptValGlb /
                                                      sqrt(pt1ErrValGlb * pt1ErrValGlb + pt2ErrValGlb * pt2ErrValGlb));
                ddxyNormalizedResiduals_global_->Fill(
                    dcurvValGlb / sqrt(pow(pt1ErrValGlb, 2) / pow(pt1, 4) + pow(pt2ErrValGlb, 2) / pow(pt2, 4)));
              }
 
            }  // end of split muons loop
          }
        }
      }
    }
  }
}

References cscGeometry, cscGeomToken_, reco::TrackBase::d0(), d0Cut_, reco::TrackBase::d0Error(), ddxyAbsoluteResiduals_global_, ddxyAbsoluteResiduals_tracker_, ddxyNormalizedResiduals_global_, ddxyNormalizedResiduals_tracker_, dtGeometry, dtGeomToken_, reco::TrackBase::dz(), dzCut_, reco::TrackBase::dzError(), dqm::impl::MonitorElement::Fill(), edm::EventSetup::getData(), reco::Muon::globalTrack(), iEvent, edm::HandleBase::isValid(), mfToken_, norchiCut_, reco::TrackBase::normalizedChi2(), reco::TrackBase::phi(), reco::TrackBase::phiError(), PixelSubdetector::PixelBarrel, pixelHitsPerLeg_, plotMuons_, funct::pow(), reco::TrackBase::pt(), HLT_FULL_cff::pt1, HLT_FULL_cff::pt2, ptCut_, reco::TrackBase::ptError(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), rpcGeometry, rpcGeomToken_, RecoMuonCosmics_cff::splitMuons, splitMuonsToken_, splitTracksToken_, mathSSE::sqrt(), theGeometry, theMagField, reco::TrackBase::theta(), reco::TrackBase::thetaError(), tkGeomToken_, and totalHitsPerLeg_.

bookHistograms()

void TrackSplittingMonitor::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 52 of file TrackSplittingMonitor.cc.

{
  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
  ibooker.setCurrentFolder(MEFolderName);
 
  // bin declarations
  int ddxyBin = conf_.getParameter<int>("ddxyBin");
  double ddxyMin = conf_.getParameter<double>("ddxyMin");
  double ddxyMax = conf_.getParameter<double>("ddxyMax");
 
  int ddzBin = conf_.getParameter<int>("ddzBin");
  double ddzMin = conf_.getParameter<double>("ddzMin");
  double ddzMax = conf_.getParameter<double>("ddzMax");
 
  int dphiBin = conf_.getParameter<int>("dphiBin");
  double dphiMin = conf_.getParameter<double>("dphiMin");
  double dphiMax = conf_.getParameter<double>("dphiMax");
 
  int dthetaBin = conf_.getParameter<int>("dthetaBin");
  double dthetaMin = conf_.getParameter<double>("dthetaMin");
  double dthetaMax = conf_.getParameter<double>("dthetaMax");
 
  int dptBin = conf_.getParameter<int>("dptBin");
  double dptMin = conf_.getParameter<double>("dptMin");
  double dptMax = conf_.getParameter<double>("dptMax");
 
  int dcurvBin = conf_.getParameter<int>("dcurvBin");
  double dcurvMin = conf_.getParameter<double>("dcurvMin");
  double dcurvMax = conf_.getParameter<double>("dcurvMax");
 
  int normBin = conf_.getParameter<int>("normBin");
  double normMin = conf_.getParameter<double>("normMin");
  double normMax = conf_.getParameter<double>("normMax");
 
  // declare histogram
  ddxyAbsoluteResiduals_tracker_ =
      ibooker.book1D("ddxyAbsoluteResiduals_tracker", "ddxyAbsoluteResiduals_tracker", ddxyBin, ddxyMin, ddxyMax);
  ddzAbsoluteResiduals_tracker_ =
      ibooker.book1D("ddzAbsoluteResiduals_tracker", "ddzAbsoluteResiduals_tracker", ddzBin, ddzMin, ddzMax);
  dphiAbsoluteResiduals_tracker_ =
      ibooker.book1D("dphiAbsoluteResiduals_tracker", "dphiAbsoluteResiduals_tracker", dphiBin, dphiMin, dphiMax);
  dthetaAbsoluteResiduals_tracker_ = ibooker.book1D(
      "dthetaAbsoluteResiduals_tracker", "dthetaAbsoluteResiduals_tracker", dthetaBin, dthetaMin, dthetaMax);
  dptAbsoluteResiduals_tracker_ =
      ibooker.book1D("dptAbsoluteResiduals_tracker", "dptAbsoluteResiduals_tracker", dptBin, dptMin, dptMax);
  dcurvAbsoluteResiduals_tracker_ =
      ibooker.book1D("dcurvAbsoluteResiduals_tracker", "dcurvAbsoluteResiduals_tracker", dcurvBin, dcurvMin, dcurvMax);
 
  ddxyNormalizedResiduals_tracker_ =
      ibooker.book1D("ddxyNormalizedResiduals_tracker", "ddxyNormalizedResiduals_tracker", normBin, normMin, normMax);
  ddzNormalizedResiduals_tracker_ =
      ibooker.book1D("ddzNormalizedResiduals_tracker", "ddzNormalizedResiduals_tracker", normBin, normMin, normMax);
  dphiNormalizedResiduals_tracker_ =
      ibooker.book1D("dphiNormalizedResiduals_tracker", "dphiNormalizedResiduals_tracker", normBin, normMin, normMax);
  dthetaNormalizedResiduals_tracker_ = ibooker.book1D(
      "dthetaNormalizedResiduals_tracker", "dthetaNormalizedResiduals_tracker", normBin, normMin, normMax);
  dptNormalizedResiduals_tracker_ =
      ibooker.book1D("dptNormalizedResiduals_tracker", "dptNormalizedResiduals_tracker", normBin, normMin, normMax);
  dcurvNormalizedResiduals_tracker_ =
      ibooker.book1D("dcurvNormalizedResiduals_tracker", "dcurvNormalizedResiduals_tracker", normBin, normMin, normMax);
 
  if (plotMuons_) {
    ddxyAbsoluteResiduals_global_ =
        ibooker.book1D("ddxyAbsoluteResiduals_global", "ddxyAbsoluteResiduals_global", ddxyBin, ddxyMin, ddxyMax);
    ddzAbsoluteResiduals_global_ =
        ibooker.book1D("ddzAbsoluteResiduals_global", "ddzAbsoluteResiduals_global", ddzBin, ddzMin, ddzMax);
    dphiAbsoluteResiduals_global_ =
        ibooker.book1D("dphiAbsoluteResiduals_global", "dphiAbsoluteResiduals_global", dphiBin, dphiMin, dphiMax);
    dthetaAbsoluteResiduals_global_ = ibooker.book1D(
        "dthetaAbsoluteResiduals_global", "dthetaAbsoluteResiduals_global", dthetaBin, dthetaMin, dthetaMax);
    dptAbsoluteResiduals_global_ =
        ibooker.book1D("dptAbsoluteResiduals_global", "dptAbsoluteResiduals_global", dptBin, dptMin, dptMax);
    dcurvAbsoluteResiduals_global_ =
        ibooker.book1D("dcurvAbsoluteResiduals_global", "dcurvAbsoluteResiduals_global", dcurvBin, dcurvMin, dcurvMax);
 
    ddxyNormalizedResiduals_global_ =
        ibooker.book1D("ddxyNormalizedResiduals_global", "ddxyNormalizedResiduals_global", normBin, normMin, normMax);
    ddzNormalizedResiduals_global_ =
        ibooker.book1D("ddzNormalizedResiduals_global", "ddzNormalizedResiduals_global", normBin, normMin, normMax);
    dphiNormalizedResiduals_global_ =
        ibooker.book1D("dphiNormalizedResiduals_global", "dphiNormalizedResiduals_global", normBin, normMin, normMax);
    dthetaNormalizedResiduals_global_ = ibooker.book1D(
        "dthetaNormalizedResiduals_global", "dthetaNormalizedResiduals_global", normBin, normMin, normMax);
    dptNormalizedResiduals_global_ =
        ibooker.book1D("dptNormalizedResiduals_global", "dptNormalizedResiduals_global", normBin, normMin, normMax);
    dcurvNormalizedResiduals_global_ =
        ibooker.book1D("dcurvNormalizedResiduals_global", "dcurvNormalizedResiduals_global", normBin, normMin, normMax);
  }
 
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta d_{xy})/#sqrt{2} [#mum]");
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta d_{z})/#sqrt{2} [#mum]");
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta #phi)/#sqrt{2} [mrad]");
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta #theta)/#sqrt{2} [mrad]");
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta pT)/#sqrt{2} [GeV]");
  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta (1/pT))/#sqrt{2} [GeV^{-1}]");
 
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta d_{xy}/#sigma(d_{xy}");
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta d_{z}/#sigma(d_{z})");
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta #phi/#sigma(d_{#phi})");
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta #theta/#sigma(d_{#theta})");
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta p_{T}/#sigma(p_{T})");
  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta 1/p_{T}/#sigma(1/p_{T})");
 
  if (plotMuons_) {
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta d_{xy})/#sqrt{2} [#mum]");
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta d_{z})/#sqrt{2} [#mum]");
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta #phi)/#sqrt{2} [mrad]");
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta #theta)/#sqrt{2} [mrad]");
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta pT)/#sqrt{2} [GeV]");
    ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta (1/pT))/#sqrt{2} [GeV^{-1}]");
 
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta d_{xy}/#sigma(d_{xy}");
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta d_{z}/#sigma(d_{z})");
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta #phi/#sigma(d_{#phi})");
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta #theta/#sigma(d_{#theta})");
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta p_{T}/#sigma(p_{T})");
    ddxyNormalizedResiduals_global_->setAxisTitle("#delta 1/p_{T}/#sigma(1/p_{T})");
  }
}

References dqm::implementation::IBooker::book1D(), conf_, dcurvAbsoluteResiduals_global_, dcurvAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::dcurvBin, TrackSplittingMonitor_cfi::dcurvMax, TrackSplittingMonitor_cfi::dcurvMin, dcurvNormalizedResiduals_global_, dcurvNormalizedResiduals_tracker_, ddxyAbsoluteResiduals_global_, ddxyAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::ddxyBin, TrackSplittingMonitor_cfi::ddxyMax, TrackSplittingMonitor_cfi::ddxyMin, ddxyNormalizedResiduals_global_, ddxyNormalizedResiduals_tracker_, ddzAbsoluteResiduals_global_, ddzAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::ddzBin, TrackSplittingMonitor_cfi::ddzMax, TrackSplittingMonitor_cfi::ddzMin, ddzNormalizedResiduals_global_, ddzNormalizedResiduals_tracker_, dphiAbsoluteResiduals_global_, dphiAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::dphiBin, TrackSplittingMonitor_cfi::dphiMax, TrackSplittingMonitor_cfi::dphiMin, dphiNormalizedResiduals_global_, dphiNormalizedResiduals_tracker_, dptAbsoluteResiduals_global_, dptAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::dptBin, TrackSplittingMonitor_cfi::dptMax, TrackSplittingMonitor_cfi::dptMin, dptNormalizedResiduals_global_, dptNormalizedResiduals_tracker_, dthetaAbsoluteResiduals_global_, dthetaAbsoluteResiduals_tracker_, TrackSplittingMonitor_cfi::dthetaBin, TrackSplittingMonitor_cfi::dthetaMax, TrackSplittingMonitor_cfi::dthetaMin, dthetaNormalizedResiduals_global_, dthetaNormalizedResiduals_tracker_, edm::ParameterSet::getParameter(), TrackSplittingMonitor_cfi::normBin, TrackSplittingMonitor_cfi::normMax, TrackSplittingMonitor_cfi::normMin, plotMuons_, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

doProfileX() [1/2]

void TrackSplittingMonitor::doProfileX ( MonitorElement *  th2m, MonitorElement *  me  )

private

doProfileX() [2/2]

void TrackSplittingMonitor::doProfileX ( TH2 *  th2, MonitorElement *  me  )

private

Member Data Documentation

conf_

edm::ParameterSet TrackSplittingMonitor::conf_

private

Definition at line 60 of file TrackSplittingMonitor.h.

Referenced by bookHistograms(), and TrackSplittingMonitor().

cscGeometry

const CSCGeometry* TrackSplittingMonitor::cscGeometry

private

Definition at line 71 of file TrackSplittingMonitor.h.

Referenced by analyze().

cscGeomToken_

const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> TrackSplittingMonitor::cscGeomToken_

private

Definition at line 65 of file TrackSplittingMonitor.h.

Referenced by analyze().

d0Cut_

double TrackSplittingMonitor::d0Cut_

private

Definition at line 82 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

dcurvAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::dcurvAbsoluteResiduals_global_

private

Definition at line 107 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dcurvAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dcurvAbsoluteResiduals_tracker_

private

Definition at line 93 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dcurvNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::dcurvNormalizedResiduals_global_

private

Definition at line 114 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dcurvNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dcurvNormalizedResiduals_tracker_

private

Definition at line 100 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

ddxyAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::ddxyAbsoluteResiduals_global_

private

Definition at line 102 of file TrackSplittingMonitor.h.

Referenced by analyze(), and bookHistograms().

ddxyAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::ddxyAbsoluteResiduals_tracker_

private

Definition at line 88 of file TrackSplittingMonitor.h.

Referenced by analyze(), and bookHistograms().

ddxyNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::ddxyNormalizedResiduals_global_

private

Definition at line 109 of file TrackSplittingMonitor.h.

Referenced by analyze(), and bookHistograms().

ddxyNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::ddxyNormalizedResiduals_tracker_

private

Definition at line 95 of file TrackSplittingMonitor.h.

Referenced by analyze(), and bookHistograms().

ddzAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::ddzAbsoluteResiduals_global_

private

Definition at line 103 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

ddzAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::ddzAbsoluteResiduals_tracker_

private

Definition at line 89 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

ddzNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::ddzNormalizedResiduals_global_

private

Definition at line 110 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

ddzNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::ddzNormalizedResiduals_tracker_

private

Definition at line 96 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dphiAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::dphiAbsoluteResiduals_global_

private

Definition at line 104 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dphiAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dphiAbsoluteResiduals_tracker_

private

Definition at line 90 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dphiNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::dphiNormalizedResiduals_global_

private

Definition at line 111 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dphiNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dphiNormalizedResiduals_tracker_

private

Definition at line 97 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dptAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::dptAbsoluteResiduals_global_

private

Definition at line 106 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dptAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dptAbsoluteResiduals_tracker_

private

Definition at line 92 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dptNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::dptNormalizedResiduals_global_

private

Definition at line 113 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dptNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dptNormalizedResiduals_tracker_

private

Definition at line 99 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dqmStore_

DQMStore* TrackSplittingMonitor::dqmStore_

private

Definition at line 59 of file TrackSplittingMonitor.h.

dtGeometry

const DTGeometry* TrackSplittingMonitor::dtGeometry

private

Definition at line 70 of file TrackSplittingMonitor.h.

Referenced by analyze().

dtGeomToken_

const edm::ESGetToken<DTGeometry, MuonGeometryRecord> TrackSplittingMonitor::dtGeomToken_

private

Definition at line 64 of file TrackSplittingMonitor.h.

Referenced by analyze().

dthetaAbsoluteResiduals_global_

MonitorElement* TrackSplittingMonitor::dthetaAbsoluteResiduals_global_

private

Definition at line 105 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dthetaAbsoluteResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dthetaAbsoluteResiduals_tracker_

private

Definition at line 91 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dthetaNormalizedResiduals_global_

MonitorElement* TrackSplittingMonitor::dthetaNormalizedResiduals_global_

private

Definition at line 112 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dthetaNormalizedResiduals_tracker_

MonitorElement* TrackSplittingMonitor::dthetaNormalizedResiduals_tracker_

private

Definition at line 98 of file TrackSplittingMonitor.h.

Referenced by bookHistograms().

dzCut_

double TrackSplittingMonitor::dzCut_

private

Definition at line 83 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

histname

std::string TrackSplittingMonitor::histname

private

Definition at line 57 of file TrackSplittingMonitor.h.

mfToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TrackSplittingMonitor::mfToken_

private

Definition at line 62 of file TrackSplittingMonitor.h.

Referenced by analyze().

norchiCut_

double TrackSplittingMonitor::norchiCut_

private

Definition at line 85 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

pixelHitsPerLeg_

int TrackSplittingMonitor::pixelHitsPerLeg_

private

Definition at line 80 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

plotMuons_

bool TrackSplittingMonitor::plotMuons_

private

Definition at line 79 of file TrackSplittingMonitor.h.

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

ptCut_

double TrackSplittingMonitor::ptCut_

private

Definition at line 84 of file TrackSplittingMonitor.h.

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

rpcGeometry

const RPCGeometry* TrackSplittingMonitor::rpcGeometry

private

Definition at line 72 of file TrackSplittingMonitor.h.

Referenced by analyze().

rpcGeomToken_

const edm::ESGetToken<RPCGeometry, MuonGeometryRecord> TrackSplittingMonitor::rpcGeomToken_

private

Definition at line 66 of file TrackSplittingMonitor.h.

Referenced by analyze().

splitMuons_

edm::InputTag TrackSplittingMonitor::splitMuons_

private

Definition at line 75 of file TrackSplittingMonitor.h.

Referenced by TrackSplittingMonitor().

splitMuonsToken_

edm::EDGetTokenT<std::vector<reco::Muon> > TrackSplittingMonitor::splitMuonsToken_

private

Definition at line 77 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

splitTracks_

edm::InputTag TrackSplittingMonitor::splitTracks_

private

Definition at line 74 of file TrackSplittingMonitor.h.

Referenced by TrackSplittingMonitor().

splitTracksToken_

edm::EDGetTokenT<std::vector<reco::Track> > TrackSplittingMonitor::splitTracksToken_

private

Definition at line 76 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().

theGeometry

const TrackerGeometry* TrackSplittingMonitor::theGeometry

private

Definition at line 68 of file TrackSplittingMonitor.h.

Referenced by analyze().

theMagField

const MagneticField* TrackSplittingMonitor::theMagField

private

Definition at line 69 of file TrackSplittingMonitor.h.

Referenced by analyze().

tkGeomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrackSplittingMonitor::tkGeomToken_

private

Definition at line 63 of file TrackSplittingMonitor.h.

Referenced by analyze().

totalHitsPerLeg_

int TrackSplittingMonitor::totalHitsPerLeg_

private

Definition at line 81 of file TrackSplittingMonitor.h.

Referenced by analyze(), and TrackSplittingMonitor().