#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::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 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_

edm::ESGetToken < EcalLaserAPDPNRatios, EcalLaserAPDPNRatiosRcd >	laserAPDPNRatiosToken_

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::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 24 of file EcalLaserCondTools.h.

Constructor & Destructor Documentation

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

Constructor

Parameters

ps	analyser configuration

Definition at line 25 of file EcalLaserCondTools.cc.

References edm::EDConsumerBase::esConsumes(), eventList_, eventListFileName_, Exception, ferr_, laserAPDPNRatiosToken_, and mode_.

     : 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")) {
   if (mode_ == "db_to_ascii_file") {
     laserAPDPNRatiosToken_ = esConsumes();
   }
 
   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";
   }
 }

EcalLaserCondTools::~EcalLaserCondTools ( )

override

Destructor

Definition at line 52 of file EcalLaserCondTools.cc.

References ferr_, and fout_.

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

Member Function Documentation

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::one::EDAnalyzerBase.

Definition at line 59 of file EcalLaserCondTools.cc.

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

                                                                                {
   if (mode_ == "ascii_file_to_db") {
     if (verb_ > 2)
       edm::LogPrint("EcalLaserCondTools") << "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";
   }
 }

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

private

Definition at line 623 of file EcalLaserCondTools.cc.

References EcalCondObjectContainer< T >::barrel(), alignCSCRings::corr, EcalCondObjectContainer< T >::endcap(), Exception, fout_, edm::EventSetup::getData(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), mps_fire::i, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, laserAPDPNRatiosToken_, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, nLmes, AlCaHLTBitMon_ParallelJobs::p, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p3, submitPVValidationJobs::t, timeToString(), EBDetId::validDetId(), EEDetId::validDetId(), and verb_.

Referenced by analyze().

                                                           {
   const auto& laserAPDPNRatios = es.getData(laserAPDPNRatiosToken_);
 
   const EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap& p = laserAPDPNRatios.getLaserMap();
   const EcalLaserAPDPNRatios::EcalLaserTimeStampMap& t = laserAPDPNRatios.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_)
     edm::LogPrint("EcalLaserCondTools") << "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);
         }
       }
     }
   }
 }

void EcalLaserCondTools::fillDb ( CorrReader & r )

private

Definition at line 259 of file EcalLaserCondTools.cc.

References c, eventList_, eventListFileName_, Exception, ntuplemaker::iline, dqmiodumpmetadata::n, nIovs_, nLmes, processIov(), EcalLaserCondTools::CorrReader::readTime(), skipIov_, submitPVValidationJobs::t, timeToString(), and verb_.

Referenced by analyze().

                                              {
   int iIov = 0;
   int processedIovs = 0;
   if (verb_ > 2)
     edm::LogPrint("EcalLaserCondTools") << "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)
         edm::LogPrint("EcalLaserCondTools") << "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)
           edm::LogPrint("EcalLaserCondTools") << "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 {
         edm::LogPrint("EcalLaserCondTools") << "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_) {
         edm::LogPrint("EcalLaserCondTools") << "Skipping IOV " << iIov << "\n";
         continue;
       } else if (processedIovs >= nIovs_ && nIovs_ >= 0) {
         edm::LogPrint("EcalLaserCondTools") << "Requested number of IOVs, " << nIovs_ << ", processed.\n";
         return;
       } else {
         processIov(r, t1, t2, t3);
         ++processedIovs;
       }
     }
 }

void EcalLaserCondTools::from_hdf_to_db ( )

Definition at line 79 of file EcalLaserCondTools.cc.

Referenced by analyze().

                                         {
   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_) {
       edm::LogPrint("EcalLaserCondTools") << " - 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_) {
       edm::LogPrint("EcalLaserCondTools") << "Number of crystals: " << nCrystals << "\n";
       edm::LogPrint("EcalLaserCondTools") << "Number of IOVs: " << nIovs << "\n";
       edm::LogPrint("EcalLaserCondTools") << "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)
       edm::LogPrint("EcalLaserCondTools") << " * 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)
       edm::LogPrint("EcalLaserCondTools") << " * 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)
       edm::LogPrint("EcalLaserCondTools") << " * 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);
 
     EcalLaserAPDPNRatios corrSet;
     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 (!std::isfinite(corr.p1) || !std::isfinite(corr.p2) || !std::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_)
         //      edm::LogPrint("EcalLaserCondTools") << "First IOV, extending starting time.\n";
         //  iovStart = db_->beginOfTime();
         //}
         timeval t;
         gettimeofday(&t, nullptr);
         if (verb_ > 1)
           edm::LogPrint("EcalLaserCondTools")
               << "[" << timeToString(t.tv_sec) << "] "
               << "Write IOV " << iIov << " starting from " << timeToString(iovStart >> 32) << "... ";
         db_->writeOneIOV(corrSet, iovStart, "EcalLaserAPDPNRatiosRcd");
       } catch (const cms::Exception& e) {
         if (verb_ > 1)
           edm::LogPrint("EcalLaserCondTools") << "Failed. ";
         edm::LogPrint("EcalLaserCondTools") << "Exception catched while writting to cond DB" << e.what() << "\n";
       }
       if (verb_ > 1)
         edm::LogPrint("EcalLaserCondTools") << "Suceeded.\n";
 
     }  // loop over IOVs
 
     H5Sclose(memspace);
     H5Fclose(file);
   }  // loop over input files
 }

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

private

Definition at line 334 of file EcalLaserCondTools.cc.

Referenced by fillDb().

                                                                                                   {
   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) {
     edm::LogPrint("EcalLaserCondTools") << "t1:" << t1 << "(" << timeToString(t1) << ") \n"
                                         << "t3: " << t3 << "(" << timeToString(t3) << ")\nt2-t1: ";
     for (int i = 0; i < EcalLaserCondTools::nLmes; ++i) {
       edm::LogPrint("EcalLaserCondTools") << t2[i] - t1 << "\t";
     }
 
     edm::LogPrint("EcalLaserCondTools") << "\n";
   }
   if (t1 < fromTime_) {
     edm::LogPrint("EcalLaserCondTools") << "Skipping IOV " << iIov << ", "
                                         << ", which is before 'fromTime'," << timeToString(fromTime_) << "("
                                         << fromTime_ << ").\n";
     return;
   }
 
   if (toTime_ != -1 && t3 < toTime_) {
     edm::LogPrint("EcalLaserCondTools") << "Skipping IOV " << iIov << ", "
                                         << ", which is beyond 'toTime'," << timeToString(toTime_) << "(" << toTime_
                                         << ").\n";
     return;
   }
 
   if (t1 == 0) {
     edm::LogPrint("EcalLaserCondTools") << "Skipping IOV with t1 = 0"
                                         << "\n";
     return;
   }
 
   EcalLaserAPDPNRatios corrSet;
 
   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)
       edm::LogPrint("EcalLaserCondTools") << "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
       edm::LogPrint("EcalLaserCondTools")
           << "Duplicate det id, for IOV " << iIov << " t1 = " << t1 << " detid = " << int(detid) << "\n";
     }
 
     if (!std::isfinite(corr.p1) || !std::isfinite(corr.p2) || !std::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) {
         edm::LogPrint("EcalLaserCondTools") << EBDetId(detid);
       } else if (detid.subdetId() == EcalEndcap) {
         edm::LogPrint("EcalLaserCondTools") << EEDetId(detid);
       } else {
         edm::LogPrint("EcalLaserCondTools") << (int)detid;
       }
       edm::LogPrint("EcalLaserCondTools") << ": "
                                           << "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_)
         edm::LogPrint("EcalLaserCondTools") << "First IOV, extending starting time.\n";
 
       iovStart = db_->beginOfTime();
     }
     timeval t;
     gettimeofday(&t, nullptr);
     if (verb_ > 1)
       edm::LogPrint("EcalLaserCondTools")
           << "[" << timeToString(t.tv_sec) << "] "
           << "Write IOV " << iIov << " starting from " << timeToString(iovStart >> 32) << "... ";
     db_->writeOneIOV(corrSet, iovStart, "EcalLaserAPDPNRatiosRcd");
   } catch (const cms::Exception& e) {
     edm::LogPrint("EcalLaserCondTools") << "Failed.\nException cathed while writting to cond DB" << e.what() << "\n";
   }
   edm::LogPrint("EcalLaserCondTools") << "Suceeded.\n";
 }

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

staticprivate

Definition at line 614 of file EcalLaserCondTools.cc.

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

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

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

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

staticprivate

Definition at line 591 of file EcalLaserCondTools.cc.

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

Referenced by processIov().

                                          {
   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();
 }