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)
	! Analyze More...

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

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

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

virtual	~EcalPedOffset ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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::vector< ModuleDescription const * > &modules, 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	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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 ( const edm::ParameterSet & ps )

Constructor.

ctor

Definition at line 38 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_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 ( )

virtual

Destructor.

dtor

Definition at line 72 of file EcalPedOffset.cc.

References 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

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

virtual

! Analyze

perform the analysis

Implements edm::EDAnalyzer.

Definition at line 104 of file EcalPedOffset.cc.

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

 {
   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->size()==0) 
   {
     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->size()==0) 
   {
     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!";
   
 }

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

virtual

BeginRun.

begin the run

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file EcalPedOffset.cc.

References ecalElectronicsMap_, edm::EventSetup::get(), cmsBatch::handle, LogDebug, and edm::ESHandle< T >::product().

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

void EcalPedOffset::endJob ( void )

virtual

EndJob.

perform the minimization and write results

Reimplemented from edm::EDAnalyzer.

Definition at line 273 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 ();       
 }

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

private

Definition at line 547 of file EcalPedOffset.cc.

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
 }

void EcalPedOffset::makePlots ( )

create the plots of the DAC pedestal trend

Definition at line 513 of file EcalPedOffset.cc.

References LogDebug, m_maxChi2OverNDFAllowed_, m_maxSlopeAllowed_, m_minSlopeAllowed_, m_pedValues, m_plotting, and indexGen::rootFile.

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

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

private

Definition at line 182 of file EcalPedOffset.cc.

References ecalMGPA::adc(), edm::DataFrameContainer::begin(), EcalElectronicsId::dccId(), ecalElectronicsMap_, edm::DataFrameContainer::end(), ecalMGPA::gainId(), EcalElectronicsMapping::getElectronicsId(), EBDetId::ic(), LogDebug, m_bestPed, m_pedValues, m_RMSmax, EcalDataFrame::MAXSAMPLES, and simplePhotonAnalyzer_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
 }

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

private

Definition at line 226 of file EcalPedOffset.cc.

References ecalMGPA::adc(), edm::DataFrameContainer::begin(), EcalElectronicsId::dccId(), ecalElectronicsMap_, edm::DataFrameContainer::end(), ecalMGPA::gainId(), EcalElectronicsMapping::getElectronicsId(), LogDebug, m_bestPed, m_pedValues, m_RMSmax, EcalDataFrame::MAXSAMPLES, and simplePhotonAnalyzer_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
 }

void EcalPedOffset::writeDb ( )

WriteDB.

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

Definition at line 296 of file EcalPedOffset.cc.

References MillePedeFileConverter_cfg::e, runTauDisplay::eid, 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, findQualityFiles::run, MonPedestalOffsetsDat::setDACG1(), MonPedestalOffsetsDat::setDACG12(), MonPedestalOffsetsDat::setDACG6(), RunTag::setGeneralTag(), MonRunTag::setGeneralTag(), LocationDef::setLocation(), RunTag::setLocationDef(), MonVersionDef::setMonitoringVersion(), MonRunIOV::setMonRunTag(), MonRunTag::setMonVersionDef(), MonRunIOV::setRunIOV(), RunTypeDef::setRunType(), RunTag::setRunTypeDef(), 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.size() != 0 )
     {
       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 (1); //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.size() != 0 ) DBconnection->insertDataSet (&DBdataset, &moniov);
       LogDebug ("EcalPedOffset") << "done." ;
     } catch (std::runtime_error &e) {
       edm::LogError ("EcalPedOffset") << e.what ();
     }
   }
 
   if ( DBconnection ) {delete DBconnection;}
 }

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

write the results into xml format

write the m_pedResults to XML files

Definition at line 461 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

const EcalElectronicsMapping* EcalPedOffset::ecalElectronicsMap_

private

Definition at line 60 of file EcalPedOffset.h.

Referenced by beginRun(), and readDACs().

edm::InputTag EcalPedOffset::m_barrelDigiCollection

private

secondary name given to collection of digis

Definition at line 66 of file EcalPedOffset.h.

Referenced by analyze().

int EcalPedOffset::m_bestPed

private

Definition at line 78 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and readDACs().

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 92 of file EcalPedOffset.h.

Referenced by writeDb().

int EcalPedOffset::m_DACmax

private

Definition at line 76 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

int EcalPedOffset::m_DACmin

private

Definition at line 75 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and endJob().

std::string EcalPedOffset::m_dbHostName

private

database host name

Definition at line 81 of file EcalPedOffset.h.

Referenced by endJob(), and writeDb().

int EcalPedOffset::m_dbHostPort

private

database

Definition at line 89 of file EcalPedOffset.h.

Referenced by writeDb().

std::string EcalPedOffset::m_dbName

private

database name

Definition at line 83 of file EcalPedOffset.h.

Referenced by writeDb().

std::string EcalPedOffset::m_dbPassword

private

database user password

Definition at line 87 of file EcalPedOffset.h.

Referenced by writeDb().

std::string EcalPedOffset::m_dbUserName

private

database user name

Definition at line 85 of file EcalPedOffset.h.

Referenced by writeDb().

edm::InputTag EcalPedOffset::m_endcapDigiCollection

private

secondary name given to collection of digis

Definition at line 67 of file EcalPedOffset.h.

Referenced by analyze().

edm::InputTag EcalPedOffset::m_headerCollection

private

name of module/plugin/producer making headers

Definition at line 68 of file EcalPedOffset.h.

Referenced by analyze().

std::string EcalPedOffset::m_location

private

Definition at line 94 of file EcalPedOffset.h.

Referenced by writeDb().

double EcalPedOffset::m_maxChi2OverNDFAllowed_

private

max chi2/ndf allowed for linearity test

Definition at line 105 of file EcalPedOffset.h.

Referenced by makePlots().

double EcalPedOffset::m_maxSlopeAllowed_

private

max slope (in magnitude) allowed for linearity test

Definition at line 101 of file EcalPedOffset.h.

Referenced by makePlots().

double EcalPedOffset::m_minSlopeAllowed_

private

min slope (in magnitude) allowed for linearity test

Definition at line 103 of file EcalPedOffset.h.

Referenced by makePlots().

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

private

Definition at line 73 of file EcalPedOffset.h.

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

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

private

Definition at line 72 of file EcalPedOffset.h.

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

std::string EcalPedOffset::m_plotting

private

the root file where to store the detail plots

Definition at line 99 of file EcalPedOffset.h.

Referenced by endJob(), and makePlots().

double EcalPedOffset::m_RMSmax

private

Definition at line 77 of file EcalPedOffset.h.

Referenced by EcalPedOffset(), and readDACs().

int EcalPedOffset::m_run

private

run number

Definition at line 96 of file EcalPedOffset.h.

Referenced by analyze(), and writeDb().

std::string EcalPedOffset::m_xmlFile

private

name of the xml file to be saved

Definition at line 70 of file EcalPedOffset.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation