#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 &)

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_

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

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	(	edm::ParameterSet &	conf,
		const MagneticField &	magf,
		const TrackerGeometry &	geom
	)

Definition at line 5 of file Phase2StripCPE.cc.

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

                                                                                                              :
    magfield_(magf),
    geom_(geom),
    tanLorentzAnglePerTesla_(conf.getParameter<double>("TanLorentzAnglePerTesla"))
 {
   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());
   }
   fillParam();
 }

Member Function Documentation

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

Definition at line 35 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 fillParam().

 {
   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);
 }

void Phase2StripCPE::fillParam ( )

private

Definition at line 48 of file Phase2StripCPE.cc.

References Surface::bounds(), constexpr, 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(), 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
    }
 }

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

overridevirtual

Implements ClusterParameterEstimator< Phase2TrackerCluster1D >.

Definition at line 20 of file Phase2StripCPE.cc.

References Phase2TrackerCluster1D::center(), Phase2TrackerCluster1D::column(), f, objects.autophobj::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 );
 }