HcalPatternSource Class Reference

#include <HcalPatternSource.h>

Inheritance diagram for HcalPatternSource:

Public Member Functions
	HcalPatternSource (const edm::ParameterSet &pset)
void	produce (edm::Event &e, const edm::EventSetup &c) override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) 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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	loadPatternFile (const std::string &filename)
void	loadPatterns (const std::string &patspec)

Private Attributes
std::vector< int >	bunches_
std::vector< HcalFiberPattern >	patterns_
int	presamples_
int	samples_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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 () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
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

Author

J. Mans - Minnesota

Definition at line 12 of file HcalPatternSource.h.

Constructor & Destructor Documentation

HcalPatternSource::HcalPatternSource

(

const edm::ParameterSet &

pset

)

Definition at line 14 of file HcalPatternSource.cc.

References edm::ParameterSet::getUntrackedParameter(), loadPatterns(), and AlCaHLTBitMon_QueryRunRegistry::string.

    : bunches_(pset.getUntrackedParameter<std::vector<int> >("Bunches", std::vector<int>())),
      presamples_(pset.getUntrackedParameter<int>("Presamples", 4)),
      samples_(pset.getUntrackedParameter<int>("Samples", 10)) {
  loadPatterns(pset.getUntrackedParameter<std::string>("Patterns"));
  produces<HBHEDigiCollection>();
  produces<HODigiCollection>();
  produces<HFDigiCollection>();
}

Member Function Documentation

void HcalPatternSource::loadPatternFile ( const std::string &  filename )

private

Definition at line 103 of file HcalPatternSource.cc.

References createPayload::block, edmScanValgrind::buffer, data, Exception, validate-o2o-wbm::f, lut2db_cfg::filename, mps_fire::i, dqmiolumiharvest::j, submitPVValidationJobs::params, HcalPatternXMLParser::parse(), writedatasetfile::parser, patterns_, SiPixelLorentzAngle_cfi::read, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by loadPatterns().

                                                                 {
  HcalPatternXMLParser parser;
  std::string buffer, element;
  std::map<std::string, std::string> params;
  std::vector<uint32_t> data;
  FILE* f = fopen(filename.c_str(), "r");
  if (f == nullptr)
    return;
  else {
    char block[4096];
    while (!feof(f)) {
      int read = fread(block, 1, 4096, f);
      buffer.append(block, block + read);
    }
    fclose(f);
  }
  if (buffer.find("<?xml") != 0) {
    throw cms::Exception("InvalidFormat") << "Not a valid XML file: " << filename;
  }
  std::string::size_type i = 0, j;
  while (buffer.find("<CFGBrick>", i) != std::string::npos) {
    i = buffer.find("<CFGBrick>", i);
    j = buffer.find("</CFGBrick>", i);
    element = "<?xml version='1.0'?>\n";
    element.append(buffer, i, j - i);
    element.append("</CFGBrick>");
    //    LogDebug("HCAL") << element;
    params.clear();
    data.clear();
    parser.parse(element, params, data);
    patterns_.push_back(HcalFiberPattern(params, data));
    i = j + 5;
  }
}

void HcalPatternSource::loadPatterns ( const std::string &  patspec )

private

Definition at line 90 of file HcalPatternSource.cc.

References submitPVResolutionJobs::files, mps_fire::i, loadPatternFile(), LogDebug, AlCaHLTBitMon_ParallelJobs::p, and patterns_.

Referenced by HcalPatternSource().

                                                             {
  wordexp_t p;
  char** files;
  wordexp(patspec.c_str(), &p, WRDE_NOCMD);  // do not run shell commands!
  files = p.we_wordv;
  for (unsigned int i = 0; i < p.we_wordc; i++) {
    LogDebug("HCAL") << "Reading pattern file '" << files[i] << "'";
    loadPatternFile(files[i]);
    LogDebug("HCAL") << "Fibers so far " << patterns_.size();
  }
  wordfree(&p);
}

void HcalPatternSource::produce ( edm::Event &  e, const edm::EventSetup &  c  )

overridevirtual

Implements edm::EDProducer.

Definition at line 24 of file HcalPatternSource.cc.

References bunches_, runTauDisplay::eid, edm::EventID::event(), edm::EventSetup::get(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcalSimParameters_cfi::ho, mps_fire::i, edm::EventBase::id(), B2GTnPMonitor_cfi::item, HcalElectronicsMap::lookup(), eostools::move(), patterns_, presamples_, edm::ESHandle< class >::product(), edm::Event::put(), and samples_.

                                                                    {
  if (e.id().event() > bunches_.size())
    return;

  edm::ESHandle<HcalElectronicsMap> item;
  es.get<HcalElectronicsMapRcd>().get(item);
  const HcalElectronicsMap* elecmap = item.product();

  auto hbhe = std::make_unique<HBHEDigiCollection>();
  auto hf = std::make_unique<HFDigiCollection>();
  auto ho = std::make_unique<HODigiCollection>();

  int bc = bunches_[e.id().event() - 1];
  for (std::vector<HcalFiberPattern>::iterator i = patterns_.begin(); i != patterns_.end(); i++) {
    std::vector<HcalQIESample> samples;
    for (int fc = 0; fc < 3; fc++) {
      samples = i->getSamples(bc, presamples_, samples_, fc);
      HcalElectronicsId eid = i->getId(fc);

      HcalDetId did(elecmap->lookup(eid));

      if (did.null()) {
        edm::LogWarning("HCAL") << "No electronics map match for id " << eid;
        continue;
      }

      switch (did.subdet()) {
        case (HcalBarrel):
        case (HcalEndcap):
          hbhe->push_back(HBHEDataFrame(did));
          hbhe->back().setSize(samples_);
          hbhe->back().setPresamples(presamples_);
          for (int i = 0; i < samples_; i++)
            hbhe->back().setSample(i, samples[i]);
          hbhe->back().setReadoutIds(eid);
          break;
        case (HcalForward):
          hf->push_back(HFDataFrame(did));
          hf->back().setSize(samples_);
          hf->back().setPresamples(presamples_);
          for (int i = 0; i < samples_; i++)
            hf->back().setSample(i, samples[i]);
          hf->back().setReadoutIds(eid);
          break;
        case (HcalOuter):
          ho->push_back(HODataFrame(did));
          ho->back().setSize(samples_);
          ho->back().setPresamples(presamples_);
          for (int i = 0; i < samples_; i++)
            ho->back().setSample(i, samples[i]);
          ho->back().setReadoutIds(eid);
          break;
        default:
          continue;
      }
    }
  }
  hbhe->sort();
  ho->sort();
  hf->sort();

  e.put(std::move(hbhe));
  e.put(std::move(ho));
  e.put(std::move(hf));
}

Member Data Documentation

std::vector<int> HcalPatternSource::bunches_

private

Definition at line 20 of file HcalPatternSource.h.

Referenced by produce().

std::vector<HcalFiberPattern> HcalPatternSource::patterns_

private

Definition at line 21 of file HcalPatternSource.h.

Referenced by loadPatternFile(), loadPatterns(), and produce().

int HcalPatternSource::presamples_

private

Definition at line 22 of file HcalPatternSource.h.

Referenced by produce().

int HcalPatternSource::samples_

private

Definition at line 22 of file HcalPatternSource.h.

Referenced by produce().