Inheritance diagram for SiStripClusters2ApproxClusters:

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

	SiStripClusters2ApproxClusters (const edm::ParameterSet &conf)

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

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

Private Attributes
edm::InputTag	beamSpot_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	clusterToken

std::string	csfLabel_

edm::ESGetToken< ClusterShapeHitFilter, CkfComponentsRecord >	csfToken_

SiStripDetInfo	detInfo_

edm::FileInPath	fileInPath_

edm::InputTag	inputClusters

unsigned int	maxNSat

edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	stripNoiseToken_

edm::ESHandle< SiStripNoises >	theNoise_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

Static Private Attributes
static constexpr double	seedCutMIPs_ = .35

static constexpr double	seedCutSN_ = 7.

static constexpr double	subclusterCutMIPs_ = .45

static constexpr double	subclusterCutSN_ = 12.

static constexpr double	subclusterWindow_ = .7

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 36 of file SiStripClusters2ApproxClusters.cc.

Constructor & Destructor Documentation

◆ SiStripClusters2ApproxClusters()

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

explicit

Definition at line 69 of file SiStripClusters2ApproxClusters.cc.

References beamSpot_, beamSpotToken_, clusterToken, csfLabel_, csfToken_, detInfo_, deDxTools::esConsumes(), fileInPath_, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), inputClusters, SiStripDetInfoFileReader::kDefaultFile, maxNSat, SiStripDetInfoFileReader::read(), AlCaHLTBitMon_QueryRunRegistry::string, stripNoiseToken_, and tkGeomToken_.

                                                                                           {
   inputClusters = conf.getParameter<edm::InputTag>("inputClusters");
   maxNSat = conf.getParameter<unsigned int>("maxSaturatedStrips");
 
   clusterToken = consumes<edmNew::DetSetVector<SiStripCluster> >(inputClusters);
 
   beamSpot_ = conf.getParameter<edm::InputTag>("beamSpot");
   beamSpotToken_ = consumes<reco::BeamSpot>(beamSpot_);
 
   tkGeomToken_ = esConsumes();
 
   fileInPath_ = edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile);
   detInfo_ = SiStripDetInfoFileReader::read(fileInPath_.fullPath());
 
   csfLabel_ = conf.getParameter<std::string>("clusterShapeHitFilterLabel");
   csfToken_ = esConsumes(edm::ESInputTag("", csfLabel_));
 
   stripNoiseToken_ = esConsumes();
   produces<SiStripApproximateClusterCollection>();
 }

Member Function Documentation

◆ fillDescriptions()

void SiStripClusters2ApproxClusters::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 158 of file SiStripClusters2ApproxClusters.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("inputClusters", edm::InputTag("siStripClusters"));
   desc.add<unsigned int>("maxSaturatedStrips", 3);
   desc.add<std::string>("clusterShapeHitFilterLabel", "ClusterShapeHitFilter");  // add CSF label
   desc.add<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot"));         // add BeamSpot tag
   descriptions.add("SiStripClusters2ApproxClusters", desc);
 }

◆ produce()

void SiStripClusters2ApproxClusters::produce	(	edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 90 of file SiStripClusters2ApproxClusters.cc.

References funct::abs(), ConfigBuilder::beamspot, align::BeamSpot, beamSpot_, beamSpotToken_, Surface::bounds(), cms::cuda::bs, reco::ceil(), ALCARECOSiPixelCalSingleMuonDQM_cff::clusterCollection, clusterToken, csfToken_, hcalRecHitTable_cff::detId, detInfo_, edm::InputTag::encode(), f, alignCSCRings::ff, edm::EventSetup::getData(), SiStripDetInfo::getNumberOfApvsAndStripLength(), maxNSat, sistrip::MeVperADCStrip, eostools::move(), hcaltpdigi_cfi::peakFilter, packedPFCandidateRefMixer_cfi::pf, mps_fire::result, seedCutMIPs_, seedCutSN_, stripNoiseToken_, sistrip::STRIPS_PER_APV, subclusterCutMIPs_, subclusterCutSN_, subclusterWindow_, GeomDet::surface(), builder_last_value_cfg::test, theNoise_, Bounds::thickness(), tkGeomToken_, Surface::toGlobal(), GeomDet::toLocal(), Bounds::width(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                          {
   const auto& clusterCollection = event.get(clusterToken);
   auto result = std::make_unique<SiStripApproximateClusterCollection>();
   result->reserve(clusterCollection.size(), clusterCollection.dataSize());
 
   auto const beamSpotHandle = event.getHandle(beamSpotToken_);
   auto const& bs = beamSpotHandle.isValid() ? *beamSpotHandle : reco::BeamSpot();
   if (not beamSpotHandle.isValid()) {
     edm::LogError("SiStripClusters2ApproxClusters")
         << "didn't find a valid beamspot with label \"" << beamSpot_.encode() << "\" -> using (0,0,0)";
   }
 
   const auto& tkGeom = &iSetup.getData(tkGeomToken_);
   const auto& theFilter = &iSetup.getData(csfToken_);
   const auto& theNoise_ = &iSetup.getData(stripNoiseToken_);
 
   for (const auto& detClusters : clusterCollection) {
     auto ff = result->beginDet(detClusters.id());
 
     unsigned int detId = detClusters.id();
     const GeomDet* det = tkGeom->idToDet(detId);
     double nApvs = detInfo_.getNumberOfApvsAndStripLength(detId).first;
     double stripLength = detInfo_.getNumberOfApvsAndStripLength(detId).second;
     double barycenter_ypos = 0.5 * stripLength;
 
     const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(det);
     float mip = 3.9 / (sistrip::MeVperADCStrip / stripDet->surface().bounds().thickness());
 
     for (const auto& cluster : detClusters) {
       const LocalPoint& lp = LocalPoint(((cluster.barycenter() * 10 / (sistrip::STRIPS_PER_APV * nApvs)) -
                                          ((stripDet->surface().bounds().width()) * 0.5f)),
                                         barycenter_ypos - (0.5f * stripLength),
                                         0.);
       const GlobalPoint& gpos = det->surface().toGlobal(lp);
       GlobalPoint beamspot(bs.position().x(), bs.position().y(), bs.position().z());
       const GlobalVector& gdir = gpos - beamspot;
       const LocalVector& ldir = det->toLocal(gdir);
 
       int hitStrips;
       float hitPredPos;
       bool usable = theFilter->getSizes(detId, cluster, lp, ldir, hitStrips, hitPredPos);
       // (almost) same logic as in StripSubClusterShapeTrajectoryFilter
       bool isTrivial = (std::abs(hitPredPos) < 2.f && hitStrips <= 2);
 
       if (!usable || isTrivial) {
         ff.push_back(SiStripApproximateCluster(cluster, maxNSat, hitPredPos, true));
       } else {
         bool peakFilter = false;
         SlidingPeakFinder pf(std::max<int>(2, std::ceil(std::abs(hitPredPos) + subclusterWindow_)));
         float mipnorm = mip / std::abs(ldir.z());
         PeakFinderTest test(mipnorm,
                             detId,
                             cluster.firstStrip(),
                             theNoise_,
                             seedCutMIPs_,
                             seedCutSN_,
                             subclusterCutMIPs_,
                             subclusterCutSN_);
         peakFilter = pf.apply(cluster.amplitudes(), test);
 
         ff.push_back(SiStripApproximateCluster(cluster, maxNSat, hitPredPos, peakFilter));
       }
     }
   }
 
   event.put(std::move(result));
 }