Phase2StripCPE Class Reference

#include <Phase2StripCPE.h>

Inheritance diagram for Phase2StripCPE:

Classes
struct	Param

Public Types
using	Phase2TrackerGeomDetUnit = PixelGeomDetUnit
using	Phase2TrackerTopology = PixelTopology
Public Types inherited from ClusterParameterEstimator< Phase2TrackerCluster1D >
typedef std::pair< LocalPoint, LocalError >	LocalValues
typedef std::vector< LocalValues >	VLocalValues

Public Member Functions
LocalVector	driftDirection (const Phase2TrackerGeomDetUnit &det) const
LocalValues	localParameters (const Phase2TrackerCluster1D &cluster, const GeomDetUnit &det) const override
	Phase2StripCPE (edm::ParameterSet &conf, const MagneticField &, const TrackerGeometry &, const SiPhase2OuterTrackerLorentzAngle &)
Public Member Functions inherited from ClusterParameterEstimator< Phase2TrackerCluster1D >
virtual void	clearParameters () const
virtual void	enterLocalParameters (unsigned int id, std::pair< int, int > &row_col, LocalValues pos_err_info) const
virtual void	enterLocalParameters (uint32_t id, uint16_t firstStrip, LocalValues pos_err_info) const
virtual LocalValues	localParameters (const Phase2TrackerCluster1D &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const
virtual LocalValues	localParameters (const Phase2TrackerCluster1D &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const
virtual VLocalValues	localParametersV (const Phase2TrackerCluster1D &cluster, const GeomDetUnit &gd) const
virtual VLocalValues	localParametersV (const Phase2TrackerCluster1D &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const
virtual	~ClusterParameterEstimator ()

Private Member Functions
void	fillParam ()

Private Attributes
const TrackerGeometry &	geom_
const SiPhase2OuterTrackerLorentzAngle &	lorentzAngleMap_
unsigned int	m_off
std::vector< Param >	m_Params
const MagneticField &	magfield_
float	tanLorentzAnglePerTesla_
bool	use_LorentzAngle_DB_

Detailed Description

Definition at line 14 of file Phase2StripCPE.h.

Member Typedef Documentation

Phase2TrackerGeomDetUnit

using Phase2StripCPE::Phase2TrackerGeomDetUnit = PixelGeomDetUnit

Definition at line 17 of file Phase2StripCPE.h.

Phase2TrackerTopology

using Phase2StripCPE::Phase2TrackerTopology = PixelTopology

Definition at line 18 of file Phase2StripCPE.h.

Constructor & Destructor Documentation

Phase2StripCPE()

Phase2StripCPE::Phase2StripCPE	(	edm::ParameterSet &	conf,
		const MagneticField &	magf,
		const TrackerGeometry &	geom,
		const SiPhase2OuterTrackerLorentzAngle &	LorentzAngle
	)

Definition at line 5 of file Phase2StripCPE.cc.

References fillParam(), edm::ParameterSet::getParameter(), and use_LorentzAngle_DB_.

    : magfield_(magf),
      geom_(geom),
      lorentzAngleMap_(LorentzAngle),
      tanLorentzAnglePerTesla_(conf.getParameter<double>("TanLorentzAnglePerTesla")) {
  use_LorentzAngle_DB_ = conf.getParameter<bool>("LorentzAngle_DB");
  fillParam();
}

Member Function Documentation

driftDirection()

LocalVector Phase2StripCPE::driftDirection

(

const Phase2TrackerGeomDetUnit &

det

)

const

Definition at line 29 of file Phase2StripCPE.cc.

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

Referenced by fillParam().

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

  float langle =
      use_LorentzAngle_DB_ ? lorentzAngleMap_.getLorentzAngle(det.geographicalId().rawId()) : tanLorentzAnglePerTesla_;

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

  return LocalVector(dir_x, dir_y, dir_z);
}

fillParam()

void Phase2StripCPE::fillParam ( )

private

Definition at line 42 of file Phase2StripCPE.cc.

References cms::cuda::assert(), Surface::bounds(), TrackerGeometry::detUnits(), shallow::drift(), driftDirection(), geom_, mps_fire::i, GeomDet::index(), createfilelist::int, LogDebug, m_off, m_Params, TrackerGeometry::offsetDU(), AlCaHLTBitMon_ParallelJobs::p, GeomDet::specificSurface(), PixelGeomDetUnit::specificTopology(), Bounds::thickness(), Calorimetry_cff::thickness, and GeomDetEnumerators::tkDetEnum.

Referenced by Phase2StripCPE().

                               {
  // in phase 2 they are all pixel topologies...
  auto const& dus = geom_.detUnits();
  m_off = dus.size();
  // skip Barrel and Foward pixels...
  for (unsigned int i = 3; i < 7; ++i) {
    LogDebug("LookingForFirstPhase2OT") << " Subdetector " << i << " GeomDetEnumerator "
                                        << GeomDetEnumerators::tkDetEnum[i] << " offset "
                                        << geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) << std::endl;
    if (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size()) {
      if (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) < m_off)
        m_off = geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]);
    }
  }
  LogDebug("LookingForFirstPhase2OT") << " Chosen offset: " << m_off;

  m_Params.resize(dus.size() - m_off);
  // very very minimal, for sure it will need to expand...
  for (auto i = m_off; i != dus.size(); ++i) {
    auto& p = m_Params[i - m_off];

    const Phase2TrackerGeomDetUnit& det = (const Phase2TrackerGeomDetUnit&)(*dus[i]);
    assert(det.index() == int(i));
    p.topology = &det.specificTopology();

    auto pitch_x = p.topology->pitch().first;
    auto pitch_y = p.topology->pitch().second;

    // see https://github.com/cms-sw/cmssw/blob/CMSSW_8_1_X/RecoLocalTracker/SiStripRecHitConverter/src/StripCPE.cc
    auto thickness = det.specificSurface().bounds().thickness();
    auto drift = driftDirection(det) * thickness;
    auto lvec = drift + LocalVector(0, 0, -thickness);
    p.coveredStrips = lvec.x() / pitch_x;  // simplifies wrt Phase0 tracker because only rectangular modules

    constexpr float o12 = 1. / 12;
    p.localErr = LocalError(o12 * pitch_x * pitch_x, 0, o12 * pitch_y * pitch_y);  // e2_xx, e2_xy, e2_yy
  }
}

localParameters()

Phase2StripCPE::LocalValues Phase2StripCPE::localParameters ( const Phase2TrackerCluster1D &  cluster, const GeomDetUnit &  det  ) const

overridevirtual

Implements ClusterParameterEstimator< Phase2TrackerCluster1D >.

Definition at line 17 of file Phase2StripCPE.cc.

References Phase2TrackerCluster1D::center(), Phase2TrackerCluster1D::column(), f, dqmMemoryStats::float, GeomDet::index(), m_off, m_Params, and AlCaHLTBitMon_ParallelJobs::p.

                                                                                              {
  auto const& p = m_Params[detunit.index() - m_off];
  auto const& topo = *p.topology;
  float ix = cluster.center() - 0.5f * p.coveredStrips;
  float iy = float(cluster.column()) + 0.5f;  // halfway the column

  LocalPoint lp(topo.localX(ix), topo.localY(iy), 0);  // x, y, z

  return std::make_pair(lp, p.localErr);
}

Member Data Documentation

geom_

const TrackerGeometry& Phase2StripCPE::geom_

private

Definition at line 40 of file Phase2StripCPE.h.

Referenced by fillParam().

lorentzAngleMap_

const SiPhase2OuterTrackerLorentzAngle& Phase2StripCPE::lorentzAngleMap_

private

Definition at line 41 of file Phase2StripCPE.h.

Referenced by driftDirection().

m_off

unsigned int Phase2StripCPE::m_off

private

Definition at line 44 of file Phase2StripCPE.h.

Referenced by fillParam(), and localParameters().

m_Params

std::vector<Param> Phase2StripCPE::m_Params

private

Definition at line 37 of file Phase2StripCPE.h.

Referenced by fillParam(), and localParameters().

magfield_

const MagneticField& Phase2StripCPE::magfield_

private

Definition at line 39 of file Phase2StripCPE.h.

Referenced by driftDirection().

tanLorentzAnglePerTesla_

float Phase2StripCPE::tanLorentzAnglePerTesla_

private

Definition at line 43 of file Phase2StripCPE.h.

Referenced by driftDirection().

use_LorentzAngle_DB_

bool Phase2StripCPE::use_LorentzAngle_DB_

private

Definition at line 46 of file Phase2StripCPE.h.

Referenced by driftDirection(), and Phase2StripCPE().