calculate the best DAC value to obtain a pedestal = 200 More...

#include <EcalPedOffset.h>

Inheritance diagram for EcalPedOffset:

Public Member Functions
void	analyze (edm::Event const &event, edm::EventSetup const &eventSetup) override
	! Analyze More...

void	beginRun (edm::Run const &, edm::EventSetup const &eventSetup) override
	BeginRun. More...

	EcalPedOffset (const edm::ParameterSet &ps)
	Constructor. More...

void	endJob (void) override
	EndJob. More...

void	makePlots ()
	create the plots of the DAC pedestal trend More...

void	writeDb ()
	WriteDB. More...

void	writeXMLFiles (std::string fileName)
	write the results into xml format More...

	~EcalPedOffset () override
	Destructor. More...

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
std::string	intToString (int num)

void	readDACs (const edm::Handle< EBDigiCollection > &pDigis, const std::map< int, int > &DACvalues)

void	readDACs (const edm::Handle< EEDigiCollection > &pDigis, const std::map< int, int > &DACvalues)

Private Attributes
const EcalElectronicsMapping *	ecalElectronicsMap_

edm::InputTag	m_barrelDigiCollection
	secondary name given to collection of digis More...

int	m_bestPed

bool	m_create_moniov

int	m_DACmax

int	m_DACmin

std::string	m_dbHostName
	database host name More...

int	m_dbHostPort
	database More...

std::string	m_dbName
	database name More...

std::string	m_dbPassword
	database user password More...

std::string	m_dbUserName
	database user name More...

edm::InputTag	m_endcapDigiCollection
	secondary name given to collection of digis More...

edm::InputTag	m_headerCollection
	name of module/plugin/producer making headers More...

std::string	m_location

double	m_maxChi2OverNDFAllowed_
	max chi2/ndf allowed for linearity test More...

double	m_maxSlopeAllowed_
	max slope (in magnitude) allowed for linearity test More...

double	m_minSlopeAllowed_
	min slope (in magnitude) allowed for linearity test More...

std::map< int, TPedResult * >	m_pedResult

std::map< int, TPedValues * >	m_pedValues

std::string	m_plotting
	the root file where to store the detail plots More...

double	m_RMSmax

int	m_run
	run number More...

std::string	m_xmlFile
	name of the xml file to be saved More...

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

calculate the best DAC value to obtain a pedestal = 200

Author: P. Govoni (pietr.nosp@m.o.go.nosp@m.voni@.nosp@m.cern.nosp@m.NOSPA.nosp@m.M.ch)

Definition at line 28 of file EcalPedOffset.h.

Constructor & Destructor Documentation

◆ EcalPedOffset()

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

Constructor.

ctor

Definition at line 38 of file EcalPedOffset.cc.

     : m_barrelDigiCollection(paramSet.getParameter<edm::InputTag>("EBdigiCollection")),
       m_endcapDigiCollection(paramSet.getParameter<edm::InputTag>("EEdigiCollection")),
       m_headerCollection(paramSet.getParameter<edm::InputTag>("headerCollection")),
       m_xmlFile(paramSet.getParameter<std::string>("xmlFile")),
       m_DACmin(paramSet.getUntrackedParameter<int>("DACmin", 0)),
       m_DACmax(paramSet.getUntrackedParameter<int>("DACmax", 256)),
       m_RMSmax(paramSet.getUntrackedParameter<double>("RMSmax", 2)),
       m_bestPed(paramSet.getUntrackedParameter<int>("bestPed", 200)),
       m_dbHostName(paramSet.getUntrackedParameter<std::string>("dbHostName", "0")),
       m_dbName(paramSet.getUntrackedParameter<std::string>("dbName", "0")),
       m_dbUserName(paramSet.getUntrackedParameter<std::string>("dbUserName")),
       m_dbPassword(paramSet.getUntrackedParameter<std::string>("dbPassword")),
       m_dbHostPort(paramSet.getUntrackedParameter<int>("dbHostPort", 1521)),
       m_create_moniov(paramSet.getUntrackedParameter<bool>("createMonIOV", false)),
       m_location(paramSet.getUntrackedParameter<std::string>("location", "H4")),
       m_run(-1),
       m_plotting(paramSet.getParameter<std::string>("plotting")),
       m_maxSlopeAllowed_(paramSet.getUntrackedParameter<double>("maxSlopeAllowed", -29)),
       m_minSlopeAllowed_(paramSet.getUntrackedParameter<double>("minSlopeAllowed", -18)),
       m_maxChi2OverNDFAllowed_(paramSet.getUntrackedParameter<double>("maxChi2OverNDF", 5)) {
   edm::LogInfo("EcalPedOffset") << " reading "
                                 << " m_DACmin: " << m_DACmin << " m_DACmax: " << m_DACmax << " m_RMSmax: " << m_RMSmax
                                 << " m_bestPed: " << m_bestPed;
 }

References m_bestPed, m_DACmax, m_DACmin, and m_RMSmax.

◆ ~EcalPedOffset()

EcalPedOffset::~EcalPedOffset ( )

override

Destructor.

dtor

Definition at line 67 of file EcalPedOffset.cc.

                               {
   for (std::map<int, TPedValues *>::iterator mapIt = m_pedValues.begin(); mapIt != m_pedValues.end(); ++mapIt)
     delete mapIt->second;
   for (std::map<int, TPedResult *>::iterator mapIt = m_pedResult.begin(); mapIt != m_pedResult.end(); ++mapIt)
     delete mapIt->second;
 }

References cms::cuda::for(), m_pedResult, and m_pedValues.

Member Function Documentation

◆ analyze()

void EcalPedOffset::analyze	(	edm::Event const &	event,
		edm::EventSetup const &	eventSetup
	)

overridevirtual

! Analyze

perform the analysis

Implements edm::EDAnalyzer.

Definition at line 88 of file EcalPedOffset.cc.

                                                                             {
   LogDebug("EcalPedOffset") << "entering analyze ...";
  
   // get the headers
   // (one header for each supermodule)
   edm::Handle<EcalRawDataCollection> DCCHeaders;
   event.getByLabel(m_headerCollection, DCCHeaders);
  
   std::map<int, int> DACvalues;
  
   if (m_run == -1)
     m_run = event.id().run();
  
   // loop over the headers
   for (EcalRawDataCollection::const_iterator headerItr = DCCHeaders->begin(); headerItr != DCCHeaders->end();
        ++headerItr) {
     EcalDCCHeaderBlock::EcalDCCEventSettings settings = headerItr->getEventSettings();
     int FEDid = 600 + headerItr->id();
     DACvalues[FEDid] = settings.ped_offset;
     LogDebug("EcalPedOffset") << "Found FED: " << FEDid << " in DCC header";
   }
  
   bool barrelDigisFound = true;
   bool endcapDigisFound = true;
   // get the barrel digis
   // (one digi for each crystal)
   Handle<EBDigiCollection> barrelDigis;
   event.getByLabel(m_barrelDigiCollection, barrelDigis);
   if (!barrelDigis.isValid()) {
     edm::LogError("EcalPedOffset") << "Error! can't get the product " << m_barrelDigiCollection
                                    << "; not reading barrel digis";
     barrelDigisFound = false;
   }
  
   if (barrelDigis->empty()) {
     edm::LogInfo("EcalPedOffset") << "Size of EBDigiCollection is zero;"
                                   << " not reading barrel digis";
     barrelDigisFound = false;
   }
  
   // get the endcap digis
   // (one digi for each crystal)
   Handle<EEDigiCollection> endcapDigis;
   event.getByLabel(m_endcapDigiCollection, endcapDigis);
   if (!endcapDigis.isValid()) {
     edm::LogError("EcalPedOffset") << "Error! can't get the product " << m_endcapDigiCollection
                                    << "; not reading endcap digis";
     endcapDigisFound = false;
   }
  
   if (endcapDigis->empty()) {
     edm::LogInfo("EcalPedOffset") << "Size of EEDigiCollection is zero;"
                                   << " not reading endcap digis";
     endcapDigisFound = false;
   }
  
   if (barrelDigisFound)
     readDACs(barrelDigis, DACvalues);
   if (endcapDigisFound)
     readDACs(endcapDigis, DACvalues);
   if (!barrelDigisFound && !endcapDigisFound)
     edm::LogError("EcalPedOffset") << "No digis found in the event!";
 }

References mtdUncalibratedRecHits_cfi::barrelDigis, edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), mtdUncalibratedRecHits_cfi::endcapDigis, LogDebug, m_barrelDigiCollection, m_endcapDigiCollection, m_headerCollection, m_run, EcalDCCHeaderBlock::EcalDCCEventSettings::ped_offset, and readDACs().

◆ beginRun()

void EcalPedOffset::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	eventSetup
	)

overridevirtual

BeginRun.

begin the run

Reimplemented from edm::EDAnalyzer.

Definition at line 77 of file EcalPedOffset.cc.

                                                                       {
   LogDebug("EcalPedOffset") << "entering beginRun...";
  
   edm::ESHandle<EcalElectronicsMapping> handle;
   eventSetup.get<EcalMappingRcd>().get(handle);
   ecalElectronicsMap_ = handle.product();
 }

References ecalElectronicsMap_, edm::EventSetup::get(), get, patZpeak::handle, and LogDebug.

◆ endJob()

void EcalPedOffset::endJob ( void )

overridevirtual

EndJob.

perform the minimization and write results

Reimplemented from edm::EDAnalyzer.

Definition at line 223 of file EcalPedOffset.cc.

                            {
   for (std::map<int, TPedValues *>::const_iterator smPeds = m_pedValues.begin(); smPeds != m_pedValues.end();
        ++smPeds) {
     m_pedResult[smPeds->first] = new TPedResult((smPeds->second)->terminate(m_DACmin, m_DACmax));
   }
   edm::LogInfo("EcalPedOffset") << " results map size " << m_pedResult.size();
   writeXMLFiles(m_xmlFile);
  
   if (m_plotting != '0')
     makePlots();
   if (m_dbHostName != '0')
     writeDb();
 }

References m_DACmax, m_DACmin, m_dbHostName, m_pedResult, m_pedValues, m_plotting, m_xmlFile, makePlots(), writeDb(), and writeXMLFiles().

Referenced by o2olib.O2ORunMgr::executeJob().

◆ intToString()

std::string EcalPedOffset::intToString ( int num )

private

Definition at line 457 of file EcalPedOffset.cc.

                                             {
   // outputs the number into the string stream and then flushes
   // the buffer (makes sure the output is put into the stream)
   std::ostringstream myStream;
   myStream << num << std::flush;
   return (myStream.str());  // returns the string form of the stringstream object
 }

References EgammaValidation_cff::num.

Referenced by writeXMLFiles().

◆ makePlots()

void EcalPedOffset::makePlots ( )

create the plots of the DAC pedestal trend

Definition at line 428 of file EcalPedOffset.cc.

                               {
   LogDebug("EcalPedOffset") << " entering makePlots ...";
  
   edm::LogInfo("EcalPedOffset") << " map size: " << m_pedValues.size();
  
   // create the ROOT file
   m_plotting += ".root";
  
   TFile *rootFile = new TFile(m_plotting.c_str(), "RECREATE");
  
   // loop over the supermodules
   for (std::map<int, TPedValues *>::const_iterator smPeds = m_pedValues.begin(); smPeds != m_pedValues.end();
        ++smPeds) {
     // make a folder in the ROOT file
     char folderName[120];
     sprintf(folderName, "FED%02d", smPeds->first);
     rootFile->mkdir(folderName);
     smPeds->second->makePlots(rootFile, folderName, m_maxSlopeAllowed_, m_minSlopeAllowed_, m_maxChi2OverNDFAllowed_);
   }
  
   rootFile->Close();
   delete rootFile;
  
   LogDebug("EcalPedOffset") << " DONE";
 }

References ALCARECODTCalibSynchCosmicsDQM_cff::folderName, LogDebug, m_maxChi2OverNDFAllowed_, m_maxSlopeAllowed_, m_minSlopeAllowed_, m_pedValues, and m_plotting.

Referenced by endJob().

◆ readDACs() [1/2]

void EcalPedOffset::readDACs	(	const edm::Handle< EBDigiCollection > &	pDigis,
		const std::map< int, int > &	DACvalues
	)

private

Definition at line 154 of file EcalPedOffset.cc.

                                                                                                           {
   std::map<int, int> DACvalues = _DACvalues;
   // loop over the digis
   for (EBDigiCollection::const_iterator itdigi = pDigis->begin(); itdigi != pDigis->end(); ++itdigi) {
     int gainId = ((EBDataFrame)(*itdigi)).sample(0).gainId();
     EBDetId detId = EBDetId(itdigi->id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
     int FEDid = 600 + elecId.dccId();
     int crystalId = detId.ic();
  
     // TODO: Behavior here
     if (DACvalues.find(FEDid) == DACvalues.end()) {
       edm::LogError("EcalPedOffset") << "Error! EB DCCid of digi does not "
                                         "match any DCCid found in DCC headers"
                                      << FEDid;
     }
  
     if (!m_pedValues.count(FEDid)) {
       LogDebug("EcalPedOffset") << "Inserting new TPedValues object for FED:" << FEDid;
       m_pedValues[FEDid] = new TPedValues(m_RMSmax, m_bestPed);
     }
  
     // loop over the samples
     for (int iSample = 0; iSample < EcalDataFrame::MAXSAMPLES; ++iSample) {
       m_pedValues[FEDid]->insert(
           gainId, crystalId, DACvalues[FEDid], ((EBDataFrame)(*itdigi)).sample(iSample).adc(), crystalId);
     }
  
   }  // end loop over digis
 }

References gpuClustering::adc, edm::DataFrameContainer::begin(), EcalElectronicsId::dccId(), ecalElectronicsMap_, edm::DataFrameContainer::end(), ecalLiteDTU::gainId(), EcalElectronicsMapping::getElectronicsId(), EBDetId::ic(), LogDebug, m_bestPed, m_pedValues, m_RMSmax, EcalDataFrame::MAXSAMPLES, and simplePhotonAnalyzer_cfi::sample.

Referenced by analyze().

◆ readDACs() [2/2]

void EcalPedOffset::readDACs	(	const edm::Handle< EEDigiCollection > &	pDigis,
		const std::map< int, int > &	DACvalues
	)

private

Definition at line 187 of file EcalPedOffset.cc.

                                                                                                           {
   std::map<int, int> DACvalues = _DACvalues;
   // loop over the digis
   for (EEDigiCollection::const_iterator itdigi = pDigis->begin(); itdigi != pDigis->end(); ++itdigi) {
     int gainId = ((EEDataFrame)(*itdigi)).sample(0).gainId();
     // int gainId = itdigi->sample(0).gainId();
     EEDetId detId = EEDetId(itdigi->id());
     EcalElectronicsId elecId = ecalElectronicsMap_->getElectronicsId(detId);
     int FEDid = 600 + elecId.dccId();
     int crystalId = 25 * elecId.towerId() + 5 * elecId.stripId() + elecId.xtalId();
     int endcapCrystalId = 100 * elecId.towerId() + 5 * (elecId.stripId() - 1) + elecId.xtalId();
  
     // TODO: Behavior here
     if (DACvalues.find(FEDid) == DACvalues.end()) {
       edm::LogError("EcalPedOffset") << "Error! EE DCCid of digi does not "
                                         "match any DCCid found in DCC headers: "
                                      << FEDid;
     }
  
     if (!m_pedValues.count(FEDid)) {
       LogDebug("EcalPedOffset") << "Inserting new TPedValues object for FED:" << FEDid;
       m_pedValues[FEDid] = new TPedValues(m_RMSmax, m_bestPed);
     }
  
     // loop over the samples
     for (int iSample = 0; iSample < EcalDataFrame::MAXSAMPLES; ++iSample) {
       m_pedValues[FEDid]->insert(
           gainId, crystalId, DACvalues[FEDid], ((EEDataFrame)(*itdigi)).sample(iSample).adc(), endcapCrystalId);
     }
  
   }  // end loop over digis
 }

References gpuClustering::adc, edm::DataFrameContainer::begin(), EcalElectronicsId::dccId(), ecalElectronicsMap_, edm::DataFrameContainer::end(), ecalLiteDTU::gainId(), EcalElectronicsMapping::getElectronicsId(), LogDebug, m_bestPed, m_pedValues, m_RMSmax, EcalDataFrame::MAXSAMPLES, and simplePhotonAnalyzer_cfi::sample.

◆ writeDb()

void EcalPedOffset::writeDb ( )

WriteDB.

write the m_pedResult in the DB FIXME divide into sub-tasks

Definition at line 241 of file EcalPedOffset.cc.

                             {
   LogDebug("EcalPedOffset") << " entering writeDb ...";
  
   // connect to the database
   EcalCondDBInterface *DBconnection;
   try {
     LogInfo("EcalPedOffset") << "Opening DB connection with TNS_ADMIN ...";
     DBconnection = new EcalCondDBInterface(m_dbName, m_dbUserName, m_dbPassword);
   } catch (std::runtime_error &e) {
     LogError("EcalPedOffset") << e.what();
     if (!m_dbHostName.empty()) {
       try {
         LogInfo("EcalPedOffset") << "Opening DB connection without TNS_ADMIN ...";
         DBconnection = new EcalCondDBInterface(m_dbHostName, m_dbName, m_dbUserName, m_dbPassword, m_dbHostPort);
       } catch (std::runtime_error &e) {
         LogError("EcalPedOffset") << e.what();
         return;
       }
     } else
       return;
   }
  
   // define the query for RunIOV to get the right place in the database
   RunTag runtag;
   LocationDef locdef;
   RunTypeDef rundef;
   locdef.setLocation(m_location);
  
   runtag.setGeneralTag("PEDESTAL-OFFSET");
   rundef.setRunType("PEDESTAL-OFFSET");
   // rundef.setRunType ("TEST");
   // runtag.setGeneralTag ("TEST");
  
   runtag.setLocationDef(locdef);
   runtag.setRunTypeDef(rundef);
  
   run_t run = m_run;  // FIXME dal config file
   // RunIOV runiov = DBconnection->fetchRunIOV (&runtag, run);
   RunIOV runiov = DBconnection->fetchRunIOV(m_location, run);
  
   // MonRunIOV
   MonVersionDef monverdef;
   monverdef.setMonitoringVersion("test01");
   MonRunTag montag;
   montag.setMonVersionDef(monverdef);
   montag.setGeneralTag("CMSSW");
  
   subrun_t subrun = 1;  // hardcoded!
  
   MonRunIOV moniov;
  
   try {
     runtag = runiov.getRunTag();
     moniov = DBconnection->fetchMonRunIOV(&runtag, &montag, run, subrun);
   } catch (std::runtime_error &e) {
     if (m_create_moniov) {
       // if not already in the DB create a new MonRunIOV
       Tm startSubRun;
       startSubRun.setToCurrentGMTime();
       // setup the MonIOV
       moniov.setRunIOV(runiov);
       moniov.setSubRunNumber(subrun);
       moniov.setSubRunStart(startSubRun);
       moniov.setMonRunTag(montag);
       LogDebug("EcalPedOffset") << " creating a new MonRunIOV";
     } else {
       edm::LogError("EcalPedOffset") << " no MonRunIOV existing in the DB";
       edm::LogError("EcalPedOffset") << " the result will not be stored into the DB";
       if (DBconnection) {
         delete DBconnection;
       }
       return;
     }
   }
  
   // create the table to be filled and the map to be inserted
   EcalLogicID ecid;
   std::map<EcalLogicID, MonPedestalOffsetsDat> DBdataset;
   MonPedestalOffsetsDat DBtable;
  
   // fill the table
  
   // loop over the super-modules
   for (std::map<int, TPedResult *>::const_iterator result = m_pedResult.begin(); result != m_pedResult.end();
        ++result) {
     // loop over the crystals
     for (int xtal = 0; xtal < 1700; ++xtal) {
       DBtable.setDACG1(result->second->m_DACvalue[2][xtal]);
       DBtable.setDACG6(result->second->m_DACvalue[1][xtal]);
       DBtable.setDACG12(result->second->m_DACvalue[0][xtal]);
       DBtable.setTaskStatus(true);  // FIXME to be set correctly
  
       // fill the table
       if (DBconnection) {
         try {
           int fedid = result->first;
           int eid = m_pedValues[fedid]->getCrystalNumber(xtal);
           // If eid is zero, that crystal was not present in digis
           if (eid == 0)
             continue;
  
           if (fedid >= 601 && fedid <= 609) {
             // Add the FEDid part in for DB
             eid = eid + 10000 * (fedid - 600);
             ecid = DBconnection->getEcalLogicID("EE_elec_crystal_number", eid);
           } else if (fedid >= 610 && fedid <= 627) {
             ecid = DBconnection->getEcalLogicID("EB_crystal_number", fedid - 610 + 19, eid);
           } else if (fedid >= 628 && fedid <= 645) {
             ecid = DBconnection->getEcalLogicID("EB_crystal_number", fedid - 628 + 1, eid);
           } else if (fedid >= 646 && fedid <= 654) {
             // Add the FEDid part in for DB
             eid = eid + 10000 * (fedid - 600);
             ecid = DBconnection->getEcalLogicID("EE_elec_crystal_number", eid);
           } else
             LogError("EcalPedOffset") << "FEDid is out of range 601-654";
  
           DBdataset[ecid] = DBtable;
         } catch (std::runtime_error &e) {
           edm::LogError("EcalPedOffset") << e.what();
         }
       }
     }  // loop over the crystals
   }    // loop over the super-modules
  
   // insert the map of tables in the database
   if (DBconnection) {
     try {
       LogDebug("EcalPedOffset") << "Inserting dataset ... " << std::flush;
       if (!DBdataset.empty())
         DBconnection->insertDataSet(&DBdataset, &moniov);
       LogDebug("EcalPedOffset") << "done.";
     } catch (std::runtime_error &e) {
       edm::LogError("EcalPedOffset") << e.what();
     }
   }
  
   if (DBconnection) {
     delete DBconnection;
   }
 }

References MillePedeFileConverter_cfg::e, runTauDisplay::eid, PixelSLinkDataInputSource_cfi::fedid, EcalCondDBInterface::fetchMonRunIOV(), EcalCondDBInterface::fetchRunIOV(), EcalCondDBInterface::getEcalLogicID(), RunIOV::getRunTag(), EcalCondDBInterface::insertDataSet(), LogDebug, m_create_moniov, m_dbHostName, m_dbHostPort, m_dbName, m_dbPassword, m_dbUserName, m_location, m_pedResult, m_pedValues, m_run, mps_fire::result, writedatasetfile::run, test_db_connect::runtag, MonPedestalOffsetsDat::setDACG1(), MonPedestalOffsetsDat::setDACG12(), MonPedestalOffsetsDat::setDACG6(), MonRunTag::setGeneralTag(), LocationDef::setLocation(), MonVersionDef::setMonitoringVersion(), MonRunIOV::setMonRunTag(), MonRunTag::setMonVersionDef(), MonRunIOV::setRunIOV(), RunTypeDef::setRunType(), MonRunIOV::setSubRunNumber(), MonRunIOV::setSubRunStart(), MonPedestalOffsetsDat::setTaskStatus(), and Tm::setToCurrentGMTime().

Referenced by endJob().

◆ writeXMLFiles()

void EcalPedOffset::writeXMLFiles ( std::string fileName )

write the results into xml format

write the m_pedResults to XML files

Definition at line 385 of file EcalPedOffset.cc.

                                                     {
   // loop over the super-modules
   for (std::map<int, TPedResult *>::const_iterator smRes = m_pedResult.begin(); smRes != m_pedResult.end(); ++smRes) {
     std::string thisSMFileName = fileName;
     // open the output stream
     thisSMFileName += "_";
     thisSMFileName += intToString(smRes->first);
     thisSMFileName += ".xml";
     std::ofstream xml_outfile;
     xml_outfile.open(thisSMFileName.c_str());
  
     // write the header file
     xml_outfile << "<offsets>" << std::endl;
     xml_outfile << "<PEDESTAL_OFFSET_RELEASE VERSION_ID = \"SM1_VER1\"> \n";
     xml_outfile << "  <RELEASE_ID>RELEASE_1</RELEASE_ID>\n";
     xml_outfile << "  <SUPERMODULE>";
     xml_outfile << smRes->first;
     xml_outfile << "</SUPERMODULE>\n";
     xml_outfile << "  <TIME_STAMP> 070705 </TIME_STAMP>" << std::endl;
  
     // loop over the crystals
     for (int xtal = 0; xtal < 1700; ++xtal) {
       int crystalNumber = m_pedValues[smRes->first]->getCrystalNumber(xtal);
       if (crystalNumber == 0)
         continue;
       xml_outfile << "  <PEDESTAL_OFFSET>\n";
       xml_outfile << "    <HIGH>" << ((smRes->second)->m_DACvalue)[0][xtal] << "</HIGH>\n";
       xml_outfile << "    <MED>" << ((smRes->second)->m_DACvalue)[1][xtal] << "</MED>\n";
       xml_outfile << "    <LOW>" << ((smRes->second)->m_DACvalue)[2][xtal] << "</LOW>\n";
       xml_outfile << "    <CRYSTAL> " << crystalNumber << " </CRYSTAL>\n";
       xml_outfile << "  </PEDESTAL_OFFSET>" << std::endl;
     }
  
     // close the open tags
     xml_outfile << " </PEDESTAL_OFFSET_RELEASE>" << std::endl;
     xml_outfile << "</offsets>" << std::endl;
     xml_outfile.close();
   }  // loop over the super-modules
 }

References MillePedeFileConverter_cfg::fileName, intToString(), m_pedResult, m_pedValues, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endJob().

Member Data Documentation

◆ ecalElectronicsMap_

const EcalElectronicsMapping* EcalPedOffset::ecalElectronicsMap_

private

Definition at line 55 of file EcalPedOffset.h.

Referenced by beginRun(), and readDACs().

◆ m_barrelDigiCollection

edm::InputTag EcalPedOffset::m_barrelDigiCollection

private

secondary name given to collection of digis

Definition at line 61 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_bestPed

int EcalPedOffset::m_bestPed

private

Definition at line 73 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and readDACs().

◆ m_create_moniov

bool EcalPedOffset::m_create_moniov

private

allow the creation of a new moniov if not existing in the DB by default it is false.

Definition at line 87 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_DACmax

int EcalPedOffset::m_DACmax

private

Definition at line 71 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

◆ m_DACmin

int EcalPedOffset::m_DACmin

private

Definition at line 70 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

◆ m_dbHostName

std::string EcalPedOffset::m_dbHostName

private

database host name

Definition at line 76 of file EcalPedOffset.h.

Referenced by endJob(), and writeDb().

◆ m_dbHostPort

int EcalPedOffset::m_dbHostPort

private

database

Definition at line 84 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbName

std::string EcalPedOffset::m_dbName

private

database name

Definition at line 78 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbPassword

std::string EcalPedOffset::m_dbPassword

private

database user password

Definition at line 82 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbUserName

std::string EcalPedOffset::m_dbUserName

private

database user name

Definition at line 80 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_endcapDigiCollection

edm::InputTag EcalPedOffset::m_endcapDigiCollection

private

secondary name given to collection of digis

Definition at line 62 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_headerCollection

edm::InputTag EcalPedOffset::m_headerCollection

private

name of module/plugin/producer making headers

Definition at line 63 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_location

std::string EcalPedOffset::m_location

private

Definition at line 89 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_maxChi2OverNDFAllowed_

double EcalPedOffset::m_maxChi2OverNDFAllowed_

private

max chi2/ndf allowed for linearity test

Definition at line 100 of file EcalPedOffset.h.

Referenced by makePlots().

◆ m_maxSlopeAllowed_

double EcalPedOffset::m_maxSlopeAllowed_

private

max slope (in magnitude) allowed for linearity test

Definition at line 96 of file EcalPedOffset.h.

Referenced by makePlots().

◆ m_minSlopeAllowed_

double EcalPedOffset::m_minSlopeAllowed_

private

min slope (in magnitude) allowed for linearity test

Definition at line 98 of file EcalPedOffset.h.

Referenced by makePlots().

◆ m_pedResult

std::map<int, TPedResult *> EcalPedOffset::m_pedResult

private

Definition at line 68 of file EcalPedOffset.h.

Referenced by endJob(), writeDb(), writeXMLFiles(), and ~EcalPedOffset().

◆ m_pedValues

std::map<int, TPedValues *> EcalPedOffset::m_pedValues

private

Definition at line 67 of file EcalPedOffset.h.

Referenced by endJob(), makePlots(), readDACs(), writeDb(), writeXMLFiles(), and ~EcalPedOffset().

◆ m_plotting

std::string EcalPedOffset::m_plotting

private

the root file where to store the detail plots

Definition at line 94 of file EcalPedOffset.h.

Referenced by endJob(), and makePlots().

◆ m_RMSmax

double EcalPedOffset::m_RMSmax

private

Definition at line 72 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and readDACs().

◆ m_run

int EcalPedOffset::m_run

private

run number

Definition at line 91 of file EcalPedOffset.h.

Referenced by analyze(), and writeDb().

◆ m_xmlFile

std::string EcalPedOffset::m_xmlFile

private

name of the xml file to be saved

Definition at line 65 of file EcalPedOffset.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalPedOffset()

◆ ~EcalPedOffset()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ intToString()

◆ makePlots()

◆ readDACs() [1/2]

◆ readDACs() [2/2]

◆ writeDb()

◆ writeXMLFiles()

Member Data Documentation

◆ ecalElectronicsMap_

◆ m_barrelDigiCollection

◆ m_bestPed

◆ m_create_moniov

◆ m_DACmax

◆ m_DACmin

◆ m_dbHostName

◆ m_dbHostPort

◆ m_dbName

◆ m_dbPassword

◆ m_dbUserName

◆ m_endcapDigiCollection

◆ m_headerCollection

◆ m_location

◆ m_maxChi2OverNDFAllowed_

◆ m_maxSlopeAllowed_

◆ m_minSlopeAllowed_

◆ m_pedResult

◆ m_pedValues

◆ m_plotting

◆ m_RMSmax

◆ m_run

◆ m_xmlFile