Int_t ws=0;
while(ret>e){
ret=fHistos[ipart-1][l-1]->GetRandom();
- if(++ws==1e5){
- Warning("GetELossRandomK",
- "Aborted reweighting; maximum loss assigned after 1e5 trials.");
+ if(++ws==1e6){
+ Warning("GetELossRandom",
+ "Aborted reweighting; maximum loss assigned after 1e6 trials.");
return e;
}
}
Int_t ws=0;
while(ret>e){
ret=fHisto->GetRandom();
- if(++ws==1e5){
+ if(++ws==1e6){
Warning("GetELossRandomK",
- "Aborted reweighting; maximum loss assigned after 1e5 trials.");
+ "Aborted reweighting; maximum loss assigned after 1e6 trials.");
return e;
}
}
{return CalcQk(fK,I0,I1);}
Double_t CalcQk(Double_t k, Double_t I0, Double_t I1) const
- {return I0*I0/2/I1/fgkConvFmToInvGeV/fgkConvFmToInvGeV*k;}
+ {return I0*I0/2/I1/fgkConvFmToInvGeV/fgkConvFmToInvGeV*k;}
+
+ Double_t CalcLk(Double_t i0, Double_t i1) const
+ {return 2.*i1/i0;}
Int_t CalcLengthMax(Double_t q) const
{Double_t l3max=fgkRMax/.5/q/fgkConvFmToInvGeV/fgkConvFmToInvGeV;