CSCWireTopology.cc

Go to the documentation of this file.

#include "Geometry/CSCGeometry/interface/CSCWireTopology.h"

#include "Geometry/CSCGeometry/src/CSCGangedWireGrouping.h"
#include "Geometry/CSCGeometry/src/CSCNonslantedWireGeometry.h"
#include "Geometry/CSCGeometry/src/CSCSlantedWireGeometry.h"

#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include "DataFormats/Math/interface/GeantUnits.h"

#include <cmath>

using namespace geant_units;
using namespace geant_units::operators;

CSCWireTopology::~CSCWireTopology() { 
  delete theWireGrouping;
  delete theWireGeometry; 
}

CSCWireTopology::CSCWireTopology( 
             const CSCWireGroupPackage& wg,
             double yOfFirstWire,
                     float wireAngleInDegrees ) :
      theWireGrouping( nullptr ), theWireGeometry( nullptr ),
      theAlignmentPinToFirstWire( wg.alignmentPinToFirstWire / 10. ) {

     // Pass consecutiveGroups and wiresInConsecutiveGroups
     // directly on to theWireGrouping ctor. These vectors
     // are transient info and can go once ctor is complete.
         
  theWireGrouping = new CSCGangedWireGrouping( wg.consecutiveGroups,
                   wg.wiresInEachGroup, wg.numberOfGroups );

  const float zeroprecision = 1.E-06; // blur zero a bit, for comparisons

  float wireAngleInRadians = convertDegToRad( wireAngleInDegrees );

  //@@ Conversion from mm to cm
  float wireSpacing = convertMmToCm( wg.wireSpacing );
  float nw =          convertMmToCm( wg.narrowWidthOfWirePlane );
  float ww =          convertMmToCm( wg.wideWidthOfWirePlane );
  float lw =          convertMmToCm( wg.lengthOfWirePlane );

   LogTrace("CSCWireTopology|CSC") <<  
       "CSCWireTopology constructing CSCWireGeometry with:\n" <<
       " wireSpacing = " << wireSpacing*10. << " (mm) " <<
       ", yOfFirstWire = " << yOfFirstWire << " (cm) " << 
     ", wireAngle = " << wireAngleInDegrees << " (deg) = " << wireAngleInRadians << " (rads)" <<
     ", extent: n, w, l = " << nw << ", " << ww << ", " << lw << " (cm)";

  if ( fabs(wireAngleInDegrees) > zeroprecision ) { 
    theWireGeometry = new CSCSlantedWireGeometry( wireSpacing, yOfFirstWire, nw, ww, lw, wireAngleInRadians );
  }
  else {
    theWireGeometry = new CSCNonslantedWireGeometry( wireSpacing, yOfFirstWire, nw, ww, lw );
  }
}

CSCWireTopology::CSCWireTopology( const CSCWireTopology& mewt ) :
  theAlignmentPinToFirstWire(mewt.theAlignmentPinToFirstWire) {
  if (mewt.theWireGrouping) theWireGrouping = mewt.theWireGrouping->clone();
  else theWireGrouping = nullptr;
  if (mewt.theWireGeometry) theWireGeometry = mewt.theWireGeometry->clone();
  else theWireGeometry = nullptr;
  
}

CSCWireTopology& CSCWireTopology::operator=( const CSCWireTopology& mewt ) {
  if ( &mewt != this ) {
    delete theWireGrouping;
    if ( mewt.theWireGrouping ) 
      theWireGrouping = mewt.theWireGrouping->clone();
    else 
      theWireGrouping = nullptr;
    
    delete theWireGeometry;
    if ( mewt.theWireGeometry ) 
      theWireGeometry = mewt.theWireGeometry->clone();
    else 
      theWireGeometry = nullptr;

    theAlignmentPinToFirstWire = mewt.theAlignmentPinToFirstWire;
    
  }
  return *this;
}

LocalPoint CSCWireTopology::localPosition( const MeasurementPoint& ) const {
  edm::LogWarning("CSC") << "CSCWireTopology: localPosition unimplemented. Don't use it." << "\n";
  return LocalPoint();
}

LocalError CSCWireTopology::localError( const MeasurementPoint&, 
                      const MeasurementError& ) const {
  edm::LogWarning("CSC") << "CSCWireTopology: localError unimplemented. Don't use it." << "\n";
  return LocalError();
}

MeasurementPoint CSCWireTopology::measurementPosition( 
                       const LocalPoint& ) const {
  edm::LogWarning("CSC") << "CSCWireTopology: measurementPosition unimplemented. Don't use it." << "\n";
  return MeasurementPoint();
}

MeasurementError CSCWireTopology::measurementError( const LocalPoint&, 
                       const LocalError& ) const {
  edm::LogWarning("CSC") << "CSCWireTopology: measurementError unimplemented. Don't use it." << "\n";
  return MeasurementError();
}

int CSCWireTopology::channel( const LocalPoint& p ) const {
  int wire = theWireGeometry->nearestWire( p );
  int group = theWireGrouping->wireGroup( wire );
  return group; // 0 means out of range or dead region
}

float CSCWireTopology::yOfWireGroup(int wireGroup, float x) const {
   float wire = middleWireOfGroup( wireGroup );
   return theWireGeometry->yOfWire( wire, x );
}
 
float CSCWireTopology::yResolution( int wireGroup ) const {
  // Return the resolution in the wire group measurement (in cm).
  // To really be resolution in local y this should be scaled by Sin
  // factors involving strip and wire angles, and also add a contribution
  // from strip measurement (which is negligible w.r.t. wire resolution.)
  // These factors can increase the number here by up to 1.26.

  // Should be improvable for smaller partial wire groups

  return wireSpacing() * theWireGrouping->numberOfWiresPerGroup( wireGroup ) / sqrt(12.);
}

std::pair<float, float> CSCWireTopology::equationOfWire( float wire ) const {
  return theWireGeometry->equationOfWire( wire );
}

float CSCWireTopology::restrictToYOfWirePlane( float y ) const {
  // Reset y to lie within bounds of wire plane at top and bottom.
  
  std::pair<float, float> ylim = theWireGeometry->yLimitsOfWirePlane();

  if ( y < ylim.first ) {
    y = ylim.first;
  }
  else if ( y > ylim.second ) {
    y = ylim.second;
  }
  return y;
}

bool CSCWireTopology::insideYOfWirePlane( float y ) const {
  // Returns true if arg falls within y limits of wire plane; false otherwise.
      
  std::pair<float, float> ylim = theWireGeometry->yLimitsOfWirePlane();

  if ( y < ylim.first ) {
    return false;
  }
  else if ( y > ylim.second ) {
    return false;
  }
  return true;
}