#include <SiStripRecHitConverterAlgorithm.h>

Classes
struct	products

Public Member Functions
void	initialize (const edm::EventSetup &)

void	run (edm::Handle< edmNew::DetSetVector< SiStripCluster > > input, products &output)

void	run (edm::Handle< edmNew::DetSetVector< SiStripCluster > > input, products &output, LocalVector trackdirection)

	SiStripRecHitConverterAlgorithm (const edm::ParameterSet &, edm::ConsumesCollector)

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

Private Types
typedef SiStripRecHit2DCollection::FastFiller	Collector

Private Member Functions
void	fillBad128StripBlocks (const uint32_t detid, bool bad128StripBlocks[6]) const

bool	isMasked (const SiStripCluster &cluster, bool bad128StripBlocks[6]) const

void	match (products &output, LocalVector trackdirection) const

bool	useModule (const uint32_t id) const

Private Attributes
edm::ESGetToken < StripClusterParameterEstimator, TkStripCPERecord >	cpeToken

bool	doMatching

bool	maskBad128StripBlocks

const SiStripRecHitMatcher *	matcher = nullptr

edm::ESGetToken < SiStripRecHitMatcher, TkStripCPERecord >	matcherToken

const StripClusterParameterEstimator *	parameterestimator = nullptr

const SiStripQuality *	quality = nullptr

edm::ESGetToken < SiStripQuality, SiStripQualityRcd >	qualityToken

const TrackerGeometry *	tracker = nullptr

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerToken

bool	useQuality

Detailed Description

Definition at line 25 of file SiStripRecHitConverterAlgorithm.h.

Member Typedef Documentation

typedef SiStripRecHit2DCollection::FastFiller SiStripRecHitConverterAlgorithm::Collector

private

Definition at line 69 of file SiStripRecHitConverterAlgorithm.h.

Constructor & Destructor Documentation

SiStripRecHitConverterAlgorithm::SiStripRecHitConverterAlgorithm	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector	iC
	)

Definition at line 16 of file SiStripRecHitConverterAlgorithm.cc.

References doMatching, edm::ConsumesCollector::esConsumes(), edm::ParameterSet::getParameter(), matcherToken, qualityToken, AlCaHarvesting_cff::SiStripQuality, and useQuality.

     : useQuality(conf.getParameter<bool>("useSiStripQuality")),
       maskBad128StripBlocks(conf.getParameter<bool>("MaskBadAPVFibers")),
       doMatching(conf.getParameter<bool>("doMatching")),
       trackerToken(iC.esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
       cpeToken(iC.esConsumes<StripClusterParameterEstimator, TkStripCPERecord>(
           conf.getParameter<edm::ESInputTag>("StripCPE"))) {
   if (doMatching) {
     matcherToken = iC.esConsumes<SiStripRecHitMatcher, TkStripCPERecord>(conf.getParameter<edm::ESInputTag>("Matcher"));
   }
   if (useQuality) {
     qualityToken =
         iC.esConsumes<SiStripQuality, SiStripQualityRcd>(conf.getParameter<edm::ESInputTag>("siStripQualityLabel"));
   }
 }

Member Function Documentation

void SiStripRecHitConverterAlgorithm::fillBad128StripBlocks	(	const uint32_t	detid,
		bool	bad128StripBlocks[6]
	)		const

private

Definition at line 279 of file SiStripRecHitConverterAlgorithm.cc.

References SiStripQuality::getBadApvs(), SiStripQuality::getBadFibers(), dqmiolumiharvest::j, maskBad128StripBlocks, and quality.

Referenced by run().

                                                                                                                  {
   if (maskBad128StripBlocks) {
     short badApvs = quality->getBadApvs(detid);
     short badFibers = quality->getBadFibers(detid);
     for (int j = 0; j < 6; j++) {
       bad128StripBlocks[j] = (badApvs & (1 << j));
     }
     for (int j = 0; j < 3; j++) {
       if (badFibers & (1 << j)) {
         bad128StripBlocks[2 * j + 0] = true;
         bad128StripBlocks[2 * j + 1] = true;
       }
     }
   }
 }

void SiStripRecHitConverterAlgorithm::fillPSetDescription ( edm::ParameterSetDescription & desc )

static

Definition at line 33 of file SiStripRecHitConverterAlgorithm.cc.

References edm::ParameterSetDescription::add().

Referenced by SiStripRecHitConverter::fillDescriptions().

                                                                                           {
   desc.add<bool>("useSiStripQuality", false);
   desc.add<bool>("MaskBadAPVFibers", false);
   desc.add<bool>("doMatching", true);
   desc.add<edm::ESInputTag>("StripCPE", edm::ESInputTag("StripCPEfromTrackAngleESProducer", "StripCPEfromTrackAngle"));
   desc.add<edm::ESInputTag>("Matcher", edm::ESInputTag("SiStripRecHitMatcherESProducer", "StandardMatcher"));
   desc.add<edm::ESInputTag>("siStripQualityLabel", edm::ESInputTag());
 }

void SiStripRecHitConverterAlgorithm::initialize ( const edm::EventSetup & es )

Definition at line 42 of file SiStripRecHitConverterAlgorithm.cc.

References cpeToken, doMatching, edm::EventSetup::getData(), matcher, matcherToken, parameterestimator, quality, qualityToken, tracker, trackerToken, and useQuality.

Referenced by SiStripRecHitConverter::produce().

                                                                         {
   tracker = &es.getData(trackerToken);
   parameterestimator = &es.getData(cpeToken);
   if (doMatching) {
     matcher = &es.getData(matcherToken);
   }
   if (useQuality) {
     quality = &es.getData(qualityToken);
   }
 }

bool SiStripRecHitConverterAlgorithm::isMasked	(	const SiStripCluster &	cluster,
		bool	bad128StripBlocks[6]
	)		const

inlineprivate

Definition at line 295 of file SiStripRecHitConverterAlgorithm.cc.

References SiStripCluster::amplitudes(), SiStripCluster::barycenter(), SiStripCluster::firstStrip(), maskBad128StripBlocks, and SiStripCluster::size().

Referenced by run().

                                                                                                                     {
   if (maskBad128StripBlocks) {
     if (bad128StripBlocks[cluster.firstStrip() >> 7]) {
       if (bad128StripBlocks[(cluster.firstStrip() + cluster.amplitudes().size()) >> 7] ||
           bad128StripBlocks[static_cast<int32_t>(cluster.barycenter() - 0.499999) >> 7]) {
         return true;
       }
     } else {
       if (bad128StripBlocks[(cluster.firstStrip() + cluster.amplitudes().size()) >> 7] &&
           bad128StripBlocks[static_cast<int32_t>(cluster.barycenter() - 0.499999) >> 7]) {
         return true;
       }
     }
   }
   return false;
 }

void SiStripRecHitConverterAlgorithm::match	(	products &	output,
		LocalVector	trackdirection
	)		const

private

Definition at line 147 of file SiStripRecHitConverterAlgorithm.cc.

References edmNew::DetSet< T >::begin(), edm::OwnVector< T, P >::begin(), edm::OwnVector< T, P >::clear(), filterCSVwithJSON::copy, edmNew::DetSet< T >::detId(), SiStripRecHitMatcher::doubleMatch(), edmNew::DetSet< T >::empty(), edmNew::DetSet< T >::end(), edm::OwnVector< T, P >::end(), TrackerGeometry::idToDet(), SiStripRecHitMatcher::match(), SiStripRecHitConverterAlgorithm::products::matched, matcher, StripSubdetector::partnerDetId(), edmNew::DetSetVector< T >::FastFiller::push_back(), edmNew::DetSetVector< T >::FastFiller::resize(), SiStripRecHitConverterAlgorithm::products::rphi, SiStripRecHitConverterAlgorithm::products::rphiUnmatched, edmNew::DetSet< T >::size(), edm::OwnVector< T, P >::size(), SiStripRecHitConverterAlgorithm::products::stereo, SiStripRecHitConverterAlgorithm::products::stereoUnmatched, and tracker.

Referenced by run().

                                                                                               {
   int nmatch = 0;
   edm::OwnVector<SiStripMatchedRecHit2D> collectorMatched;  // gp/FIXME: avoid this
 
   // Remember the ends of the collections, as we will use them a lot
   SiStripRecHit2DCollection::const_iterator edStereoDet = output.stereo->end();
   SiStripRecHit2DCollection::const_iterator edRPhiDet = output.rphi->end();
 
   // two work vectors for bookeeping clusters used by the stereo part of the matched hits
   std::vector<SiStripRecHit2D::ClusterRef::key_type> matchedSteroClusters;
 
   for (SiStripRecHit2DCollection::const_iterator itRPhiDet = output.rphi->begin(); itRPhiDet != edRPhiDet;
        ++itRPhiDet) {
     edmNew::DetSet<SiStripRecHit2D> rphiHits = *itRPhiDet;
     StripSubdetector specDetId(rphiHits.detId());
     uint32_t partnerId = specDetId.partnerDetId();
 
     // if not part of a glued pair
     if (partnerId == 0) {
       // I must copy these as unmatched
       if (!rphiHits.empty()) {
         SiStripRecHit2DCollection::FastFiller filler(*output.rphiUnmatched, rphiHits.detId());
         filler.resize(rphiHits.size());
         std::copy(rphiHits.begin(), rphiHits.end(), filler.begin());
       }
       continue;
     }
 
     SiStripRecHit2DCollection::const_iterator itStereoDet = output.stereo->find(partnerId);
 
     // if the partner is not found (which probably can happen if it's empty)
     if (itStereoDet == edStereoDet) {
       // I must copy these as unmatched
       if (!rphiHits.empty()) {
         SiStripRecHit2DCollection::FastFiller filler(*output.rphiUnmatched, rphiHits.detId());
         filler.resize(rphiHits.size());
         std::copy(rphiHits.begin(), rphiHits.end(), filler.begin());
       }
       continue;
     }
 
     edmNew::DetSet<SiStripRecHit2D> stereoHits = *itStereoDet;
 
     // Get ready for making glued hits
     const GluedGeomDet* gluedDet = (const GluedGeomDet*)tracker->idToDet(DetId(specDetId.glued()));
     typedef SiStripMatchedRecHit2DCollection::FastFiller Collector;
     Collector collector(*output.matched, specDetId.glued());
 
     // Prepare also the list for unmatched rphi hits
     SiStripRecHit2DCollection::FastFiller fillerRphiUnm(*output.rphiUnmatched, rphiHits.detId());
 
     // a list of clusters used by the matched part of the stereo hits in this detector
     matchedSteroClusters.clear();  // at the beginning, empty
 
 #ifdef DOUBLE_MATCH
     CollectorHelper chelper(collector, collectorMatched, fillerRphiUnm, matchedSteroClusters);
     matcher->doubleMatch(
         rphiHits.begin(), rphiHits.end(), stereoHits.begin(), stereoHits.end(), gluedDet, trackdirection, chelper);
     nmatch += chelper.nmatch;
 #else
     // Make simple collection of this (gp:FIXME: why do we need it?)
     SiStripRecHitMatcher::SimpleHitCollection stereoSimpleHits;
     // gp:FIXME: use std::transform
     stereoSimpleHits.reserve(stereoHits.size());
     for (edmNew::DetSet<SiStripRecHit2D>::const_iterator it = stereoHits.begin(), ed = stereoHits.end(); it != ed;
          ++it) {
       stereoSimpleHits.push_back(&*it);
     }
 
     for (edmNew::DetSet<SiStripRecHit2D>::const_iterator it = rphiHits.begin(), ed = rphiHits.end(); it != ed; ++it) {
       matcher->match(
           &(*it), stereoSimpleHits.begin(), stereoSimpleHits.end(), collectorMatched, gluedDet, trackdirection);
       if (collectorMatched.size() > 0) {
         nmatch += collectorMatched.size();
         for (edm::OwnVector<SiStripMatchedRecHit2D>::const_iterator itm = collectorMatched.begin(),
                                                                     edm = collectorMatched.end();
              itm != edm;
              ++itm) {
           collector.push_back(*itm);
           // mark the stereo hit cluster as used, so that the hit won't go in the unmatched stereo ones
           matchedSteroClusters.push_back(itm->stereoClusterRef().key());
         }
         collectorMatched.clear();
       } else {
         // store a copy of this rphi hit as an unmatched rphi hit
         fillerRphiUnm.push_back(*it);
       }
     }
 
 #endif
 
     // discard matched hits if the collection is empty
     if (collector.empty())
       collector.abort();
 
     // discard unmatched rphi hits if there are none
     if (fillerRphiUnm.empty())
       fillerRphiUnm.abort();
 
     // now look for unmatched stereo hits
     SiStripRecHit2DCollection::FastFiller fillerStereoUnm(*output.stereoUnmatched, stereoHits.detId());
     std::sort(matchedSteroClusters.begin(), matchedSteroClusters.end());
     for (edmNew::DetSet<SiStripRecHit2D>::const_iterator it = stereoHits.begin(), ed = stereoHits.end(); it != ed;
          ++it) {
       if (!std::binary_search(matchedSteroClusters.begin(), matchedSteroClusters.end(), it->cluster().key())) {
         fillerStereoUnm.push_back(*it);
       }
     }
     if (fillerStereoUnm.empty())
       fillerStereoUnm.abort();
   }
 
   for (SiStripRecHit2DCollection::const_iterator itStereoDet = output.stereo->begin(); itStereoDet != edStereoDet;
        ++itStereoDet) {
     edmNew::DetSet<SiStripRecHit2D> stereoHits = *itStereoDet;
     StripSubdetector specDetId(stereoHits.detId());
     uint32_t partnerId = specDetId.partnerDetId();
     if (partnerId == 0)
       continue;
     SiStripRecHit2DCollection::const_iterator itRPhiDet = output.rphi->find(partnerId);
     if (itRPhiDet == edRPhiDet) {
       if (!stereoHits.empty()) {
         SiStripRecHit2DCollection::FastFiller filler(*output.stereoUnmatched, stereoHits.detId());
         filler.resize(stereoHits.size());
         std::copy(stereoHits.begin(), stereoHits.end(), filler.begin());
       }
     }
   }
 
   edm::LogInfo("SiStripRecHitConverter") << "found\n" << nmatch << "  matched RecHits\n";
 }

void SiStripRecHitConverterAlgorithm::run	(	edm::Handle< edmNew::DetSetVector< SiStripCluster > >	input,
		products &	output
	)

Definition at line 53 of file SiStripRecHitConverterAlgorithm.cc.

References input.

Referenced by SiStripRecHitConverter::produce().

                                                                                                                 {
   run(input, output, LocalVector(0., 0., 0.));
 }

void SiStripRecHitConverterAlgorithm::run	(	edm::Handle< edmNew::DetSetVector< SiStripCluster > >	input,
		products &	output,
		LocalVector	trackdirection
	)

Definition at line 57 of file SiStripRecHitConverterAlgorithm.cc.

References edmNew::DetSetVector< T >::FastFiller::abort(), doMatching, relval_steps::DS, edmNew::DetSetVector< T >::FastFiller::empty(), fillBad128StripBlocks(), TrackerGeometry::idToDetUnit(), isMasked(), StripClusterParameterEstimator::localParameters(), match(), parameterestimator, edmNew::DetSetVector< T >::FastFiller::push_back(), SiStripRecHitConverterAlgorithm::products::rphi, SiStripRecHitConverterAlgorithm::products::stereo, StripSubdetector::stereo(), tracker, and useModule().

                                                                       {
   for (auto const& DS : *inputhandle) {
     auto id = DS.id();
     if (!useModule(id))
       continue;
 
     Collector collector = StripSubdetector(id).stereo() ? Collector(*output.stereo, id) : Collector(*output.rphi, id);
 
     bool bad128StripBlocks[6];
     fillBad128StripBlocks(id, bad128StripBlocks);
 
     GeomDetUnit const& du = *(tracker->idToDetUnit(id));
     for (auto const& cluster : DS) {
       if (isMasked(cluster, bad128StripBlocks))
         continue;
 
       StripClusterParameterEstimator::LocalValues parameters = parameterestimator->localParameters(cluster, du);
       collector.push_back(
           SiStripRecHit2D(parameters.first, parameters.second, du, DS.makeRefTo(inputhandle, &cluster)));
     }
 
     if (collector.empty())
       collector.abort();
   }
   if (doMatching) {
     match(output, trackdirection);
   }
 }

bool SiStripRecHitConverterAlgorithm::useModule ( const uint32_t id ) const

inlineprivate

Definition at line 312 of file SiStripRecHitConverterAlgorithm.cc.

References TrackerGeometry::idToDetUnit(), SiStripQuality::IsModuleUsable(), quality, tracker, and useQuality.

Referenced by run().

                                                                               {
   const StripGeomDetUnit* stripdet = (const StripGeomDetUnit*)tracker->idToDetUnit(id);
   if (stripdet == nullptr)
     edm::LogWarning("SiStripRecHitConverter") << "Detid=" << id << " not found";
   return stripdet != nullptr && (!useQuality || quality->IsModuleUsable(id));
 }