-
Notifications
You must be signed in to change notification settings - Fork 44
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
missing MPGD hits #774
Comments
@ShujieL Can you post some more diagnostics here since you already started looking into it? |
@wdconinc, With my 20k 10GeV pi+ sample at 70<theta<110 degrees, I got the following number of events with at least one hit at given radius range. The events missed MPGD seem to have no preference in theta and phi. Hit collection name / rmin / rmax [mm] / # of events |
Can you rerun one specific event with the track reconstruction in debug log level and analyze whether the hit even gets turned into a measurement? We know that for slanted tracks through the 3mm thick MPGD layers there is a non-negligible chance that there is a geant4 step that ends in the middle of the gas volume. That leads to two tracker hits then and each is going to be a bit displaced from the actual location of the center plane crossing. Can you rule that out based on debug output in the tracker to measurement converter? |
those numbers are from RecHits which are converted from sim hits after digitization, before the tracker measurement conversion. Also, MPGDBarrelRawHits and MPGDBarrelRecHits have the same number of entries. So the hits are very likely missing in sim hits (geant4). |
Interesting. Please post more as you find out more. |
I will check the impact of missing MPGD hits in tracking. Meanwhile, I'd suggest
|
@ShujieL Thanks for already starting to look into this. I can look into the acceptance a bit more and I will also bring this up during today's MPGD meeting. One thing to keep in mind is that we know there are acceptance gaps in the outer MPGD barrel by design (follows DIRC bar geometry), so having less Sim hits in the outer MPGD barrel compared to the inner MPGD barrel should be expected. |
It seems that the hits are missing from the sim hits to raw hits transition, i.e. digitization. See https://indico.bnl.gov/event/24735/contributions/96326/attachments/57062/97925/missing%20hits.pdf |
Because the SiliconTrackerDigi is used for the MPGDs, the only reasons why sim hits may not turn into sim hits are:
|
With a threshold of 0.25 keV, that is the main reason for cutting MPGD hits here. Below is a plot of the energy deposition for simulated hits, in x-axis units of GeV, so keV if you ignore the 10^-6. A cut of 0.25 keV cuts a sizeable fraction of these events. @mposik1983 Are these thresholds correct? |
@wdconinc for now the threshold values should stay as they are. From the discussion at the tracking meeting, there should be some follow up next week about a new threshold estimate |
Yeah, I didn't know you were going to present on this exact issue today :-) |
Environment: (where does this bug occur, have you tried other environments)
I generated some 10 GeV pi+ events at theta between 70 to 110 degrees with epic 24.08.1, and noticed that more than 10% events has no hits on inner MPGD.
I also checked a simulation campaign file and saw similar behaviors
EPIC/RECO/24.07.0/epic_craterlake/SINGLE/pi+/10GeV/45to135deg/pi+_10GeV_45to135deg.0000.eicrecon.tree.edm4eic.root'
Steps to reproduce: (give a step by step account of how to trigger the bug)
Expected Result: (what do you expect when you execute the steps above)
At 10 GeV and central angle, almost every single pi+ particle should trigger inner MPGD.
Actual Result: (what do you get when you execute the steps above)
Less than 90% events have inner MPGD hits. This number is even lower for outer MPGD
The text was updated successfully, but these errors were encountered: