#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 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	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 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 () noexcept

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

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

override

Definition at line 152 of file StoreESCondition.cc.

152 {}

Member Function Documentation

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

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

void StoreESCondition::fillHeader ( char * header )

private

Definition at line 190 of file StoreESCondition.cc.

References EnviromentSettings::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);
 }

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, ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, ESDetId::IX_MAX, ESDetId::IX_MIN, 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;
 }

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

Definition at line 211 of file StoreESCondition.cc.

References contentValuesFiles::fullPath.

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

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

Definition at line 327 of file StoreESCondition.cc.

References contentValuesFiles::fullPath.

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

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

Definition at line 353 of file StoreESCondition.cc.

References esgain_, contentValuesFiles::fullPath, mps_fire::i, 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;
 }

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

Definition at line 251 of file StoreESCondition.cc.

References contentValuesFiles::fullPath.

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

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

Definition at line 268 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, mps_fire::i, ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, ESDetId::IX_MAX, ESDetId::IX_MIN, ESDetId::IY_MAX, ESDetId::IY_MIN, ESDetId::IZ_NUM, ESPedestal::mean, 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;
 }

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

Definition at line 314 of file StoreESCondition.cc.

References contentValuesFiles::fullPath.

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

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

Definition at line 201 of file StoreESCondition.cc.

References contentValuesFiles::fullPath.

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

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

Definition at line 338 of file StoreESCondition.cc.

References contentValuesFiles::fullPath, isotrackApplyRegressor::k, and w.

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

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

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

Definition at line 156 of file StoreESCondition.cc.

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

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

void StoreESCondition::writeToLogFileResults ( char * arg )

Definition at line 177 of file StoreESCondition.cc.

References fillHeader(), 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
 }