#include <BesTofDigitizerEcV4.hh>

Public Member Functions
	StripStruct ()

void	initial ()

void	print ()

void	avalanche ()

void	calFirstHit ()

long int	calNextN (int num)

long int	multiply (double rdm)

double	calSigma ()

void	setPar (double V, double threshold, bool saturationFlag=true)

double	getAlpha (double E)

double	getEta (double E)

double	getV (double E)

	StripStruct ()

void	initial ()

void	print ()

void	avalanche ()

void	calFirstHit ()

long int	calNextN (int num)

long int	multiply (double rdm)

double	calSigma ()

void	setPar (double V, double threshold, bool saturationFlag=true)

double	getAlpha (double E)

double	getEta (double E)

double	getV (double E)

Public Attributes
Param	m_param

vector< HitStruct >	hitStructCol

int	strip

int	trkIndex

double	tStart

double	tPropagate_sphi

double	tPropagate_xphi

double	tThreshold

double	charge

double	E

double	alpha

double	eta

long int	threshold

double	v_drift

bool	saturationFlag

Detailed Description

Definition at line 101 of file InstallArea/include/TofSim/TofSim/BesTofDigitizerEcV4.hh.

Constructor & Destructor Documentation

◆ StripStruct() [1/2]

BesTofDigitizerEcV4::StripStruct::StripStruct ( )

Definition at line 657 of file BesTofDigitizerEcV4.cc.

{
    initial();
}

◆ StripStruct() [2/2]

BesTofDigitizerEcV4::StripStruct::StripStruct ( )

Member Function Documentation

◆ avalanche() [1/2]

void BesTofDigitizerEcV4::StripStruct::avalanche ( )

Definition at line 458 of file BesTofDigitizerEcV4.cc.

{
    //process each hit
    for(unsigned int i=0; i<hitStructCol.size(); i++)
    {
        hitStructCol[i].ava_pos.clear();
        hitStructCol[i].ava_num.clear();
 
        hitStructCol[i].ava_pos[0] = hitStructCol[i].calAvaLength();
        hitStructCol[i].ava_num[0] = hitStructCol[i].ions;
        //cout<<"i= "<<i<<"  gap= "<<hitStructCol[i].gap<<"  initial pos= "<<hitStructCol[i].ava_pos[0]<<endl;
        for(int j=1; j<m_param.nstep; j++)
        {
            hitStructCol[i].ava_pos[j] = hitStructCol[i].ava_pos[j-1] + m_param.stepWidth;
            if(saturationFlag==true && hitStructCol[i].ava_num[j-1]>1.5e+07) //saturation e+07~e+08, ~2pC, Reather limit
            {
                hitStructCol[i].ava_num[j] = hitStructCol[i].ava_num[j-1];
            }
            else
            {
                hitStructCol[i].ava_num[j] = calNextN(hitStructCol[i].ava_num[j-1]);
            }
            if(hitStructCol[i].ava_pos[j]>m_param.gapWidth) break;
        }
    }
 
    //decide threshold and charge
    bool over_threshold = false;
    long int sum = 0;
    for(int i=0; i<m_param.nstep; i++)
    {
        for(unsigned int j=0; j<hitStructCol.size(); j++)
        {
            if(i<hitStructCol[j].ava_pos.size() && hitStructCol[j].ava_pos[i]<m_param.gapWidth) 
            {
                sum += hitStructCol[j].ava_num[i];
            }
        }
        //cout<<"sum= "<<sum<<"  avaSize= "<<hitStructCol.size()<<endl;
 
        if(over_threshold==false)
        {
            if(sum>threshold)
            {
                over_threshold = true;
                tThreshold = (m_param.gapWidth)/(m_param.nstep)/v_drift*(i+1);
            }
        }
    }
 
    charge = sum*(m_param.E_weight)*(v_drift)*(m_param.eCharge)*(m_param.gapWidth)/(m_param.nstep)/v_drift;
}

Referenced by BesTofDigitizerEcV4::Digitize().

◆ avalanche() [2/2]

void BesTofDigitizerEcV4::StripStruct::avalanche ( )

◆ calFirstHit() [1/2]

void BesTofDigitizerEcV4::StripStruct::calFirstHit ( )

Definition at line 410 of file BesTofDigitizerEcV4.cc.

{
    for(unsigned int i=0; i<hitStructCol.size(); i++)
    {
        if(hitStructCol[i].glbTime<tStart)
        {
            tStart = hitStructCol[i].glbTime;
            trkIndex = hitStructCol[i].trkIndex;
            hitStructCol[i].calTPropagate();
            tPropagate_sphi = hitStructCol[i].tPropagate_sphi;
            tPropagate_xphi = hitStructCol[i].tPropagate_xphi;
        }
    }
}

Referenced by BesTofDigitizerEcV4::Digitize().

◆ calFirstHit() [2/2]

void BesTofDigitizerEcV4::StripStruct::calFirstHit ( )

◆ calNextN() [1/2]

long int BesTofDigitizerEcV4::StripStruct::calNextN ( int num )

Definition at line 512 of file BesTofDigitizerEcV4.cc.

{
    if(num<150)
    {
        long int nextN = 0;
        double rdm;
        for(int i=0; i<num; i++)
        {
            rdm = G4UniformRand();
            nextN += multiply(rdm);
        }
        return nextN;
    }
    else
    {
        double nbar = exp((alpha-eta)*m_param.stepWidth);
        double sigma = calSigma();
        double mean = num*nbar;
        double resolution = G4RandGauss::shoot(0,(sqrt(num)*sigma));
        long int nextN = mean+resolution;
        return nextN;
    }
}

◆ calNextN() [2/2]

long int BesTofDigitizerEcV4::StripStruct::calNextN ( int num )

◆ calSigma() [1/2]

double BesTofDigitizerEcV4::StripStruct::calSigma ( )

Definition at line 552 of file BesTofDigitizerEcV4.cc.

{
    double nbar = exp((alpha-eta)*m_param.stepWidth);
    double k = eta/alpha;
    double sigma = sqrt((1+k)/(1-k)*nbar*(nbar-1));
    return sigma;
}

◆ calSigma() [2/2]

double BesTofDigitizerEcV4::StripStruct::calSigma ( )

◆ getAlpha() [1/2]

double BesTofDigitizerEcV4::StripStruct::getAlpha ( double E )

Definition at line 727 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);
    return alphaVal;
}

◆ getAlpha() [2/2]

double BesTofDigitizerEcV4::StripStruct::getAlpha ( double E )

◆ getEta() [1/2]

double BesTofDigitizerEcV4::StripStruct::getEta ( double E )

Definition at line 789 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);
    return etaVal;
}

◆ getEta() [2/2]

double BesTofDigitizerEcV4::StripStruct::getEta ( double E )

◆ getV() [1/2]

double BesTofDigitizerEcV4::StripStruct::getV ( double E )

Definition at line 851 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);
    return vVal;
}

◆ getV() [2/2]

double BesTofDigitizerEcV4::StripStruct::getV ( double E )

◆ initial() [1/2]

void BesTofDigitizerEcV4::StripStruct::initial ( )

Definition at line 662 of file BesTofDigitizerEcV4.cc.

{
    //properties to get
    strip = -999.0;
    trkIndex = -999.0;
    tStart = 99999.0;
    tPropagate_sphi = -999.0;
    tPropagate_xphi = -999.0;
    tThreshold = -999.0;
    charge = -999.0;
 
    //parameters to tune
    E = 106;
    alpha = 144800./1000; //-999.0; /mm^-1
    eta = 5013./1000; //-999.0; /mm^-1
    threshold = 1.5e+08; //Correspond to induced charge of 15 fC
    v_drift = 210.9e-3; // mm/ns
 
    hitStructCol.clear();
}

◆ initial() [2/2]

void BesTofDigitizerEcV4::StripStruct::initial ( )

◆ multiply() [1/2]

long int BesTofDigitizerEcV4::StripStruct::multiply ( double rdm )

Definition at line 536 of file BesTofDigitizerEcV4.cc.

{
    double nbar = exp((alpha-eta)*m_param.stepWidth);
    double k = eta/alpha;
    double rdm_border = k*(nbar-1)/(nbar-k);
    if(rdm<rdm_border)
    {
        return 0;
    }
    else
    {
        long int number = 1.+1./log((nbar-1.)/(nbar-k))*log((nbar-k)*(rdm-1)/(k-1)/nbar);
        return number;
    }
}

◆ multiply() [2/2]

long int BesTofDigitizerEcV4::StripStruct::multiply ( double rdm )

◆ print() [1/2]

void BesTofDigitizerEcV4::StripStruct::print ( )

Definition at line 963 of file BesTofDigitizerEcV4.cc.

{
    cout<<"Strip information:  "<<endl;
    cout<<"  strip= "<<strip
        <<"  trkIndex= "<<trkIndex
        <<"  tStart= "<<tStart
        <<"  tPropagate_sphi= "<<tPropagate_sphi
        <<"  tPropagate_xphi= "<<tPropagate_xphi
        <<"  tThreshold "<<tThreshold
        <<"  charge= "<<charge
        <<"  E= "<<E
        <<"  alpha= "<<alpha
        <<"  eta= "<<eta
        <<"  threshold= "<<threshold
        <<"  v_drift= "<<v_drift
        <<endl;
}

◆ print() [2/2]

void BesTofDigitizerEcV4::StripStruct::print ( )

◆ setPar() [1/2]

void BesTofDigitizerEcV4::StripStruct::setPar	(	double	V,
		double	threshold,
		bool	saturationFlag = `true`
	)

Definition at line 683 of file BesTofDigitizerEcV4.cc.

{
    threshold = threshold_n;
    E = E_V/1000*2/6/(m_param.gapWidth/10);  //kV/cm
    alpha = getAlpha(E); //mm^-1
    eta = getEta(E); //mm^-1
    v_drift = getV(E); // mm/ns
 
    saturationFlag = saturationFlag_n;
}