Message: Problem with G4ExtrudedSolid on CentOS 7 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problem with G4ExtrudedSolid on CentOS 7 

Forum: Geometry
Date: 31 Oct, 2017
From: Emanuele Leonardi <Emanuele Leonardi>

In the NA62 simulation I am using a G4ExtrudedSolid to create a rail inside the experimental tube.

We are now undergoing the migration to CentOS 7. Compiling the code with gcc 6.2.0 and G4.10.3.p01 on this OS the simulation crashes issuing the following messages:

 -------- WWWW ------- UException-START -------- WWWW -------

 *** ExceptionHandler is not defined ***
 *** Exception : GeomSolids1001
       issued by : UTriangularFacet::UTriangularFacet()
 Length of sides of facet are too small.
 P[0] = (1.58556e-319, 9.50122e-317, 8.72718e-320)
 P[1] = (2.22586e-319, 9.50122e-317, 8.72718e-320)
 P[2] = (2.22586e-319, 9.50046e-317, 8.72718e-320)
 Side lengths = P[0]->P[1]0
 Side lengths = P[0]->P[2]0
 Side lengths = P[1]->P[2]0
 *** This is just a warning message ***
 -------- WWWW -------- UException-END --------- WWWW -------

 -------- WWWW ------- UException-START -------- WWWW -------

 *** ExceptionHandler is not defined ***
 *** Exception : GeomSolids1002
       issued by : UTessellatedSolid::AddFacet()
 Attempt to add facet not properly defined.
 *** This is just a warning message ***
 -------- WWWW -------- UException-END --------- WWWW -------

 *********************************************************************
 FACET TYPE       = TriangularFacet
 ABSOLUTE VECTORS =
 P[0]      = (1.58556e-319, 9.50122e-317, 8.72718e-320)
 P[1]      = (2.22586e-319, 9.50122e-317, 8.72718e-320)
 P[2]      = (2.22586e-319, 9.50046e-317, 8.72718e-320)
 *********************************************************************

 -------- EEEE ------- UException-START -------- EEEE -------

 *** ExceptionHandler is not defined ***
 *** Exception : GeomSolids0003
       issued by : UExtrudedSolid::UExtrudedSolid()
 Making facets failed - btRail
 *** Fatal Exception ***
 -------- EEEE ------- UException-END --------- EEEE -------

 *** GException: Aborting execution ***
 Aborted (core dumped)

To investigate the problem, I extracted the code used to create the rail and inserted it in an empty world using the standard B1 example from the G4 package.

The code used to create the rail is the following:

  // NA62 Rail
  G4int nV = 8;
  std::vector<G4TwoVector> v(nV);
  v[0] = G4TwoVector( 395.00*mm,-760.25*mm);
  v[1] = G4TwoVector( 555.00*mm,-760.25*mm);
  v[2] = G4TwoVector( 555.00*mm,-779.25*mm);
  v[3] = G4TwoVector( 490.00*mm,-779.25*mm);
  v[4] = G4TwoVector( 490.00*mm,-825.53*mm);
  v[5] = G4TwoVector( 460.00*mm,-842.61*mm);
  v[6] = G4TwoVector( 460.00*mm,-779.25*mm);
  v[7] = G4TwoVector( 395.00*mm,-779.25*mm);

  G4ExtrudedSolid* btRailSolid = new G4ExtrudedSolid("btRail",v,1.*m,G4TwoVector(0.,0.),1.,G4TwoVector(0.,0.),1.);
  G4LogicalVolume* btRailLV = new G4LogicalVolume(btRailSolid, nist->FindOrBuildMaterial("G4_Fe"),"BlueTubeRail",0,0,0);
  new G4PVPlacement(0,G4ThreeVector(0.,0.,0.),btRailLV,"BlueTubeRail",logicWorld,false,0,checkOverlaps);

Compiling the modified example with the "old" configuration, i.e.

    source /cvmfs/sft.cern.ch/lcg/releases/LCG_87/gcc/4.9.3/x86_64-slc6/setup.sh
    source /cvmfs/geant4.cern.ch/geant4/10.3.p01/x86_64-slc6-gcc49-opt/share/Geant4-10.3.1/geant4make/geant4make.sh

the rail is correctly created and placed in the world, while with the "new" configuration under CentOS 7, i.e.

    source /cvmfs/sft.cern.ch/lcg/contrib/gcc/6.2.0/x86_64-centos7/setup.sh ""
    source /cvmfs/sft.cern.ch/lcg/releases/LCG_89/Geant4/10.03.p01/x86_64-centos7-gcc62-opt/share/Geant4-10.3.1/geant4make/geant4make.sh

I get the crash.

Any idea of what could be wrong with this G4ExtrudedSolid?

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Problem with G4ExtrudedSolid on CentOS 7   (Evgueni Tcherniaev - 10 Nov, 2017)
 Add Message Add Message
to: "Problem with G4ExtrudedSolid on CentOS 7"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]