Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

unknown particle type 50 being generated through interactions #187

Closed
rjones30 opened this issue Mar 21, 2021 · 5 comments
Closed

unknown particle type 50 being generated through interactions #187

rjones30 opened this issue Mar 21, 2021 · 5 comments

Comments

@rjones30
Copy link

The following error shows up in simulations of events containing certain types of primary particles. This message is going to stderr in hdgeant4.
ParticleMass: Error: Unknown particle type 50, returning HUGE_VAL...
It happens in bursts of order 10 at a time, as if the particle were being generated in bunches. In Geant3, this is the particle type for optical photons, eg. Cerenkov photons or scintillation. I am not aware of any explicit simulation of scintillation light in hdgeant4, so it is probably some residue of G3 code that made it into G4 for the dirc.

@rjones30
Copy link
Author

The only possible place this could happen is in GlueXSensitiveDetectorFCAL.cc:
int pmass = track->GetDynamicParticle()->GetMass();
This raises the question, how does a Cerenkov photon get inside the FCAL? This is a Cerenkov detector, but AFAIK there is no explicit Cerenkov generation implemented for the FCAL in this simulation. The run number is 71000 so there is a DIRC in the geometry, but it should not leak like this unless there is an overlap problem somewhere. Also, why does it only happen for certain kinds of particles in the particle gun sim samples?

  • Lambda
  • Klong
  • Kshort
  • eta
  • and more....

but not

  • proton
  • antiproton
  • pi0
  • neutron
  • and more...

@jrstevenjlab
Copy link

jrstevenjlab commented Mar 21, 2021

I suppose one possibility if you see this in the FCAL is the Plexiglas material. At one point for the DIRC we were considering a plexiglass window instead of fused silica and added optical properties to it. I see now that they are still in the material definition which is used for the LED monitoring panels in front of the FCAL but not used anywhere in the DIRC

https://github.com/JeffersonLab/hdds/blob/de899d6a375b3ef696eaecb5c2f5b119b757b54b/Material_HDDS.xml#L679

If you remove the optical properties, it you shouldn't get Cherenkov photons even if the CKOV option is on, right? Why it's different for certain kinds of particles... I have no idea.

@rjones30
Copy link
Author

rjones30 commented Mar 22, 2021 via email

@markdalton
Copy link

The FCAL LED monitoring system is directly upstream of the FCAL blocks. The are 4 large sheets which each cover a single quadrant of the detector. The LEDs themselves are at the edge of these sheets and the sheets act to transport the light from the edge to all of the blocks.

@rjones30
Copy link
Author

rjones30 commented Sep 7, 2021

Cerenkov light can be produced in the light guides at the back of the fcal. This is normal, but explicit following of CKOV photons should be suppressed in the light guides. I have not seen this warning in a while, closing this issue until it pops up again. -Richard

@rjones30 rjones30 closed this as completed Sep 7, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants