Comments (7)
Hi Diego, if I understand this warning right, it's a warning that you've also propagated from NEST itself for where we are outside the fit-regime of validity with NRs. (By the way we can probably increase this to oh, about 330 keV now with newer AmBe data I believe. Or so the NEST docs tell me!)
Can you explain how the microclustering will affect this particular line? Because this warning is thrown based on only input arguments (energy and interaction only.) So how will you envision the MicroSeparation will impact this error?
from epix.
To clarify, it seems you're generating quanta internally from some energy you fill in. Can you link here where in the code that is so I Can take a look at the related step prior to this error?
from epix.
Hi Sophia, the step for quanta generation is indeed just calling NEST with some energy and drift field value for a given cluster. It happens in this file, while the actual warning I quoted is here and relies on the NEST limits that we just hardcoded (maybe there's a better way to do this....).
The relation with MicroSeparation
arises because the clustering step happens before we fetch the efield at the interaction site, call NEST, etc. MicroSeparation
decides if two energy deposits were close enough to tag them as a single one and, once this is done, this function takes care of merging the actual interactions into a cluster. The resulting vertex has weighted position according to the energy of the merged steps, recoil type according to time/energy, and energy deposit as the sum of the individual interactions. If MicroSeparation
is too big, we will merge more events into a single cluster and the resulting total energy will be larger, therefore higher chances to reach the NEST limits afterwards.
from epix.
Hi Diego, thanks for the explanation. I think this makes sense to me.
The only thing to be careful of is that just because it hits NEST limits does not mean that the energy is unphysical... The cryoneutrons root file you provided: does that have a max cutoff energy or is it monoenergetic? (I'm guessing since it is called tpc_cryoneutrons_200.root
you have some 200 keV energy cutoff or something but let me know if I'm wrong.)
from epix.
Quick reply on the second thing: sorry, the filename is a bit stupid since this was just a small sample test. The 200 stands for the number of simulated primaries (neutrons). The initial energy of these neutrons (with starting position confined to the inner and outer cryostat volumes) follows a flat distribution from 0 to 10 MeV.
(Update: tpc
stands for the fact that only LXe records entries in this simulation, and not PMTs or other sensitive volumes -- please ignore).
from epix.
I see. So really we could physically have an energy deposition of >300 keV, given the initial energy spectrum. I don't know of any constraints there but maybe the kinematics don't allow it. The only constraint I'm aware of is that we don't have data fit in this region.. meaning, you want to split up the events into smaller events where we have data (and then I assume sum these quanta back together into a final "event".) Let me know if I'm still missing something.
from epix.
This is indeed the question. We wonder if our MicroSeparation
is too big, since we get such warnings so often, or, on the contrary, it is fine and we need to accept that many interactions are beyond the NEST fitted data; but in that case we should decide how to handle these interactions (removing them sounds messy).
meaning, you want to split up the events into smaller events where we have data (and then I assume sum these quanta back together into a final "event".) Let me know if I'm still missing something.
Since summing the quanta from interactions with energy E1 and E2 would differ from the quanta obtained for an interaction with E=E1+E2, we keep trying to find a reasonable (~ NEST-being-happy-with-it) microseparation before quanta generation.
I don't really have an asnwer, but this issue is just part of the recoil classification and clustering review that I am trying to conduct with @HenningSE (and you, if interested!), who looked compared more in-depth the NEST methods with ours. Do you think this would deserve a call?
from epix.
Related Issues (20)
- Chunk for large files
- Add excitons to instructions HOT 1
- Classification of gamma interactions
- Epix is untested
- maxTime in NEST Kr83m model HOT 4
- Keep track of removed interactions/events
- fast_simulator - Macro-clustering issue HOT 1
- Significant fraction of time is just spent in loading text objects
- EPIX doesn't make events except initial one after v0.1.2 HOT 2
- Problem with the cluster indexing after DBSCAN HOT 3
- Incorrect input parameters for NEST HOT 3
- Slower code after fixing clustering index issue (#45)
- Kr83m deltaT_ns
- `culstering._write_result` fails in numba `v0.56.0` HOT 1
- Unify quanta generation classes
- Calculation of average time per cluster can lead to numerical errors HOT 1
- Numba - Awkward Array Version Problem HOT 2
- Alpha events are always misclassified as NR
- module 'awkward' has no attribute 'to_pandas' since awkward version 2.0 HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from epix.