I created a simple project to simulate the expected non-uniform response in a plastic scintillator tile read out by two SiPMs attached to wavelength shifting bars. I created a setup with the panel and WLS bars with the appropriate optical surfaces, and did a sanity check in a visualization. So far, so good. The scintillation light is generated, is partially attenuated in the panel, and is trapped in the WLS bars.
To check the result more quantitatively, I added some scoring meshes to the ends of the bars where I want the SiPMs to be. I’m still learning in regards to using the scorers, so I started by prototyping with the commands interactively. Now if I place a scorer, draw it, and run a new event, the visualization changes drastically. Specifically, it seems like the optical surface properties of my volumes change. The scintillation light is no longer reflecting inside the panel, but passes straight through it as if there were no wrapping.
For instance, I open a gui session and run
/gps/energy 3 GeV
I then add a scorer and run again:
/score/mesh/boxSize 7.5 7.5 0.5 mm
/score/mesh/nBin 1 1 1
/score/mesh/translate/xyz 27.005025253169418 28.065685424949237 0 cm
/score/mesh/rotate/rotateZ 45 deg
/score/mesh/rotate/rotateX 90 deg
/score/quantity/flatSurfaceCurrent sipm1current 1 0 0
/score/filter/particle phot opticalphoton
On the left is the run before the scoring mesh is added, on the right is after. This is the only change I have made:
Can anyone help me understand why this is happening? I can probably get the information I want out of the simulation without the scoring, but I would rather learn how to use the built-in tool for future projects.
Thanks for any help!