#include <EcalLaserCondTools.h>

Inheritance diagram for EcalLaserCondTools:

Classes
class	CorrReader

class	FileReader

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

	EcalLaserCondTools (const edm::ParameterSet &)

void	from_hdf_to_db ()

	~EcalLaserCondTools () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

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

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 Member Functions
void	dbToAscii (const edm::EventSetup &es)

void	fillDb (CorrReader &r)

void	processIov (CorrReader &r, int t1, int t2[nLmes], int t3)

Static Private Member Functions
static std::string	timeToString (time_t t)

static std::string	toNth (int n)

Private Attributes
edm::Service< cond::service::PoolDBOutputService >	db_

FILE *	eventList_

std::string	eventListFileName_

FILE *	ferr_

std::vector< std::string >	fnames_

FILE *	fout_

int	fromTime_

double	maxP_

double	minP_

std::string	mode_

int	nIovs_

int	skipIov_

int	toTime_

int	verb_

Static Private Attributes
static constexpr int	nLmes = 92

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 23 of file EcalLaserCondTools.h.

Constructor & Destructor Documentation

◆ EcalLaserCondTools()

EcalLaserCondTools::EcalLaserCondTools ( const edm::ParameterSet & ps )

Constructor

Parameters

ps	analyser configuration

Definition at line 26 of file EcalLaserCondTools.cc.

     : fout_(nullptr),
       eventList_(nullptr),
       eventListFileName_(ps.getParameter<std::string>("eventListFile")),
       verb_(ps.getParameter<int>("verbosity")),
       mode_(ps.getParameter<std::string>("mode")),
       fnames_(ps.getParameter<std::vector<std::string> >("inputFiles")),
       skipIov_(ps.getParameter<int>("skipIov")),
       nIovs_(ps.getParameter<int>("nIovs")),
       fromTime_(ps.getParameter<int>("fromTime")),
       toTime_(ps.getParameter<int>("toTime")),
       minP_(ps.getParameter<double>("transparencyMin")),
       maxP_(ps.getParameter<double>("transparencyMax")) {
   ferr_ = fopen("corr_errors.txt", "w");
   fprintf(ferr_, "#t1\tdetid\tp1\tp2\tp3");
  
   if (!eventListFileName_.empty()) {
     eventList_ = fopen(eventListFileName_.c_str(), "r");
     if (eventList_ == nullptr)
       throw cms::Exception("User") << "Failed to open file " << eventListFileName_ << "\n";
   }
 }

References eventList_, eventListFileName_, Exception, and ferr_.

◆ ~EcalLaserCondTools()

EcalLaserCondTools::~EcalLaserCondTools ( )

override

Destructor

Definition at line 49 of file EcalLaserCondTools.cc.

                                         {
   if (ferr_)
     fclose(ferr_);
   if (fout_)
     fclose(fout_);
 }

References ferr_, and fout_.

Member Function Documentation

◆ analyze()

void EcalLaserCondTools::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

overridevirtual

Called by CMSSW event loop

Parameters

evt	the event
es	events setup

Implements edm::EDAnalyzer.

Definition at line 56 of file EcalLaserCondTools.cc.

                                                                                {
   static bool done = false;
  
   if (done && (mode_ == "ascii_file_to_db" || mode_ == "hdf_file_to_db")) {
     return;
   }
  
   if (mode_ == "ascii_file_to_db") {
     if (verb_ > 2)
       std::cout << "ascii_file_to_db mode\n";
  
     if (!db_.isAvailable()) {
       throw cms::Exception("CondDBAccess") << "Failed to connect to PoolDBOutputService\n";
     }
     FileReader corrReader(fnames_);
     corrReader.setVerbosity(verb_);
     fillDb(corrReader);
   } else if (mode_ == "hdf_file_to_db") {
     from_hdf_to_db();
   } else if (mode_ == "db_to_ascii_file") {
     dbToAscii(es);
   } else {
     cms::Exception("InvalidParam") << "Value of parameter mode is not valid. Expecting ascii_file_to_db or read";
   }
 }

References gather_cfg::cout, db_, dbToAscii(), fileCollector::done, Exception, fillDb(), fnames_, from_hdf_to_db(), edm::Service< T >::isAvailable(), mode_, EcalLaserCondTools::CorrReader::setVerbosity(), and verb_.

◆ dbToAscii()

void EcalLaserCondTools::dbToAscii ( const edm::EventSetup & es )

private

Definition at line 620 of file EcalLaserCondTools.cc.

                                                           {
   edm::ESHandle<EcalLaserAPDPNRatios> hCorr;
   es.get<EcalLaserAPDPNRatiosRcd>().get(hCorr);
  
   const EcalLaserAPDPNRatios* corr = hCorr.product();
  
   const EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap& p = corr->getLaserMap();
   const EcalLaserAPDPNRatios::EcalLaserTimeStampMap& t = corr->getTimeMap();
  
   if (t.size() != EcalLaserCondTools::nLmes)
     throw cms::Exception("LasCor") << "Unexpected number time parameter triplets\n";
  
   if (fout_ == nullptr) {
     fout_ = fopen("corr_dump.txt", "w");
     if (fout_ == nullptr)
       throw cms::Exception("LasCor") << "Failed to create file corr_dump.txt\n";
   }
  
   unsigned t1 = t[0].t1.unixTime();
   unsigned t3 = t[0].t3.unixTime();
   fprintf(fout_, "T %d\t%d", t1, t3);
  
   if (verb_)
     std::cout << "Processing IOV " << t1 << " - " << t3 << "(" << timeToString(t1) << " - " << timeToString(t3) << "\n";
  
   for (unsigned i = 0; i < t.size(); ++i) {
     if (t[i].t1.unixTime() != t1 || t[i].t3.unixTime() != t3) {
       throw cms::Exception("LasCor") << "Inconsitency in t1, t3: "
                                      << "t1(lme 1) =" << t1 << ", t1(lme " << (i + 1) << ") = " << t[i].t1.unixTime()
                                      << ", t3(lme 1) =" << t3 << ", t3(lme " << (i + 1) << ") = " << t[i].t3.unixTime()
                                      << "\n";
     }
     fprintf(fout_, "\t%d", t[i].t2.unixTime());
   }
   fputc('\n', fout_);
   fflush(fout_);
  
   for (int ieta = -EBDetId::MAX_IETA; ieta <= EBDetId::MAX_IETA; ++ieta) {
     if (ieta == 0)
       continue;
     for (int iphi = EBDetId::MIN_IPHI; iphi <= EBDetId::MAX_IPHI; ++iphi) {
       if (EBDetId::validDetId(ieta, iphi)) {
         EBDetId detId(ieta, iphi);
         EcalLaserAPDPNRatios::EcalLaserAPDPNpair corr = p.barrel(detId.hashedIndex());
         fprintf(fout_, "P %d\t%f\t%f\t%f\n", (int)detId, corr.p1, corr.p2, corr.p3);
       }
     }
   }
  
   for (int iZ = 1; iZ >= -1; --iZ) {
     for (int iX = EEDetId::IX_MIN; iX <= EEDetId::IX_MAX; ++iX) {
       for (int iY = EEDetId::IY_MIN; iY <= EEDetId::IY_MAX; ++iY) {
         if (EEDetId::validDetId(iX, iY, iZ)) {
           EEDetId detId(iX, iY, iZ);
           EcalLaserAPDPNRatios::EcalLaserAPDPNpair corr = p.endcap(detId.hashedIndex());
           fprintf(fout_, "P %d\t%f\t%f\t%f\n", (int)detId, corr.p1, corr.p2, corr.p3);
         }
       }
     }
   }
 }

References alignCSCRings::corr, gather_cfg::cout, Exception, fout_, edm::EventSetup::get(), get, EBDetId::hashedIndex(), EEDetId::hashedIndex(), mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, nLmes, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), submitPVValidationJobs::t, RandomServiceHelper::t1, RandomServiceHelper::t2, RandomServiceHelper::t3, timeToString(), EBDetId::validDetId(), EEDetId::validDetId(), and verb_.

Referenced by analyze().

◆ fillDb()

void EcalLaserCondTools::fillDb ( CorrReader & r )

private

Definition at line 261 of file EcalLaserCondTools.cc.

                                              {
   int iIov = 0;
   int processedIovs = 0;
   if (verb_ > 2)
     std::cout << "Starting filling DB...\n";
   int t1 = 0;
   int t3 = 0;
   int t2[nLmes];
  
   int prevT1 = 0;
   int prevT3 = 0;
   int prevT = 0;
   int t = 0;
   if (eventList_) {
     int iline = 0;
     while (!feof(eventList_)) {
       //skips comment lines:
       char c;
       while (fscanf(eventList_, " %1[#]%*[^\n]\n", &c) == 1)
         ++iline;
  
       int n = fscanf(eventList_, "%*d %*d %*d %d%*[^\n]\n", &t);
       if (verb_ > 1)
         std::cout << "Event time: t = " << t << ", " << timeToString(t) << "\n";
       ++iline;
       if (n != 1)
         throw cms::Exception("User") << "Syntax error in event list file " << eventListFileName_ << " at line " << iline
                                      << " " << n << " "
                                      << ".\n";
       if (t < prevT)
         throw cms::Exception("User") << "Events in the event list file " << eventListFileName_
                                      << " are not ordered in increased time as required! See line " << iline << "\n";
       if (t == 0)
         throw cms::Exception("User") << "Found an unexpected t = 0 time value "
                                         "in the event list file"
                                      << eventListFileName_ << " at line " << iline << "!\n";
       //Look for IOV:
       bool iovFound = true;
       if (t <= t3) {  //IOV already inserted for previous event.
         if (verb_ > 1)
           std::cout << "Event in same IOV than previous one.\n";
         continue;
       }
  
       while ((iovFound = r.readTime(t1, t2, t3)) && t3 < t) /*NOP*/
         ;
  
       if (iovFound) {
         if (t1 < prevT1 && t3 < prevT3)
           throw cms::Exception("User")
               << "IOVs in the correction ascii file are not ordered in increased time as required!\n";
         else if (t1 < prevT1 || t3 < prevT3)
           throw cms::Exception("User") << "Found interleaved IOVs in the correction ascii file!\n";
         processIov(r, t1, t2, t3);
       } else {
         std::cout << "Warning: event beyond last IOV t3. Event time: " << timeToString(t)
                   << ". Last IOV t3: " << timeToString(t3) << "\n";
       }
     }
   } else
     while (r.readTime(t1, t2, t3)) {
       ++iIov;
       if (iIov <= skipIov_) {
         std::cout << "Skipping IOV " << iIov << "\n";
         continue;
       } else if (processedIovs >= nIovs_) {
         std::cout << "Requested number of IOVs, " << nIovs_ << ", processed.\n";
         return;
       } else {
         processIov(r, t1, t2, t3);
         ++processedIovs;
       }
     }
 }

References c, gather_cfg::cout, eventList_, eventListFileName_, Exception, ntuplemaker::iline, dqmiodumpmetadata::n, nIovs_, nLmes, processIov(), alignCSCRings::r, skipIov_, submitPVValidationJobs::t, RandomServiceHelper::t1, RandomServiceHelper::t2, RandomServiceHelper::t3, timeToString(), and verb_.

Referenced by analyze().

◆ from_hdf_to_db()

void EcalLaserCondTools::from_hdf_to_db ( )

Definition at line 82 of file EcalLaserCondTools.cc.

                                         {
   cond::Time_t iovStart = 0;
  
   hid_t file, space, memspace;
   hid_t dset_rawid, dset_t2, dset;
  
   hsize_t dims[2] = {};
  
   for (unsigned int ifile = 0; ifile < fnames_.size(); ++ifile) {
     if (verb_) {
       std::cout << " - converting file: " << fnames_[ifile] << "\n";
     }
  
     file = H5Fopen(fnames_[ifile].c_str(), H5F_ACC_RDONLY, H5P_DEFAULT);
  
     dset_rawid = H5Dopen(file, "cmssw_id", H5P_DEFAULT);
     space = H5Dget_space(dset_rawid);
     H5Sget_simple_extent_dims(space, dims, nullptr);
  
     unsigned int nCrystals = dims[0];
     int rawid[nCrystals];
     herr_t status;
  
     status = H5Dread(dset_rawid, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rawid);
     if (status < 0)
       throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
     H5Dclose(dset_rawid);
     H5Sclose(space);
  
     dset_t2 = H5Dopen(file, "t2", H5P_DEFAULT);
     space = H5Dget_space(dset_t2);
     H5Sget_simple_extent_dims(space, dims, nullptr);
  
     unsigned int nIovs = dims[0];
     unsigned int nLME = dims[1];
  
     if (verb_) {
       std::cout << "Number of crystals: " << nCrystals << "\n";
       std::cout << "Number of IOVs: " << nIovs << "\n";
       std::cout << "Number of Monitoring regions: " << nLME << "\n";
     }
  
     int t1[nIovs], t3[nIovs], t2[nIovs][nLME];
  
     // -- reading data (cmsswid, t2, t1, t3, p2, p1, p3
     if (verb_ > 1)
       std::cout << " * reading t2 table "
                 << "\n";
     status = H5Dread(dset_t2, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, t2[0]);
     if (status < 0)
       throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
     H5Dclose(dset_t2);
     //H5Sclose(space);
  
     if (verb_ > 1)
       std::cout << " * reading t1 table "
                 << "\n";
     dset = H5Dopen(file, "t1", H5P_DEFAULT);
     status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, t1);
     if (status < 0)
       throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
     H5Dclose(dset);
  
     if (verb_ > 1)
       std::cout << " * reading t3 table "
                 << "\n";
     dset = H5Dopen(file, "t3", H5P_DEFAULT);
     status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, t3);
     if (status < 0)
       throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
     H5Dclose(dset);
  
     assert(EcalLaserCondTools::nLmes == nLME);
  
     // read crystal info IOV by IOV (otherwise too large)
     float p1[nCrystals], p2[nCrystals], p3[nCrystals];
     hsize_t iov_dim[1] = {nCrystals};
     memspace = H5Screate_simple(1, iov_dim, nullptr);
  
     auto corrSet = std::make_unique<EcalLaserAPDPNRatios>();
     for (unsigned int iIov = skipIov_; iIov < nIovs && iIov < unsigned(nIovs_); ++iIov) {
       EcalLaserAPDPNRatios::EcalLaserTimeStamp t;
       iovStart = uint64_t(t1[iIov]) << 32;
       for (size_t iLme = 0; iLme < EcalLaserCondTools::nLmes; ++iLme) {
         t.t1 = edm::Timestamp(uint64_t(t1[iIov]) << 32);
         t.t2 = edm::Timestamp(uint64_t(t2[iIov][iLme]) << 32);
         t.t3 = edm::Timestamp(uint64_t(t3[iIov]) << 32);
         corrSet->setTime(iLme, t);
       }
  
       hsize_t offset[2] = {iIov, 0};      // shift rows: iIov, columns: 0
       hsize_t count[2] = {1, nCrystals};  // 1 row, nXtal columns
  
       dset = H5Dopen(file, "p1", H5P_DEFAULT);
       space = H5Dget_space(dset);
       status = H5Sselect_hyperslab(space, H5S_SELECT_SET, offset, nullptr, count, nullptr);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
       status = H5Dread(dset, H5T_NATIVE_FLOAT, memspace, space, H5P_DEFAULT, p1);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
       H5Dclose(dset);
       //H5Sclose(space);
  
       dset = H5Dopen(file, "p2", H5P_DEFAULT);
       space = H5Dget_space(dset);
       status = H5Sselect_hyperslab(space, H5S_SELECT_SET, offset, nullptr, count, nullptr);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
       status = H5Dread(dset, H5T_NATIVE_FLOAT, memspace, space, H5P_DEFAULT, p2);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
       H5Dclose(dset);
       //      H5Sclose(space);
  
       dset = H5Dopen(file, "p3", H5P_DEFAULT);
       space = H5Dget_space(dset);
       status = H5Sselect_hyperslab(space, H5S_SELECT_SET, offset, nullptr, count, nullptr);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
  
       status = H5Dread(dset, H5T_NATIVE_FLOAT, memspace, space, H5P_DEFAULT, p3);
       if (status < 0)
         throw cms::Exception("EcalLaserCondTool:HDF") << "Error while reading HD file.";
       H5Dclose(dset);
       H5Sclose(space);
  
       for (size_t iXtal = 0; iXtal < nCrystals; ++iXtal) {
         DetId detid = rawid[iXtal];
  
         EcalLaserAPDPNRatios::EcalLaserAPDPNpair corr = EcalLaserAPDPNRatios::EcalLaserAPDPNpair();
         corr.p1 = p1[iXtal];
         corr.p2 = p2[iXtal];
         corr.p3 = p3[iXtal];
  
         if (!isfinite(corr.p1) || !isfinite(corr.p2) || !isfinite(corr.p3) || corr.p1 < minP_ || corr.p1 > maxP_ ||
             corr.p2 < minP_ || corr.p2 > maxP_ || corr.p3 < minP_ || corr.p3 > maxP_) {
           fprintf(ferr_, "%d %d %f %f %f\n", t1[iIov], (int)detid, corr.p1, corr.p2, corr.p3);
           corr.p1 = corr.p2 = corr.p3 = 1;
         }
         corrSet->setValue((int)detid, corr);
       }
  
       try {
         //Write correction set in DB (one IOV):
         //if (db_->isNewTagRequest("EcalLaserAPDPNRatiosRcd")) {
         //  if (verb_)
         //    std::cout << "First IOV, extending starting time.\n";
         //  iovStart = db_->beginOfTime();
         //}
         timeval t;
         gettimeofday(&t, nullptr);
         if (verb_ > 1)
           std::cout << "[" << timeToString(t.tv_sec) << "] "
                     << "Write IOV " << iIov << " starting from " << timeToString(iovStart >> 32) << "... ";
         db_->writeOne(corrSet.get(), iovStart, "EcalLaserAPDPNRatiosRcd");
       } catch (const cms::Exception& e) {
         if (verb_ > 1)
           std::cout << "Failed. ";
         std::cout << "Exception catched while writting to cond DB" << e.what() << "\n";
       }
       if (verb_ > 1)
         std::cout << "Suceeded.\n";
  
     }  // loop over IOVs
  
     H5Sclose(memspace);
     H5Fclose(file);
   }  // loop over input files
 }

References cms::cuda::assert(), alignCSCRings::corr, submitPVResolutionJobs::count, gather_cfg::cout, db_, MillePedeFileConverter_cfg::e, Exception, ferr_, geometryDiff::file, fnames_, compare_using_db::ifile, maxP_, minP_, ecaldqm::nCrystals(), nIovs_, nLmes, hltrates_dqm_sourceclient-live_cfg::offset, p1, p2, p3, skipIov_, mps_update::status, submitPVValidationJobs::t, RandomServiceHelper::t1, RandomServiceHelper::t2, RandomServiceHelper::t3, timeToString(), verb_, and cond::service::PoolDBOutputService::writeOne().

Referenced by analyze().

◆ processIov()

void EcalLaserCondTools::processIov	(	CorrReader &	r,
		int	t1,
		int	t2[nLmes],
		int	t3
	)

private

Definition at line 336 of file EcalLaserCondTools.cc.

                                                                                                   {
   static int iIov = 0;
   ++iIov;
  
   //   FILE* fdebug = 0;
   //   if(i==) fdebug = fopen("debug.txt", "w");
  
   //  if(iIov <= skipIov_) { std::cout << "Skipping IOV " << iIov << "\n"; return; }
  
   cond::Time_t iovStart = 0;
  
   if (verb_ > 1) {
     std::cout << "t1:" << t1 << "(" << timeToString(t1) << ") \n"
               << "t3: " << t3 << "(" << timeToString(t3) << ")\nt2-t1: ";
     for (int i = 0; i < EcalLaserCondTools::nLmes; ++i)
       std::cout << t2[i] - t1 << "\t";
  
     std::cout << "\n";
   }
   if (t1 < fromTime_) {
     std::cout << "Skipping IOV " << iIov << ", "
               << ", which is before 'fromTime'," << timeToString(fromTime_) << "(" << fromTime_ << ").\n";
     return;
   }
  
   if (toTime_ != -1 && t3 < toTime_) {
     std::cout << "Skipping IOV " << iIov << ", "
               << ", which is beyond 'toTime'," << timeToString(toTime_) << "(" << toTime_ << ").\n";
     return;
   }
  
   if (t1 == 0) {
     std::cout << "Skipping IOV with t1 = 0"
               << "\n";
     return;
   }
  
   auto corrSet = std::make_unique<EcalLaserAPDPNRatios>();
  
   EcalLaserAPDPNRatios::EcalLaserTimeStamp t;
   iovStart = uint64_t(t1) << 32;
   for (size_t i = 0; i < EcalLaserCondTools::nLmes; ++i) {
     t.t1 = edm::Timestamp(uint64_t(t1) << 32);
     t.t2 = edm::Timestamp(uint64_t(t2[i]) << 32);
     t.t3 = edm::Timestamp(uint64_t(t3) << 32);
     corrSet->setTime(i, t);
   }
  
   constexpr int ncrystals = 75848;
   std::set<int> detidList;
   for (int i = 0; i < ncrystals; ++i) {
     DetId detid;
     //EcalLaserAPDPNRatios::EcalLaserAPDPNpair corr = {0, 0, 0};
     EcalLaserAPDPNRatios::EcalLaserAPDPNpair corr = EcalLaserAPDPNRatios::EcalLaserAPDPNpair();
     if (verb_ > 2)
       std::cout << "Reading " << toNth(i + 1) << " crystal\n";
     if (!r.readPs(detid, corr)) {
       throw cms::Exception("LasCor") << "Failed to read " << toNth(i + 1) << " crystal correction.\n";
     }
  
     std::pair<std::set<int>::iterator, bool> res = detidList.insert(int(detid));
  
     if (!res.second) {  //detid already processed
       std::cout << "Duplicate det id, for IOV " << iIov << " t1 = " << t1 << " detid = " << int(detid) << "\n";
     }
  
     if (!isfinite(corr.p1) || !isfinite(corr.p2) || !isfinite(corr.p3) || corr.p1 < minP_ || corr.p1 > maxP_ ||
         corr.p2 < minP_ || corr.p2 > maxP_ || corr.p3 < minP_ || corr.p3 > maxP_) {
       fprintf(ferr_, "%d %d %f %f %f\n", t1, (int)detid, corr.p1, corr.p2, corr.p3);
       corr.p1 = corr.p2 = corr.p3 = 1;
     }
  
     if (verb_ > 2) {
       if (detid.subdetId() == EcalBarrel) {
         std::cout << EBDetId(detid);
       } else if (detid.subdetId() == EcalEndcap) {
         std::cout << EEDetId(detid);
       } else {
         std::cout << (int)detid;
       }
       std::cout << ": "
                 << "p1 = " << corr.p1 << "\t"
                 << "p2 = " << corr.p2 << "\t"
                 << "p3 = " << corr.p3 << "\n";
     }
  
     corrSet->setValue((int)detid, corr);
   }
  
   try {
     //Write correction set in DB (one IOV):
     if (db_->isNewTagRequest("EcalLaserAPDPNRatiosRcd")) {
       if (verb_)
         std::cout << "First IOV, extending starting time.\n";
  
       iovStart = db_->beginOfTime();
     }
     timeval t;
     gettimeofday(&t, nullptr);
     if (verb_ > 1)
       std::cout << "[" << timeToString(t.tv_sec) << "] "
                 << "Write IOV " << iIov << " starting from " << timeToString(iovStart >> 32) << "... ";
     db_->writeOne(corrSet.get(), iovStart, "EcalLaserAPDPNRatiosRcd");
   } catch (const cms::Exception& e) {
     std::cout << "Failed.\nException cathed while writting to cond DB" << e.what() << "\n";
   }
   std::cout << "Suceeded.\n";
 }

References cond::service::PoolDBOutputService::beginOfTime(), alignCSCRings::corr, gather_cfg::cout, db_, MillePedeFileConverter_cfg::e, EcalBarrel, EcalEndcap, Exception, ferr_, fromTime_, mps_fire::i, createfilelist::int, cond::service::PoolDBOutputService::isNewTagRequest(), maxP_, minP_, nLmes, alignCSCRings::r, DetId::subdetId(), submitPVValidationJobs::t, RandomServiceHelper::t1, RandomServiceHelper::t2, RandomServiceHelper::t3, timeToString(), toNth(), toTime_, verb_, and cond::service::PoolDBOutputService::writeOne().

Referenced by fillDb().

◆ timeToString()

std::string EcalLaserCondTools::timeToString ( time_t t )

staticprivate

Definition at line 611 of file EcalLaserCondTools.cc.

                                                    {
   char buf[256];
   struct tm lt;
   localtime_r(&t, &lt);
   strftime(buf, sizeof(buf), "%F %R:%S", &lt);
   buf[sizeof(buf) - 1] = 0;
   return std::string(buf);
 }

References visDQMUpload::buf, AlCaHLTBitMon_QueryRunRegistry::string, and submitPVValidationJobs::t.

Referenced by dbToAscii(), fillDb(), from_hdf_to_db(), and processIov().

◆ toNth()

std::string EcalLaserCondTools::toNth ( int n )

staticprivate

Definition at line 588 of file EcalLaserCondTools.cc.

                                          {
   std::stringstream s;
   s << n;
   if (n % 100 < 10 || n % 100 > 20) {
     switch (n % 10) {
       case 1:
         s << "st";
         break;
       case 2:
         s << "nd";
         break;
       case 3:
         s << "rd";
         break;
       default:
         s << "th";
     }
   } else {
     s << "th";
   }
   return s.str();
 }

References dqmiodumpmetadata::n, and alignCSCRings::s.

Referenced by processIov().

Member Data Documentation

◆ db_

edm::Service<cond::service::PoolDBOutputService> EcalLaserCondTools::db_

private

Definition at line 92 of file EcalLaserCondTools.h.

Referenced by analyze(), from_hdf_to_db(), and processIov().

◆ eventList_

FILE* EcalLaserCondTools::eventList_

private

Definition at line 87 of file EcalLaserCondTools.h.

Referenced by EcalLaserCondTools(), and fillDb().

◆ eventListFileName_

std::string EcalLaserCondTools::eventListFileName_

private

Definition at line 88 of file EcalLaserCondTools.h.

Referenced by EcalLaserCondTools(), and fillDb().

◆ ferr_

FILE* EcalLaserCondTools::ferr_

private

Definition at line 98 of file EcalLaserCondTools.h.

Referenced by EcalLaserCondTools(), from_hdf_to_db(), processIov(), and ~EcalLaserCondTools().

◆ fnames_

std::vector<std::string> EcalLaserCondTools::fnames_

private

Definition at line 91 of file EcalLaserCondTools.h.

Referenced by analyze(), from_hdf_to_db(), EcalLaserCondTools::FileReader::readPs(), and EcalLaserCondTools::FileReader::readTime().

◆ fout_

FILE* EcalLaserCondTools::fout_

private

Definition at line 86 of file EcalLaserCondTools.h.

Referenced by dbToAscii(), and ~EcalLaserCondTools().

◆ fromTime_

int EcalLaserCondTools::fromTime_

private

Definition at line 95 of file EcalLaserCondTools.h.

Referenced by processIov().

◆ maxP_

double EcalLaserCondTools::maxP_

private

Definition at line 97 of file EcalLaserCondTools.h.

Referenced by from_hdf_to_db(), and processIov().

◆ minP_

double EcalLaserCondTools::minP_

private

Definition at line 97 of file EcalLaserCondTools.h.

Referenced by from_hdf_to_db(), and processIov().

◆ mode_

std::string EcalLaserCondTools::mode_

private

Definition at line 90 of file EcalLaserCondTools.h.

Referenced by analyze().

◆ nIovs_

int EcalLaserCondTools::nIovs_

private

Definition at line 94 of file EcalLaserCondTools.h.

Referenced by fillDb(), and from_hdf_to_db().

◆ nLmes

constexpr int EcalLaserCondTools::nLmes = 92

staticconstexprprivate

Number of extended laser monitoring regions

Definition at line 28 of file EcalLaserCondTools.h.

Referenced by dbToAscii(), fillDb(), from_hdf_to_db(), processIov(), and EcalLaserCondTools::FileReader::readTime().

◆ skipIov_

int EcalLaserCondTools::skipIov_

private

Definition at line 93 of file EcalLaserCondTools.h.

Referenced by fillDb(), and from_hdf_to_db().

◆ toTime_

int EcalLaserCondTools::toTime_

private

Definition at line 96 of file EcalLaserCondTools.h.

Referenced by processIov().

◆ verb_

int EcalLaserCondTools::verb_

private

Definition at line 89 of file EcalLaserCondTools.h.

Referenced by analyze(), dbToAscii(), fillDb(), from_hdf_to_db(), processIov(), EcalLaserCondTools::FileReader::readPs(), and EcalLaserCondTools::FileReader::readTime().

Classes

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalLaserCondTools()

◆ ~EcalLaserCondTools()

Member Function Documentation

◆ analyze()

◆ dbToAscii()

◆ fillDb()

◆ from_hdf_to_db()

◆ processIov()

◆ timeToString()

◆ toNth()

Member Data Documentation

◆ db_

◆ eventList_

◆ eventListFileName_

◆ ferr_

◆ fnames_

◆ fout_

◆ fromTime_

◆ maxP_

◆ minP_

◆ mode_

◆ nIovs_

◆ nLmes

◆ skipIov_

◆ toTime_

◆ verb_