SiStripFineDelayTOF Class Reference

#include <SiStripFineDelayTOF.h>

List of all members.

Static Public Member Functions
static double	timeOfFlight (bool cosmics, bool field, double *trackParameters, double *hit, double *phit, bool onDisk)
static void	trackParameters (const reco::Track &tk, double *trackParameters)
Private Member Functions
	SiStripFineDelayTOF ()
virtual	~SiStripFineDelayTOF ()
Static Private Member Functions
static double	getPhi (double *trackParameters, double *hit, bool onDisk)
static double	timeOfFlightBeam (double *hit, double *phit)
static double	timeOfFlightBeamB (double *trackParameters, double *hit, double *phit, bool onDisk)
static double	timeOfFlightCosmic (double *hit, double *phit)
static double	timeOfFlightCosmicB (double *trackParameters, double *hit, double *phit, bool onDisk)
static double	x (double *trackParameters, double phi)
static double	y (double *trackParameters, double phi)
static double	z (double *trackParameters, double phi)

---

Detailed Description

Definition at line 10 of file SiStripFineDelayTOF.h.

---

Constructor & Destructor Documentation

SiStripFineDelayTOF::SiStripFineDelayTOF

(

)

[private]

virtual SiStripFineDelayTOF::~SiStripFineDelayTOF

(

)

[private, virtual]

---

Member Function Documentation

double SiStripFineDelayTOF::getPhi	(	double *	trackParameters,
		double *	hit,
		bool	onDisk
	)		[static, private]

Definition at line 119 of file SiStripFineDelayTOF.cc.

References funct::cos(), phi, Pi, funct::sin(), funct::tan(), x(), and y().

Referenced by timeOfFlightCosmicB().

{
 if(onDisk) //use z coordinate to find phi
   return trackParameters[2] + trackParameters[0]*tan(trackParameters[1])*(hit[2]-trackParameters[4]);
 else // use x,y coordinate to find phi
 {
   double phi =0;
   if(trackParameters[2]>-2.35 && trackParameters[2]<-0.78) {
     // first use y
     phi = acos(-trackParameters[0]*(hit[1]-(trackParameters[3]-1/trackParameters[0])*(-cos(trackParameters[2]))));
     // use x to resolve the ambiguity
     if(fabs(x(trackParameters,phi)-hit[0])>fabs(x(trackParameters,-phi)-hit[0])) phi *= -1.;
   } else {
     // first use x
     phi = asin(trackParameters[0]*(hit[0]-(trackParameters[3]-1/trackParameters[0])*sin(trackParameters[2])));
     // use y to resolve the ambiguity
     if((fabs(y(trackParameters,phi)-hit[1])>fabs(y(trackParameters,TMath::Pi()-phi)-hit[1]))) 
       phi = phi>0 ? TMath::Pi()-phi: -TMath::Pi()-phi;
   }
   return phi;
 }
 return 0.;
}

double SiStripFineDelayTOF::timeOfFlight	(	bool	cosmics,
		bool	field,
		double *	trackParameters,
		double *	hit,
		double *	phit,
		bool	onDisk
	)		[static]

Definition at line 10 of file SiStripFineDelayTOF.cc.

References timeOfFlightBeam(), timeOfFlightBeamB(), timeOfFlightCosmic(), and timeOfFlightCosmicB().

{

  // case 1: cosmics with no field.
  if(cosmics && !field) {
   return timeOfFlightCosmic(hit,phit);
  }
  // case 2: cosmics with field.
  else if(cosmics && field) {
   return timeOfFlightCosmicB(trackParameters,hit,phit,onDisk);
  }
  // case 3: beam with no field.
  else if(!cosmics && !field) {
   return timeOfFlightBeam(hit,phit);
  }
  // case 4: beam with field.
  else {
   return timeOfFlightBeamB(trackParameters,hit,phit,onDisk);
  }
}

double SiStripFineDelayTOF::timeOfFlightBeam	(	double *	hit,
		double *	phit
	)		[static, private]

Definition at line 86 of file SiStripFineDelayTOF.cc.

References trackerHits::c.

Referenced by timeOfFlight().

{
  // constants
  const double c = 30; // cm/ns
  TVector3 r0(hit[0],hit[1],hit[2]);
  return r0.Mag()/c;
}

double SiStripFineDelayTOF::timeOfFlightBeamB	(	double *	trackParameters,
		double *	hit,
		double *	phit,
		bool	onDisk
	)		[static, private]

Definition at line 94 of file SiStripFineDelayTOF.cc.

References trackerHits::c, funct::cos(), funct::sin(), mathSSE::sqrt(), and theta().

Referenced by timeOfFlight().

{
  // constants
  const double c = 30; // cm/ns
  // track parameters
  double& theta = trackParameters[1];
  // returns the time of flight from the origin
  return onDisk? fabs(hit[2]/(c*cos(theta))) : fabs(sqrt(hit[0]*hit[0]+hit[1]*hit[1])/(c*sin(theta)));
}

double SiStripFineDelayTOF::timeOfFlightCosmic	(	double *	hit,
		double *	phit
	)		[static, private]

Definition at line 31 of file SiStripFineDelayTOF.cc.

References trackerHits::c, L1TEmulatorMonitor_cff::p, and mathSSE::sqrt().

Referenced by timeOfFlight().

{
  // constants
  const double c = 30; // cm/ns
#ifndef TIF_COSMIC_SETUP
  const double Rmu = 385; // cm
  const double zmu = 560; // cm
  // estimate the time for crossing the barrel
  TVector3 r0(hit[0],hit[1],hit[2]);
  TVector3 pr(phit[0],phit[1],phit[2]);
  TVector2 r0_xy = r0.XYvector();
  TVector2 pr_xy = pr.Unit().XYvector();
  double t_barrel = ((r0_xy*pr_xy)+sqrt( (r0_xy*pr_xy)*(r0_xy*pr_xy) - r0_xy.Mod2() + Rmu*Rmu  ))/c*pr.Mag()/pr.Perp();
  // estimate the time for crossing endcaps
  double t_endcap = fabs(((phit[2]/fabs(phit[2])*zmu)+hit[2])/c*pr.Mag()/pr.Pz());
  // take the minimum
  return t_barrel<t_endcap ? t_barrel : t_endcap;
#else
  const double y_trigger = 100; //cm
  double p = sqrt(phit[0]*phit[0]+phit[1]*phit[1]+phit[2]*phit[2]);
//  LogDebug("TOF") << "momentum:" << phit[0] << " " << phit[1] << " " << phit[2];
//  LogDebug("TOF") << "p/py=" << p/phit[1] << "  Y0,Y,dY = " << y_trigger << " " << hit[1] << " " << (y_trigger-hit[1]); 
//  LogDebug("TOF") << "d, t=d/c : " << ((y_trigger-hit[1])*(p/phit[1])) << " " << ((y_trigger-hit[1])*(p/phit[1])/c); 
  return fabs((y_trigger-hit[1])*(p/phit[1])/c);
#endif
}

double SiStripFineDelayTOF::timeOfFlightCosmicB	(	double *	trackParameters,
		double *	hit,
		double *	phit,
		bool	onDisk
	)		[static, private]

Definition at line 58 of file SiStripFineDelayTOF.cc.

References trackerHits::c, getPhi(), phi, funct::sin(), funct::tan(), and theta().

Referenced by timeOfFlight().

{
  // constants
  const double Rmu = 385; // cm
  const double zmu = 560; // cm
  const double c = 30; // cm/ns
  // track parameters
  double& kappa = trackParameters[0];
  double& theta = trackParameters[1];
  double& phi_0 = trackParameters[2];
  double& d_0   = trackParameters[3];
  double& z_0   = trackParameters[4];
  double invkappa = 1/kappa;
  // computes the value of the track parameter that correspond to the hit, relative to phi_0
  //double phi = kappa*tan(theta)*(hit[2]-z_0);
  double phi = getPhi(trackParameters,hit,onDisk) - phi_0;
  // computes the value of the track parameter that correspond to the muon system, relative to phi_0
  // phi_mu = phi_mu0 - phi_0
  double phi_mu_b = (kappa>0 ? -1 : 1 ) * acos((Rmu*Rmu-d_0*d_0-2*invkappa*invkappa+2*invkappa*d_0)/(2*invkappa*d_0-2*invkappa*invkappa));
  double phi_mu_e = kappa*tan(theta)*((phit[2]<0 ? 1 : -1)*zmu-z_0);
  // estimate the time for crossing the barrel
  double t_barrel = fabs(invkappa/(c*sin(theta))*(phi-phi_mu_b));
  // estimate the time for crossing endcaps
  double t_endcap = fabs(invkappa/(c*sin(theta))*(phi-phi_mu_e));
  // take the minimum
  return t_barrel<t_endcap ? t_barrel : t_endcap;
}

void SiStripFineDelayTOF::trackParameters	(	const reco::Track &	tk,
		double *	trackParameters
	)		[static]

Definition at line 143 of file SiStripFineDelayTOF.cc.

References reco::Track::innerPosition(), LogDebug, reco::Track::outerMomentum(), reco::Track::outerPosition(), Pi, position, and CosmicsPD_Skims::radius.

Referenced by SiStripFineDelayHit::detId().

{
      math::XYZPoint position = tk.outerPosition();
      math::XYZVector momentum = tk.outerMomentum();
      LogDebug("trackParameters") << "outer position: " << position.x() << " " << position.y() << " " << position.z();
      LogDebug("trackParameters") << "outer momentum: " << momentum.x() << " " << momentum.y() << " " << momentum.z();
      math::XYZVector fieldDirection(0,0,1);
      // computes the center of curvature
      math::XYZVector radius = momentum.Cross(fieldDirection).Unit();
      position -= radius / trackParameters[0];
      LogDebug("trackParameters") << "center of curvature: " << position.x() << " " << position.y() << " " << position.z();
      // the transverse IP
      trackParameters[3] = position.rho() - fabs(1/trackParameters[0]);
      if(trackParameters[0]>0) trackParameters[3] *= -1.;
      // phi_0
      double phi_out = trackParameters[2];
      double phi_0 = position.phi() - TMath::Pi()/2;
      if(trackParameters[0]<0) phi_0 -= TMath::Pi();
      if(phi_0<-2*TMath::Pi()) phi_0 += 2*TMath::Pi();
      if(phi_0>2*TMath::Pi()) phi_0 -= 2*TMath::Pi();
      if(phi_0 > 0) phi_0 -= 2*TMath::Pi();
      LogDebug("trackParameters") << "phi_0: " << phi_0;
      trackParameters[2] = phi_0;
      // z_0
      trackParameters[4] = tk.outerPosition().z() - (phi_out-phi_0)/TMath::Tan(trackParameters[1])/trackParameters[0];
      LogDebug("trackParameters") << "z_0: " << tk.outerPosition().z() << " " << trackParameters[4] << " " << tk.innerPosition().z();  
}

double SiStripFineDelayTOF::x	(	double *	trackParameters,
		double	phi
	)		[static, private]

Definition at line 104 of file SiStripFineDelayTOF.cc.

References funct::sin().

Referenced by getPhi().

{
  return trackParameters[3]*sin(trackParameters[2]) + (1/trackParameters[0])*(sin(phi)-sin(trackParameters[2]));
}

double SiStripFineDelayTOF::y	(	double *	trackParameters,
		double	phi
	)		[static, private]

Definition at line 109 of file SiStripFineDelayTOF.cc.

References funct::cos().

Referenced by getPhi().

{
  return -trackParameters[3]*cos(trackParameters[2]) - (1/trackParameters[0])*(cos(phi)-cos(trackParameters[2]));
}

double SiStripFineDelayTOF::z	(	double *	trackParameters,
		double	phi
	)		[static, private]

Definition at line 114 of file SiStripFineDelayTOF.cc.

References funct::tan().

{
  return trackParameters[4] + (phi-trackParameters[2])/(trackParameters[0]*tan(trackParameters[1]));
}