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 noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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(), ALCARECOPPSCalTrackBasedSel_cff::detid, fp_, edm::FileInPath::fullPath(), GainsMap, SiStripGainFromAsciiFile::ModuleGain::hard_reset(), timingPdfMaker::infile, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, mps_splice::line, getGTfromDQMFile::obj, isotrackApplyRegressor::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().