| Summary: | Hanging up in infinity loop while tracking | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | Alexander.Dietz |
| Component: | geometry | Assignee: | Gabriele Cosmo <Gabriele.Cosmo> |
| Status: | CLOSED INVALID | ||
| Severity: | major | ||
| Priority: | P2 | ||
| Version: | 3.1 | ||
| Hardware: | PC | ||
| OS: | Linux | ||
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. |
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?