#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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::ESGetToken< HcalTopology, HcalRecNumberingRecord >	htopoToken_

edm::InputTag	inputLabel_

edm::ESGetToken < HcalChannelPropertiesVec, HcalChannelPropertiesRecord >	propertiesToken_

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

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

Implementation: [Notes on implementation]

Definition at line 52 of file HFPreReconstructor.cc.

Member Typedef Documentation

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

private

Definition at line 60 of file HFPreReconstructor.cc.

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

private

Definition at line 61 of file HFPreReconstructor.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 93 of file HFPreReconstructor.cc.

References htopoToken_, inputLabel_, propertiesToken_, 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>();
 
   // ES tokens
   htopoToken_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   propertiesToken_ = esConsumes<HcalChannelPropertiesVec, HcalChannelPropertiesRecord>();
 }

HFPreReconstructor::~HFPreReconstructor ( )

override

Definition at line 111 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 205 of file HFPreReconstructor.cc.

205 {}

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

static

Definition at line 268 of file HFPreReconstructor.cc.

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

                                                                                     {
   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 147 of file HFPreReconstructor.cc.

References HcalChannelProperties::channelCoder, dropZSmarkedPassed_, HcalRecoParam::firstSample(), forceSOI_, edm::Event::getByToken(), edm::EventSetup::getData(), HcalForward, htopoToken_, HFQIE10Info::id(), QIE10DataFrame::id(), info(), HcalChannelProperties::paramTs, QIE10DataFrame::presamples(), propertiesToken_, qie10Infos_, DetId::rawId(), reco_, HFPreRecAlgo::reconstruct(), HcalChannelProperties::shape, 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 calibrations
   const HcalTopology& htopo = eventSetup.getData(htopoToken_);
   const HcalChannelPropertiesVec& prop = eventSetup.getData(propertiesToken_);
 
   // 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;
 
       // Look up the channel properties. This lookup is O(1).
       const HcalChannelProperties& properties(prop.at(htopo.detId2denseId(cell)));
 
       // ADC decoding tool
       const HcalCoderDb coder(*properties.channelCoder, *properties.shape);
 
       int tsToUse = forceSOI_;
       if (tsToUse < 0) {
         if (tsFromDB_) {
           tsToUse = properties.paramTs->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, properties);
       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 cms::cuda::assert(), fillInfos(), first, mps_fire::i, eostools::move(), submitPVResolutionJobs::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 119 of file HFPreReconstructor.cc.

References HLT_FULL_cff::depth, first, mps_fire::i, gpuClustering::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;
 }