#include <CSCLayerGeometry.h>

Inheritance diagram for CSCLayerGeometry:

Public Member Functions
int	channel (int strip) const

Bounds *	clone () const override

	CSCLayerGeometry (const CSCGeometry *geom, int iChamberType, const TrapezoidalPlaneBounds &bounds, int nstrips, float stripOffset, float stripPhiPitch, float whereStripsMeet, float extentOfStripPlane, float yCentreOfStripPlane, const CSCWireGroupPackage &wg, float wireAngleInDegrees, double yOfFirstWire, float hThickness)

	CSCLayerGeometry (const CSCLayerGeometry &)

bool	inside (const Local2DPoint &) const override

bool	inside (const Local3DPoint &) const override
	Determine if the point is inside the bounds. More...

bool	inside (const Local3DPoint &, const LocalError &, float scale=1.f) const override

LocalPoint	intersectionOfStripAndWire (float s, int w) const

LocalPoint	intersectionOfTwoLines (std::pair< float, float > p1, std::pair< float, float > p2) const

float	lengthOfWireGroup (int wireGroup) const

LocalPoint	localCenterOfWireGroup (int wireGroup) const

LocalError	localError (int strip, float sigmaStrip, float sigmaWire) const

float	middleWireOfGroup (int wireGroup) const

int	nearestStrip (const LocalPoint &lp) const

int	nearestWire (const LocalPoint &lp) const

int	numberOfStrips () const

int	numberOfWireGroups () const

int	numberOfWires () const

int	numberOfWiresPerGroup (int wireGroup) const

CSCLayerGeometry &	operator= (const CSCLayerGeometry &)

std::pair< LocalPoint, float >	possibleRecHitPosition (float s, int w1, int w2) const

void	setTopology (CSCStripTopology *topology)

int	stagger () const

float	strip (const LocalPoint &lp) const

float	stripAngle (int strip) const

float	stripOffset (void) const

float	stripPhiPitch () const

float	stripPitch () const

float	stripPitch (const LocalPoint &lp) const

LocalPoint	stripWireGroupIntersection (int strip, int wireGroup) const

LocalPoint	stripWireIntersection (int strip, float wire) const

const CSCStripTopology *	topology () const

float	wireAngle () const

int	wireGroup (int wire) const

float	wirePitch () const

const CSCWireTopology *	wireTopology () const

float	xOfStrip (int strip, float y=0.) const

std::pair< float, float >	yLimitsOfStripPlane () const

float	yOfWire (float wire, float x=0.) const

float	yOfWireGroup (int wireGroup, float x=0.) const

float	yResolution (int wireGroup=1) const

	~CSCLayerGeometry () override

Public Member Functions inherited from TrapezoidalPlaneBounds
virtual bool	inside (const Local2DPoint &p) const

bool	inside (const Local2DPoint &p, const LocalError &err, float scale) const override
	Determine if a 2D point is inside the bounds, taking error into account. More...

virtual bool	inside (const Local2DPoint &p, const LocalError &err, float scale=1.f) const
	Determine if a 2D point is inside the bounds, taking error into account. More...

virtual bool	inside (const Local2DPoint &p, float tollerance) const
	Determine if a point is inside the bounds, taking a tollerance into account. More...

virtual bool	inside (const Local3DPoint &) const=0
	Determine if the point is inside the bounds. More...

virtual bool	inside (const Local3DPoint &, const LocalError &, float scale=1.f) const=0
	Determine if a point is inside the bounds, taking error into account. More...

float	length () const override

virtual const std::array< const float, 4 >	parameters () const

float	significanceInside (const Local3DPoint &, const LocalError &) const override

float	thickness () const override

	TrapezoidalPlaneBounds (float be, float te, float a, float t)

float	width () const override

float	widthAtHalfLength () const override

virtual int	yAxisOrientation () const

Public Member Functions inherited from Bounds
void	computeSpan (Surface const &plane)

virtual bool	inside (const Local2DPoint &p, float tollerance) const
	Determine if a point is inside the bounds, taking a tollerance into account. More...

std::pair< float, float > const &	phiSpan () const

std::pair< float, float > const &	rSpan () const

std::pair< float, float > const &	zSpan () const

virtual	~Bounds ()

Private Attributes
float	apothem

int	chamberType

float	hBottomEdge

float	hTopEdge

const std::string	myName

CSCStripTopology *	theStripTopology

CSCWireTopology *	theWireTopology

Friends
std::ostream &	operator<< (std::ostream &, const CSCLayerGeometry &)

Detailed Description

Encapsulates the geometrical details of a CSCLayer in a WireTopology for the wires and in a StripTopology for the strips. Note that it does not have the capability of calculating global values, so all values are in local coordinates.

Author: Tim Cox

Definition at line 25 of file CSCLayerGeometry.h.

Constructor & Destructor Documentation

◆ CSCLayerGeometry() [1/2]

CSCLayerGeometry::CSCLayerGeometry	(	const CSCGeometry *	geom,
		int	iChamberType,
		const TrapezoidalPlaneBounds &	bounds,
		int	nstrips,
		float	stripOffset,
		float	stripPhiPitch,
		float	whereStripsMeet,
		float	extentOfStripPlane,
		float	yCentreOfStripPlane,
		const CSCWireGroupPackage &	wg,
		float	wireAngleInDegrees,
		double	yOfFirstWire,
		float	hThickness
	)

Ctor from basic trapezoidal Chamber parameters.

Parameters

geom	The pointer to the actual CSCGeometry we're building.
iChamberType	The index 1-9 for station/ring combination.
TrapezoidalPlaneBounds	describing geometry of face.
nstrips	No. of strips in cathode plane of a Layer.
stripOffset	Alternate strip planes are relatively shifted by +/-0.25 strip widths.
stripPhiPitch	Delta-phi width of strips (they're fan-shaped) in radians
whereStripsMeet	radial distance from projected intersection of strips to centre of strip plane
extentOfStripPlane	height of strip plane (along its long symmetry axis)
yCentreOfStripPlane	local y of symmetry centre of strip plane (before any offset rotation)
wg	CSCWireGroupPackage encapsulating wire group info.
wireAngleInDegrees	angle of wires w.r.t local x axis.
yOfFirstWire	local y coordinate of first (lowest) wire in wire plane - nearest narrow edge.
hThickness	half-thickness of chamber layer in cm (i.e. half the gas gap).

Definition at line 23 of file CSCLayerGeometry.cc.

     : TrapezoidalPlaneBounds(
           bounds.widthAtHalfLength() - bounds.width() / 2., bounds.width() / 2., bounds.length() / 2., hThickness),
       theWireTopology(nullptr),
       theStripTopology(nullptr),
       hBottomEdge(bounds.widthAtHalfLength() - bounds.width() / 2.),
       hTopEdge(bounds.width() / 2.),
       apothem(bounds.length() / 2.),
       myName("CSCLayerGeometry"),
       chamberType(iChamberType) {
   LogTrace("CSCLayerGeometry|CSC") << myName << ": being constructed, this=" << this;
  
   // Ganged strips in ME1A?
   bool gangedME1A = (iChamberType == 1 && geom->gangedStrips());
  
   CSCStripTopology* aStripTopology = new CSCUngangedStripTopology(
       nstrips, stripPhiPitch, extentOfStripPlane, whereStripsMeet, stripOffset, yCentreOfStripPlane);
  
   if (gangedME1A) {
     theStripTopology = new CSCGangedStripTopology(*aStripTopology, 16);
     delete aStripTopology;
   } else {
     theStripTopology = aStripTopology;
   }
  
   if (!geom->realWireGeometry()) {
     // Approximate ORCA_8_8_0 and earlier calculated geometry...
     float wangler = convertDegToRad(wireAngleInDegrees);
     float wireCos = cos(wangler);
     float wireSin = sin(wangler);
     float y2 = apothem * wireCos + hBottomEdge * fabs(wireSin);
     float wireSpacing = convertMmToCm(wg.wireSpacing);
     float wireOffset = -y2 + wireSpacing / 2.;
     yOfFirstWire = wireOffset / wireCos;
   }
  
   theWireTopology = new CSCWireTopology(wg, yOfFirstWire, wireAngleInDegrees);
   LogTrace("CSCLayerGeometry|CSC") << myName << ": constructed: " << *this;
 }

References apothem, angle_units::operators::convertDegToRad(), geant_units::operators::convertMmToCm(), funct::cos(), relativeConstraints::geom, hBottomEdge, LogTrace, myName, funct::sin(), stripOffset(), stripPhiPitch(), theStripTopology, theWireTopology, CSCWireGroupPackage::wireSpacing, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by clone().

◆ CSCLayerGeometry() [2/2]

CSCLayerGeometry::CSCLayerGeometry ( const CSCLayerGeometry & melg )

Definition at line 75 of file CSCLayerGeometry.cc.

     : TrapezoidalPlaneBounds(melg.hBottomEdge, melg.hTopEdge, melg.apothem, 0.5 * melg.thickness()),
       theWireTopology(nullptr),
       theStripTopology(nullptr),
       hBottomEdge(melg.hBottomEdge),
       hTopEdge(melg.hTopEdge),
       apothem(melg.apothem),
       chamberType(melg.chamberType) {
   // CSCStripTopology is abstract, so need clone()
   if (melg.theStripTopology)
     theStripTopology = melg.theStripTopology->clone();
   // CSCWireTopology is concrete, so direct copy
   if (melg.theWireTopology)
     theWireTopology = new CSCWireTopology(*(melg.theWireTopology));
 }

References CSCStripTopology::clone(), theStripTopology, and theWireTopology.

◆ ~CSCLayerGeometry()

CSCLayerGeometry::~CSCLayerGeometry ( )

override

Definition at line 112 of file CSCLayerGeometry.cc.

                                     {
   LogTrace("CSCLayerGeometry|CSC") << myName << ": being destroyed, this=" << this
                                    << "\nDeleting theStripTopology=" << theStripTopology
                                    << " and theWireTopology=" << theWireTopology;
   delete theStripTopology;
   delete theWireTopology;
 }

References LogTrace, myName, theStripTopology, and theWireTopology.

Member Function Documentation

◆ channel()

int CSCLayerGeometry::channel ( int strip ) const

inline

Electronics channel corresponding to a given strip ...sometimes there will be more than one strip OR'ed into a channel

Definition at line 112 of file CSCLayerGeometry.h.

112 { return theStripTopology->channel(strip); }

References OffsetRadialStripTopology::channel(), strip(), and theStripTopology.

Referenced by MuonTruth::analyze(), CSCHitAssociator::associateCSCHitId(), CSCHitAssociator::associateHitId(), and CSCStripElectronicsSim::readoutElement().

◆ clone()

Bounds* CSCLayerGeometry::clone ( void ) const

inlineoverridevirtual

Utility method to handle proper copying of the class

Reimplemented from TrapezoidalPlaneBounds.

Definition at line 287 of file CSCLayerGeometry.h.

287 { return new CSCLayerGeometry(*this); }

References CSCLayerGeometry().

◆ inside() [1/3]

bool CSCLayerGeometry::inside ( const Local2DPoint & lp ) const

overridevirtual

Reimplemented from TrapezoidalPlaneBounds.

Definition at line 277 of file CSCLayerGeometry.cc.

                                                           {
   LocalPoint lp2(lp.x(), lp.y(), 0.);
   return inside(lp2);
 }

References inside(), PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

◆ inside() [2/3]

bool CSCLayerGeometry::inside ( const Local3DPoint & ) const

overridevirtual

Determine if the point is inside the bounds.

Reimplemented from TrapezoidalPlaneBounds.

Definition at line 262 of file CSCLayerGeometry.cc.

                                                           {
   bool result = false;
   const float epsilon = 1.e-06;
   if (fabs(lp.z()) < thickness() / 2.) {  // thickness of TPB is that of gas layer
     std::pair<float, float> ylims = yLimitsOfStripPlane();
     if ((lp.y() > ylims.first) && (lp.y() < ylims.second)) {
       // 'strip' returns float value between 0. and float(Nstrips) and value outside
       // is set to 0. or float(Nstrips)... add a conservative precision of 'epsilon'
       if ((theStripTopology->strip(lp) > epsilon) && (theStripTopology->strip(lp) < (numberOfStrips() - epsilon)))
         result = true;
     }
   }
   return result;
 }

References geometryDiff::epsilon, numberOfStrips(), mps_fire::result, OffsetRadialStripTopology::strip(), theStripTopology, TrapezoidalPlaneBounds::thickness(), PV3DBase< T, PVType, FrameType >::y(), yLimitsOfStripPlane(), and PV3DBase< T, PVType, FrameType >::z().

◆ inside() [3/3]

bool CSCLayerGeometry::inside	(	const Local3DPoint &	lp,
		const LocalError &	le,
		float	scale = `1.f`
	)		const

overridevirtual

Is a supplied LocalPoint inside the strip region?

This is a more reliable fiducial cut for CSCs than the 'Bounds' of the GeomDet(Unit) since those ranges are looser than the sensitive gas region. There are three versions, to parallel those of the TrapezoidalPlaneBounds which a naive user might otherwise employ.

Reimplemented from TrapezoidalPlaneBounds.

Definition at line 282 of file CSCLayerGeometry.cc.

                                                                                              {
   // Effectively consider that the LocalError components extend the area which is acceptable.
   // Form a little box centered on the point, with x, y diameters defined by the errors
   // and require that ALL four corners of the box fall outside the strip region for failure
  
   // Note that LocalError is 2-dim x,y and doesn't supply a z error
   float deltaX = scale * sqrt(le.xx());
   float deltaY = scale * sqrt(le.yy());
  
   LocalPoint lp1(lp.x() - deltaX, lp.y() - deltaY, lp.z());
   LocalPoint lp2(lp.x() - deltaX, lp.y() + deltaY, lp.z());
   LocalPoint lp3(lp.x() + deltaX, lp.y() + deltaY, lp.z());
   LocalPoint lp4(lp.x() + deltaX, lp.y() - deltaY, lp.z());
  
   return (inside(lp1) || inside(lp2) || inside(lp3) || inside(lp4));
 }

References L1EGammaCrystalsEmulatorProducer_cfi::scale, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by CSCWireHitSim::getIonizationClusters(), CSCMake2DRecHit::hitFromStripAndWire(), inside(), CSCMake2DRecHit::isHitInFiducial(), and MuonSimHitProducer::produce().

◆ intersectionOfStripAndWire()

LocalPoint CSCLayerGeometry::intersectionOfStripAndWire	(	float	s,
		int	w
	)		const

Return 2-dim point at which a strip and a wire intersects.

Input arguments: a (float) strip number, and an (int) wire.
Output: LocalPoint which is at their intersection, or at extreme y of wire plane, as appropriate. (If y is adjusted, x is also adjusted to keep it on same strip.)

Definition at line 199 of file CSCLayerGeometry.cc.

                                                                             {
   std::pair<float, float> pw = theWireTopology->equationOfWire(static_cast<float>(w));
   std::pair<float, float> ps = theStripTopology->equationOfStrip(s);
   LocalPoint sw = intersectionOfTwoLines(ps, pw);
  
   // If point falls outside wire plane, at extremes in local y,
   // replace its y by that of appropriate edge of wire plane
   if (!(theWireTopology->insideYOfWirePlane(sw.y()))) {
     float y = theWireTopology->restrictToYOfWirePlane(sw.y());
     // and adjust x to match new y
     float x = sw.x() + (y - sw.y()) * tan(theStripTopology->stripAngle(s));
     sw = LocalPoint(x, y);
   }
  
   return sw;
 }

References CSCStripTopology::equationOfStrip(), CSCWireTopology::equationOfWire(), CSCWireTopology::insideYOfWirePlane(), intersectionOfTwoLines(), CSCWireTopology::restrictToYOfWirePlane(), alignCSCRings::s, OffsetRadialStripTopology::stripAngle(), funct::tan(), theStripTopology, theWireTopology, w, x, PV3DBase< T, PVType, FrameType >::x(), y, and PV3DBase< T, PVType, FrameType >::y().

Referenced by possibleRecHitPosition().

◆ intersectionOfTwoLines()

LocalPoint CSCLayerGeometry::intersectionOfTwoLines	(	std::pair< float, float >	p1,
		std::pair< float, float >	p2
	)		const

Return the point of intersection of two straight lines (in 2-dim).

Input arguments are pair(m1,c1) and pair(m2,c2) where m=slope, c=intercept (y=mx+c).
BEWARE! Do not call with m1 = m2 ! No trapping !

Definition at line 216 of file CSCLayerGeometry.cc.

                                                                                                             {
   // Calculate the point of intersection of two straight lines (in 2-dim)
   // input arguments are pair(m1,c1) and pair(m2,c2) where m=slope, c=intercept (y=mx+c)
   // BEWARE! Do not call with m1 = m2 ! No trapping !
  
   float m1 = p1.first;
   float c1 = p1.second;
   float m2 = p2.first;
   float c2 = p2.second;
   float x = (c2 - c1) / (m1 - m2);
   float y = (m1 * c2 - m2 * c1) / (m1 - m2);
   return LocalPoint(x, y);
 }

References alignmentValidation::c1, p1, p2, x, and y.

Referenced by intersectionOfStripAndWire().

◆ lengthOfWireGroup()

float CSCLayerGeometry::lengthOfWireGroup ( int wireGroup ) const

Length of a wire group (center wire, across chamber face)

Definition at line 179 of file CSCLayerGeometry.cc.

                                                              {
   // Return length of 'wire' in the middle of the wire group
   float w = middleWireOfGroup(wireGroup);
   std::vector<float> store = theWireTopology->wireValues(w);
   return store[2];
 }

References middleWireOfGroup(), theWireTopology, w, wireGroup(), and CSCWireTopology::wireValues().

◆ localCenterOfWireGroup()

LocalPoint CSCLayerGeometry::localCenterOfWireGroup ( int wireGroup ) const

Local coordinates of center of a wire group. \WARNING Used to be centerOfWireGroup in ORCA but that version now returns GlobalPoint.

Definition at line 163 of file CSCLayerGeometry.cc.

                                                                        {
   // It can use CSCWireTopology::yOfWireGroup for y,
   // But x requires mixing with 'extent' of wire plane
  
   // If the wires are NOT tilted, default to simple calculation...
   if (fabs(wireAngle()) < 1.E-6) {
     float y = yOfWireGroup(wireGroup);
     return LocalPoint(0., y);
   } else {
     // w is "wire" at the center of the wire group
     float w = middleWireOfGroup(wireGroup);
     std::vector<float> store = theWireTopology->wireValues(w);
     return LocalPoint(store[0], store[1]);
   }
 }

References middleWireOfGroup(), theWireTopology, w, wireAngle(), wireGroup(), CSCWireTopology::wireValues(), y, and yOfWireGroup().

Referenced by CSCAnodeLCTAnalyzer::getWGEta(), and ValidateGeometry::validateCSCLayerGeometry().

◆ localError()

LocalError CSCLayerGeometry::localError	(	int	strip,
		float	sigmaStrip,
		float	sigmaWire
	)		const

Transform strip and wire errors to local x, y frame. Need to supply (central) strip of the hit. The sigma's are in distance units.

Definition at line 230 of file CSCLayerGeometry.cc.

                                                                                           {
   // Input sigmas are expected to be in _distance units_
   // - uncertainty in strip measurement (typically from Gatti fit, value is in local x units)
   // - uncertainty in wire measurement (along direction perpendicular to wires)
  
   float wangle = this->wireAngle();
   float strangle = this->stripAngle(strip);
  
   float sinAngdif = sin(strangle - wangle);
   float sinAngdif2 = sinAngdif * sinAngdif;
  
   float du = sigmaStrip / sin(strangle);  // sigmaStrip is just x-component of strip error
   float dv = sigmaWire;
  
   // The notation is
   // wsins = wire resol  * sin(strip angle)
   // wcoss = wire resol  * cos(strip angle)
   // ssinw = strip resol * sin(wire angle)
   // scosw = strip resol * cos(wire angle)
  
   float wsins = dv * sin(strangle);
   float wcoss = dv * cos(strangle);
   float ssinw = du * sin(wangle);
   float scosw = du * cos(wangle);
  
   float dx2 = (scosw * scosw + wcoss * wcoss) / sinAngdif2;
   float dy2 = (ssinw * ssinw + wsins * wsins) / sinAngdif2;
   float dxy = (scosw * ssinw + wcoss * wsins) / sinAngdif2;
  
   return LocalError(dx2, dxy, dy2);
 }

References funct::cos(), PVValHelper::dxy, funct::sin(), stripAngle(), and wireAngle().

Referenced by CSCMake2DRecHit::hitFromStripAndWire().

◆ middleWireOfGroup()

float CSCLayerGeometry::middleWireOfGroup ( int wireGroup ) const

inline

Middle of wire-group. This is the central wire no. for a group with an odd no. of wires. This is a pseudo-wire no. for a group with an even no. of wires. Accordingly, it is non-integer.

Definition at line 191 of file CSCLayerGeometry.h.

191 { return theWireTopology->middleWireOfGroup(wireGroup); }

References CSCWireTopology::middleWireOfGroup(), theWireTopology, and wireGroup().

Referenced by CSCTriggerPrimitivesReader::getGlobalPosition(), CSCMake2DRecHit::hitFromStripAndWire(), lengthOfWireGroup(), localCenterOfWireGroup(), and stripWireGroupIntersection().

◆ nearestStrip()

int CSCLayerGeometry::nearestStrip ( const LocalPoint & lp ) const

inline

Strip nearest a given local point

Definition at line 96 of file CSCLayerGeometry.h.

96 { return theStripTopology->nearestStrip(lp); }

References CSCRadialStripTopology::nearestStrip(), and theStripTopology.

Referenced by MuonCSCChamberResidual::addResidual(), CSCPairResidualsConstraint::calculatePhi(), and CSCEfficiency::recHitSegment_Efficiencies().

◆ nearestWire()

int CSCLayerGeometry::nearestWire ( const LocalPoint & lp ) const

inline

Wire nearest a given local point

Definition at line 101 of file CSCLayerGeometry.h.

101 { return theWireTopology->nearestWire(lp); }

References CSCWireTopology::nearestWire(), and theWireTopology.

Referenced by CSCAnodeLCTAnalyzer::nearestWG().

◆ numberOfStrips()

int CSCLayerGeometry::numberOfStrips ( ) const

inline

How many strips in layer

Definition at line 66 of file CSCLayerGeometry.h.

66 { return theStripTopology->nstrips(); }

References CSCRadialStripTopology::nstrips(), and theStripTopology.

Referenced by CSCRecoBadChannelsAnalyzer::analyze(), CSCSectorReceiverLUT::calcGlobalPhiME(), CSCSectorReceiverLUT::getGlobalEtaValue(), CSCStripElectronicsSim::initParameters(), inside(), operator<<(), CSCCathodeLCTProcessor::run(), CSCGeometryValidate::validateCSCLayerGeometry(), and ValidateGeometry::validateCSCLayerGeometry().

◆ numberOfWireGroups()

int CSCLayerGeometry::numberOfWireGroups ( ) const

inline

How many wire groups in Layer

Definition at line 76 of file CSCLayerGeometry.h.

76 { return theWireTopology->numberOfWireGroups(); }

References CSCWireTopology::numberOfWireGroups(), and theWireTopology.

Referenced by CSCSectorReceiverLUT::getGlobalEtaValue(), CSCWireElectronicsSim::initParameters(), operator<<(), CSCAnodeLCTProcessor::run(), CSCWireElectronicsSim::timeOfFlightCalibration(), and ValidateGeometry::validateCSCLayerGeometry().

◆ numberOfWires()

int CSCLayerGeometry::numberOfWires ( ) const

inline

How many wires in layer

Definition at line 71 of file CSCLayerGeometry.h.

71 { return theWireTopology->numberOfWires(); }

References CSCWireTopology::numberOfWires(), and theWireTopology.

Referenced by operator<<(), and CSCGeometryValidate::validateCSCLayerGeometry().

◆ numberOfWiresPerGroup()

int CSCLayerGeometry::numberOfWiresPerGroup ( int wireGroup ) const

inline

How many wires in a wiregroup

Definition at line 81 of file CSCLayerGeometry.h.

81 { return theWireTopology->numberOfWiresPerGroup(wireGroup); }

References CSCWireTopology::numberOfWiresPerGroup(), theWireTopology, and wireGroup().

Referenced by CSCMake2DRecHit::hitFromStripAndWire().

◆ operator=()

CSCLayerGeometry & CSCLayerGeometry::operator= ( const CSCLayerGeometry & melg )

Definition at line 91 of file CSCLayerGeometry.cc.

                                                                           {
   if (&melg != this) {
     delete theStripTopology;
     if (melg.theStripTopology)
       theStripTopology = melg.theStripTopology->clone();
     else
       theStripTopology = nullptr;
  
     delete theWireTopology;
     if (melg.theWireTopology)
       theWireTopology = new CSCWireTopology(*(melg.theWireTopology));
     else
       theWireTopology = nullptr;
  
     hBottomEdge = melg.hBottomEdge;
     hTopEdge = melg.hTopEdge;
     apothem = melg.apothem;
   }
   return *this;
 }

References apothem, CSCStripTopology::clone(), hBottomEdge, hTopEdge, theStripTopology, and theWireTopology.

◆ possibleRecHitPosition()

std::pair< LocalPoint, float > CSCLayerGeometry::possibleRecHitPosition	(	float	s,
		int	w1,
		int	w2
	)		const

Return estimate of the 2-dim point of intersection of a strip and a cluster of wires.

Input arguments: a (float) strip number, and the wires which delimit a cluster of wires. The wires are expected to be real wire numbers, and not wire-group numbers.

Returned: pair, with members:
first: LocalPoint which is midway along "the" strip between the wire limits,
or the chamber edges, as appropriate. <bf> second: length of the strip between the wires (or edges as appropriate).

Definition at line 186 of file CSCLayerGeometry.cc.

                                                                                                  {
   LocalPoint sw1 = intersectionOfStripAndWire(s, w1);
   LocalPoint sw2 = intersectionOfStripAndWire(s, w2);
  
   // Average the two points
   LocalPoint midpt((sw1.x() + sw2.x()) / 2., (sw1.y() + sw2.y()) / 2);
  
   // Length of strip crossing this group of wires
   float length = sqrt((sw1.x() - sw2.x()) * (sw1.x() - sw2.x()) + (sw1.y() - sw2.y()) * (sw1.y() - sw2.y()));
  
   return std::pair<LocalPoint, float>(midpt, length);
 }

References intersectionOfStripAndWire(), TrapezoidalPlaneBounds::length(), alignCSCRings::s, mathSSE::sqrt(), w2, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by CSCMake2DRecHit::hitFromStripAndWire().

◆ setTopology()

void CSCLayerGeometry::setTopology ( CSCStripTopology * topology )

This class takes ownership of the pointer, and will destroy it

Definition at line 299 of file CSCLayerGeometry.cc.

                                                                 {
   delete theStripTopology;
   theStripTopology = newTopology;
 }

References theStripTopology.

◆ stagger()

int CSCLayerGeometry::stagger ( ) const

inline

Return +1 or -1 for a stripOffset of +0.25 or -0.25 respectively. Requested by trigger people.

Definition at line 124 of file CSCLayerGeometry.h.

124 { return static_cast<int>(4.1 * stripOffset()); }

References stripOffset().

Referenced by CSCCathodeLCTAnalyzer::preselectDigis(), and CSCCathodeLCTProcessor::run().

◆ strip()

float CSCLayerGeometry::strip ( const LocalPoint & lp ) const

inline

Strip in which a given LocalPoint lies. This is a float which represents the fractional strip position within the detector.
Returns zero if the LocalPoint falls at the extreme low edge of the detector or BELOW, and float(nstrips) if it falls at the extreme high edge or ABOVE.

Definition at line 183 of file CSCLayerGeometry.h.

183 { return theStripTopology->strip(lp); }

References OffsetRadialStripTopology::strip(), and theStripTopology.

Referenced by channel(), CSCRecHit2DValidation::plotResolution(), stripAngle(), stripWireGroupIntersection(), stripWireIntersection(), and xOfStrip().

◆ stripAngle()

float CSCLayerGeometry::stripAngle ( int strip ) const

The angle (in radians) of a strip wrt local x-axis.

Definition at line 147 of file CSCLayerGeometry.cc.

                                                   {
   // Cleverly subtly change meaning of stripAngle once more.
   // In TrapezoidalStripTopology it is angle measured
   // counter-clockwise from y axis.
   // In APTST and RST it is angle measured
   // clockwise from y axis.
   // Output of this function is measured counter-clockwise
   // from x-axis, so it is a conventional 2-dim azimuthal angle
   // in the (x,y) local coordinates
  
   // We want angle at centre of strip (strip N covers
   // *float* range N-1 to N-epsilon)
  
   return 0.5_pi - theStripTopology->stripAngle(strip - 0.5);
 }

References strip(), OffsetRadialStripTopology::stripAngle(), and theStripTopology.

Referenced by MuonCSCChamberResidual::addResidual(), CSCPairResidualsConstraint::calculatePhi(), CSCMake2DRecHit::hitFromStripAndWire(), localError(), CSCEfficiency::recHitSegment_Efficiencies(), and stripWireIntersection().

◆ stripOffset()

float CSCLayerGeometry::stripOffset ( void ) const

inline

Offset of strips from symmetrical distribution about local y axis as a fraction of a strip (0 default, but usually +0.25 or -0.25)

Definition at line 118 of file CSCLayerGeometry.h.

118 { return theStripTopology->stripOffset(); }

References OffsetRadialStripTopology::stripOffset(), and theStripTopology.

Referenced by CSCLayerGeometry(), stagger(), and CSCGeometryValidate::validateCSCLayerGeometry().

◆ stripPhiPitch()

float CSCLayerGeometry::stripPhiPitch ( ) const

inline

The phi width of the strips (radians)

Definition at line 157 of file CSCLayerGeometry.h.

157 { return theStripTopology->phiPitch(); }

References CSCRadialStripTopology::phiPitch(), and theStripTopology.

Referenced by CSCLayerGeometry().

◆ stripPitch() [1/2]

float CSCLayerGeometry::stripPitch ( ) const

inline

The width of the strips (in middle)

Definition at line 162 of file CSCLayerGeometry.h.

                            {
     //    return theStripTopology->pitch(); }
     return stripPitch(LocalPoint(0., 0.));
   }

Referenced by CSCMake2DRecHit::hitFromStripAndWire(), and operator<<().

◆ stripPitch() [2/2]

float CSCLayerGeometry::stripPitch ( const LocalPoint & lp ) const

inline

The width of the strip at a given local point

Definition at line 170 of file CSCLayerGeometry.h.

170 { return theStripTopology->localPitch(lp); }

References CSCRadialStripTopology::localPitch(), and theStripTopology.

◆ stripWireGroupIntersection()

LocalPoint CSCLayerGeometry::stripWireGroupIntersection	(	int	strip,
		int	wireGroup
	)		const

Local point at which strip and centre of wire group intersect

Definition at line 141 of file CSCLayerGeometry.cc.

                                                                                       {
   // middleWire is only an actual wire for a group with an odd no. of wires
   float middleWire = middleWireOfGroup(wireGroup);
   return stripWireIntersection(strip, middleWire);
 }

References middleWireOfGroup(), strip(), stripWireIntersection(), and wireGroup().

Referenced by CSCSectorReceiverLUT::getGlobalEtaValue().

◆ stripWireIntersection()

LocalPoint CSCLayerGeometry::stripWireIntersection	(	int	strip,
		float	wire
	)		const

Local point at which strip and wire intersect

Definition at line 120 of file CSCLayerGeometry.cc.

                                                                               {
   // This allows _float_ wire no. so that we can calculate the
   // intersection of a strip with the mid point of a wire group
   // containing an even no. of wires (which is not an actual wire),
   // as well as for a group containing an odd no. of wires.
  
   // Equation of wire and strip as straight lines in local xy
   // y = mx + c where m = tan(angle w.r.t. x axis)
   // At the intersection x = -(cs-cw)/(ms-mw)
   // At y=0, 0 = ms * xOfStrip(strip) + cs => cs = -ms*xOfStrip
   // At x=0, yOfWire(wire) = 0 + cw => cw = yOfWire
  
   float ms = tan(stripAngle(strip));
   float mw = tan(wireAngle());
   float xs = xOfStrip(strip);
   float xi = (ms * xs + yOfWire(wire)) / (ms - mw);
   float yi = ms * (xi - xs);
  
   return LocalPoint(xi, yi);
 }

References strip(), stripAngle(), funct::tan(), wireAngle(), hybridSuperClusters_cfi::xi, xOfStrip(), and yOfWire().

Referenced by stripWireGroupIntersection().

◆ topology()

const CSCStripTopology* CSCLayerGeometry::topology ( ) const

inline

'The' Topology (i.e. Strip Topology) owned by this MELG

Definition at line 272 of file CSCLayerGeometry.h.

272 { return theStripTopology; }

References theStripTopology.

Referenced by FWRecoGeometryESProducer::addCSCGeometry(), CSCCathodeLCTAnalyzer::getStripPhi(), CSCCathodeLCTAnalyzer::nearestHS(), CSCLayer::topology(), CSCChamberSpecs::topology(), CSCGeometryValidate::validateCSCLayerGeometry(), and ValidateGeometry::validateCSCLayerGeometry().

◆ wireAngle()

float CSCLayerGeometry::wireAngle ( ) const

inline

The angle (in radians) of (any) wire wrt local x-axis.

Definition at line 134 of file CSCLayerGeometry.h.

134 { return theWireTopology->wireAngle(); }

References theWireTopology, and CSCWireTopology::wireAngle().

Referenced by localCenterOfWireGroup(), localError(), operator<<(), and stripWireIntersection().

◆ wireGroup()

int CSCLayerGeometry::wireGroup ( int wire ) const

inline

Wire group containing a given wire

Definition at line 106 of file CSCLayerGeometry.h.

106 { return theWireTopology->wireGroup(wire); }

References theWireTopology, and CSCWireTopology::wireGroup().

Referenced by lengthOfWireGroup(), localCenterOfWireGroup(), middleWireOfGroup(), CSCAnodeLCTAnalyzer::nearestWG(), numberOfWiresPerGroup(), CSCWireElectronicsSim::readoutElement(), stripWireGroupIntersection(), yOfWireGroup(), and yResolution().

◆ wirePitch()

float CSCLayerGeometry::wirePitch ( ) const

inline

The distance (in cm) between anode wires

Definition at line 139 of file CSCLayerGeometry.h.

139 { return theWireTopology->wirePitch(); }

References theWireTopology, and CSCWireTopology::wirePitch().

Referenced by operator<<(), and CSCChamberSpecs::wireSpacing().

◆ wireTopology()

const CSCWireTopology* CSCLayerGeometry::wireTopology ( ) const

inline

The Wire Topology owned by this MELG

Definition at line 282 of file CSCLayerGeometry.h.

282 { return theWireTopology; }

References theWireTopology.

Referenced by FWRecoGeometryESProducer::addCSCGeometry(), TrackDetectorAssociator::getTAMuonChamberMatches(), CSCGeometryValidate::validateCSCLayerGeometry(), and ValidateGeometry::validateCSCLayerGeometry().

◆ xOfStrip()

float CSCLayerGeometry::xOfStrip	(	int	strip,
		float	y = `0.`
	)		const

inline

The local x-position of the center of the strip.

Definition at line 175 of file CSCLayerGeometry.h.

175 { return theStripTopology->xOfStrip(strip, y); }

References strip(), theStripTopology, CSCRadialStripTopology::xOfStrip(), and y.

Referenced by CSCLayer::centerOfStrip(), MuonShowerDigiFiller::fill(), CSCStripDigiValidation::plotResolution(), stripWireIntersection(), and ValidateGeometry::validateCSCLayerGeometry().

◆ yLimitsOfStripPlane()

std::pair<float, float> CSCLayerGeometry::yLimitsOfStripPlane ( ) const

inline

Local y limits of the strip plane

Definition at line 218 of file CSCLayerGeometry.h.

218 { return theStripTopology->yLimitsOfStripPlane(); }

References theStripTopology, and CSCStripTopology::yLimitsOfStripPlane().

Referenced by inside().

◆ yOfWire()

float CSCLayerGeometry::yOfWire	(	float	wire,
		float	x = `0.`
	)		const

inline

Local y of a given wire 'number' (float) at given x

Definition at line 196 of file CSCLayerGeometry.h.

196 { return theWireTopology->yOfWire(wire, x); }

References theWireTopology, x, and CSCWireTopology::yOfWire().

Referenced by stripWireIntersection().

◆ yOfWireGroup()

float CSCLayerGeometry::yOfWireGroup	(	int	wireGroup,
		float	x = `0.`
	)		const

inline

Local y of a given wire group at given x

Definition at line 201 of file CSCLayerGeometry.h.

201 { return theWireTopology->yOfWireGroup(wireGroup, x); }

References theWireTopology, wireGroup(), x, and CSCWireTopology::yOfWireGroup().

Referenced by CSCEfficiency::fillWG_info(), localCenterOfWireGroup(), and CSCWireDigiValidation::plotResolution().

◆ yResolution()

float CSCLayerGeometry::yResolution ( int wireGroup = 1 ) const

inline

The measurement resolution from wire groups (in cm.) This approximates the measurement resolution in the local y direction but may be too small by a factor of up to 1.26 due to stripAngle contributions which are neglected here. The last wiregroup may have more wires than others. The other wiregroups, including the first, are the same. One day the wiregroups will be matched to the hardware by using the DDD.

Definition at line 152 of file CSCLayerGeometry.h.

152 { return theWireTopology->yResolution(wireGroup); }

References theWireTopology, wireGroup(), and CSCWireTopology::yResolution().

Referenced by CSCMake2DRecHit::hitFromStripAndWire().

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<<	(	std::ostream &	stream,
		const CSCLayerGeometry &	lg
	)

friend

Output operator for members of class.

Definition at line 304 of file CSCLayerGeometry.cc.

                                                                        {
   stream << "LayerGeometry " << std::endl
          << "------------- " << std::endl
          << "numberOfStrips               " << lg.numberOfStrips() << std::endl
          << "numberOfWires                " << lg.numberOfWires() << std::endl
          << "numberOfWireGroups           " << lg.numberOfWireGroups() << std::endl
          << "wireAngle  (rad)             " << lg.wireAngle()
          << std::endl
          //         << "wireAngle  (deg)      " << lg.theWireAngle << std::endl
          //         << "sin(wireAngle)        " << lg.theWireSin << std::endl
          //         << "cos(wireAngle)        " << lg.theWireCos << std::endl
          << "wirePitch                    " << lg.wirePitch() << std::endl
          << "stripPitch                   " << lg.stripPitch()
          << std::endl
          //         << "numberOfWiresPerGroup " << lg.theNumberOfWiresPerGroup << std::endl
          //         << "numberOfWiresInLastGroup " << lg.theNumberOfWiresInLastGroup << std::endl
          //         << "wireOffset            " << lg.theWireOffset << std::endl
          //         << "whereStripsMeet       " << lg.whereStripsMeet << std::endl;
          << "hBottomEdge                  " << lg.hBottomEdge << std::endl
          << "hTopEdge                     " << lg.hTopEdge << std::endl
          << "apothem                      " << lg.apothem << std::endl
          << "length (should be 2xapothem) " << lg.length() << std::endl
          << "thickness                    " << lg.thickness() << std::endl;
   return stream;
 }

Member Data Documentation

◆ apothem

float CSCLayerGeometry::apothem

private

Definition at line 312 of file CSCLayerGeometry.h.

Referenced by CSCLayerGeometry(), operator<<(), and operator=().

◆ chamberType

int CSCLayerGeometry::chamberType

private

Definition at line 315 of file CSCLayerGeometry.h.

◆ hBottomEdge

float CSCLayerGeometry::hBottomEdge

private

Definition at line 310 of file CSCLayerGeometry.h.

Referenced by CSCLayerGeometry(), operator<<(), and operator=().

◆ hTopEdge

float CSCLayerGeometry::hTopEdge

private

Definition at line 311 of file CSCLayerGeometry.h.

Referenced by operator<<(), and operator=().

◆ myName

const std::string CSCLayerGeometry::myName

private

Definition at line 314 of file CSCLayerGeometry.h.

Referenced by CSCLayerGeometry(), and ~CSCLayerGeometry().

◆ theStripTopology

CSCStripTopology* CSCLayerGeometry::theStripTopology

private

Definition at line 305 of file CSCLayerGeometry.h.

Referenced by channel(), CSCLayerGeometry(), inside(), intersectionOfStripAndWire(), nearestStrip(), numberOfStrips(), operator=(), setTopology(), strip(), stripAngle(), stripOffset(), stripPhiPitch(), stripPitch(), topology(), xOfStrip(), yLimitsOfStripPlane(), and ~CSCLayerGeometry().

◆ theWireTopology

CSCWireTopology* CSCLayerGeometry::theWireTopology

private

Definition at line 299 of file CSCLayerGeometry.h.

Referenced by CSCLayerGeometry(), intersectionOfStripAndWire(), lengthOfWireGroup(), localCenterOfWireGroup(), middleWireOfGroup(), nearestWire(), numberOfWireGroups(), numberOfWires(), numberOfWiresPerGroup(), operator=(), wireAngle(), wireGroup(), wirePitch(), wireTopology(), yOfWire(), yOfWireGroup(), yResolution(), and ~CSCLayerGeometry().

Public Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ CSCLayerGeometry() [1/2]

◆ CSCLayerGeometry() [2/2]

◆ ~CSCLayerGeometry()

Member Function Documentation

◆ channel()

◆ clone()

◆ inside() [1/3]

◆ inside() [2/3]

◆ inside() [3/3]

◆ intersectionOfStripAndWire()

◆ intersectionOfTwoLines()

◆ lengthOfWireGroup()

◆ localCenterOfWireGroup()

◆ localError()

◆ middleWireOfGroup()

◆ nearestStrip()

◆ nearestWire()

◆ numberOfStrips()

◆ numberOfWireGroups()

◆ numberOfWires()

◆ numberOfWiresPerGroup()

◆ operator=()

◆ possibleRecHitPosition()

◆ setTopology()

◆ stagger()

◆ strip()

◆ stripAngle()

◆ stripOffset()

◆ stripPhiPitch()

◆ stripPitch() [1/2]

◆ stripPitch() [2/2]

◆ stripWireGroupIntersection()

◆ stripWireIntersection()

◆ topology()

◆ wireAngle()

◆ wireGroup()

◆ wirePitch()

◆ wireTopology()

◆ xOfStrip()

◆ yLimitsOfStripPlane()

◆ yOfWire()

◆ yOfWireGroup()

◆ yResolution()

Friends And Related Function Documentation

◆ operator<<

Member Data Documentation

◆ apothem

◆ chamberType

◆ hBottomEdge

◆ hTopEdge

◆ myName

◆ theStripTopology

◆ theWireTopology