Support #11998
Provide a better way to simulate non-standard physics inside the nucleus
0%
Description
The "non-standard" in the subject of the issue refers to all processes that are not natively implemented by the supported generators (focus on GENIE).
This issue arises from the intention for MicroBooNE to simulate the phenomenology of a new process (boosted dark matter) elaborated by a theorist.While this specific process may be worked around pretty easily, since it turns into just a single proton that can be simulated with
SingleGen
and it is not unreasonable to neglect nucleus interactions, other processes don't have that luxury.For example, the procedure actually used by Jeremy Hewes for neutron-antineutron oscillation events is incredibly cumbersome.
The simulated bits are decay products from a
The procedure Jeremy follows is:
- run a custom GENIE-based executable outputting a ROOT file
- run a GENIE executable dumping the ROOT format into a text file
- run LArSoft with a source module1 parsing that text file to produce a HEPEVT file
Possible better solutions may be:
- have a LArSoft module similar to
GENIEGen
that can run an arbitrary GENIE executable and perform the output conversion as needed - have a class like nutools's
GENIEHelper
, with the additional flexibility to allow GENIE commands to be executed - have a LArSoft module able to read a HEPEVT event and run GENIE's FSI on it:
- read a HEPEVT file
- fill GENIE with that information
- let GENIE run FSI
- collect the results in the same way as the standard
GENIEGen
module does
Note that solution 3 does not answer Jeremy Hewes's use case since GENIE is also used to define the initial state based on nucleus characteristics.
1 That module is fairly unsupported and required some adaption (performed by MicroBooNE) to work for Jeremy's use case.
History
#1 Updated by Robert Hatcher almost 5 years ago
Note: My (R. Hatcher) near term plan -- to satisfy immediate DUNE Near Det Task Force needs -- involves writing a module to read standard GENIE gntp.*.ghep.root files which would in this case cut out the whole "squishing" into a text file (and thus losing some auxiliary information not represented in a HEPEVT format).
#2 Updated by Katherine Lato over 4 years ago
- Assignee set to Robert Hatcher
Via email, from Robert, that module is already written. It currently lives in the DUNE NearDetector Task Force repository, but the plan is to migrate it to NuTools.
#3 Updated by Katherine Lato about 4 years ago
- Status changed from New to Assigned