CEA Codes forum

Dear professor,

When Plotwin read particles from an ASICII file informat 1, I am not sure how does the dp/p change into W.
If dp/p=0, W=W0. But when dp/p<-0.96, it shows negative energy. Others, the output W seems not exactly equal the value calculated with dp(p=beta*gamma*m0c). There must be some misunderstandings.

Dear Zhaoliangchao,

The correct formula is dW=(dp/p)*beta2*gamma*m0c
In your case, with a particle whose dp/p=0.96 you are very far from the dp/p<<1 condition, so you must answer "Yes" to the question asked.

Regards,

Didier