Dear Koichi, a simplification of my MC is a Paraboloid inside a Cylinder to simulate a mirror for optical photons. All checks by "geantinos" indicate that there are no overlaps. But when I run: /geometry/test/recursive_test I have a lot of warnings of the same type: GeomTest Error: SolidProblem DistanceToOut(p,v) = kInfinity for point inside Solid name = solidMirror1 Local position = 17 -15.98 -0.05 cm When I run my macro (testmirror.in)I have: *** G4Exception : OutOfRange issued by : G4MaterialPropertyVector::GetProperty() Attempt to retrieve property below range ! but the MC ends ok. In my complete MC, when it runs, I also have the first warning. I wrote to Gabriele Cosmo, who advised me to open a bug. In attach: testmirror.tar.gz that cointains: TestMirror.cc GNUmakefile testmirror.in src.tar.gz include.tar.gz Thanks in advance Best regards Sergio Lo Meo
Can you please specify what version of Geant4 you're using ? Also, there is no attachment uploaded to this report. It is therefore impossible for us to make any investigation. Can you please provide the missing information ? The warning you get from G4MaterialPropertyVector is likely uncorrelated to a geometry problem and eventually associated to the optical physics setup of the application...
Created attachment 130 [details] tar.gz that contains main.cc macro src.tar and include.tar >Can you please specify what version of Geant4 you're using ? 4.9.3 >Also, there is no attachment uploaded to this report. It is therefore >impossible for us to make any investigation Sorry, I don't find any "upload" button in my bugzilla account. May I send to you an email? >The warning you get from G4MaterialPropertyVector is likely uncorrelated to a >geometry problem and eventually associated to the optical physics setup of the >application. The example that I made for bugzilla is a simplification of my complete setup. Perhaps, it is better to upload (or email) the tar.gz of my complete setup. Is it possible?
Created attachment 131 [details] tar.gz cointains a complete MC. During run I have: G4Exception : Notification issued by : G4Paraboloid::DistanceToIn(p,v) Point p is inside! This MC runs ok with al lot of warnings: *** G4Exception : Notification issued by : G4Paraboloid::DistanceToIn(p,v) Point p is inside! *** This is just a warning message. WARNING - G4Paraboloid::DistanceToIn(p,v) p = (30.76042335,0.5739114165,-21.63727978) mm v = (0.7561243347,-0.2456864795,0.6065592669) mm Setup: A black box (made of Plexi and inside Air)that contains: 1) Lead Glass (for Cerenkov) 2) Dummy Detector (air) 3) A Mirror inside a Cylinder 4) Detector (little box) at the bottom of Lead Glass Outside of BlackBox there is a box (Plexi) simulating a BuildUp material testcomplete.in is the macro. 1000 gamma impinging LeadGlass, then Cerenkov Optical Photons to the Mirror.
Dear Sergio, Unfortunately, I'm not able to run your complete application. I have a segmentation fault comming from HistoManager::ScoreNewTrack (). Could you check if the attachment is complete and if the attached application is working ? Best Regards, Tatiana Nikitina
Created attachment 132 [details] Complete tar.gz (after a second check and make clean) It's ok on my PC. No segmentation fault. I have made a “make clean” and then “make” Then, I run with: /home/g4/bin/Linux-g++/TestComplete /home/lomeo/TestComplete/testcomplete.in In attach the complete tar.gz This is the final summary: many warnings as: *** G4Exception : Notification issued by : G4Paraboloid::DistanceToIn(p,v) Point p is inside! *** This is just a warning message. Total Number of Optical Photons (Cerenkov) in Target: 0 Total Number of Optical Photons (Scintillation): 0 RunAction: End of run actions are started HistoManager: End of run actions are started then the summary: ======================================================== Beam particle gamma Beam Energy(MeV) 4.24546 Number of events 1000 Average energy deposit (MeV) 0.6236 RMS(MeV) 0.9309 Average number of steps 238.8 Average number of gamma 0.72 Average number of e- 48.71 Average number of e+ 0.036 Average number of neutrons 0 Average number of protons 0 Average number of antiprotons 0 Average number of pi+ & pi- 0 Average number of pi0 0 Average number of kaons 0 Average number of muons 0 Average number of deuterons+tritons 0 Average number of He3+alpha 0 Average number of ions 0 Average number of forward neutrons 0 Average number of reflected neutrons 0 Average number of leaked neutrons 0 Average number of proton leak 0 Average number of pion leak 0 Average number of optical photos (T) 31.34 ======================================================== and it’s ok, no segmentation.
Dear Sergio, Thank you for this new version, it is working for me. I will investigate these warning messages from G4Paraboloid. I have tested your application on the possible overlaps using G4PVPlacement(...,last paramter=true). There is No overlaps detected. You can have a look in Application Developers Guide, Chapter 4, to see more details about different tools to detect overlaps. In any case I will also investigate if the problem reported by 'geometry/test' is a real problem of the Solid or Not. Best Regards, Tatiana
Created attachment 133 [details] Fix for G4Paraboloid
Dear Sergio, The raison for Warnings from DistanceToIn(p,v) is found. There was an error in G4Paraboloid::DistanceToIn(p). In attachment you will find the fix for geant4.9.3 version. Could you try it out and tell if it is works for you also? Best Regards, Tatiana
Created attachment 134 [details] Fix for 9.3 version This second version is a fix for error messages from /geometry/test/recursive_test
Dear Sergio, The problem reported by /geometry/test/recursive_test was fixed. I have included fix for geant4.9.3 version. This fix will be also included in next geant4 release. Thank you for your help, Best Regards, Tatiana
