Inheritance diagram for CTPPSInterpolatedOpticalFunctionsESSource:

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

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection >	produce (const CTPPSInterpolatedOpticsRcd &)

	~CTPPSInterpolatedOpticalFunctionsESSource () override

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)

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

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

Private Attributes
float	currentCrossingAngle_

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection >	currentData_

bool	currentDataValid_

edm::ESGetToken< LHCInfoPerFill, LHCInfoPerFillRcd >	lhcInfoPerFillToken_

edm::ESGetToken< LHCInfoPerLS, LHCInfoPerLSRcd >	lhcInfoPerLSToken_

edm::ESGetToken< LHCInfo, LHCInfoRcd >	lhcInfoToken_

edm::ESGetToken< LHCOpticalFunctionsSetCollection, CTPPSOpticsRcd >	opticsToken_

const bool	useNewLHCInfo_

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)

Detailed Description

Definition at line 20 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

Constructor & Destructor Documentation

◆ CTPPSInterpolatedOpticalFunctionsESSource()

CTPPSInterpolatedOpticalFunctionsESSource::CTPPSInterpolatedOpticalFunctionsESSource ( const edm::ParameterSet & iConfig )

Definition at line 43 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References gpuPixelDoublets::cc, edm::ParameterSet::getParameter(), lhcInfoPerFillToken_, lhcInfoPerLSToken_, lhcInfoToken_, opticsToken_, edm::ESProducer::setWhatProduced(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : currentCrossingAngle_(LHCInfoCombined::crossingAngleInvalid),
       currentDataValid_(false),
       useNewLHCInfo_(iConfig.getParameter<bool>("useNewLHCInfo")) {
   auto cc = setWhatProduced(this, iConfig.getParameter<std::string>("opticsLabel"));
   opticsToken_ = cc.consumes(edm::ESInputTag("", iConfig.getParameter<std::string>("opticsLabel")));
   lhcInfoToken_ = cc.consumes(edm::ESInputTag("", iConfig.getParameter<std::string>("lhcInfoLabel")));
   lhcInfoPerLSToken_ = cc.consumes(edm::ESInputTag("", iConfig.getParameter<std::string>("lhcInfoPerLSLabel")));
   lhcInfoPerFillToken_ = cc.consumes(edm::ESInputTag("", iConfig.getParameter<std::string>("lhcInfoPerFillLabel")));
 }

◆ ~CTPPSInterpolatedOpticalFunctionsESSource()

CTPPSInterpolatedOpticalFunctionsESSource::~CTPPSInterpolatedOpticalFunctionsESSource ( )

inlineoverride

Definition at line 23 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

23 {}

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 56 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

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

                                                                                                            {
   edm::ParameterSetDescription desc;
 
   desc.add<std::string>("lhcInfoLabel", "")->setComment("label of the LHCInfo record");
   desc.add<std::string>("lhcInfoPerFillLabel", "")->setComment("label of the LHCInfoPerFill record");
   desc.add<std::string>("lhcInfoPerLSLabel", "")->setComment("label of the LHCInfoPerLS record");
   desc.add<std::string>("opticsLabel", "")->setComment("label of the optics records");
   desc.add<bool>("useNewLHCInfo", false)->setComment("flag whether to use new LHCInfoPer* records or old LHCInfo");
 
   descriptions.add("ctppsInterpolatedOpticalFunctionsESSourceDefault", desc);
 }

◆ produce()

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection > CTPPSInterpolatedOpticalFunctionsESSource::produce ( const CTPPSInterpolatedOpticsRcd & iRecord )

Definition at line 70 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References LHCInfoCombined::createLHCInfoCombined(), LHCInfoCombined::crossingAngleInvalid, currentCrossingAngle_, currentData_, currentDataValid_, MillePedeFileConverter_cfg::e, Exception, dqmdumpme::first, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), LHCInterpolatedOpticalFunctionsSet::initializeSplines(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, lhcInfoPerFillToken_, lhcInfoPerLSToken_, lhcInfoToken_, LHCOpticalFunctionsSet::m_fcn_values, LHCOpticalFunctionsSet::m_xi_values, LHCOpticalFunctionsSet::m_z, eostools::move(), GetRecoTauVFromDQM_MC_cff::next, LHCOpticalFunctionsSet::nFunctions, opticsToken_, pi, useNewLHCInfo_, and protons_cff::xi.

                                                {
   // get the input data
   LHCOpticalFunctionsSetCollection const &ofColl = iRecord.get(opticsToken_);
 
   auto lhcInfoCombined = LHCInfoCombined::createLHCInfoCombined<
       CTPPSInterpolatedOpticsRcd,
       edm::mpl::Vector<CTPPSOpticsRcd, LHCInfoRcd, LHCInfoPerFillRcd, LHCInfoPerLSRcd>>(
       iRecord, lhcInfoPerLSToken_, lhcInfoPerFillToken_, lhcInfoToken_, useNewLHCInfo_);
 
   // is there anything to do?
   if (currentDataValid_ && lhcInfoCombined.crossingAngle() == currentCrossingAngle_)
     return currentData_;
 
   // is crossing angle reasonable (LHCInfo is correctly filled in DB)?
   if (lhcInfoCombined.isCrossingAngleInvalid()) {
     edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource")
         << "Invalid crossing angle, no optical functions produced.";
 
     currentDataValid_ = false;
     currentCrossingAngle_ = LHCInfoCombined::crossingAngleInvalid;
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
 
     return currentData_;
   }
 
   // set new crossing angle
   currentCrossingAngle_ = lhcInfoCombined.crossingAngle();
   edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource")
       << "Crossing angle has changed to " << currentCrossingAngle_ << ".";
 
   // is input optics available ?
   if (ofColl.empty()) {
     edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource")
         << "No input optics available, no optical functions produced.";
 
     currentDataValid_ = false;
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
 
     return currentData_;
   }
 
   // regular case with single-xangle input
   if (ofColl.size() == 1) {
     const auto &it = ofColl.begin();
 
     // does the input xangle correspond to the actual one?
     if (fabs(currentCrossingAngle_ - it->first) > 1e-6)
       throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource")
           << "Cannot interpolate: input given only for xangle " << it->first << " while interpolation requested for "
           << currentCrossingAngle_ << ".";
 
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
     for (const auto &rp_p : it->second) {
       const auto rpId = rp_p.first;
       LHCInterpolatedOpticalFunctionsSet iof(rp_p.second);
       iof.initializeSplines();
       currentData_->emplace(rpId, std::move(iof));
     }
 
     currentDataValid_ = true;
   }
 
   // regular case with multi-xangle input
   if (ofColl.size() > 1) {
     // find the closest xangle points for interpolation
     auto it1 = ofColl.begin();
     auto it2 = std::next(it1);
 
     if (currentCrossingAngle_ > it1->first) {
       for (; it1 != ofColl.end(); ++it1) {
         it2 = std::next(it1);
 
         if (it2 == ofColl.end()) {
           it2 = it1;
           it1 = std::prev(it1);
           break;
         }
 
         if (it1->first <= currentCrossingAngle_ && currentCrossingAngle_ < it2->first)
           break;
       }
     }
 
     const auto &xangle1 = it1->first;
     const auto &xangle2 = it2->first;
 
     const auto &ofs1 = it1->second;
     const auto &ofs2 = it2->second;
 
     // do the interpoaltion RP by RP
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
     for (const auto &rp_p : ofs1) {
       const auto rpId = rp_p.first;
       const auto &rp_it2 = ofs2.find(rpId);
       if (rp_it2 == ofs2.end())
         throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "RP mismatch between ofs1 and ofs2.";
 
       const auto &of1 = rp_p.second;
       const auto &of2 = rp_it2->second;
 
       const size_t num_xi_vals1 = of1.getXiValues().size();
       const size_t num_xi_vals2 = of2.getXiValues().size();
 
       if (num_xi_vals1 != num_xi_vals2)
         throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "Size mismatch between ofs1 and ofs2.";
 
       const size_t num_xi_vals = num_xi_vals1;
 
       LHCInterpolatedOpticalFunctionsSet iof;
       iof.m_z = of1.getScoringPlaneZ();
       iof.m_fcn_values.resize(LHCInterpolatedOpticalFunctionsSet::nFunctions);
       iof.m_xi_values.resize(num_xi_vals);
 
       for (size_t fi = 0; fi < of1.getFcnValues().size(); ++fi) {
         iof.m_fcn_values[fi].resize(num_xi_vals);
 
         for (size_t pi = 0; pi < num_xi_vals; ++pi) {
           double xi = of1.getXiValues()[pi];
           double xi_control = of2.getXiValues()[pi];
 
           if (fabs(xi - xi_control) > 1e-6)
             throw cms::Exception("CTPPSInterpolatedOpticalFunctionsESSource") << "Xi mismatch between ofs1 and ofs2.";
 
           iof.m_xi_values[pi] = xi;
 
           double v1 = of1.getFcnValues()[fi][pi];
           double v2 = of2.getFcnValues()[fi][pi];
           iof.m_fcn_values[fi][pi] = v1 + (v2 - v1) / (xangle2 - xangle1) * (currentCrossingAngle_ - xangle1);
         }
       }
 
       iof.initializeSplines();
 
       currentDataValid_ = true;
       currentData_->emplace(rpId, std::move(iof));
     }
   }
 
   return currentData_;
 }