#include <StoreESCondition.h>

Inheritance diagram for StoreESCondition:

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

void	endJob () override

std::shared_ptr< ESChannelStatus >	readESChannelStatusFromFile (const char *)

std::shared_ptr< ESEEIntercalibConstants >	readESEEIntercalibConstantsFromFile (const char *)

std::shared_ptr< ESGain >	readESGainFromFile (const char *)

std::shared_ptr< ESIntercalibConstants >	readESIntercalibConstantsFromFile (const char *)

std::shared_ptr< ESMissingEnergyCalibration >	readESMissingEnergyFromFile (const char *)

std::shared_ptr< ESPedestals >	readESPedestalsFromFile (const char *)

std::shared_ptr< ESRecHitRatioCuts >	readESRecHitRatioCutsFromFile (const char *)

std::shared_ptr< ESThresholds >	readESThresholdsFromFile (const char *)

std::shared_ptr< ESTimeSampleWeights >	readESTimeSampleWeightsFromFile (const char *)

	StoreESCondition (const edm::ParameterSet &iConfig)

void	writeToLogFile (std::string, std::string, unsigned long long)

void	writeToLogFileResults (char *)

	~StoreESCondition () 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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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 Member Functions
void	fillHeader (char *)

std::string	to_string (char value[])

Private Attributes
unsigned int	esgain_

std::vector< std::string >	inpFileName_

std::vector< std::string >	inpFileNameEE_

std::string	logfile_

std::vector< std::string >	objectName_

std::string	prog_name_

std::vector< unsigned long long >	since_

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 31 of file StoreESCondition.h.

Constructor & Destructor Documentation

◆ StoreESCondition()

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

explicit

Definition at line 17 of file StoreESCondition.cc.

References esgain_, edm::ParameterSet::getParameter(), inpFileName_, logfile_, objectName_, prog_name_, since_, AlCaHLTBitMon_QueryRunRegistry::string, and convertXMLtoSQLite_cfg::toPut.

                                                                  {
   prog_name_ = "StoreESCondition";
 
   logfile_ = iConfig.getParameter<std::string>("logfile");
 
   esgain_ = iConfig.getParameter<unsigned int>("gain");
 
   typedef std::vector<edm::ParameterSet> Parameters;
   Parameters toPut = iConfig.getParameter<Parameters>("toPut");
   for (Parameters::iterator itToPut = toPut.begin(); itToPut != toPut.end(); ++itToPut) {
     inpFileName_.push_back(itToPut->getUntrackedParameter<std::string>("inputFile"));
     objectName_.push_back(itToPut->getUntrackedParameter<std::string>("conditionType"));
     since_.push_back(itToPut->getUntrackedParameter<unsigned int>("since"));
   }
 }

◆ ~StoreESCondition()

StoreESCondition::~StoreESCondition ( )

override

Definition at line 152 of file StoreESCondition.cc.

152 {}

Member Function Documentation

◆ analyze()

void StoreESCondition::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	evtSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 154 of file StoreESCondition.cc.

154 {}

◆ endJob()

void StoreESCondition::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 33 of file StoreESCondition.cc.

References cond::service::PoolDBOutputService::appendOneIOV(), cond::service::PoolDBOutputService::beginOfTime(), cond::service::PoolDBOutputService::createOneIOV(), mps_fire::i, inpFileName_, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), objectName_, readESChannelStatusFromFile(), readESEEIntercalibConstantsFromFile(), readESGainFromFile(), readESIntercalibConstantsFromFile(), readESMissingEnergyFromFile(), readESPedestalsFromFile(), readESRecHitRatioCutsFromFile(), readESThresholdsFromFile(), readESTimeSampleWeightsFromFile(), since_, AlCaHLTBitMon_QueryRunRegistry::string, and writeToLogFileResults().

                               {
   edm::Service<cond::service::PoolDBOutputService> mydbservice;
   if (!mydbservice.isAvailable()) {
     edm::LogError("StoreESCondition") << "PoolDBOutputService is unavailable"
                                       << "\n";
     return;
   }
 
   bool toAppend = false;
   // copy a string to the char *
   std::string message = "finished OK\n";
   size_t messageSize = message.size() + 1;
   char* messChar = new char[messageSize];
   strncpy(messChar, message.c_str(), messageSize);
 
   for (unsigned int i = 0; i < objectName_.size(); i++) {
     cond::Time_t newTime;
 
     if (mydbservice->isNewTagRequest(objectName_[i] + std::string("Rcd"))) {
       // This is the first object for this tag.
       // Append mode should be off.
       // newTime is the end of this new objects IOV.
       newTime = mydbservice->beginOfTime();
     } else {
       // There should already be an object in the DB for this tag.
       // Append IOV mode should be on.
       // newTime is the beginning of this new objects IOV.
       toAppend = true;
       newTime = (cond::Time_t)since_[i];
     }
     edm::LogInfo("StoreESCondition") << "Reading " << objectName_[i] << " from file and writing to DB with newTime "
                                      << newTime << "\n";
     if (objectName_[i] == "ESChannelStatus") {
       edm::LogInfo("StoreESCondition") << " ESChannelStatus file " << inpFileName_[i] << "\n";
       const auto mycali = readESChannelStatusFromFile(inpFileName_[i].c_str());
       edm::LogInfo("StoreESCondition") << " ESChannelStatus file read "
                                        << "\n";
       if (!toAppend) {
         edm::LogInfo("StoreESCondition") << " before create "
                                          << "\n";
         mydbservice->createOneIOV<ESChannelStatus>(*mycali, newTime, "ESChannelStatusRcd");
         edm::LogInfo("StoreESCondition") << " after create "
                                          << "\n";
       } else {
         edm::LogInfo("StoreESCondition") << " before append "
                                          << "\n";
         mydbservice->appendOneIOV<ESChannelStatus>(*mycali, newTime, "ESChannelStatusRcd");
         edm::LogInfo("StoreESCondition") << " after append "
                                          << "\n";
       }
     } else if (objectName_[i] == "ESIntercalibConstants") {
       const auto myintercalib = readESIntercalibConstantsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESIntercalibConstants>(*myintercalib, newTime, "ESIntercalibConstantsRcd");
       } else {
         mydbservice->appendOneIOV<ESIntercalibConstants>(*myintercalib, newTime, "ESIntercalibConstantsRcd");
       }
     } else if (objectName_[i] == "ESTimeSampleWeights") {
       const auto myintercalib = readESTimeSampleWeightsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESTimeSampleWeights>(*myintercalib, newTime, "ESTimeSampleWeightsRcd");
       } else {
         mydbservice->appendOneIOV<ESTimeSampleWeights>(*myintercalib, newTime, "ESTimeSampleWeightsRcd");
       }
     } else if (objectName_[i] == "ESGain") {
       const auto myintercalib = readESGainFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESGain>(*myintercalib, newTime, "ESGainRcd");
       } else {
         mydbservice->appendOneIOV<ESGain>(*myintercalib, newTime, "ESGainRcd");
       }
     } else if (objectName_[i] == "ESMissingEnergyCalibration") {
       const auto myintercalib = readESMissingEnergyFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESMissingEnergyCalibration>(*myintercalib, newTime, "ESMissingEnergyCalibrationRcd");
       } else {
         mydbservice->appendOneIOV<ESMissingEnergyCalibration>(*myintercalib, newTime, "ESMissingEnergyCalibrationRcd");
       }
     } else if (objectName_[i] == "ESRecHitRatioCuts") {
       const auto myintercalib = readESRecHitRatioCutsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESRecHitRatioCuts>(*myintercalib, newTime, "ESRecHitRatioCutsRcd");
       } else {
         mydbservice->appendOneIOV<ESRecHitRatioCuts>(*myintercalib, newTime, "ESRecHitRatioCutsRcd");
       }
     } else if (objectName_[i] == "ESThresholds") {
       const auto myintercalib = readESThresholdsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESThresholds>(*myintercalib, newTime, "ESThresholdsRcd");
       } else {
         mydbservice->appendOneIOV<ESThresholds>(*myintercalib, newTime, "ESThresholdsRcd");
       }
     } else if (objectName_[i] == "ESPedestals") {
       const auto myintercalib = readESPedestalsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESPedestals>(*myintercalib, newTime, "ESPedestalsRcd");
       } else {
         mydbservice->appendOneIOV<ESPedestals>(*myintercalib, newTime, "ESPedestalsRcd");
       }
     } else if (objectName_[i] == "ESEEIntercalibConstants") {
       const auto myintercalib = readESEEIntercalibConstantsFromFile(inpFileName_[i].c_str());
       if (!toAppend) {
         mydbservice->createOneIOV<ESEEIntercalibConstants>(*myintercalib, newTime, "ESEEIntercalibConstantsRcd");
       } else {
         mydbservice->appendOneIOV<ESEEIntercalibConstants>(*myintercalib, newTime, "ESEEIntercalibConstantsRcd");
       }
     } else {
       edm::LogError("StoreESCondition") << "Object " << objectName_[i] << " is not supported by this program.";
     }
     // if more records write here else if ....
 
     writeToLogFileResults(messChar);
 
     edm::LogInfo("StoreESCondition") << "Finished endJob";
   }
 
   delete[] messChar;
 }

◆ fillHeader()

void StoreESCondition::fillHeader ( char * header )

private

Definition at line 190 of file StoreESCondition.cc.

References RecoTauValidation_cfi::header, hcalRecHitTable_cff::time, and SiStripO2O_cfg_template::user.

Referenced by writeToLogFile(), and writeToLogFileResults().

                                               {
   time_t rawtime;
   struct tm* timeinfo;
   time(&rawtime);
   timeinfo = localtime(&rawtime);
   char user[50];
   strcpy(user, getlogin());
   strcpy(header, asctime(timeinfo));
   strcpy(header, user);
 }

◆ readESChannelStatusFromFile()

std::shared_ptr< ESChannelStatus > StoreESCondition::readESChannelStatusFromFile ( const char * inputFile )

Definition at line 381 of file StoreESCondition.cc.

References beamvalidation::exit(), contentValuesFiles::fullPath, mps_fire::i, makeListRunsInFiles::inputFile, ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ESDetId::IX_MAX, ESDetId::IX_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), ESDetId::IY_MAX, ESDetId::IY_MIN, AlCaHLTBitMon_ParallelJobs::p, mps_update::status, ESDetId::validDetId(), x, y, and z.

Referenced by endJob().

                                                                                                   {
   int z[1000], p[1000], x[1000], y[1000], nsensors;
   std::ifstream statusFile(edm::FileInPath(inputFile).fullPath().c_str());
   statusFile >> nsensors;
   edm::LogInfo("StoreESCondition") << " nsensors " << nsensors << "\n";
   if (nsensors >= 1000) {
     edm::LogInfo("StoreESCondition") << " *** value too high, modify the method!***"
                                      << "\n";
     exit(-1);
   }
   for (int i = 0; i < nsensors; ++i) {
     statusFile >> z[i] >> p[i] >> x[i] >> y[i];
   }
   auto ecalStatus = std::make_shared<ESChannelStatus>();
   int Nbstatus = 0, Nbstrip = 0;
   for (int istrip = ESDetId::ISTRIP_MIN; istrip <= ESDetId::ISTRIP_MAX; istrip++) {
     for (int ix = ESDetId::IX_MIN; ix <= ESDetId::IX_MAX; ix++) {
       for (int iy = ESDetId::IY_MIN; iy <= ESDetId::IY_MAX; iy++) {
         for (int iplane = 1; iplane <= 2; iplane++) {
           for (int izeta = -1; izeta <= 1; izeta = izeta + 2) {
             //      if (izeta==0) continue;
             //      try {
 
             //ESDetId Plane iplane Zside izeta
             //      if(!ESDetId::validDetId(istrip,ix,iy,iplane,izeta)) cout << " Unvalid DetId" << endl;
             //      else {
             if (ESDetId::validDetId(istrip, ix, iy, iplane, izeta)) {
               ESDetId anESId(istrip, ix, iy, iplane, izeta);
               int status = 0;
               //          std::ifstream statusFile(edm::FileInPath(inputFile).fullPath().c_str());
               Nbstrip++;
               for (int i = 0; i < nsensors; ++i) {
                 if (izeta == z[i] && iplane == p[i] && ix == x[i] && iy == y[i])
                   status = 1;
               }
               if (status == 1) {
                 Nbstatus++;
                 if (istrip == 1)
                   edm::LogInfo("StoreESCondition") << " Bad channel ix " << ix << " iy " << iy << " iplane " << iplane
                                                    << " iz " << izeta << "\n";  // print only once
               }
               ecalStatus->setValue(anESId, status);
               //          statusFile.close();
             }  // valid DetId
             //      catch ( cms::Exception &e ) { }
           }  // loop over z
         }  //  loop over plane
       }  //   loop over y
     }  //    loop over x
   }  //     loop over strips
   edm::LogInfo("StoreESCondition") << " Nb of strips " << Nbstrip << " Number of bad channels " << Nbstatus << "\n";
   statusFile.close();
 
   // overwrite the statuses which are in the file
   return ecalStatus;
 }

◆ readESEEIntercalibConstantsFromFile()

std::shared_ptr< ESEEIntercalibConstants > StoreESCondition::readESEEIntercalibConstantsFromFile ( const char * inputFile )

Definition at line 211 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, and makeListRunsInFiles::inputFile.

Referenced by endJob().

                                                                                                                   {
   std::ifstream ESEEIntercalibFile(edm::FileInPath(inputFile).fullPath().c_str());
   float gammaLow0, alphaLow0, gammaHigh0, alphaHigh0, gammaLow1, alphaLow1, gammaHigh1, alphaHigh1, gammaLow2,
       alphaLow2, gammaHigh2, alphaHigh2, gammaLow3, alphaLow3, gammaHigh3, alphaHigh3;
   ESEEIntercalibFile >> gammaLow0;
   ESEEIntercalibFile >> alphaLow0;
   ESEEIntercalibFile >> gammaHigh0;
   ESEEIntercalibFile >> alphaHigh0;
   ESEEIntercalibFile >> gammaLow1;
   ESEEIntercalibFile >> alphaLow1;
   ESEEIntercalibFile >> gammaHigh1;
   ESEEIntercalibFile >> alphaHigh1;
   ESEEIntercalibFile >> gammaLow2;
   ESEEIntercalibFile >> alphaLow2;
   ESEEIntercalibFile >> gammaHigh2;
   ESEEIntercalibFile >> alphaHigh2;
   ESEEIntercalibFile >> gammaLow3;
   ESEEIntercalibFile >> alphaLow3;
   ESEEIntercalibFile >> gammaHigh3;
   ESEEIntercalibFile >> alphaHigh3;
   auto eseeIntercalibConstants = std::make_shared<ESEEIntercalibConstants>(gammaLow0,
                                                                            alphaLow0,
                                                                            gammaHigh0,
                                                                            alphaHigh0,
                                                                            gammaLow1,
                                                                            alphaLow1,
                                                                            gammaHigh1,
                                                                            alphaHigh1,
                                                                            gammaLow2,
                                                                            alphaLow2,
                                                                            gammaHigh2,
                                                                            alphaHigh2,
                                                                            gammaLow3,
                                                                            alphaLow3,
                                                                            gammaHigh3,
                                                                            alphaHigh3);
 
   return eseeIntercalibConstants;
 }

◆ readESGainFromFile()

std::shared_ptr< ESGain > StoreESCondition::readESGainFromFile ( const char * inputFile )

Definition at line 327 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, Hcal_Conditions_forGlobalTag_cff::gain, and makeListRunsInFiles::inputFile.

Referenced by endJob().

                                                                                 {
   std::ifstream amplFile(edm::FileInPath(inputFile).fullPath().c_str());
 
   int gain;
   amplFile >> gain;
   edm::LogInfo("StoreESCondition") << "gain : " << gain << "\n";
 
   auto esGain = std::make_shared<ESGain>(gain);  // 1: LG, 2: HG
   return esGain;
 }

◆ readESIntercalibConstantsFromFile()

std::shared_ptr< ESIntercalibConstants > StoreESCondition::readESIntercalibConstantsFromFile ( const char * inputFile )

Definition at line 353 of file StoreESCondition.cc.

References esgain_, contentValuesFiles::fullPath, mps_fire::i, makeListRunsInFiles::inputFile, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), DetId::rawId(), and ESDetId::validDetId().

Referenced by endJob().

                                                                                                               {
   auto ical = std::make_shared<ESIntercalibConstants>();
 
   std::ifstream mipFile(edm::FileInPath(inputFile).fullPath().c_str());
 
   for (int i = 0; i < 137216; ++i) {
     int iz, ip, ix, iy, is;
     double mip;
     mipFile >> iz >> ip >> ix >> iy >> is >> mip;
     //if (mip <20 || mip> 70) cout<<iz<<" "<<ip<<" "<<ix<<" "<<iy<<" "<<is<<" "<<mip<<endl; // HG
     // LG : HG MIP/6/1.14
     //mip = mip/6/1.14;
     // LG : HG MIP/6
     if (esgain_ == 1)
       mip = mip / 6.;  // LG
     if (mip < 20 || mip > 70)
       edm::LogInfo("StoreESCondition") << iz << " " << ip << " " << ix << " " << iy << " " << is << " " << mip
                                        << "\n";  // LG
 
     if (ESDetId::validDetId(is, ix, iy, ip, iz)) {
       ESDetId esId(is, ix, iy, ip, iz);
       ical->setValue(esId.rawId(), mip);
     }
   }
 
   return ical;
 }

◆ readESMissingEnergyFromFile()

std::shared_ptr< ESMissingEnergyCalibration > StoreESCondition::readESMissingEnergyFromFile ( const char * inputFile )

Definition at line 251 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, and makeListRunsInFiles::inputFile.

Referenced by endJob().

                                                                                                              {
   std::ifstream ESMissingEnergyFile(edm::FileInPath(inputFile).fullPath().c_str());
   float ConstAEta0, ConstBEta0, ConstAEta1, ConstBEta1, ConstAEta2, ConstBEta2, ConstAEta3, ConstBEta3;
   ESMissingEnergyFile >> ConstAEta0;
   ESMissingEnergyFile >> ConstBEta0;
   ESMissingEnergyFile >> ConstAEta1;
   ESMissingEnergyFile >> ConstBEta1;
   ESMissingEnergyFile >> ConstAEta2;
   ESMissingEnergyFile >> ConstBEta2;
   ESMissingEnergyFile >> ConstAEta3;
   ESMissingEnergyFile >> ConstBEta3;
   auto esMissingEnergy = std::make_shared<ESMissingEnergyCalibration>(
       ConstAEta0, ConstBEta0, ConstAEta1, ConstBEta1, ConstAEta2, ConstBEta2, ConstAEta3, ConstBEta3);
 
   return esMissingEnergy;
 }

◆ readESPedestalsFromFile()

std::shared_ptr< ESPedestals > StoreESCondition::readESPedestalsFromFile ( const char * inputFile )

Definition at line 268 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, mps_fire::i, makeListRunsInFiles::inputFile, ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ESDetId::IX_MAX, ESDetId::IX_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), ESDetId::IY_MAX, ESDetId::IY_MIN, ESDetId::IZ_NUM, ESPedestal::mean, QIE10Task_cfi::ped, ESDetId::PLANE_MAX, ESDetId::PLANE_MIN, DetId::rawId(), ESPedestal::rms, ESDetId::validDetId(), and ecaldqm::zside().

Referenced by endJob().

                                                                                           {
   auto esPedestals = std::make_shared<ESPedestals>();
 
   // int ped[2][2][40][40][32];
   // for (int i=0; i<2; ++i)
   //  for (int j=0; j<2; ++j)
   //    for (int k=0; k<40; ++k)
   //    for (int m=0; m<40; ++m)
   //      for (int n=0; n<32; ++n)
   //        ped[i][j][k][m][n] = 0;
 
   int ped[ESDetId::IZ_NUM][ESDetId::PLANE_MAX][ESDetId::IX_MAX][ESDetId::IY_MAX][ESDetId::ISTRIP_MAX] = {};
 
   int iz, ip, ix, iy, is, ped_, zside;
   std::ifstream pedestalFile(edm::FileInPath(inputFile).fullPath().c_str());
 
   for (int i = 0; i < 137216; ++i) {
     pedestalFile >> iz >> ip >> ix >> iy >> is >> ped_;
 
     zside = (iz == -1) ? 1 : 0;
     ped[zside][ip - 1][ix - 1][iy - 1][is - 1] = ped_;
   }
 
   for (int iz = -1; iz <= 1; ++iz) {
     if (iz == 0)
       continue;
     zside = (iz == -1) ? 1 : 0;
 
     for (int iplane = ESDetId::PLANE_MIN; iplane <= ESDetId::PLANE_MAX; ++iplane)
       for (int ix = ESDetId::IX_MIN; ix <= ESDetId::IX_MAX; ++ix)
         for (int iy = ESDetId::IY_MIN; iy <= ESDetId::IY_MAX; ++iy)
           for (int istrip = ESDetId::ISTRIP_MIN; istrip <= ESDetId::ISTRIP_MAX; ++istrip) {
             ESPedestals::Item ESitem;
             ESitem.mean = ped[zside][iplane - 1][ix - 1][iy - 1][istrip - 1];
             ESitem.rms = 3;  // LG : 3, HG : 6
 
             if (ESDetId::validDetId(istrip, ix, iy, iplane, iz)) {
               ESDetId esId(istrip, ix, iy, iplane, iz);
               esPedestals->insert(std::make_pair(esId.rawId(), ESitem));
             }
           }
   }
 
   return esPedestals;
 }

◆ readESRecHitRatioCutsFromFile()

std::shared_ptr< ESRecHitRatioCuts > StoreESCondition::readESRecHitRatioCutsFromFile ( const char * inputFile )

Definition at line 314 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, and makeListRunsInFiles::inputFile.

Referenced by endJob().

                                                                                                       {
   std::ifstream ESRecHitRatioCutsFile(edm::FileInPath(inputFile).fullPath().c_str());
 
   float r12Low, r23Low, r12High, r23High;
   ESRecHitRatioCutsFile >> r12Low;
   ESRecHitRatioCutsFile >> r23Low;
   ESRecHitRatioCutsFile >> r12High;
   ESRecHitRatioCutsFile >> r23High;
   auto esRecHitRatioCuts = std::make_shared<ESRecHitRatioCuts>(r12Low, r23Low, r12High, r23High);
 
   return esRecHitRatioCuts;
 }

◆ readESThresholdsFromFile()

std::shared_ptr< ESThresholds > StoreESCondition::readESThresholdsFromFile ( const char * inputFile )

Definition at line 201 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, and makeListRunsInFiles::inputFile.

Referenced by endJob().

                                                                                             {
   std::ifstream ESThresholdsFile(edm::FileInPath(inputFile).fullPath().c_str());
   float ts2, zs;  //2nd time sample, ZS threshold
   ESThresholdsFile >> ts2;
   ESThresholdsFile >> zs;
   auto esThresholds = std::make_shared<ESThresholds>(ts2, zs);
 
   return esThresholds;
 }

◆ readESTimeSampleWeightsFromFile()

std::shared_ptr< ESTimeSampleWeights > StoreESCondition::readESTimeSampleWeightsFromFile ( const char * inputFile )

Definition at line 338 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, makeListRunsInFiles::inputFile, isotrackApplyRegressor::k, w(), and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::ww.

Referenced by endJob().

                                                                                                           {
   std::ifstream amplFile(edm::FileInPath(inputFile).fullPath().c_str());
 
   float w[3];
   for (int k = 0; k < 3; ++k) {
     float ww;
     amplFile >> ww;
     w[k] = ww;
     edm::LogInfo("StoreESCondition") << "weight : " << k << " " << w[k] << "\n";
   }
 
   auto esWeights = std::make_shared<ESTimeSampleWeights>(w[0], w[1], w[2]);
   return esWeights;
 }

◆ to_string()

std::string StoreESCondition::to_string ( char value[] )

inlineprivate

Definition at line 63 of file StoreESCondition.h.

References relativeConstraints::value.

                                     {
     std::ostringstream streamOut;
     streamOut << value;
     return streamOut.str();
   }

◆ writeToLogFile()

void StoreESCondition::writeToLogFile	(	std::string	,
		std::string	,
		unsigned long	long
	)

Definition at line 156 of file StoreESCondition.cc.

References fillHeader(), RecoTauValidation_cfi::header, logfile_, produceTPGParameters_beamv6_transparency_spikekill_2016_script::outFile, and SiStripO2O_cfg_template::since.

                                                                                   {
   FILE* outFile;  // output log file for appending
   outFile = fopen(logfile_.c_str(), "a");
   if (!outFile) {
     edm::LogError("StoreESCondition") << "*** Can not open file: " << logfile_;
     return;
   }
   char header[256];
   fillHeader(header);
   char appendMode[10];
   if (since != 0)
     strcpy(appendMode, "append");
   else
     strcpy(appendMode, "create");
 
   //fprintf(outFile, "%s %s condition from file %s written into DB for SM %d (mapped to SM %d) in %s mode (since run %u)\n",
   //header, a.c_str(), b.c_str(),  appendMode, (unsigned int)since);
 
   fclose(outFile);  // close out file
 }

◆ writeToLogFileResults()

void StoreESCondition::writeToLogFileResults ( char * arg )

Definition at line 177 of file StoreESCondition.cc.

References fillHeader(), RecoTauValidation_cfi::header, logfile_, and produceTPGParameters_beamv6_transparency_spikekill_2016_script::outFile.

Referenced by endJob().

                                                       {
   FILE* outFile;  // output log file for appending
   outFile = fopen(logfile_.c_str(), "a");
   if (!outFile) {
     edm::LogError("StoreESCondition") << "*** Can not open file: " << logfile_;
     return;
   }
   char header[256];
   fillHeader(header);
   fprintf(outFile, "%s %s\n", header, arg);
   fclose(outFile);  // close out file
 }

Member Data Documentation

◆ esgain_

unsigned int StoreESCondition::esgain_

private

Definition at line 61 of file StoreESCondition.h.

Referenced by readESIntercalibConstantsFromFile(), and StoreESCondition().

◆ inpFileName_

std::vector<std::string> StoreESCondition::inpFileName_

private

Definition at line 55 of file StoreESCondition.h.

Referenced by endJob(), and StoreESCondition().

◆ inpFileNameEE_

std::vector<std::string> StoreESCondition::inpFileNameEE_

private

Definition at line 56 of file StoreESCondition.h.

◆ logfile_

std::string StoreESCondition::logfile_

private

Definition at line 59 of file StoreESCondition.h.

Referenced by StoreESCondition(), writeToLogFile(), and writeToLogFileResults().

◆ objectName_

std::vector<std::string> StoreESCondition::objectName_

private

Definition at line 54 of file StoreESCondition.h.

Referenced by endJob(), and StoreESCondition().

◆ prog_name_

std::string StoreESCondition::prog_name_

private

Definition at line 57 of file StoreESCondition.h.

Referenced by StoreESCondition().

◆ since_

std::vector<unsigned long long> StoreESCondition::since_

private

Definition at line 58 of file StoreESCondition.h.

Referenced by endJob(), and StoreESCondition().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ StoreESCondition()

◆ ~StoreESCondition()

Member Function Documentation

◆ analyze()

◆ endJob()

◆ fillHeader()

◆ readESChannelStatusFromFile()

◆ readESEEIntercalibConstantsFromFile()

◆ readESGainFromFile()

◆ readESIntercalibConstantsFromFile()

◆ readESMissingEnergyFromFile()

◆ readESPedestalsFromFile()

◆ readESRecHitRatioCutsFromFile()

◆ readESThresholdsFromFile()

◆ readESTimeSampleWeightsFromFile()

◆ to_string()

◆ writeToLogFile()

◆ writeToLogFileResults()

Member Data Documentation

◆ esgain_

◆ inpFileName_

◆ inpFileNameEE_

◆ logfile_

◆ objectName_

◆ prog_name_

◆ since_