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

	ESProducer (ESProducer &&)=delete

ESResolverIndex const *	getTokenIndices (unsigned int iIndex) const

ESRecordIndex const *	getTokenRecordIndices (unsigned int iIndex) const

bool	hasMayConsumes () const noexcept

size_t	numberOfTokenIndices (unsigned int iIndex) const

ESProducer &	operator= (const ESProducer &)=delete

ESProducer &	operator= (ESProducer &&)=delete

SerialTaskQueueChain &	queue ()

template<typename Record >
std::optional< std::vector< ESResolverIndex > >	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord) const

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &) final

	~ESProducer () noexcept(false) override

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

	ESProductResolverFactoryProducer (const ESProductResolverFactoryProducer &)=delete

const ESProductResolverFactoryProducer &	operator= (const ESProductResolverFactoryProducer &)=delete

	~ESProductResolverFactoryProducer () noexcept(false) override

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

const ComponentDescription &	description () const

	ESProductResolverProvider ()

	ESProductResolverProvider (const ESProductResolverProvider &)=delete

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

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

bool	isUsingRecord (const EventSetupRecordKey &key) const

KeyedResolvers &	keyedResolvers (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)

const ESProductResolverProvider &	operator= (const ESProductResolverProvider &)=delete

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~ESProductResolverProvider () 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::ESProductResolverProvider
static void	prevalidate (ConfigurationDescriptions &)

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

Private Attributes
double	m_beamEnergy

std::unique_ptr< CLHEP::HepRandomEngine >	m_engine

unsigned int	m_generateEveryNEvents

std::string	m_label

std::vector< BinData >	xangleBetaStarBins

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

Protected Types inherited from edm::eventsetup::ESProductResolverProvider
using	KeyedResolversVector = std::vector< std::pair< DataKey, std::shared_ptr< ESProductResolver > >>

Protected Member Functions inherited from edm::ESProducer
ESConsumesInfo *	consumesInfoPushBackNew ()

unsigned int	consumesInfoSize () const

template<typename CallbackT , typename TList , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>> iCallback, const TList , const TRecord iRecord, const es::Label &iLabel)

template<typename CallbackT , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>>, const eventsetup::produce::Null , const TRecord , const es::Label &)

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

template<typename T >
auto	setWhatProduced (T iThis, const char 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 TDecorator >
auto	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const TDecorator &iDec, const es::Label &iLabel={})

template<typename TFunc >
auto	setWhatProduced (TFunc &&func, const es::Label &iLabel={})

template<typename TReturn , typename TRecord , typename TFunc , typename TDecorator >
ESConsumesCollectorT< TRecord >	setWhatProduced (TFunc &&func, TDecorator &&iDec, const es::Label &iLabel={})

void	usesResources (std::vector< std::string > const &)

Protected Member Functions inherited from edm::ESProductResolverFactoryProducer
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::ESProductResolverFactoryBase > iFactory, const std::string &iLabel=std::string())

KeyedResolversVector	registerResolvers (const EventSetupRecordKey &, unsigned int iovIndex) override

Protected Member Functions inherited from edm::eventsetup::ESProductResolverProvider
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.

References simBeamSpotPI::betaStar, DummyCfis::c, Exception, edm::ParameterSet::getParameter(), m_label, edm::ESProducer::setWhatProduced(), AlCaHLTBitMon_QueryRunRegistry::string, x, profile_base_cff::xangle, xangleBetaStarBins, profile_base_cff::xangleBetaStarHistogramFile, profile_base_cff::xangleBetaStarHistogramObject, and y.

     : m_label(conf.getParameter<std::string>("label")),
 
       m_generateEveryNEvents(conf.getParameter<unsigned int>("generateEveryNEvents")),
 
       m_beamEnergy(conf.getParameter<double>("beamEnergy")),
 
       m_engine(new CLHEP::HepJamesRandom(conf.getParameter<unsigned int>("seed"))) {
   // get input beta* vs. xangle histogram
   const auto &xangleBetaStarHistogramFile = conf.getParameter<std::string>("xangleBetaStarHistogramFile");
   const auto &xangleBetaStarHistogramObject = conf.getParameter<std::string>("xangleBetaStarHistogramObject");
 
   edm::FileInPath fip(xangleBetaStarHistogramFile);
   std::unique_ptr<TFile> f_in(TFile::Open(fip.fullPath().c_str()));
   if (!f_in)
     throw cms::Exception("PPS") << "Cannot open input file '" << xangleBetaStarHistogramFile << "'.";
 
   TH2D *h_xangle_beta_star = (TH2D *)f_in->Get(xangleBetaStarHistogramObject.c_str());
   if (!h_xangle_beta_star)
     throw cms::Exception("PPS") << "Cannot load input object '" << xangleBetaStarHistogramObject << "'.";
 
   // parse histogram
   double sum = 0.;
   for (int x = 1; x <= h_xangle_beta_star->GetNbinsX(); ++x) {
     for (int y = 1; y <= h_xangle_beta_star->GetNbinsY(); ++y)
       sum += h_xangle_beta_star->GetBinContent(x, y);
   }
 
   double cw = 0;
   for (int x = 1; x <= h_xangle_beta_star->GetNbinsX(); ++x) {
     for (int y = 1; y <= h_xangle_beta_star->GetNbinsY(); ++y) {
       const double c = h_xangle_beta_star->GetBinContent(x, y);
       const double xangle = h_xangle_beta_star->GetXaxis()->GetBinCenter(x);
       const double betaStar = h_xangle_beta_star->GetYaxis()->GetBinCenter(y);
 
       if (c > 0.) {
         const double rc = c / sum;
         xangleBetaStarBins.push_back({cw, cw + rc, xangle, betaStar});
         cw += rc;
       }
     }
   }
 
   setWhatProduced(this, m_label);
   findingRecord<LHCInfoRcd>();
 }

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 105 of file CTPPSLHCInfoRandomXangleESSource.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                   {
   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>("xangleBetaStarHistogramFile", "")->setComment("ROOT file with xangle distribution");
   desc.add<std::string>("xangleBetaStarHistogramObject", "")->setComment("xangle distribution object in the ROOT file");
 
   desc.add<double>("beamEnergy", 0.)->setComment("beam energy");
 
   descriptions.add("ctppsLHCInfoRandomXangleESSource", desc);
 }

◆ produce()

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

Definition at line 134 of file CTPPSLHCInfoRandomXangleESSource.cc.

References simBeamSpotPI::betaStar, ztail::d, m_beamEnergy, m_engine, eostools::move(), edm::es::products(), profile_base_cff::xangle, and xangleBetaStarBins.

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

◆ setIntervalFor()

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

overrideprivatevirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 124 of file CTPPSLHCInfoRandomXangleESSource.cc.

References edm::EventID::event(), edm::IOVSyncValue::eventID(), edm::EventID::luminosityBlock(), m_generateEveryNEvents, and edm::EventID::run().

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