magneticfield::AutoParametrizedMagneticFieldProducer Class Reference

Inheritance diagram for magneticfield::AutoParametrizedMagneticFieldProducer:

Public Member Functions
	AutoParametrizedMagneticFieldProducer (const edm::ParameterSet &)
int	closerNominaCurrent (float current)
std::unique_ptr< MagneticField >	produce (const IdealMagneticFieldRecord &)
	~AutoParametrizedMagneticFieldProducer ()
Public Member Functions inherited from edm::ESProducer
	ESProducer ()
virtual	~ESProducer ()(false)
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
virtual void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval)
	overrides DataProxyProvider method More...
virtual	~ESProxyFactoryProducer () noexcept(false)
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()
const ComponentDescription &	description () const
bool	isUsingRecord (const EventSetupRecordKey &) const
const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const
void	resetProxies (const EventSetupRecordKey &iRecordType)
void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)
void	setAppendToDataLabel (const edm::ParameterSet &)
void	setDescription (const ComponentDescription &iDescription)
std::set< EventSetupRecordKey >	usingRecords () const
virtual	~DataProxyProvider () noexcept(false)

Public Attributes
std::vector< int >	nominalCurrents
edm::ParameterSet	pset

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair< DataKey, edm::propagate_const< std::shared_ptr< DataProxy > > > >	KeyedProxies
typedef std::vector< EventSetupRecordKey >	Keys
typedef std::map< EventSetupRecordKey, KeyedProxies >	RecordProxies
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::ESProducer
template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())
template<typename T >
void	setWhatProduced (T *iThis, const char *iLabel)
template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)
template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())
template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel=es::Label())
template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())
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 eventsetup::EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())
virtual void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
	override DataProxyProvider method More...
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...
void	invalidateProxies (const EventSetupRecordKey &iRecordKey)
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &)

Detailed Description

Definition at line 29 of file AutoParametrizedMagneticFieldProducer.cc.

Constructor & Destructor Documentation

AutoParametrizedMagneticFieldProducer::AutoParametrizedMagneticFieldProducer

(

const edm::ParameterSet &

iConfig

)

Definition at line 46 of file AutoParametrizedMagneticFieldProducer.cc.

References edm::ParameterSet::getUntrackedParameter(), nominalCurrents, pset, edm::ESProducer::setWhatProduced(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                           : pset(iConfig) {
  setWhatProduced(this, pset.getUntrackedParameter<std::string>("label",""));
  nominalCurrents={-1, 0,9558,14416,16819,18268,19262};
  //  nominalLabels  ={"3.8T","0T","2T", "3T", "3.5T", "3.8T", "4T"};
}

magneticfield::AutoParametrizedMagneticFieldProducer::~AutoParametrizedMagneticFieldProducer ( )

inline

Definition at line 33 of file AutoParametrizedMagneticFieldProducer.cc.

{}

Member Function Documentation

int AutoParametrizedMagneticFieldProducer::closerNominaCurrent

(

float

current

)

Definition at line 98 of file AutoParametrizedMagneticFieldProducer.cc.

References DEFINE_FWK_EVENTSETUP_MODULE, mps_fire::i, createfilelist::int, and nominalCurrents.

Referenced by produce().

                                                                            {
  int i=0;
  for(;i<(int)nominalCurrents.size()-1;i++) {
    if(2*current < nominalCurrents[i]+nominalCurrents[i+1] )
      return nominalCurrents[i];
  }
  return nominalCurrents[i];
}

std::unique_ptr< MagneticField > AutoParametrizedMagneticFieldProducer::produce

(

const IdealMagneticFieldRecord &

iRecord

)

Definition at line 53 of file AutoParametrizedMagneticFieldProducer.cc.

References closerNominaCurrent(), MillePedeFileConverter_cfg::e, Exception, ParametrizedMagneticFieldFactory::get(), edm::ParameterSet::getParameter(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), RunInfo::m_avg_current, python.rootplot.argparse::message, metProducer_cfi::parameters, pset, Scenarios_cff::scale, and jetIDSelector_cfi::version.

{
  string version = pset.getParameter<string>("version");

  // Get value of the current from condition DB
  float current = pset.getParameter<int>("valueOverride");
  string message;
  if (current < 0) {
    ESHandle<RunInfo> rInfo;
    iRecord.getRecord<RunInfoRcd>().get(rInfo);
    current = rInfo->m_avg_current;
    message = " (from RunInfo DB)";
  } else {
    message = " (from valueOverride card)";
  }
  float cnc= closerNominaCurrent(current);

  edm::LogInfo("MagneticField|AutoParametrizedMagneticField") << "Current: " << current << message << "; using map for: " << cnc;

  vector<double> parameters;
  
  if (cnc==0) {
    version = "Uniform";
    parameters.push_back(0);
  }

  else if (version=="Parabolic"){
    parameters.push_back(3.8114);       //c1
    parameters.push_back(-3.94991e-06); //b0
    parameters.push_back(7.53701e-06);  //b1
    parameters.push_back(2.43878e-11);   //a
    if (cnc!=18268){ // Linear scaling for B!= 3.8T; note that just c1, b0 and b1 have to be scaled to get linear scaling
      double scale=double(cnc)/double(18268);
      parameters[0]*=scale;
      parameters[1]*=scale;
      parameters[2]*=scale;
    }
  } else {  
    //Other parametrizations are not relevant here and not supported
    throw cms::Exception("InvalidParameter") << "version " << version << " is not supported";
  }

  return ParametrizedMagneticFieldFactory::get(version, parameters);
}

Member Data Documentation

std::vector<int> magneticfield::AutoParametrizedMagneticFieldProducer::nominalCurrents

Definition at line 39 of file AutoParametrizedMagneticFieldProducer.cc.

Referenced by AutoParametrizedMagneticFieldProducer(), and closerNominaCurrent().

edm::ParameterSet magneticfield::AutoParametrizedMagneticFieldProducer::pset

Definition at line 38 of file AutoParametrizedMagneticFieldProducer.cc.

Referenced by AutoParametrizedMagneticFieldProducer(), and produce().