#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 alpha_n, double eta_n, double drift_n, double threshold, bool saturationFlag=true)

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

Constructor & Destructor Documentation

◆ StripStruct()

BesTofDigitizerEcV4::StripStruct::StripStruct ( )

Definition at line 660 of file BesTofDigitizerEcV4.cc.

{
    initial();
}

Member Function Documentation

◆ avalanche()

void BesTofDigitizerEcV4::StripStruct::avalanche ( )

Definition at line 461 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().

◆ calFirstHit()

void BesTofDigitizerEcV4::StripStruct::calFirstHit ( )

Definition at line 413 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().

◆ calNextN()

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

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

Referenced by avalanche().

◆ calSigma()

double BesTofDigitizerEcV4::StripStruct::calSigma ( )

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

Referenced by calNextN().

◆ initial()

void BesTofDigitizerEcV4::StripStruct::initial ( )

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

◆ multiply()

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

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

Referenced by calNextN().

◆ print()

void BesTofDigitizerEcV4::StripStruct::print ( )

Definition at line 981 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
        <<"  alpha= "<<alpha
        <<"  eta= "<<eta
        <<"  threshold= "<<threshold
        <<"  v_drift= "<<v_drift
        <<endl;
}

◆ setPar()

void BesTofDigitizerEcV4::StripStruct::setPar	(	double	alpha_n,
		double	eta_n,
		double	drift_n,
		double	threshold,
		bool	saturationFlag = true )

Definition at line 686 of file BesTofDigitizerEcV4.cc.

{
    alpha = alpha_n;
    eta = eta_n;
    v_drift = drift_n;
    threshold = threshold_n;
 
    saturationFlag = saturationFlag_n;
}

Referenced by BesTofDigitizerEcV4::Digitize().

Member Data Documentation

◆ alpha

double BesTofDigitizerEcV4::StripStruct::alpha

Definition at line 133 of file BesTofDigitizerEcV4.hh.

Referenced by calNextN(), calSigma(), initial(), multiply(), print(), and setPar().

◆ charge

double BesTofDigitizerEcV4::StripStruct::charge

Definition at line 129 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ E

double BesTofDigitizerEcV4::StripStruct::E

Definition at line 132 of file BesTofDigitizerEcV4.hh.

Referenced by initial().

◆ eta

double BesTofDigitizerEcV4::StripStruct::eta

Definition at line 134 of file BesTofDigitizerEcV4.hh.

Referenced by calNextN(), calSigma(), initial(), multiply(), print(), and setPar().

◆ hitStructCol

vector<HitStruct> BesTofDigitizerEcV4::StripStruct::hitStructCol

Definition at line 120 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), calFirstHit(), BesTofDigitizerEcV4::Digitize(), and initial().

◆ m_param

Param BesTofDigitizerEcV4::StripStruct::m_param

Definition at line 119 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), calNextN(), calSigma(), BesTofDigitizerEcV4::Digitize(), and multiply().

◆ saturationFlag

bool BesTofDigitizerEcV4::StripStruct::saturationFlag

Definition at line 137 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), and setPar().

◆ strip

int BesTofDigitizerEcV4::StripStruct::strip

Definition at line 123 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ threshold

long int BesTofDigitizerEcV4::StripStruct::threshold

Definition at line 135 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), initial(), print(), and setPar().

◆ tPropagate_sphi

double BesTofDigitizerEcV4::StripStruct::tPropagate_sphi

Definition at line 126 of file BesTofDigitizerEcV4.hh.

Referenced by calFirstHit(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ tPropagate_xphi

double BesTofDigitizerEcV4::StripStruct::tPropagate_xphi

Definition at line 127 of file BesTofDigitizerEcV4.hh.

Referenced by calFirstHit(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ trkIndex

int BesTofDigitizerEcV4::StripStruct::trkIndex

Definition at line 124 of file BesTofDigitizerEcV4.hh.

Referenced by calFirstHit(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ tStart

double BesTofDigitizerEcV4::StripStruct::tStart

Definition at line 125 of file BesTofDigitizerEcV4.hh.

Referenced by calFirstHit(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ tThreshold

double BesTofDigitizerEcV4::StripStruct::tThreshold

Definition at line 128 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), BesTofDigitizerEcV4::Digitize(), initial(), and print().

◆ v_drift

double BesTofDigitizerEcV4::StripStruct::v_drift

Definition at line 136 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), initial(), print(), and setPar().

The documentation for this struct was generated from the following files:

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ StripStruct()

Member Function Documentation

◆ avalanche()

◆ calFirstHit()

◆ calNextN()

◆ calSigma()

◆ initial()

◆ multiply()

◆ print()

◆ setPar()

Member Data Documentation

◆ alpha

◆ charge

◆ E

◆ eta

◆ hitStructCol

◆ m_param

◆ saturationFlag

◆ strip

◆ threshold

◆ tPropagate_sphi

◆ tPropagate_xphi

◆ trkIndex

◆ tStart

◆ tThreshold

◆ v_drift