ShallowClustersProducer Class Reference

Inheritance diagram for ShallowClustersProducer:

Classes
struct	moduleVars
struct	NearDigis

Public Member Functions
	ShallowClustersProducer (const edm::ParameterSet &)
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
std::string	Prefix
SiStripClusterInfo	siStripClusterInfo_
edm::InputTag	theClustersLabel
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	theClustersToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > >	theDigisToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

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

Definition at line 18 of file ShallowClustersProducer.cc.

Constructor & Destructor Documentation

ShallowClustersProducer()

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

explicit

Definition at line 52 of file ShallowClustersProducer.cc.

    : Prefix(iConfig.getParameter<std::string>("Prefix")),
      topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>()),
      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_FULL_cff::InputTag, Prefix, theClustersToken_, and theDigisToken_.

Member Function Documentation

produce()

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

overrideprivate

Definition at line 97 of file ShallowClustersProducer.cc.

                                                                                   {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(topoToken_);
  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::getHandle(), 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_, topoToken_, and ApeEstimator_cff::width.

Member Data Documentation

Prefix

std::string ShallowClustersProducer::Prefix

private

Definition at line 24 of file ShallowClustersProducer.cc.

Referenced by produce(), and ShallowClustersProducer().

siStripClusterInfo_

SiStripClusterInfo ShallowClustersProducer::siStripClusterInfo_

private

Definition at line 46 of file ShallowClustersProducer.cc.

Referenced by produce().

theClustersLabel

edm::InputTag ShallowClustersProducer::theClustersLabel

private

Definition at line 23 of file ShallowClustersProducer.cc.

theClustersToken_

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

private

Definition at line 44 of file ShallowClustersProducer.cc.

Referenced by produce(), and ShallowClustersProducer().

theDigisToken_

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

private

Definition at line 45 of file ShallowClustersProducer.cc.

Referenced by produce(), and ShallowClustersProducer().

topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> ShallowClustersProducer::topoToken_

private

Definition at line 43 of file ShallowClustersProducer.cc.

Referenced by produce().