HFPreReconstructor Class Reference

#include <RecoLocalCalo/HcalRecProducers/src/HFPreReconstructor.cc>

Inheritance diagram for HFPreReconstructor:

Public Member Functions
	HFPreReconstructor (const edm::ParameterSet &)
	~HFPreReconstructor () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Types
typedef std::pair< HcalDetId, int >	PmtAnodeId
typedef std::pair< PmtAnodeId, const HFQIE10Info * >	QIE10InfoWithId

Private Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override
void	fillInfos (const edm::Event &e, const edm::EventSetup &eventSetup)
void	produce (edm::Event &, const edm::EventSetup &) override
unsigned	sortDataByPmt ()

Private Attributes
bool	dropZSmarkedPassed_
int	forceSOI_
edm::InputTag	inputLabel_
std::unique_ptr< HcalRecoParams >	paramTS_
std::vector< HFQIE10Info >	qie10Infos_
HFPreRecAlgo	reco_
int	soiShift_
std::vector< QIE10InfoWithId >	sortedQIE10Infos_
edm::EDGetTokenT< QIE10DigiCollection >	tok_hfQIE10_
bool	tsFromDB_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: Phase 1 HF reco with QIE 10 and split-anode readout

Implementation: [Notes on implementation]

Definition at line 47 of file HFPreReconstructor.cc.

Member Typedef Documentation

PmtAnodeId

typedef std::pair<HcalDetId, int> HFPreReconstructor::PmtAnodeId

private

Definition at line 55 of file HFPreReconstructor.cc.

QIE10InfoWithId

typedef std::pair<PmtAnodeId, const HFQIE10Info*> HFPreReconstructor::QIE10InfoWithId

private

Definition at line 56 of file HFPreReconstructor.cc.

Constructor & Destructor Documentation

HFPreReconstructor()

HFPreReconstructor::HFPreReconstructor ( const edm::ParameterSet &  conf )

explicit

Definition at line 85 of file HFPreReconstructor.cc.

    : inputLabel_(conf.getParameter<edm::InputTag>("digiLabel")),
      forceSOI_(conf.getParameter<int>("forceSOI")),
      soiShift_(conf.getParameter<int>("soiShift")),
      dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")),
      tsFromDB_(conf.getParameter<bool>("tsFromDB")),
      reco_(conf.getParameter<bool>("sumAllTimeSlices")) {
  // Describe consumed data
  tok_hfQIE10_ = consumes<QIE10DigiCollection>(inputLabel_);
 
  // Register the product
  produces<HFPreRecHitCollection>();
}

References inputLabel_, and tok_hfQIE10_.

~HFPreReconstructor()

HFPreReconstructor::~HFPreReconstructor ( )

override

Definition at line 99 of file HFPreReconstructor.cc.

                                        {
  // do anything here that needs to be done at destruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

beginRun()

void HFPreReconstructor::beginRun ( const edm::Run &  r, const edm::EventSetup &  es  )

overrideprivate

Definition at line 199 of file HFPreReconstructor.cc.

                                                                          {
  if (tsFromDB_) {
    edm::ESHandle<HcalRecoParams> p;
    es.get<HcalRecoParamsRcd>().get(p);
    paramTS_ = std::make_unique<HcalRecoParams>(*p.product());
  }
}

References edm::EventSetup::get(), get, AlCaHLTBitMon_ParallelJobs::p, paramTS_, and tsFromDB_.

fillDescriptions()

void HFPreReconstructor::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 268 of file HFPreReconstructor.cc.

                                                                                    {
  edm::ParameterSetDescription desc;
 
  desc.add<edm::InputTag>("digiLabel");
  desc.add<int>("forceSOI", -1);
  desc.add<int>("soiShift", 0);
  desc.add<bool>("dropZSmarkedPassed");
  desc.add<bool>("tsFromDB");
  desc.add<bool>("sumAllTimeSlices");
 
  descriptions.addDefault(desc);
}

References edm::ParameterSetDescription::add(), and edm::ConfigurationDescriptions::addDefault().

fillInfos()

void HFPreReconstructor::fillInfos ( const edm::Event &  e, const edm::EventSetup &  eventSetup  )

private

Definition at line 135 of file HFPreReconstructor.cc.

                                                                                     {
  using namespace edm;
 
  // Clear the collection we want to fill in this method
  qie10Infos_.clear();
 
  // Get the Hcal topology if needed
  ESHandle<HcalTopology> htopo;
  if (tsFromDB_) {
    eventSetup.get<HcalRecNumberingRecord>().get(htopo);
    paramTS_->setTopo(htopo.product());
  }
 
  // Get the calibrations
  ESHandle<HcalDbService> conditions;
  eventSetup.get<HcalDbRecord>().get(conditions);
 
  // Get the input collection
  Handle<QIE10DigiCollection> digi;
  e.getByToken(tok_hfQIE10_, digi);
 
  const unsigned inputSize = digi->size();
  if (inputSize) {
    // Process the digis and fill out the HFQIE10Info vector
    qie10Infos_.reserve(inputSize);
 
    for (QIE10DigiCollection::const_iterator it = digi->begin(); it != digi->end(); ++it) {
      const QIE10DataFrame& frame(*it);
      const HcalDetId cell(frame.id());
 
      // Protection against calibration channels which are not
      // in the database but can still come in the QIE10DataFrame
      // in the laser calibs, etc.
      if (cell.subdet() != HcalSubdetector::HcalForward)
        continue;
 
      // Check zero suppression
      if (dropZSmarkedPassed_)
        if (frame.zsMarkAndPass())
          continue;
 
      const HcalCalibrations& calibrations(conditions->getHcalCalibrations(cell));
      const HcalQIECoder* channelCoder = conditions->getHcalCoder(cell);
      const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
      const HcalCoderDb coder(*channelCoder, *shape);
 
      int tsToUse = forceSOI_;
      if (tsToUse < 0) {
        if (tsFromDB_) {
          const HcalRecoParam* param_ts = paramTS_->getValues(cell.rawId());
          tsToUse = param_ts->firstSample();
        } else
          // Get the "sample of interest" from the data frame itself
          tsToUse = frame.presamples();
      }
 
      // Reconstruct the charge, energy, etc
      const HFQIE10Info& info = reco_.reconstruct(frame, tsToUse + soiShift_, coder, calibrations);
      if (info.id().rawId())
        qie10Infos_.push_back(info);
    }
  }
}

References edm::DataFrameContainer::begin(), AlignmentProducer_cff::calibrations, dropZSmarkedPassed_, MillePedeFileConverter_cfg::e, edm::DataFrameContainer::end(), options_cfi::eventSetup, HcalRecoParam::firstSample(), forceSOI_, amptDefault_cfi::frame, get, HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalForward, info(), paramTS_, edm::ESHandle< T >::product(), qie10Infos_, reco_, HFPreRecAlgo::reconstruct(), edm::DataFrameContainer::size(), soiShift_, tok_hfQIE10_, and tsFromDB_.

Referenced by produce().

produce()

void HFPreReconstructor::produce ( edm::Event &  e, const edm::EventSetup &  eventSetup  )

overrideprivate

Definition at line 208 of file HFPreReconstructor.cc.

                                                                             {
  // Process the input data
  fillInfos(e, eventSetup);
 
  // Create a new output collection
  std::unique_ptr<HFPreRecHitCollection> out(std::make_unique<HFPreRecHitCollection>());
 
  // Fill the output collection
  const unsigned pmtCount = sortDataByPmt();
  if (pmtCount) {
    out->reserve(pmtCount);
    const unsigned sz = sortedQIE10Infos_.size();
    HcalDetId previousBaseId(sortedQIE10Infos_[0].first.first);
    unsigned nFound = 1;
    for (unsigned i = 1; i <= sz; ++i) {
      bool appendData = i == sz;
      if (i < sz) {
        const HcalDetId baseId(sortedQIE10Infos_[i].first.first);
        if (baseId == previousBaseId)
          ++nFound;
        else {
          appendData = true;
          previousBaseId = baseId;
        }
      }
 
      if (appendData) {
        // If we have found more than two QIE10 with the same base id,
        // there is a bug somewhere in the dataframe. We can't do
        // anything useful about it here. Once we make sure that
        // everything works as expected, this assertion can be removed.
        assert(nFound <= 2);
 
        const HFQIE10Info* first = nullptr;
        const HFQIE10Info* second = sortedQIE10Infos_[i - 1].second;
 
        if (nFound >= 2)
          first = sortedQIE10Infos_[i - 2].second;
        else if (sortedQIE10Infos_[i - 1].first.second < 3) {
          // Only one QIE10 readout found for this PMT.
          // Arrange for depth 1 and 2 to be "first".
          first = second;
          second = nullptr;
        }
 
        out->push_back(HFPreRecHit(sortedQIE10Infos_[i - nFound].first.first, first, second));
 
        // Reset the QIE find count for this base id
        nFound = 1;
      }
    }
 
    assert(out->size() == pmtCount);
  }
 
  // Add the output collection to the event record
  e.put(std::move(out));
}

References cms::cuda::assert(), MillePedeFileConverter_cfg::e, options_cfi::eventSetup, fillInfos(), dqmdumpme::first, mps_fire::i, eostools::move(), MillePedeFileConverter_cfg::out, edm::second(), sortDataByPmt(), and sortedQIE10Infos_.

sortDataByPmt()

unsigned HFPreReconstructor::sortDataByPmt ( )

private

Definition at line 107 of file HFPreReconstructor.cc.

                                           {
  sortedQIE10Infos_.clear();
  unsigned pmtCount = 0;
  const unsigned sz = qie10Infos_.size();
  if (sz) {
    // Perform sorting
    sortedQIE10Infos_.reserve(sz);
    const HFQIE10Info* info = &qie10Infos_[0];
    for (unsigned i = 0; i < sz; ++i) {
      const HcalDetId id(info[i].id());
      sortedQIE10Infos_.push_back(QIE10InfoWithId(PmtAnodeId(id.baseDetId(), id.depth()), info + i));
    }
    std::sort(sortedQIE10Infos_.begin(), sortedQIE10Infos_.end());
 
    // Count the PMTs
    HcalDetId previousBaseId(sortedQIE10Infos_[0].first.first);
    pmtCount = 1;
    for (unsigned i = 1; i < sz; ++i) {
      const HcalDetId baseId(sortedQIE10Infos_[i].first.first);
      if (baseId != previousBaseId) {
        previousBaseId = baseId;
        ++pmtCount;
      }
    }
  }
  return pmtCount;
}

References LEDCalibrationChannels::depth, dqmdumpme::first, mps_fire::i, triggerObjects_cff::id, info(), qie10Infos_, and sortedQIE10Infos_.

Referenced by produce().

Member Data Documentation

dropZSmarkedPassed_

bool HFPreReconstructor::dropZSmarkedPassed_

private

Definition at line 65 of file HFPreReconstructor.cc.

Referenced by fillInfos().

forceSOI_

int HFPreReconstructor::forceSOI_

private

Definition at line 63 of file HFPreReconstructor.cc.

Referenced by fillInfos().

inputLabel_

edm::InputTag HFPreReconstructor::inputLabel_

private

Definition at line 62 of file HFPreReconstructor.cc.

Referenced by HFPreReconstructor().

paramTS_

std::unique_ptr<HcalRecoParams> HFPreReconstructor::paramTS_

private

Definition at line 73 of file HFPreReconstructor.cc.

Referenced by beginRun(), and fillInfos().

qie10Infos_

std::vector<HFQIE10Info> HFPreReconstructor::qie10Infos_

private

Definition at line 71 of file HFPreReconstructor.cc.

Referenced by fillInfos(), and sortDataByPmt().

reco_

HFPreRecAlgo HFPreReconstructor::reco_

private

Definition at line 69 of file HFPreReconstructor.cc.

Referenced by fillInfos().

soiShift_

int HFPreReconstructor::soiShift_

private

Definition at line 64 of file HFPreReconstructor.cc.

Referenced by fillInfos().

sortedQIE10Infos_

std::vector<QIE10InfoWithId> HFPreReconstructor::sortedQIE10Infos_

private

Definition at line 72 of file HFPreReconstructor.cc.

Referenced by produce(), and sortDataByPmt().

tok_hfQIE10_

edm::EDGetTokenT<QIE10DigiCollection> HFPreReconstructor::tok_hfQIE10_

private

Definition at line 70 of file HFPreReconstructor.cc.

Referenced by fillInfos(), and HFPreReconstructor().

tsFromDB_

bool HFPreReconstructor::tsFromDB_

private

Definition at line 66 of file HFPreReconstructor.cc.

Referenced by beginRun(), and fillInfos().