#include <SiPixelClusterSource.h>

Inheritance diagram for SiPixelClusterSource:

Public Types
typedef edmNew::DetSet< SiPixelCluster >::const_iterator	ClusterIterator

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

virtual void	bookMEs (DQMStore::IBooker &, const edm::EventSetup &iSetup)

virtual void	buildStructure (edm::EventSetup const &)

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

	SiPixelClusterSource (const edm::ParameterSet &conf)

	~SiPixelClusterSource () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Public Attributes
std::string	topFolderName_

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

void	getrococcupancye (DetId detId, const edmNew::DetSetVector< SiPixelCluster > &clustColl, const TrackerTopology const pTT, edm::ESHandle< TrackerGeometry > pDD, MonitorElement meinput)

Private Attributes
int	bigEventSize

MonitorElement *	bigFpixClusterEventRate

bool	bladeOn

edm::ParameterSet	conf_

edm::InputTag	digisrc_

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

bool	diskOn

int	eventNo

bool	firstRun

bool	isPIB

bool	isUpgrade

bool	ladOn

bool	layOn

int	lumSec

std::vector< MonitorElement * >	meClPosDiskmz

std::vector< MonitorElement * >	meClPosDiskpz

std::vector< MonitorElement * >	meClPosLayer

MonitorElement *	meClusBarrelProf

MonitorElement *	meClusFpixMProf

MonitorElement *	meClusFpixPProf

std::vector< MonitorElement * >	meZeroRocBPIX

MonitorElement *	meZeroRocFPIX

bool	modOn

int	nBigEvents

int	nLumiSecs

int	noOfDisks

int	noOfLayers

bool	phiOn

bool	reducedSet

bool	ringOn

bool	saveFile

bool	slowDown

bool	smileyOn

edm::InputTag	src_

edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	srcToken_

std::map< uint32_t, SiPixelClusterModule * >	thePixelStructure

bool	twoDimOn

Detailed Description

Definition at line 59 of file SiPixelClusterSource.h.

Member Typedef Documentation

typedef edmNew::DetSet<SiPixelCluster>::const_iterator SiPixelClusterSource::ClusterIterator

Definition at line 64 of file SiPixelClusterSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 of file SiPixelClusterSource.cc.

References conf_, digisrcToken_, firstRun, edm::ParameterSet::getParameter(), srcToken_, AlCaHLTBitMon_QueryRunRegistry::string, and topFolderName_.

     : conf_(iConfig),
       src_(conf_.getParameter<edm::InputTag>("src")),
       digisrc_(conf_.getParameter<edm::InputTag>("digisrc")),
       saveFile(conf_.getUntrackedParameter<bool>("saveFile", false)),
       isPIB(conf_.getUntrackedParameter<bool>("isPIB", false)),
       slowDown(conf_.getUntrackedParameter<bool>("slowDown", false)),
       modOn(conf_.getUntrackedParameter<bool>("modOn", true)),
       twoDimOn(conf_.getUntrackedParameter<bool>("twoDimOn", true)),
       reducedSet(conf_.getUntrackedParameter<bool>("reducedSet", false)),
       ladOn(conf_.getUntrackedParameter<bool>("ladOn", false)),
       layOn(conf_.getUntrackedParameter<bool>("layOn", false)),
       phiOn(conf_.getUntrackedParameter<bool>("phiOn", false)),
       ringOn(conf_.getUntrackedParameter<bool>("ringOn", false)),
       bladeOn(conf_.getUntrackedParameter<bool>("bladeOn", false)),
       diskOn(conf_.getUntrackedParameter<bool>("diskOn", false)),
       smileyOn(conf_.getUntrackedParameter<bool>("smileyOn", false)),
       bigEventSize(conf_.getUntrackedParameter<int>("bigEventSize", 100)),
       isUpgrade(conf_.getUntrackedParameter<bool>("isUpgrade", false)),
       noOfLayers(0),
       noOfDisks(0) {
   LogInfo("PixelDQM") << "SiPixelClusterSource::SiPixelClusterSource: Got DQM BackEnd interface" << endl;
 
   // set Token(-s)
   srcToken_ = consumes<edmNew::DetSetVector<SiPixelCluster>>(conf_.getParameter<edm::InputTag>("src"));
   digisrcToken_ = consumes<edm::DetSetVector<PixelDigi>>(conf_.getParameter<edm::InputTag>("digisrc"));
   firstRun = true;
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
 }

SiPixelClusterSource::~SiPixelClusterSource ( )

override

Definition at line 78 of file SiPixelClusterSource.cc.

References thePixelStructure.

                                             {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   LogInfo("PixelDQM") << "SiPixelClusterSource::~SiPixelClusterSource: Destructor" << endl;
 
   std::map<uint32_t, SiPixelClusterModule *>::iterator struct_iter;
   for (struct_iter = thePixelStructure.begin(); struct_iter != thePixelStructure.end(); struct_iter++) {
     delete struct_iter->second;
     struct_iter->second = nullptr;
   }
 }

Member Function Documentation

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

override

Definition at line 213 of file SiPixelClusterSource.cc.

References bigEventSize, bigFpixClusterEventRate, bladeOn, TrackerGeometry::dets(), digisrcToken_, diskOn, eventNo, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), getrococcupancy(), getrococcupancye(), mps_fire::i, input, createfilelist::int, isUpgrade, ladOn, layOn, edm::EventBase::luminosityBlock(), meClPosDiskmz, meClPosDiskpz, meClPosLayer, meClusBarrelProf, meClusFpixMProf, meClusFpixPProf, meZeroRocBPIX, meZeroRocFPIX, modOn, noOfDisks, noOfLayers, edm::EventBase::orbitNumber(), phiOn, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::Handle< T >::product(), edm::ESHandle< T >::product(), reducedSet, ringOn, slowDown, smileyOn, srcToken_, DetId::subdetId(), thePixelStructure, trackingTruthProducer_cfi::tracker, and twoDimOn.

                                                                                       {
   eventNo++;
 
   if (meClPosLayer.at(0) && meClPosLayer.at(0)->getEntries() > 150000) {
     for (int i = 0; i < noOfLayers; i++) {
       meClPosLayer.at(i)->Reset();
     }
     for (int i = 0; i < noOfDisks; i++) {
       meClPosDiskpz.at(i)->Reset();
       meClPosDiskmz.at(i)->Reset();
     }
   }
 
   // get input data
   edm::Handle<edmNew::DetSetVector<SiPixelCluster>> input;
   iEvent.getByToken(srcToken_, input);
   auto const &clustColl = *(input.product());
 
   edm::ESHandle<TrackerGeometry> pDD;
   iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
   const TrackerGeometry *tracker = &(*pDD);
 
   edm::Handle<edm::DetSetVector<PixelDigi>> digiinput;
   iEvent.getByToken(digisrcToken_, digiinput);
   const edm::DetSetVector<PixelDigi> diginp = *(digiinput.product());
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *pTT = tTopoHandle.product();
 
   int lumiSection = (int)iEvent.luminosityBlock();
   int nEventFpixClusters = 0;
 
   int nEventsBarrel = 0;
   int nEventsFPIXm = 0;
   int nEventsFPIXp = 0;
 
   std::map<uint32_t, SiPixelClusterModule *>::iterator struct_iter;
   for (struct_iter = thePixelStructure.begin(); struct_iter != thePixelStructure.end(); struct_iter++) {
     int numberOfFpixClusters = (*struct_iter)
                                    .second->fill(*input,
                                                  tracker,
                                                  &nEventsBarrel,
                                                  &nEventsFPIXp,
                                                  &nEventsFPIXm,
                                                  meClPosLayer,
                                                  meClPosDiskpz,
                                                  meClPosDiskmz,
                                                  modOn,
                                                  ladOn,
                                                  layOn,
                                                  phiOn,
                                                  bladeOn,
                                                  diskOn,
                                                  ringOn,
                                                  twoDimOn,
                                                  reducedSet,
                                                  smileyOn,
                                                  isUpgrade);
     nEventFpixClusters = nEventFpixClusters + numberOfFpixClusters;
   }
 
   if (nEventFpixClusters > bigEventSize) {
     if (bigFpixClusterEventRate) {
       bigFpixClusterEventRate->Fill(lumiSection, 1. / 23.);
     }
   }
 
   float trendVar = iEvent.orbitNumber() / 262144.0;  // lumisection : seconds - matches strip trend plot
 
   meClusBarrelProf->Fill(trendVar, nEventsBarrel);
   meClusFpixPProf->Fill(trendVar, nEventsFPIXp);
   meClusFpixMProf->Fill(trendVar, nEventsFPIXm);
 
   // std::cout<<"nEventFpixClusters: "<<nEventFpixClusters<<" , nLumiSecs:
   // "<<nLumiSecs<<" , nBigEvents: "<<nBigEvents<<std::endl;
 
   for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) {
     DetId detId = (*it)->geographicalId();
 
     // fill barrel
     if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
       getrococcupancy(detId, diginp, pTT, meZeroRocBPIX);
     }
 
     // fill endcap
     if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
       getrococcupancye(detId, clustColl, pTT, pDD, meZeroRocFPIX);
     }
   }
 
   // slow down...
   if (slowDown)
     usleep(10000);
 }

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

override

Definition at line 109 of file SiPixelClusterSource.cc.

References DQMStore::IBooker::book2D(), bookMEs(), DQMStore::IBooker::bookProfile(), digisrc_, MonitorElement::getTH1(), mps_fire::i, edm::InputTag::label(), meClPosDiskmz, meClPosDiskpz, meClPosLayer, meClusBarrelProf, meClusFpixMProf, meClusFpixPProf, meZeroRocBPIX, meZeroRocFPIX, noOfDisks, noOfLayers, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::ymin.

                                                                                                                {
   bookMEs(iBooker, iSetup);
   // Book occupancy maps in global coordinates for all clusters:
   iBooker.setCurrentFolder(topFolderName_ + "/Clusters/OffTrack");
   // bpix
   std::stringstream ss1, ss2;
   for (int i = 1; i <= noOfLayers; i++) {
     ss1.str(std::string());
     ss1 << "position_siPixelClusters_Layer_" << i;
     ss2.str(std::string());
     ss2 << "Clusters Layer" << i << ";Global Z (cm);Global #phi";
     meClPosLayer.push_back(iBooker.book2D(ss1.str(), ss2.str(), 200, -30., 30., 128, -3.2, 3.2));
   }
   for (int i = 1; i <= noOfDisks; i++) {
     ss1.str(std::string());
     ss1 << "position_siPixelClusters_pz_Disk_" << i;
     ss2.str(std::string());
     ss2 << "Clusters +Z Disk" << i << ";Global X (cm);Global Y (cm)";
     meClPosDiskpz.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
     ss1.str(std::string());
     ss1 << "position_siPixelClusters_mz_Disk_" << i;
     ss2.str(std::string());
     ss2 << "Clusters -Z Disk" << i << ";Global X (cm);Global Y (cm)";
     meClPosDiskmz.push_back(iBooker.book2D(ss1.str(), ss2.str(), 80, -20., 20., 80, -20., 20.));
   }
 
   // Book trend cluster plots for barrel and endcap. Lumisections for offline
   // and second for online - taken from strips
   iBooker.setCurrentFolder(topFolderName_ + "/Barrel");
   ss1.str(std::string());
   ss1 << "totalNumberOfClustersProfile_siPixelClusters_Barrel";
   ss2.str(std::string());
   ss2 << "Total number of barrel clusters profile;Lumisection;";
   meClusBarrelProf = iBooker.bookProfile(ss1.str(), ss2.str(), 2400, 0., 150, 0, 0, "");
   meClusBarrelProf->getTH1()->SetCanExtend(TH1::kAllAxes);
 
   iBooker.setCurrentFolder(topFolderName_ + "/Endcap");
 
   ss1.str(std::string());
   ss1 << "totalNumberOfClustersProfile_siPixelClusters_FPIX+";
   ss2.str(std::string());
   ss2 << "Total number of FPIX+ clusters profile;Lumisection;";
   meClusFpixPProf = iBooker.bookProfile(ss1.str(), ss2.str(), 2400, 0., 150, 0, 0, "");
   meClusFpixPProf->getTH1()->SetCanExtend(TH1::kAllAxes);
 
   ss1.str(std::string());
   ss1 << "totalNumberOfClustersProfile_siPixelClusters_FPIX-";
   ss2.str(std::string());
   ss2 << "Total number of FPIX- clusters profile;Lumisection;";
   meClusFpixMProf = iBooker.bookProfile(ss1.str(), ss2.str(), 2400, 0., 150, 0, 0, "");
   meClusFpixMProf->getTH1()->SetCanExtend(TH1::kAllAxes);
 
   iBooker.setCurrentFolder(topFolderName_ + "/Barrel");
   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;
     }
     if (i == 4) {
       ybins = 130;
       ymin = -32.5;
       ymax = 32.5;
     }
     ss1.str(std::string());
     ss1 << "pix_bar Occ_roc_online_" + digisrc_.label() + "_layer_" << i;
     ss2.str(std::string());
     ss2 << "Pixel Barrel Occupancy, ROC level (Online): Layer " << i;
     meZeroRocBPIX.push_back(iBooker.book2D(ss1.str(), ss2.str(), 72, -4.5, 4.5, ybins, ymin, ymax));
     meZeroRocBPIX.at(i - 1)->setAxisTitle("ROC / Module", 1);
     meZeroRocBPIX.at(i - 1)->setAxisTitle("ROC / Ladder", 2);
   }
 
   iBooker.setCurrentFolder(topFolderName_ + "/Endcap");
   meZeroRocFPIX = iBooker.book2D(
       "ROC_endcap_occupancy", "Pixel Endcap Occupancy, ROC level (Online)", 72, -4.5, 4.5, 288, -12.5, 12.5);
   meZeroRocFPIX->setBinLabel(1, "Disk-2 Pnl2", 1);
   meZeroRocFPIX->setBinLabel(9, "Disk-2 Pnl1", 1);
   meZeroRocFPIX->setBinLabel(19, "Disk-1 Pnl2", 1);
   meZeroRocFPIX->setBinLabel(27, "Disk-1 Pnl1", 1);
   meZeroRocFPIX->setBinLabel(41, "Disk+1 Pnl1", 1);
   meZeroRocFPIX->setBinLabel(49, "Disk+1 Pnl2", 1);
   meZeroRocFPIX->setBinLabel(59, "Disk+2 Pnl1", 1);
   meZeroRocFPIX->setBinLabel(67, "Disk+2 Pnl2", 1);
   meZeroRocFPIX->setAxisTitle("Blades in Inner (>0) / Outer(<) Halves", 2);
   meZeroRocFPIX->setAxisTitle("ROC occupancy", 3);
 }

void SiPixelClusterSource::bookMEs	(	DQMStore::IBooker &	iBooker,
		const edm::EventSetup &	iSetup
	)

virtual

Create folder tree and book histograms

Definition at line 384 of file SiPixelClusterSource.cc.

References bigEventSize, bigFpixClusterEventRate, bladeOn, DQMStore::IBooker::book1D(), conf_, diskOn, Exception, isPIB, isUpgrade, ladOn, layOn, LogDebug, modOn, phiOn, reducedSet, ringOn, DQMStore::IBooker::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), smileyOn, thePixelStructure, runGCPTkAlMap::title, topFolderName_, and twoDimOn.

Referenced by bookHistograms().

                                                                                         {
   // Get DQM interface
   iBooker.setCurrentFolder(topFolderName_);
   char title[256];
   snprintf(title,
            256,
            "Rate of events with >%i FPIX clusters;LumiSection;Rate of large "
            "FPIX events per LS [Hz]",
            bigEventSize);
   bigFpixClusterEventRate = iBooker.book1D("bigFpixClusterEventRate", title, 5000, 0., 5000.);
 
   std::map<uint32_t, SiPixelClusterModule *>::iterator struct_iter;
 
   SiPixelFolderOrganizer theSiPixelFolder(false);
 
   for (struct_iter = thePixelStructure.begin(); struct_iter != thePixelStructure.end(); struct_iter++) {
     if (modOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 0, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 0, twoDimOn, reducedSet, isUpgrade);
       } else {
         if (!isPIB)
           throw cms::Exception("LogicError") << "[SiPixelClusterSource::bookMEs] Creation of DQM folder "
                                                 "failed";
       }
     }
     if (ladOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 1, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 1, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH LADDER-FOLDER\n";
       }
     }
     if (layOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 2, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 2, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH LAYER-FOLDER\n";
       }
     }
     if (phiOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 3, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 3, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH PHI-FOLDER\n";
       }
     }
     if (bladeOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 4, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 4, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH BLADE-FOLDER\n";
       }
     }
     if (diskOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 5, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 5, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH DISK-FOLDER\n";
       }
     }
     if (ringOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 6, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 6, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH RING-FOLDER\n";
       }
     }
     if (smileyOn) {
       if (theSiPixelFolder.setModuleFolder(iBooker, (*struct_iter).first, 7, isUpgrade)) {
         (*struct_iter).second->book(conf_, iSetup, iBooker, 7, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH BARREL-FOLDER\n";
       }
     }
   }
 }

void SiPixelClusterSource::buildStructure ( edm::EventSetup const & )

virtual

Definition at line 312 of file SiPixelClusterSource.cc.

References PixelEndcapName::bladeName(), TrackerGeometry::dets(), TrackerGeometry::detTypes(), PixelEndcapName::diskName(), edm::EventSetup::get(), PixelEndcapName::halfCylinder(), TrackerGeometry::idToDetUnit(), isPIB, isUpgrade, PixelBarrelName::layerName(), LogDebug, RecoTauDiscriminantConfiguration::mask, hgcalPlots::ncols, noOfDisks, noOfLayers, PixelEndcapName::pannelName(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapName::plaquetteName(), edm::ESHandle< T >::product(), DetId::rawId(), PixelGeomDetUnit::specificTopology(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), and thePixelStructure.

Referenced by dqmBeginRun().

                                                                      {
   LogInfo("PixelDQM") << " SiPixelClusterSource::buildStructure";
   edm::ESHandle<TrackerGeometry> pDD;
   iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *pTT = tTopoHandle.product();
 
   LogVerbatim("PixelDQM") << " *** Geometry node for TrackerGeom is  " << &(*pDD) << std::endl;
   LogVerbatim("PixelDQM") << " *** I have " << pDD->dets().size() << " detectors" << std::endl;
   LogVerbatim("PixelDQM") << " *** I have " << pDD->detTypes().size() << " types" << std::endl;
 
   for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) {
     if (dynamic_cast<PixelGeomDetUnit const *>((*it)) != nullptr) {
       DetId detId = (*it)->geographicalId();
       const GeomDetUnit *geoUnit = pDD->idToDetUnit(detId);
       const PixelGeomDetUnit *pixDet = dynamic_cast<const PixelGeomDetUnit *>(geoUnit);
       int nrows = (pixDet->specificTopology()).nrows();
       int ncols = (pixDet->specificTopology()).ncolumns();
 
       if ((detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) ||
           (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap))) {
         uint32_t id = detId();
         if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
           if (isPIB)
             continue;
           LogDebug("PixelDQM") << " ---> Adding Barrel Module " << detId.rawId() << endl;
           int layer = PixelBarrelName(DetId(id), pTT, isUpgrade).layerName();
           if (layer > noOfLayers)
             noOfLayers = layer;
           SiPixelClusterModule *theModule = new SiPixelClusterModule(id, ncols, nrows);
           thePixelStructure.insert(pair<uint32_t, SiPixelClusterModule *>(id, theModule));
         } else if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
           LogDebug("PixelDQM") << " ---> Adding Endcap Module " << detId.rawId() << endl;
           PixelEndcapName::HalfCylinder side = PixelEndcapName(DetId(id), pTT, isUpgrade).halfCylinder();
           int disk = PixelEndcapName(DetId(id), pTT, isUpgrade).diskName();
           if (disk > noOfDisks)
             noOfDisks = disk;
           int blade = PixelEndcapName(DetId(id), pTT, isUpgrade).bladeName();
           int panel = PixelEndcapName(DetId(id), pTT, isUpgrade).pannelName();
           int module = PixelEndcapName(DetId(id), pTT, isUpgrade).plaquetteName();
           char sside[80];
           sprintf(sside, "HalfCylinder_%i", side);
           char sdisk[80];
           sprintf(sdisk, "Disk_%i", disk);
           char sblade[80];
           sprintf(sblade, "Blade_%02i", blade);
           char spanel[80];
           sprintf(spanel, "Panel_%i", panel);
           char smodule[80];
           sprintf(smodule, "Module_%i", module);
           std::string side_str = sside;
           std::string disk_str = sdisk;
           bool mask = side_str.find("HalfCylinder_1") != string::npos ||
                       side_str.find("HalfCylinder_2") != string::npos ||
                       side_str.find("HalfCylinder_4") != string::npos || disk_str.find("Disk_2") != string::npos;
           // clutch to take all of FPIX, but no BPIX:
           mask = false;
           if (isPIB && mask)
             continue;
           SiPixelClusterModule *theModule = new SiPixelClusterModule(id, ncols, nrows);
           thePixelStructure.insert(pair<uint32_t, SiPixelClusterModule *>(id, theModule));
         }
       }
     }
   }
   LogInfo("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
 }

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

override

Definition at line 90 of file SiPixelClusterSource.cc.

References bladeOn, buildStructure(), diskOn, eventNo, firstRun, ladOn, layOn, lumSec, modOn, nBigEvents, nLumiSecs, phiOn, ringOn, smileyOn, and twoDimOn.

                                                                                    {
   LogInfo("PixelDQM") << " SiPixelClusterSource::beginJob - Initialisation ... " << std::endl;
   LogInfo("PixelDQM") << "Mod/Lad/Lay/Phi " << modOn << "/" << ladOn << "/" << layOn << "/" << phiOn << std::endl;
   LogInfo("PixelDQM") << "Blade/Disk/Ring" << bladeOn << "/" << diskOn << "/" << ringOn << std::endl;
   LogInfo("PixelDQM") << "2DIM IS " << twoDimOn << "\n";
   LogInfo("PixelDQM") << "Smiley (Cluster sizeY vs. Cluster eta) is " << smileyOn << "\n";
 
   if (firstRun) {
     eventNo = 0;
     lumSec = 0;
     nLumiSecs = 0;
     nBigEvents = 0;
     // Build map
     buildStructure(iSetup);
 
     firstRun = false;
   }
 }

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

private

Definition at line 462 of file SiPixelClusterSource.cc.

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

Referenced by analyze().

                                                                                        {
   edm::DetSetVector<PixelDigi>::const_iterator ipxsearch = diginp.find(detId);
   if (ipxsearch != diginp.end()) {
     // Look at digis now
     edm::DetSet<PixelDigi>::const_iterator pxdi;
     for (pxdi = ipxsearch->begin(); pxdi != ipxsearch->end(); pxdi++) {
       bool isHalfModule = PixelBarrelName(DetId(detId), tTopo, isUpgrade).isHalfModule();
       int DBlayer = PixelBarrelName(DetId(detId), tTopo, isUpgrade).layerName();
       int DBmodule = PixelBarrelName(DetId(detId), tTopo, isUpgrade).moduleName();
       int DBladder = PixelBarrelName(DetId(detId), tTopo, isUpgrade).ladderName();
       int DBshell = PixelBarrelName(DetId(detId), tTopo, isUpgrade).shell();
 
       // add sign to the modules
       if (DBshell == 1 || DBshell == 2) {
         DBmodule = -DBmodule;
       }
       if (DBshell == 1 || DBshell == 3) {
         DBladder = -DBladder;
       }
 
       int col = pxdi->column();
       int row = pxdi->row();
 
       float modsign = (float)DBmodule / (abs((float)DBmodule));
       float ladsign = (float)DBladder / (abs((float)DBladder));
       float rocx = ((float)col / (52. * 8.)) * modsign + ((float)DBmodule - (modsign)*0.5);
       float rocy = ((float)row / (80. * 2.)) * ladsign + ((float)DBladder - (ladsign)*0.5);
 
       // do the flip where need
       bool flip = false;
       if ((DBladder % 2 == 0) && (!isHalfModule)) {
         flip = true;
       }
       if ((flip) && (DBladder > 0)) {
         if ((((float)DBladder - (ladsign)*0.5) <= rocy) && (rocy < (float)DBladder)) {
           rocy = rocy + ladsign * 0.5;
         } else if ((((float)DBladder) <= rocy) && (rocy < ((float)DBladder + (ladsign)*0.5))) {
           rocy = rocy - ladsign * 0.5;
         }
       }
 
       // tweak border effect for negative modules/ladders
       if (modsign < 0) {
         rocx = rocx - 0.0001;
       }
       if (ladsign < 0) {
         rocy = rocy - 0.0001;
       } else {
         rocy = rocy + 0.0001;
       }
       if (abs(DBladder) == 1) {
         rocy = rocy + ladsign * 0.5;
       }  // take care of the half module
       meinput[DBlayer - 1]->Fill(rocx, rocy);
     }  // end of looping over pxdi
   }
 }

void SiPixelClusterSource::getrococcupancye	(	DetId	detId,
		const edmNew::DetSetVector< SiPixelCluster > &	clustColl,
		const TrackerTopology *const	pTT,
		edm::ESHandle< TrackerGeometry >	pDD,
		MonitorElement *	meinput
	)

private

Definition at line 523 of file SiPixelClusterSource.cc.

References funct::abs(), edmNew::DetSetVector< T >::begin(), PixelEndcapName::bladeName(), DEFINE_FWK_MODULE, PixelEndcapName::diskName(), edmNew::DetSetVector< T >::end(), edmNew::DetSetVector< T >::find(), MonitorElement::getBinContent(), PixelEndcapName::halfCylinder(), mps_fire::i, TrackerGeometry::idToDetUnit(), isUpgrade, Topology::localPosition(), PixelEndcapName::pannelName(), PixelEndcapName::plaquetteName(), MonitorElement::setBinContent(), PixelGeomDetUnit::specificTopology(), GeomDet::surface(), Surface::toGlobal(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                                      {
   edmNew::DetSetVector<SiPixelCluster>::const_iterator ipxsearch = clustColl.find(detId);
   if (ipxsearch != clustColl.end()) {
     // Look at clusters now
     edmNew::DetSet<SiPixelCluster>::const_iterator pxclust;
     for (pxclust = ipxsearch->begin(); pxclust != ipxsearch->end(); pxclust++) {
       const GeomDetUnit *geoUnit = pDD->idToDetUnit(detId);
       const PixelGeomDetUnit *pixDet = dynamic_cast<const PixelGeomDetUnit *>(geoUnit);
       const PixelTopology *topol = &(pixDet->specificTopology());
       LocalPoint clustlp = topol->localPosition(MeasurementPoint(pxclust->x(), pxclust->y()));
       GlobalPoint clustgp = geoUnit->surface().toGlobal(clustlp);
 
       float xclust = pxclust->x();
       float yclust = pxclust->y();
       float z = clustgp.z();
 
       int pxfside = PixelEndcapName(detId, pTT, isUpgrade).halfCylinder();
       int pxfpanel = PixelEndcapName(detId, pTT, isUpgrade).pannelName();
       int pxfmodule = PixelEndcapName(detId, pTT, isUpgrade).plaquetteName();
       int pxfdisk = PixelEndcapName(detId, pTT, isUpgrade).diskName();
       int pxfblade = PixelEndcapName(detId, pTT, isUpgrade).bladeName();
 
       if ((pxfside == 1) || (pxfside == 3)) {
         pxfblade = -1. * pxfblade;
       }
 
       if (z < 0.) {
         pxfdisk = -1. * pxfdisk;
       }
 
       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);
     }
   }
 }