#include <SiPixelRecHitSource.h>

Inheritance diagram for SiPixelRecHitSource:

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

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

virtual void	bookMEs (DQMStore::IBooker &)

virtual void	buildStructure (edm::EventSetup const &)

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

	SiPixelRecHitSource (const edm::ParameterSet &conf)

	~SiPixelRecHitSource ()

Public Member Functions inherited from thread_unsafe::DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

	DQMEDAnalyzer (void)

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

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 50 of file SiPixelRecHitSource.cc.

References firstRun.

                                                                        :
   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;
 }

SiPixelRecHitSource::~SiPixelRecHitSource ( )

Definition at line 72 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 = 0;
   }
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 110 of file SiPixelRecHitSource.cc.

References edmNew::DetSet< T >::begin(), bladeOn, diskOn, edmNew::DetSet< T >::end(), eventNo, edm::Event::getByToken(), ladOn, layOn, match(), modOn, phiOn, rechit_count, reducedSet, ringOn, 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 &	,
		edm::EventSetup const &
	)

overridevirtual

Implements thread_unsafe::DQMEDAnalyzer.

Definition at line 103 of file SiPixelRecHitSource.cc.

References bookMEs().

                                                                                                         {
   bookMEs(iBooker);
 }

void SiPixelRecHitSource::bookMEs ( DQMStore::IBooker & iBooker )

virtual

Create folder tree and book histograms

Definition at line 251 of file SiPixelRecHitSource.cc.

References bladeOn, conf_, diskOn, edm::hlt::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,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,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,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,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,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,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,6,twoDimOn, reducedSet, isUpgrade);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH RING-FOLDER\n";
       }
     }
 
   }
 
 }

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

virtual

Definition at line 172 of file SiPixelRecHitSource.cc.

References PixelEndcapNameUpgrade::bladeName(), PixelEndcapName::bladeName(), PixelEndcapNameUpgrade::diskName(), PixelEndcapName::diskName(), edm::EventSetup::get(), PixelEndcapNameUpgrade::halfCylinder(), PixelEndcapName::halfCylinder(), isPIB, isUpgrade, LogDebug, PixelEndcapNameUpgrade::pannelName(), PixelEndcapName::pannelName(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapNameUpgrade::plaquetteName(), PixelEndcapName::plaquetteName(), DetId::rawId(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), and thePixelStructure.

Referenced by dqmBeginRun().

                                                                    {
 
   LogInfo ("PixelDQM") <<" SiPixelRecHitSource::buildStructure" ;
   edm::ESHandle<TrackerGeometry> pDD;
   iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
 
   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))!=0){
 
       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)) && (!isUpgrade)) {
       LogDebug ("PixelDQM") << " ---> Adding Endcap Module " <<  detId.rawId() << endl;
       PixelEndcapName::HalfCylinder side = PixelEndcapName(DetId(id)).halfCylinder();
       int disk   = PixelEndcapName(DetId(id)).diskName();
       int blade  = PixelEndcapName(DetId(id)).bladeName();
       int panel  = PixelEndcapName(DetId(id)).pannelName();
       int module = PixelEndcapName(DetId(id)).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));
     }   else if( (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) && (isUpgrade)) {
       LogDebug ("PixelDQM") << " ---> Adding Endcap Module " <<  detId.rawId() << endl;
       PixelEndcapNameUpgrade::HalfCylinder side = PixelEndcapNameUpgrade(DetId(id)).halfCylinder();
       int disk   = PixelEndcapNameUpgrade(DetId(id)).diskName();
       int blade  = PixelEndcapNameUpgrade(DetId(id)).bladeName();
       int panel  = PixelEndcapNameUpgrade(DetId(id)).pannelName();
       int module = PixelEndcapNameUpgrade(DetId(id)).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));
     }//endif(isUpgrade)
       }
     }       
   }
 
   LogInfo ("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
 }

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

virtual

Reimplemented from thread_unsafe::DQMEDAnalyzer.

Definition at line 85 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;
   }
 }