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

Inheritance diagram for SiPixelCalibDigiProducer:

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

	~SiPixelCalibDigiProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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)

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

virtual void	setPattern ()

virtual bool	store ()

Private Attributes
edm::ESGetToken < SiPixelFedCablingMap, SiPixelFedCablingMapRcd >	cablingMapToken_

edm::ESHandle < SiPixelCalibConfiguration >	calib_

edm::ESGetToken < SiPixelCalibConfiguration, SiPixelCalibConfigurationRcd >	calibToken_

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

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_

bool	use_realeventnumber_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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 69 of file SiPixelCalibDigiProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 47 of file SiPixelCalibDigiProducer.cc.

References cablingMapToken_, calibToken_, includeErrors_, src_, tPixelDigi, and trackerGeomToken_.

     : 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>>();
 
   calibToken_ = esConsumes<SiPixelCalibConfiguration, SiPixelCalibConfigurationRcd>();
   trackerGeomToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
   cablingMapToken_ = esConsumes<SiPixelFedCablingMap, SiPixelFedCablingMapRcd>();
 }

SiPixelCalibDigiProducer::~SiPixelCalibDigiProducer ( )

override

Definition at line 70 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 120 of file SiPixelCalibDigiProducer.cc.

References 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 309 of file SiPixelCalibDigiProducer.cc.

References cuy::col, control_pattern_size_, currentpair_, currentpattern_, sipixelobjects::ElectronicIndex::dcol, sipixelobjects::LocalPixel::DcolPxid::dcol, detid_to_fedid_, error_data_, mps_fire::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_.empty())
     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 172 of file SiPixelCalibDigiProducer.cc.

References error_data_, intermediate_data_, and number_of_pixels_per_pattern_.

Referenced by produce(), and BeautifulSoup.Tag::setString().

                                      {
   //  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 96 of file SiPixelCalibDigiProducer.cc.

References edm::DetSetVector< T >::begin(), calib_, 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 141 of file SiPixelCalibDigiProducer.cc.

References calib_, checkFED(), checkPixel(), cuy::col, 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
	)

overrideprivate

Definition at line 253 of file SiPixelCalibDigiProducer.cc.

References cablingMapToken_, calib_, calibToken_, checkPixel(), clear(), control_pattern_size_, edm::DetSet< T >::data, error_data_, edm::EventID::event(), fill(), edm::EventSetup::getHandle(), edm::EventBase::id(), iEventCounter_, includeErrors_, intermediate_data_, eostools::move(), pattern_repeat_, edm::Event::put(), setPattern(), store(), theCablingMap_, theGeometry_, trackerGeomToken_, and use_realeventnumber_.

                                                                                     {
   //  edm::LogInfo("SiPixelCalibDigiProducer") <<"in produce() " << std::endl;
   using namespace edm;
   calib_ = iSetup.getHandle(calibToken_);
   theGeometry_ = iSetup.getHandle(trackerGeomToken_);
   theCablingMap_ = iSetup.getHandle(cablingMapToken_);
   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;
   auto pOut = std::make_unique<edm::DetSetVector<SiPixelCalibDigi>>();
   auto pErr = std::make_unique<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(std::move(pOut));
   if (includeErrors_)
     iEvent.put(std::move(pErr));
 }

void SiPixelCalibDigiProducer::setPattern ( )

privatevirtual

Definition at line 199 of file SiPixelCalibDigiProducer.cc.

References calib_, currentpattern_, iEventCounter_, pattern_repeat_, and hgcalPerformanceValidation::val.

Referenced by checkPixel(), and produce().

                                           {
   //  edm::LogInfo("SiPixelCalibProducer") << "in setPattern()" << std::endl;
   uint32_t patternnumber = (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 82 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 true;
   } else if (iEventCounter_ == calib_->expectedTotalEvents())
     return true;
   else
     return false;
   return true;
 }