Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
Geant4 Version:11.3.2
Operating System:ubuntu 22.04
Compiler/Version:gcc/g++13.1.0
CMake Version:3.22.1
The yield of scintillation is not same
I am simulating a scintillation example : this is my code
the problem is on the prat : scintillationyield1/2/3
you can see that the total scintillationyield is 19/MeV , and the particle shot to air is alpha(5.49MeV) , this mean the “myPMT2” should produce 104 photons.
But if I annotating the code :
myMPT2->AddConstProperty(“SCINTILLATIONYIELD1”, 0.15);
myMPT2->AddConstProperty(“SCINTILLATIONYIELD2”, 0.65);
myMPT2->AddConstProperty(“SCINTILLATIONYIELD3”, 0.2);
the photon yield is 104 ; if not ,the yield is only 28-32 ; when I change the value of scintillationyield1/2/3 ,the yield still is 28-32
Could you send your entire construction file? If I had to warrant a guess I would say there is an issue with the vectors you input into SCINTILLATIONCOMPONENT1/2/3.
1 Like
OpNoviceDetectorConstruction.cc (26.4 KB)this is my code, the main part is in the “construct“ ,if you need the more information , i can send my entire code to your email, thanks very much
I was mistaken, I used your scintillationcomponent vectors in my code and got approximately the expected value 104 optical photons.
wow，that’s good，but i tried ，in my code ，If I don’t comment on those “scintillationyield1，2，3” codes, the result I get is still around 30.
May I have a look at your code.
my email is 1924925272@qq.com
thanks very much
.
i tried again
this is the result when i do not comment the scintillationyield1/2/3
this is the result when i comment the scintillationyield1/2/3
Where are you generating your primary particles? If you are not randomizing your seed, could you repetitively be seeing a rare case of some kind of escape, such as edge effects i.e. electron leaving your scorer?
i do not think my particle get out of my scorer.
If you can accept it, may I send you the entire code?
