#include <SiPixelClusterSource.h>

Inheritance diagram for SiPixelClusterSource:

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

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 60 of file SiPixelClusterSource.h.

Member Typedef Documentation

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

Definition at line 65 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 80 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 189 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 &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 116 of file SiPixelClusterSource.cc.

References DQMStore::IBooker::book2D(), bookMEs(), DQMStore::IBooker::bookProfile(), digisrc_, MonitorElement::getTH1(), mps_fire::i, diffTwoXMLs::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 351 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, fftjetcommon_cfi::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 284 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, trackingPlots::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();
         SiPixelClusterModule* theModule = new SiPixelClusterModule(id, ncols, nrows);
         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;
       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;
       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 95 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 429 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 475 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);
       
     }
   }
   
 }