#include <DTTime2DriftParametrization.h>

Classes
struct	drift_distance
	Structure used to return output values. More...

Public Member Functions
bool	computeDriftDistance_mean (double time, double alpha, double by, double bz, short interpolate, drift_distance *dx) const

bool	computeDriftDistance_mode (double time, double alpha, double by, double bz, short interpolate, drift_distance *dx) const

	DTTime2DriftParametrization ()
	Constructor. More...

virtual	~DTTime2DriftParametrization ()
	Destructor. More...

Private Member Functions
unsigned short	MB_DT_Check_boundaries (double distime, double alpha, double by, double bz, short ifl) const

double	MB_DT_delta_x (double, double *) const

double	MB_DT_dist (double, double *) const

unsigned short	MB_DT_drift_distance (double time, double alpha, double by, double bz, drift_distance *DX, short interpolate) const
	Calculate drift distance and spread. More...

void	MB_DT_Get_grid_points (double alpha, double by, double bz, unsigned short p_alpha, unsigned short p_By, unsigned short p_Bz, unsigned short q_alpha, unsigned short q_By, unsigned short q_Bz) const

void	MB_DT_Get_grid_values (double Var, unsigned short pi, unsigned short pj, short Initial, unsigned short N, const double *Values) const

double	MB_DT_MLInterpolation (double al, double by, double bz, double f) const

double	MB_DT_sigma_t_m (double, double *) const

double	MB_DT_sigma_t_p (double, double *) const

Static Private Attributes
static const double	alpha_value [N_alpha]

static const double	By_value [N_By]

static const double	Bz_value [N_Bz]

static const double	fun_sigma_t [N_alpha][N_By][N_Bz][N_Sigma_t]

static const double	fun_x [N_alpha][N_By][N_Bz][N_Par_x]

static const double	velocitySkew

Detailed Description

This class access the cell parametrization by Pablo Garcia-Abia and Jesus Puerta (CIEMAT). The two public methods act as an interface for the parametrization converting the reference frame and the units accordingly to CMSSW convention.

Author: Pablo Garcia-Abia and Jesus Puerta (CIEMAT, Madrid)

Definition at line 21 of file DTTime2DriftParametrization.h.

Constructor & Destructor Documentation

DTTime2DriftParametrization::DTTime2DriftParametrization ( )

Constructor.

Definition at line 18 of file DTTime2DriftParametrization.cc.

18 {}

DTTime2DriftParametrization::~DTTime2DriftParametrization ( )

virtual

Destructor.

Definition at line 20 of file DTTime2DriftParametrization.cc.

20 {}

Member Function Documentation

bool DTTime2DriftParametrization::computeDriftDistance_mean	(	double	time,
		double	alpha,
		double	by,
		double	bz,
		short	interpolate,
		drift_distance *	dx
	)		const

Compute the drift distance and relative errors (cm). The drift distance is computed as the MEAN value of the distance distribution for this particular time. This method takes care of the conversions of units and reference frame used by CMSSW and the parametrization. time is in ns, alpha (=atan(x/-z)) in rad, by and bz are the magnetic field in the layer RF. False is returned if the parametrization fails

Definition at line 61 of file DTTime2DriftParametrization.cc.

References DTTime2DriftParametrization::drift_distance::delta_x, archive::flag, max(), MB_DT_drift_distance(), Geom::pi(), mathSSE::sqrt(), DTTime2DriftParametrization::drift_distance::v_drift, DTTime2DriftParametrization::drift_distance::x_drift, DTTime2DriftParametrization::drift_distance::x_width_m, and DTTime2DriftParametrization::drift_distance::x_width_p.

                                                           {
   // NOTE: This method takes care of convertion of the units in order to use the CMSSW conventions
 
   // The convention used by the parametrization for the reference frame are different
   // with respect to CMSSW
   // X_par = X_CMSSW; Y_par=Z_CMSSW; Z_par = -Y_CMSSW 
   float By_par = bz;  // Bnorm
   float Bz_par = -by; // Bwire
  
   // Express alpha in degrees
   float alpha_par = alpha * 180./Geom::pi();
 
   //--------------------------------------------------------------------
   // Calculate the drift distance and the resolution from the parametrization
   unsigned short flag =
     MB_DT_drift_distance(time, alpha_par, By_par, Bz_par, dx, interpolate);
 
   if(flag!=1)
     return false;
 
   // Convert from mm (used by the parametrization) to cm (used by CMSSW)
   dx->v_drift = dx->v_drift/10.;
   dx->x_drift = dx->x_drift/10.;
   dx->delta_x = dx->delta_x/10.;
   dx->x_width_m = dx->x_width_m/10.;
   dx->x_width_p = dx->x_width_p/10.;
 
   // Correct drift time for the difference between mode and mean
   dx->x_drift = std::max(0.,dx->x_drift - (dx->x_width_p-dx->x_width_m)*sqrt(2./Geom::pi()));
 
   return true;
 
  }

bool DTTime2DriftParametrization::computeDriftDistance_mode	(	double	time,
		double	alpha,
		double	by,
		double	bz,
		short	interpolate,
		drift_distance *	dx
	)		const

Compute the drift distance and relative errors (cm). The drift distance is computed as the MODE value of the distance distribution for this particular time. This method takes care of the conversions of units and reference frame used by CMSSW and the parametrization. time is in ns, alpha (=atan(x/-z)) in rad, by and bz are the magnetic field in the layer RF. False is returned if the parametrization fails

Definition at line 24 of file DTTime2DriftParametrization.cc.

References DTTime2DriftParametrization::drift_distance::delta_x, archive::flag, MB_DT_drift_distance(), Geom::pi(), DTTime2DriftParametrization::drift_distance::v_drift, DTTime2DriftParametrization::drift_distance::x_drift, DTTime2DriftParametrization::drift_distance::x_width_m, and DTTime2DriftParametrization::drift_distance::x_width_p.

                                                           {
   // NOTE: This method takes care of convertion of the units in order to use the CMSSW conventions
 
   // The convention used by the parametrization for the reference frame are different
   // with respect to CMSSW
   // X_par = X_CMSSW; Y_par=Z_CMSSW; Z_par = -Y_CMSSW 
   float By_par = bz;  // Bnorm
   float Bz_par = -by; // Bwire
  
   // Express alpha in degrees 
   float alpha_par = alpha * 180./Geom::pi();
 
   //--------------------------------------------------------------------
   // Calculate the drift distance and the resolution from the parametrization
   unsigned short flag =
     MB_DT_drift_distance(time, alpha_par, By_par, Bz_par, dx, interpolate);
 
   if(flag!=1)
     return false;
 
   // Convert from mm (used by the parametrization) to cm (used by CMSSW)
   dx->v_drift = dx->v_drift/10.;
   dx->x_drift = dx->x_drift/10.;
   dx->delta_x = dx->delta_x/10.;
   dx->x_width_m = dx->x_width_m/10.;
   dx->x_width_p = dx->x_width_p/10.;
 
   return true;
 
  }

unsigned short DTTime2DriftParametrization::MB_DT_Check_boundaries	(	double	distime,
		double	alpha,
		double	by,
		double	bz,
		short	ifl
	)		const

private

Referenced by MB_DT_drift_distance().

double DTTime2DriftParametrization::MB_DT_delta_x	(	double	t,
		double *	par
	)		const

private

Definition at line 240 of file DTTime2DriftParametrization.cc.

References edmStreamStallGrapher::t, and reco::t2.

Referenced by MB_DT_drift_distance().

                                                                               {
       double t2 = t*t;
 
 #ifdef MB_DT_DEBUG
       printf("MB_DT_delta_x:  time = %f\n",t);
 #endif
 
 //
 // Now distances are in microns. Divide by 1000 to get mm.
 //
       if ( t <= par[12] ) { return ( -0.001 * (par[1] + par[2] *t + par[3] *t2) ); }
       if ( t <= par[13] ) { return ( -0.001 * (par[4] + par[5] *t + par[6] *t2 + par[7]*t2*t + par[8]*t2*t2) ); }
                             return ( -0.001 * (par[9] + par[10]*t + par[11]*t2) );
 }

double DTTime2DriftParametrization::MB_DT_dist	(	double	time,
		double *	par
	)		const

private

Definition at line 256 of file DTTime2DriftParametrization.cc.

Referenced by MB_DT_drift_distance().

                                                                               {
       return ( time * par[0] );        // par[0] is the drift velocity, 'dist' is the linear part of the drift distance
 }

unsigned short DTTime2DriftParametrization::MB_DT_drift_distance	(	double	time,
		double	alpha,
		double	by,
		double	bz,
		drift_distance *	DX,
		short	interpolate
	)		const

private

Calculate drift distance and spread.

Definition at line 100 of file DTTime2DriftParametrization.cc.

References alpha, alpha_value, By_value, Bz_value, DTTime2DriftParametrization::drift_distance::delta_x, fun_sigma_t, fun_x, i, j, MB_DT_Check_boundaries(), MB_DT_delta_x(), MB_DT_dist(), MB_DT_Get_grid_points(), MB_DT_MLInterpolation(), MB_DT_sigma_t_m(), MB_DT_sigma_t_p(), N_Par_x, N_Sigma_t, DTTime2DriftParametrization::drift_distance::v_drift, DTTime2DriftParametrization::drift_distance::x_drift, DTTime2DriftParametrization::drift_distance::x_width_m, and DTTime2DriftParametrization::drift_distance::x_width_p.

Referenced by computeDriftDistance_mean(), and computeDriftDistance_mode().

                                                            {
       unsigned short i, j, n_func, ial, iby, ibz;
       unsigned short i_alpha, i_By, i_Bz;
       unsigned short j_alpha, j_By, j_Bz;
 
       double OffSet;
       double par_x[N_Par_x];
       double par_sigma_t[N_Sigma_t];
 
       double V_al[3], V_by[3], V_bz[3];
 
       double DXV_v_drift  [N_Func];
       double DXV_x_drift  [N_Func];
       double DXV_delta_x  [N_Func];
       double DXV_x_width_m[N_Func];
       double DXV_x_width_p[N_Func];
 
       DX->v_drift   = -1;
       DX->x_drift   = -1;
       DX->delta_x   = -1;
       DX->x_width_m = -1;
       DX->x_width_p = -1;
 
 /* Check 'interpolate' and initialize DXV */
 
       switch(interpolate) {
          case 1:  n_func = N_Func;
                   for ( j=0 ; j<N_Func ; j++ ) {
                      DXV_v_drift[j]   = -1;
                      DXV_x_drift[j]   = -1;
                      DXV_delta_x[j]   = -1;
                      DXV_x_width_m[j] = -1;
                      DXV_x_width_p[j] = -1;
                   }
                   break ;
 
          case 0:  n_func = 1;
                   break ;
 
          default: printf ("*** ERROR, MB_DT_drift_distance:  invalid interpolate value = %d\n",interpolate);
                   return (0);
       }
 
 #ifdef MB_DT_DEBUG
 
 /* Dump input values */
 
       printf ("\nMB_DT_drift_distance:  Function called with values:\n\n");
       printf ("MB_DT_drift_distance:    time  = %f\n",time );
       printf ("MB_DT_drift_distance:    alpha = %f\n",alpha);
       printf ("MB_DT_drift_distance:    by    = %f\n",by   );
       printf ("MB_DT_drift_distance:    bz    = %f\n",bz   );
 
 #endif
 
 /* Take into account the symmetries of the parametrisations */
 
       by = fabs(by);         //  f-1 (By) = f-1 (-By)
 
       if ( bz < 0 ) {        //  f-1 (alpha,Bz) = f-1 (-alpha, -Bz)
          bz    = -bz;
          alpha = -alpha;
       }
 
 /* Check boundaries of the variables and take the closest values */
 
       MB_DT_Check_boundaries (time, alpha, by, bz, -1) ;
 
       MB_DT_Get_grid_points (alpha, by, bz, &i_alpha, &i_By, &i_Bz, &j_alpha, &j_By, &j_Bz) ;
 
 #ifdef MB_DT_DEBUG
       printf("MB_DT_drift_distance:\n");
       printf("MB_DT_drift_distance:  i_alpha j_alpha alpha_value's %d %d %.0f %.0f\n",i_alpha,j_alpha,alpha_value[i_alpha],alpha_value[j_alpha]);
       printf("MB_DT_drift_distance:  i_By    j_By    By_value's    %d %d %.2f %.2f\n",i_By   ,j_By   ,   By_value[i_By]   ,   By_value[j_By]);
       printf("MB_DT_drift_distance:  i_Bz    j_Bz    Bz_value's    %d %d %.2f %.2f\n",i_Bz   ,j_Bz   ,   Bz_value[i_Bz]   ,   Bz_value[j_Bz]);
 #endif
 
 /* Get the parametrisations for the different grid points */
 
       for ( j=0 ; j<n_func ; j++ ) {
 
          ial = (j&4) ? j_alpha : i_alpha;
          iby = (j&2) ? j_By    : i_By   ;
          ibz = (j&1) ? j_Bz    : i_Bz   ;
 
          for ( i=0 ; i<N_Par_x   ; i++ ) par_x[i]       =       fun_x[ial][iby][ibz][i];
          for ( i=0 ; i<N_Sigma_t ; i++ ) par_sigma_t[i] = fun_sigma_t[ial][iby][ibz][i];
 
          OffSet = par_x[N_Par_x-1];
 
          DXV_v_drift[j]   = par_x[0];                                                       // drift velocity
          DXV_delta_x[j]   = MB_DT_delta_x (time, par_x);                                    // deviation from linearity
          DXV_x_drift[j]   = MB_DT_dist (time, par_x) + DXV_delta_x[j] - OffSet*par_x[0];    // drift distance to the wire (signed)
          DXV_x_width_m[j] = MB_DT_sigma_t_m (DXV_x_drift[j], par_sigma_t) * DXV_v_drift[j]; // distance width ~ vd * sigma_t(x)  (minus, left)
          DXV_x_width_p[j] = MB_DT_sigma_t_p (DXV_x_drift[j], par_sigma_t) * DXV_v_drift[j]; // distance width ~ vd * sigma_t(x)  (plus, right)
 
 #ifdef MB_DT_DEBUG
          printf("MB_DT_drift_LOOP: OffSet = %f\n",OffSet);
          printf("MB_DT_drift_LOOP: time, v_drift, x_drift, lin, x_width_m, x_width_p = %f %f %f %f %f %f\n",
                                    time, DXV_v_drift[j], DXV_x_drift[j], DXV_delta_x[j], DXV_x_width_m[j], DXV_x_width_p[j]) ;
 #endif
       }
 
 /* Return interpolated or grid values */
 
       if ( interpolate == 0 ) {
 
          DX->v_drift   = DXV_v_drift[0];
          DX->delta_x   = DXV_delta_x[0];
          DX->x_drift   = DXV_x_drift[0];
          DX->x_width_m = DXV_x_width_m[0];
          DX->x_width_p = DXV_x_width_p[0];
       }
       else {
 
          V_al[0] = alpha ; V_al[1] = alpha_value[i_alpha] ; V_al[2] = alpha_value[j_alpha];
          V_by[0] = by    ; V_by[1] =    By_value[i_By]    ; V_by[2] =    By_value[j_By];
          V_bz[0] = bz    ; V_bz[1] =    Bz_value[i_Bz]    ; V_bz[2] =    Bz_value[j_Bz];
 
          DX->v_drift   = MB_DT_MLInterpolation (V_al, V_by, V_bz, DXV_v_drift  );
          DX->delta_x   = MB_DT_MLInterpolation (V_al, V_by, V_bz, DXV_delta_x  );
          DX->x_drift   = MB_DT_MLInterpolation (V_al, V_by, V_bz, DXV_x_drift  );
          DX->x_width_m = MB_DT_MLInterpolation (V_al, V_by, V_bz, DXV_x_width_m);
          DX->x_width_p = MB_DT_MLInterpolation (V_al, V_by, V_bz, DXV_x_width_p);
       }
 
 #ifdef MB_DT_DEBUG
       printf("MB_DT_drift_distance: time, v_drift, x_drift, lin, x_width_m, x_width_p = %f %f %f %f %f %f\n",
                                     time, DX->v_drift, DX->x_drift, DX->delta_x, DX->x_width_m, DX->x_width_p) ;
 #endif
 
       return (1);
 }

void DTTime2DriftParametrization::MB_DT_Get_grid_points	(	double	alpha,
		double	by,
		double	bz,
		unsigned short *	p_alpha,
		unsigned short *	p_By,
		unsigned short *	p_Bz,
		unsigned short *	q_alpha,
		unsigned short *	q_By,
		unsigned short *	q_Bz
	)		const

private

Referenced by MB_DT_drift_distance().

void DTTime2DriftParametrization::MB_DT_Get_grid_values	(	double	Var,
		unsigned short *	pi,
		unsigned short *	pj,
		short	Initial,
		unsigned short	N,
		const double *	Values
	)		const