EXP1803

Виталик!

Нам понадобится Гадаст, а там энергетическое разрешение применяется к пойнтам, а надо к сумме по всем пойнтам в заданном кубике.

Кроме того, надо применить то, что мы делали в стильбенах:

    Float_t edepSigma = sqrt(pow(fEdepErrorA,2) + pow(fEdepErrorB*TMath::Sqrt(edep/1000.),2) + pow(fEdepErrorC*edep,2));
    edep = gRandom->Gaus(edep, edepSigma);

    /*
    fEdepErrorB and fLYErrorB are expressed in per cent divided by square root of edep or LY,
    expressed in MeVlike it is usually done in detector description. e.g. for Edep=0.001 (i.e. 1 Mev)
    sigma is equal0.04*0.001Gev = 0.04 MeV
    */

Сейчас в ГАДАСТ неправильно:

  Float_t disp = a*a*edep + b*b*edep*edep;
    edep = gRandom->Gaus(edep*LC, disp);