| Summary: | DistanceToOut(p,v) = kInfinity for Paraboloid Mirror | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | Sergio Lo Meo <sergio.lomeo> |
| Component: | geometry/solids | Assignee: | tatiana.nikitina |
| Status: | RESOLVED FIXED | ||
| Severity: | normal | CC: | sergio.lomeo |
| Priority: | P3 | ||
| Version: | 9.4 | ||
| Hardware: | PC | ||
| OS: | Linux | ||
| Attachments: |
tar.gz that contains main.cc macro src.tar and include.tar
tar.gz cointains a complete MC. During run I have: G4Exception : Notification issued by : G4Paraboloid::DistanceToIn(p,v) Point p is inside! Complete tar.gz (after a second check and make clean) It's ok on my PC. No segmentation fault. Fix for G4Paraboloid Fix for 9.3 version |
||
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 |
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