CSCRecHitDBuilder Class Reference

#include <CSCRecHitDBuilder.h>

Public Member Functions
void	build (const CSCStripDigiCollection *stripds, const CSCWireDigiCollection *wireds, CSCRecHit2DCollection &oc)
	CSCRecHitDBuilder (const edm::ParameterSet &ps)
const CSCLayer *	getLayer (const CSCDetId &detId)
void	setConditions (CSCRecoConditions *reco)
void	setGeometry (const CSCGeometry *geom)
	~CSCRecHitDBuilder ()

Private Attributes
const CSCGeometry *	geom_
CSCHitFromStripOnly *	hitsFromStripOnly_
CSCHitFromWireOnly *	hitsFromWireOnly_
CSCMake2DRecHit *	make2DHits_
bool	makePseudo2DHits
CSCRecoConditions *	recoConditions_
int	stripWireDeltaT
bool	useCalib

Detailed Description

Algorithm to build 2-D RecHit from wire and strip digis in endcap muon CSCs by implementing a 'build' function required by CSCRecHitDProducer.

The builder goes through many stages before building 2-D hits: 1) It finds wire clusters and form wire hits which it stores in CSCWireHit. 2) From these wire hits, it builds pseudo-wire segments to clean up the wire hit collection from noisy hits. Only the hits falling on the segment or far away from existing segments are retained. 1) It then finds strip cluster and hits which it stores in CSCStripHit. 2) Similary to the wire hits, segments are build using the strip hits. Because of the trapezoidal geometry of the strips, all strip hits falling close to the pseudo-strip segments are retained.

Author

Stoyan Stoynev - NU

Definition at line 38 of file CSCRecHitDBuilder.h.

Constructor & Destructor Documentation

CSCRecHitDBuilder()

CSCRecHitDBuilder::CSCRecHitDBuilder ( const edm::ParameterSet &  ps )

explicit

Configure the algorithm via ctor. Receives ParameterSet percolated down from EDProducer which owns this Builder.

Definition at line 31 of file CSCRecHitDBuilder.cc.

                                                              : geom_(nullptr) {
  // Receives ParameterSet percolated down from EDProducer
 
  useCalib = ps.getParameter<bool>("CSCUseCalibrations");
  stripWireDeltaT = ps.getParameter<int>("CSCstripWireDeltaTime");
 
  hitsFromStripOnly_ = new CSCHitFromStripOnly(ps);
  hitsFromWireOnly_ = new CSCHitFromWireOnly(ps);
  make2DHits_ = new CSCMake2DRecHit(ps);
}

References edm::ParameterSet::getParameter(), hitsFromStripOnly_, hitsFromWireOnly_, make2DHits_, stripWireDeltaT, and useCalib.

~CSCRecHitDBuilder()

CSCRecHitDBuilder::~CSCRecHitDBuilder

(

)

Definition at line 42 of file CSCRecHitDBuilder.cc.

                                      {
  delete hitsFromStripOnly_;
  delete hitsFromWireOnly_;
  delete make2DHits_;
}

References hitsFromStripOnly_, hitsFromWireOnly_, and make2DHits_.

Member Function Documentation

build()

void CSCRecHitDBuilder::build	(	const CSCStripDigiCollection *	stripds,
		const CSCWireDigiCollection *	wireds,
		CSCRecHit2DCollection &	oc
	)

Find digis in each CSCLayer, build strip and wire proto-hits in each layer from which pseudo-segments are build to select hits. Then, strip/wire hits are combined to form 2-D hits, whereas remaining "good" strip and wire only hits are also stored
into output collection.

Definition at line 48 of file CSCRecHitDBuilder.cc.

                                                         {
  LogTrace("CSCRecHitDBuilder") << "[CSCRecHitDBuilder] build entered";
 
  if (!geom_)
    throw cms::Exception("MissingGeometry") << "[CSCRecHitDBuilder::getLayer] Missing geometry" << std::endl;
 
  //  create 2-D hits by looking at superposition of strip and wire hit in a layer
  //
  // N.B.  I've sorted the hits from layer 1-6 always, so can test if there are "holes",
  // that is layers without hits for a given chamber.
 
  int layer_idx = 0;
  CSCDetId old_id;
 
  for (CSCStripDigiCollection::DigiRangeIterator it = stripdc->begin(); it != stripdc->end(); ++it) {
    const CSCDetId& id = (*it).first;
    const CSCLayer* layer = getLayer(id);
    const CSCStripDigiCollection::Range& rstripd = (*it).second;
 
    // Skip if no strip digis in this layer
    if (rstripd.second == rstripd.first)
      continue;
 
    const CSCDetId& sDetId = id;
 
    // This is used to test for gaps in layers and needs to be initialized here
    if (layer_idx == 0) {
      old_id = sDetId;
    }
 
    CSCDetId compId = sDetId;
    CSCWireDigiCollection::Range rwired = wiredc->get(sDetId);
    // Skip if no wire digis in this layer
    // But for ME11, real wire digis are labelled as belonging to ME1b, so that's where ME1a must look
    // (We try ME1a - above - anyway, because simulated wire digis are labelled as ME1a.)
    if (rwired.second == rwired.first) {
      if (sDetId.station() != 1 || sDetId.ring() != 4) {
        continue;  // not ME1a, skip to next layer
      }
      // So if ME1a has no wire digis (always the case for data) make the
      // wire digi ID point to ME1b. This is what is compared to the
      // strip digi ID below (and not used anywhere else).
      // Later, rechits use the strip digi ID for construction.
 
      // It is ME1a but no wire digis there, so try ME1b...
      int endcap = sDetId.endcap();
      int chamber = sDetId.chamber();
      int layer = sDetId.layer();
      CSCDetId idw(endcap, 1, 1, chamber, layer);  // Set idw to same layer in ME1b
      compId = idw;
      rwired = wiredc->get(compId);
    }
 
    // Fill bad channel bitsets for this layer
    recoConditions_->fillBadChannelWords(id);
 
    // Build strip hits for this layer
    std::vector<CSCStripHit> const& cscStripHit = hitsFromStripOnly_->runStrip(id, layer, rstripd);
 
    if (cscStripHit.empty())
      continue;
 
    // now build collection of wire only hits !
    std::vector<CSCWireHit> const& cscWireHit = hitsFromWireOnly_->runWire(compId, layer, rwired);
 
    // Build 2D hit for all possible strip-wire pairs
    // overlapping within this layer
 
    LogTrace("CSCRecHitBuilder") << "[CSCRecHitDBuilder] found " << cscStripHit.size() << " strip and "
                                 << cscWireHit.size() << " wire hits in layer " << sDetId;
 
    // Vector to store rechit within layer
    std::vector<CSCRecHit2D> hitsInLayer;
    unsigned int hits_in_layer = 0;
 
    for (auto const& s_hit : cscStripHit) {
      for (auto const& w_hit : cscWireHit) {
        CSCRecHit2D rechit = make2DHits_->hitFromStripAndWire(sDetId, layer, w_hit, s_hit);
        // Store rechit as a Local Point:
        LocalPoint rhitlocal = rechit.localPosition();
        float yreco = rhitlocal.y();
        bool isInFiducial = false;
        //in me1/1 chambers the strip cut region is at local y = 30 cm, +-5 cm area around it proved to be a suitabla region for omiting the check
        if ((sDetId.station() == 1) && (sDetId.ring() == 1 || sDetId.ring() == 4) && (fabs(yreco + 30.) < 5.)) {
          isInFiducial = true;
        } else {
          isInFiducial = make2DHits_->isHitInFiducial(layer, rechit);
        }
        if (isInFiducial) {
          hitsInLayer.push_back(rechit);
          hits_in_layer++;
        }
      }
    }
 
    LogTrace("CSCRecHitDBuilder") << "[CSCRecHitDBuilder] " << hits_in_layer << " rechits found in layer " << sDetId;
 
    // output vector of 2D rechits to collection
    if (hits_in_layer > 0) {
      oc.put(sDetId, hitsInLayer.begin(), hitsInLayer.end());
    }
    layer_idx++;
    old_id = sDetId;
  }
 
  LogTrace("CSCRecHitDBuilder") << "[CSCRecHitDBuilder] " << oc.size() << " 2d rechits created in this event.";
}

References relativeConstraints::chamber, CSCDetId::chamber(), CSCDetId::endcap(), makeMuonMisalignmentScenario::endcap, Exception, CSCRecoConditions::fillBadChannelWords(), geom_, getLayer(), CSCMake2DRecHit::hitFromStripAndWire(), hitsFromStripOnly_, hitsFromWireOnly_, triggerObjects_cff::id, CSCMake2DRecHit::isHitInFiducial(), CSCDetId::layer(), phase1PixelTopology::layer, CSCRecHit2D::localPosition(), LogTrace, make2DHits_, recoConditions_, CSCDetId::ring(), CSCHitFromStripOnly::runStrip(), CSCHitFromWireOnly::runWire(), CSCDetId::station(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by CSCRecHitDProducer::produce().

getLayer()

const CSCLayer * CSCRecHitDBuilder::getLayer

(

const CSCDetId &

detId

)

Definition at line 158 of file CSCRecHitDBuilder.cc.

{ return geom_->layer(detId); }

References geom_, and CSCGeometry::layer().

Referenced by build().

setConditions()

void CSCRecHitDBuilder::setConditions

(

CSCRecoConditions *

reco

)

Pass conditions downstream

Definition at line 160 of file CSCRecHitDBuilder.cc.

                                                             {
  recoConditions_ = reco;
  hitsFromStripOnly_->setConditions(reco);
  hitsFromWireOnly_->setConditions(reco);
  make2DHits_->setConditions(reco);
}

References hitsFromStripOnly_, hitsFromWireOnly_, make2DHits_, dt_dqm_sourceclient_common_cff::reco, recoConditions_, CSCHitFromStripOnly::setConditions(), CSCMake2DRecHit::setConditions(), and CSCHitFromWireOnly::setConditions().

Referenced by CSCRecHitDProducer::CSCRecHitDProducer().

setGeometry()

void CSCRecHitDBuilder::setGeometry ( const CSCGeometry *  geom )

inline

Cache pointer to geometry so it can be passed downstream

Definition at line 61 of file CSCRecHitDBuilder.h.

{ geom_ = geom; }

References relativeConstraints::geom, and geom_.

Referenced by CSCRecHitDProducer::produce().

Member Data Documentation

geom_

const CSCGeometry* CSCRecHitDBuilder::geom_

private

Definition at line 88 of file CSCRecHitDBuilder.h.

Referenced by build(), getLayer(), and setGeometry().

hitsFromStripOnly_

CSCHitFromStripOnly* CSCRecHitDBuilder::hitsFromStripOnly_

private

The Program first constructs proto wire/strip hits which it stores in a special collection. Proto strip/wire segments are then build from these hits and allow to clean up up the list of hits.

Definition at line 80 of file CSCRecHitDBuilder.h.

Referenced by build(), CSCRecHitDBuilder(), setConditions(), and ~CSCRecHitDBuilder().

hitsFromWireOnly_

CSCHitFromWireOnly* CSCRecHitDBuilder::hitsFromWireOnly_

private

Definition at line 81 of file CSCRecHitDBuilder.h.

Referenced by build(), CSCRecHitDBuilder(), setConditions(), and ~CSCRecHitDBuilder().

make2DHits_

CSCMake2DRecHit* CSCRecHitDBuilder::make2DHits_

private

Definition at line 83 of file CSCRecHitDBuilder.h.

Referenced by build(), CSCRecHitDBuilder(), setConditions(), and ~CSCRecHitDBuilder().

makePseudo2DHits

bool CSCRecHitDBuilder::makePseudo2DHits

private

Definition at line 73 of file CSCRecHitDBuilder.h.

recoConditions_

CSCRecoConditions* CSCRecHitDBuilder::recoConditions_

private

Definition at line 92 of file CSCRecHitDBuilder.h.

Referenced by build(), and setConditions().

stripWireDeltaT

int CSCRecHitDBuilder::stripWireDeltaT

private

Definition at line 72 of file CSCRecHitDBuilder.h.

Referenced by CSCRecHitDBuilder().

useCalib

bool CSCRecHitDBuilder::useCalib

private

Definition at line 71 of file CSCRecHitDBuilder.h.

Referenced by CSCRecHitDBuilder().