SiPixelTrackResidualSource Class Reference

#include <SiPixelTrackResidualSource.h>

Inheritance diagram for SiPixelTrackResidualSource:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &r, edm::EventSetup const &iSetup) override
void	getrococcupancy (DetId detId, const edm::DetSetVector< PixelDigi > &diginp, const TrackerTopology *const tTopo, std::vector< MonitorElement *> meinput)
	SiPixelTrackResidualSource (const edm::ParameterSet &)
void	triplets (double x1, double y1, double z1, double x2, double y2, double z2, double x3, double y3, double z3, double ptsig, double &dc, double &dz, double kap)
	~SiPixelTrackResidualSource () 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
	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

Public Attributes
std::string	topFolderName_

Private Member Functions
void	getepixrococcupancyofftrk (DetId detId, const TrackerTopology *const tTopo, float xclust, float yclust, float z, MonitorElement *meinput)
void	getepixrococcupancyontrk (const TrackerTopology *const tTopo, TransientTrackingRecHit::ConstRecHitPointer hit, float xclust, float yclust, float z, MonitorElement *meinput)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
bool	bladeOn
edm::InputTag	clustersrc_
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	clustersrcToken_
bool	debug_
edm::InputTag	digisrc_
edm::EDGetTokenT< edm::DetSetVector< PixelDigi > >	digisrcToken_
bool	diskOn
bool	firstRun
edm::EDGetTokenT< reco::TrackCollection >	generalTracksToken_
bool	isUpgrade
bool	ladOn
bool	layOn
MonitorElement *	meClChargeNotOnTrack_all
MonitorElement *	meClChargeNotOnTrack_bpix
std::vector< MonitorElement * >	meClChargeNotOnTrack_diskms
std::vector< MonitorElement * >	meClChargeNotOnTrack_diskps
MonitorElement *	meClChargeNotOnTrack_fpix
std::vector< MonitorElement * >	meClChargeNotOnTrack_layers
MonitorElement *	meClChargeOnTrack_all
MonitorElement *	meClChargeOnTrack_bpix
std::vector< MonitorElement * >	meClChargeOnTrack_diskms
std::vector< MonitorElement * >	meClChargeOnTrack_diskps
MonitorElement *	meClChargeOnTrack_fpix
std::vector< MonitorElement * >	meClChargeOnTrack_layers
std::vector< MonitorElement * >	meClPosDisksmzNotOnTrack
std::vector< MonitorElement * >	meClPosDisksmzOnTrack
std::vector< MonitorElement * >	meClPosDiskspzNotOnTrack
std::vector< MonitorElement * >	meClPosDiskspzOnTrack
std::vector< MonitorElement * >	meClPosLayersLadVsModOnTrack
std::vector< MonitorElement * >	meClPosLayersNotOnTrack
std::vector< MonitorElement * >	meClPosLayersOnTrack
MonitorElement *	meClSizeNotOnTrack_all
MonitorElement *	meClSizeNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeNotOnTrack_diskps
MonitorElement *	meClSizeNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeNotOnTrack_layers
MonitorElement *	meClSizeOnTrack_all
MonitorElement *	meClSizeOnTrack_bpix
std::vector< MonitorElement * >	meClSizeOnTrack_diskms
std::vector< MonitorElement * >	meClSizeOnTrack_diskps
MonitorElement *	meClSizeOnTrack_fpix
std::vector< MonitorElement * >	meClSizeOnTrack_layers
MonitorElement *	meClSizeXNotOnTrack_all
MonitorElement *	meClSizeXNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeXNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeXNotOnTrack_diskps
MonitorElement *	meClSizeXNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeXNotOnTrack_layers
MonitorElement *	meClSizeXOnTrack_all
MonitorElement *	meClSizeXOnTrack_bpix
std::vector< MonitorElement * >	meClSizeXOnTrack_diskms
std::vector< MonitorElement * >	meClSizeXOnTrack_diskps
MonitorElement *	meClSizeXOnTrack_fpix
std::vector< MonitorElement * >	meClSizeXOnTrack_layers
MonitorElement *	meClSizeYNotOnTrack_all
MonitorElement *	meClSizeYNotOnTrack_bpix
std::vector< MonitorElement * >	meClSizeYNotOnTrack_diskms
std::vector< MonitorElement * >	meClSizeYNotOnTrack_diskps
MonitorElement *	meClSizeYNotOnTrack_fpix
std::vector< MonitorElement * >	meClSizeYNotOnTrack_layers
MonitorElement *	meClSizeYOnTrack_all
MonitorElement *	meClSizeYOnTrack_bpix
std::vector< MonitorElement * >	meClSizeYOnTrack_diskms
std::vector< MonitorElement * >	meClSizeYOnTrack_diskps
MonitorElement *	meClSizeYOnTrack_fpix
std::vector< MonitorElement * >	meClSizeYOnTrack_layers
MonitorElement *	meHitProbability
MonitorElement *	meNClustersNotOnTrack_all
MonitorElement *	meNClustersNotOnTrack_bpix
std::vector< MonitorElement * >	meNClustersNotOnTrack_diskms
std::vector< MonitorElement * >	meNClustersNotOnTrack_diskps
MonitorElement *	meNClustersNotOnTrack_fpix
std::vector< MonitorElement * >	meNClustersNotOnTrack_layers
MonitorElement *	meNClustersOnTrack_all
MonitorElement *	meNClustersOnTrack_bpix
std::vector< MonitorElement * >	meNClustersOnTrack_diskms
std::vector< MonitorElement * >	meNClustersOnTrack_diskps
MonitorElement *	meNClustersOnTrack_fpix
std::vector< MonitorElement * >	meNClustersOnTrack_layers
MonitorElement *	meNofClustersNotOnTrack_
MonitorElement *	meNofClustersOnTrack_
std::vector< MonitorElement * >	meNofClustersvsPhiOnTrack_diskms
std::vector< MonitorElement * >	meNofClustersvsPhiOnTrack_diskps
std::vector< MonitorElement * >	meNofClustersvsPhiOnTrack_layers
MonitorElement *	meNofTracks_
MonitorElement *	meNofTracksInPixVol_
std::vector< MonitorElement * >	meResidualXSummedLay
std::vector< MonitorElement * >	meResidualYSummedLay
MonitorElement *	meRocBladevsDiskEndcapOffTrk
MonitorElement *	meRocBladevsDiskEndcapOnTrk
MonitorElement *	meSubdetResidualX [3]
MonitorElement *	meSubdetResidualY [3]
std::vector< MonitorElement * >	meZeroRocLadvsModOffTrackBarrel
std::vector< MonitorElement * >	meZeroRocLadvsModOnTrackBarrel
bool	modOn
int	NLowProb
int	noOfDisks
int	noOfLayers
int	NTotal
edm::EDGetTokenT< reco::VertexCollection >	offlinePrimaryVerticesToken_
bool	phiOn
edm::ParameterSet	pSet_
double	ptminres_
bool	reducedSet
bool	ringOn
edm::InputTag	src_
std::map< uint32_t, SiPixelTrackResidualModule * >	theSiPixelStructure
edm::EDGetTokenT< TrajTrackAssociationCollection >	trackAssociationToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomTokenBeginRun_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoTokenBeginRun_
edm::InputTag	tracksrc_
edm::EDGetTokenT< std::vector< Trajectory > >	tracksrcToken_
edm::EDGetTokenT< std::vector< reco::Track > >	trackToken_
edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	transientTrackBuilderToken_
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	transientTrackingRecHitBuilderToken_
std::string	ttrhbuilder_
std::string	vtxsrc_

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

Definition at line 46 of file SiPixelTrackResidualSource.h.

Constructor & Destructor Documentation

SiPixelTrackResidualSource()

SiPixelTrackResidualSource::SiPixelTrackResidualSource ( const edm::ParameterSet &  pSet )

explicit

Definition at line 48 of file SiPixelTrackResidualSource.cc.

References beamSpotToken_, bladeOn, clustersrc_, clustersrcToken_, debug_, digisrc_, digisrcToken_, diskOn, firstRun, generalTracksToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ladOn, layOn, modOn, NLowProb, NTotal, offlinePrimaryVerticesToken_, phiOn, pSet_, ptminres_, ringOn, src_, AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, trackAssociationToken_, trackerGeomToken_, trackerGeomTokenBeginRun_, trackerTopoToken_, trackerTopoTokenBeginRun_, tracksrc_, tracksrcToken_, trackToken_, transientTrackBuilderToken_, transientTrackingRecHitBuilderToken_, ttrhbuilder_, and vtxsrc_.

    : pSet_(pSet),
      modOn(pSet.getUntrackedParameter<bool>("modOn", true)),
      reducedSet(pSet.getUntrackedParameter<bool>("reducedSet", true)),
      ladOn(pSet.getUntrackedParameter<bool>("ladOn", false)),
      layOn(pSet.getUntrackedParameter<bool>("layOn", false)),
      phiOn(pSet.getUntrackedParameter<bool>("phiOn", false)),
      ringOn(pSet.getUntrackedParameter<bool>("ringOn", false)),
      bladeOn(pSet.getUntrackedParameter<bool>("bladeOn", false)),
      diskOn(pSet.getUntrackedParameter<bool>("diskOn", false)),
      isUpgrade(pSet.getUntrackedParameter<bool>("isUpgrade", false)),
      noOfLayers(0),
      noOfDisks(0) {
  pSet_ = pSet;
  debug_ = pSet_.getUntrackedParameter<bool>("debug", false);
  src_ = pSet_.getParameter<edm::InputTag>("src");
  clustersrc_ = pSet_.getParameter<edm::InputTag>("clustersrc");
  tracksrc_ = pSet_.getParameter<edm::InputTag>("trajectoryInput");
  ttrhbuilder_ = pSet_.getParameter<std::string>("TTRHBuilder");
  ptminres_ = pSet.getUntrackedParameter<double>("PtMinRes", 4.0);
  beamSpotToken_ = consumes<reco::BeamSpot>(std::string("offlineBeamSpot"));
  vtxsrc_ = pSet_.getUntrackedParameter<std::string>("vtxsrc", "offlinePrimaryVertices");
  offlinePrimaryVerticesToken_ =
      consumes<reco::VertexCollection>(vtxsrc_);  // consumes<reco::VertexCollection>(std::string("hiSelectedVertex"));
                                                  // //"offlinePrimaryVertices"));
  generalTracksToken_ = consumes<reco::TrackCollection>(pSet_.getParameter<edm::InputTag>("tracksrc"));
  tracksrcToken_ = consumes<std::vector<Trajectory>>(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
  trackToken_ = consumes<std::vector<reco::Track>>(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
  trackAssociationToken_ =
      consumes<TrajTrackAssociationCollection>(pSet_.getParameter<edm::InputTag>("trajectoryInput"));
  clustersrcToken_ = consumes<edmNew::DetSetVector<SiPixelCluster>>(pSet_.getParameter<edm::InputTag>("clustersrc"));
  digisrc_ = pSet_.getParameter<edm::InputTag>("digisrc");
  digisrcToken_ = consumes<edm::DetSetVector<PixelDigi>>(pSet_.getParameter<edm::InputTag>("digisrc"));

  trackerTopoToken_ = esConsumes<TrackerTopology, TrackerTopologyRcd>();
  trackerGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
  transientTrackBuilderToken_ =
      esConsumes<TransientTrackBuilder, TransientTrackRecord>(edm::ESInputTag("", "TransientTrackBuilder"));
  transientTrackingRecHitBuilderToken_ =
      esConsumes<TransientTrackingRecHitBuilder, TransientRecHitRecord>(edm::ESInputTag("", ttrhbuilder_));
  trackerTopoTokenBeginRun_ = esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>();
  trackerGeomTokenBeginRun_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>();

  LogInfo("PixelDQM") << "SiPixelTrackResidualSource constructor" << endl;
  LogInfo("PixelDQM") << "Mod/Lad/Lay/Phi " << modOn << "/" << ladOn << "/" << layOn << "/" << phiOn << std::endl;
  LogInfo("PixelDQM") << "Blade/Disk/Ring" << bladeOn << "/" << diskOn << "/" << ringOn << std::endl;

  topFolderName_ = pSet_.getParameter<std::string>("TopFolderName");

  firstRun = true;
  NTotal = 0;
  NLowProb = 0;
}

~SiPixelTrackResidualSource()

SiPixelTrackResidualSource::~SiPixelTrackResidualSource ( )

override

Definition at line 102 of file SiPixelTrackResidualSource.cc.

References theSiPixelStructure.

                                                        {
  LogInfo("PixelDQM") << "SiPixelTrackResidualSource destructor" << endl;

  std::map<uint32_t, SiPixelTrackResidualModule *>::iterator struct_iter;
  for (struct_iter = theSiPixelStructure.begin(); struct_iter != theSiPixelStructure.end(); struct_iter++) {
    delete struct_iter->second;
    struct_iter->second = nullptr;
  }
}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 823 of file SiPixelTrackResidualSource.cc.

References funct::abs(), Matriplex::atan2(), Reference_intrackfit_cff::barrel, beamSpotToken_, edmNew::DetSetVector< T >::begin(), bladeOn, TransientTrackingRecHitBuilder::build(), TransientTrackBuilder::build(), ALCARECOTkAlJpsiMuMu_cff::charge, SiPixelRecHit::clusterProbability(), clustersrcToken_, gather_cfg::cout, d0, debug_, hcalRecHitTable_cff::detId, TrackerGeometry::dets(), digisrcToken_, PixelEndcapName::diskName(), diskOn, PVValHelper::dz, edmNew::DetSetVector< T >::end(), makeMuonMisalignmentScenario::endcap, edm::HandleBase::failedToGet(), dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), generalTracksToken_, getepixrococcupancyofftrk(), getepixrococcupancyontrk(), edm::EventSetup::getHandle(), getrococcupancy(), corrVsCorr::gX, corrVsCorr::gY, SiPixelRecHit::hasFilledProb(), mps_fire::i, TrackerGeometry::idToDet(), iEvent, reco::TransientTrack::initialFreeState(), reco::TransientTrack::innermostMeasurementState(), createfilelist::int, edm::Ref< C, T, F >::isNonnull(), isUpgrade, sistrip::SpyUtilities::isValid(), TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, SummaryClient_cfi::labels, PVValHelper::ladder, PixelBarrelName::ladderName(), ladOn, PixelBarrelNameUpgrade::layerName(), PixelBarrelName::layerName(), layOn, TrajectoryStateOnSurface::localParameters(), Topology::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), match(), meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskms, meClChargeNotOnTrack_diskps, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layers, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskms, meClChargeOnTrack_diskps, meClChargeOnTrack_fpix, meClChargeOnTrack_layers, meClPosDisksmzNotOnTrack, meClPosDisksmzOnTrack, meClPosDiskspzNotOnTrack, meClPosDiskspzOnTrack, meClPosLayersLadVsModOnTrack, meClPosLayersNotOnTrack, meClPosLayersOnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskms, meClSizeNotOnTrack_diskps, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layers, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskms, meClSizeOnTrack_diskps, meClSizeOnTrack_fpix, meClSizeOnTrack_layers, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskms, meClSizeXNotOnTrack_diskps, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layers, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskms, meClSizeXOnTrack_diskps, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layers, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskms, meClSizeYNotOnTrack_diskps, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layers, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskms, meClSizeYOnTrack_diskps, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layers, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskms, meNClustersNotOnTrack_diskps, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layers, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskms, meNClustersOnTrack_diskps, meNClustersOnTrack_fpix, meNClustersOnTrack_layers, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofClustersvsPhiOnTrack_diskms, meNofClustersvsPhiOnTrack_diskps, meNofClustersvsPhiOnTrack_layers, meNofTracks_, meNofTracksInPixVol_, meResidualXSummedLay, meResidualYSummedLay, meRocBladevsDiskEndcapOffTrk, meRocBladevsDiskEndcapOnTrk, meZeroRocLadvsModOffTrackBarrel, meZeroRocLadvsModOnTrackBarrel, modOn, callgraph::module, PixelBarrelName::moduleName(), LocalTrajectoryParameters::momentum(), reco::Vertex::ndof(), NLowProb, noOfDisks, noOfLayers, NTotal, offlinePrimaryVerticesToken_, phi, PV3DBase< T, PVType, FrameType >::phi(), phiOn, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, funct::pow(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), DiDispStaMuonMonitor_cfi::pt, ptminres_, TrackerTopology::pxbLayer(), reducedSet, ringOn, PixelBarrelName::shell(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), DetId::subdetId(), GeomDet::surface(), funct::tan(), groupFilesInBlocks::temp, theSiPixelStructure, Surface::toGlobal(), trackAssociationToken_, DetId::Tracker, trackerGeomToken_, trackerTopoToken_, DiMuonV_cfg::tracks, tracksrcToken_, trackToken_, transientTrackBuilderToken_, transientTrackingRecHitBuilderToken_, FreeTrajectoryState::transverseCurvature(), triplets(), parallelization::uint, AlignmentTracksFromVertexSelector_cfi::vertices, x, PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x0(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, y, PV2DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y0(), testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, z, PV3DBase< T, PVType, FrameType >::z(), reco::BeamSpot::z0(), ppsModifySingularModes_cfi::z1, and ppsModifySingularModes_cfi::z2.

                                                                                            {
  // Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(trackerTopoToken_);
  const TrackerTopology *tTopo = tTopoHandle.product();

  // retrieve TrackerGeometry again and MagneticField for use in transforming
  // a TrackCandidate's P(ersistent)TrajectoryStateoOnDet (PTSoD) to a
  // TrajectoryStateOnSurface (TSoS)
  edm::ESHandle<TrackerGeometry> TG = iSetup.getHandle(trackerGeomToken_);
  const TrackerGeometry *theTrackerGeometry = TG.product();

  // analytic triplet method to calculate the track residuals in the pixe barrel
  // detector

  //--------------------------------------------------------------------
  // beam spot:
  //
  edm::Handle<reco::BeamSpot> rbs;
  // iEvent.getByLabel( "offlineBeamSpot", rbs );
  iEvent.getByToken(beamSpotToken_, rbs);
  math::XYZPoint bsP = math::XYZPoint(0, 0, 0);
  if (!rbs.failedToGet() && rbs.isValid()) {
    bsP = math::XYZPoint(rbs->x0(), rbs->y0(), rbs->z0());
  }

  //--------------------------------------------------------------------
  // primary vertices:
  //
  edm::Handle<reco::VertexCollection> vertices;
  // iEvent.getByLabel("offlinePrimaryVertices", vertices );
  iEvent.getByToken(offlinePrimaryVerticesToken_, vertices);

  if (vertices.failedToGet())
    return;
  if (!vertices.isValid())
    return;

  math::XYZPoint vtxN = math::XYZPoint(0, 0, 0);
  math::XYZPoint vtxP = math::XYZPoint(0, 0, 0);

  double bestNdof = 0.0;
  double maxSumPt = 0.0;
  reco::Vertex bestPvx;
  for (reco::VertexCollection::const_iterator iVertex = vertices->begin(); iVertex != vertices->end(); ++iVertex) {
    if (iVertex->ndof() > bestNdof) {
      bestNdof = iVertex->ndof();
      vtxN = math::XYZPoint(iVertex->x(), iVertex->y(), iVertex->z());
    }  // ndof
    if (iVertex->p4().pt() > maxSumPt) {
      maxSumPt = iVertex->p4().pt();
      vtxP = math::XYZPoint(iVertex->x(), iVertex->y(), iVertex->z());
      bestPvx = *iVertex;
    }  // sumpt

  }  // vertex

  if (maxSumPt < 1.0)
    vtxP = vtxN;

  //---------------------------------------------
  // get Tracks
  //
  edm::Handle<reco::TrackCollection> TracksForRes;
  // iEvent.getByLabel( "generalTracks", TracksForRes );
  iEvent.getByToken(generalTracksToken_, TracksForRes);

  if (debug_) {
    edm::EDConsumerBase::Labels labels;
    labelsForToken(generalTracksToken_, labels);
    std::cout << "Track for Res from " << labels.module << std::endl;
  }

  //
  // transient track builder, needs B-field from data base (global tag in .py)
  //
  edm::ESHandle<TransientTrackBuilder> theB = iSetup.getHandle(transientTrackBuilderToken_);

  // get the TransienTrackingRecHitBuilder needed for extracting the global
  // position of the hits in the pixel
  edm::ESHandle<TransientTrackingRecHitBuilder> theTrackerRecHitBuilder =
      iSetup.getHandle(transientTrackingRecHitBuilderToken_);

  // check that tracks are valid
  if (TracksForRes.failedToGet())
    return;
  if (!TracksForRes.isValid())
    return;

  //----------------------------------------------------------------------------
  // Residuals:
  //
  for (reco::TrackCollection::const_iterator iTrack = TracksForRes->begin(); iTrack != TracksForRes->end(); ++iTrack) {
    // count
    // Calculate minimal track pt before curling
    // cpt = cqRB = 0.3*R[m]*B[T] = 1.14*R[m] for B=3.8T
    // D = 2R = 2*pt/1.14
    // calo: D = 1.3 m => pt = 0.74 GeV/c
    double pt = iTrack->pt();
    if (pt < 0.75)
      continue;  // curls up
    if (abs(iTrack->dxy(vtxP)) > 5 * iTrack->dxyError())
      continue;  // not prompt

    double charge = iTrack->charge();

    reco::TransientTrack tTrack = theB->build(*iTrack);
    // get curvature of the track, needed for the residuals
    double kap = tTrack.initialFreeState().transverseCurvature();
    // needed for the TransienTrackingRecHitBuilder
    TrajectoryStateOnSurface initialTSOS = tTrack.innermostMeasurementState();
    if (iTrack->extra().isNonnull() && iTrack->extra().isAvailable()) {
      double x1 = 0;
      double y1 = 0;
      double z1 = 0;
      double x2 = 0;
      double y2 = 0;
      double z2 = 0;
      double x3 = 0;
      double y3 = 0;
      double z3 = 0;
      int n1 = 0;
      int n2 = 0;
      int n3 = 0;

      // for saving the pixel barrel hits
      vector<TransientTrackingRecHit::RecHitPointer> GoodPixBarrelHits;
      // looping through the RecHits of the track
      for (trackingRecHit_iterator irecHit = iTrack->recHitsBegin(); irecHit != iTrack->recHitsEnd(); ++irecHit) {
        if ((*irecHit)->isValid()) {
          DetId detId = (*irecHit)->geographicalId();
          // enum Detector { Tracker=1, Muon=2, Ecal=3, Hcal=4, Calo=5 };
          if (detId.det() != 1) {
            if (debug_) {
              cout << "rec hit ID = " << detId.det() << " not in tracker!?!?\n";
            }
            continue;
          }
          uint32_t subDet = detId.subdetId();

          // enum SubDetector{ PixelBarrel=1, PixelEndcap=2 };
          // enum SubDetector{ TIB=3, TID=4, TOB=5, TEC=6 };

          TransientTrackingRecHit::RecHitPointer trecHit = theTrackerRecHitBuilder->build(&*(*irecHit), initialTSOS);

          double gX = trecHit->globalPosition().x();
          double gY = trecHit->globalPosition().y();
          double gZ = trecHit->globalPosition().z();

          if (subDet == PixelSubdetector::PixelBarrel) {
            int ilay = tTopo->pxbLayer(detId);

            if (ilay == 1) {
              n1++;
              x1 = gX;
              y1 = gY;
              z1 = gZ;

              GoodPixBarrelHits.push_back((trecHit));
            }  // PXB1
            if (ilay == 2) {
              n2++;
              x2 = gX;
              y2 = gY;
              z2 = gZ;

              GoodPixBarrelHits.push_back((trecHit));

            }  // PXB2
            if (ilay == 3) {
              n3++;
              x3 = gX;
              y3 = gY;
              z3 = gZ;
              GoodPixBarrelHits.push_back((trecHit));
            }
          }  // PXB

        }  // valid
      }  // loop rechits

      // CS extra plots

      if (n1 + n2 + n3 == 3 && n1 * n2 * n3 > 0) {
        for (unsigned int i = 0; i < GoodPixBarrelHits.size(); i++) {
          if (GoodPixBarrelHits[i]->isValid()) {
            DetId detId = GoodPixBarrelHits[i]->geographicalId().rawId();
            int ilay = tTopo->pxbLayer(detId);
            if (pt > ptminres_) {
              double dca2 = 0.0, dz2 = 0.0;
              double ptsig = pt;
              if (charge < 0.)
                ptsig = -pt;
              // Filling the histograms in modules

              MeasurementPoint Test;
              MeasurementPoint Test2;
              Test = MeasurementPoint(0, 0);
              Test2 = MeasurementPoint(0, 0);
              Measurement2DVector residual;

              if (ilay == 1) {
                triplets(x2, y2, z2, x1, y1, z1, x3, y3, z3, ptsig, dca2, dz2, kap);

                Test = MeasurementPoint(dca2 * 1E4, dz2 * 1E4);
                residual = Test - Test2;
              }

              if (ilay == 2) {
                triplets(x1, y1, z1, x2, y2, z2, x3, y3, z3, ptsig, dca2, dz2, kap);

                Test = MeasurementPoint(dca2 * 1E4, dz2 * 1E4);
                residual = Test - Test2;
              }

              if (ilay == 3) {
                triplets(x1, y1, z1, x3, y3, z3, x2, y2, z2, ptsig, dca2, dz2, kap);

                Test = MeasurementPoint(dca2 * 1E4, dz2 * 1E4);
                residual = Test - Test2;
              }
              // fill the residual histograms

              std::map<uint32_t, SiPixelTrackResidualModule *>::iterator pxd = theSiPixelStructure.find(detId);
              if (pxd != theSiPixelStructure.end())
                (*pxd).second->fill(residual, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);

              if (ladOn && pxd != theSiPixelStructure.end()) {
                meResidualXSummedLay.at(PixelBarrelNameUpgrade((*pxd).first).layerName() - 1)->Fill(residual.x());
                meResidualYSummedLay.at(PixelBarrelNameUpgrade((*pxd).first).layerName() - 1)->Fill(residual.y());
              }

            }  // three hits
          }  // is valid
        }  // rechits loop
      }  // pt 4
    }

  }  //-----Tracks
  // get tracks
  edm::Handle<std::vector<reco::Track>> trackCollectionHandle;
  // iEvent.getByLabel(tracksrc_,trackCollectionHandle);
  iEvent.getByToken(trackToken_, trackCollectionHandle);
  auto const &trackColl = *(trackCollectionHandle.product());

  // get clusters
  edm::Handle<edmNew::DetSetVector<SiPixelCluster>> clusterColl;
  // iEvent.getByLabel( clustersrc_, clusterColl );
  iEvent.getByToken(clustersrcToken_, clusterColl);
  auto const &clustColl = *(clusterColl.product());

  // get digis
  edm::Handle<edm::DetSetVector<PixelDigi>> digiinput;
  iEvent.getByToken(digisrcToken_, digiinput);
  edm::DetSetVector<PixelDigi> const &diginp = *(digiinput.product());

  std::set<SiPixelCluster> clusterSet;
  TrajectoryStateCombiner tsoscomb;
  int tracks = 0, pixeltracks = 0, bpixtracks = 0, fpixtracks = 0;
  int trackclusters = 0, barreltrackclusters = 0, endcaptrackclusters = 0;
  int otherclusters = 0, barrelotherclusters = 0, endcapotherclusters = 0;

  // get trajectories
  edm::Handle<std::vector<Trajectory>> trajCollectionHandle;
  iEvent.getByToken(tracksrcToken_, trajCollectionHandle);

  if (debug_) {
    edm::EDConsumerBase::Labels labels;
    labelsForToken(tracksrcToken_, labels);
    std::cout << "Trajectories for Res from " << labels.module << std::endl;
  }

  if (trajCollectionHandle.isValid()) {
    auto const &trajColl = *(trajCollectionHandle.product());
    // get the map
    edm::Handle<TrajTrackAssociationCollection> match;
    iEvent.getByToken(trackAssociationToken_, match);
    auto const &ttac = *(match.product());

    if (debug_) {
      std::cout << "Trajectories\t : " << trajColl.size() << std::endl;
      std::cout << "recoTracks  \t : " << trackColl.size() << std::endl;
      std::cout << "Map entries \t : " << ttac.size() << std::endl;
    }

    // Loop over map entries
    for (TrajTrackAssociationCollection::const_iterator it = ttac.begin(); it != ttac.end(); ++it) {
      const edm::Ref<std::vector<Trajectory>> traj_iterator = it->key;
      // Trajectory Map, extract Trajectory for this track
      reco::TrackRef trackref = it->val;
      tracks++;

      bool isBpixtrack = false, isFpixtrack = false, crossesPixVol = false;

      // find out whether track crosses pixel fiducial volume (for cosmic
      // tracks)

      double d0 = (*trackref).d0(), dz = (*trackref).dz();

      if (abs(d0) < 15 && abs(dz) < 50)
        crossesPixVol = true;

      const std::vector<TrajectoryMeasurement> &tmeasColl = traj_iterator->measurements();
      std::vector<TrajectoryMeasurement>::const_iterator tmeasIt;
      // loop on measurements to find out whether there are bpix and/or fpix
      // hits
      for (tmeasIt = tmeasColl.begin(); tmeasIt != tmeasColl.end(); tmeasIt++) {
        if (!tmeasIt->updatedState().isValid())
          continue;
        TransientTrackingRecHit::ConstRecHitPointer testhit = tmeasIt->recHit();
        if (!testhit->isValid() || testhit->geographicalId().det() != DetId::Tracker)
          continue;
        uint testSubDetID = (testhit->geographicalId().subdetId());
        if (testSubDetID == PixelSubdetector::PixelBarrel)
          isBpixtrack = true;
        if (testSubDetID == PixelSubdetector::PixelEndcap)
          isFpixtrack = true;
      }  // end loop on measurements
      if (isBpixtrack) {
        bpixtracks++;
        if (debug_)
          std::cout << "bpixtrack\n";
      }
      if (isFpixtrack) {
        fpixtracks++;
        if (debug_)
          std::cout << "fpixtrack\n";
      }
      if (isBpixtrack || isFpixtrack) {
        pixeltracks++;

        if (crossesPixVol)
          meNofTracksInPixVol_->Fill(0, 1);

        const std::vector<TrajectoryMeasurement> &tmeasColl = traj_iterator->measurements();
        for (std::vector<TrajectoryMeasurement>::const_iterator tmeasIt = tmeasColl.begin(); tmeasIt != tmeasColl.end();
             tmeasIt++) {
          if (!tmeasIt->updatedState().isValid())
            continue;

          TrajectoryStateOnSurface tsos = tsoscomb(tmeasIt->forwardPredictedState(), tmeasIt->backwardPredictedState());
          if (!tsos.isValid())
            continue;  // Happens rarely, due to singular matrix or similar

          TransientTrackingRecHit::ConstRecHitPointer hit = tmeasIt->recHit();
          if (!hit->isValid() || hit->geographicalId().det() != DetId::Tracker) {
            continue;
          } else {
            //    //residual
            const DetId &hit_detId = hit->geographicalId();
            // uint IntRawDetID = (hit_detId.rawId());
            uint IntSubDetID = (hit_detId.subdetId());

            if (IntSubDetID == 0)
              continue;  // don't look at SiStrip hits!

            // get the enclosed persistent hit
            const TrackingRecHit *persistentHit = hit->hit();
            // check if it's not null, and if it's a valid pixel hit
            if ((persistentHit != nullptr) && (typeid(*persistentHit) == typeid(SiPixelRecHit))) {
              // tell the C++ compiler that the hit is a pixel hit
              const SiPixelRecHit *pixhit = static_cast<const SiPixelRecHit *>(hit->hit());
              // Hit probability:
              float hit_prob = -1.;
              if (pixhit->hasFilledProb()) {
                hit_prob = pixhit->clusterProbability(0);
                // std::cout<<"HITPROB= "<<hit_prob<<std::endl;
                if (hit_prob < pow(10., -15.))
                  NLowProb++;
                NTotal++;
                if (NTotal > 0)
                  meHitProbability->Fill(float(NLowProb / NTotal));
              }

              // get the edm::Ref to the cluster
              edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> const &clust = (*pixhit).cluster();

              //  check if the ref is not null
              if (clust.isNonnull()) {
                // define tracker and pixel geometry and topology
                const TrackerGeometry &theTracker(*theTrackerGeometry);
                const PixelGeomDetUnit *theGeomDet =
                    static_cast<const PixelGeomDetUnit *>(theTracker.idToDet(hit_detId));

                // test if PixelGeomDetUnit exists
                if (theGeomDet == nullptr) {
                  if (debug_)
                    std::cout << "NO THEGEOMDET\n";
                  continue;
                }

                const PixelTopology *topol = &(theGeomDet->specificTopology());
                // fill histograms for clusters on tracks
                // correct SiPixelTrackResidualModule
                std::map<uint32_t, SiPixelTrackResidualModule *>::iterator pxd =
                    theSiPixelStructure.find((*hit).geographicalId().rawId());

                // CHARGE CORRECTION (for track impact angle)
                // calculate alpha and beta from cluster position
                LocalTrajectoryParameters ltp = tsos.localParameters();
                LocalVector localDir = ltp.momentum() / ltp.momentum().mag();

                float clust_alpha = atan2(localDir.z(), localDir.x());
                float clust_beta = atan2(localDir.z(), localDir.y());
                double corrCharge = clust->charge() *
                                    sqrt(1.0 / (1.0 / pow(tan(clust_alpha), 2) + 1.0 / pow(tan(clust_beta), 2) + 1.0)) /
                                    1000.;

                if (pxd != theSiPixelStructure.end())
                  (*pxd).second->fill(
                      (*clust), true, corrCharge, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);

                trackclusters++;
                // CORR CHARGE
                meClChargeOnTrack_all->Fill(corrCharge);
                meClSizeOnTrack_all->Fill((*clust).size());
                meClSizeXOnTrack_all->Fill((*clust).sizeX());
                meClSizeYOnTrack_all->Fill((*clust).sizeY());
                clusterSet.insert(*clust);

                // find cluster global position (rphi, z)
                // get cluster center of gravity (of charge)
                float xcenter = clust->x();
                float ycenter = clust->y();
                // get the cluster position in local coordinates (cm)
                LocalPoint clustlp = topol->localPosition(MeasurementPoint(xcenter, ycenter));
                // get the cluster position in global coordinates (cm)
                GlobalPoint clustgp = theGeomDet->surface().toGlobal(clustlp);

                // find location of hit (barrel or endcap, same for cluster)
                bool barrel =
                    DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
                bool endcap =
                    DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
                if (barrel) {
                  barreltrackclusters++;

                  DetId detId = (*hit).geographicalId();
                  if (detId >= 302055684 && detId <= 352477708) {
                    getrococcupancy(detId, diginp, tTopo, meZeroRocLadvsModOnTrackBarrel);
                  }

                  // CORR CHARGE
                  meClChargeOnTrack_bpix->Fill(corrCharge);
                  meClSizeOnTrack_bpix->Fill((*clust).size());
                  meClSizeXOnTrack_bpix->Fill((*clust).sizeX());
                  meClSizeYOnTrack_bpix->Fill((*clust).sizeY());
                  int DBlayer;
                  DBlayer = PixelBarrelName(DetId((*hit).geographicalId()), tTopo, isUpgrade).layerName();
                  float phi = clustgp.phi();
                  float z = clustgp.z();

                  PixelBarrelName pbn(DetId((*hit).geographicalId()), tTopo, isUpgrade);
                  int ladder = pbn.ladderName();
                  int module = pbn.moduleName();

                  PixelBarrelName::Shell sh = pbn.shell();  // enum
                  int ladderSigned = ladder;
                  int moduleSigned = module;
                  // Shell { mO = 1, mI = 2 , pO =3 , pI =4 };
                  int shell = int(sh);
                  // change the module sign for z<0
                  if (shell == 1 || shell == 2) {
                    moduleSigned = -module;
                  }
                  // change ladeer sign for Outer )x<0)
                  if (shell == 1 || shell == 3) {
                    ladderSigned = -ladder;
                  }

                  for (int i = 0; i < noOfLayers; i++) {
                    if (DBlayer == i + 1) {
                      meClPosLayersOnTrack.at(i)->Fill(z, phi);
                      meClPosLayersLadVsModOnTrack.at(i)->Fill(moduleSigned, ladderSigned);
                      meClChargeOnTrack_layers.at(i)->Fill(corrCharge);
                      meClSizeOnTrack_layers.at(i)->Fill((*clust).size());
                      meClSizeXOnTrack_layers.at(i)->Fill((*clust).sizeX());
                      meClSizeYOnTrack_layers.at(i)->Fill((*clust).sizeY());
                      meNofClustersvsPhiOnTrack_layers.at(i)->Fill(phi);
                    }
                  }
                }
                if (endcap) {
                  endcaptrackclusters++;
                  // CORR CHARGE
                  meClChargeOnTrack_fpix->Fill(corrCharge);
                  meClSizeOnTrack_fpix->Fill((*clust).size());
                  meClSizeXOnTrack_fpix->Fill((*clust).sizeX());
                  meClSizeYOnTrack_fpix->Fill((*clust).sizeY());
                  int DBdisk = 0;
                  DBdisk = PixelEndcapName(DetId((*hit).geographicalId()), tTopo, isUpgrade).diskName();
                  float x = clustgp.x();
                  float y = clustgp.y();
                  float z = clustgp.z();
                  float phi = clustgp.phi();

                  float xclust = clust->x();
                  float yclust = clust->y();

                  getepixrococcupancyontrk(tTopo, hit, xclust, yclust, z, meRocBladevsDiskEndcapOnTrk);

                  if (z > 0) {
                    for (int i = 0; i < noOfDisks; i++) {
                      if (DBdisk == i + 1) {
                        meClPosDiskspzOnTrack.at(i)->Fill(x, y);
                        meClChargeOnTrack_diskps.at(i)->Fill(corrCharge);
                        meClSizeOnTrack_diskps.at(i)->Fill((*clust).size());
                        meClSizeXOnTrack_diskps.at(i)->Fill((*clust).sizeX());
                        meClSizeYOnTrack_diskps.at(i)->Fill((*clust).sizeY());
                        meNofClustersvsPhiOnTrack_diskps.at(i)->Fill(phi);
                      }
                    }
                  } else {
                    for (int i = 0; i < noOfDisks; i++) {
                      if (DBdisk == i + 1) {
                        meClPosDisksmzOnTrack.at(i)->Fill(x, y);
                        meClChargeOnTrack_diskms.at(i)->Fill(corrCharge);
                        meClSizeOnTrack_diskms.at(i)->Fill((*clust).size());
                        meClSizeXOnTrack_diskms.at(i)->Fill((*clust).sizeX());
                        meClSizeYOnTrack_diskms.at(i)->Fill((*clust).sizeY());
                        meNofClustersvsPhiOnTrack_diskms.at(i)->Fill(phi);
                      }
                    }
                  }
                }

              }  // end if (cluster exists)

            }  // end if (persistent hit exists and is pixel hit)

          }  // end of else

        }  // end for (all traj measurements of pixeltrack)
      }  // end if (is pixeltrack)
      else {
        if (debug_)
          std::cout << "no pixeltrack:\n";
        if (crossesPixVol)
          meNofTracksInPixVol_->Fill(1, 1);
      }

    }  // end loop on map entries

  }  // end valid trajectory:

  // find clusters that are NOT on track
  // edmNew::DetSet<SiPixelCluster>::const_iterator  di;
  if (debug_)
    std::cout << "clusters not on track: (size " << clustColl.size() << ") ";

  for (TrackerGeometry::DetContainer::const_iterator it = TG->dets().begin(); it != TG->dets().end(); it++) {
    // if(dynamic_cast<PixelGeomDetUnit const *>((*it))!=0){
    DetId detId = (*it)->geographicalId();
    if (detId >= 302055684 && detId <= 352477708) {  // make sure it's a Pixel module WITHOUT using
                                                     // dynamic_cast!
      int nofclOnTrack = 0, nofclOffTrack = 0;
      float z = 0.;
      edmNew::DetSetVector<SiPixelCluster>::const_iterator isearch = clustColl.find(detId);
      if (isearch != clustColl.end()) {  // Not an empty iterator
        edmNew::DetSet<SiPixelCluster>::const_iterator di;
        for (di = isearch->begin(); di != isearch->end(); di++) {
          unsigned int temp = clusterSet.size();
          clusterSet.insert(*di);
          // check if cluster is off track
          if (clusterSet.size() > temp) {
            otherclusters++;
            nofclOffTrack++;
            // fill histograms for clusters off tracks
            // correct SiPixelTrackResidualModule
            bool barrel = DetId(detId).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
            if (barrel) {
              getrococcupancy(detId, diginp, tTopo, meZeroRocLadvsModOffTrackBarrel);
            }

            std::map<uint32_t, SiPixelTrackResidualModule *>::iterator pxd =
                theSiPixelStructure.find((*it)->geographicalId().rawId());

            if (pxd != theSiPixelStructure.end())
              (*pxd).second->fill((*di), false, -1., reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);

            meClSizeNotOnTrack_all->Fill((*di).size());
            meClSizeXNotOnTrack_all->Fill((*di).sizeX());
            meClSizeYNotOnTrack_all->Fill((*di).sizeY());
            meClChargeNotOnTrack_all->Fill((*di).charge() / 1000);

            // find cluster global position (rphi, z) get cluster
            // define tracker and pixel geometry and topology
            const TrackerGeometry &theTracker(*theTrackerGeometry);
            const PixelGeomDetUnit *theGeomDet = static_cast<const PixelGeomDetUnit *>(theTracker.idToDet(detId));
            // test if PixelGeomDetUnit exists
            if (theGeomDet == nullptr) {
              if (debug_)
                std::cout << "NO THEGEOMDET\n";
              continue;
            }
            const PixelTopology *topol = &(theGeomDet->specificTopology());

            // center of gravity (of charge)
            float xcenter = di->x();
            float ycenter = di->y();
            // get the cluster position in local coordinates (cm)
            LocalPoint clustlp = topol->localPosition(MeasurementPoint(xcenter, ycenter));
            // get the cluster position in global coordinates (cm)
            GlobalPoint clustgp = theGeomDet->surface().toGlobal(clustlp);


            // barrel
            if (DetId(detId).subdetId() == 1) {
              meClSizeNotOnTrack_bpix->Fill((*di).size());
              meClSizeXNotOnTrack_bpix->Fill((*di).sizeX());
              meClSizeYNotOnTrack_bpix->Fill((*di).sizeY());
              meClChargeNotOnTrack_bpix->Fill((*di).charge() / 1000);
              barrelotherclusters++;
              int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
              float phi = clustgp.phi();
              // float r = clustgp.perp();
              z = clustgp.z();
              for (int i = 0; i < noOfLayers; i++) {
                if (DBlayer == i + 1) {
                  meClPosLayersNotOnTrack.at(i)->Fill(z, phi);
                  meClSizeNotOnTrack_layers.at(i)->Fill((*di).size());
                  meClSizeXNotOnTrack_layers.at(i)->Fill((*di).sizeX());
                  meClSizeYNotOnTrack_layers.at(i)->Fill((*di).sizeY());
                  meClChargeNotOnTrack_layers.at(i)->Fill((*di).charge() / 1000);
                }
              }
            }
            // endcap
            if (DetId(detId).subdetId() == 2) {
              meClSizeNotOnTrack_fpix->Fill((*di).size());
              meClSizeXNotOnTrack_fpix->Fill((*di).sizeX());
              meClSizeYNotOnTrack_fpix->Fill((*di).sizeY());
              meClChargeNotOnTrack_fpix->Fill((*di).charge() / 1000);
              endcapotherclusters++;
              int DBdisk = PixelEndcapName(DetId(detId), tTopo, isUpgrade).diskName();
              float x = clustgp.x();
              float y = clustgp.y();
              z = clustgp.z();

              getepixrococcupancyofftrk(DetId(detId), tTopo, xcenter, ycenter, z, meRocBladevsDiskEndcapOffTrk);

              if (z > 0) {
                for (int i = 0; i < noOfDisks; i++) {
                  if (DBdisk == i + 1) {
                    meClPosDiskspzNotOnTrack.at(i)->Fill(x, y);
                    meClSizeNotOnTrack_diskps.at(i)->Fill((*di).size());
                    meClSizeXNotOnTrack_diskps.at(i)->Fill((*di).sizeX());
                    meClSizeYNotOnTrack_diskps.at(i)->Fill((*di).sizeY());
                    meClChargeNotOnTrack_diskps.at(i)->Fill((*di).charge() / 1000);
                  }
                }
              } else {
                for (int i = 0; i < noOfDisks; i++) {
                  if (DBdisk == i + 1) {
                    meClPosDisksmzNotOnTrack.at(i)->Fill(x, y);
                    meClSizeNotOnTrack_diskms.at(i)->Fill((*di).size());
                    meClSizeXNotOnTrack_diskms.at(i)->Fill((*di).sizeX());
                    meClSizeYNotOnTrack_diskms.at(i)->Fill((*di).sizeY());
                    meClChargeNotOnTrack_diskms.at(i)->Fill((*di).charge() / 1000);
                  }
                }
              }
            }
          }  // end "if cluster off track"
          else {
            nofclOnTrack++;
            if (z == 0) {
              // find cluster global position (rphi, z) get cluster
              // define tracker and pixel geometry and topology
              const TrackerGeometry &theTracker(*theTrackerGeometry);
              const PixelGeomDetUnit *theGeomDet = static_cast<const PixelGeomDetUnit *>(theTracker.idToDet(detId));
              // test if PixelGeomDetUnit exists
              if (theGeomDet == nullptr) {
                if (debug_)
                  std::cout << "NO THEGEOMDET\n";
                continue;
              }
              const PixelTopology *topol = &(theGeomDet->specificTopology());
              // center of gravity (of charge)
              float xcenter = di->x();
              float ycenter = di->y();
              // get the cluster position in local coordinates (cm)
              LocalPoint clustlp = topol->localPosition(MeasurementPoint(xcenter, ycenter));
              // get the cluster position in global coordinates (cm)
              GlobalPoint clustgp = theGeomDet->surface().toGlobal(clustlp);
              z = clustgp.z();
            }
          }
        }
      }
      //++ fill the number of clusters on a module
      std::map<uint32_t, SiPixelTrackResidualModule *>::iterator pxd =
          theSiPixelStructure.find((*it)->geographicalId().rawId());
      if (pxd != theSiPixelStructure.end())
        (*pxd).second->nfill(
            nofclOnTrack, nofclOffTrack, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
      if (nofclOnTrack != 0)
        meNClustersOnTrack_all->Fill(nofclOnTrack);
      if (nofclOffTrack != 0)
        meNClustersNotOnTrack_all->Fill(nofclOffTrack);
      // barrel
      if (DetId(detId).subdetId() == 1) {
        if (nofclOnTrack != 0)
          meNClustersOnTrack_bpix->Fill(nofclOnTrack);
        if (nofclOffTrack != 0)
          meNClustersNotOnTrack_bpix->Fill(nofclOffTrack);
        int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
        for (int i = 0; i < noOfLayers; i++) {
          if (DBlayer == i + 1) {
            if (nofclOnTrack != 0)
              meNClustersOnTrack_layers.at(i)->Fill(nofclOnTrack);
            if (nofclOffTrack != 0)
              meNClustersNotOnTrack_layers.at(i)->Fill(nofclOffTrack);
          }
        }
      }  // end barrel
      // endcap
      if (DetId(detId).subdetId() == 2) {
        int DBdisk = PixelEndcapName(DetId(detId)).diskName();
        // z = clustgp.z();
        if (nofclOnTrack != 0)
          meNClustersOnTrack_fpix->Fill(nofclOnTrack);
        if (nofclOffTrack != 0)
          meNClustersNotOnTrack_fpix->Fill(nofclOffTrack);
        if (z > 0) {
          for (int i = 0; i < noOfDisks; i++) {
            if (DBdisk == i + 1) {
              if (nofclOnTrack != 0)
                meNClustersOnTrack_diskps.at(i)->Fill(nofclOnTrack);
              if (nofclOffTrack != 0)
                meNClustersNotOnTrack_diskps.at(i)->Fill(nofclOffTrack);
            }
          }
        }
        if (z < 0) {
          for (int i = 0; i < noOfDisks; i++) {
            if (DBdisk == i + 1) {
              if (nofclOnTrack != 0)
                meNClustersOnTrack_diskms.at(i)->Fill(nofclOnTrack);
              if (nofclOffTrack != 0)
                meNClustersNotOnTrack_diskms.at(i)->Fill(nofclOffTrack);
            }
          }
        }
      }

    }  // end if it's a Pixel module
  }  // end for loop over tracker detector geometry modules

  if (trackclusters > 0)
    (meNofClustersOnTrack_)->Fill(0, trackclusters);
  if (barreltrackclusters > 0)
    (meNofClustersOnTrack_)->Fill(1, barreltrackclusters);
  if (endcaptrackclusters > 0)
    (meNofClustersOnTrack_)->Fill(2, endcaptrackclusters);
  if (otherclusters > 0)
    (meNofClustersNotOnTrack_)->Fill(0, otherclusters);
  if (barrelotherclusters > 0)
    (meNofClustersNotOnTrack_)->Fill(1, barrelotherclusters);
  if (endcapotherclusters > 0)
    (meNofClustersNotOnTrack_)->Fill(2, endcapotherclusters);
  if (tracks > 0)
    (meNofTracks_)->Fill(0, tracks);
  if (pixeltracks > 0)
    (meNofTracks_)->Fill(1, pixeltracks);
  if (bpixtracks > 0)
    (meNofTracks_)->Fill(2, bpixtracks);
  if (fpixtracks > 0)
    (meNofTracks_)->Fill(3, fpixtracks);
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 149 of file SiPixelTrackResidualSource.cc.

References bladeOn, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), clustersrc_, debug_, digisrc_, diskOn, Exception, firstRun, edm::EventSetup::getHandle(), mps_fire::i, isUpgrade, edm::InputTag::label(), ladOn, layOn, meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskms, meClChargeNotOnTrack_diskps, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layers, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskms, meClChargeOnTrack_diskps, meClChargeOnTrack_fpix, meClChargeOnTrack_layers, meClPosDisksmzNotOnTrack, meClPosDisksmzOnTrack, meClPosDiskspzNotOnTrack, meClPosDiskspzOnTrack, meClPosLayersLadVsModOnTrack, meClPosLayersNotOnTrack, meClPosLayersOnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskms, meClSizeNotOnTrack_diskps, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layers, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskms, meClSizeOnTrack_diskps, meClSizeOnTrack_fpix, meClSizeOnTrack_layers, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskms, meClSizeXNotOnTrack_diskps, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layers, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskms, meClSizeXOnTrack_diskps, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layers, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskms, meClSizeYNotOnTrack_diskps, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layers, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskms, meClSizeYOnTrack_diskps, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layers, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskms, meNClustersNotOnTrack_diskps, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layers, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskms, meNClustersOnTrack_diskps, meNClustersOnTrack_fpix, meNClustersOnTrack_layers, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofClustersvsPhiOnTrack_diskms, meNofClustersvsPhiOnTrack_diskps, meNofClustersvsPhiOnTrack_layers, meNofTracks_, meNofTracksInPixVol_, meResidualXSummedLay, meResidualYSummedLay, meRocBladevsDiskEndcapOffTrk, meRocBladevsDiskEndcapOnTrk, meSubdetResidualX, meSubdetResidualY, meZeroRocLadvsModOffTrackBarrel, meZeroRocLadvsModOnTrackBarrel, modOn, hfnoseParametersInitialization_cfi::nameX, noOfDisks, noOfLayers, phiOn, edm::ESHandle< T >::product(), pSet_, reducedSet, ringOn, alignCSCRings::s, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theSiPixelStructure, topFolderName_, trackerTopoTokenBeginRun_, tracksrc_, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.

                                                                             {
  // book residual histograms in theSiPixelFolder - one (x,y) pair of histograms
  // per det
  SiPixelFolderOrganizer theSiPixelFolder(false);

  edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(trackerTopoTokenBeginRun_);
  const TrackerTopology *pTT = tTopoHandle.product();

  std::stringstream nameX, titleX, nameY, titleY;

  if (ladOn) {
    iBooker.setCurrentFolder(topFolderName_ + "/Barrel");
    for (int i = 1; i <= noOfLayers; i++) {
      nameX.str(std::string());
      nameX << "siPixelTrackResidualsX_SummedLayer_" << i;
      titleX.str(std::string());
      titleX << "Layer" << i << "Hit-to-Track Residual in r-phi";
      meResidualXSummedLay.push_back(iBooker.book1D(nameX.str(), titleX.str(), 100, -150, 150));
      meResidualXSummedLay.at(i - 1)->setAxisTitle("hit-to-track residual in r-phi (um)", 1);
      nameY.str(std::string());
      nameY << "siPixelTrackResidualsY_SummedLayer_" << i;
      titleY.str(std::string());
      titleY << "Layer" << i << "Hit-to-Track Residual in Z";
      meResidualYSummedLay.push_back(iBooker.book1D(nameY.str(), titleY.str(), 100, -300, 300));
      meResidualYSummedLay.at(i - 1)->setAxisTitle("hit-to-track residual in z (um)", 1);
    }
  }

  for (std::map<uint32_t, SiPixelTrackResidualModule *>::iterator pxd = theSiPixelStructure.begin();
       pxd != theSiPixelStructure.end();
       pxd++) {
    if (modOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 0, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 0, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Folder Creation Failed! ";
    }
    if (ladOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 1, isUpgrade)) {
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 1, isUpgrade);
      } else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource ladder Folder Creation Failed! ";
    }
    if (layOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 2, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 2, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource layer Folder Creation Failed! ";
    }
    if (phiOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 3, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 3, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource phi Folder Creation Failed! ";
    }
    if (bladeOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 4, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 4, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Blade Folder Creation Failed! ";
    }
    if (diskOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 5, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 5, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Disk Folder Creation Failed! ";
    }
    if (ringOn) {
      if (theSiPixelFolder.setModuleFolder(iBooker, (*pxd).first, 6, isUpgrade))
        (*pxd).second->book(pSet_, pTT, iBooker, reducedSet, 6, isUpgrade);
      else
        throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Ring Folder Creation Failed! ";
    }
  }

  //   edm::InputTag tracksrc =
  //   pSet_.getParameter<edm::InputTag>("trajectoryInput"); edm::InputTag
  //   clustersrc = pSet_.getParameter<edm::InputTag>("clustersrc");

  // number of tracks
  iBooker.setCurrentFolder(topFolderName_ + "/Tracks");
  meNofTracks_ = iBooker.book1D("ntracks_" + tracksrc_.label(), "Number of Tracks", 4, 0, 4);
  meNofTracks_->setAxisTitle("Number of Tracks", 1);
  meNofTracks_->setBinLabel(1, "All");
  meNofTracks_->setBinLabel(2, "Pixel");
  meNofTracks_->setBinLabel(3, "BPix");
  meNofTracks_->setBinLabel(4, "FPix");

  // number of tracks in pixel fiducial volume
  iBooker.setCurrentFolder(topFolderName_ + "/Tracks");
  meNofTracksInPixVol_ = iBooker.book1D(
      "ntracksInPixVol_" + tracksrc_.label(), "Number of Tracks crossing Pixel fiducial Volume", 2, 0, 2);
  meNofTracksInPixVol_->setAxisTitle("Number of Tracks", 1);
  meNofTracksInPixVol_->setBinLabel(1, "With Hits");
  meNofTracksInPixVol_->setBinLabel(2, "Without Hits");

  // number of clusters (associated to track / not associated)
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OnTrack");
  meNofClustersOnTrack_ =
      iBooker.book1D("nclusters_" + clustersrc_.label() + "_tot", "Number of Clusters (on track)", 3, 0, 3);
  meNofClustersOnTrack_->setAxisTitle("Number of Clusters on Track", 1);
  meNofClustersOnTrack_->setBinLabel(1, "All");
  meNofClustersOnTrack_->setBinLabel(2, "BPix");
  meNofClustersOnTrack_->setBinLabel(3, "FPix");
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OffTrack");
  meNofClustersNotOnTrack_ =
      iBooker.book1D("nclusters_" + clustersrc_.label() + "_tot", "Number of Clusters (off track)", 3, 0, 3);
  meNofClustersNotOnTrack_->setAxisTitle("Number of Clusters off Track", 1);
  meNofClustersNotOnTrack_->setBinLabel(1, "All");
  meNofClustersNotOnTrack_->setBinLabel(2, "BPix");
  meNofClustersNotOnTrack_->setBinLabel(3, "FPix");

  // cluster charge and size
  // charge
  // on track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OnTrack");
  std::stringstream ss1, ss2;
  meClChargeOnTrack_all = iBooker.book1D("charge_" + clustersrc_.label(), "Charge (on track)", 500, 0., 500.);
  meClChargeOnTrack_all->setAxisTitle("Charge size (in ke)", 1);
  meClChargeOnTrack_bpix =
      iBooker.book1D("charge_" + clustersrc_.label() + "_Barrel", "Charge (on track, barrel)", 500, 0., 500.);
  meClChargeOnTrack_bpix->setAxisTitle("Charge size (in ke)", 1);
  meClChargeOnTrack_fpix =
      iBooker.book1D("charge_" + clustersrc_.label() + "_Endcap", "Charge (on track, endcap)", 500, 0., 500.);
  meClChargeOnTrack_fpix->setAxisTitle("Charge size (in ke)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Charge (on track, layer" << i << ")";
    meClChargeOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeOnTrack_layers.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Charge (on track, diskp" << i << ")";
    meClChargeOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeOnTrack_diskps.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "Charge (on track, diskm" << i << ")";
    meClChargeOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeOnTrack_diskms.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }
  // off track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OffTrack");
  meClChargeNotOnTrack_all = iBooker.book1D("charge_" + clustersrc_.label(), "Charge (off track)", 500, 0., 500.);
  meClChargeNotOnTrack_all->setAxisTitle("Charge size (in ke)", 1);
  meClChargeNotOnTrack_bpix =
      iBooker.book1D("charge_" + clustersrc_.label() + "_Barrel", "Charge (off track, barrel)", 500, 0., 500.);
  meClChargeNotOnTrack_bpix->setAxisTitle("Charge size (in ke)", 1);
  meClChargeNotOnTrack_fpix =
      iBooker.book1D("charge_" + clustersrc_.label() + "_Endcap", "Charge (off track, endcap)", 500, 0., 500.);
  meClChargeNotOnTrack_fpix->setAxisTitle("Charge size (in ke)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Charge (off track, layer" << i << ")";
    meClChargeNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeNotOnTrack_layers.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Charge (off track, diskp" << i << ")";
    meClChargeNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeNotOnTrack_diskps.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "charge_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "Charge (off track, diskm" << i << ")";
    meClChargeNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClChargeNotOnTrack_diskms.at(i - 1)->setAxisTitle("Charge size (in ke)", 1);
  }

  // size
  // on track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OnTrack");
  meClSizeOnTrack_all = iBooker.book1D("size_" + clustersrc_.label(), "Size (on track)", 100, 0., 100.);
  meClSizeOnTrack_all->setAxisTitle("Cluster size (in pixels)", 1);
  meClSizeOnTrack_bpix =
      iBooker.book1D("size_" + clustersrc_.label() + "_Barrel", "Size (on track, barrel)", 100, 0., 100.);
  meClSizeOnTrack_bpix->setAxisTitle("Cluster size (in pixels)", 1);
  meClSizeOnTrack_fpix =
      iBooker.book1D("size_" + clustersrc_.label() + "_Endcap", "Size (on track, endcap)", 100, 0., 100.);
  meClSizeOnTrack_fpix->setAxisTitle("Cluster size (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Size (on track, layer" << i << ")";
    meClSizeOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Size (on track, diskp" << i << ")";
    meClSizeOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Disk_m1" << i;
    ss2.str(std::string());
    ss2 << "Size (on track, diskm" << i << ")";
    meClSizeOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  meClSizeXOnTrack_all = iBooker.book1D("sizeX_" + clustersrc_.label(), "SizeX (on track)", 100, 0., 100.);
  meClSizeXOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)", 1);
  meClSizeXOnTrack_bpix =
      iBooker.book1D("sizeX_" + clustersrc_.label() + "_Barrel", "SizeX (on track, barrel)", 100, 0., 100.);
  meClSizeXOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)", 1);
  meClSizeXOnTrack_fpix =
      iBooker.book1D("sizeX_" + clustersrc_.label() + "_Endcap", "SizeX (on track, endcap)", 100, 0., 100.);
  meClSizeXOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "SizeX (on track, layer" << i << ")";
    meClSizeXOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "SizeX (on track, diskp" << i << ")";
    meClSizeXOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "SizeX (on track, diskm" << i << ")";
    meClSizeXOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  meClSizeYOnTrack_all = iBooker.book1D("sizeY_" + clustersrc_.label(), "SizeY (on track)", 100, 0., 100.);
  meClSizeYOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)", 1);
  meClSizeYOnTrack_bpix =
      iBooker.book1D("sizeY_" + clustersrc_.label() + "_Barrel", "SizeY (on track, barrel)", 100, 0., 100.);
  meClSizeYOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)", 1);
  meClSizeYOnTrack_fpix =
      iBooker.book1D("sizeY_" + clustersrc_.label() + "_Endcap", "SizeY (on track, endcap)", 100, 0., 100.);
  meClSizeYOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "SizeY (on track, layer" << i << ")";
    meClSizeYOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "SizeY (on track, diskp" << i << ")";
    meClSizeYOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "SizeY (on track, diskm" << i << ")";
    meClSizeYOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  // off track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OffTrack");
  meClSizeNotOnTrack_all = iBooker.book1D("size_" + clustersrc_.label(), "Size (off track)", 100, 0., 100.);
  meClSizeNotOnTrack_all->setAxisTitle("Cluster size (in pixels)", 1);
  meClSizeNotOnTrack_bpix =
      iBooker.book1D("size_" + clustersrc_.label() + "_Barrel", "Size (off track, barrel)", 100, 0., 100.);
  meClSizeNotOnTrack_bpix->setAxisTitle("Cluster size (in pixels)", 1);
  meClSizeNotOnTrack_fpix =
      iBooker.book1D("size_" + clustersrc_.label() + "_Endcap", "Size (off track, endcap)", 100, 0., 100.);
  meClSizeNotOnTrack_fpix->setAxisTitle("Cluster size (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Size (off track, layer" << i << ")";
    meClSizeNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeNotOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }

  for (int i = 1; i <= noOfLayers; i++) {
    int ybins = -1;
    float ymin = 0.;
    float ymax = 0.;
    if (i == 1) {
      ybins = 42;
      ymin = -10.5;
      ymax = 10.5;
    }
    if (i == 2) {
      ybins = 66;
      ymin = -16.5;
      ymax = 16.5;
    }
    if (i == 3) {
      ybins = 90;
      ymin = -22.5;
      ymax = 22.5;
    }
    ss1.str(std::string());
    ss1 << "pix_bar Occ_roc_offtrack" + digisrc_.label() + "_layer_" << i;
    ss2.str(std::string());
    ss2 << "Pixel Barrel Occupancy, ROC level (Off Track): Layer " << i;
    meZeroRocLadvsModOffTrackBarrel.push_back(iBooker.book2D(ss1.str(), ss2.str(), 72, -4.5, 4.5, ybins, ymin, ymax));
    meZeroRocLadvsModOffTrackBarrel.at(i - 1)->setAxisTitle("ROC / Module", 1);
    meZeroRocLadvsModOffTrackBarrel.at(i - 1)->setAxisTitle("ROC / Ladder", 2);
  }

  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Size (off track, diskp" << i << ")";
    meClSizeNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeNotOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "size_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "Size (off track, diskm" << i << ")";
    meClSizeNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeNotOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  meClSizeXNotOnTrack_all = iBooker.book1D("sizeX_" + clustersrc_.label(), "SizeX (off track)", 100, 0., 100.);
  meClSizeXNotOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)", 1);
  meClSizeXNotOnTrack_bpix =
      iBooker.book1D("sizeX_" + clustersrc_.label() + "_Barrel", "SizeX (off track, barrel)", 100, 0., 100.);
  meClSizeXNotOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)", 1);
  meClSizeXNotOnTrack_fpix =
      iBooker.book1D("sizeX_" + clustersrc_.label() + "_Endcap", "SizeX (off track, endcap)", 100, 0., 100.);
  meClSizeXNotOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "SizeX (off track, layer" << i << ")";
    meClSizeXNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXNotOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "SizeX (off track, diskp" << i << ")";
    meClSizeXNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXNotOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeX_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "SizeX (off track, diskm" << i << ")";
    meClSizeXNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeXNotOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  meClSizeYNotOnTrack_all = iBooker.book1D("sizeY_" + clustersrc_.label(), "SizeY (off track)", 100, 0., 100.);
  meClSizeYNotOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)", 1);
  meClSizeYNotOnTrack_bpix =
      iBooker.book1D("sizeY_" + clustersrc_.label() + "_Barrel", "SizeY (off track, barrel)", 100, 0., 100.);
  meClSizeYNotOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)", 1);
  meClSizeYNotOnTrack_fpix =
      iBooker.book1D("sizeY_" + clustersrc_.label() + "_Endcap", "SizeY (off track, endcap)", 100, 0., 100.);
  meClSizeYNotOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "SizeY (off track, layer" << i << ")";
    meClSizeYNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYNotOnTrack_layers.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "SizeY (off track, diskp" << i << ")";
    meClSizeYNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYNotOnTrack_diskps.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "sizeY_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "SizeY (off track, diskm" << i << ")";
    meClSizeYNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meClSizeYNotOnTrack_diskms.at(i - 1)->setAxisTitle("Cluster size (in pixels)", 1);
  }

  // cluster global position
  // on track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OnTrack");
  // bpix
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Clusters Layer" << i << " (on track)";
    meClPosLayersOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 200, -30., 30., 128, -3.2, 3.2));
    meClPosLayersOnTrack.at(i - 1)->setAxisTitle("Global Z (cm)", 1);
    meClPosLayersOnTrack.at(i - 1)->setAxisTitle("Global #phi", 2);

    int ybins = -1;
    float ymin = 0.;
    float ymax = 0.;
    if (i == 1) {
      ybins = 23;
      ymin = -11.5;
      ymax = 11.5;
    }
    if (i == 2) {
      ybins = 33;
      ymin = -17.5;
      ymax = 17.5;
    }
    if (i == 3) {
      ybins = 45;
      ymin = -24.5;
      ymax = 24.5;
    }
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_LadvsMod_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Clusters Layer" << i << "_LadvsMod (on track)";
    meClPosLayersLadVsModOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 11, -5.5, 5.5, ybins, ymin, ymax));
    meClPosLayersLadVsModOnTrack.at(i - 1)->setAxisTitle("z-module", 1);
    meClPosLayersLadVsModOnTrack.at(i - 1)->setAxisTitle("Ladder", 2);
  }

  for (int i = 1; i <= noOfLayers; i++) {
    int ybins = -1;
    float ymin = 0.;
    float ymax = 0.;
    if (i == 1) {
      ybins = 42;
      ymin = -10.5;
      ymax = 10.5;
    }
    if (i == 2) {
      ybins = 66;
      ymin = -16.5;
      ymax = 16.5;
    }
    if (i == 3) {
      ybins = 90;
      ymin = -22.5;
      ymax = 22.5;
    }
    ss1.str(std::string());
    ss1 << "pix_bar Occ_roc_ontrack" + digisrc_.label() + "_layer_" << i;
    ss2.str(std::string());
    ss2 << "Pixel Barrel Occupancy, ROC level (On Track): Layer " << i;
    meZeroRocLadvsModOnTrackBarrel.push_back(iBooker.book2D(ss1.str(), ss2.str(), 72, -4.5, 4.5, ybins, ymin, ymax));
    meZeroRocLadvsModOnTrackBarrel.at(i - 1)->setAxisTitle("ROC / Module", 1);
    meZeroRocLadvsModOnTrackBarrel.at(i - 1)->setAxisTitle("ROC / Ladder", 2);
  }

  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "nclustersvsPhi_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "nclusters (on track, layer" << i << ")";
    meNofClustersvsPhiOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 1400., -3.5, 3.5));
    meNofClustersvsPhiOnTrack_layers.at(i - 1)->setAxisTitle("Global #Phi", 1);
    meNofClustersvsPhiOnTrack_layers.at(i - 1)->setAxisTitle("Number of Clusters/Layer on Track", 2);
  }

  // fpix
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_pz_Disk_" << i;
    ss2.str(std::string());
    ss2 << "Clusters +Z Disk" << i << " (on track)";
    meClPosDiskspzOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
    meClPosDiskspzOnTrack.at(i - 1)->setAxisTitle("Global X (cm)", 1);
    meClPosDiskspzOnTrack.at(i - 1)->setAxisTitle("Global Y (cm)", 2);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_mz_Disk_" << i;
    ss2.str(std::string());
    ss2 << "Clusters -Z Disk" << i << " (on track)";
    meClPosDisksmzOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
    meClPosDisksmzOnTrack.at(i - 1)->setAxisTitle("Global X (cm)", 1);
    meClPosDisksmzOnTrack.at(i - 1)->setAxisTitle("Global Y (cm)", 2);
  }
  meNClustersOnTrack_all =
      iBooker.book1D("nclusters_" + clustersrc_.label(), "Number of Clusters (on Track)", 50, 0., 50.);
  meNClustersOnTrack_all->setAxisTitle("Number of Clusters", 1);
  meNClustersOnTrack_bpix = iBooker.book1D(
      "nclusters_" + clustersrc_.label() + "_Barrel", "Number of Clusters (on track, barrel)", 50, 0., 50.);
  meNClustersOnTrack_bpix->setAxisTitle("Number of Clusters", 1);
  meNClustersOnTrack_fpix = iBooker.book1D(
      "nclusters_" + clustersrc_.label() + "_Endcap", "Number of Clusters (on track, endcap)", 50, 0., 50.);
  meNClustersOnTrack_fpix->setAxisTitle("Number of Clusters", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (on track, layer" << i << ")";
    meNClustersOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meNClustersOnTrack_layers.at(i - 1)->setAxisTitle("Number of Clusters", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (on track, diskp" << i << ")";
    meNClustersOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 50, 0., 50.));
    meNClustersOnTrack_diskps.at(i - 1)->setAxisTitle("Number of Clusters", 1);

    ss1.str(std::string());
    ss1 << "nclustersvsPhi_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "nclusters (on track, diskp" << i << ")";
    meNofClustersvsPhiOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 1400., -3.5, 3.5));
    meNofClustersvsPhiOnTrack_diskps.at(i - 1)->setAxisTitle("Global #Phi", 1);
    meNofClustersvsPhiOnTrack_diskps.at(i - 1)->setAxisTitle("Number of Clusters/Disk on Track", 2);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (on track, diskm" << i << ")";
    meNClustersOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meNClustersOnTrack_diskms.at(i - 1)->setAxisTitle("Number of Clusters", 1);

    ss1.str(std::string());
    ss1 << "nclustersvsPhi_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "nclusters (on track, diskm" << i << ")";
    meNofClustersvsPhiOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 1400., -3.5, 3.5));
    meNofClustersvsPhiOnTrack_diskms.at(i - 1)->setAxisTitle("Global #Phi", 1);
    meNofClustersvsPhiOnTrack_diskms.at(i - 1)->setAxisTitle("Number of Clusters/Disk on Track", 2);
  }

  meRocBladevsDiskEndcapOnTrk = iBooker.book2D(
      "ROC_endcap_occupancy_ontrk", "Pixel Endcap Occupancy, ROC level (On Track)", 72, -4.5, 4.5, 288, -12.5, 12.5);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(1, "Disk-2 Pnl2", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(9, "Disk-2 Pnl1", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(19, "Disk-1 Pnl2", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(27, "Disk-1 Pnl1", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(41, "Disk+1 Pnl1", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(49, "Disk+1 Pnl2", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(59, "Disk+2 Pnl1", 1);
  meRocBladevsDiskEndcapOnTrk->setBinLabel(67, "Disk+2 Pnl2", 1);
  meRocBladevsDiskEndcapOnTrk->setAxisTitle("Blades in Inner (>0) / Outer(<) Halves", 2);
  meRocBladevsDiskEndcapOnTrk->setAxisTitle("ROC occupancy", 3);

  // not on track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OffTrack");
  // bpix
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Clusters Layer" << i << " (off track)";
    meClPosLayersNotOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 200, -30., 30., 128, -3.2, 3.2));
    meClPosLayersNotOnTrack.at(i - 1)->setAxisTitle("Global Z (cm)", 1);
    meClPosLayersNotOnTrack.at(i - 1)->setAxisTitle("Global #phi", 2);
  }
  // fpix
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_pz_Disk_" << i;
    ss2.str(std::string());
    ss2 << "Clusters +Z Disk" << i << " (off track)";
    meClPosDiskspzNotOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
    meClPosDiskspzNotOnTrack.at(i - 1)->setAxisTitle("Global X (cm)", 1);
    meClPosDiskspzNotOnTrack.at(i - 1)->setAxisTitle("Global Y (cm)", 2);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "position_" + clustersrc_.label() + "_mz_Disk_" << i;
    ss2.str(std::string());
    ss2 << "Clusters -Z Disk" << i << " (off track)";
    meClPosDisksmzNotOnTrack.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
    meClPosDisksmzNotOnTrack.at(i - 1)->setAxisTitle("Global X (cm)", 1);
    meClPosDisksmzNotOnTrack.at(i - 1)->setAxisTitle("Global Y (cm)", 2);
  }
  meNClustersNotOnTrack_all =
      iBooker.book1D("nclusters_" + clustersrc_.label(), "Number of Clusters (off Track)", 50, 0., 50.);
  meNClustersNotOnTrack_all->setAxisTitle("Number of Clusters", 1);
  meNClustersNotOnTrack_bpix = iBooker.book1D(
      "nclusters_" + clustersrc_.label() + "_Barrel", "Number of Clusters (off track, barrel)", 50, 0., 50.);
  meNClustersNotOnTrack_bpix->setAxisTitle("Number of Clusters", 1);
  meNClustersNotOnTrack_fpix = iBooker.book1D(
      "nclusters_" + clustersrc_.label() + "_Endcap", "Number of Clusters (off track, endcap)", 50, 0., 50.);
  meNClustersNotOnTrack_fpix->setAxisTitle("Number of Clusters", 1);
  for (int i = 1; i <= noOfLayers; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Layer_" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (off track, layer" << i << ")";
    meNClustersNotOnTrack_layers.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meNClustersNotOnTrack_layers.at(i - 1)->setAxisTitle("Number of Clusters", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Disk_p" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (off track, diskp" << i << ")";
    meNClustersNotOnTrack_diskps.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meNClustersNotOnTrack_diskps.at(i - 1)->setAxisTitle("Number of Clusters", 1);
  }
  for (int i = 1; i <= noOfDisks; i++) {
    ss1.str(std::string());
    ss1 << "nclusters_" + clustersrc_.label() + "_Disk_m" << i;
    ss2.str(std::string());
    ss2 << "Number of Clusters (off track, diskm" << i << ")";
    meNClustersNotOnTrack_diskms.push_back(iBooker.book1D(ss1.str(), ss2.str(), 500, 0., 500.));
    meNClustersNotOnTrack_diskms.at(i - 1)->setAxisTitle("Number of Clusters", 1);
  }

  meRocBladevsDiskEndcapOffTrk = iBooker.book2D(
      "ROC_endcap_occupancy_offtrk", "Pixel Endcap Occupancy, ROC level (Off Track)", 72, -4.5, 4.5, 288, -12.5, 12.5);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(1, "Disk-2 Pnl2", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(9, "Disk-2 Pnl1", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(19, "Disk-1 Pnl2", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(27, "Disk-1 Pnl1", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(41, "Disk+1 Pnl1", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(49, "Disk+1 Pnl2", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(59, "Disk+2 Pnl1", 1);
  meRocBladevsDiskEndcapOffTrk->setBinLabel(67, "Disk+2 Pnl2", 1);
  meRocBladevsDiskEndcapOffTrk->setAxisTitle("Blades in Inner (>0) / Outer(<) Halves", 2);
  meRocBladevsDiskEndcapOffTrk->setAxisTitle("ROC occupancy", 3);

  // HitProbability
  // on track
  iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OnTrack");
  meHitProbability = iBooker.book1D(
      "FractionLowProb", "Fraction of hits with low probability;FractionLowProb;#HitsOnTrack", 100, 0., 1.);

  if (debug_) {
    // book summary residual histograms in a debugging folder - one (x,y) pair
    // of histograms per subdetector
    iBooker.setCurrentFolder("debugging");
    char hisID[80];
    for (int s = 0; s < 3; s++) {
      sprintf(hisID, "residual_x_subdet_%i", s);
      meSubdetResidualX[s] = iBooker.book1D(hisID, "Pixel Hit-to-Track Residual in X", 500, -5., 5.);

      sprintf(hisID, "residual_y_subdet_%i", s);
      meSubdetResidualY[s] = iBooker.book1D(hisID, "Pixel Hit-to-Track Residual in Y", 500, -5., 5.);
    }
  }

  firstRun = false;
}

dqmBeginRun()

void SiPixelTrackResidualSource::dqmBeginRun ( const edm::Run &  r, edm::EventSetup const &  iSetup  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 112 of file SiPixelTrackResidualSource.cc.

References debug_, TrackerGeometry::dets(), TrackerGeometry::detsPXB(), TrackerGeometry::detsPXF(), PixelEndcapName::diskName(), edm::EventSetup::getHandle(), isUpgrade, PixelBarrelName::layerName(), callgraph::module, noOfDisks, noOfLayers, edm::ESHandle< T >::product(), theSiPixelStructure, trackerGeomTokenBeginRun_, and trackerTopoTokenBeginRun_.

                                                                                         {
  LogInfo("PixelDQM") << "SiPixelTrackResidualSource beginRun()" << endl;
  // retrieve TrackerGeometry for pixel dets
  edm::ESHandle<TrackerGeometry> TG = iSetup.getHandle(trackerGeomTokenBeginRun_);
  if (debug_)
    LogVerbatim("PixelDQM") << "TrackerGeometry " << &(*TG) << " size is " << TG->dets().size() << endl;
  edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(trackerTopoTokenBeginRun_);
  const TrackerTopology *pTT = tTopoHandle.product();

  // build theSiPixelStructure with the pixel barrel and endcap dets from
  // TrackerGeometry
  for (TrackerGeometry::DetContainer::const_iterator pxb = TG->detsPXB().begin(); pxb != TG->detsPXB().end(); pxb++) {
    if (dynamic_cast<PixelGeomDetUnit const *>((*pxb)) != nullptr) {
      SiPixelTrackResidualModule *module = new SiPixelTrackResidualModule((*pxb)->geographicalId().rawId());
      theSiPixelStructure.insert(
          pair<uint32_t, SiPixelTrackResidualModule *>((*pxb)->geographicalId().rawId(), module));
      // int DBlayer = PixelBarrelNameWrapper(pSet_,
      // DetId((*pxb)->geographicalId())).layerName();
      int DBlayer = PixelBarrelName(DetId((*pxb)->geographicalId()), pTT, isUpgrade).layerName();
      if (noOfLayers < DBlayer)
        noOfLayers = DBlayer;
    }
  }
  for (TrackerGeometry::DetContainer::const_iterator pxf = TG->detsPXF().begin(); pxf != TG->detsPXF().end(); pxf++) {
    if (dynamic_cast<PixelGeomDetUnit const *>((*pxf)) != nullptr) {
      SiPixelTrackResidualModule *module = new SiPixelTrackResidualModule((*pxf)->geographicalId().rawId());
      theSiPixelStructure.insert(
          pair<uint32_t, SiPixelTrackResidualModule *>((*pxf)->geographicalId().rawId(), module));
      int DBdisk;
      DBdisk = PixelEndcapName(DetId((*pxf)->geographicalId()), pTT, isUpgrade).diskName();
      if (noOfDisks < DBdisk)
        noOfDisks = DBdisk;
    }
  }
  LogInfo("PixelDQM") << "SiPixelStructure size is " << theSiPixelStructure.size() << endl;
}

getepixrococcupancyofftrk()

void SiPixelTrackResidualSource::getepixrococcupancyofftrk ( DetId  detId, const TrackerTopology *const  tTopo, float  xclust, float  yclust, float  z, MonitorElement *  meinput  )

private

Definition at line 1827 of file SiPixelTrackResidualSource.cc.

References funct::abs(), PXFDetId::blade(), hcalRecHitTable_cff::detId, PXFDetId::disk(), dqm::impl::MonitorElement::getBinContent(), mps_fire::i, dqmiolumiharvest::j, PXFDetId::module(), PXFDetId::panel(), dqm::impl::MonitorElement::setBinContent(), and z.

Referenced by analyze().

                                                                                                             {
  PXFDetId pxfid = PXFDetId(detId);

  // int pxfside   = PixelEndcapName(detId,tTopo,isUpgrade).halfCylinder();
  int pxfpanel = pxfid.panel();
  int pxfmodule = pxfid.module();
  int pxfdisk = pxfid.disk();
  int pxfblade_off = pxfid.blade();

  if (z < 0.) {
    pxfdisk = -1. * pxfdisk;
  }

  int pxfblade = -99;
  if (pxfblade_off <= 6 && pxfblade_off >= 1) {
    pxfblade = 7 - pxfblade_off;
  } else if (pxfblade_off <= 18 && pxfblade_off >= 7) {
    pxfblade = 6 - pxfblade_off;
  } else if (pxfblade_off <= 24 && pxfblade_off >= 19) {
    pxfblade = 31 - pxfblade_off;
  }

  int clu_sdpx = ((pxfdisk > 0) ? 1 : -1) * (2 * (abs(pxfdisk) - 1) + pxfpanel);
  int binselx = (pxfpanel == 1 && (pxfmodule == 1 || pxfmodule == 4))
                    ? (pxfmodule == 1)
                    : ((pxfpanel == 1 && xclust < 80.0) || (pxfpanel == 2 && xclust >= 80.0));
  int nperpan = 2 * pxfmodule + pxfpanel - 1 + binselx;
  int clu_roc_binx =
      ((pxfdisk > 0) ? nperpan : 9 - nperpan) + (clu_sdpx + 4) * 8 - 2 * ((abs(pxfdisk) == 1) ? pxfdisk : 0);

  int clu_roc_biny = -99.;
  int nrocly = pxfmodule + pxfpanel;
  for (int i = 0; i < nrocly; i++) {
    int j = (pxfdisk < 0) ? i : nrocly - 1 - i;
    if (yclust >= (j * 52.0) && yclust < ((j + 1) * 52.0))
      clu_roc_biny = 6 - nrocly + 2 * i + ((pxfblade > 0) ? pxfblade - 1 : pxfblade + 12) * 12 + 1;
  }
  if (pxfblade > 0) {
    clu_roc_biny = clu_roc_biny + 144;
  }

  meinput->setBinContent(clu_roc_binx, clu_roc_biny, meinput->getBinContent(clu_roc_binx, clu_roc_biny) + 1);
  meinput->setBinContent(clu_roc_binx, clu_roc_biny + 1, meinput->getBinContent(clu_roc_binx, clu_roc_biny + 1) + 1);
}

getepixrococcupancyontrk()

void SiPixelTrackResidualSource::getepixrococcupancyontrk ( const TrackerTopology *const  tTopo, TransientTrackingRecHit::ConstRecHitPointer  hit, float  xclust, float  yclust, float  z, MonitorElement *  meinput  )

private

Definition at line 1773 of file SiPixelTrackResidualSource.cc.

References funct::abs(), dqm::impl::MonitorElement::getBinContent(), mps_fire::i, dqmiolumiharvest::j, TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), dqm::impl::MonitorElement::setBinContent(), and z.

Referenced by analyze().

                                                                                   {
  int pxfpanel = tTopo->pxfPanel((*hit).geographicalId());
  int pxfmodule = tTopo->pxfModule((*hit).geographicalId());
  int pxfdisk = tTopo->pxfDisk((*hit).geographicalId());
  int pxfblade_off = tTopo->pxfBlade((*hit).geographicalId());

  // translate to online conventions
  // each EPIX disk is split in 2 half-disk - each one consists of 12 blades
  // in offline: blades num 0->24; here translate numbering to online convetion
  // positive blades pointing to beam;  negative pointing away
  // each blade has two panels: each consisting of an array of ROC plaquettes
  // front (rear) pannel: 3 (4) plaquettes
  // number of ROCs in each plaquette depends on the position on the panel
  // each ROC has 80x52 pixel cells
  if (z < 0.) {
    pxfdisk = -1. * pxfdisk;
  }
  int pxfblade = -99;
  if (pxfblade_off <= 6 && pxfblade_off >= 1) {
    pxfblade = 7 - pxfblade_off;
  } else if (pxfblade_off <= 18 && pxfblade_off >= 7) {
    pxfblade = 6 - pxfblade_off;
  } else if (pxfblade_off <= 24 && pxfblade_off >= 19) {
    pxfblade = 31 - pxfblade_off;
  }

  int clu_sdpx = ((pxfdisk > 0) ? 1 : -1) * (2 * (abs(pxfdisk) - 1) + pxfpanel);
  int binselx = (pxfpanel == 1 && (pxfmodule == 1 || pxfmodule == 4))
                    ? (pxfmodule == 1)
                    : ((pxfpanel == 1 && xclust < 80.0) || (pxfpanel == 2 && xclust >= 80.0));
  int nperpan = 2 * pxfmodule + pxfpanel - 1 + binselx;
  int clu_roc_binx =
      ((pxfdisk > 0) ? nperpan : 9 - nperpan) + (clu_sdpx + 4) * 8 - 2 * ((abs(pxfdisk) == 1) ? pxfdisk : 0);

  int clu_roc_biny = -99.;
  int nrocly = pxfmodule + pxfpanel;
  for (int i = 0; i < nrocly; i++) {
    int j = (pxfdisk < 0) ? i : nrocly - 1 - i;
    if (yclust >= (j * 52.0) && yclust < ((j + 1) * 52.0))
      clu_roc_biny = 6 - nrocly + 2 * i + ((pxfblade > 0) ? pxfblade - 1 : pxfblade + 12) * 12 + 1;
  }
  if (pxfblade > 0) {
    clu_roc_biny = clu_roc_biny + 144;
  }

  meinput->setBinContent(clu_roc_binx, clu_roc_biny, meinput->getBinContent(clu_roc_binx, clu_roc_biny) + 1);
  meinput->setBinContent(clu_roc_binx, clu_roc_biny + 1, meinput->getBinContent(clu_roc_binx, clu_roc_biny + 1) + 1);
}

getrococcupancy()

void SiPixelTrackResidualSource::getrococcupancy	(	DetId	detId,
		const edm::DetSetVector< PixelDigi > &	diginp,
		const TrackerTopology *const	tTopo,
		std::vector< MonitorElement *>	meinput
	)

Definition at line 1596 of file SiPixelTrackResidualSource.cc.

References funct::abs(), edm::DetSetVector< T >::begin(), cuy::col, hcalRecHitTable_cff::detId, edm::DetSetVector< T >::end(), edm::DetSetVector< T >::find(), HLT_2024v14_cff::flip, ALCARECOEcalPhiSym_cff::float, PixelBarrelName::isHalfModule(), isUpgrade, PixelBarrelName::ladderName(), PixelBarrelName::layerName(), PixelBarrelName::moduleName(), and PixelBarrelName::shell().

Referenced by analyze().

                                                                                      {
  edm::DetSetVector<PixelDigi>::const_iterator ipxsearch = diginp.find(detId);
  if (ipxsearch != diginp.end()) {
    // Look at digis now
    edm::DetSet<PixelDigi>::const_iterator pxdi;
    for (pxdi = ipxsearch->begin(); pxdi != ipxsearch->end(); pxdi++) {
      bool isHalfModule = PixelBarrelName(DetId(detId), tTopo, isUpgrade).isHalfModule();
      int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
      int DBmodule = PixelBarrelName(DetId(detId), tTopo, isUpgrade).moduleName();
      int DBladder = PixelBarrelName(DetId(detId), tTopo, isUpgrade).ladderName();
      int DBshell = PixelBarrelName(DetId(detId), tTopo, isUpgrade).shell();

      // add sign to the modules
      if (DBshell == 1 || DBshell == 2) {
        DBmodule = -DBmodule;
      }
      if (DBshell == 1 || DBshell == 3) {
        DBladder = -DBladder;
      }

      int col = pxdi->column();
      int row = pxdi->row();

      float modsign = (float)DBmodule / (abs((float)DBmodule));
      float ladsign = (float)DBladder / (abs((float)DBladder));
      float rocx = ((float)col / (52. * 8.)) * modsign + ((float)DBmodule - (modsign) * 0.5);
      float rocy = ((float)row / (80. * 2.)) * ladsign + ((float)DBladder - (ladsign) * 0.5);

      // do the flip where need
      bool flip = false;
      if ((DBladder % 2 == 0) && (!isHalfModule)) {
        flip = true;
      }
      if ((flip) && (DBladder > 0)) {
        if ((((float)DBladder - (ladsign) * 0.5) <= rocy) && (rocy < (float)DBladder)) {
          rocy = rocy + ladsign * 0.5;
        } else if ((((float)DBladder) <= rocy) && (rocy < ((float)DBladder + (ladsign) * 0.5))) {
          rocy = rocy - ladsign * 0.5;
        }
      }

      // tweak border effect for negative modules/ladders
      if (modsign < 0) {
        rocx = rocx - 0.0001;
      }
      if (ladsign < 0) {
        rocy = rocy - 0.0001;
      } else {
        rocy = rocy + 0.0001;
      }
      if (abs(DBladder) == 1) {
        rocy = rocy + ladsign * 0.5;
      }  // take care of the half module

      if (DBlayer == 1) {
        meinput.at(0)->Fill(rocx, rocy);
      }
      if (DBlayer == 2) {
        meinput.at(1)->Fill(rocx, rocy);
      }
      if (DBlayer == 3) {
        meinput.at(2)->Fill(rocx, rocy);
      }
    }  // end of looping over pxdi
  }
}

triplets()

void SiPixelTrackResidualSource::triplets	(	double	x1,
		double	y1,
		double	z1,
		double	x2,
		double	y2,
		double	z2,
		double	x3,
		double	y3,
		double	z3,
		double	ptsig,
		double &	dc,
		double &	dz,
		double	kap
	)

Definition at line 1666 of file SiPixelTrackResidualSource.cc.

References gather_cfg::cout, dttmaxenums::L, EgammaValidation_cff::num, diffTwoXMLs::r1, rho, mathSSE::sqrt(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, geometryCSVtoXML::xx, testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, geometryCSVtoXML::yy, ppsModifySingularModes_cfi::z1, and ppsModifySingularModes_cfi::z2.

Referenced by analyze().

                                                      {
  // Define some constants
  using namespace std;

  // Curvature kap from global Track

  // inverse of the curvature is the radius in the transverse plane
  double rho = 1 / kap;
  // Check that the hits are in the correct layers
  double r1 = sqrt(x1 * x1 + y1 * y1);
  double r3 = sqrt(x3 * x3 + y3 * y3);

  if (r3 - r1 < 2.0)
    cout << "warn r1 = " << r1 << ", r3 = " << r3 << endl;

  // Calculate the centre of the helix in xy-projection with radius rho from the
  // track.
  // start with a line (sekante) connecting the two points (x1,y1) and (x3,y3)
  // vec_L = vec_x3-vec_x1 with L being the length of that vector.
  double L = sqrt((x3 - x1) * (x3 - x1) + (y3 - y1) * (y3 - y1));
  // lam is the line from the middel point of vec_q towards the center of the
  // circle X0,Y0
  // we already have kap and rho = 1/kap
  double lam = sqrt(rho * rho - L * L / 4);

  // There are two solutions, the sign of kap gives the information
  // which of them is correct.
  //
  if (kap > 0)
    lam = -lam;

  //
  // ( X0, Y0 ) is the centre of the circle that describes the helix in
  // xy-projection.
  //
  double x0 = 0.5 * (x1 + x3) + lam / L * (-y1 + y3);
  double y0 = 0.5 * (y1 + y3) + lam / L * (x1 - x3);

  // Calculate the dipangle in z direction (needed later for z residual) :
  // Starting from the heliz equation whihc has to hold for both points z1,z3
  double num = (y3 - y0) * (x1 - x0) - (x3 - x0) * (y1 - y0);
  double den = (x1 - x0) * (x3 - x0) + (y1 - y0) * (y3 - y0);
  double tandip = kap * (z3 - z1) / atan(num / den);

  // angle from first hit to dca point:
  //
  double dphi = atan(((x1 - x0) * y0 - (y1 - y0) * x0) / ((x1 - x0) * x0 + (y1 - y0) * y0));
  // z position of the track based on the middle of the circle
  // track equation for the z component
  double uz0 = z1 + tandip * dphi * rho;

  // RESIDUAL IN R-PHI
  // Calculate distance dca2 from point (x2,y2) to the circle which is given by
  // the distance of the point to the middlepoint dcM = sqrt((x0-x2)^2+(y0-y2))
  // and rho dca = rho +- dcM
  if (kap > 0)
    dca2 = rho - sqrt((x0 - x2) * (x0 - x2) + (y0 - y2) * (y0 - y2));
  else
    dca2 = rho + sqrt((-x0 + x2) * (-x0 + x2) + (-y0 + y2) * (-y0 + y2));

  // RESIDUAL IN Z
  double xx = 0;
  double yy = 0;
  // sign of kappa determines the calculation
  // xx and yy are the new coordinates starting from x2, y2 that are on the
  // track itself vec_X2+-dca2*vec(X0-X2)/|(X0-X2)|
  if (kap < 0) {
    xx = x2 + (dca2 * ((x0 - x2)) / sqrt((x0 - x2) * (x0 - x2) + (y0 - y2) * (y0 - y2)));
    yy = y2 + (dca2 * ((y0 - y2)) / sqrt((x0 - x2) * (x0 - x2) + (y0 - y2) * (y0 - y2)));
  } else if (kap >= 0) {
    xx = x2 - (dca2 * ((x0 - x2)) / sqrt((x0 - x2) * (x0 - x2) + (y0 - y2) * (y0 - y2)));
    yy = y2 - (dca2 * ((y0 - y2)) / sqrt((x0 - x2) * (x0 - x2) + (y0 - y2) * (y0 - y2)));
  }

  // to get residual in z start with calculating the new uz2 position if one has
  // moved to xx, yy on the track. First calculate the change in phi2 with
  // respect to the center X0, Y0
  double dphi2 = atan(((xx - x0) * y0 - (yy - y0) * x0) / ((xx - x0) * x0 + (yy - y0) * y0));
  // Solve track equation for this new z depending on the dip angle of the track
  // (see above calculated based on X1, X3 and X0, use uz0 as reference point
  // again.
  double uz2 = uz0 - dphi2 * tandip * rho;

  // subtract new z position from the old one
  dz2 = z2 - uz2;

  // if we are interested in the arclength this is unsigned though
  //  double cosphi2 = (x2*xx+y2*yy)/(sqrt(x2*x2+y2*y2)*sqrt(xx*xx+yy*yy));
  // double arcdca2=sqrt(x2*x2+y2*y2)*acos(cosphi2);
}

Member Data Documentation

beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> SiPixelTrackResidualSource::beamSpotToken_

private

Definition at line 80 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

bladeOn

bool SiPixelTrackResidualSource::bladeOn

private

Definition at line 104 of file SiPixelTrackResidualSource.h.

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

clustersrc_

edm::InputTag SiPixelTrackResidualSource::clustersrc_

private

Definition at line 77 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

clustersrcToken_

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > SiPixelTrackResidualSource::clustersrcToken_

private

Definition at line 86 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

debug_

bool SiPixelTrackResidualSource::debug_

private

Definition at line 98 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and SiPixelTrackResidualSource().

digisrc_

edm::InputTag SiPixelTrackResidualSource::digisrc_

private

Definition at line 88 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

digisrcToken_

edm::EDGetTokenT<edm::DetSetVector<PixelDigi> > SiPixelTrackResidualSource::digisrcToken_

private

Definition at line 89 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

diskOn

bool SiPixelTrackResidualSource::diskOn

private

Definition at line 104 of file SiPixelTrackResidualSource.h.

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

firstRun

bool SiPixelTrackResidualSource::firstRun

private

Definition at line 107 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

generalTracksToken_

edm::EDGetTokenT<reco::TrackCollection> SiPixelTrackResidualSource::generalTracksToken_

private

Definition at line 82 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

isUpgrade

bool SiPixelTrackResidualSource::isUpgrade

private

Definition at line 105 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and getrococcupancy().

ladOn

bool SiPixelTrackResidualSource::ladOn

private

Definition at line 102 of file SiPixelTrackResidualSource.h.

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

layOn

bool SiPixelTrackResidualSource::layOn

private

Definition at line 102 of file SiPixelTrackResidualSource.h.

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

meClChargeNotOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_all

private

Definition at line 132 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeNotOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_bpix

private

Definition at line 133 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeNotOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeNotOnTrack_diskms

private

Definition at line 137 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeNotOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeNotOnTrack_diskps

private

Definition at line 136 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeNotOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClChargeNotOnTrack_fpix

private

Definition at line 134 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeNotOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeNotOnTrack_layers

private

Definition at line 135 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_all

private

Definition at line 126 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_bpix

private

Definition at line 127 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeOnTrack_diskms

private

Definition at line 131 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeOnTrack_diskps

private

Definition at line 130 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClChargeOnTrack_fpix

private

Definition at line 128 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClChargeOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClChargeOnTrack_layers

private

Definition at line 129 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosDisksmzNotOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosDisksmzNotOnTrack

private

Definition at line 196 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosDisksmzOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosDisksmzOnTrack

private

Definition at line 194 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosDiskspzNotOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosDiskspzNotOnTrack

private

Definition at line 195 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosDiskspzOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosDiskspzOnTrack

private

Definition at line 193 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosLayersLadVsModOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosLayersLadVsModOnTrack

private

Definition at line 191 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosLayersNotOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosLayersNotOnTrack

private

Definition at line 192 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClPosLayersOnTrack

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClPosLayersOnTrack

private

Definition at line 190 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_all

private

Definition at line 144 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_bpix

private

Definition at line 145 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeNotOnTrack_diskms

private

Definition at line 149 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeNotOnTrack_diskps

private

Definition at line 148 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeNotOnTrack_fpix

private

Definition at line 146 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeNotOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeNotOnTrack_layers

private

Definition at line 147 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_all

private

Definition at line 138 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_bpix

private

Definition at line 139 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeOnTrack_diskms

private

Definition at line 143 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeOnTrack_diskps

private

Definition at line 142 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeOnTrack_fpix

private

Definition at line 140 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeOnTrack_layers

private

Definition at line 141 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_all

private

Definition at line 156 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_bpix

private

Definition at line 157 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXNotOnTrack_diskms

private

Definition at line 161 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXNotOnTrack_diskps

private

Definition at line 160 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeXNotOnTrack_fpix

private

Definition at line 158 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXNotOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXNotOnTrack_layers

private

Definition at line 159 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_all

private

Definition at line 150 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_bpix

private

Definition at line 151 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXOnTrack_diskms

private

Definition at line 155 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXOnTrack_diskps

private

Definition at line 154 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeXOnTrack_fpix

private

Definition at line 152 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeXOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeXOnTrack_layers

private

Definition at line 153 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_all

private

Definition at line 168 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_bpix

private

Definition at line 169 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYNotOnTrack_diskms

private

Definition at line 173 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYNotOnTrack_diskps

private

Definition at line 172 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeYNotOnTrack_fpix

private

Definition at line 170 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYNotOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYNotOnTrack_layers

private

Definition at line 171 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_all

private

Definition at line 162 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_bpix

private

Definition at line 163 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYOnTrack_diskms

private

Definition at line 167 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYOnTrack_diskps

private

Definition at line 166 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meClSizeYOnTrack_fpix

private

Definition at line 164 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meClSizeYOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meClSizeYOnTrack_layers

private

Definition at line 165 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meHitProbability

MonitorElement* SiPixelTrackResidualSource::meHitProbability

private

Definition at line 201 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_all

private

Definition at line 182 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_bpix

private

Definition at line 183 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersNotOnTrack_diskms

private

Definition at line 187 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersNotOnTrack_diskps

private

Definition at line 186 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meNClustersNotOnTrack_fpix

private

Definition at line 184 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersNotOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersNotOnTrack_layers

private

Definition at line 185 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_all

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_all

private

Definition at line 176 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_bpix

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_bpix

private

Definition at line 177 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersOnTrack_diskms

private

Definition at line 181 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersOnTrack_diskps

private

Definition at line 180 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_fpix

MonitorElement* SiPixelTrackResidualSource::meNClustersOnTrack_fpix

private

Definition at line 178 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNClustersOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNClustersOnTrack_layers

private

Definition at line 179 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofClustersNotOnTrack_

MonitorElement* SiPixelTrackResidualSource::meNofClustersNotOnTrack_

private

Definition at line 125 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofClustersOnTrack_

MonitorElement* SiPixelTrackResidualSource::meNofClustersOnTrack_

private

Definition at line 121 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofClustersvsPhiOnTrack_diskms

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNofClustersvsPhiOnTrack_diskms

private

Definition at line 124 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofClustersvsPhiOnTrack_diskps

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNofClustersvsPhiOnTrack_diskps

private

Definition at line 123 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofClustersvsPhiOnTrack_layers

std::vector<MonitorElement *> SiPixelTrackResidualSource::meNofClustersvsPhiOnTrack_layers

private

Definition at line 122 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofTracks_

MonitorElement* SiPixelTrackResidualSource::meNofTracks_

private

Definition at line 119 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meNofTracksInPixVol_

MonitorElement* SiPixelTrackResidualSource::meNofTracksInPixVol_

private

Definition at line 120 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meResidualXSummedLay

std::vector<MonitorElement *> SiPixelTrackResidualSource::meResidualXSummedLay

private

Definition at line 116 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meResidualYSummedLay

std::vector<MonitorElement *> SiPixelTrackResidualSource::meResidualYSummedLay

private

Definition at line 117 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meRocBladevsDiskEndcapOffTrk

MonitorElement* SiPixelTrackResidualSource::meRocBladevsDiskEndcapOffTrk

private

Definition at line 203 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meRocBladevsDiskEndcapOnTrk

MonitorElement* SiPixelTrackResidualSource::meRocBladevsDiskEndcapOnTrk

private

Definition at line 202 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meSubdetResidualX

MonitorElement* SiPixelTrackResidualSource::meSubdetResidualX[3]

private

Definition at line 113 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms().

meSubdetResidualY

MonitorElement* SiPixelTrackResidualSource::meSubdetResidualY[3]

private

Definition at line 114 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms().

meZeroRocLadvsModOffTrackBarrel

std::vector<MonitorElement *> SiPixelTrackResidualSource::meZeroRocLadvsModOffTrackBarrel

private

Definition at line 199 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

meZeroRocLadvsModOnTrackBarrel

std::vector<MonitorElement *> SiPixelTrackResidualSource::meZeroRocLadvsModOnTrackBarrel

private

Definition at line 198 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

modOn

bool SiPixelTrackResidualSource::modOn

private

Definition at line 99 of file SiPixelTrackResidualSource.h.

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

NLowProb

int SiPixelTrackResidualSource::NLowProb

private

Definition at line 109 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

noOfDisks

int SiPixelTrackResidualSource::noOfDisks

private

Definition at line 215 of file SiPixelTrackResidualSource.h.

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

noOfLayers

int SiPixelTrackResidualSource::noOfLayers

private

Definition at line 214 of file SiPixelTrackResidualSource.h.

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

NTotal

int SiPixelTrackResidualSource::NTotal

private

Definition at line 108 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

offlinePrimaryVerticesToken_

edm::EDGetTokenT<reco::VertexCollection> SiPixelTrackResidualSource::offlinePrimaryVerticesToken_

private

Definition at line 81 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

phiOn

bool SiPixelTrackResidualSource::phiOn

private

Definition at line 102 of file SiPixelTrackResidualSource.h.

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

pSet_

edm::ParameterSet SiPixelTrackResidualSource::pSet_

private

Definition at line 75 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

ptminres_

double SiPixelTrackResidualSource::ptminres_

private

Definition at line 106 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

reducedSet

bool SiPixelTrackResidualSource::reducedSet

private

Definition at line 100 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and bookHistograms().

ringOn

bool SiPixelTrackResidualSource::ringOn

private

Definition at line 104 of file SiPixelTrackResidualSource.h.

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

src_

edm::InputTag SiPixelTrackResidualSource::src_

private

Definition at line 76 of file SiPixelTrackResidualSource.h.

Referenced by SiPixelTrackResidualSource().

theSiPixelStructure

std::map<uint32_t, SiPixelTrackResidualModule *> SiPixelTrackResidualSource::theSiPixelStructure

private

Definition at line 111 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and ~SiPixelTrackResidualSource().

topFolderName_

std::string SiPixelTrackResidualSource::topFolderName_

Definition at line 72 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

trackAssociationToken_

edm::EDGetTokenT<TrajTrackAssociationCollection> SiPixelTrackResidualSource::trackAssociationToken_

private

Definition at line 85 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

trackerGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelTrackResidualSource::trackerGeomToken_

private

Definition at line 92 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

trackerGeomTokenBeginRun_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SiPixelTrackResidualSource::trackerGeomTokenBeginRun_

private

Definition at line 96 of file SiPixelTrackResidualSource.h.

Referenced by dqmBeginRun(), and SiPixelTrackResidualSource().

trackerTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelTrackResidualSource::trackerTopoToken_

private

Definition at line 91 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

trackerTopoTokenBeginRun_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiPixelTrackResidualSource::trackerTopoTokenBeginRun_

private

Definition at line 95 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), dqmBeginRun(), and SiPixelTrackResidualSource().

tracksrc_

edm::InputTag SiPixelTrackResidualSource::tracksrc_

private

Definition at line 78 of file SiPixelTrackResidualSource.h.

Referenced by bookHistograms(), and SiPixelTrackResidualSource().

tracksrcToken_

edm::EDGetTokenT<std::vector<Trajectory> > SiPixelTrackResidualSource::tracksrcToken_

private

Definition at line 83 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

trackToken_

edm::EDGetTokenT<std::vector<reco::Track> > SiPixelTrackResidualSource::trackToken_

private

Definition at line 84 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

transientTrackBuilderToken_

edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> SiPixelTrackResidualSource::transientTrackBuilderToken_

private

Definition at line 93 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

transientTrackingRecHitBuilderToken_

edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> SiPixelTrackResidualSource::transientTrackingRecHitBuilderToken_

private

Definition at line 94 of file SiPixelTrackResidualSource.h.

Referenced by analyze(), and SiPixelTrackResidualSource().

ttrhbuilder_

std::string SiPixelTrackResidualSource::ttrhbuilder_

private

Definition at line 79 of file SiPixelTrackResidualSource.h.

Referenced by SiPixelTrackResidualSource().

vtxsrc_

std::string SiPixelTrackResidualSource::vtxsrc_

private

Definition at line 87 of file SiPixelTrackResidualSource.h.

Referenced by SiPixelTrackResidualSource().