#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 659 of file BesTofDigitizerEcV4.cc.

{
    initial();
}

Member Function Documentation

◆ avalanche()

void BesTofDigitizerEcV4::StripStruct::avalanche ( )

Definition at line 460 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 412 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 514 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;
    }
}

◆ calSigma()

double BesTofDigitizerEcV4::StripStruct::calSigma ( )

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

◆ initial()

void BesTofDigitizerEcV4::StripStruct::initial ( )

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

◆ print()

void BesTofDigitizerEcV4::StripStruct::print ( )

Definition at line 980 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 685 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.

◆ charge

double BesTofDigitizerEcV4::StripStruct::charge

Definition at line 129 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize().

◆ E

double BesTofDigitizerEcV4::StripStruct::E

Definition at line 132 of file BesTofDigitizerEcV4.hh.

◆ eta

double BesTofDigitizerEcV4::StripStruct::eta

Definition at line 134 of file BesTofDigitizerEcV4.hh.

◆ hitStructCol

vector<HitStruct> BesTofDigitizerEcV4::StripStruct::hitStructCol

Definition at line 120 of file BesTofDigitizerEcV4.hh.

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

◆ m_param

Param BesTofDigitizerEcV4::StripStruct::m_param

Definition at line 119 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize().

◆ saturationFlag

bool BesTofDigitizerEcV4::StripStruct::saturationFlag

Definition at line 137 of file BesTofDigitizerEcV4.hh.

◆ strip

int BesTofDigitizerEcV4::StripStruct::strip

Definition at line 123 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize().

◆ threshold

long int BesTofDigitizerEcV4::StripStruct::threshold

Definition at line 135 of file BesTofDigitizerEcV4.hh.

◆ tPropagate_sphi

double BesTofDigitizerEcV4::StripStruct::tPropagate_sphi

Definition at line 126 of file BesTofDigitizerEcV4.hh.

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

◆ tPropagate_xphi

double BesTofDigitizerEcV4::StripStruct::tPropagate_xphi

Definition at line 127 of file BesTofDigitizerEcV4.hh.

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

◆ trkIndex

int BesTofDigitizerEcV4::StripStruct::trkIndex

Definition at line 124 of file BesTofDigitizerEcV4.hh.

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

◆ tStart

double BesTofDigitizerEcV4::StripStruct::tStart

Definition at line 125 of file BesTofDigitizerEcV4.hh.

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

◆ tThreshold

double BesTofDigitizerEcV4::StripStruct::tThreshold

Definition at line 128 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize().

◆ v_drift

double BesTofDigitizerEcV4::StripStruct::v_drift

Definition at line 136 of file BesTofDigitizerEcV4.hh.

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