fw Namespace Reference

Functions
double	estimate_field (const reco::Track &track, bool highQuality=false)
double	estimate_field (double vx1, double vy1, double vx2, double vy2, double px, double py)

Function Documentation

double fw::estimate_field	(	const reco::Track &	track,
		bool	highQuality = false
	)

Definition at line 13 of file estimate_field.cc.

References reco::Track::extra(), reco::Track::innerPosition(), edm::Ref< C, T, F >::isAvailable(), reco::Track::outerMomentum(), reco::Track::outerPosition(), funct::pow(), reco::TrackBase::pt(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::vx(), and reco::TrackBase::vy().

Referenced by FWTrackProxyBuilder::build(), FWMuonBuilder::calculateField(), FWPFTrackUtils::setupLegoTrack(), and FWPFTrackUtils::setupTrack().

{
   if ( !track.extra().isAvailable() ) return -1;
   // We have 3 states so we can form 3 pairs to check for curvature
   // In practice we want to avoid estimates over lost of material and
   // in regions with small field. Therefore we should look at the pair
   // based on POCA and closest state to it. Than we will estimate the
   // filed inside the solenoid. 
   double estimate = -1;
   if ( pow(track.outerPosition().x()-track.vx(),2)+pow(track.outerPosition().y()-track.vy(),2) <
    pow(track.innerPosition().x()-track.vx(),2)+pow(track.innerPosition().y()-track.vy(),2) )
     {
       estimate = estimate_field( track.outerPosition().x(), 
                  track.outerPosition().y(),
                  track.vx(),
                  track.vy(),
                  track.px(),
                  track.py() );
     }
   else
     {
       estimate = estimate_field( track.innerPosition().x(), 
                  track.innerPosition().y(),
                  track.vx(),
                  track.vy(),
                  track.px(),
                  track.py() );
     }
   
   if ( highQuality ){
     if ( fabs( track.outerMomentum().rho()/track.pt()-1 )<0.1 && 
      track.innerPosition().rho()<129 &&
      track.outerPosition().rho()<129 ){
       /*printf("outer-poca: \tPOCA Pt: %0.2f, \tinner Rho: %0.1f, \touter Rho: %0.1f, \tchange in pt: %0.1f%%, \testimate: %0.2f\n",
          track.pt(), track.innerPosition().rho(), track.outerPosition().rho(), 
          fabs( track.outerMomentum().rho()/track.pt()-1 )*100, estimate);*/
       return estimate;
     } else {
       return -1;
     }
   } else
     return estimate;
}

double fw::estimate_field	(	double	vx1,
		double	vy1,
		double	vx2,
		double	vy2,
		double	px,
		double	py
	)

Definition at line 58 of file estimate_field.cc.

References mathSSE::sqrt().

{
  math::XYZVector displacement(vx2-vx1, vy2-vy1, 0);
  math::XYZVector transverseMomentum(px, py, 0);
  double cosAlpha = transverseMomentum.Dot(displacement)/transverseMomentum.r()/displacement.r();
  return 200*sqrt(1-cosAlpha*cosAlpha)/0.2998*transverseMomentum.r()/displacement.r();
}