Provides LHCInfo data necessary for CTPPS reconstruction (and direct simulation). More...

Inheritance diagram for CTPPSLHCInfoRandomXangleESSource:

Classes
struct	BinData

Public Member Functions
	CTPPSLHCInfoRandomXangleESSource (const edm::ParameterSet &)

edm::ESProducts< std::unique_ptr< LHCInfo > >	produce (const LHCInfoRcd &)

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

	ESProducer (const ESProducer &)=delete

ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const

ESProducer const &	operator= (const ESProducer &)=delete

template<typename Record >
void	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final

	~ESProducer () noexcept(false) override

Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()

	ESProxyFactoryProducer (const ESProxyFactoryProducer &)=delete

const ESProxyFactoryProducer &	operator= (const ESProxyFactoryProducer &)=delete

	~ESProxyFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
void	createKeyedProxies (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

	DataProxyProvider ()

	DataProxyProvider (const DataProxyProvider &)=delete

const ComponentDescription &	description () const

void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const

virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)

bool	isUsingRecord (const EventSetupRecordKey &key) const

KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)

const DataProxyProvider &	operator= (const DataProxyProvider &)=delete

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~DataProxyProvider () noexcept(false)

Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
bool	concurrentFinder () const

const eventsetup::ComponentDescription &	descriptionForFinder () const

	EventSetupRecordIntervalFinder ()

	EventSetupRecordIntervalFinder (const EventSetupRecordIntervalFinder &)=delete

std::set< eventsetup::EventSetupRecordKey >	findingForRecords () const

const ValidityInterval &	findIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)

bool	nonconcurrentAndIOVNeedsUpdate (const eventsetup::EventSetupRecordKey &key, const IOVSyncValue &syncValue) const

const EventSetupRecordIntervalFinder &	operator= (const EventSetupRecordIntervalFinder &)=delete

void	resetInterval (const eventsetup::EventSetupRecordKey &)

void	setDescriptionForFinder (const eventsetup::ComponentDescription &iDescription)

virtual	~EventSetupRecordIntervalFinder () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &)

Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
void	setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &) override

Private Attributes
std::vector< BinData >	binData

double	m_beamEnergy

double	m_betaStar

std::unique_ptr< CLHEP::HepRandomEngine >	m_engine

unsigned int	m_generateEveryNEvents

std::string	m_label

Additional Inherited Members
Protected Types inherited from edm::ESProxyFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey

Protected Types inherited from edm::eventsetup::DataProxyProvider
using	KeyedProxiesVector = std::vector< std::pair< DataKey, std::shared_ptr< DataProxy > >>

Protected Member Functions inherited from edm::ESProducer
template<typename T >
auto	setWhatProduced (T iThis, const char iLabel)

template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})

template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const es::Label &iLabel={})

template<typename T , typename TReturn , typename TRecord , typename TArg >
ESConsumesCollectorT< TRecord >	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel={})

Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())

KeyedProxiesVector	registerProxies (const EventSetupRecordKey &, unsigned int iovIndex) override

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &key)

Protected Member Functions inherited from edm::EventSetupRecordIntervalFinder
template<class T >
void	findingRecord ()

void	findingRecordWithKey (const eventsetup::EventSetupRecordKey &)

Detailed Description

Provides LHCInfo data necessary for CTPPS reconstruction (and direct simulation).

Definition at line 26 of file CTPPSLHCInfoRandomXangleESSource.cc.

Constructor & Destructor Documentation

◆ CTPPSLHCInfoRandomXangleESSource()

CTPPSLHCInfoRandomXangleESSource::CTPPSLHCInfoRandomXangleESSource ( const edm::ParameterSet & conf )

Definition at line 56 of file CTPPSLHCInfoRandomXangleESSource.cc.

     : m_label(conf.getParameter<std::string>("label")),
  
       m_generateEveryNEvents(conf.getParameter<unsigned int>("generateEveryNEvents")),
  
       m_beamEnergy(conf.getParameter<double>("beamEnergy")),
       m_betaStar(conf.getParameter<double>("betaStar")),
  
       m_engine(new CLHEP::HepJamesRandom(conf.getParameter<unsigned int>("seed"))) {
   const auto &xangleHistogramFile = conf.getParameter<std::string>("xangleHistogramFile");
   const auto &xangleHistogramObject = conf.getParameter<std::string>("xangleHistogramObject");
  
   TFile *f_in = TFile::Open(xangleHistogramFile.c_str());
   TH1D *h_xangle = (TH1D *)f_in->Get(xangleHistogramObject.c_str());
  
   double s = 0.;
   for (int bi = 1; bi <= h_xangle->GetNbinsX(); ++bi)
     s += h_xangle->GetBinContent(bi);
  
   double cw = 0.;
   for (int bi = 1; bi <= h_xangle->GetNbinsX(); ++bi) {
     double xangle = h_xangle->GetBinCenter(bi);
     double w = h_xangle->GetBinContent(bi) / s;
  
     binData.push_back({cw, cw + w, xangle});
  
     cw += w;
   }
  
   delete f_in;
  
   setWhatProduced(this, m_label);
   findingRecord<LHCInfoRcd>();
 }

References binData, edm::ParameterSet::getParameter(), m_label, alignCSCRings::s, edm::ESProducer::setWhatProduced(), AlCaHLTBitMon_QueryRunRegistry::string, w, and OpticalFunctionsConfig_cfi::xangle.

Member Function Documentation

◆ fillDescriptions()

void CTPPSLHCInfoRandomXangleESSource::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 93 of file CTPPSLHCInfoRandomXangleESSource.cc.

                                                                                                   {
   edm::ParameterSetDescription desc;
  
   desc.add<std::string>("label", "")->setComment("label of the LHCInfo record");
  
   desc.add<unsigned int>("seed", 1)->setComment("random seed");
  
   desc.add<unsigned int>("generateEveryNEvents", 1)->setComment("how often to generate new xangle");
  
   desc.add<std::string>("xangleHistogramFile", "")->setComment("ROOT file with xangle distribution");
   desc.add<std::string>("xangleHistogramObject", "")->setComment("xangle distribution object in the ROOT file");
  
   desc.add<double>("beamEnergy", 0.)->setComment("beam energy");
   desc.add<double>("betaStar", 0.)->setComment("beta*");
  
   descriptions.add("ctppsLHCInfoRandomXangleESSource", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterDescriptionNode::setComment(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ produce()

edm::ESProducts< std::unique_ptr< LHCInfo > > CTPPSLHCInfoRandomXangleESSource::produce ( const LHCInfoRcd & )

Definition at line 123 of file CTPPSLHCInfoRandomXangleESSource.cc.

                                                                                                   {
   auto output = std::make_unique<LHCInfo>();
  
   const double u = CLHEP::RandFlat::shoot(m_engine.get(), 0., 1.);
  
   double xangle = 0.;
   for (const auto &d : binData) {
     if (d.min <= u && u <= d.max) {
       xangle = d.xangle;
       break;
     }
   }
  
   output->setEnergy(m_beamEnergy);
   output->setCrossingAngle(xangle);
   output->setBetaStar(m_betaStar);
  
   return edm::es::products(std::move(output));
 }

References binData, ztail::d, m_beamEnergy, m_betaStar, m_engine, eostools::move(), convertSQLitetoXML_cfg::output, edm::es::products(), and OpticalFunctionsConfig_cfi::xangle.

◆ setIntervalFor()

void CTPPSLHCInfoRandomXangleESSource::setIntervalFor	(	const edm::eventsetup::EventSetupRecordKey &	key,
		const edm::IOVSyncValue &	iosv,
		edm::ValidityInterval &	oValidity
	)

overrideprivatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 113 of file CTPPSLHCInfoRandomXangleESSource.cc.

                                                                                       {
   edm::EventID beginEvent = iosv.eventID();
   edm::EventID endEvent(beginEvent.run(), beginEvent.luminosityBlock(), beginEvent.event() + m_generateEveryNEvents);
   oValidity = edm::ValidityInterval(edm::IOVSyncValue(beginEvent), edm::IOVSyncValue(endEvent));
 }