r += dr;
} // radial steps
//
-// My length definition (is exact for hard sphere)
+// My length definition (is exact for hard disk)
return (2. * rw / w);
}
//
Double_t b = x[0];
Double_t rA = par[0];
-
- Double_t saa = (TMath::Pi() - 2. * TMath::ASin(b/ 2./ rA)) * rA * rA - b * TMath::Sqrt(rA * rA - b * b/ 4.);
+//
+// Attention: area of transverse reaction zone in hard-sphere approximation !
+ Double_t saa = (TMath::Pi() - 2. * TMath::ASin(b/ 2./ rA)) * rA * rA
+ - b * TMath::Sqrt(rA * rA - b * b/ 4.);
Double_t taa = fWStaa->Eval(b);
return (taa/saa);
// Return number of binary collisions normalized to 1 at b=0
//
- return fWSN->Eval(b)/fWSN->Eval(0.);
+ return fWSN->Eval(b)/fWSN->Eval(0.001);
}