#include <SiPixelClusterSource.h>

Inheritance diagram for SiPixelClusterSource:

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

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

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 DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Public Attributes
std::string	topFolderName_

Private Member Functions
void	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

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomTokenBeginRun_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	trackerTopoTokenBeginRun_

bool	twoDimOn

Additional Inherited Members
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 35 of file SiPixelClusterSource.h.

Member Typedef Documentation

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

Definition at line 40 of file SiPixelClusterSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 45 of file SiPixelClusterSource.cc.

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

     : 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"));
 
   trackerTopoToken_ = esConsumes<TrackerTopology, TrackerTopologyRcd>();
   trackerGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
   trackerTopoTokenBeginRun_ = esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>();
   trackerGeomTokenBeginRun_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>();
 
   firstRun = true;
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
 }

SiPixelClusterSource::~SiPixelClusterSource ( )

override

Definition at line 81 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
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 216 of file SiPixelClusterSource.cc.

References bigEventSize, bigFpixClusterEventRate, bladeOn, digisrcToken_, diskOn, eventNo, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventSetup::getHandle(), getrococcupancy(), getrococcupancye(), mps_fire::i, input, 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< class >::product(), reducedSet, ringOn, slowDown, smileyOn, srcToken_, DetId::subdetId(), thePixelStructure, trackerGeomToken_, trackerTopoToken_, 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.getHandle(trackerGeomToken_);
   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.getHandle(trackerTopoToken_);
   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,
                                                  pTT,
                                                  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 &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 112 of file SiPixelClusterSource.cc.

References dqm::implementation::IBooker::book2D(), bookMEs(), dqm::implementation::IBooker::bookProfile(), digisrc_, dqm::impl::MonitorElement::getTH1(), mps_fire::i, edm::InputTag::label(), meClPosDiskmz, meClPosDiskpz, meClPosLayer, meClusBarrelProf, meClusFpixMProf, meClusFpixPProf, meZeroRocBPIX, meZeroRocFPIX, noOfDisks, noOfLayers, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, SiStrip_OfflineMonitoring_cff::ymax, and SiStrip_OfflineMonitoring_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, dqm::implementation::IBooker::book1D(), conf_, diskOn, Exception, edm::EventSetup::getHandle(), isPIB, isUpgrade, ladOn, layOn, LogDebug, modOn, phiOn, edm::ESHandle< class >::product(), reducedSet, ringOn, dqm::implementation::NavigatorBase::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), smileyOn, thePixelStructure, runGCPTkAlMap::title, topFolderName_, trackerTopoTokenBeginRun_, 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);
 
   edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(trackerTopoTokenBeginRun_);
   const TrackerTopology *pTT = tTopoHandle.product();
 
   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_, pTT, 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_, pTT, 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_, pTT, 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_, pTT, 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_, pTT, 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_, pTT, 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_, pTT, 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_, pTT, iBooker, 7, twoDimOn, reducedSet, isUpgrade);
       } else {
         LogDebug("PixelDQM") << "PROBLEM WITH BARREL-FOLDER\n";
       }
     }
   }
 }

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

virtual

Definition at line 314 of file SiPixelClusterSource.cc.

References PixelEndcapName::bladeName(), PixelEndcapName::diskName(), edm::EventSetup::getHandle(), PixelEndcapName::halfCylinder(), isPIB, isUpgrade, phase1PixelTopology::layer, PixelBarrelName::layerName(), LogDebug, callgraph::module, noOfDisks, noOfLayers, PixelEndcapName::pannelName(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapName::plaquetteName(), edm::ESHandle< class >::product(), DetId::rawId(), PixelGeomDetUnit::specificTopology(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), thePixelStructure, trackerGeomTokenBeginRun_, and trackerTopoTokenBeginRun_.

Referenced by dqmBeginRun().

                                                                      {
   LogInfo("PixelDQM") << " SiPixelClusterSource::buildStructure";
   edm::ESHandle<TrackerGeometry> pDD = iSetup.getHandle(trackerGeomTokenBeginRun_);
 
   edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(trackerTopoTokenBeginRun_);
   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 &	r,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 93 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 465 of file SiPixelClusterSource.cc.

References funct::abs(), edm::DetSetVector< T >::begin(), cuy::col, edm::DetSetVector< T >::end(), edm::DetSetVector< T >::find(), 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 526 of file SiPixelClusterSource.cc.

References funct::abs(), edmNew::DetSetVector< T >::begin(), PixelEndcapName::bladeName(), PixelEndcapName::diskName(), edmNew::DetSetVector< T >::end(), edmNew::DetSetVector< T >::find(), dqm::impl::MonitorElement::getBinContent(), PixelEndcapName::halfCylinder(), mps_fire::i, isUpgrade, dqmiolumiharvest::j, Topology::localPosition(), PixelEndcapName::pannelName(), PixelEndcapName::plaquetteName(), dqm::impl::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);
     }
   }
 }