ShallowClustersProducer Class Reference

#include <ShallowClustersProducer.h>

Inheritance diagram for ShallowClustersProducer:

Classes
struct	moduleVars
struct	NearDigis

Public Member Functions
	ShallowClustersProducer (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
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 ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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 &	itemsToGetFrom (BranchType iType) 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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	Prefix
SiStripClusterInfo	siStripClusterInfo_
edm::InputTag	theClustersLabel
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	theClustersToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > >	theDigisToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 17 of file ShallowClustersProducer.h.

Constructor & Destructor Documentation

ShallowClustersProducer()

ShallowClustersProducer::ShallowClustersProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 12 of file ShallowClustersProducer.cc.

    : Prefix(iConfig.getParameter<std::string>("Prefix")), siStripClusterInfo_(consumesCollector()) {
  produces<std::vector<unsigned>>(Prefix + "number");
  produces<std::vector<unsigned>>(Prefix + "width");
  produces<std::vector<float>>(Prefix + "variance");
  produces<std::vector<float>>(Prefix + "barystrip");
  produces<std::vector<float>>(Prefix + "middlestrip");
  produces<std::vector<unsigned>>(Prefix + "charge");
  produces<std::vector<float>>(Prefix + "noise");
  produces<std::vector<float>>(Prefix + "ston");
  produces<std::vector<unsigned>>(Prefix + "seedstrip");
  produces<std::vector<unsigned>>(Prefix + "seedindex");
  produces<std::vector<unsigned>>(Prefix + "seedcharge");
  produces<std::vector<float>>(Prefix + "seednoise");
  produces<std::vector<float>>(Prefix + "seedgain");
  produces<std::vector<unsigned>>(Prefix + "qualityisbad");
 
  produces<std::vector<float>>(Prefix + "rawchargeC");
  produces<std::vector<float>>(Prefix + "rawchargeL");
  produces<std::vector<float>>(Prefix + "rawchargeR");
  produces<std::vector<float>>(Prefix + "rawchargeLL");
  produces<std::vector<float>>(Prefix + "rawchargeRR");
  produces<std::vector<float>>(Prefix + "eta");
  produces<std::vector<float>>(Prefix + "foldedeta");
  produces<std::vector<float>>(Prefix + "etaX");
  produces<std::vector<float>>(Prefix + "etaasymm");
  produces<std::vector<float>>(Prefix + "outsideasymm");
  produces<std::vector<float>>(Prefix + "neweta");
  produces<std::vector<float>>(Prefix + "newetaerr");
 
  produces<std::vector<unsigned>>(Prefix + "detid");
  produces<std::vector<int>>(Prefix + "subdetid");
  produces<std::vector<int>>(Prefix + "module");
  produces<std::vector<int>>(Prefix + "side");
  produces<std::vector<int>>(Prefix + "layerwheel");
  produces<std::vector<int>>(Prefix + "stringringrod");
  produces<std::vector<int>>(Prefix + "petal");
  produces<std::vector<int>>(Prefix + "stereo");
 
  theClustersToken_ = consumes<edmNew::DetSetVector<SiStripCluster>>(iConfig.getParameter<edm::InputTag>("Clusters"));
  theDigisToken_ = consumes<edm::DetSetVector<SiStripProcessedRawDigi>>(edm::InputTag("siStripProcessedRawDigis", ""));
}

References edm::ParameterSet::getParameter(), HLT_2018_cff::InputTag, Prefix, theClustersToken_, and theDigisToken_.

Member Function Documentation

produce()

void ShallowClustersProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 55 of file ShallowClustersProducer.cc.

                                                                                   {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();
 
  siStripClusterInfo_.initEvent(iSetup);
 
  auto number = std::make_unique<std::vector<unsigned>>(7, 0);
  auto width = std::make_unique<std::vector<unsigned>>();
  auto variance = std::make_unique<std::vector<float>>();
  auto barystrip = std::make_unique<std::vector<float>>();
  auto middlestrip = std::make_unique<std::vector<float>>();
  auto charge = std::make_unique<std::vector<unsigned>>();
  auto noise = std::make_unique<std::vector<float>>();
  auto ston = std::make_unique<std::vector<float>>();
  auto seedstrip = std::make_unique<std::vector<unsigned>>();
  auto seedindex = std::make_unique<std::vector<unsigned>>();
  auto seedcharge = std::make_unique<std::vector<unsigned>>();
  auto seednoise = std::make_unique<std::vector<float>>();
  auto seedgain = std::make_unique<std::vector<float>>();
  auto qualityisbad = std::make_unique<std::vector<unsigned>>();
 
  auto rawchargeC = std::make_unique<std::vector<float>>();
  auto rawchargeL = std::make_unique<std::vector<float>>();
  auto rawchargeR = std::make_unique<std::vector<float>>();
  auto rawchargeLL = std::make_unique<std::vector<float>>();
  auto rawchargeRR = std::make_unique<std::vector<float>>();
  auto etaX = std::make_unique<std::vector<float>>();
  auto eta = std::make_unique<std::vector<float>>();
  auto foldedeta = std::make_unique<std::vector<float>>();
  auto etaasymm = std::make_unique<std::vector<float>>();
  auto outsideasymm = std::make_unique<std::vector<float>>();
  auto neweta = std::make_unique<std::vector<float>>();
  auto newetaerr = std::make_unique<std::vector<float>>();
 
  auto detid = std::make_unique<std::vector<unsigned>>();
  auto subdetid = std::make_unique<std::vector<int>>();
  auto side = std::make_unique<std::vector<int>>();
  auto module = std::make_unique<std::vector<int>>();
  auto layerwheel = std::make_unique<std::vector<int>>();
  auto stringringrod = std::make_unique<std::vector<int>>();
  auto petal = std::make_unique<std::vector<int>>();
  auto stereo = std::make_unique<std::vector<int>>();
 
  edm::Handle<edmNew::DetSetVector<SiStripCluster>> clusters;
  //  iEvent.getByLabel(theClustersLabel, clusters);
  iEvent.getByToken(theClustersToken_, clusters);
 
  edm::Handle<edm::DetSetVector<SiStripProcessedRawDigi>> rawProcessedDigis;
  //  iEvent.getByLabel("siStripProcessedRawDigis", "", rawProcessedDigis);
  iEvent.getByToken(theDigisToken_, rawProcessedDigis);
 
  edmNew::DetSetVector<SiStripCluster>::const_iterator itClusters = clusters->begin();
  for (; itClusters != clusters->end(); ++itClusters) {
    uint32_t id = itClusters->id();
    const moduleVars moduleV(id, tTopo);
    for (edmNew::DetSet<SiStripCluster>::const_iterator cluster = itClusters->begin(); cluster != itClusters->end();
         ++cluster) {
      siStripClusterInfo_.setCluster(*cluster, id);
      const SiStripClusterInfo& info = siStripClusterInfo_;
      const NearDigis digis = rawProcessedDigis.isValid() ? NearDigis(info, *rawProcessedDigis) : NearDigis(info);
 
      (number->at(0))++;
      (number->at(moduleV.subdetid))++;
      width->push_back(cluster->amplitudes().size());
      barystrip->push_back(cluster->barycenter());
      variance->push_back(info.variance());
      middlestrip->push_back(info.firstStrip() + info.width() / 2.0);
      charge->push_back(info.charge());
      noise->push_back(info.noiseRescaledByGain());
      ston->push_back(info.signalOverNoise());
      seedstrip->push_back(info.maxStrip());
      seedindex->push_back(info.maxIndex());
      seedcharge->push_back(info.maxCharge());
      seednoise->push_back(info.stripNoisesRescaledByGain().at(info.maxIndex()));
      seedgain->push_back(info.stripGains().at(info.maxIndex()));
      qualityisbad->push_back(info.IsAnythingBad());
 
      rawchargeC->push_back(digis.max);
      rawchargeL->push_back(digis.left);
      rawchargeR->push_back(digis.right);
      rawchargeLL->push_back(digis.Lleft);
      rawchargeRR->push_back(digis.Rright);
      etaX->push_back(digis.etaX());
      eta->push_back(digis.eta());
      etaasymm->push_back(digis.etaasymm());
      outsideasymm->push_back(digis.outsideasymm());
      neweta->push_back((digis.last - digis.first) / info.charge());
      newetaerr->push_back((sqrt(digis.last + digis.first)) / pow(info.charge(), 1.5));
 
      detid->push_back(id);
      subdetid->push_back(moduleV.subdetid);
      side->push_back(moduleV.side);
      module->push_back(moduleV.module);
      layerwheel->push_back(moduleV.layerwheel);
      stringringrod->push_back(moduleV.stringringrod);
      petal->push_back(moduleV.petal);
      stereo->push_back(moduleV.stereo);
    }
  }
 
  iEvent.put(std::move(number), Prefix + "number");
  iEvent.put(std::move(width), Prefix + "width");
  iEvent.put(std::move(variance), Prefix + "variance");
  iEvent.put(std::move(barystrip), Prefix + "barystrip");
  iEvent.put(std::move(middlestrip), Prefix + "middlestrip");
  iEvent.put(std::move(charge), Prefix + "charge");
  iEvent.put(std::move(noise), Prefix + "noise");
  iEvent.put(std::move(ston), Prefix + "ston");
  iEvent.put(std::move(seedstrip), Prefix + "seedstrip");
  iEvent.put(std::move(seedindex), Prefix + "seedindex");
  iEvent.put(std::move(seedcharge), Prefix + "seedcharge");
  iEvent.put(std::move(seednoise), Prefix + "seednoise");
  iEvent.put(std::move(seedgain), Prefix + "seedgain");
  iEvent.put(std::move(qualityisbad), Prefix + "qualityisbad");
 
  iEvent.put(std::move(rawchargeC), Prefix + "rawchargeC");
  iEvent.put(std::move(rawchargeL), Prefix + "rawchargeL");
  iEvent.put(std::move(rawchargeR), Prefix + "rawchargeR");
  iEvent.put(std::move(rawchargeLL), Prefix + "rawchargeLL");
  iEvent.put(std::move(rawchargeRR), Prefix + "rawchargeRR");
  iEvent.put(std::move(etaX), Prefix + "etaX");
  iEvent.put(std::move(eta), Prefix + "eta");
  iEvent.put(std::move(foldedeta), Prefix + "foldedeta");
  iEvent.put(std::move(etaasymm), Prefix + "etaasymm");
  iEvent.put(std::move(outsideasymm), Prefix + "outsideasymm");
  iEvent.put(std::move(neweta), Prefix + "neweta");
  iEvent.put(std::move(newetaerr), Prefix + "newetaerr");
 
  iEvent.put(std::move(detid), Prefix + "detid");
  iEvent.put(std::move(subdetid), Prefix + "subdetid");
  iEvent.put(std::move(module), Prefix + "module");
  iEvent.put(std::move(side), Prefix + "side");
  iEvent.put(std::move(layerwheel), Prefix + "layerwheel");
  iEvent.put(std::move(stringringrod), Prefix + "stringringrod");
  iEvent.put(std::move(petal), Prefix + "petal");
  iEvent.put(std::move(stereo), Prefix + "stereo");
}

References edmNew::DetSetVector< T >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, bsc_activity_cfg::clusters, edmNew::DetSetVector< T >::end(), ShallowClustersProducer::NearDigis::eta(), PVValHelper::eta, ShallowClustersProducer::NearDigis::etaasymm(), ShallowClustersProducer::NearDigis::etaX(), ShallowClustersProducer::NearDigis::first, edm::EventSetup::get(), get, edmNew::DetSetVector< T >::id(), iEvent, info(), SiStripClusterInfo::initEvent(), edm::HandleBase::isValid(), ShallowClustersProducer::NearDigis::last, ShallowClustersProducer::moduleVars::layerwheel, ShallowClustersProducer::NearDigis::left, ShallowClustersProducer::NearDigis::Lleft, ShallowClustersProducer::NearDigis::max, ShallowClustersProducer::moduleVars::module, eostools::move(), hgcalDigitizer_cfi::noise, contentValuesFiles::number, ShallowClustersProducer::NearDigis::outsideasymm(), ShallowClustersProducer::moduleVars::petal, funct::pow(), Prefix, edm::ESHandle< T >::product(), ShallowClustersProducer::NearDigis::right, ShallowClustersProducer::NearDigis::Rright, SiStripClusterInfo::setCluster(), ShallowClustersProducer::moduleVars::side, siStripClusterInfo_, mathSSE::sqrt(), ShallowClustersProducer::moduleVars::stereo, ShallowClustersProducer::moduleVars::stringringrod, ShallowClustersProducer::moduleVars::subdetid, theClustersToken_, theDigisToken_, and ApeEstimator_cff::width.

Member Data Documentation

Prefix

std::string ShallowClustersProducer::Prefix

private

Definition at line 23 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().

siStripClusterInfo_

SiStripClusterInfo ShallowClustersProducer::siStripClusterInfo_

private

Definition at line 44 of file ShallowClustersProducer.h.

Referenced by produce().

theClustersLabel

edm::InputTag ShallowClustersProducer::theClustersLabel

private

Definition at line 22 of file ShallowClustersProducer.h.

theClustersToken_

edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > ShallowClustersProducer::theClustersToken_

private

Definition at line 42 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().

theDigisToken_

edm::EDGetTokenT<edm::DetSetVector<SiStripProcessedRawDigi> > ShallowClustersProducer::theDigisToken_

private

Definition at line 43 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().