#include <BesTofDigitizerEcV4.hh>

Public Member Functions
	Param ()
void	print ()
void	setPar (int nstep, double E_weight, double E)
double	getAlpha (double E)
double	getEta (double E)
double	getV (double E)
	Param ()
void	print ()
void	setPar (int nstep, double E_weight, double E)
double	getAlpha (double E)
double	getEta (double E)
double	getV (double E)
	Param ()
void	print ()
void	setPar (int nstep, double E_weight, double E)
double	getAlpha (double E)
double	getEta (double E)
double	getV (double E)

Public Attributes
BesTofGeoParameter *	tofPara
int	nstrip
int	nmodule
double	strip_x [12]
double	strip_z
double	strip_gap
int	deadChannel [72][12]
double	ngap
double	gapWidth
int	nstep
double	stepWidth
double	E_weight
double	eCharge
double	E
double	alpha
double	eta
double	v_drift

Detailed Description

Definition at line 45 of file InstallArea/x86_64-el9-gcc13-dbg/include/TofSim/BesTofDigitizerEcV4.hh.

Constructor & Destructor Documentation

◆ Param() [1/3]

BesTofDigitizerEcV4::Param::Param ( )

Definition at line 645 of file BesTofDigitizerEcV4.cc.

                                {
  // parameters fixed
  tofPara = BesTofGeoParameter::GetInstance();
  nstrip  = 12;
  nmodule = 72;
  std::stringstream ss;
  for ( int i = 0; i < nstrip; i++ )
  {
    ss.str( "" );
    ss << "strip_x[" << i << "]";
    strip_x[i] = tofPara->Get( ss.str().c_str() ); // mm
  }
  strip_z   = tofPara->Get( "strip_z" );
  strip_gap = tofPara->Get( "strip_gap" );
 
  ngap      = 12;
  gapWidth  = 0.22; // mm
  nstep     = 200;
  stepWidth = gapWidth / nstep;
  E_weight  = 1. / ( 6. * 0.22 + ( 5. * 0.4 + 2 * 0.55 ) / 3.7 ); // V/mm
  eCharge   = 1.60217733e-7;                                      // pC
  tofPara->Get_deadChannel( deadChannel );
 
  // print();
}

◆ Param() [2/3]

BesTofDigitizerEcV4::Param::Param ( )

◆ Param() [3/3]

BesTofDigitizerEcV4::Param::Param ( )

Member Function Documentation

◆ getAlpha() [1/3]

double BesTofDigitizerEcV4::Param::getAlpha ( double E )

Definition at line 671 of file BesTofDigitizerEcV4.cc.

                                                    {
  // electric field: kV/cm; alpha: mm-1
  // kV/cm
  double e[22] = { 65,  70,  75,  80,  85,  90,  95,  100, 102, 104, 106,
                   108, 110, 112, 114, 116, 118, 120, 125, 130, 135, 140 };
 
  // mm-1
  double alpha[22] = {
      383.5 / 10, 471 / 10,  564.5 / 10, 663.6 / 10, 777.1 / 10, 877 / 10,
      990.8 / 10, 1106 / 10, 1154 / 10,  1199 / 10,  1253 / 10,  1296 / 10,
      1344 / 10,  1396 / 10, 1448 / 10,  1502 / 10,  1545 / 10,  1597 / 10,
      1726 / 10,  1858 / 10, 1992 / 10,  2124 / 10,
  };
 
  TSpline3* sp_alpha = new TSpline3( "sp_alpha", e, alpha, 22 );
  double    alphaVal = sp_alpha->Eval( E );
  delete sp_alpha;
  return alphaVal;
}

Referenced by setPar().

◆ getAlpha() [2/3]

double BesTofDigitizerEcV4::Param::getAlpha ( double E )

◆ getAlpha() [3/3]

double BesTofDigitizerEcV4::Param::getAlpha ( double E )

◆ getEta() [1/3]

double BesTofDigitizerEcV4::Param::getEta ( double E )

Definition at line 691 of file BesTofDigitizerEcV4.cc.

                                                  {
  // electric field: kV/cm; eta: mm-1
  // kV/cm
  double e[22] = { 65,  70,  75,  80,  85,  90,  95,  100, 102, 104, 106,
                   108, 110, 112, 114, 116, 118, 120, 125, 130, 135, 140 };
 
  // mm-1
  double eta[22] = { 132.6 / 10, 117.2 / 10, 102.6 / 10, 88.26 / 10, 79.81 / 10, 74.0 / 10,
                     66.7 / 10,  62.7 / 10,  61.4 / 10,  57.4 / 10,  55.45 / 10, 54.35 / 10,
                     52.48 / 10, 51.3 / 10,  50.1 / 10,  48.3 / 10,  48.28 / 10, 46.00 / 10,
                     44.08 / 10, 41.67 / 10, 39.97 / 10, 38.04 / 10 };
 
  TSpline3* sp_eta = new TSpline3( "sp_eta", e, eta, 22 );
  double    etaVal = sp_eta->Eval( E );
  delete sp_eta;
  return etaVal;
}

Referenced by setPar().

◆ getEta() [2/3]

double BesTofDigitizerEcV4::Param::getEta ( double E )

◆ getEta() [3/3]

double BesTofDigitizerEcV4::Param::getEta ( double E )

◆ getV() [1/3]

double BesTofDigitizerEcV4::Param::getV ( double E )

Definition at line 709 of file BesTofDigitizerEcV4.cc.

                                                {
  // electric field: kV/cm; velocity: mm/ns
  // kV/cm
  double e[22] = { 65,  70,  75,  80,  85,  90,  95,  100, 102, 104, 106,
                   108, 110, 112, 114, 116, 118, 120, 125, 130, 135, 140 };
 
  // mm/ns
  double v[22] = { 130.2 / 1000, 138.5 / 1000, 146.7 / 1000, 155.0 / 1000, 163.3 / 1000,
                   171.4 / 1000, 179.7 / 1000, 187.7 / 1000, 191.2 / 1000, 194.5 / 1000,
                   197.9 / 1000, 201.2 / 1000, 204.5 / 1000, 207.6 / 1000, 210.9 / 1000,
                   214.4 / 1000, 217.5 / 1000, 220.9 / 1000, 228.8 / 1000, 237.0 / 1000,
                   244.7 / 1000, 252.9 / 1000 };
 
  TSpline3* sp_v = new TSpline3( "sp_v", e, v, 22 );
  double    vVal = sp_v->Eval( E );
  delete sp_v;
  return vVal;
}

Referenced by setPar().

◆ getV() [2/3]

double BesTofDigitizerEcV4::Param::getV ( double E )

◆ getV() [3/3]

double BesTofDigitizerEcV4::Param::getV ( double E )

◆ print() [1/3]

void BesTofDigitizerEcV4::Param::print ( )

Definition at line 728 of file BesTofDigitizerEcV4.cc.

                                     {
  cout << "Fixed parameters: " << endl;
  for ( int i = 0; i < nstrip; i++ ) { cout << "  strip_x[" << i << "]= " << strip_x[i]; }
  for ( int i = 0; i < nmodule; i++ )
  {
    for ( int j = 0; j < nstrip; j++ )
    {
      if ( deadChannel[i][j] != -999 )
      { cout << "  deadChannel[" << i << "][" << j << "]= " << deadChannel[i][j]; }
    }
  }
 
  cout << "  strip_z= " << strip_z << "  strip_gap= " << strip_gap << "  ngap= " << ngap
       << "  gapWidth= " << gapWidth << "  nstep= " << nstep << "  stepWidth= " << stepWidth
       << "  E_weight= " << E_weight << "  eCharge= " << eCharge << "  E= " << E
       << "  alpha= " << alpha << "  eta= " << eta << "  v_drift= " << v_drift << endl;
}

◆ print() [2/3]

void BesTofDigitizerEcV4::Param::print ( )

◆ print() [3/3]

void BesTofDigitizerEcV4::Param::print ( )

◆ setPar() [1/3]

void BesTofDigitizerEcV4::Param::setPar	(	int	nstep,
		double	E_weight,
		double	E )

Definition at line 634 of file BesTofDigitizerEcV4.cc.

                                                                                  {
  if ( nstep_n > 0 ) nstep = nstep_n;
  if ( E_weight_n > 0 ) E_weight = E_weight_n;
 
  E            = E_n;
  double E_eff = E / 1000 * 2 / 6 / ( gapWidth / 10 ); // kV/cm
  alpha        = getAlpha( E_eff );                    // mm^-1
  eta          = getEta( E_eff );                      // mm^-1
  v_drift      = getV( E_eff );                        // mm/ns
}