Problem 1176

Summary: Binary Cascade model error while running an hadrontherapy revised example
Product: Examples/Advanced Reporter: Adele <adele.rimoldi>
Component: hadrontherapyAssignee: cirrone
Status: RESOLVED WORKSFORME    
Severity: normal CC: Gunter.Folger
Priority: P1    
Version: 9.4   
Hardware: Apple   
OS: Mac OS X   

Description Adele 2011-02-22 11:22:51 CET
this error happens when running an hadrontherapy revised example using geant4 9.4b01.
the program bombs after a successful run of about 700 events in a subdetector which material is :

//Cristallino (Materiale)
G4Material* CristallinoMat = new G4Material("CristallinoMat", density= 1.067*g/cm3, nel=2);
CristallinoMat->AddMaterial(H2O,       60*perCent);
CristallinoMat->AddMaterial(Proteine,  40*perCent);

//Proteine
G4Material* Proteine = new G4Material("Proteine", density=1*g/cm3, nel=5);
Proteine->AddElement(H, 50*perCent);
Proteine->AddElement(C, 28*perCent);
Proteine->AddElement(O, 13*perCent);
Proteine->AddElement(N,  8*perCent);
Proteine->AddElement(S,  1*perCent);

the same happened for another material:
//Retina (Materiale)
G4Material* RetinaMat = new G4Material("RetinaMat", density= 1.0174*g/cm3, nel=4);
RetinaMat->AddMaterial(H2O, 80*perCent);
RetinaMat->AddMaterial(NAA, 10*perCent);
RetinaMat->AddMaterial(Cho,  5*perCent);
RetinaMat->AddMaterial(Creatina, 5*perCent);


where:
//Colina
G4Material* Cho = new G4Material("Colina", density=1*g/cm3, nel=4);
Cho->AddElement(H, natoms=14);
Cho->AddElement(C, natoms=5);
Cho->AddElement(O, natoms=1);
Cho->AddElement(N, natoms=1);

//Creatina
G4Material* Creatina = new G4Material("Creatina", density=1*g/cm3, nel=4);
Creatina->AddElement(H, natoms=9);
Creatina->AddElement(C, natoms=4);
Creatina->AddElement(O, natoms=2);
Creatina->AddElement(N, natoms=3);


G4Fragment::CalculateExcitationEnergy(): Excitation Energy = -17.977786 MeV for A = 2 and Z= 2
In src/G4Fragment.cc, line 234:
===> G4Fragment::G4Fragment Excitation Energy < 0.0!
Call for Binary Cascade
Target element Hydrogen  Z= 1  A= 1
Unrecoverable error in the method ApplyYourself of hInelastic
TrackID= 1  ParentID= 0  proton
Ekin(GeV)= 0.036859948;  direction= (0.99949149,0.023409148,-0.021651078)
Position(mm)= (6.2444002,2.3479921,-10.126127);  material CristallinoMat
PhysicalVolume  <Encefalo>

*** G4Exception : 007
    issued by : G4HadronicProcess
PostStepDoIt failed.
*** Fatal Exception *** core dump ***

*** G4Exception: Aborting execution ***
Comment 1 Gabriele Cosmo 2011-02-22 16:07:38 CET
Can you please try to run your application on the final release 9.4 and see if the problem still persists ?
Comment 2 Adele 2011-02-23 10:23:54 CET
Hi,
  I have installed g4 9.4 on my Mac  and my job runs successfully. Apparently the problem has gone.
Adele