---

StripCPEfromTrackAngle Class Reference

#include <RecoLocalTracker/SiStripRecHitConverter/interface/StripCPEfromTrackAngle.h>

Inheritance diagram for StripCPEfromTrackAngle:

List of all members.

Public Member Functions
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const LocalTrajectoryParameters &ltp) const
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const
	StripCPEfromTrackAngle (edm::ParameterSet &conf, const MagneticField *mag, const TrackerGeometry *geom, const SiStripLorentzAngle *LorentzAngle)

---

Detailed Description

Definition at line 12 of file StripCPEfromTrackAngle.h.

---

Constructor & Destructor Documentation

StripCPEfromTrackAngle::StripCPEfromTrackAngle	(	edm::ParameterSet &	conf,
		const MagneticField *	mag,
		const TrackerGeometry *	geom,
		const SiStripLorentzAngle *	LorentzAngle
	)			[inline]

Definition at line 16 of file StripCPEfromTrackAngle.h.

    :StripCPE(conf,mag, geom, LorentzAngle ){}

---

Member Function Documentation

StripClusterParameterEstimator::LocalValues StripCPEfromTrackAngle::localParameters	(	const SiStripCluster &	cl,
		const LocalTrajectoryParameters &	ltp
	)			const

Definition at line 6 of file StripCPEfromTrackAngle.cc.

References SiStripCluster::amplitudes(), SiStripCluster::barycenter(), StripCPE::Param::drift, funct::exp(), SiStripCluster::geographicalId(), StripTopology::localError(), StripTopology::localPosition(), m1, PV3DBase< T, PVType, FrameType >::mag(), StripCPE::Param::maxLength, Topology::measurementPosition(), min, LocalTrajectoryParameters::momentum(), StripCPE::Param::nstrips, p, p1, p2, StripCPE::param(), LocalTrajectoryParameters::position(), HLT_VtxMuL3::result, StripCPE::Param::thickness, StripCPE::Param::topology, PV3DBase< T, PVType, FrameType >::x(), PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                                      {
  //
  // get the det from the geometry
  //

  StripCPE::Param const & p = param(DetId(cl.geographicalId()));
  

 
  // ionisation length
   
  //float ionLen = std::min( trackDir.mag(), maxLength);

  //  const float par1=38.07;
  // const float par2=0.3184; 
  // const float par3=0.09828; 
  // const float P1 = par1 * p.thickness; 
  // const float P2 = par2; 
  // const float P3 = par3;


  //  drift*=(thickness/2);

  //calculate error form track angle


  LocalPoint  middlepoint = ltp.position();
  LocalVector trackDir = ltp.momentum()/ltp.momentum().mag();
  LocalPoint  position = p.topology->localPosition(cl.barycenter());



  if(trackDir.z()*p.drift.z() > 0.) trackDir *= -1.;


  if(trackDir.z() !=0.) {
    trackDir *= fabs(p.thickness/trackDir.z());
  } else {
    trackDir *= p.maxLength/trackDir.mag();
  }


      
  if(p.drift.z() == 0.) {
    //  if(drift.z() == 0.||cl.amplitudes().size()==1) {
    edm::LogError("StripCPE") <<"No drift towards anodes !!!";
    LocalError eresult = p.topology->localError(cl.barycenter(),1/12.);
    //  LocalPoint  result=LocalPoint(position.x()-drift.x()/2,position.y()-drift.y()/2,0);
    return std::make_pair(position-p.drift*(p.thickness/2),eresult);
  }      




  // covered length along U
      
  LocalVector middleOfProjection = 0.5*(trackDir + p.drift);

  LocalPoint middlePointOnStrips = middlepoint + 0.5*p.drift;

  LocalPoint p1 = LocalPoint(middlePointOnStrips.x() + middleOfProjection.x()
                             ,middlePointOnStrips.y() + middleOfProjection.y());
  LocalPoint p2 = LocalPoint(middlePointOnStrips.x() - middleOfProjection.x()
                             ,middlePointOnStrips.y() - middleOfProjection.y());

  MeasurementPoint m1 = p.topology->measurementPosition(p1);
  MeasurementPoint m2 = p.topology->measurementPosition(p2);
  float u1 = m1.x();
  float u2 = m2.x();
  float uerr2;
  float uProj = std::min( float(fabs( u1 - u2)), float(p.nstrips));
  if((cl.amplitudes().size() - (uProj+2.5)) > 1)uerr2=cl.amplitudes().size() * cl.amplitudes().size() * (1./12.);
  else{
    //   float par1=38.07;
    //   float par2=0.3184; 
    //   float par3=0.09828; 
    //   float P1 = par1 * thickness; 
    //   float P2 = par2; 
    //   float P3 = par3;
    //   float uerr;
    
    //   uerr =(uProj-P1)*(uProj-P1)*(P2-P3)/(P1*P1)+P3;
    
    const float P1=-0.339;
    const float P2=0.90;
    const float P3=0.279;
    
    float uerr=P1*uProj*exp(-uProj*P2)+P3;
    uerr2=uerr*uerr;
  }

  MeasurementError merror=MeasurementError( uerr2, 0., 1./12.);
  LocalPoint result=LocalPoint(position.x()-0.5*p.drift.x(),position.y()-0.5*p.drift.y(),0);
  MeasurementPoint mpoint=p.topology->measurementPosition(result);
  LocalError eresult=p.topology->localError(mpoint,merror);
  return std::make_pair(result,eresult);
}

StripClusterParameterEstimator::LocalValues StripCPEfromTrackAngle::localParameters	(	const SiStripCluster &	cl,
		const GeomDetUnit &	det,
		const LocalTrajectoryParameters &	ltp
	)			const [inline, virtual]

Reimplemented from ClusterParameterEstimator< SiStripCluster >.

Definition at line 25 of file StripCPEfromTrackAngle.h.

                                                                                                           {
    return localParameters(cl,ltp);
  }; 

---

The documentation for this class was generated from the following files:

• RecoLocalTracker/SiStripRecHitConverter/interface/StripCPEfromTrackAngle.h
• RecoLocalTracker/SiStripRecHitConverter/src/StripCPEfromTrackAngle.cc

---