cms::SiPixelCondObjBuilder Class Reference

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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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 ( 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");
  }

cms::SiPixelCondObjBuilder::~SiPixelCondObjBuilder ( )

overridedefault

Member Function Documentation

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 calmap_, cppFunctionSkipper::exception, fromFile_, edm::EventSetup::getData(), mps_fire::i, edm::EventBase::id(), dqmiolumiharvest::j, meanGain_, meanPed_, PixelTopology::ncolumns(), PixelTopology::nrows(), numberOfModules_, cms::SiPixelCondObjBuilder::CalParameters::p0, cms::SiPixelCondObjBuilder::CalParameters::p1, pddToken_, PixelIndices::pixelToChannelROC(), sistrip::SpyUtilities::range(), recordName_, rmsGain_, rmsPed_, edm::EventID::run(), submitPVValidationJobs::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;
    }
  }

void SiPixelCondObjBuilder::beginJob ( void  )

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;
      }
    }
  }

bool SiPixelCondObjBuilder::loadFromFile

(

)

Definition at line 284 of file SiPixelCondObjBuilder.cc.

References calmap_, officialStyle::chan, fileName_, mps_fire::i, recoMuon::in, geometryCSVtoXML::line, mergeVDriftHistosByStation::name, cms::SiPixelCondObjBuilder::CalParameters::p0, cms::SiPixelCondObjBuilder::CalParameters::p1, 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

const bool cms::SiPixelCondObjBuilder::appendMode_

private

Definition at line 55 of file SiPixelCondObjBuilder.cc.

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

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

private

Definition at line 69 of file SiPixelCondObjBuilder.cc.

Referenced by loadFromFile().

const bool cms::SiPixelCondObjBuilder::fromFile_

private

Definition at line 68 of file SiPixelCondObjBuilder.cc.

Referenced by analyze(), and beginJob().

const double cms::SiPixelCondObjBuilder::meanGain_

private

Definition at line 63 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const double cms::SiPixelCondObjBuilder::meanPed_

private

Definition at line 61 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const int cms::SiPixelCondObjBuilder::numberOfModules_

private

Definition at line 67 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

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

private

Definition at line 56 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

PixelIndices* cms::SiPixelCondObjBuilder::pIndexConverter_

private

Definition at line 79 of file SiPixelCondObjBuilder.cc.

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

private

Definition at line 59 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const double cms::SiPixelCondObjBuilder::rmsGain_

private

Definition at line 64 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const double cms::SiPixelCondObjBuilder::rmsPed_

private

Definition at line 62 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const double cms::SiPixelCondObjBuilder::secondRocRowGainOffset_

private

Definition at line 65 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

const double cms::SiPixelCondObjBuilder::secondRocRowPedOffset_

private

Definition at line 66 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

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

private

Definition at line 57 of file SiPixelCondObjBuilder.cc.

Referenced by analyze().

SiPixelGainCalibrationService cms::SiPixelCondObjBuilder::SiPixelGainCalibrationService_

private

Definition at line 58 of file SiPixelCondObjBuilder.cc.