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

ESProxyIndex 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 const &	operator= (const ESProducer &)=delete

SerialTaskQueueChain &	queue ()

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

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)

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

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

Private Attributes
float	currentCrossingAngle_

std::shared_ptr< LHCInterpolatedOpticalFunctionsSetCollection >	currentData_

bool	currentDataValid_

edm::ESGetToken< LHCInfo, LHCInfoRcd >	lhcInfoToken_

edm::ESGetToken< LHCOpticalFunctionsSetCollection, CTPPSOpticsRcd >	opticsToken_

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 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 TArg >
auto	setWhatProduced (T iThis, TReturn(T ::iMethod)(const TRecord &), const TArg &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::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)

Detailed Description

Definition at line 19 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

Constructor & Destructor Documentation

◆ CTPPSInterpolatedOpticalFunctionsESSource()

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

Definition at line 39 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

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

     : currentCrossingAngle_(-1.), currentDataValid_(false) {
   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")));
 }

◆ ~CTPPSInterpolatedOpticalFunctionsESSource()

CTPPSInterpolatedOpticalFunctionsESSource::~CTPPSInterpolatedOpticalFunctionsESSource ( )

inlineoverride

Definition at line 22 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

22 {}

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 48 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>("opticsLabel", "")->setComment("label of the optics records");
 
   descriptions.add("ctppsInterpolatedOpticalFunctionsESSource", desc);
 }

◆ produce()

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

Definition at line 59 of file CTPPSInterpolatedOpticalFunctionsESSource.cc.

References LHCInfo::crossingAngle(), currentCrossingAngle_, currentData_, currentDataValid_, MillePedeFileConverter_cfg::e, Exception, first, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), LHCInterpolatedOpticalFunctionsSet::initializeSplines(), lhcInfoToken_, LHCOpticalFunctionsSet::m_fcn_values, LHCOpticalFunctionsSet::m_xi_values, LHCOpticalFunctionsSet::m_z, eostools::move(), GetRecoTauVFromDQM_MC_cff::next, LHCOpticalFunctionsSet::nFunctions, opticsToken_, pi, and protons_cff::xi.

                                                {
   // get the input data
   LHCOpticalFunctionsSetCollection const &ofColl = iRecord.get(opticsToken_);
 
   LHCInfo const &lhcInfo = iRecord.get(lhcInfoToken_);
 
   // is there anything to do?
   if (currentDataValid_ && lhcInfo.crossingAngle() == currentCrossingAngle_)
     return currentData_;
 
   // is crossing angle reasonable (LHCInfo is correctly filled in DB)?
   if (lhcInfo.crossingAngle() == 0.) {
     edm::LogInfo("CTPPSInterpolatedOpticalFunctionsESSource")
         << "Invalid crossing angle, no optical functions produced.";
 
     currentDataValid_ = false;
     currentCrossingAngle_ = -1;
     currentData_ = std::make_shared<LHCInterpolatedOpticalFunctionsSetCollection>();
 
     return currentData_;
   }
 
   // set new crossing angle
   currentCrossingAngle_ = lhcInfo.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_;
 }