#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)

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 44 of file BesTofDigitizerEcV4.hh.

Constructor & Destructor Documentation

◆ Param()

BesTofDigitizerEcV4::Param::Param ( )

Definition at line 708 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();
}

Member Function Documentation

◆ getAlpha()

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

Definition at line 735 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;
}

◆ getEta()

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

Definition at line 798 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;
}

◆ getV()

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

Definition at line 861 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;
}

◆ print()

void BesTofDigitizerEcV4::Param::print ( )

Definition at line 924 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;
}

Referenced by BesTofDigitizerEcV4::initial().

◆ setPar()

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

Definition at line 696 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
}