#include <SiStripFineDelayTLA.h>

Public Member Functions
std::vector< std::pair < std::pair< DetId, LocalPoint > ,float > >	findtrackangle (const TrajectorySeed &seed, const reco::Track &theT)

std::vector< std::pair < std::pair< DetId, LocalPoint > ,float > >	findtrackangle (const reco::Track &theT)

std::vector< std::pair < std::pair< DetId, LocalPoint > ,float > >	findtrackangle (const std::vector< Trajectory > &traj)

std::vector< std::pair < std::pair< DetId, LocalPoint > ,float > >	findtrackangle (const Trajectory &traj)

void	init (const edm::Event &e, const edm::EventSetup &c)

	SiStripFineDelayTLA (const edm::ParameterSet &conf)

virtual	~SiStripFineDelayTLA ()

Private Member Functions
double	computeAngleCorr (const LocalVector &v, double pitch, double thickness)

Private Attributes
edm::ParameterSet	conf_

SimpleTrackRefitter *	refitter_

const TrackerGeometry *	tracker

Detailed Description

Definition at line 21 of file SiStripFineDelayTLA.h.

Constructor & Destructor Documentation

SiStripFineDelayTLA::SiStripFineDelayTLA ( const edm::ParameterSet & conf )

explicit

Definition at line 27 of file SiStripFineDelayTLA.cc.

References refitter_.

                                                                     : 
   conf_(conf)
 {
   refitter_ = new SimpleTrackRefitter(conf);
 }

SiStripFineDelayTLA::~SiStripFineDelayTLA ( )

virtual

Definition at line 44 of file SiStripFineDelayTLA.cc.

References refitter_.

 {  
   delete refitter_;
 }  

Member Function Documentation

double SiStripFineDelayTLA::computeAngleCorr	(	const LocalVector &	v,
		double	pitch,
		double	thickness
	)

private

Definition at line 128 of file SiStripFineDelayTLA.cc.

References dttmaxenums::L, LogDebug, PV3DBase< T, PVType, FrameType >::mag(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by findtrackangle().

 {
   double v_xy = sqrt(v.x()*v.x()+v.y()*v.y());
   double L = fabs(thickness*v_xy/v.z());
   double Lmax = fabs(pitch/v.x()*v_xy);
   if(L<Lmax) {
     LogDebug("SiStripFineDelayTLA ") << L << " vs " << Lmax 
         << " Signal contained in strip. Correction is " << v.z()/v.mag();
     return v.z()/v.mag();
   } else {
     LogDebug("SiStripFineDelayTLA ") << L << " vs " << Lmax
        << " Signal not contained in strip. Correction is " << thickness/pitch*v.x()/v_xy*v.z()/v.mag()
        << " instead of " << v.z()/v.mag();
     return thickness/pitch*v.x()/v_xy*v.z()/v.mag();
   }
 }

std::vector< std::pair< std::pair< DetId, LocalPoint >,float > > SiStripFineDelayTLA::findtrackangle	(	const TrajectorySeed &	seed,
		const reco::Track &	theT
	)

Definition at line 55 of file SiStripFineDelayTLA.cc.

References refitter_, and SimpleTrackRefitter::refitTrack().

Referenced by SiStripFineDelayHit::detId(), and findtrackangle().

                                                                                                                                             {
   vector<Trajectory> traj = refitter_->refitTrack(seed,theT);
   return findtrackangle(traj);
 }       

std::vector< std::pair< std::pair< DetId, LocalPoint >,float > > SiStripFineDelayTLA::findtrackangle ( const reco::Track & theT )

Definition at line 49 of file SiStripFineDelayTLA.cc.

References findtrackangle(), refitter_, and SimpleTrackRefitter::refitTrack().

 {
   std::vector<Trajectory> traj = refitter_->refitTrack(theT);
   return findtrackangle(traj);
 }

std::vector< std::pair< std::pair< DetId, LocalPoint >,float > > SiStripFineDelayTLA::findtrackangle ( const std::vector< Trajectory > & traj )

Definition at line 60 of file SiStripFineDelayTLA.cc.

References findtrackangle().

 {
   if (trajVec.size()) {
   return findtrackangle(trajVec.front()); }
   std::vector<std::pair< std::pair<DetId, LocalPoint> ,float> > hitangleassociation;
   return hitangleassociation;
 }

std::vector< std::pair< std::pair< DetId, LocalPoint >,float > > SiStripFineDelayTLA::findtrackangle ( const Trajectory & traj )

Definition at line 68 of file SiStripFineDelayTLA.cc.

References angle(), computeAngleCorr(), TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), TrajectoryStateOnSurface::localDirection(), BaseSiTrackerRecHit2DLocalPos::localPosition(), TrajectoryStateOnSurface::localPosition(), Trajectory::measurements(), GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), GluedGeomDet::stereoDet(), SiStripMatchedRecHit2D::stereoHit(), GeomDet::toGlobal(), GeomDet::toLocal(), GeomDetUnit::topology(), tracker, and PV3DBase< T, PVType, FrameType >::z().

 {
   std::vector<std::pair< std::pair<DetId, LocalPoint> ,float> >hitangleassociation;
   std::vector<TrajectoryMeasurement> TMeas=traj.measurements();
   std::vector<TrajectoryMeasurement>::iterator itm;
   for (itm=TMeas.begin();itm!=TMeas.end();itm++){
     TrajectoryStateOnSurface tsos=itm->updatedState();
     const TransientTrackingRecHit::ConstRecHitPointer thit=itm->recHit();
     const SiStripMatchedRecHit2D* matchedhit=dynamic_cast<const SiStripMatchedRecHit2D*>((*thit).hit());
     const SiStripRecHit2D* hit=dynamic_cast<const SiStripRecHit2D*>((*thit).hit());
     LocalVector trackdirection=tsos.localDirection();
     if(matchedhit){//if matched hit...
     GluedGeomDet * gdet=(GluedGeomDet *)tracker->idToDet(matchedhit->geographicalId());
     GlobalVector gtrkdir=gdet->toGlobal(trackdirection);
     // trackdirection on mono det
     // this the pointer to the mono hit of a matched hit 
     const SiStripRecHit2D *monohit=matchedhit->monoHit();
     const GeomDetUnit * monodet=gdet->monoDet();
     LocalVector monotkdir=monodet->toLocal(gtrkdir);
     if(monotkdir.z()!=0){
       // the local angle (mono)
           float localpitch = ((StripTopology*)(&monodet->topology()))->localPitch(tsos.localPosition());
           float thickness = ((((((monohit->geographicalId())>>25)&0x7f)==0xd)||
                          ((((monohit->geographicalId())>>25)&0x7f)==0xe))&&
                        ((((monohit->geographicalId())>>5)&0x7)>4)) ? 0.0500 : 0.0320;
           float angle = computeAngleCorr(monotkdir, localpitch, thickness);
       hitangleassociation.push_back(make_pair(make_pair(monohit->geographicalId(),monohit->localPosition()), angle)); 
       // trackdirection on stereo det
       // this the pointer to the stereo hit of a matched hit 
       const SiStripRecHit2D *stereohit=matchedhit->stereoHit();
       const GeomDetUnit * stereodet=gdet->stereoDet(); 
       LocalVector stereotkdir=stereodet->toLocal(gtrkdir);
       if(stereotkdir.z()!=0){
         // the local angle (stereo)
             float localpitch = ((StripTopology*)(&stereodet->topology()))->localPitch(tsos.localPosition());
             float thickness = ((((((stereohit->geographicalId())>>25)&0x7f)==0xd)||
                            ((((stereohit->geographicalId())>>25)&0x7f)==0xe))&&
                          ((((stereohit->geographicalId())>>5)&0x7)>4)) ? 0.0500 : 0.0320;
             float angle = computeAngleCorr(stereotkdir, localpitch, thickness);
         hitangleassociation.push_back(make_pair(make_pair(stereohit->geographicalId(),stereohit->localPosition()), angle)); 
       }
     }
     }
     else if(hit){
         GeomDetUnit * det=(GeomDetUnit *)tracker->idToDet(hit->geographicalId());
     //  hit= pointer to the rechit
     if(trackdirection.z()!=0){
           // the local angle (single hit)
           float localpitch = ((StripTopology*)(&det->topology()))->localPitch(tsos.localPosition());
           float thickness = ((((((hit->geographicalId())>>25)&0x7f)==0xd)||
                          ((((hit->geographicalId())>>25)&0x7f)==0xe))&&
                        ((((hit->geographicalId())>>5)&0x7)>4)) ? 0.0500 : 0.0320;
           float angle = computeAngleCorr(trackdirection, localpitch, thickness);
       hitangleassociation.push_back(make_pair(make_pair(hit->geographicalId(),hit->localPosition()), angle)); 
     }
     }
   }
   return hitangleassociation;
 }

void SiStripFineDelayTLA::init	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

Definition at line 33 of file SiStripFineDelayTLA.cc.

References edm::EventSetup::get(), refitter_, SimpleTrackRefitter::setServices(), and tracker.

Referenced by SiStripFineDelayHit::produce().

 {
   // get geometry
   edm::ESHandle<TrackerGeometry> estracker;
   es.get<TrackerDigiGeometryRecord>().get(estracker);
   tracker=&(*estracker);
   // the refitter
   refitter_->setServices(es);
 }