SiStripGainFromAsciiFile Class Reference

#include <SiStripGainFromAsciiFile.h>

Inheritance diagram for SiStripGainFromAsciiFile:

Classes
struct	ModuleGain

Public Member Functions
	SiStripGainFromAsciiFile (const edm::ParameterSet &)
	~SiStripGainFromAsciiFile () override
Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)
	~ConditionDBWriter () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () 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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::unique_ptr< SiStripApvGain >	getNewObject () override

Private Attributes
std::string	Asciifilename_
edm::FileInPath	fp_
std::unordered_map< unsigned int, ModuleGain >	GainsMap
float	referenceValue_

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 ConditionDBWriter< SiStripApvGain >
static void	fillPSetDescription (edm::ParameterSetDescription &desc)
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 ConditionDBWriter< SiStripApvGain >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...
void	storeOnDbNow ()
cond::Time_t	timeOfLastIOV ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 13 of file SiStripGainFromAsciiFile.h.

Constructor & Destructor Documentation

SiStripGainFromAsciiFile()

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

explicit

Definition at line 13 of file SiStripGainFromAsciiFile.cc.

References Asciifilename_, fp_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SiStripDetInfoFileReader::kDefaultFile, referenceValue_, and AlCaHLTBitMon_QueryRunRegistry::string.

    : ConditionDBWriter<SiStripApvGain>(iConfig) {
  Asciifilename_ = iConfig.getParameter<std::string>("InputFileName");
  referenceValue_ = iConfig.getParameter<double>("referenceValue");
  fp_ = iConfig.getUntrackedParameter<edm::FileInPath>("file", edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile));
}

~SiStripGainFromAsciiFile()

SiStripGainFromAsciiFile::~SiStripGainFromAsciiFile ( )

override

Definition at line 20 of file SiStripGainFromAsciiFile.cc.

                                                    {
  edm::LogInfo("SiStripGainFromAsciiFile::~SiStripGainFromAsciiFile");
}

Member Function Documentation

getNewObject()

std::unique_ptr< SiStripApvGain > SiStripGainFromAsciiFile::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 24 of file SiStripGainFromAsciiFile.cc.

References SiStripGainFromAsciiFile::ModuleGain::apv, Asciifilename_, cms::cuda::assert(), DetId::det(), fp_, edm::FileInPath::fullPath(), GainsMap, SiStripGainFromAsciiFile::ModuleGain::hard_reset(), timingPdfMaker::infile, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, mps_splice::line, getGTfromDQMFile::obj, FastTimerService_cff::range, SiStripDetInfoFileReader::read(), referenceValue_, contentValuesCheck::ss, and DetId::Tracker.

                                                                     {
  auto obj = std::make_unique<SiStripApvGain>();

  std::stringstream ss;
  ss.str("");
  ss << "[SiStripGainFromAsciiFile::getNewObject]\n Reading Ascii File\n";
  FILE* infile = fopen(Asciifilename_.c_str(), "r");
  char line[4096];
  if (infile) {
    while (fgets(line, 4096, infile) != nullptr) {
      uint32_t detid;
      ModuleGain MG;
      MG.apv[0] = 0.0;
      MG.apv[1] = 0.0;
      MG.apv[2] = 0.0;
      MG.apv[3] = 0.0;
      MG.apv[4] = 0.0;
      MG.apv[5] = 0.0;
      char* saveptr;
      char* pch = strtok_r(line, " ", &saveptr);
      int Arg = 0;
      while (pch != nullptr) {
        if (Arg == 0) {
          sscanf(pch, "%d", &detid);
        } else if (Arg <= 6) {
          sscanf(pch, "%f", &(MG.apv[Arg - 1]));
        } else {
          //nothing to do here
        }
        pch = strtok_r(nullptr, " ", &saveptr);
        Arg++;
      }
      ss << detid << " " << MG.apv[0] << " " << MG.apv[1] << " " << MG.apv[2] << " " << MG.apv[3] << " " << MG.apv[4]
         << " " << MG.apv[5] << std::endl;
      GainsMap.insert(std::pair<unsigned int, ModuleGain>(detid, MG));
    }
    fclose(infile);
    edm::LogInfo("SiStripGainFromAsciiFile") << ss.str();
  } else {
    edm::LogError("SiStripGainFromAsciiFile")
        << " [SiStripGainFromAsciiFile::getNewObject] Error opening file " << Asciifilename_ << std::endl;
    assert(0);
  }

  const auto detInfo = SiStripDetInfoFileReader::read(fp_.fullPath());

  ss.str("");
  ss << "[SiStripGainFromAsciiFile::getNewObject]\n Filling SiStripApvGain object";
  short nApvPair;
  for (const auto it : detInfo.getAllDetIds()) {
    ModuleGain MG;
    if (DetId(it).det() != DetId::Tracker)
      continue;

    nApvPair = detInfo.getNumberOfApvsAndStripLength(it).first / 2;

    ss << "Looking at detid " << it << " nApvPairs  " << nApvPair << std::endl;
    auto iter = GainsMap.find(it);
    if (iter != GainsMap.end()) {
      MG = iter->second;
      ss << " " << MG.apv[0] << " " << MG.apv[1] << " " << MG.apv[2] << " " << MG.apv[3] << " " << MG.apv[4] << " "
         << MG.apv[5] << std::endl;
    } else {
      ss << "Hard reset for detid " << it << std::endl;
      MG.hard_reset(referenceValue_);
    }

    std::vector<float> DetGainsVector;

    if (nApvPair == 2) {
      DetGainsVector.push_back(MG.apv[0] / referenceValue_);
      DetGainsVector.push_back(MG.apv[1] / referenceValue_);
      DetGainsVector.push_back(MG.apv[2] / referenceValue_);
      DetGainsVector.push_back(MG.apv[3] / referenceValue_);
    } else if (nApvPair == 3) {
      DetGainsVector.push_back(MG.apv[0] / referenceValue_);
      DetGainsVector.push_back(MG.apv[1] / referenceValue_);
      DetGainsVector.push_back(MG.apv[2] / referenceValue_);
      DetGainsVector.push_back(MG.apv[3] / referenceValue_);
      DetGainsVector.push_back(MG.apv[4] / referenceValue_);
      DetGainsVector.push_back(MG.apv[5] / referenceValue_);
    } else {
      edm::LogError("SiStripGainFromAsciiFile") << " SiStripGainFromAsciiFile::getNewObject] ERROR for detid " << it
                                                << " not expected number of APV pairs " << nApvPair << std::endl;
    }

    SiStripApvGain::Range range(DetGainsVector.begin(), DetGainsVector.end());
    if (!obj->put(it, range)) {
      edm::LogError("SiStripGainFromAsciiFile")
          << " [SiStripGainFromAsciiFile::getNewObject] detid already exists" << std::endl;
      ss << " [SiStripGainFromAsciiFile::getNewObject] detid already exists" << std::endl;
    }
  }
  edm::LogInfo("SiStripGainFromAsciiFile") << ss.str();

  return obj;
}

Member Data Documentation

Asciifilename_

std::string SiStripGainFromAsciiFile::Asciifilename_

private

Definition at line 36 of file SiStripGainFromAsciiFile.h.

Referenced by getNewObject(), and SiStripGainFromAsciiFile().

fp_

edm::FileInPath SiStripGainFromAsciiFile::fp_

private

Definition at line 38 of file SiStripGainFromAsciiFile.h.

Referenced by getNewObject(), and SiStripGainFromAsciiFile().

GainsMap

std::unordered_map<unsigned int, ModuleGain> SiStripGainFromAsciiFile::GainsMap

private

Definition at line 40 of file SiStripGainFromAsciiFile.h.

Referenced by getNewObject().

referenceValue_

float SiStripGainFromAsciiFile::referenceValue_

private

Definition at line 37 of file SiStripGainFromAsciiFile.h.

Referenced by getNewObject(), and SiStripGainFromAsciiFile().