#include <ShallowRechitClustersProducer.h>

Inheritance diagram for ShallowRechitClustersProducer:

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

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

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

Private Attributes
const edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	clusters_token_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

std::string	Prefix

std::vector< edm::EDGetTokenT< SiStripRecHit2DCollection > >	rec_hits_tokens_

std::string	Suffix

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

Definition at line 13 of file ShallowRechitClustersProducer.h.

Constructor & Destructor Documentation

◆ ShallowRechitClustersProducer()

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

explicit

Definition at line 11 of file ShallowRechitClustersProducer.cc.

References edm::ParameterSet::getParameter(), Prefix, rec_hits_tokens_, and Suffix.

     : Suffix(iConfig.getParameter<std::string>("Suffix")),
       Prefix(iConfig.getParameter<std::string>("Prefix")),
       geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
       clusters_token_(consumes<edmNew::DetSetVector<SiStripCluster>>(iConfig.getParameter<edm::InputTag>("Clusters"))) {
   std::vector<edm::InputTag> rec_hits_tags = iConfig.getParameter<std::vector<edm::InputTag>>("InputTags");
   for (const auto& itag : rec_hits_tags) {
     rec_hits_tokens_.push_back(consumes<SiStripRecHit2DCollection>(itag));
   }
 
   produces<std::vector<float>>(Prefix + "strip" + Suffix);
   produces<std::vector<float>>(Prefix + "merr" + Suffix);
   produces<std::vector<float>>(Prefix + "localx" + Suffix);
   produces<std::vector<float>>(Prefix + "localy" + Suffix);
   produces<std::vector<float>>(Prefix + "localxerr" + Suffix);
   produces<std::vector<float>>(Prefix + "localyerr" + Suffix);
   produces<std::vector<float>>(Prefix + "globalx" + Suffix);
   produces<std::vector<float>>(Prefix + "globaly" + Suffix);
   produces<std::vector<float>>(Prefix + "globalz" + Suffix);
 }

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 32 of file ShallowRechitClustersProducer.cc.

References clusters_token_, Exception, geomToken_, edm::EventSetup::getHandle(), mps_fire::i, iEvent, shallow::make_cluster_map(), Topology::measurementError(), eostools::move(), Prefix, rec_hits_tokens_, FastTrackerRecHitMaskProducer_cfi::recHits, findQualityFiles::size, StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), digitizers_cfi::strip, StripTopology::strip(), Suffix, GeomDet::toGlobal(), x, hit::x, y, hit::y, and z.

                                                                                          {
   shallow::CLUSTERMAP clustermap = shallow::make_cluster_map(iEvent, clusters_token_);
 
   int size = clustermap.size();
   auto strip = std::make_unique<std::vector<float>>(size, -10000);
   auto merr = std::make_unique<std::vector<float>>(size, -10000);
   auto localx = std::make_unique<std::vector<float>>(size, -10000);
   auto localy = std::make_unique<std::vector<float>>(size, -10000);
   auto localxerr = std::make_unique<std::vector<float>>(size, -1);
   auto localyerr = std::make_unique<std::vector<float>>(size, -1);
   auto globalx = std::make_unique<std::vector<float>>(size, -10000);
   auto globaly = std::make_unique<std::vector<float>>(size, -10000);
   auto globalz = std::make_unique<std::vector<float>>(size, -10000);
 
   edm::ESHandle<TrackerGeometry> theTrackerGeometry = iSetup.getHandle(geomToken_);
 
   for (auto recHit_token : rec_hits_tokens_) {
     edm::Handle<SiStripRecHit2DCollection> recHits;
     iEvent.getByToken(recHit_token, recHits);
     for (auto const& ds : *recHits) {
       for (auto const& hit : ds) {
         shallow::CLUSTERMAP::iterator cluster =
             clustermap.find(std::make_pair(hit.geographicalId().rawId(), hit.cluster()->firstStrip()));
         if (cluster != clustermap.end()) {
           const StripGeomDetUnit* theStripDet =
               dynamic_cast<const StripGeomDetUnit*>(theTrackerGeometry->idToDet(hit.geographicalId()));
           unsigned int i = cluster->second;
           strip->at(i) = theStripDet->specificTopology().strip(hit.localPosition());
           merr->at(i) = sqrt(
               theStripDet->specificTopology().measurementError(hit.localPosition(), hit.localPositionError()).uu());
           localx->at(i) = hit.localPosition().x();
           localy->at(i) = hit.localPosition().y();
           localxerr->at(i) = sqrt(hit.localPositionError().xx());
           localyerr->at(i) = sqrt(hit.localPositionError().yy());
           globalx->at(i) = theStripDet->toGlobal(hit.localPosition()).x();
           globaly->at(i) = theStripDet->toGlobal(hit.localPosition()).y();
           globalz->at(i) = theStripDet->toGlobal(hit.localPosition()).z();
         } else {
           throw cms::Exception("cluster not found");
         }
       }
     }
   }
 
   iEvent.put(std::move(strip), Prefix + "strip" + Suffix);
   iEvent.put(std::move(merr), Prefix + "merr" + Suffix);
   iEvent.put(std::move(localx), Prefix + "localx" + Suffix);
   iEvent.put(std::move(localy), Prefix + "localy" + Suffix);
   iEvent.put(std::move(localxerr), Prefix + "localxerr" + Suffix);
   iEvent.put(std::move(localyerr), Prefix + "localyerr" + Suffix);
   iEvent.put(std::move(globalx), Prefix + "globalx" + Suffix);
   iEvent.put(std::move(globaly), Prefix + "globaly" + Suffix);
   iEvent.put(std::move(globalz), Prefix + "globalz" + Suffix);
 }