#include <Phase2StripCPE.h>

Inheritance diagram for Phase2StripCPE:

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

	Phase2StripCPE ()

	Phase2StripCPE (edm::ParameterSet &conf, const MagneticField &)

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 ()

Protected Attributes
const MagneticField *	magfield_

double	tanLorentzAnglePerTesla_

bool	use_LorentzAngle_DB_

Detailed Description

Definition at line 13 of file Phase2StripCPE.h.

Member Typedef Documentation

using Phase2StripCPE::Phase2TrackerGeomDetUnit = PixelGeomDetUnit

Definition at line 18 of file Phase2StripCPE.h.

using Phase2StripCPE::Phase2TrackerTopology = PixelTopology

Definition at line 19 of file Phase2StripCPE.h.

Constructor & Destructor Documentation

Phase2StripCPE::Phase2StripCPE ( )

inline

Definition at line 23 of file Phase2StripCPE.h.

23 {};

Phase2StripCPE::Phase2StripCPE	(	edm::ParameterSet &	conf,
		const MagneticField &	magf
	)

Definition at line 6 of file Phase2StripCPE.cc.

References Exception, edm::ParameterSet::getParameter(), magfield_, tanLorentzAnglePerTesla_, and use_LorentzAngle_DB_.

 {
   magfield_ = &magf;
   use_LorentzAngle_DB_ = conf.getParameter<bool>("LorentzAngle_DB");
   if (use_LorentzAngle_DB_) {
     throw cms::Exception("Lorentz Angle from DB not implemented yet");
     tanLorentzAnglePerTesla_ = 0;
     // old code: LorentzAngleMap_.getLorentzAngle(det->geographicalId().rawId());
   } else {
     tanLorentzAnglePerTesla_ = conf.getParameter<double>("TanLorentzAnglePerTesla");
   }
 }

Member Function Documentation

LocalVector Phase2StripCPE::driftDirection ( const Phase2TrackerGeomDetUnit & det ) const

Definition at line 47 of file Phase2StripCPE.cc.

References MagneticField::inTesla(), magfield_, GloballyPositioned< T >::position(), GeomDet::surface(), tanLorentzAnglePerTesla_, toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by localParameters().

 {
   LocalVector lbfield = (det.surface()).toLocal(magfield_->inTesla(det.surface().position()));  
 
   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);
 }

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

virtual

Implements ClusterParameterEstimator< Phase2TrackerCluster1D >.

Definition at line 20 of file Phase2StripCPE.cc.

References Surface::bounds(), Phase2TrackerCluster1D::center(), Phase2TrackerCluster1D::column(), shallow::drift(), driftDirection(), PixelTopology::localX(), PixelTopology::localY(), PixelTopology::pitch(), funct::pow(), GeomDet::specificSurface(), PixelGeomDetUnit::specificTopology(), Bounds::thickness(), and PV3DBase< T, PVType, FrameType >::x().

 {
   const Phase2TrackerGeomDetUnit & det = (const Phase2TrackerGeomDetUnit &) detunit;
   const Phase2TrackerTopology * topo = &det.specificTopology();
 
   float pitch_x = topo->pitch().first;
   float pitch_y = topo->pitch().second;
 
   // see https://github.com/cms-sw/cmssw/blob/CMSSW_8_1_X/RecoLocalTracker/SiStripRecHitConverter/src/StripCPE.cc
 
   float thickness = det.specificSurface().bounds().thickness();
   LocalVector drift = driftDirection(det) * thickness;
   LocalVector lvec = drift + LocalVector(0,0,-thickness);
   float coveredStrips = lvec.x() / pitch_x; // simplifies wrt Phase0 tracker because only rectangular modules
 
   float ix = cluster.center() - 0.5 * coveredStrips;
   float iy = cluster.column()+0.5; // halfway the column
 
   LocalPoint lp( topo->localX(ix), topo->localY(iy), 0 );          // x, y, z
   LocalError le( pow(pitch_x, 2) / 12, 0, pow(pitch_y, 2) / 12);   // e2_xx, e2_xy, e2_yy
 
   return std::make_pair( lp, le );
 }