When I try to simulate electrons coming from the x-direction the simulation is hang up in an infinity loop. A anti-muon-neutrino is generated which is trapped between two volumes: 296 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cmBLE1Transportation 297 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cm1PI2Transportation 298 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cmBLE1Transportation 299 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cm1PI2Transportation 300 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cmBLE1Transportation 301 -31.6 cm -7.5 cm -9.1 cm 48.6 MeV 0 eV 0 fm 7.46 cm1PI2Transportation Here are the specifications of these volumes. The materials are used also earlier (without any problem) so that the definition of the materials should be alright. solidBLE1 = new G4Box("BLE1",45.00000*cm,45.00000*cm,42.50000*cm); logicBLE1 = new G4LogicalVolume(solidBLE1, matPb,"BLE1"); physiBLE1= new G4PVPlacement(0,G4ThreeVector( ), "BLE1",logicBLE1,physiSTHL,false,0 ); solid1PI2 = new G4Tubs("1PI2",0.00000*cm,2.00000*cm,10.00000*cm,0,6.2831853); logic1PI2 = new G4LogicalVolume(solid1PI2, matAir,"1PI2"); physi1PI2= new G4PVPlacement(rotY,G4ThreeVector(-35.00000*cm,-7.50000*cm,-9.60000*cm ), "1PI2",logic1PI2,physiBLE1,false,0 ); Sa can you help me? Do you need more informations? What can I do to avoid this infinite loop?
The information provided is not sufficient to debug the problem (what's the world-volume? What's "rotY"? ...). A description of the whole geometry is required, possibly with a testcase where the problem can be easily reproduced.
It was verified the problem is due to bugs in the user geometry. Overlapping volumes were detected.
