#include <SiPixelRecHitSource.h>

Inheritance diagram for SiPixelRecHitSource:

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

virtual void	beginJob ()

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

virtual void	bookMEs ()

virtual void	buildStructure (edm::EventSetup const &)

virtual void	endJob ()

	SiPixelRecHitSource (const edm::ParameterSet &conf)

	~SiPixelRecHitSource ()

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

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	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_

DQMStore *	theDMBE

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 53 of file SiPixelRecHitSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 49 of file SiPixelRecHitSource.cc.

References cppFunctionSkipper::operator, and theDMBE.

                                                                        :
   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) )
 {
    theDMBE = edm::Service<DQMStore>().operator->();
    LogInfo ("PixelDQM") << "SiPixelRecHitSource::SiPixelRecHitSource: Got DQM BackEnd interface"<<endl;
 }

SiPixelRecHitSource::~SiPixelRecHitSource ( )

Definition at line 70 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 122 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.getByLabel( src_, recHitColl ); //original
   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;
       //cout << "hh" << endl;
       (*struct_iter).second->fill(rechit_x, rechit_y, sizeX, sizeY, lerr_x, lerr_y, 
                                   modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn, 
                       twoDimOn, reducedSet);
       //cout << "ii" << endl;
     
     }
       }
       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::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 83 of file SiPixelRecHitSource.cc.

References firstRun.

                                   {
   firstRun = true;
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 88 of file SiPixelRecHitSource.cc.

References bladeOn, bookMEs(), 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
     bookMEs();
     firstRun = false;
   }
 }

void SiPixelRecHitSource::bookMEs ( )

virtual

Create folder tree and book histograms

Definition at line 251 of file SiPixelRecHitSource.cc.

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

Referenced by beginRun().

                                  {
   
   std::map<uint32_t,SiPixelRecHitModule*>::iterator struct_iter;
     
   SiPixelFolderOrganizer theSiPixelFolder;
   
   for(struct_iter = thePixelStructure.begin(); struct_iter != thePixelStructure.end(); struct_iter++){
     
     if(modOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first)){
     (*struct_iter).second->book( conf_,0,twoDimOn, reducedSet);
       } else {
     if(!isPIB) throw cms::Exception("LogicError")
       << "[SiPixelDigiSource::bookMEs] Creation of DQM folder failed";
       }
     }
     if(ladOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,1)){
     (*struct_iter).second->book( conf_,1,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH LADDER-FOLDER\n";
       }
     }
     if(layOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,2)){
     (*struct_iter).second->book( conf_,2,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH LAYER-FOLDER\n";
       }
     }
     if(phiOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,3)){
     (*struct_iter).second->book( conf_,3,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH PHI-FOLDER\n";
       }
     }
     if(bladeOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,4)){
     (*struct_iter).second->book( conf_,4,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH BLADE-FOLDER\n";
       }
     }
     if(diskOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,5)){
     (*struct_iter).second->book( conf_,5,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH DISK-FOLDER\n";
       }
     }
     if(ringOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,6)){
     (*struct_iter).second->book( conf_,6,twoDimOn, reducedSet);
     } else {
     LogDebug ("PixelDQM") << "PROBLEM WITH RING-FOLDER\n";
       }
     }
 
   }
 
 }

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

virtual

Definition at line 187 of file SiPixelRecHitSource.cc.

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

Referenced by beginRun().

                                                                    {
 
   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*>((*it))!=0){
 
       DetId detId = (*it)->geographicalId();
       // const GeomDetUnit      * geoUnit = pDD->idToDetUnit( detId );
       //const PixelGeomDetUnit * pixDet  = dynamic_cast<const PixelGeomDetUnit*>(geoUnit);
 
           
       
           // SiPixelRecHitModule *theModule = new SiPixelRecHitModule(id, rechit_x, rechit_y, x_res, y_res, x_pull, y_pull);
     
     
             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)).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));
           }
         }
     }       
   }
 
   LogInfo ("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
 }

void SiPixelRecHitSource::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 108 of file SiPixelRecHitSource.cc.

References conf_, edm::ParameterSet::getParameter(), download_sqlite_cfg::outputFile, DQMStore::save(), saveFile, AlCaHLTBitMon_QueryRunRegistry::string, and theDMBE.

                                     {
 
 
   if(saveFile){
     LogInfo ("PixelDQM") << " SiPixelRecHitSource::endJob - Saving Root File " << std::endl;
     std::string outputFile = conf_.getParameter<std::string>("outputFile");
     theDMBE->save( outputFile );
   }
 
 }