#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

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

private

Definition at line 55 of file HFPreReconstructor.cc.

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

private

Definition at line 56 of file HFPreReconstructor.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 85 of file HFPreReconstructor.cc.

References inputLabel_, and tok_hfQIE10_.

     : 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>();
 }

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

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

overrideprivate

Definition at line 199 of file HFPreReconstructor.cc.

References edm::EventSetup::get(), AlCaHLTBitMon_ParallelJobs::p, paramTS_, edm::ESHandle< T >::product(), and tsFromDB_.

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

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

static

Definition at line 268 of file HFPreReconstructor.cc.

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

                                                                                     {
   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);
 }

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

private

Definition at line 135 of file HFPreReconstructor.cc.

References edm::DataFrameContainer::begin(), AlignmentProducer_cff::calibrations, dropZSmarkedPassed_, edm::DataFrameContainer::end(), HcalRecoParam::firstSample(), forceSOI_, amptDefault_cfi::frame, edm::EventSetup::get(), edm::Event::getByToken(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalForward, QIE10DataFrame::id(), HFQIE10Info::id(), info(), paramTS_, QIE10DataFrame::presamples(), edm::ESHandle< T >::product(), qie10Infos_, DetId::rawId(), reco_, HFPreRecAlgo::reconstruct(), edm::DataFrameContainer::size(), soiShift_, tok_hfQIE10_, tsFromDB_, and QIE10DataFrame::zsMarkAndPass().

Referenced by produce().

                                                                                      {
   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);
     }
   }
 }

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

overrideprivate

Definition at line 208 of file HFPreReconstructor.cc.

References fillInfos(), dqmdumpme::first, mps_fire::i, eostools::move(), MillePedeFileConverter_cfg::out, edm::Event::put(), edm::second(), sortDataByPmt(), and sortedQIE10Infos_.

                                                                              {
   // 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));
 }

unsigned HFPreReconstructor::sortDataByPmt ( )

private

Definition at line 107 of file HFPreReconstructor.cc.

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

Referenced by produce().

                                            {
   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;
 }