BesTofDigitizerEcV4::StripStruct Struct Reference

#include <BesTofDigitizerEcV4.hh>

List of all members.

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)
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 BesTofDigitizerEcV4.hh.

---

Constructor & Destructor Documentation

BesTofDigitizerEcV4::StripStruct::StripStruct

(

)

Definition at line 657 of file BesTofDigitizerEcV4.cc.

References initial().

{
    initial();
}

---

Member Function Documentation

void BesTofDigitizerEcV4::StripStruct::avalanche

(

)

Definition at line 458 of file BesTofDigitizerEcV4.cc.

References calNextN(), BesTofDigitizerEcV4::Param::gapWidth, hitStructCol, genRecEmupikp::i, ganga-rec::j, m_param, BesTofDigitizerEcV4::Param::nstep, saturationFlag, and BesTofDigitizerEcV4::Param::stepWidth.

Referenced by BesTofDigitizerEcV4::Digitize().

{
    //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;
}

void BesTofDigitizerEcV4::StripStruct::calFirstHit

(

)

Definition at line 410 of file BesTofDigitizerEcV4.cc.

References hitStructCol, genRecEmupikp::i, tPropagate_sphi, tPropagate_xphi, trkIndex, and tStart.

Referenced by BesTofDigitizerEcV4::Digitize().

{
    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;
        }
    }
}

long int BesTofDigitizerEcV4::StripStruct::calNextN

(

int

num

)

Definition at line 512 of file BesTofDigitizerEcV4.cc.

References alpha, calSigma(), eta, exp(), genRecEmupikp::i, m_param, multiply(), and BesTofDigitizerEcV4::Param::stepWidth.

Referenced by avalanche().

{
    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;
    }
}

double BesTofDigitizerEcV4::StripStruct::calSigma

(

)

Definition at line 552 of file BesTofDigitizerEcV4.cc.

References alpha, eta, exp(), m_param, and BesTofDigitizerEcV4::Param::stepWidth.

Referenced by calNextN().

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

double BesTofDigitizerEcV4::StripStruct::getAlpha

(

double

E

)

Definition at line 727 of file BesTofDigitizerEcV4.cc.

References alpha.

Referenced by setPar().

{
    //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;
}

double BesTofDigitizerEcV4::StripStruct::getEta

(

double

E

)

Definition at line 789 of file BesTofDigitizerEcV4.cc.

References eta.

Referenced by setPar().

{
    //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;
}

double BesTofDigitizerEcV4::StripStruct::getV

(

double

E

)

Definition at line 851 of file BesTofDigitizerEcV4.cc.

References v.

Referenced by setPar().

{
    //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;
}

void BesTofDigitizerEcV4::StripStruct::initial

(

)

Definition at line 662 of file BesTofDigitizerEcV4.cc.

References alpha, E, eta, hitStructCol, strip, threshold, tPropagate_sphi, tPropagate_xphi, trkIndex, tStart, tThreshold, and v_drift.

Referenced by StripStruct().

{
    //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();
}

long int BesTofDigitizerEcV4::StripStruct::multiply

(

double

rdm

)

Definition at line 536 of file BesTofDigitizerEcV4.cc.

References alpha, eta, exp(), m_param, and BesTofDigitizerEcV4::Param::stepWidth.

Referenced by calNextN().

{
    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;
    }
}

void BesTofDigitizerEcV4::StripStruct::print

(

)

Definition at line 963 of file BesTofDigitizerEcV4.cc.

References alpha, E, eta, strip, threshold, tPropagate_sphi, tPropagate_xphi, trkIndex, tStart, tThreshold, and v_drift.

{
    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;
}

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

Definition at line 683 of file BesTofDigitizerEcV4.cc.

References alpha, E, eta, BesTofDigitizerEcV4::Param::gapWidth, getAlpha(), getEta(), getV(), m_param, saturationFlag, threshold, and v_drift.

Referenced by BesTofDigitizerEcV4::Digitize().

{
    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;
}

---

Member Data Documentation

double BesTofDigitizerEcV4::StripStruct::alpha

Definition at line 131 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::charge

Definition at line 127 of file BesTofDigitizerEcV4.hh.

Referenced by BesTofDigitizerEcV4::Digitize().

double BesTofDigitizerEcV4::StripStruct::E

Definition at line 130 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::eta

Definition at line 132 of file BesTofDigitizerEcV4.hh.

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

vector<HitStruct> BesTofDigitizerEcV4::StripStruct::hitStructCol

Definition at line 118 of file BesTofDigitizerEcV4.hh.

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

Param BesTofDigitizerEcV4::StripStruct::m_param

Definition at line 117 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), calNextN(), calSigma(), multiply(), and setPar().

bool BesTofDigitizerEcV4::StripStruct::saturationFlag

Definition at line 135 of file BesTofDigitizerEcV4.hh.

Referenced by avalanche(), and setPar().

int BesTofDigitizerEcV4::StripStruct::strip

Definition at line 121 of file BesTofDigitizerEcV4.hh.

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

long int BesTofDigitizerEcV4::StripStruct::threshold

Definition at line 133 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::tPropagate_sphi

Definition at line 124 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::tPropagate_xphi

Definition at line 125 of file BesTofDigitizerEcV4.hh.

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

int BesTofDigitizerEcV4::StripStruct::trkIndex

Definition at line 122 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::tStart

Definition at line 123 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::tThreshold

Definition at line 126 of file BesTofDigitizerEcV4.hh.

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

double BesTofDigitizerEcV4::StripStruct::v_drift

Definition at line 134 of file BesTofDigitizerEcV4.hh.

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

---