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

Inheritance diagram for HFPreReconstructor:

Public Member Functions
	HFPreReconstructor (const edm::ParameterSet &)

	~HFPreReconstructor ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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 ()

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, 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	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Types
typedef std::pair< HcalDetId, int >	PmtAnodeId

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

Private Member Functions
virtual void	beginRun (const edm::Run &, const edm::EventSetup &) override

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

virtual void	produce (edm::Event &, const edm::EventSetup &) override

unsigned	sortDataByPmt ()

Private Attributes
bool	dropZSmarkedPassed_

edm::InputTag	inputLabel_

std::unique_ptr< HcalRecoParams >	paramTS_

std::vector< HFQIE10Info >	qie10Infos_

HFPreRecAlgo	reco_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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 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)

Detailed Description

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

Implementation: [Notes on implementation]

Definition at line 48 of file HFPreReconstructor.cc.

Member Typedef Documentation

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

private

Definition at line 57 of file HFPreReconstructor.cc.

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

private

Definition at line 58 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")),
       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 ( )

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
	)

overrideprivatevirtual

Definition at line 212 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 293 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<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 146 of file HFPreReconstructor.cc.

References edm::DataFrameContainer::begin(), AlignmentProducer_cff::calibrations, dropZSmarkedPassed_, edm::DataFrameContainer::end(), HcalRecoParam::firstSample(), edm::EventSetup::get(), edm::Event::getByToken(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalForward, QIE10DataFrame::id(), paramTS_, QIE10DataFrame::presamples(), edm::ESHandle< T >::product(), qie10Infos_, reco_, HFPreRecAlgo::reconstruct(), edm::DataFrameContainer::size(), 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);
 
             // Get the "sample of interest" from the data frame itself
             int tsToUse = frame.presamples();
             if (tsFromDB_)
             {
                 const HcalRecoParam* param_ts = paramTS_->getValues(cell.rawId());
                 tsToUse = param_ts->firstSample();
             }
 
             // Reconstruct the charge, energy, etc
             qie10Infos_.push_back(reco_.reconstruct(frame, tsToUse, coder, calibrations));
         }
     }
 }

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

overrideprivatevirtual

Definition at line 224 of file HFPreReconstructor.cc.

References fillInfos(), plotBeamSpotDB::first, 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 = 0;
                 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 = 0;
                 }
 
                 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 112 of file HFPreReconstructor.cc.

References particleFlowClusterECALTimeSelected_cfi::depth, plotBeamSpotDB::first, i, hcalTTPDigis_cfi::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;
 }