Inheritance diagram for cms::SiPixelCondObjBuilder:

Classes
class	CalParameters

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

void	beginJob () override

bool	loadFromFile ()

	SiPixelCondObjBuilder (const edm::ParameterSet &iConfig)

	~SiPixelCondObjBuilder () override=default

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const bool	appendMode_

std::map< int, CalParameters, std::less< int > >	calmap_

const std::string	fileName_

const bool	fromFile_

const double	meanGain_

const double	meanPed_

const int	numberOfModules_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	pddToken_

PixelIndices *	pIndexConverter_

const std::string	recordName_

const double	rmsGain_

const double	rmsPed_

const double	secondRocRowGainOffset_

const double	secondRocRowPedOffset_

std::unique_ptr< SiPixelGainCalibration >	SiPixelGainCalibration_

SiPixelGainCalibrationService	SiPixelGainCalibrationService_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 45 of file SiPixelCondObjBuilder.cc.

Constructor & Destructor Documentation

◆ SiPixelCondObjBuilder()

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

explicit

Definition at line 84 of file SiPixelCondObjBuilder.cc.

       : appendMode_(iConfig.getUntrackedParameter<bool>("appendMode", true)),
         pddToken_(esConsumes()),
         SiPixelGainCalibration_(nullptr),
         SiPixelGainCalibrationService_(iConfig, consumesCollector()),
         recordName_(iConfig.getParameter<std::string>("record")),
         meanPed_(iConfig.getParameter<double>("meanPed")),
         rmsPed_(iConfig.getParameter<double>("rmsPed")),
         meanGain_(iConfig.getParameter<double>("meanGain")),
         rmsGain_(iConfig.getParameter<double>("rmsGain")),
         secondRocRowGainOffset_(iConfig.getParameter<double>("secondRocRowGainOffset")),
         secondRocRowPedOffset_(iConfig.getParameter<double>("secondRocRowPedOffset")),
         numberOfModules_(iConfig.getParameter<int>("numberOfModules")),
         fromFile_(iConfig.getParameter<bool>("fromFile")),
         fileName_(iConfig.getParameter<std::string>("fileName")) {
     ::putenv((char*)"CORAL_AUTH_USER=me");
     ::putenv((char*)"CORAL_AUTH_PASSWORD=test");
   }

◆ ~SiPixelCondObjBuilder()

cms::SiPixelCondObjBuilder::~SiPixelCondObjBuilder ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 103 of file SiPixelCondObjBuilder.cc.

References cond::service::PoolDBOutputService::appendOneIOV(), cond::service::PoolDBOutputService::beginOfTime(), calmap_, cond::service::PoolDBOutputService::createOneIOV(), cond::service::PoolDBOutputService::currentTime(), cppFunctionSkipper::exception, fromFile_, PedestalClient_cfi::gain, edm::EventSetup::getData(), mps_fire::i, iEvent, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, meanGain_, meanPed_, hgcalPlots::ncols, PixelTopology::ncolumns(), PixelTopology::nrows(), numberOfModules_, cms::SiPixelCondObjBuilder::CalParameters::p0, cms::SiPixelCondObjBuilder::CalParameters::p1, pddToken_, PixelIndices::pixelToChannelROC(), FastTimerService_cff::range, nano_mu_digi_cff::rawId, recordName_, rmsGain_, rmsPed_, writedatasetfile::run, secondRocRowGainOffset_, secondRocRowPedOffset_, SiPixelGainCalibration_, PixelGeomDetUnit::specificTopology(), PixelIndices::transformToROC(), and cms::Exception::what().

                                                                                          {
     using namespace edm;
     unsigned int run = iEvent.id().run();
     int nmodules = 0;
     uint32_t nchannels = 0;
     //    int mycol = 415;
     //    int myrow = 159;
 
     edm::LogInfo("SiPixelCondObjBuilder")
         << "... creating dummy SiPixelGainCalibration Data for Run " << run << "\n " << std::endl;
     //
     // Instantiate Gain calibration offset and define pedestal/gain range
     //
     float mingain = 0;
     float maxgain = 10;
     float minped = 0;
     float maxped = 255;
     SiPixelGainCalibration_ = std::make_unique<SiPixelGainCalibration>(minped, maxped, mingain, maxgain);
 
     const TrackerGeometry* pDD = &iSetup.getData(pddToken_);
     edm::LogInfo("SiPixelCondObjBuilder") << " There are " << pDD->dets().size() << " detectors" << std::endl;
 
     for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) {
       if (dynamic_cast<PixelGeomDetUnit const*>((*it)) != nullptr) {
         uint32_t detid = ((*it)->geographicalId()).rawId();
 
         // Stop if module limit reached
         nmodules++;
         if (nmodules > numberOfModules_)
           break;
 
         const PixelGeomDetUnit* pixDet = dynamic_cast<const PixelGeomDetUnit*>((*it));
         const PixelTopology& topol = pixDet->specificTopology();
         // Get the module sizes.
         int nrows = topol.nrows();     // rows in x
         int ncols = topol.ncolumns();  // cols in y
         //edm::LogPrint("SiPixelCondObjBuilder") << " ---> PIXEL DETID " << detid << " Cols " << ncols << " Rows " << nrows << std::endl;
 
         PixelIndices pIndexConverter(ncols, nrows);
 
         std::vector<char> theSiPixelGainCalibration;
 
         // Loop over columns and rows of this DetID
         for (int i = 0; i < ncols; i++) {
           for (int j = 0; j < nrows; j++) {
             nchannels++;
 
             float ped = 0.0, gain = 0.0;
 
             if (fromFile_) {
               // Use calibration from a file
               int chipIndex = 0, colROC = 0, rowROC = 0;
 
               pIndexConverter.transformToROC(i, j, chipIndex, colROC, rowROC);
               int chanROC = PixelIndices::pixelToChannelROC(rowROC, colROC);  // use ROC coordinates
               //         float pp0=0, pp1=0;
               std::map<int, CalParameters, std::less<int> >::const_iterator it = calmap_.find(chanROC);
               CalParameters theCalParameters = (*it).second;
               ped = theCalParameters.p0;
               gain = theCalParameters.p1;
 
             } else {
               if (rmsPed_ > 0) {
                 ped = CLHEP::RandGauss::shoot(meanPed_, rmsPed_);
                 while (minped > ped || maxped < ped)
                   ped = CLHEP::RandGauss::shoot(meanPed_, rmsPed_);
 
               } else
                 ped = meanPed_;
               if (rmsGain_ > 0) {
                 gain = CLHEP::RandGauss::shoot(meanGain_, rmsGain_);
                 while (mingain > gain || maxgain < gain)
                   gain = CLHEP::RandGauss::shoot(meanGain_, rmsGain_);
               } else
                 gain = meanGain_;
             }
 
             //if in the second row of rocs (i.e. a 2xN plaquette) add an offset (if desired) for testing
             if (j >= 80) {
               ped += secondRocRowPedOffset_;
               gain += secondRocRowGainOffset_;
 
               if (gain > maxgain)
                 gain = maxgain;
               else if (gain < mingain)
                 gain = mingain;
 
               if (ped > maxped)
                 ped = maxped;
               else if (ped < minped)
                 ped = minped;
             }
 
             //     if(i==mycol && j==myrow) {
             //       edm::LogPrint("SiPixelCondObjBuilder") << "       Col "<<i<<" Row "<<j<<" Ped "<<ped<<" Gain "<<gain<<std::endl;
             //     }
 
             //     gain =  2.8;
             //     ped  = 28.2;
 
             // Insert data in the container
             SiPixelGainCalibration_->setData(ped, gain, theSiPixelGainCalibration);
           }
         }
 
         SiPixelGainCalibration::Range range(theSiPixelGainCalibration.begin(), theSiPixelGainCalibration.end());
         if (!SiPixelGainCalibration_->put(detid, range, ncols))
           edm::LogError("SiPixelCondObjBuilder")
               << "[SiPixelCondObjBuilder::analyze] detid already exists" << std::endl;
       }
     }
     edm::LogPrint("SiPixelCondObjBuilder") << " ---> PIXEL Modules  " << nmodules << std::endl;
     edm::LogPrint("SiPixelCondObjBuilder") << " ---> PIXEL Channels " << nchannels << std::endl;
 
     //   // Try to read object
     //    int mynmodules =0;
     //    for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){
     //      if( dynamic_cast<PixelGeomDetUnit const*>((*it))!=0){
     //        uint32_t mydetid=((*it)->geographicalId()).rawId();
     //        mynmodules++;
     //        if( mynmodules > numberOfModules_) break;
     //        SiPixelGainCalibration::Range myrange = SiPixelGainCalibration_->getRange(mydetid);
     //        float mypedestal = SiPixelGainCalibration_->getPed (mycol,myrow,myrange,416);
     //        float mygain     = SiPixelGainCalibration_->getGain(mycol,myrow,myrange,416);
     //        //edm::LogPrint("SiPixelCondObjBuilder")<<" PEDESTAL "<< mypedestal<<" GAIN "<<mygain<<std::endl;
     //      }
     //    }
     // Write into DB
     edm::LogInfo(" --- writeing to DB!");
     edm::Service<cond::service::PoolDBOutputService> mydbservice;
     if (!mydbservice.isAvailable()) {
       edm::LogPrint("SiPixelCondObjBuilder") << "Didn't get DB" << std::endl;
       edm::LogError("db service unavailable");
       return;
     } else {
       edm::LogInfo("DB service OK");
     }
 
     try {
       //     size_t callbackToken = mydbservice->callbackToken("SiPixelGainCalibration");
       //     edm::LogInfo("SiPixelCondObjBuilder")<<"CallbackToken SiPixelGainCalibration "
       //         <<callbackToken<<std::endl;
       //       unsigned long long tillTime;
       //     if ( appendMode_)
       //     tillTime = mydbservice->currentTime();
       //       else
       //     tillTime = mydbservice->endOfTime();
       //     edm::LogInfo("SiPixelCondObjBuilder")<<"[SiPixelCondObjBuilder::analyze] tillTime = "
       //         <<tillTime<<std::endl;
       //     mydbservice->newValidityForNewPayload<SiPixelGainCalibration>(
       //           SiPixelGainCalibration_, tillTime , callbackToken);
 
       if (mydbservice->isNewTagRequest(recordName_)) {
         mydbservice->createOneIOV<SiPixelGainCalibration>(
             *SiPixelGainCalibration_, mydbservice->beginOfTime(), recordName_);
       } else {
         mydbservice->appendOneIOV<SiPixelGainCalibration>(
             *SiPixelGainCalibration_, mydbservice->currentTime(), recordName_);
       }
       edm::LogInfo(" --- all OK");
     } catch (const cond::Exception& er) {
       edm::LogPrint("SiPixelCondObjBuilder") << "Database exception!   " << er.what() << std::endl;
       edm::LogError("SiPixelCondObjBuilder") << er.what() << std::endl;
     } catch (const std::exception& er) {
       edm::LogPrint("SiPixelCondObjBuilder") << "Standard exception!   " << er.what() << std::endl;
       edm::LogError("SiPixelCondObjBuilder") << "caught std::exception " << er.what() << std::endl;
     } catch (...) {
       edm::LogError("SiPixelCondObjBuilder") << "Funny error" << std::endl;
     }
   }

◆ beginJob()

void SiPixelCondObjBuilder::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 275 of file SiPixelCondObjBuilder.cc.

References calmap_, fromFile_, and loadFromFile().

                                        {
     if (fromFile_) {
       if (loadFromFile()) {
         edm::LogInfo("SiPixelCondObjBuilder") << " Calibration loaded: Map size " << calmap_.size() << " max "
                                               << calmap_.max_size() << " " << calmap_.empty() << std::endl;
       }
     }
   }

◆ loadFromFile()

bool SiPixelCondObjBuilder::loadFromFile ( )

Definition at line 284 of file SiPixelCondObjBuilder.cc.

References calmap_, officialStyle::chan, fileName_, RemoveAddSevLevel::flag, mps_fire::i, recoMuon::in, runGCPTkAlMap::in_file, mps_splice::line, Skims_PA_cff::name, cms::SiPixelCondObjBuilder::CalParameters::p0, cms::SiPixelCondObjBuilder::CalParameters::p1, pfMETCorrectionType0_cfi::par0, pfMETCorrectionType0_cfi::par1, PixelIndices::pixelToChannelROC(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by beginJob(), and VarParsing.VarParsing::parseArguments().

                                            {
     float par0, par1;  //,par2,par3;
     int colid, rowid;  //rocid
     std::string name;
 
     std::ifstream in_file;                          // data file pointer
     in_file.open(fileName_.c_str(), std::ios::in);  // in C++
     if (in_file.bad()) {
       edm::LogError("SiPixelCondObjBuilder") << "Input file not found" << std::endl;
     }
     if (in_file.eof() != 0) {
       edm::LogError("SiPixelCondObjBuilder") << in_file.eof() << " " << in_file.gcount() << " " << in_file.fail() << " "
                                              << in_file.good() << " end of file " << std::endl;
       return false;
     }
     //Load file header
     char line[500];
     for (int i = 0; i < 3; i++) {
       in_file.getline(line, 500, '\n');
       edm::LogInfo("SiPixelCondObjBuilder") << line << std::endl;
     }
     //Loading calibration constants from file, loop on pixels
     for (int i = 0; i < (52 * 80); i++) {
       in_file >> par0 >> par1 >> name >> colid >> rowid;
 
       edm::LogPrint("SiPixelCondObjBuilder")
           << " Col " << colid << " Row " << rowid << " P0 " << par0 << " P1 " << par1 << std::endl;
 
       CalParameters onePix;
       onePix.p0 = par0;
       onePix.p1 = par1;
 
       // Convert ROC pixel index to channel
       int chan = PixelIndices::pixelToChannelROC(rowid, colid);
       calmap_.insert(std::pair<int, CalParameters>(chan, onePix));
     }
 
     bool flag;
     if (calmap_.size() == 4160) {
       flag = true;
     } else {
       flag = false;
     }
     return flag;
   }

Member Data Documentation

◆ appendMode_

const bool cms::SiPixelCondObjBuilder::appendMode_

private

Definition at line 55 of file SiPixelCondObjBuilder.cc.

◆ calmap_

std::map<int, CalParameters, std::less<int> > cms::SiPixelCondObjBuilder::calmap_

private

Definition at line 78 of file SiPixelCondObjBuilder.cc.

Referenced by analyze(), beginJob(), and loadFromFile().

◆ fileName_

const std::string cms::SiPixelCondObjBuilder::fileName_

private

Definition at line 69 of file SiPixelCondObjBuilder.cc.

Referenced by loadFromFile().

◆ fromFile_

const bool cms::SiPixelCondObjBuilder::fromFile_

private

Definition at line 68 of file SiPixelCondObjBuilder.cc.

Referenced by analyze(), and beginJob().

◆ meanGain_

const double cms::SiPixelCondObjBuilder::meanGain_

private

Definition at line 63 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ meanPed_

const double cms::SiPixelCondObjBuilder::meanPed_

private

Definition at line 61 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ numberOfModules_

const int cms::SiPixelCondObjBuilder::numberOfModules_

private

Definition at line 67 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ pddToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> cms::SiPixelCondObjBuilder::pddToken_

private

Definition at line 56 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ pIndexConverter_

PixelIndices* cms::SiPixelCondObjBuilder::pIndexConverter_

private

Definition at line 79 of file SiPixelCondObjBuilder.cc.

◆ recordName_

const std::string cms::SiPixelCondObjBuilder::recordName_

private

Definition at line 59 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ rmsGain_

const double cms::SiPixelCondObjBuilder::rmsGain_

private

Definition at line 64 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ rmsPed_

const double cms::SiPixelCondObjBuilder::rmsPed_

private

Definition at line 62 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ secondRocRowGainOffset_

const double cms::SiPixelCondObjBuilder::secondRocRowGainOffset_

private

Definition at line 65 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ secondRocRowPedOffset_

const double cms::SiPixelCondObjBuilder::secondRocRowPedOffset_

private

Definition at line 66 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ SiPixelGainCalibration_

std::unique_ptr<SiPixelGainCalibration> cms::SiPixelCondObjBuilder::SiPixelGainCalibration_

private

Definition at line 57 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

◆ SiPixelGainCalibrationService_

SiPixelGainCalibrationService cms::SiPixelCondObjBuilder::SiPixelGainCalibrationService_

private

Definition at line 58 of file SiPixelCondObjBuilder.cc.

Classes

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiPixelCondObjBuilder()

◆ ~SiPixelCondObjBuilder()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ loadFromFile()

Member Data Documentation

◆ appendMode_

◆ calmap_

◆ fileName_

◆ fromFile_

◆ meanGain_

◆ meanPed_

◆ numberOfModules_

◆ pddToken_

◆ pIndexConverter_

◆ recordName_

◆ rmsGain_

◆ rmsPed_

◆ secondRocRowGainOffset_

◆ secondRocRowPedOffset_

◆ SiPixelGainCalibration_

◆ SiPixelGainCalibrationService_