modGains Class Reference

Inheritance diagram for modGains:

Public Member Functions
	modGains (const edm::ParameterSet &)
	~modGains () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	beginRun (edm::Run const &iEvent, edm::EventSetup const &) override
void	endRun (edm::Run const &iEvent, edm::EventSetup const &) override

Private Attributes
std::string	fileCorr
std::string	fileIn
std::string	fileOut
std::string	s_operation
double	val
bool	vectorop

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 21 of file modGains.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file modGains.cc.

References MessageLogger_cfi::cerr, Exception, fileCorr, fileIn, fileOut, edm::ParameterSet::getUntrackedParameter(), s_operation, AlCaHLTBitMon_QueryRunRegistry::string, val, and vectorop.

                                                 : vectorop(false) {
  s_operation = iConfig.getUntrackedParameter<std::string>("Operation");
  fileIn      = iConfig.getUntrackedParameter<std::string>("FileIn");
  fileOut     = iConfig.getUntrackedParameter<std::string>("FileOut");
  fileCorr    = iConfig.getUntrackedParameter<std::string>("FileCorr");
  val         = iConfig.getUntrackedParameter<double>("ScalarFactor");

  if ( (std::strcmp(s_operation.c_str(),"add")==0) || 
       (std::strcmp(s_operation.c_str(),"sub")==0) || 
       (std::strcmp(s_operation.c_str(),"mult")==0) || 
       (std::strcmp(s_operation.c_str(),"div")==0) ) { // vector operation
    vectorop = true;
  }  else if ((std::strcmp(s_operation.c_str(),"sadd")==0) || 
          (std::strcmp(s_operation.c_str(),"ssub")==0) || 
          (std::strcmp(s_operation.c_str(),"smult")==0) || 
          (std::strcmp(s_operation.c_str(),"sdiv")==0)) {// scalar operation
    std::cerr << "Scalar operation: using val=" << val << std::endl;
  } else {
    throw cms::Exception("Unknown", "modGains") << "Unknown operator. "
                        << s_operation <<" Stopping.\n";
  }
}

modGains::~modGains ( )

override

Definition at line 60 of file modGains.cc.

{ }

Member Function Documentation

void modGains::analyze ( edm::Event const &  , edm::EventSetup const &  iSetup  )

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 62 of file modGains.cc.

References HcalCondObjectContainer< Item >::addValues(), MessageLogger_cfi::cerr, DEFINE_FWK_MODULE, HcalDbASCIIIO::dumpObject(), HcalCondObjectContainer< Item >::exists(), fileCorr, fileIn, fileOut, edm::EventSetup::get(), HcalDbASCIIIO::getObject(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), mps_fire::i, s_operation, val, and vectorop.

Referenced by beginRun().

                                                                   {

  edm::ESHandle<HcalTopology> htopo;
  iSetup.get<HcalRecNumberingRecord>().get(htopo);
  HcalTopology topo = (*htopo);

  // get base conditions
  std::cerr << fileIn << std::endl;
  std::ifstream inStream  (fileIn.c_str());
  HcalGains gainsIn(&topo);;
  HcalDbASCIIIO::getObject (inStream, &gainsIn);
  inStream.close();

  HcalRespCorrs corrsIn(&topo);;
  if (vectorop) {
    std::ifstream inCorr     (fileCorr.c_str());
    HcalDbASCIIIO::getObject (inCorr, &corrsIn);
    inCorr.close();
  }

  HcalGains gainsOut(&topo);;
  std::vector<DetId> channels = gainsIn.getAllChannels ();
  std::cerr << "size = " << channels.size() << std::endl;
  for (unsigned i = 0; i < channels.size(); i++) {
    DetId id = channels[i];

    if (vectorop) { // vector operation
      if ((std::strcmp(s_operation.c_str(),"mult")==0)||(std::strcmp(s_operation.c_str(),"div")==0)) val = 1.0; // mult,div
      if ((std::strcmp(s_operation.c_str(),"add")==0)||(std::strcmp(s_operation.c_str(),"sub")==0)) val = 0.0; // add,sub
      if (corrsIn.exists(id))  {
    val = corrsIn.getValues(id)->getValue();
      }
      if (i%100 == 0)
    std::cerr << "Vector operation, " << i << "th channel: using val=" << val << std::endl;
    }
    
    HcalGain* p_item = nullptr;
    if ((std::strcmp(s_operation.c_str(),"add")==0) || (std::strcmp(s_operation.c_str(),"sadd")==0))
      p_item = new HcalGain(id, gainsIn.getValues(id)->getValue(0) + val, gainsIn.getValues(id)->getValue(1) + val, 
                gainsIn.getValues(id)->getValue(2) + val, gainsIn.getValues(id)->getValue(3) + val);

    if ((std::strcmp(s_operation.c_str(),"sub")==0) || (std::strcmp(s_operation.c_str(),"ssub")==0))
      p_item = new HcalGain(id, gainsIn.getValues(id)->getValue(0) - val, gainsIn.getValues(id)->getValue(1) - val, 
                gainsIn.getValues(id)->getValue(2) - val, gainsIn.getValues(id)->getValue(3) - val);
    
    if ((std::strcmp(s_operation.c_str(),"mult")==0) || (std::strcmp(s_operation.c_str(),"smult")==0))
      p_item = new HcalGain(id, gainsIn.getValues(id)->getValue(0) * val, gainsIn.getValues(id)->getValue(1) * val, 
                gainsIn.getValues(id)->getValue(2) * val, gainsIn.getValues(id)->getValue(3) * val);

    if ((std::strcmp(s_operation.c_str(),"div")==0) || (std::strcmp(s_operation.c_str(),"sdiv")==0))
      p_item = new HcalGain(id, gainsIn.getValues(id)->getValue(0) / val, gainsIn.getValues(id)->getValue(1) / val, 
                gainsIn.getValues(id)->getValue(2) / val, gainsIn.getValues(id)->getValue(3) / val);


    // for all
    if (p_item)  gainsOut.addValues(*p_item);
    //    std::cerr << i << std::endl;
  }
  // write out
  std::ofstream outStream (fileOut.c_str());
  HcalDbASCIIIO::dumpObject (outStream, gainsOut);
  outStream.close();
}

void modGains::beginRun ( edm::Run const &  iEvent, edm::EventSetup const &    )

inlineoverrideprivate

Definition at line 28 of file modGains.cc.

References analyze().

{}

void modGains::endRun ( edm::Run const &  iEvent, edm::EventSetup const &    )

inlineoverrideprivate

Definition at line 30 of file modGains.cc.

{}

Member Data Documentation

std::string modGains::fileCorr

private

Definition at line 32 of file modGains.cc.

Referenced by analyze(), and modGains().

std::string modGains::fileIn

private

Definition at line 32 of file modGains.cc.

Referenced by analyze(), and modGains().

std::string modGains::fileOut

private

Definition at line 32 of file modGains.cc.

Referenced by analyze(), and modGains().

std::string modGains::s_operation

private

Definition at line 31 of file modGains.cc.

Referenced by analyze(), and modGains().

double modGains::val

private

Definition at line 33 of file modGains.cc.

Referenced by value.Value::__eq__(), value.Value::__iadd__(), value.Value::__idiv__(), value.Value::__isub__(), value.Value::__str__(), analyze(), modGains(), and value.Value::relerr().

bool modGains::vectorop

private

Definition at line 34 of file modGains.cc.

Referenced by analyze(), and modGains().