StripCPE Class Reference

#include <StripCPE.h>

Inheritance diagram for StripCPE:

Classes
struct	AlgoParam
struct	Param

Public Member Functions
LocalVector	driftDirection (const StripGeomDetUnit *det) const override
AlgoParam	getAlgoParam (const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const
virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const
virtual LocalValues	localParameters (const SiStripCluster &, const GeomDetUnit &) const
virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const
virtual void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const override
virtual StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, AlgoParam const &ap) const
	StripCPE (edm::ParameterSet &conf, const MagneticField &, const TrackerGeometry &, const SiStripLorentzAngle &, const SiStripBackPlaneCorrection &, const SiStripConfObject &, const SiStripLatency &)
Public Member Functions inherited from StripClusterParameterEstimator
virtual void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const
virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const
virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const
virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd) const
virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const
virtual	~StripClusterParameterEstimator ()

Protected Member Functions
Param const &	param (const GeomDetUnit &det) const

Protected Attributes
const SiStripBackPlaneCorrection &	BackPlaneCorrectionMap_
const TrackerGeometry &	geom_
const SiStripLorentzAngle &	LorentzAngleMap_
const MagneticField &	magfield_
const bool	peakMode_
std::vector< float >	xtalk1
std::vector< float >	xtalk2

Private Types
typedef std::vector< Param >	Params

Private Member Functions
void	fillParams ()

Private Attributes
unsigned int	m_off
Params	m_Params

Additional Inherited Members
Public Types inherited from StripClusterParameterEstimator
using	AClusters = DynArray< SiStripCluster const * >
using	ALocalValues = DynArray< LocalValues >
using	LocalValues = std::pair< LocalPoint, LocalError >
typedef std::vector< LocalValues >	VLocalValues

Detailed Description

Definition at line 17 of file StripCPE.h.

Member Typedef Documentation

Params

typedef std::vector<Param> StripCPE::Params

private

Definition at line 85 of file StripCPE.h.

Constructor & Destructor Documentation

StripCPE()

StripCPE::StripCPE	(	edm::ParameterSet &	conf,
		const MagneticField &	mag,
		const TrackerGeometry &	geom,
		const SiStripLorentzAngle &	LorentzAngle,
		const SiStripBackPlaneCorrection &	BackPlaneCorrection,
		const SiStripConfObject &	confObj,
		const SiStripLatency &	latency
	)

Definition at line 9 of file StripCPE.cc.

References fillParams(), SiStripConfObject::get(), IB1, IB2, SiStripConfObject::isParameter(), LogMessageMonitor_cff::modules, OB1, OB2, peakMode_, edm::second(), findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, W1A, W1B, W2A, W2B, W3A, W3B, W4, W5, W6, W7, xtalk1, and xtalk2.

    : peakMode_(latency.singleReadOutMode() == 1),
      geom_(geom),
      magfield_(mag),
      LorentzAngleMap_(LorentzAngle),
      BackPlaneCorrectionMap_(BackPlaneCorrection) {
  typedef std::map<std::string, SiStripModuleGeometry> map_t;
  map_t modules;
  modules["IB1"] = SiStripModuleGeometry::IB1;
  modules["IB2"] = SiStripModuleGeometry::IB2;
  modules["OB1"] = SiStripModuleGeometry::OB1;
  modules["OB2"] = SiStripModuleGeometry::OB2;
  modules["W1A"] = SiStripModuleGeometry::W1A;
  modules["W2A"] = SiStripModuleGeometry::W2A;
  modules["W3A"] = SiStripModuleGeometry::W3A;
  modules["W1B"] = SiStripModuleGeometry::W1B;
  modules["W2B"] = SiStripModuleGeometry::W2B;
  modules["W3B"] = SiStripModuleGeometry::W3B;
  modules["W4"] = SiStripModuleGeometry::W4;
  modules["W5"] = SiStripModuleGeometry::W5;
  modules["W6"] = SiStripModuleGeometry::W6;
  modules["W7"] = SiStripModuleGeometry::W7;

  const unsigned size =
      static_cast<unsigned int>(
          max_element(modules.begin(), modules.end(), [&](auto& arg1, auto& arg2) { return arg1.second < arg2.second; })
              ->second) +
      1;
  xtalk1.resize(size);
  xtalk2.resize(size);

  for (map_t::const_iterator it = modules.begin(); it != modules.end(); it++) {
    const std::string modeS(peakMode_ ? "Peak" : "Deco"), xtalk1S("xtalk1_" + it->first + modeS),
        xtalk2S("xtalk2_" + it->first + modeS);

    if (!confObj.isParameter(xtalk1S))
      throw cms::Exception("SiStripConfObject does not contain: ") << xtalk1S;
    if (!confObj.isParameter(xtalk2S))
      throw cms::Exception("SiStripConfObject does not contain: ") << xtalk2S;

    xtalk1[static_cast<int>(it->second)] = confObj.get<double>(xtalk1S);
    xtalk2[static_cast<int>(it->second)] = confObj.get<double>(xtalk2S);
  }

  fillParams();
}

Member Function Documentation

driftDirection()

LocalVector StripCPE::driftDirection ( const StripGeomDetUnit *  det ) const

overridevirtual

Implements StripClusterParameterEstimator.

Definition at line 77 of file StripCPE.cc.

References GeomDet::geographicalId(), SiStripLorentzAngle::getLorentzAngle(), MagneticField::inTesla(), LorentzAngleMap_, magfield_, GloballyPositioned< T >::position(), DetId::rawId(), GeomDet::surface(), toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by fillParams().

                                                                      {
  LocalVector lbfield = (det->surface()).toLocal(magfield_.inTesla(det->surface().position()));

  float tanLorentzAnglePerTesla = LorentzAngleMap_.getLorentzAngle(det->geographicalId().rawId());

  float dir_x = -tanLorentzAnglePerTesla * lbfield.y();
  float dir_y = tanLorentzAnglePerTesla * lbfield.x();
  float dir_z = 1.f;  // E field always in z direction

  return LocalVector(dir_x, dir_y, dir_z);
}

fillParams()

void StripCPE::fillParams ( )

private

Definition at line 89 of file StripCPE.cc.

References TkRadialStripTopology::angularWidth(), cms::cuda::assert(), BackPlaneCorrectionMap_, Surface::bounds(), TrackerGeometry::detUnits(), driftDirection(), f, GeomDet::geographicalId(), geom_, SiStripBackPlaneCorrection::getBackPlaneCorrection(), mps_fire::i, GeomDet::index(), createfilelist::int, GeomDetType::isTrackerStrip(), TkRadialStripTopology::localPitch(), LogDebug, m_off, m_Params, TrackerGeometry::offsetDU(), AlCaHLTBitMon_ParallelJobs::p, funct::pow(), DetId::rawId(), GeomDet::specificSurface(), StripGeomDetType::specificTopology(), StripGeomDetUnit::specificType(), mathSSE::sqrt(), TkRadialStripTopology::stripLength(), GeomDetEnumerators::tkDetEnum, StripGeomDetUnit::topology(), and StripGeomDetUnit::type().

Referenced by StripCPE().

                          {
  auto const& dus = geom_.detUnits();
  m_off = dus.size();
  for (unsigned int i = 1; i < 7; ++i) {
    LogDebug("LookingForFirstStrip")
        << "Subdetector " << i << " GeomDetEnumerator " << GeomDetEnumerators::tkDetEnum[i] << " offset "
        << geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) << " is it strip? "
        << (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size()
                ? dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip()
                : false);
    if (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() &&
        dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip()) {
      if (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) < m_off)
        m_off = geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]);
    }
  }
  LogDebug("LookingForFirstStrip") << " Chosen offset: " << m_off;
  m_Params.resize(dus.size() - m_off);
  for (auto i = m_off; i != dus.size(); ++i) {
    auto& p = m_Params[i - m_off];
    const StripGeomDetUnit* stripdet = (const StripGeomDetUnit*)(dus[i]);
    assert(stripdet->index() == int(i));
    //    assert(stripdet->geographicalId().subdetId()>1); // not pixel..
    assert(stripdet->type().isTrackerStrip());  // not pixel

    const Bounds& bounds = stripdet->specificSurface().bounds();
    p.maxLength = std::sqrt(std::pow(bounds.length(), 2.f) + std::pow(bounds.width(), 2.f));
    p.thickness = bounds.thickness();
    p.invThickness = 1.f / p.thickness;
    p.drift = driftDirection(stripdet) * p.thickness;
    p.topology = (StripTopology*)(&stripdet->topology());
    p.nstrips = p.topology->nstrips();
    p.moduleGeom = SiStripDetId(stripdet->geographicalId()).moduleGeometry();
    p.backplanecorrection = BackPlaneCorrectionMap_.getBackPlaneCorrection(stripdet->geographicalId().rawId());

    const TkRadialStripTopology* rtop =
        dynamic_cast<const TkRadialStripTopology*>(&stripdet->specificType().specificTopology());
    p.pitch_rel_err2 =
        (rtop)
            ? pow(0.5f * rtop->angularWidth() * rtop->stripLength() / rtop->localPitch(LocalPoint(0, 0, 0)), 2.f) / 12.f
            : 0.f;
  }
}

getAlgoParam()

AlgoParam StripCPE::getAlgoParam ( const GeomDetUnit &  det, const LocalTrajectoryParameters &  ltp  ) const

inline

Definition at line 57 of file StripCPE.h.

References funct::abs(), alignCSCRings::corr, LocalTrajectoryParameters::directionNotNormalized(), f, GeomDet::geographicalId(), AlCaHLTBitMon_ParallelJobs::p, param(), LocalTrajectoryParameters::position(), and HLT_2022v14_cff::track.

Referenced by StripCPEfromTrackAngle::localParameters(), TkStripMeasurementDet::recHits(), and TkStripMeasurementDet::simpleRecHits().

                                                                                             {
    StripCPE::Param const& p = param(det);
    SiStripDetId::SubDetector loc = SiStripDetId(det.geographicalId()).subDetector();

    LocalVector track = ltp.directionNotNormalized();
    track *= -p.thickness;

    const float fullProjection = p.coveredStrips(track + p.drift, ltp.position());

    auto const corr =
        -0.5f * (1.f - p.backplanecorrection) * fullProjection + 0.5f * p.coveredStrips(track, ltp.position());

    return AlgoParam{p, ltp, loc, std::abs(fullProjection), corr};
  }

localParameters() [1/6]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 40 of file StripClusterParameterEstimator.h.

                                                                              {
    return localParameters(cluster, gd);
  }

localParameters() [2/6]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 37 of file StripClusterParameterEstimator.h.

                                                                                       {
    return std::make_pair(LocalPoint(), LocalError());
  }

localParameters() [3/6]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 45 of file StripClusterParameterEstimator.h.

                                                                                  {
    return localParameters(cluster, gd, tsos.localParameters());
  }

localParameters() [4/6]

virtual void StripClusterParameterEstimator::localParameters

inline

Definition at line 32 of file StripClusterParameterEstimator.h.

                                                                           {}

localParameters() [5/6]

StripClusterParameterEstimator::LocalValues StripCPE::localParameters ( const SiStripCluster &  cl, const GeomDetUnit &  det  ) const

overridevirtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 62 of file StripCPE.cc.

References SiStripCluster::barycenter(), f, AlCaHLTBitMon_ParallelJobs::p, param(), and digitizers_cfi::strip.

Referenced by TkStripMeasurementDet::hit().

                                                                                                    {
  StripCPE::Param const& p = param(det);
  const float barycenter = cluster.barycenter();
  const float fullProjection =
      p.coveredStrips(p.drift + LocalVector(0, 0, -p.thickness), p.topology->localPosition(barycenter));
  const float strip = barycenter - 0.5f * (1.f - p.backplanecorrection) * fullProjection;

  return std::make_pair(p.topology->localPosition(strip), p.topology->localError(strip, 1.f / 12.f));
}

localParameters() [6/6]

virtual StripClusterParameterEstimator::LocalValues StripCPE::localParameters ( const SiStripCluster &  cl, AlgoParam const &  ap  ) const

inlinevirtual

Reimplemented in StripCPEfromTrackAngle.

Definition at line 52 of file StripCPE.h.

                                                                                                 {
    return std::make_pair(LocalPoint(), LocalError());
  }

param()

Param const& StripCPE::param ( const GeomDetUnit &  det ) const

inlineprotected

Definition at line 81 of file StripCPE.h.

References GeomDet::index(), m_off, and m_Params.

Referenced by getAlgoParam(), StripCPEgeometric::localParameters(), and localParameters().

{ return m_Params[det.index() - m_off]; }

Member Data Documentation

BackPlaneCorrectionMap_

const SiStripBackPlaneCorrection& StripCPE::BackPlaneCorrectionMap_

protected

Definition at line 77 of file StripCPE.h.

Referenced by fillParams().

geom_

const TrackerGeometry& StripCPE::geom_

protected

Definition at line 74 of file StripCPE.h.

Referenced by fillParams().

LorentzAngleMap_

const SiStripLorentzAngle& StripCPE::LorentzAngleMap_

protected

Definition at line 76 of file StripCPE.h.

Referenced by driftDirection().

m_off

unsigned int StripCPE::m_off

private

Definition at line 87 of file StripCPE.h.

Referenced by fillParams(), and param().

m_Params

Params StripCPE::m_Params

private

Definition at line 86 of file StripCPE.h.

Referenced by fillParams(), and param().

magfield_

const MagneticField& StripCPE::magfield_

protected

Definition at line 75 of file StripCPE.h.

Referenced by driftDirection().

peakMode_

const bool StripCPE::peakMode_

protected

Definition at line 73 of file StripCPE.h.

Referenced by StripCPE().

xtalk1

std::vector<float> StripCPE::xtalk1

protected

Definition at line 78 of file StripCPE.h.

Referenced by StripCPEgeometric::localParameters(), and StripCPE().

xtalk2

std::vector<float> StripCPE::xtalk2

protected

Definition at line 79 of file StripCPE.h.

Referenced by StripCPE().