#include <SiPixelRecHitsValid.h>

Inheritance diagram for SiPixelRecHitsValid:

Public Member Functions
	SiPixelRecHitsValid (const edm::ParameterSet &conf)

	~SiPixelRecHitsValid () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	fillBarrel (const SiPixelRecHit &, const PSimHit &, DetId, const PixelGeomDetUnit , const TrackerTopology tTopo)

void	fillForward (const SiPixelRecHit &, const PSimHit &, DetId, const PixelGeomDetUnit , const TrackerTopology tTopo)

Private Attributes
MonitorElement *	clustChargeDisk1Plaquettes [7]

MonitorElement *	clustChargeDisk2Plaquettes [7]

MonitorElement *	clustChargeLayer1Modules [8]

MonitorElement *	clustChargeLayer2Modules [8]

MonitorElement *	clustChargeLayer3Modules [8]

MonitorElement *	clustXSizeDisk1Plaquettes [7]

MonitorElement *	clustXSizeDisk2Plaquettes [7]

MonitorElement *	clustXSizeLayer [3]

MonitorElement *	clustYSizeDisk1Plaquettes [7]

MonitorElement *	clustYSizeDisk2Plaquettes [7]

MonitorElement *	clustYSizeModule [8]

MonitorElement *	recHitBunchB

MonitorElement *	recHitBunchF

MonitorElement *	recHitEventB

MonitorElement *	recHitEventF

MonitorElement *	recHitNsimHitDisk1

MonitorElement *	recHitNsimHitDisk2

MonitorElement *	recHitNsimHitLayer [3]

MonitorElement *	recHitXFullModules

MonitorElement *	recHitXHalfModules

MonitorElement *	recHitXPlaquetteSize1

MonitorElement *	recHitXPlaquetteSize2

MonitorElement *	recHitXPullAllB

MonitorElement *	recHitXPullAllF

MonitorElement *	recHitXPullDisk1Plaquettes [7]

MonitorElement *	recHitXPullDisk2Plaquettes [7]

MonitorElement *	recHitXPullFlippedLadderLayers [3]

MonitorElement *	recHitXPullNonFlippedLadderLayers [3]

MonitorElement *	recHitXResAllB

MonitorElement *	recHitXResAllF

MonitorElement *	recHitXResDisk1Plaquettes [7]

MonitorElement *	recHitXResDisk2Plaquettes [7]

MonitorElement *	recHitXResFlippedLadderLayers [3]

MonitorElement *	recHitXResNonFlippedLadderLayers [3]

MonitorElement *	recHitYAllModules

MonitorElement *	recHitYPlaquetteSize2

MonitorElement *	recHitYPlaquetteSize3

MonitorElement *	recHitYPlaquetteSize4

MonitorElement *	recHitYPlaquetteSize5

MonitorElement *	recHitYPullAllB

MonitorElement *	recHitYPullAllF

MonitorElement *	recHitYPullDisk1Plaquettes [7]

MonitorElement *	recHitYPullDisk2Plaquettes [7]

MonitorElement *	recHitYPullLayer1Modules [8]

MonitorElement *	recHitYPullLayer2Modules [8]

MonitorElement *	recHitYPullLayer3Modules [8]

MonitorElement *	recHitYResAllB

MonitorElement *	recHitYResAllF

MonitorElement *	recHitYResDisk1Plaquettes [7]

MonitorElement *	recHitYResDisk2Plaquettes [7]

MonitorElement *	recHitYResLayer1Modules [8]

MonitorElement *	recHitYResLayer2Modules [8]

MonitorElement *	recHitYResLayer3Modules [8]

edm::EDGetTokenT< SiPixelRecHitCollection >	siPixelRecHitCollectionToken_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tGeomEsToken

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoEsToken

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

File: SiPixelRecHitsValid.h

Author: Jason Shaev, JHU Created: 6/7/06

Definition at line 29 of file SiPixelRecHitsValid.h.

Constructor & Destructor Documentation

◆ SiPixelRecHitsValid()

SiPixelRecHitsValid::SiPixelRecHitsValid ( const edm::ParameterSet & conf )

Definition at line 46 of file SiPixelRecHitsValid.cc.

     : tTopoEsToken(esConsumes()),
       tGeomEsToken(esConsumes()),
       trackerHitAssociatorConfig_(ps, consumesCollector()),
       siPixelRecHitCollectionToken_(consumes<SiPixelRecHitCollection>(ps.getParameter<edm::InputTag>("src"))) {}

◆ ~SiPixelRecHitsValid()

SiPixelRecHitsValid::~SiPixelRecHitsValid ( )

overridedefault

Member Function Documentation

◆ analyze()

void SiPixelRecHitsValid::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 275 of file SiPixelRecHitsValid.cc.

References TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), gather_cfg::cout, TrackerGeometry::dets(), MillePedeFileConverter_cfg::e, edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), fillBarrel(), fillForward(), edmNew::DetSetVector< T >::find(), edm::EventSetup::getData(), mps_fire::i, TrackerGeometry::idToDet(), visualization-live-secondInstance_cfg::m, muonTagProbeFilters_cff::matched, TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), recHitNsimHitDisk1, recHitNsimHitDisk2, recHitNsimHitLayer, siPixelRecHitCollectionToken_, mathSSE::sqrt(), DetId::subdetId(), tGeomEsToken, trackerHitAssociatorConfig_, tTopoEsToken, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                             {
   //Retrieve tracker topology from geometry
   const TrackerTopology* tTopo = &es.getData(tTopoEsToken);
 
   edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
   if (e.id().event() % 1000 == 0)
     std::cout << " Run = " << e.id().run() << " Event = " << e.id().event() << std::endl;
 
   //Get RecHits
   edm::Handle<SiPixelRecHitCollection> recHitColl;
   e.getByToken(siPixelRecHitCollectionToken_, recHitColl);
 
   //Get event setup
   const TrackerGeometry* theTracker = &es.getData(tGeomEsToken);
 
   TrackerHitAssociator associate(e, trackerHitAssociatorConfig_);
 
   //iterate over detunits
   //for (TrackerGeometry::DetContainer::const_iterator it = geom->dets().begin(); it != geom->dets().end(); it++) {
   for (const auto& it : theTracker->dets()) {
     DetId detId = it->geographicalId();
     unsigned int subid = detId.subdetId();
 
     if (!((subid == 1) || (subid == 2)))
       continue;
 
     const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(theTracker->idToDet(detId));
 
     SiPixelRecHitCollection::const_iterator pixeldet = recHitColl->find(detId);
     if (pixeldet == recHitColl->end())
       continue;
     SiPixelRecHitCollection::DetSet pixelrechitRange = *pixeldet;
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
     SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
     std::vector<PSimHit> matched;
 
     //----Loop over rechits for this detId
     for (; pixeliter != pixelrechitRangeIteratorEnd; pixeliter++) {
       matched.clear();
       matched = associate.associateHit(*pixeliter);
 
       if (!matched.empty()) {
         float closest = 9999.9;
         std::vector<PSimHit>::const_iterator closestit = matched.begin();
         LocalPoint lp = pixeliter->localPosition();
         float rechit_x = lp.x();
         float rechit_y = lp.y();
 
         //loop over sim hits and fill closet
         for (std::vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); m++) {
           float sim_x1 = (*m).entryPoint().x();
           float sim_x2 = (*m).exitPoint().x();
           float sim_xpos = 0.5 * (sim_x1 + sim_x2);
 
           float sim_y1 = (*m).entryPoint().y();
           float sim_y2 = (*m).exitPoint().y();
           float sim_ypos = 0.5 * (sim_y1 + sim_y2);
 
           float x_res = fabs(sim_xpos - rechit_x);
           float y_res = fabs(sim_ypos - rechit_y);
 
           float dist = sqrt(x_res * x_res + y_res * y_res);
 
           if (dist < closest) {
             closest = dist;
             closestit = m;
           }
         }  // end sim hit loop
 
         if (subid == 1) {  //<----------barrel
           fillBarrel(*pixeliter, *closestit, detId, theGeomDet, tTopo);
         }                  // end barrel
         if (subid == 2) {  // <-------forward
           fillForward(*pixeliter, *closestit, detId, theGeomDet, tTopo);
         }
 
       }  // end matched emtpy
 
       int NsimHit = matched.size();
       if (subid == 1) {  //<----------barrel
         for (unsigned int i = 0; i < 3; i++)
           if (tTopo->pxbLayer(detId) == i + 1)
             recHitNsimHitLayer[i]->Fill(NsimHit);
       }                  // end barrel
       if (subid == 2) {  // <-------forward
         if (tTopo->pxfDisk(detId) == 1)
           recHitNsimHitDisk1->Fill(NsimHit);
         else
           recHitNsimHitDisk2->Fill(NsimHit);
       }
     }  // <-----end rechit loop
   }    // <------ end detunit loop
 }

◆ bookHistograms()

void SiPixelRecHitsValid::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const edm::Run &	run,
		const edm::EventSetup &	es
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 54 of file SiPixelRecHitsValid.cc.

References dqm::implementation::IBooker::book1D(), clustChargeDisk1Plaquettes, clustChargeDisk2Plaquettes, clustChargeLayer1Modules, clustChargeLayer2Modules, clustChargeLayer3Modules, clustXSizeDisk1Plaquettes, clustXSizeDisk2Plaquettes, clustXSizeLayer, clustYSizeDisk1Plaquettes, clustYSizeDisk2Plaquettes, clustYSizeModule, timingPdfMaker::histo, mps_fire::i, recHitBunchB, recHitBunchF, recHitEventB, recHitEventF, recHitNsimHitDisk1, recHitNsimHitDisk2, recHitNsimHitLayer, recHitXFullModules, recHitXHalfModules, recHitXPlaquetteSize1, recHitXPlaquetteSize2, recHitXPullAllB, recHitXPullAllF, recHitXPullDisk1Plaquettes, recHitXPullDisk2Plaquettes, recHitXPullFlippedLadderLayers, recHitXPullNonFlippedLadderLayers, recHitXResAllB, recHitXResAllF, recHitXResDisk1Plaquettes, recHitXResDisk2Plaquettes, recHitXResFlippedLadderLayers, recHitXResNonFlippedLadderLayers, recHitYAllModules, recHitYPlaquetteSize2, recHitYPlaquetteSize3, recHitYPlaquetteSize4, recHitYPlaquetteSize5, recHitYPullAllB, recHitYPullAllF, recHitYPullDisk1Plaquettes, recHitYPullDisk2Plaquettes, recHitYPullLayer1Modules, recHitYPullLayer2Modules, recHitYPullLayer3Modules, recHitYResAllB, recHitYResAllF, recHitYResDisk1Plaquettes, recHitYResDisk2Plaquettes, recHitYResLayer1Modules, recHitYResLayer2Modules, recHitYResLayer3Modules, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                              {
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/clustBPIX");
 
   Char_t histo[200];
 
   // ---------------------------------------------------------------
   // All histograms that depend on plaquette number have 7 indexes.
   // The first 4 (0-3) correspond to Panel 1 plaquettes 1-4.
   // The last 3 (4-6) correspond to Panel 2 plaquettes 1-3.
   // ---------------------------------------------------------------
 
   //Cluster y-size by module number for barrel
   for (int i = 0; i < 8; i++) {
     sprintf(histo, "Clust_y_size_Module%d", i + 1);
     clustYSizeModule[i] = ibooker.book1D(histo, "Cluster y-size by Module", 20, 0.5, 20.5);
   }  // end for
 
   //Cluster x-size by layer for barrel
   for (int i = 0; i < 3; i++) {
     sprintf(histo, "Clust_x_size_Layer%d", i + 1);
     clustXSizeLayer[i] = ibooker.book1D(histo, "Cluster x-size by Layer", 20, 0.5, 20.5);
   }  // end for
 
   //Cluster charge by module for 3 layers of barrel
   for (int i = 0; i < 8; i++) {
     //Cluster charge by module for Layer1
     sprintf(histo, "Clust_charge_Layer1_Module%d", i + 1);
     clustChargeLayer1Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 1 by Module", 50, 0., 200000.);
 
     //Cluster charge by module for Layer2
     sprintf(histo, "Clust_charge_Layer2_Module%d", i + 1);
     clustChargeLayer2Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 2 by Module", 50, 0., 200000.);
 
     //Cluster charge by module for Layer3
     sprintf(histo, "Clust_charge_Layer3_Module%d", i + 1);
     clustChargeLayer3Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 3 by Module", 50, 0., 200000.);
   }  // end for
 
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/clustFPIX");
   //Cluster x-size, y-size, and charge by plaquette for Disks in Forward
   for (int i = 0; i < 7; i++) {
     //Cluster x-size for Disk1 by Plaquette
     sprintf(histo, "Clust_x_size_Disk1_Plaquette%d", i + 1);
     clustXSizeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster X-size for Disk1 by Plaquette", 20, 0.5, 20.5);
 
     //Cluster x-size for Disk2 by Plaquette
     sprintf(histo, "Clust_x_size_Disk2_Plaquette%d", i + 1);
     clustXSizeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster X-size for Disk2 by Plaquette", 20, 0.5, 20.5);
 
     //Cluster y-size for Disk1 by Plaquette
     sprintf(histo, "Clust_y_size_Disk1_Plaquette%d", i + 1);
     clustYSizeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster Y-size for Disk1 by Plaquette", 20, 0.5, 20.5);
 
     //Cluster y-size for Disk2 by Plaquette
     sprintf(histo, "Clust_y_size_Disk2_Plaquette%d", i + 1);
     clustYSizeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster Y-size for Disk2 by Plaquette", 20, 0.5, 20.5);
 
     //Cluster charge for Disk1 by Plaquette
     sprintf(histo, "Clust_charge_Disk1_Plaquette%d", i + 1);
     clustChargeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster charge for Disk1 by Plaquette", 50, 0., 200000.);
 
     //Cluster charge for Disk2 by Plaquette
     sprintf(histo, "Clust_charge_Disk2_Plaquette%d", i + 1);
     clustChargeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster charge for Disk2 by Plaquette", 50, 0., 200000.);
   }  // end for
 
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitBPIX");
   //RecHit Bunch crossing all barrel hits
   recHitBunchB = ibooker.book1D("RecHit_Bunch_Barrel", "RecHit Bunch Crossing, Barrel", 20, -10., 10.);
 
   //RecHit Event, in-time bunch, all barrel hits
   recHitEventB = ibooker.book1D("RecHit_Event_Barrel", "RecHit Event (in-time bunch), Barrel", 100, 0., 100.);
 
   //RecHit X Resolution all barrel hits
   recHitXResAllB = ibooker.book1D("RecHit_xres_b_All", "RecHit X Res All Modules in Barrel", 100, -200., 200.);
 
   //RecHit Y Resolution all barrel hits
   recHitYResAllB = ibooker.book1D("RecHit_yres_b_All", "RecHit Y Res All Modules in Barrel", 100, -200., 200.);
 
   //RecHit X distribution for full modules for barrel
   recHitXFullModules = ibooker.book1D("RecHit_x_FullModules", "RecHit X distribution for full modules", 100, -2., 2.);
 
   //RecHit X distribution for half modules for barrel
   recHitXHalfModules = ibooker.book1D("RecHit_x_HalfModules", "RecHit X distribution for half modules", 100, -1., 1.);
 
   //RecHit Y distribution all modules for barrel
   recHitYAllModules = ibooker.book1D("RecHit_y_AllModules", "RecHit Y distribution for all modules", 100, -4., 4.);
 
   //RecHit X resolution for flipped and unflipped ladders by layer for barrel
   for (int i = 0; i < 3; i++) {
     //RecHit no. of matched simHits all ladders by layer
     sprintf(histo, "RecHit_NsimHit_Layer%d", i + 1);
     recHitNsimHitLayer[i] = ibooker.book1D(histo, "RecHit Number of simHits by Layer", 30, 0., 30.);
 
     //RecHit X resolution for flipped ladders by layer
     sprintf(histo, "RecHit_XRes_FlippedLadder_Layer%d", i + 1);
     recHitXResFlippedLadderLayers[i] = ibooker.book1D(histo, "RecHit XRes Flipped Ladders by Layer", 100, -200., 200.);
 
     //RecHit X resolution for unflipped ladders by layer
     sprintf(histo, "RecHit_XRes_UnFlippedLadder_Layer%d", i + 1);
     recHitXResNonFlippedLadderLayers[i] =
         ibooker.book1D(histo, "RecHit XRes NonFlipped Ladders by Layer", 100, -200., 200.);
   }  // end for
 
   //RecHit Y resolutions for layers by module for barrel
   for (int i = 0; i < 8; i++) {
     //Rec Hit Y resolution by module for Layer1
     sprintf(histo, "RecHit_YRes_Layer1_Module%d", i + 1);
     recHitYResLayer1Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer1 by module", 100, -200., 200.);
 
     //RecHit Y resolution by module for Layer2
     sprintf(histo, "RecHit_YRes_Layer2_Module%d", i + 1);
     recHitYResLayer2Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer2 by module", 100, -200., 200.);
 
     //RecHit Y resolution by module for Layer3
     sprintf(histo, "RecHit_YRes_Layer3_Module%d", i + 1);
     recHitYResLayer3Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer3 by module", 100, -200., 200.);
   }  // end for
 
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitFPIX");
   //RecHit Bunch crossing all plaquettes
   recHitBunchF = ibooker.book1D("RecHit_Bunch_Forward", "RecHit Bunch Crossing, Forward", 20, -10., 10.);
 
   //RecHit Event, in-time bunch, all plaquettes
   recHitEventF = ibooker.book1D("RecHit_Event_Forward", "RecHit Event (in-time bunch), Forward", 100, 0., 100.);
 
   //RecHit No. of simHits, by disk
   recHitNsimHitDisk1 = ibooker.book1D("RecHit_NsimHit_Disk1", "RecHit Number of simHits, Disk1", 30, 0., 30.);
   recHitNsimHitDisk2 = ibooker.book1D("RecHit_NsimHit_Disk2", "RecHit Number of simHits, Disk2", 30, 0., 30.);
 
   //RecHit X resolution all plaquettes
   recHitXResAllF = ibooker.book1D("RecHit_xres_f_All", "RecHit X Res All in Forward", 100, -200., 200.);
 
   //RecHit Y resolution all plaquettes
   recHitYResAllF = ibooker.book1D("RecHit_yres_f_All", "RecHit Y Res All in Forward", 100, -200., 200.);
 
   //RecHit X distribution for plaquette with x-size 1 in forward
   recHitXPlaquetteSize1 =
       ibooker.book1D("RecHit_x_Plaquette_xsize1", "RecHit X Distribution for plaquette x-size1", 100, -2., 2.);
 
   //RecHit X distribution for plaquette with x-size 2 in forward
   recHitXPlaquetteSize2 =
       ibooker.book1D("RecHit_x_Plaquette_xsize2", "RecHit X Distribution for plaquette x-size2", 100, -2., 2.);
 
   //RecHit Y distribution for plaquette with y-size 2 in forward
   recHitYPlaquetteSize2 =
       ibooker.book1D("RecHit_y_Plaquette_ysize2", "RecHit Y Distribution for plaquette y-size2", 100, -4., 4.);
 
   //RecHit Y distribution for plaquette with y-size 3 in forward
   recHitYPlaquetteSize3 =
       ibooker.book1D("RecHit_y_Plaquette_ysize3", "RecHit Y Distribution for plaquette y-size3", 100, -4., 4.);
 
   //RecHit Y distribution for plaquette with y-size 4 in forward
   recHitYPlaquetteSize4 =
       ibooker.book1D("RecHit_y_Plaquette_ysize4", "RecHit Y Distribution for plaquette y-size4", 100, -4., 4.);
 
   //RecHit Y distribution for plaquette with y-size 5 in forward
   recHitYPlaquetteSize5 =
       ibooker.book1D("RecHit_y_Plaquette_ysize5", "RecHit Y Distribution for plaquette y-size5", 100, -4., 4.);
 
   //X and Y resolutions for both disks by plaquette in forward
   for (int i = 0; i < 7; i++) {
     //X resolution for Disk1 by plaquette
     sprintf(histo, "RecHit_XRes_Disk1_Plaquette%d", i + 1);
     recHitXResDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit XRes Disk1 by plaquette", 100, -200., 200.);
     //X resolution for Disk2 by plaquette
     sprintf(histo, "RecHit_XRes_Disk2_Plaquette%d", i + 1);
     recHitXResDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit XRes Disk2 by plaquette", 100, -200., 200.);
 
     //Y resolution for Disk1 by plaquette
     sprintf(histo, "RecHit_YRes_Disk1_Plaquette%d", i + 1);
     recHitYResDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit YRes Disk1 by plaquette", 100, -200., 200.);
     //Y resolution for Disk2 by plaquette
     sprintf(histo, "RecHit_YRes_Disk2_Plaquette%d", i + 1);
     recHitYResDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit YRes Disk2 by plaquette", 100, -200., 200.);
   }
 
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitPullsBPIX");
   recHitXPullAllB = ibooker.book1D("RecHit_xres_b_All", "RecHit X Pull All Modules in Barrel", 100, -10.0, 10.0);
   recHitYPullAllB = ibooker.book1D("RecHit_yres_b_All", "RecHit Y Pull All Modules in Barrel", 100, -10.0, 10.0);
 
   for (int i = 0; i < 3; i++) {
     sprintf(histo, "RecHit_XPull_FlippedLadder_Layer%d", i + 1);
     recHitXPullFlippedLadderLayers[i] =
         ibooker.book1D(histo, "RecHit XPull Flipped Ladders by Layer", 100, -10.0, 10.0);
 
     sprintf(histo, "RecHit_XPull_UnFlippedLadder_Layer%d", i + 1);
     recHitXPullNonFlippedLadderLayers[i] =
         ibooker.book1D(histo, "RecHit XPull NonFlipped Ladders by Layer", 100, -10.0, 10.0);
   }
 
   for (int i = 0; i < 8; i++) {
     sprintf(histo, "RecHit_YPull_Layer1_Module%d", i + 1);
     recHitYPullLayer1Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer1 by module", 100, -10.0, 10.0);
 
     sprintf(histo, "RecHit_YPull_Layer2_Module%d", i + 1);
     recHitYPullLayer2Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer2 by module", 100, -10.0, 10.0);
 
     sprintf(histo, "RecHit_YPull_Layer3_Module%d", i + 1);
     recHitYPullLayer3Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer3 by module", 100, -10.0, 10.0);
   }
 
   ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitPullsFPIX");
   recHitXPullAllF = ibooker.book1D("RecHit_XPull_f_All", "RecHit X Pull All in Forward", 100, -10.0, 10.0);
 
   recHitYPullAllF = ibooker.book1D("RecHit_YPull_f_All", "RecHit Y Pull All in Forward", 100, -10.0, 10.0);
 
   for (int i = 0; i < 7; i++) {
     sprintf(histo, "RecHit_XPull_Disk1_Plaquette%d", i + 1);
     recHitXPullDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit XPull Disk1 by plaquette", 100, -10.0, 10.0);
     sprintf(histo, "RecHit_XPull_Disk2_Plaquette%d", i + 1);
     recHitXPullDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit XPull Disk2 by plaquette", 100, -10.0, 10.0);
 
     sprintf(histo, "RecHit_YPull_Disk1_Plaquette%d", i + 1);
     recHitYPullDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit YPull Disk1 by plaquette", 100, -10.0, 10.0);
 
     sprintf(histo, "RecHit_YPull_Disk2_Plaquette%d", i + 1);
     recHitYPullDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit YPull Disk2 by plaquette", 100, -10.0, 10.0);
   }
 }

◆ fillBarrel()

void SiPixelRecHitsValid::fillBarrel	(	const SiPixelRecHit &	recHit,
		const PSimHit &	simHit,
		DetId	detId,
		const PixelGeomDetUnit *	theGeomDet,
		const TrackerTopology *	tTopo
	)

private

Definition at line 370 of file SiPixelRecHitsValid.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, clustChargeLayer1Modules, clustChargeLayer2Modules, clustChargeLayer3Modules, clustXSizeLayer, clustYSizeModule, cmtomicron, dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), mps_fire::i, PixelTopology::nrows(), PV3DBase< T, PVType, FrameType >::perp(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), rpcPointValidation_cfi::recHit, recHitBunchB, recHitEventB, recHitXFullModules, recHitXHalfModules, recHitXPullAllB, recHitXPullFlippedLadderLayers, recHitXPullNonFlippedLadderLayers, recHitXResAllB, recHitXResFlippedLadderLayers, recHitXResNonFlippedLadderLayers, recHitYAllModules, recHitYPullAllB, recHitYPullLayer1Modules, recHitYPullLayer2Modules, recHitYPullLayer3Modules, recHitYResAllB, recHitYResLayer1Modules, recHitYResLayer2Modules, recHitYResLayer3Modules, postprocess-scan-build::rows, rpcPointValidation_cfi::simHit, RecoTauValidation_cfi::sizeX, RecoTauValidation_cfi::sizeY, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by analyze().

                                                                    {
   const float cmtomicron = 10000.0;
 
   int bunch = simHit.eventId().bunchCrossing();
   int event = simHit.eventId().event();
 
   recHitBunchB->Fill(bunch);
   if (bunch == 0)
     recHitEventB->Fill(event);
 
   LocalPoint lp = recHit.localPosition();
   float lp_y = lp.y();
   float lp_x = lp.x();
 
   LocalError lerr = recHit.localPositionError();
   float lerr_x = sqrt(lerr.xx());
   float lerr_y = sqrt(lerr.yy());
 
   recHitYAllModules->Fill(lp_y);
 
   float sim_x1 = simHit.entryPoint().x();
   float sim_x2 = simHit.exitPoint().x();
   float sim_xpos = 0.5 * (sim_x1 + sim_x2);
   float res_x = (lp.x() - sim_xpos) * cmtomicron;
 
   recHitXResAllB->Fill(res_x);
 
   float sim_y1 = simHit.entryPoint().y();
   float sim_y2 = simHit.exitPoint().y();
   float sim_ypos = 0.5 * (sim_y1 + sim_y2);
   float res_y = (lp.y() - sim_ypos) * cmtomicron;
 
   recHitYResAllB->Fill(res_y);
 
   float pull_x = (lp_x - sim_xpos) / lerr_x;
   float pull_y = (lp_y - sim_ypos) / lerr_y;
 
   recHitXPullAllB->Fill(pull_x);
   recHitYPullAllB->Fill(pull_y);
 
   int rows = theGeomDet->specificTopology().nrows();
 
   if (rows == 160) {
     recHitXFullModules->Fill(lp_x);
   } else if (rows == 80) {
     recHitXHalfModules->Fill(lp_x);
   }
 
   float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 0.)).perp();
   float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 1.)).perp();
 
   if (tmp2 < tmp1) {  // flipped
     for (unsigned int i = 0; i < 3; i++) {
       if (tTopo->pxbLayer(detId) == i + 1) {
         recHitXResFlippedLadderLayers[i]->Fill(res_x);
         recHitXPullFlippedLadderLayers[i]->Fill(pull_x);
       }
     }
   } else {
     for (unsigned int i = 0; i < 3; i++) {
       if (tTopo->pxbLayer(detId) == i + 1) {
         recHitXResNonFlippedLadderLayers[i]->Fill(res_x);
         recHitXPullNonFlippedLadderLayers[i]->Fill(pull_x);
       }
     }
   }
 
   //get cluster
   SiPixelRecHit::ClusterRef const& clust = recHit.cluster();
 
   // fill module dependent info
   for (unsigned int i = 0; i < 8; i++) {
     if (tTopo->pxbModule(detId) == i + 1) {
       int sizeY = (*clust).sizeY();
       clustYSizeModule[i]->Fill(sizeY);
 
       if (tTopo->pxbLayer(detId) == 1) {
         float charge = (*clust).charge();
         clustChargeLayer1Modules[i]->Fill(charge);
         recHitYResLayer1Modules[i]->Fill(res_y);
         recHitYPullLayer1Modules[i]->Fill(pull_y);
       } else if (tTopo->pxbLayer(detId) == 2) {
         float charge = (*clust).charge();
         clustChargeLayer2Modules[i]->Fill(charge);
         recHitYResLayer2Modules[i]->Fill(res_y);
         recHitYPullLayer2Modules[i]->Fill(pull_y);
       } else if (tTopo->pxbLayer(detId) == 3) {
         float charge = (*clust).charge();
         clustChargeLayer3Modules[i]->Fill(charge);
         recHitYResLayer3Modules[i]->Fill(res_y);
         recHitYPullLayer3Modules[i]->Fill(pull_y);
       }
     }
   }
   int sizeX = (*clust).sizeX();
   if (tTopo->pxbLayer(detId) == 1)
     clustXSizeLayer[0]->Fill(sizeX);
   if (tTopo->pxbLayer(detId) == 2)
     clustXSizeLayer[1]->Fill(sizeX);
   if (tTopo->pxbLayer(detId) == 3)
     clustXSizeLayer[2]->Fill(sizeX);
 }

◆ fillForward()

void SiPixelRecHitsValid::fillForward	(	const SiPixelRecHit &	recHit,
		const PSimHit &	simHit,
		DetId	detId,
		const PixelGeomDetUnit *	theGeomDet,
		const TrackerTopology *	tTopo
	)

private

Definition at line 477 of file SiPixelRecHitsValid.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, clustChargeDisk1Plaquettes, clustChargeDisk2Plaquettes, clustXSizeDisk1Plaquettes, clustXSizeDisk2Plaquettes, clustYSizeDisk1Plaquettes, clustYSizeDisk2Plaquettes, cmtomicron, dqm::impl::MonitorElement::Fill(), mps_fire::i, PixelTopology::ncolumns(), PixelTopology::nrows(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), rpcPointValidation_cfi::recHit, recHitBunchF, recHitEventF, recHitXPlaquetteSize1, recHitXPlaquetteSize2, recHitXPullAllF, recHitXPullDisk1Plaquettes, recHitXPullDisk2Plaquettes, recHitXResAllF, recHitXResDisk1Plaquettes, recHitXResDisk2Plaquettes, recHitYPlaquetteSize2, recHitYPlaquetteSize3, recHitYPlaquetteSize4, recHitYPlaquetteSize5, recHitYPullAllF, recHitYPullDisk1Plaquettes, recHitYPullDisk2Plaquettes, recHitYResDisk1Plaquettes, recHitYResDisk2Plaquettes, postprocess-scan-build::rows, rpcPointValidation_cfi::simHit, RecoTauValidation_cfi::sizeX, RecoTauValidation_cfi::sizeY, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by analyze().

                                                                     {
   int rows = theGeomDet->specificTopology().nrows();
   int cols = theGeomDet->specificTopology().ncolumns();
 
   const float cmtomicron = 10000.0;
 
   int bunch = simHit.eventId().bunchCrossing();
   int event = simHit.eventId().event();
 
   recHitBunchF->Fill(bunch);
   if (bunch == 0)
     recHitEventF->Fill(event);
 
   LocalPoint lp = recHit.localPosition();
   float lp_x = lp.x();
   float lp_y = lp.y();
 
   LocalError lerr = recHit.localPositionError();
   float lerr_x = sqrt(lerr.xx());
   float lerr_y = sqrt(lerr.yy());
 
   float sim_x1 = simHit.entryPoint().x();
   float sim_x2 = simHit.exitPoint().x();
   float sim_xpos = 0.5 * (sim_x1 + sim_x2);
 
   float sim_y1 = simHit.entryPoint().y();
   float sim_y2 = simHit.exitPoint().y();
   float sim_ypos = 0.5 * (sim_y1 + sim_y2);
 
   float pull_x = (lp_x - sim_xpos) / lerr_x;
   float pull_y = (lp_y - sim_ypos) / lerr_y;
 
   if (rows == 80) {
     recHitXPlaquetteSize1->Fill(lp_x);
   } else if (rows == 160) {
     recHitXPlaquetteSize2->Fill(lp_x);
   }
 
   if (cols == 104) {
     recHitYPlaquetteSize2->Fill(lp_y);
   } else if (cols == 156) {
     recHitYPlaquetteSize3->Fill(lp_y);
   } else if (cols == 208) {
     recHitYPlaquetteSize4->Fill(lp_y);
   } else if (cols == 260) {
     recHitYPlaquetteSize5->Fill(lp_y);
   }
 
   float res_x = (lp.x() - sim_xpos) * cmtomicron;
 
   recHitXResAllF->Fill(res_x);
   recHitXPullAllF->Fill(pull_x);
 
   float res_y = (lp.y() - sim_ypos) * cmtomicron;
 
   recHitYPullAllF->Fill(pull_y);
 
   // get cluster
   SiPixelRecHit::ClusterRef const& clust = recHit.cluster();
 
   // fill plaquette dependent info
   for (unsigned int i = 0; i < 7; i++) {
     if (tTopo->pxfModule(detId) == i + 1) {
       if (tTopo->pxfDisk(detId) == 1) {
         int sizeX = (*clust).sizeX();
         clustXSizeDisk1Plaquettes[i]->Fill(sizeX);
 
         int sizeY = (*clust).sizeY();
         clustYSizeDisk1Plaquettes[i]->Fill(sizeY);
 
         float charge = (*clust).charge();
         clustChargeDisk1Plaquettes[i]->Fill(charge);
 
         recHitXResDisk1Plaquettes[i]->Fill(res_x);
         recHitYResDisk1Plaquettes[i]->Fill(res_y);
 
         recHitXPullDisk1Plaquettes[i]->Fill(pull_x);
         recHitYPullDisk1Plaquettes[i]->Fill(pull_y);
       } else {
         int sizeX = (*clust).sizeX();
         clustXSizeDisk2Plaquettes[i]->Fill(sizeX);
 
         int sizeY = (*clust).sizeY();
         clustYSizeDisk2Plaquettes[i]->Fill(sizeY);
 
         float charge = (*clust).charge();
         clustChargeDisk2Plaquettes[i]->Fill(charge);
 
         recHitXResDisk2Plaquettes[i]->Fill(res_x);
         recHitYResDisk2Plaquettes[i]->Fill(res_y);
 
         recHitXPullDisk2Plaquettes[i]->Fill(pull_x);
         recHitYPullDisk2Plaquettes[i]->Fill(pull_y);
 
       }  // end else
     }    // end if module
     else if (tTopo->pxfPanel(detId) == 2 && (tTopo->pxfModule(detId) + 4) == i + 1) {
       if (tTopo->pxfDisk(detId) == 1) {
         int sizeX = (*clust).sizeX();
         clustXSizeDisk1Plaquettes[i]->Fill(sizeX);
 
         int sizeY = (*clust).sizeY();
         clustYSizeDisk1Plaquettes[i]->Fill(sizeY);
 
         float charge = (*clust).charge();
         clustChargeDisk1Plaquettes[i]->Fill(charge);
 
         recHitXResDisk1Plaquettes[i]->Fill(res_x);
         recHitYResDisk1Plaquettes[i]->Fill(res_y);
 
         recHitXPullDisk1Plaquettes[i]->Fill(pull_x);
         recHitYPullDisk1Plaquettes[i]->Fill(pull_y);
       } else {
         int sizeX = (*clust).sizeX();
         clustXSizeDisk2Plaquettes[i]->Fill(sizeX);
 
         int sizeY = (*clust).sizeY();
         clustYSizeDisk2Plaquettes[i]->Fill(sizeY);
 
         float charge = (*clust).charge();
         clustChargeDisk2Plaquettes[i]->Fill(charge);
 
         recHitXResDisk2Plaquettes[i]->Fill(res_x);
         recHitYResDisk2Plaquettes[i]->Fill(res_y);
 
         recHitXPullDisk2Plaquettes[i]->Fill(pull_x);
         recHitYPullDisk2Plaquettes[i]->Fill(pull_y);
 
       }  // end else
     }    // end else
   }      // end for
 }