#include <SiStripClusterizer.h>

Inheritance diagram for SiStripClusterizer:

Public Member Functions
void	beginRun (const edm::Run &run, const edm::EventSetup &es) override

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

	SiStripClusterizer (const edm::ParameterSet &conf)

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef edm::EDGetTokenT < edm::DetSetVector < SiStripDigi > >	token_t

typedef std::vector< token_t >	token_v

Private Member Functions
template<class T >
bool	findInput (const edm::EDGetTokenT< T > &, edm::Handle< T > &, const edm::Event &)

void	refineCluster (const edm::Handle< edm::DetSetVector< SiStripDigi > > &input, std::auto_ptr< edmNew::DetSetVector< SiStripCluster > > &output)

Private Attributes
std::auto_ptr < StripClusterizerAlgorithm >	algorithm

edm::ParameterSet	confClusterizer_

bool	doRefineCluster_

const std::vector< edm::InputTag >	inputTags

token_v	inputTokens

float	occupancyThreshold_

edm::ESHandle< SiStripQuality >	quality_

edm::ESHandle< SiStripDetCabling >	SiStripDetCabling_

unsigned	widthThreshold_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Definition at line 15 of file SiStripClusterizer.h.

Member Typedef Documentation

typedef edm::EDGetTokenT< edm::DetSetVector<SiStripDigi> > SiStripClusterizer::token_t

private

Definition at line 34 of file SiStripClusterizer.h.

typedef std::vector<token_t> SiStripClusterizer::token_v

private

Definition at line 35 of file SiStripClusterizer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file SiStripClusterizer.cc.

References confClusterizer_, doRefineCluster_, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), inputTags, inputTokens, occupancyThreshold_, GlobalPosition_Frontier_DevDB_cff::tag, edm::vector_transform(), and widthThreshold_.

   : confClusterizer_(conf.getParameter<edm::ParameterSet>("Clusterizer")),
     inputTags( conf.getParameter<std::vector<edm::InputTag> >("DigiProducersList") ),
     algorithm( StripClusterizerAlgorithmFactory::create(conf.getParameter<edm::ParameterSet>("Clusterizer")) ) {
   produces< edmNew::DetSetVector<SiStripCluster> > ();
   inputTokens = edm::vector_transform( inputTags, [this](edm::InputTag const & tag) { return consumes< edm::DetSetVector<SiStripDigi> >(tag);} );
   doRefineCluster_ = confClusterizer_.existsAs<bool>("doRefineCluster") ? confClusterizer_.getParameter<bool>("doRefineCluster") : false;
   occupancyThreshold_ = confClusterizer_.existsAs<double>("occupancyThreshold") ? confClusterizer_.getParameter<double>("occupancyThreshold") : 0.05;
   widthThreshold_ = confClusterizer_.existsAs<unsigned>("widthThreshold") ? confClusterizer_.getParameter<unsigned>("widthThreshold") : 4;
 }

Member Function Documentation

void SiStripClusterizer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

overridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 84 of file SiStripClusterizer.cc.

References doRefineCluster_, edm::EventSetup::get(), quality_, and SiStripDetCabling_.

                                                                              {
   if (doRefineCluster_) {
     es.get<SiStripDetCablingRcd>().get( SiStripDetCabling_);
     es.get<SiStripQualityRcd>().get("", quality_);
   }
 }

template<class T >

bool SiStripClusterizer::findInput	(	const edm::EDGetTokenT< T > &	tag,
		edm::Handle< T > &	handle,
		const edm::Event &	e
	)

inlineprivate

Definition at line 51 of file SiStripClusterizer.cc.

References edm::Event::getByToken(), and edm::HandleBase::isValid().

Referenced by produce().

                                                                                  {
     e.getByToken( tag, handle);
     return handle.isValid();
 }

void SiStripClusterizer::produce	(	edm::Event &	event,
		const edm::EventSetup &	es
	)

virtual

Implements edm::stream::EDProducerBase.

Definition at line 23 of file SiStripClusterizer.cc.

References algorithm, doRefineCluster_, findInput(), inputTokens, LogDebug, convertSQLitetoXML_cfg::output, and refineCluster().

                                                    {
 
   std::auto_ptr< edmNew::DetSetVector<SiStripCluster> > output(new edmNew::DetSetVector<SiStripCluster>());
   output->reserve(10000,4*10000);
 
   edm::Handle< edm::DetSetVector<SiStripDigi> >     inputOld;  
 //   edm::Handle< edmNew::DetSetVector<SiStripDigi> >  inputNew;  
 
   algorithm->initialize(es);  
 
   BOOST_FOREACH( const edm::EDGetTokenT< edm::DetSetVector<SiStripDigi> >& token, inputTokens) {
     if(      findInput( token, inputOld, event) ) {
       algorithm->clusterize(*inputOld, *output);
       if (doRefineCluster_) refineCluster(inputOld, output);
     } 
 //     else if( findInput( tag, inputNew, event) ) algorithm->clusterize(*inputNew, *output);
     else edm::LogError("Input Not Found") << "[SiStripClusterizer::produce] ";// << tag;
   }
 
   LogDebug("Output") << output->dataSize() << " clusters from " 
              << output->size()     << " modules";
   output->shrink_to_fit();
   event.put(output);
 }

void SiStripClusterizer::refineCluster	(	const edm::Handle< edm::DetSetVector< SiStripDigi > > &	input,
		std::auto_ptr< edmNew::DetSetVector< SiStripCluster > > &	output
	)

private

Definition at line 57 of file SiStripClusterizer.cc.

References input, occupancyThreshold_, convertSQLitetoXML_cfg::output, quality_, SiStripDetCabling_, edm::DetSetVector< T >::size(), and widthThreshold_.

Referenced by produce().

                                                                      {
   if (input->size() == 0) return;
 
   // Flag merge-prone clusters for relaxed CPE errors
   // Criterion is sensor occupancy and cluster width exceeding thresholds
 
   for (edmNew::DetSetVector<SiStripCluster>::const_iterator det=output->begin(); det!=output->end(); det++) {
     uint32_t detId = det->id();
     // Find the number of good strips in this sensor
     int nchannideal = SiStripDetCabling_->nApvPairs(detId) * 2 * 128;
     int nchannreal = 0;
     for(int strip = 0; strip < nchannideal; ++strip)
       if(!quality_->IsStripBad(detId,strip)) ++nchannreal;
 
     edm::DetSetVector<SiStripDigi>::const_iterator digis = input->find(detId);
     if (digis != input->end()) {
       int ndigi = digis->size();
       for (edmNew::DetSet<SiStripCluster>::iterator clust = det->begin(); clust != det->end(); clust++) {
     if (ndigi > occupancyThreshold_*nchannreal && clust->amplitudes().size() >= widthThreshold_) clust->setMerged(true);
     else clust->setMerged(false);
       }
       // std::cout << "Sensor:strips_occStrips_clust " << nchannreal << " " << ndigi << " " << det->size() << std::endl;
     }
   }  // traverse sensors
 }