#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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 >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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.

     : 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;
 }

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_.

◆ ~SiPixelTrackResidualSource()

SiPixelTrackResidualSource::~SiPixelTrackResidualSource ( )

override

Definition at line 102 of file SiPixelTrackResidualSource.cc.

                                                         {
   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;
   }
 }

References theSiPixelStructure.

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.

                                                                                             {
   // 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;
  
   int kk = -1;
   //----------------------------------------------------------------------------
   // Residuals:
   //
   for (reco::TrackCollection::const_iterator iTrack = TracksForRes->begin(); iTrack != TracksForRes->end(); ++iTrack) {
     // count
     kk++;
     // 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);
 }

References funct::abs(), Reference_intrackfit_cff::barrel, beamSpotToken_, edmNew::DetSetVector< T >::begin(), bladeOn, TransientTrackingRecHitBuilder::build(), TransientTrackBuilder::build(), ALCARECOTkAlJpsiMuMu_cff::charge, SiPixelRecHit::clusterProbability(), clustersrcToken_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, gather_cfg::cout, d0, debug_, DetId::det(), 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(), edm::Ref< C, T, F >::key(), GetRecoTauVFromDQM_MC_cff::kk, 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, 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(), DetId::rawId(), reducedSet, ringOn, PixelBarrelName::shell(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), DetId::subdetId(), GeomDet::surface(), funct::tan(), groupFilesInBlocks::temp, theSiPixelStructure, Surface::toGlobal(), trackAssociationToken_, pfDisplacedTrackerVertex_cfi::trackColl, DetId::Tracker, trackerGeomToken_, trackerTopoToken_, tracks, tracksrcToken_, trackToken_, transientTrackBuilderToken_, transientTrackingRecHitBuilderToken_, FreeTrajectoryState::transverseCurvature(), triplets(), parallelization::uint, pwdgSkimBPark_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(), and testProducerWithPsetDescEmpty_cfi::z2.

◆ 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.

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

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.

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 112 of file SiPixelTrackResidualSource.cc.

                                                                                          {
   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;
 }

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

◆ getepixrococcupancyofftrk()

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

private

Definition at line 1829 of file SiPixelTrackResidualSource.cc.

                                                                                                              {
   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);
 }

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

Referenced by analyze().

◆ getepixrococcupancyontrk()

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

private

Definition at line 1775 of file SiPixelTrackResidualSource.cc.

                                                                                    {
   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);
 }

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().

◆ getrococcupancy()

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

Definition at line 1598 of file SiPixelTrackResidualSource.cc.

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

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

Referenced by analyze().

◆ 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 1668 of file SiPixelTrackResidualSource.cc.

                                                       {
   // 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);
 }