In G4WilsonRadius.cc file,GetWilsonRMSRadius(double A) fuction want to calculates The root-mean-square radius.But, this fuction return the radius of a nucleus. According to Tripathi, Cucinnota and Wilson,Nucl Instrum Meth Phys Res B, 152, pp425-431, 1999, RMS RADIUS = 0.84*A**(1./3.)+0.55. I am not sure that this is a bug,please check and send me an e-mail. thank you.
Although I am not the author of the code, by inspecting the code it seems fine to me. In fact, the method you refer to, G4WilsonRadius::GetWilsonRMSRadius, does correctly compute the RMS of the radius - not the radius. Perhaps you got confused by the "bad" choice of the name of the local variable used inside the method: it is named "radius", instead of using a better-looking name as, for instance, "rmsRadius", but the name does not play any role! The name of the method, instead, correctly refers to the RMS of the radius.
(In reply to Alberto.Ribon from comment #1) > Although I am not the author of the code, by inspecting the code it seems > fine to me. In fact, the method you refer to, > G4WilsonRadius::GetWilsonRMSRadius, does correctly compute the RMS of the > radius - not the radius. > Perhaps you got confused by the "bad" choice of the name of the local > variable used inside the method: it is named "radius", instead of using a > better-looking name as, for instance, "rmsRadius", but the name does not > play any role! > The name of the method, instead, correctly refers to the RMS of the radius. I'm not sure that RMS RADIUS = 0.84*A**(1./3.)+0.55 or RMS RADIUS = 1.29*0.84*A**(1./3.)+0.55? Maybe I misunderstand the formula in the reference. and thanks for your reply.
(In reply to Alberto.Ribon from comment #1) > Although I am not the author of the code, by inspecting the code it seems > fine to me. In fact, the method you refer to, > G4WilsonRadius::GetWilsonRMSRadius, does correctly compute the RMS of the > radius - not the radius. > Perhaps you got confused by the "bad" choice of the name of the local > variable used inside the method: it is named "radius", instead of using a > better-looking name as, for instance, "rmsRadius", but the name does not > play any role! > The name of the method, instead, correctly refers to the RMS of the radius. I'm not sure that RMS RADIUS = 0.84*A**(1./3.)+0.55 or RMS RADIUS = 1.29*(0.84*A**(1./3.)+0.55)? Maybe I misunderstand the formula in the reference. and thanks for your reply.
