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	endRun (edm::Run const &, edm::EventSetup const &) override
	EndRun. 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::one::EDAnalyzer< edm::one::WatchRuns >
	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
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_

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

const edm::EDGetTokenT< EBDigiCollection >	m_ebDigiToken

const edm::EDGetTokenT< EEDigiCollection >	m_eeDigiToken

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

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

std::string	m_location

const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	m_mappingToken

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

const edm::EDGetTokenT< EcalRawDataCollection >	m_rawDataToken

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

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 31 of file EcalPedOffset.cc.

References m_bestPed, m_DACmax, m_DACmin, and m_RMSmax.

     : m_barrelDigiCollection(paramSet.getParameter<edm::InputTag>("EBdigiCollection")),
       m_endcapDigiCollection(paramSet.getParameter<edm::InputTag>("EEdigiCollection")),
       m_headerCollection(paramSet.getParameter<edm::InputTag>("headerCollection")),
       m_rawDataToken(consumes<EcalRawDataCollection>(m_headerCollection)),
       m_ebDigiToken(consumes<EBDigiCollection>(m_barrelDigiCollection)),
       m_eeDigiToken(consumes<EEDigiCollection>(m_endcapDigiCollection)),
       m_mappingToken(esConsumes()),
       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;
 }

◆ ~EcalPedOffset()

EcalPedOffset::~EcalPedOffset ( )

override

Destructor.

dtor

Definition at line 64 of file EcalPedOffset.cc.

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

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

Member Function Documentation

◆ analyze()

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

overridevirtual

! Analyze

perform the analysis

Implements edm::one::EDAnalyzerBase.

Definition at line 90 of file EcalPedOffset.cc.

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

                                                                             {
   LogDebug("EcalPedOffset") << "entering analyze ...";
 
   // get the headers
   // (one header for each supermodule)
   edm::Handle<EcalRawDataCollection> DCCHeaders;
   event.getByToken(m_rawDataToken, 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.getByToken(m_ebDigiToken, 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.getByToken(m_eeDigiToken, 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!";
 }

◆ beginRun()

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

override

BeginRun.

begin the run

Definition at line 74 of file EcalPedOffset.cc.

References ecalElectronicsMap_, options_cfi::eventSetup, LogDebug, and m_mappingToken.

                                                                       {
   LogDebug("EcalPedOffset") << "entering beginRun...";
 
   // get the electronics map
   const auto mappingHandle = eventSetup.getHandle(m_mappingToken);
   ecalElectronicsMap_ = mappingHandle.product();
 }

◆ endJob()

void EcalPedOffset::endJob ( void )

overridevirtual

EndJob.

perform the minimization and write results

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 225 of file EcalPedOffset.cc.

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

Referenced by o2olib.O2ORunMgr::executeJob().

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

◆ endRun()

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

override

EndRun.

end the run

Definition at line 85 of file EcalPedOffset.cc.

85 {}

◆ intToString()

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

private

Definition at line 459 of file EcalPedOffset.cc.

References EgammaValidation_cff::num.

Referenced by writeXMLFiles().

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

◆ makePlots()

void EcalPedOffset::makePlots ( )

create the plots of the DAC pedestal trend

Definition at line 430 of file EcalPedOffset.cc.

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

Referenced by endJob().

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

◆ readDACs() [1/2]

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

private

Definition at line 156 of file EcalPedOffset.cc.

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

Referenced by analyze().

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

◆ readDACs() [2/2]

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

private

Definition at line 189 of file EcalPedOffset.cc.

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

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

◆ writeDb()

void EcalPedOffset::writeDb ( )

WriteDB.

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

Definition at line 243 of file EcalPedOffset.cc.

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

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

◆ writeXMLFiles()

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

write the results into xml format

write the m_pedResults to XML files

Definition at line 387 of file EcalPedOffset.cc.

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

Referenced by endJob().

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

Member Data Documentation

◆ ecalElectronicsMap_

const EcalElectronicsMapping* EcalPedOffset::ecalElectronicsMap_

private

Definition at line 58 of file EcalPedOffset.h.

Referenced by beginRun(), and readDACs().

◆ m_barrelDigiCollection

const edm::InputTag EcalPedOffset::m_barrelDigiCollection

private

secondary name given to collection of digis

Definition at line 64 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_bestPed

int EcalPedOffset::m_bestPed

private

Definition at line 81 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 95 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_DACmax

int EcalPedOffset::m_DACmax

private

Definition at line 79 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

◆ m_DACmin

int EcalPedOffset::m_DACmin

private

Definition at line 78 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

◆ m_dbHostName

std::string EcalPedOffset::m_dbHostName

private

database host name

Definition at line 84 of file EcalPedOffset.h.

Referenced by endJob(), and writeDb().

◆ m_dbHostPort

int EcalPedOffset::m_dbHostPort

private

database

Definition at line 92 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbName

std::string EcalPedOffset::m_dbName

private

database name

Definition at line 86 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbPassword

std::string EcalPedOffset::m_dbPassword

private

database user password

Definition at line 90 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_dbUserName

std::string EcalPedOffset::m_dbUserName

private

database user name

Definition at line 88 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_ebDigiToken

const edm::EDGetTokenT<EBDigiCollection> EcalPedOffset::m_ebDigiToken

private

Definition at line 69 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_eeDigiToken

const edm::EDGetTokenT<EEDigiCollection> EcalPedOffset::m_eeDigiToken

private

Definition at line 70 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_endcapDigiCollection

const edm::InputTag EcalPedOffset::m_endcapDigiCollection

private

secondary name given to collection of digis

Definition at line 65 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_headerCollection

const edm::InputTag EcalPedOffset::m_headerCollection

private

name of module/plugin/producer making headers

Definition at line 66 of file EcalPedOffset.h.

◆ m_location

std::string EcalPedOffset::m_location

private

Definition at line 97 of file EcalPedOffset.h.

Referenced by writeDb().

◆ m_mappingToken

const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalPedOffset::m_mappingToken

private

Definition at line 71 of file EcalPedOffset.h.

Referenced by beginRun().

◆ m_maxChi2OverNDFAllowed_

double EcalPedOffset::m_maxChi2OverNDFAllowed_

private

max chi2/ndf allowed for linearity test

Definition at line 108 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 104 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 106 of file EcalPedOffset.h.

Referenced by makePlots().

◆ m_pedResult

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

private

Definition at line 76 of file EcalPedOffset.h.

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

◆ m_pedValues

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

private

Definition at line 75 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 102 of file EcalPedOffset.h.

Referenced by endJob(), and makePlots().

◆ m_rawDataToken

const edm::EDGetTokenT<EcalRawDataCollection> EcalPedOffset::m_rawDataToken

private

Definition at line 68 of file EcalPedOffset.h.

Referenced by analyze().

◆ m_RMSmax

double EcalPedOffset::m_RMSmax

private

Definition at line 80 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and readDACs().

◆ m_run

int EcalPedOffset::m_run

private

run number

Definition at line 99 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 73 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()

◆ endRun()

◆ 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_ebDigiToken

◆ m_eeDigiToken

◆ m_endcapDigiCollection

◆ m_headerCollection

◆ m_location

◆ m_mappingToken

◆ m_maxChi2OverNDFAllowed_

◆ m_maxSlopeAllowed_

◆ m_minSlopeAllowed_

◆ m_pedResult

◆ m_pedValues

◆ m_plotting

◆ m_rawDataToken

◆ m_RMSmax

◆ m_run

◆ m_xmlFile