#include <SiPixelRecHitSource.h>

Inheritance diagram for SiPixelRecHitSource:

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

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, const edm::EventSetup &) 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

	SiPixelRecHitSource (const edm::ParameterSet &conf)

	~SiPixelRecHitSource () 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	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

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 Attributes
bool	bladeOn

edm::ParameterSet	conf_

bool	diskOn

int	eventNo

bool	firstRun

bool	isPIB

bool	isUpgrade

bool	ladOn

bool	layOn

bool	modOn

bool	phiOn

std::map< uint32_t, int >	rechit_count

bool	reducedSet

bool	ringOn

bool	saveFile

bool	slowDown

edm::EDGetTokenT< SiPixelRecHitCollection >	src_

std::map< uint32_t, SiPixelRecHitModule * >	thePixelStructure

bool	twoDimOn

Additional Inherited Members
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

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

Description: header file for Pixel Monitor Rec Hits

Usage: see description

Definition at line 54 of file SiPixelRecHitSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 51 of file SiPixelRecHitSource.cc.

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

                                                                        :
   conf_(iConfig),
   src_( consumes<SiPixelRecHitCollection>( conf_.getParameter<edm::InputTag>( "src" ))),
   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) ), 
   isUpgrade( conf_.getUntrackedParameter<bool>("isUpgrade",false) )
 {
   firstRun = true;
   LogInfo ("PixelDQM") << "SiPixelRecHitSource::SiPixelRecHitSource: Got DQM BackEnd interface"<<endl;
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
 }

SiPixelRecHitSource::~SiPixelRecHitSource ( )

override

Definition at line 74 of file SiPixelRecHitSource.cc.

References thePixelStructure.

 {
    // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   LogInfo ("PixelDQM") << "SiPixelRecHitSource::~SiPixelRecHitSource: Destructor"<<endl;
   std::map<uint32_t,SiPixelRecHitModule*>::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 SiPixelRecHitSource::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 112 of file SiPixelRecHitSource.cc.

References edmNew::DetSet< T >::begin(), bladeOn, diskOn, edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), eventNo, edmNew::DetSetVector< T >::find(), edm::Event::getByToken(), ladOn, layOn, match(), modOn, phiOn, rechit_count, reducedSet, ringOn, RecoTauValidation_cfi::sizeX, RecoTauValidation_cfi::sizeY, slowDown, mathSSE::sqrt(), src_, thePixelStructure, twoDimOn, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

 {
   eventNo++;
   //cout << eventNo << endl;
   // get input data
   edm::Handle<SiPixelRecHitCollection>  recHitColl;
   iEvent.getByToken( src_, recHitColl );
 
   std::map<uint32_t,SiPixelRecHitModule*>::iterator struct_iter;
   for (struct_iter = thePixelStructure.begin() ; struct_iter != thePixelStructure.end() ; struct_iter++) {
     uint32_t TheID = (*struct_iter).first;
 
     SiPixelRecHitCollection::const_iterator match = recHitColl->find(TheID);
     
       // if( pixelrechitRangeIteratorBegin == pixelrechitRangeIteratorEnd) {cout << "oops" << endl;}
       float rechit_x = 0;
       float rechit_y = 0;
       int rechit_count = 0;
 
       if (match != recHitColl->end()) {
     SiPixelRecHitCollection::DetSet pixelrechitRange = *match;
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
     SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
     SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
 
        for ( ; pixeliter != pixelrechitRangeIteratorEnd; pixeliter++) 
      {
       
 
       rechit_count++;
       //cout << TheID << endl;
       SiPixelRecHit::ClusterRef const& clust = pixeliter->cluster();
       int sizeX = (*clust).sizeX();
       //cout << sizeX << endl;
       int sizeY = (*clust).sizeY();
       //cout << sizeY << endl;
       LocalPoint lp = pixeliter->localPosition();
       rechit_x = lp.x();
       rechit_y = lp.y();
       
       LocalError lerr = pixeliter->localPositionError();
       float lerr_x = sqrt(lerr.xx());
       float lerr_y = sqrt(lerr.yy());
       //std::cout << "errors " << lerr_x << " " << lerr_y << std::endl;
       (*struct_iter).second->fill(rechit_x, rechit_y, sizeX, sizeY, lerr_x, lerr_y, 
                                   modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn, 
                       twoDimOn, reducedSet);
     
     }
       }
       if(rechit_count > 0) (*struct_iter).second->nfill(rechit_count, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
     
   }
 
   // slow down...
   if(slowDown) usleep(10000);
   
 }

void SiPixelRecHitSource::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 105 of file SiPixelRecHitSource.cc.

References bookMEs().

                                                                                                               {
   bookMEs(iBooker, iSetup);
 }

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

virtual

Create folder tree and book histograms

Definition at line 243 of file SiPixelRecHitSource.cc.

References bladeOn, conf_, DEFINE_FWK_MODULE, diskOn, Exception, isPIB, isUpgrade, ladOn, layOn, LogDebug, modOn, phiOn, reducedSet, ringOn, SiPixelFolderOrganizer::setModuleFolder(), thePixelStructure, and twoDimOn.

Referenced by bookHistograms().

                                                                                        {
   
   std::map<uint32_t,SiPixelRecHitModule*>::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_,iBooker, iSetup, 0, twoDimOn, reducedSet, isUpgrade);
       } else {
     if(!isPIB) throw cms::Exception("LogicError")
       << "[SiPixelDigiSource::bookMEs] Creation of DQM folder failed";
       }
     }
     if(ladOn){
       if(theSiPixelFolder.setModuleFolder(iBooker,(*struct_iter).first,1,isUpgrade)){
     (*struct_iter).second->book( conf_,iBooker, iSetup, 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_,iBooker, iSetup, 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_,iBooker, iSetup, 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_,iBooker, iSetup, 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_,iBooker, iSetup, 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_,iBooker, iSetup, 6, twoDimOn, reducedSet, isUpgrade);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH RING-FOLDER\n";
       }
     }
 
   }
 
 }

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

virtual

Definition at line 174 of file SiPixelRecHitSource.cc.

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

Referenced by dqmBeginRun().

                                                                    {
 
   LogInfo ("PixelDQM") <<" SiPixelRecHitSource::buildStructure" ;
 
 
   edm::ESHandle<TrackerGeometry> pDD;
   edm::ESHandle<TrackerTopology> tTopoHandle;
 
   iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
   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();
       
       if((detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) ||
      (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap))){ 
 
     uint32_t id = detId();
     SiPixelRecHitModule* theModule = new SiPixelRecHitModule(id);
 
     if(detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
       if(isPIB) continue;
       LogDebug ("PixelDQM") << " ---> Adding Barrel Module " <<  detId.rawId() << endl;
       thePixelStructure.insert(pair<uint32_t,SiPixelRecHitModule*> (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();
       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;
       if(isPIB && mask) continue;
     
       thePixelStructure.insert(pair<uint32_t,SiPixelRecHitModule*> (id,theModule));
     }
 
       }
     }
   }//FOR_LOOP
 
   LogInfo ("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
 }

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

override

Definition at line 87 of file SiPixelRecHitSource.cc.

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

                                                                                  {
 
   LogInfo ("PixelDQM") << " SiPixelRecHitSource::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";
   
   if(firstRun){
     eventNo = 0;
     // Build map
     buildStructure(iSetup);
     // Book Monitoring Elements
     firstRun = false;
   }
 }