#include <SiPixelRawDataErrorSource.h>

Inheritance diagram for SiPixelRawDataErrorSource:

Public Types
typedef edm::DetSet < SiPixelRawDataError > ::const_iterator	ErrorIterator

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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 ()

	SiPixelRawDataErrorSource (const edm::ParameterSet &conf)

	~SiPixelRawDataErrorSource ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Attributes
bool	bladeOn

MonitorElement *	byLumiErrors

edm::ParameterSet	conf_

MonitorElement *	errorRate

int	eventNo

bool	firstRun

bool	isPIB

bool	ladOn

bool	modOn

bool	reducedSet

bool	saveFile

bool	slowDown

edm::InputTag	src_

DQMStore *	theDMBE

std::map< uint32_t, SiPixelRawDataErrorModule * >	theFEDStructure

std::map< uint32_t, SiPixelRawDataErrorModule * >	thePixelStructure

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Description: Produces histograms for error information generated at the raw2digi stage for the pixel tracker.

Usage: Takes a DetSetVector<SiPixelRawDataError> as input, and uses it to populate a folder hierarchy (organized by detId) with histograms containing information about the errors. Uses SiPixelRawDataErrorModule class to book and fill individual folders. Note that this source is different than other DQM sources in the creation of an unphysical detId folder (detId=0xffffffff) to hold information about errors for which there is no detId available (except the dummy detId given to it at raw2digi).

Definition at line 53 of file SiPixelRawDataErrorSource.h.

Member Typedef Documentation

typedef edm::DetSet<SiPixelRawDataError>::const_iterator SiPixelRawDataErrorSource::ErrorIterator

Definition at line 58 of file SiPixelRawDataErrorSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 49 of file SiPixelRawDataErrorSource.cc.

References cmsCodeRules.cppFunctionSkipper::operator, and theDMBE.

                                                                                    :
   conf_(iConfig),
   src_( conf_.getParameter<edm::InputTag>( "src" ) ),
   saveFile( conf_.getUntrackedParameter<bool>("saveFile",false) ),
   isPIB( conf_.getUntrackedParameter<bool>("isPIB",false) ),
   slowDown( conf_.getUntrackedParameter<bool>("slowDown",false) ),
   reducedSet( conf_.getUntrackedParameter<bool>("reducedSet",false) ),
   modOn( conf_.getUntrackedParameter<bool>("modOn",true) ),
   ladOn( conf_.getUntrackedParameter<bool>("ladOn",false) ), 
   bladeOn( conf_.getUntrackedParameter<bool>("bladeOn",false) )
 {
    theDMBE = edm::Service<DQMStore>().operator->();
    LogInfo ("PixelDQM") << "SiPixelRawDataErrorSource::SiPixelRawDataErrorSource: Got DQM BackEnd interface"<<endl;
 }

SiPixelRawDataErrorSource::~SiPixelRawDataErrorSource ( )

Definition at line 65 of file SiPixelRawDataErrorSource.cc.

 {
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
   LogInfo ("PixelDQM") << "SiPixelRawDataErrorSource::~SiPixelRawDataErrorSource: Destructor"<<endl;
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 107 of file SiPixelRawDataErrorSource.cc.

References bladeOn, eventNo, MonitorElement::Fill(), DQMStore::get(), MonitorElement::getBinContent(), edm::Event::getByLabel(), collect_tpl::input, edm::HandleBase::isValid(), ladOn, edm::EventBase::luminosityBlock(), modOn, cmsCodeRules.cppFunctionSkipper::operator, MonitorElement::setBinContent(), slowDown, src_, theDMBE, theFEDStructure, and thePixelStructure.

 {
   eventNo++;
   //std::cout<<"Event number: "<<eventNo<<std::endl;
   // get input data
   edm::Handle< edm::DetSetVector<SiPixelRawDataError> >  input;
   iEvent.getByLabel( src_, input );
   if (!input.isValid()) return; 
 
   // Get DQM interface
   DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
    
   //float iOrbitSec = iEvent.orbitNumber()/11223.;
   //int bx = iEvent.bunchCrossing();
   //long long tbx = (long long)iEvent.orbitNumber() * 3564 + bx;
   int lumiSection = (int)iEvent.luminosityBlock();
   
   int nEventBPIXModuleErrors = 0; int nEventFPIXModuleErrors = 0; int nEventBPIXFEDErrors = 0; int nEventFPIXFEDErrors = 0;
   int nErrors = 0;
 
   std::map<uint32_t,SiPixelRawDataErrorModule*>::iterator struct_iter;
   std::map<uint32_t,SiPixelRawDataErrorModule*>::iterator struct_iter2;
   for (struct_iter = thePixelStructure.begin() ; struct_iter != thePixelStructure.end() ; struct_iter++) {
     
     int numberOfModuleErrors = (*struct_iter).second->fill(*input, modOn, ladOn, bladeOn);
     if((*struct_iter).first >= 302055684 && (*struct_iter).first <= 302197792 ) nEventBPIXModuleErrors = nEventBPIXModuleErrors + numberOfModuleErrors;
     if((*struct_iter).first >= 343999748 && (*struct_iter).first <= 352477708 ) nEventFPIXModuleErrors = nEventFPIXModuleErrors + numberOfModuleErrors;
     nErrors = nErrors + numberOfModuleErrors;
   }
   for (struct_iter2 = theFEDStructure.begin() ; struct_iter2 != theFEDStructure.end() ; struct_iter2++) {
     
     int numberOfFEDErrors = (*struct_iter2).second->fillFED(*input);
     if((*struct_iter2).first <= 31) nEventBPIXFEDErrors = nEventBPIXFEDErrors + numberOfFEDErrors; // (*struct_iter2).first >= 0, since (*struct_iter2).first is unsigned
     if((*struct_iter2).first >= 32 && (*struct_iter2).first <= 39) nEventFPIXFEDErrors = nEventFPIXFEDErrors + numberOfFEDErrors;    
     nErrors = nErrors + numberOfFEDErrors;
   }
   MonitorElement* me = theDMBE->get("Pixel/AdditionalPixelErrors/byLumiErrors");
   if(me){
     me->setBinContent(0,eventNo);
     if(nEventBPIXModuleErrors+nEventBPIXFEDErrors>0) me->setBinContent(1,nEventBPIXModuleErrors+nEventBPIXFEDErrors);
     if(nEventFPIXModuleErrors+nEventFPIXFEDErrors>0) me->setBinContent(2,nEventFPIXModuleErrors+nEventFPIXFEDErrors);
   }  
   
   // Rate of errors per lumi section:
   MonitorElement* me1 = theDMBE->get("Pixel/AdditionalPixelErrors/errorRate");
   if(me1){
     int nLumiErrors = int(me1->getBinContent(lumiSection+1)) + nErrors;
     me1->Fill(lumiSection, nLumiErrors);
   }
   
     
 
   // slow down...
   if(slowDown) usleep(100000);
   
 }

void SiPixelRawDataErrorSource::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 73 of file SiPixelRawDataErrorSource.cc.

References firstRun.

                                         {
   firstRun = true;
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 77 of file SiPixelRawDataErrorSource.cc.

References bladeOn, bookMEs(), buildStructure(), eventNo, firstRun, ladOn, and modOn.

                                                                                     {
 
   LogInfo ("PixelDQM") << " SiPixelRawDataErrorSource::beginRun - Initialisation ... " << std::endl;
   LogInfo ("PixelDQM") << "Mod/Lad/Blade " << modOn << "/" << ladOn << "/" << bladeOn << std::endl;
 
   if(firstRun){
     eventNo = 0;
     // Build map
     buildStructure(iSetup);
     // Book Monitoring Elements
     bookMEs();
     
     firstRun = false;
   }
 }

void SiPixelRawDataErrorSource::bookMEs ( )

virtual

Create folder tree and book histograms

Create folder tree for errors without detId and book histograms

Definition at line 239 of file SiPixelRawDataErrorSource.cc.

References bladeOn, DQMStore::book1D(), byLumiErrors, conf_, errorRate, edm::hlt::Exception, isPIB, ladOn, LogDebug, modOn, cmsCodeRules.cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), SiPixelFolderOrganizer::setFedFolder(), MonitorElement::setLumiFlag(), SiPixelFolderOrganizer::setModuleFolder(), theDMBE, theFEDStructure, thePixelStructure, and indexGen::title.

Referenced by beginRun().

                                        {
   
   // Get DQM interface
   DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
   theDMBE->setCurrentFolder("Pixel/AdditionalPixelErrors");
   char title[80]; sprintf(title, "By-LumiSection Error counters");
   byLumiErrors = theDMBE->book1D("byLumiErrors",title,2,0.,2.);
   byLumiErrors->setLumiFlag();
   char title1[80]; sprintf(title1, "Errors per LumiSection;LumiSection;NErrors");
   errorRate = theDMBE->book1D("errorRate",title1,5000,0.,5000.);
   
   std::map<uint32_t,SiPixelRawDataErrorModule*>::iterator struct_iter;
   std::map<uint32_t,SiPixelRawDataErrorModule*>::iterator struct_iter2;
   
   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 );
       }
       else {
         //std::cout<<"PIB! not booking histograms for non-PIB modules!"<<std::endl;
         if(!isPIB) throw cms::Exception("LogicError")
                        << "[SiPixelRawDataErrorSource::bookMEs] Creation of DQM folder failed";
       }
     }
     
     if(ladOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,1)) {
         (*struct_iter).second->book( conf_, 1 );
       }
       else {
         LogDebug ("PixelDQM") << "PROBLEM WITH LADDER-FOLDER\n";
       }
     }
     
     if(bladeOn){
       if(theSiPixelFolder.setModuleFolder((*struct_iter).first,4)) {
         (*struct_iter).second->book( conf_, 4 );
       }
       else {
         LogDebug ("PixelDQM") << "PROBLEM WITH BLADE-FOLDER\n";
       }
     }
     
   }//for loop
 
   for(struct_iter2 = theFEDStructure.begin(); struct_iter2 != theFEDStructure.end(); struct_iter2++){
     if(theSiPixelFolder.setFedFolder((*struct_iter2).first)) {
       (*struct_iter2).second->bookFED( conf_ );
     }
     else {
       throw cms::Exception("LogicError")
     << "[SiPixelRawDataErrorSource::bookMEs] Creation of DQM folder failed";
     }
 
   }
 
 }

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

virtual

Definition at line 167 of file SiPixelRawDataErrorSource.cc.

References PixelEndcapName::bladeName(), PixelEndcapName::diskName(), edm::EventSetup::get(), PixelEndcapName::halfCylinder(), isPIB, LogDebug, FEDNumbering::MAXSiPixelFEDID, FEDNumbering::MINSiPixelFEDID, PixelEndcapName::pannelName(), GeomDetEnumerators::PixelBarrel, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, GeomDetEnumerators::PixelEndcap, PixelEndcapName::plaquetteName(), DetId::rawId(), PixelGeomDetUnit::specificTopology(), DetId::subdetId(), theFEDStructure, and thePixelStructure.

Referenced by beginRun().

                                                                          {
 
   LogInfo ("PixelDQM") <<" SiPixelRawDataErrorSource::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( ((*it)->subDetector()==GeomDetEnumerators::PixelBarrel) || ((*it)->subDetector()==GeomDetEnumerators::PixelEndcap) ){
       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)) {
         if(isPIB) continue;
         LogDebug ("PixelDQM") << " ---> Adding Barrel Module " <<  detId.rawId() << endl;
     uint32_t id = detId();
     SiPixelRawDataErrorModule* theModule = new SiPixelRawDataErrorModule(id, ncols, nrows);
     thePixelStructure.insert(pair<uint32_t,SiPixelRawDataErrorModule*> (id,theModule));
 
       } else if(detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
     LogDebug ("PixelDQM") << " ---> Adding Endcap Module " <<  detId.rawId() << endl;
     uint32_t id = detId();
     SiPixelRawDataErrorModule* theModule = new SiPixelRawDataErrorModule(id, ncols, nrows);
     
         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;
     // clutch to take all of FPIX, but no BPIX:
     mask = false;
     if(isPIB && mask) continue;
         
     thePixelStructure.insert(pair<uint32_t,SiPixelRawDataErrorModule*> (id,theModule));
       }
     }
   }
   LogDebug ("PixelDQM") << " ---> Adding Module for Additional Errors " << endl;
   pair<int,int> fedIds (FEDNumbering::MINSiPixelFEDID, FEDNumbering::MAXSiPixelFEDID);
   fedIds.first = 0;
   fedIds.second = 39;
   for (int fedId = fedIds.first; fedId <= fedIds.second; fedId++) {
     //std::cout<<"Adding FED module: "<<fedId<<std::endl;
     uint32_t id = static_cast<uint32_t> (fedId);
     SiPixelRawDataErrorModule* theModule = new SiPixelRawDataErrorModule(id);
     theFEDStructure.insert(pair<uint32_t,SiPixelRawDataErrorModule*> (id,theModule));
   }
   
   LogInfo ("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
 }

void SiPixelRawDataErrorSource::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 94 of file SiPixelRawDataErrorSource.cc.

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

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