#include <CalibTracker/SiPixelGainCalibration/src/SiPixelCalibDigiProducer.cc>

Inheritance diagram for SiPixelCalibDigiProducer:

Public Member Functions
	SiPixelCalibDigiProducer (const edm::ParameterSet &iConfig)

	~SiPixelCalibDigiProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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 Types
typedef std::pair< uint32_t, std::pair< short, short > >	pixelstruct

Private Member Functions
virtual bool	checkFED (uint32_t detid)

virtual bool	checkPixel (uint32_t detid, short row, short col)

virtual void	clear ()

virtual void	fill (edm::Event &iEvent, const edm::EventSetup &iSetup)

virtual void	fillPixel (uint32_t detid, short row, short col, short ipoint, short adc)

virtual void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

virtual void	setPattern ()

virtual bool	store ()

Private Attributes
edm::ESHandle < SiPixelCalibConfiguration >	calib_

edm::ParameterSet	conf_

bool	control_pattern_size_

std::pair< short, short >	currentpair_

std::vector< std::pair< short, short > >	currentpattern_

std::map< uint32_t, uint32_t >	detid_to_fedid_

std::vector< pixelstruct >	detPixelMap_

std::map< pixelstruct, SiPixelCalibDigiError >	error_data_

int	errorType

uint32_t	iEventCounter_

bool	ignore_non_pattern_

bool	includeErrors_

std::string	instance_

std::map< pixelstruct, SiPixelCalibDigi >	intermediate_data_

std::string	label_

uint32_t	number_of_pixels_per_pattern_

uint32_t	pattern_repeat_

edm::InputTag	src_

edm::ESHandle < SiPixelFedCablingMap >	theCablingMap_

edm::ESHandle< TrackerGeometry >	theGeometry_

edm::EDGetTokenT < edm::DetSetVector< PixelDigi > >	tPixelDigi

bool	use_realeventnumber_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 55 of file SiPixelCalibDigiProducer.h.

Member Typedef Documentation

typedef std::pair<uint32_t,std::pair<short,short> > SiPixelCalibDigiProducer::pixelstruct

private

Definition at line 70 of file SiPixelCalibDigiProducer.h.

Constructor & Destructor Documentation

SiPixelCalibDigiProducer::SiPixelCalibDigiProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 54 of file SiPixelCalibDigiProducer.cc.

References includeErrors_, src_, and tPixelDigi.

                                                                                 :
   src_(iConfig.getParameter<edm::InputTag>("src")),
   iEventCounter_(0),
   ignore_non_pattern_(iConfig.getParameter<bool>("ignoreNonPattern")),
   control_pattern_size_(iConfig.getParameter<bool>("checkPatternEachEvent")),
   includeErrors_(iConfig.getUntrackedParameter<bool>("includeErrors",false)),
   errorType(iConfig.getUntrackedParameter<int>("errorTypeNumber",1)),
   conf_(iConfig),
   number_of_pixels_per_pattern_(0),
   use_realeventnumber_(iConfig.getParameter<bool>("useRealEventNumber"))
 
 {
   tPixelDigi = consumes<edm::DetSetVector<PixelDigi>> (src_);
    //register your products
   produces< edm::DetSetVector<SiPixelCalibDigi> >();
   if(includeErrors_)
     produces< edm::DetSetVector<SiPixelCalibDigiError> > ();
 
 }

SiPixelCalibDigiProducer::~SiPixelCalibDigiProducer ( )

Definition at line 75 of file SiPixelCalibDigiProducer.cc.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

bool SiPixelCalibDigiProducer::checkFED ( uint32_t detid )

privatevirtual

Definition at line 134 of file SiPixelCalibDigiProducer.cc.

References cond::rpcobgas::detid, detid_to_fedid_, edm::ESHandle< class >::product(), and theCablingMap_.

Referenced by fillPixel().

                                                      {
   //  edm::LogInfo("SiPixelCalibProducer") << "in checkFED" << std::endl;
 
   if(detid_to_fedid_[detid])
     return true;
   for(int fedid=0; fedid<=40; ++fedid){
     //    edm::LogInfo("SiPixelCalibProducer") << " looking at fedid " << fedid << std::endl;
     SiPixelFrameConverter converter(theCablingMap_.product(),fedid);
     if(converter.hasDetUnit(detid)){
       detid_to_fedid_[detid]=fedid;
       edm::LogInfo("SiPixelCalibDigiProducer") << "matched detid " << detid << " to fed " << detid_to_fedid_[detid] << std::endl;
       return true;
     }
   }
   return false;
 }

bool SiPixelCalibDigiProducer::checkPixel	(	uint32_t	detid,
		short	row,
		short	col
	)

privatevirtual

Definition at line 325 of file SiPixelCalibDigiProducer.cc.

References cuy::col, control_pattern_size_, currentpair_, currentpattern_, sipixelobjects::ElectronicIndex::dcol, sipixelobjects::LocalPixel::DcolPxid::dcol, cond::rpcobgas::detid, detid_to_fedid_, error_data_, i, iEventCounter_, ignore_non_pattern_, includeErrors_, edm::ESHandle< class >::product(), sipixelobjects::ElectronicIndex::pxid, sipixelobjects::LocalPixel::DcolPxid::pxid, sipixelobjects::LocalPixel::rocCol(), sipixelobjects::LocalPixel::rocRow(), setPattern(), store(), and theCablingMap_.

Referenced by fillPixel(), and produce().

                                                                              {
 
   if(!control_pattern_size_ && !store())
     return true;
   
   if( !ignore_non_pattern_ )
     return true;
   
   
   edm::LogInfo("SiPixelCalibDigiProducer") << "Event" << iEventCounter_ << ",now in checkpixel() " << std::endl;
   if(currentpattern_.size()==0)
     setPattern();
   //  uint32_t iroc;
   uint32_t fedid = detid_to_fedid_[detid];
 
   SiPixelFrameConverter formatter(theCablingMap_.product(),fedid);
   sipixelobjects::DetectorIndex detector ={detid, row, col};
   sipixelobjects::ElectronicIndex cabling;
   
   formatter.toCabling(cabling,detector);
   // cabling should now contain cabling.roc and cabling.dcol  and cabling.pxid
 
   // however, the coordinates now need to be converted from dcl, pxid to the row,col coordinates used in the calibration info
   sipixelobjects::LocalPixel::DcolPxid loc;
   loc.dcol = cabling.dcol;
   loc.pxid = cabling.pxid;
   sipixelobjects::LocalPixel locpixel(loc);
   currentpair_.first = locpixel.rocRow();
   currentpair_.second = locpixel.rocCol();
 
   for(uint32_t i=0; i<currentpattern_.size(); ++i){
     //    edm::LogInfo("SiPixelCalibDigiProducer") << "found pair " << currentpair_.first << "," << currentpair_.second << " calib " << currentpattern_[i].first << ","<< currentpattern_[i].second << " input " << row << "," << col << std::endl;
     if(currentpair_==currentpattern_[i]){
       return true;
     }
   }
   std::ostringstream errorlog;
   errorlog <<  "DETID "<<detid<<", row, col (offline)="<<row<<","<<col<<" row, col (ROC) ="<<currentpair_.first<<","<< currentpair_.second<< " found no match in list of patterns: " ;
   for(uint32_t i=0; i<currentpattern_.size(); ++i){
     if(i!=0 && i!=currentpattern_.size()-1)
       errorlog<<" ";
     errorlog<<"(";
     errorlog<<currentpattern_[i].first;
     errorlog<<",";
     errorlog<<currentpattern_[i].second;
     errorlog<<")";
   }
   edm::LogError("ERROR") << errorlog.str() << std::endl;
   if(includeErrors_){// book the error
     
     pixelstruct temppixelworker;
     temppixelworker.first=detid;
     temppixelworker.second.first=row;
     temppixelworker.second.second=col;
     std::map<pixelstruct,SiPixelCalibDigiError>::const_iterator ierr = error_data_.find(temppixelworker);
     if(ierr== error_data_.end()){
       SiPixelCalibDigiError temperr(row,col,1);
       error_data_[temppixelworker]=temperr;
     }
   }
     
   return false;
 }

void SiPixelCalibDigiProducer::clear ( void )

privatevirtual

Definition at line 187 of file SiPixelCalibDigiProducer.cc.

References error_data_, intermediate_data_, and number_of_pixels_per_pattern_.

Referenced by Vispa.Views.WidgetView.WidgetView::closeEvent(), Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), Vispa.Share.FindAlgorithm.FindAlgorithm::findUsingFindDialog(), produce(), Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), Vispa.Views.WidgetView.WidgetView::setDataObjects(), BeautifulSoup.Tag::setString(), Vispa.Views.TreeView.TreeView::updateContent(), Vispa.Views.TableView.TableView::updateContent(), Vispa.Views.BoxDecayView.BoxDecayView::updateContent(), and Vispa.Views.PropertyView.PropertyView::updateContent().

                                {
   //  edm::LogInfo("SiPixelCalibProducer") << "in clear() " << std::endl;
   // this is where we empty the containers so they can be re-filled
   // the idea: the detPixelMap_ container shrinks/expands as a function 
   // of the number of pixels looked at at one particular time... 
   // unlike the intermediate_data_ container which only expands when 
   // detPixelMap_ becomes bigger than intermedate_data_
   
   // shrink the detPixelMap_
   uint32_t tempsize = intermediate_data_.size();
   if(tempsize>number_of_pixels_per_pattern_){
     edm::LogError("SiPixelCalibDigiProducer") << "Number of pixels in pattern is now: " << tempsize<< ", size is was " << number_of_pixels_per_pattern_ << std::endl;
     number_of_pixels_per_pattern_=tempsize;
   }
 
   intermediate_data_.erase(intermediate_data_.begin(),intermediate_data_.end());
   intermediate_data_.clear();
 
   // and erase the error bits
   error_data_.erase(error_data_.begin(), error_data_.end());
   error_data_.clear();
 }

void SiPixelCalibDigiProducer::fill	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Definition at line 109 of file SiPixelCalibDigiProducer.cc.

References edm::DetSetVector< T >::begin(), calib_, cond::rpcobgas::detid, edm::DetSetVector< T >::end(), fillPixel(), edm::Event::getByToken(), iEventCounter_, and tPixelDigi.

Referenced by produce().

 {
 
   // figure out which calibration point we're on now..
   short icalibpoint = calib_->vcalIndexForEvent(iEventCounter_);
   edm::Handle< edm::DetSetVector<PixelDigi> > pixelDigis;
   iEvent.getByToken( tPixelDigi, pixelDigis );
   
   edm::LogInfo("SiPixelCalibProducer") << "in fill(), calibpoint " << icalibpoint <<" ndigis " << pixelDigis->size() <<  std::endl;
     // loop over the data and store things
   edm::DetSetVector<PixelDigi>::const_iterator digiIter;
   for(digiIter=pixelDigis->begin(); digiIter!=pixelDigis->end(); ++digiIter){// ITERATOR OVER DET IDs
     uint32_t detid = digiIter->id;
     edm::DetSet<PixelDigi>::const_iterator ipix; // ITERATOR OVER DIGI DATA  
 
     for(ipix = digiIter->data.begin(); ipix!=digiIter->end(); ++ipix){
       // fill in the appropriate location of the temporary data container
       fillPixel(detid,ipix->row(),ipix->column(),icalibpoint,ipix->adc());
     }
   }
 }

void SiPixelCalibDigiProducer::fillPixel	(	uint32_t	detid,
		short	row,
		short	col,
		short	ipoint,
		short	adc
	)

privatevirtual

Definition at line 154 of file SiPixelCalibDigiProducer.cc.

References calib_, checkFED(), checkPixel(), cuy::col, cond::rpcobgas::detid, intermediate_data_, and SiPixelCalibDigi::setrowcol().

Referenced by fill().

                                                                                                      {
   //  edm::LogInfo("SiPixelCalibProducer") << " in fillpixel()" << std::endl;
  
   //  edm::LogInfo("SiPixelCalibProducer") << "in fillPixel " << detid << " " << row << " " << col << " " << ipoint << " " << adc << std::endl;
   if(!checkFED(detid)){
     edm::LogError("SiPixelCalibDigiProducer") << " was unable to match detid " << detid << " to a FED!" << std::endl;
     return;
   }
   if(!checkPixel(detid,row,col)){
     return;
   }
   // now the check if the pixel exists and fill
   //
   pixelstruct temppixelworker;
   temppixelworker.first=detid;
   temppixelworker.second.first=row;
   temppixelworker.second.second=col;  
   std::map<pixelstruct,SiPixelCalibDigi>::const_iterator ipix = intermediate_data_.find(temppixelworker);
   
   if(ipix == intermediate_data_.end()){
     SiPixelCalibDigi tempdigi(calib_->nVCal());
     tempdigi.setrowcol(row,col);
     intermediate_data_[temppixelworker]=tempdigi;
   }
 
   intermediate_data_[temppixelworker].fill(ipoint,adc);
   return;
 
 }

void SiPixelCalibDigiProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 271 of file SiPixelCalibDigiProducer.cc.

References calib_, checkPixel(), clear(), control_pattern_size_, edm::DetSet< T >::data, cond::rpcobgas::detid, error_data_, edm::EventID::event(), fill(), edm::EventSetup::get(), edm::EventBase::id(), iEventCounter_, includeErrors_, intermediate_data_, pattern_repeat_, edm::Event::put(), setPattern(), store(), theCablingMap_, theGeometry_, and use_realeventnumber_.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

 {
   //  edm::LogInfo("SiPixelCalibDigiProducer") <<"in produce() " << std::endl;
   using namespace edm;
   iSetup.get<SiPixelCalibConfigurationRcd>().get(calib_);
   iSetup.get<TrackerDigiGeometryRecord>().get( theGeometry_ );
   iSetup.get<SiPixelFedCablingMapRcd>().get(theCablingMap_);
   pattern_repeat_=calib_->getNTriggers()*calib_->nVCal();
   if(use_realeventnumber_){
     iEventCounter_= iEvent.id().event()-1;
   }
   else
     iEventCounter_++;
   if(iEventCounter_%pattern_repeat_==1)
     setPattern();
   
   //  edm::LogInfo("SiPixelCalibDigiProducer") << "now starting fill..." << std::endl;
   fill(iEvent,iSetup); // fill method where the actual looping over the digis is done.
   //  edm::LogInfo("SiPixelCalibDigiProducer") << "done filling..." << std::endl;
   std::auto_ptr<edm::DetSetVector<SiPixelCalibDigi> > pOut(new edm::DetSetVector<SiPixelCalibDigi>);
   std::auto_ptr<edm::DetSetVector<SiPixelCalibDigiError> > pErr (new edm::DetSetVector<SiPixelCalibDigiError> );
   
   // copy things over into pOut if necessary (this is only once per pattern)
   if(store()){
     //    edm::LogInfo("SiPixelCalibDigiProducer") << "in loop" << std::endl;
     for(std::map<pixelstruct,SiPixelCalibDigi>::const_iterator idet=intermediate_data_.begin(); idet!=intermediate_data_.end();++idet){
       uint32_t detid=idet->first.first;
       if(!control_pattern_size_){
     if(! checkPixel(idet->first.first,idet->first.second.first,idet->first.second.second))
       continue;
       }
 
       
       SiPixelCalibDigi tempdigi=idet->second;
       edm::DetSet<SiPixelCalibDigi> & detSet = pOut->find_or_insert(detid);
       detSet.data.push_back(tempdigi);
     }
     if(includeErrors_){
       for(std::map<pixelstruct,SiPixelCalibDigiError>::const_iterator ierr=error_data_.begin(); ierr!=error_data_.end();++ierr){
     uint32_t detid=ierr->first.first;
     SiPixelCalibDigiError temperror = ierr->second;
     edm::DetSet<SiPixelCalibDigiError> & errSet = pErr->find_or_insert(detid);
     errSet.data.push_back(temperror);
       }
     }
     edm::LogInfo("INFO") << "now filling event " << iEventCounter_ << " as pixel pattern changes every " <<  pattern_repeat_ << " events..." << std::endl;
     clear();
   }
   iEvent.put(pOut);
   if(includeErrors_)
     iEvent.put(pErr);
 }

void SiPixelCalibDigiProducer::setPattern ( )

privatevirtual

Definition at line 214 of file SiPixelCalibDigiProducer.cc.

References funct::abs(), calib_, currentpattern_, iEventCounter_, chain::pattern, and pattern_repeat_.

Referenced by checkPixel(), and produce().

                                     {
   //  edm::LogInfo("SiPixelCalibProducer") << "in setPattern()" << std::endl;
   uint32_t patternnumber = abs(iEventCounter_-1)/pattern_repeat_;
   uint32_t rowpatternnumber = patternnumber/calib_->nColumnPatterns();
   uint32_t colpatternnumber = patternnumber%calib_->nColumnPatterns();
   edm::LogInfo("SiPixelCalibDigiProducer") << " rowpatternnumbers = " << rowpatternnumber << " " << colpatternnumber << " " << patternnumber << std::endl;
   // update currentpattern_
   std::vector<short> calibcols = calib_->getColumnPattern();
   std::vector<short> calibrows = calib_->getRowPattern();
   std::vector<short> temprowvals(0);
   std::vector<short> tempcolvals(0);
   uint32_t nminuscol=0;
   uint32_t nminusrow=0;
   uint32_t npatterns=0;
   for(uint32_t icol=0; icol<calibcols.size(); icol++){
     if(calibcols[icol]==-1){
       nminuscol++;
     }
     else if(nminuscol==colpatternnumber){
       //edm::LogInfo("SiPixelCalibProducer") << "col " << calibcols[icol] << std::endl;
       short val=calibcols[icol];
       tempcolvals.push_back(val);
     }
     else if (nminuscol> colpatternnumber)
       break;
   }
   for(uint32_t irow=0; irow<calibrows.size(); irow++){
     // edm::LogInfo("SiPixelCalibProducer") << "row " << irow <<" "<< nminusrow<<" "  << calibrows[irow] << std::endl;
     if(calibrows[irow]==-1)
       nminusrow++;
     else if(nminusrow==rowpatternnumber){
       short val=calibrows[irow];
       temprowvals.push_back(val);
     }
     else if(nminusrow>rowpatternnumber)
       break;
   }
   //now clean up the currentpattern_;
   while(currentpattern_.size()>temprowvals.size()*tempcolvals.size()){
     currentpattern_.erase(currentpattern_.end());
   }
   for(uint32_t irow=0; irow<temprowvals.size(); irow++){
     for(uint32_t icol=0; icol<tempcolvals.size(); icol++){
       std::pair<short,short> pattern(temprowvals[irow],tempcolvals[icol]);
       npatterns++;
       if(npatterns>currentpattern_.size())
     currentpattern_.push_back(pattern);
       else
     currentpattern_[npatterns-1]=pattern;
     }
   }
 }

bool SiPixelCalibDigiProducer::store ( )

privatevirtual

Definition at line 92 of file SiPixelCalibDigiProducer.cc.

References calib_, iEventCounter_, and pattern_repeat_.

Referenced by checkPixel(), and produce().

 {
   //  std::cout << "in store() " << std::endl;
   if(iEventCounter_%pattern_repeat_==0){
     //    std::cout << "now at event " << iEventCounter_ <<" where we save the calibration information into the CMSSW digi";
     return 1;
   }
   else if(iEventCounter_==calib_->expectedTotalEvents())
     return 1;
   else
     return 0;
   return 1;
 }