Пример #1
0
//_____________________________________________________________________________
Gammagammaleptonpair::Gammagammaleptonpair(const inputParameters& inputParametersInstance, randomGenerator* randy, beamBeamSystem& bbsystem)
: eventChannel(inputParametersInstance, randy, bbsystem)
, _GGlepInputpidtest(inputParametersInstance.prodParticleType())
, _GGlepInputnumw(inputParametersInstance.nmbWBins())
, _GGlepInputnumy(inputParametersInstance.nmbRapidityBins())
, _GGlepInputGamma_em(inputParametersInstance.beamLorentzGamma())
{
    //Initialize randomgenerator with our seed.
    _randy->Rndom();
    //Let us read in the luminosity tables
    read();
    //Now we will calculate the crosssection
    twoLeptonCrossSection();
    //If it is a tauon, calculate its tables
    if(inputParametersInstance.prodParticleId()==starlightConstants::TAUONDECAY) calculateTable();
}