Charged fusion product diagnostic forward model and velocity-space weight function

[alvin-garcia]

FIDASIM community,

It is my pleasure to announce the implementation of CFPD signals and weight functions into FIDASIM.

IDL procedures run_tests.pro and test_cfpd.pro in FIDASIM/test provide an example of how to prepare the inputs. The user must insert an "orbits table" into *_geometry.h5 and set calc_cfpd = 1 to output the count rate and weight function in the output file *_cfpd.h5. The table process happens in the back end of test_cfpd.pro and dependencies are located in FIDASIM/lib/idl. In order to prepare the orbits table using our IDL code, the user must supply a GEQDSK file and CFPD diagnostic geometry. Also, you can use plot_cfpd_orbits.pro to visualize sightline. Note: A new CFPD tab will be inserted into our documentation website with better instructions soon.

Figure 1: Orbits produced by FIDASIM test cases

Figure 2: 3 MeV proton signal produced by test cases

Below is a recent paper and poster with more details related to this PR:

1. Heidbrink, W. W., Garcia, A. V., et al. "Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics." Plasma Physics and Controlled Fusion (accepted)
2. Development of weight functions for 3 MeV proton diagnostics in FIDASIM

@lstagner ,

Please review this branch and note the following

• Implementing Python equivalents of IDL prefida code is deferred to a later date, so don't worry about that for now.
• docs will be updated in a forthcoming PR
• The NPA-like definitions of the CFPD geometry are not really used by the main program in FIDASIM. However, I threw them in anyways in case they are needed in the future.

[alvin-garcia]

Hi @lstagner ,

Please merge this PR if you are done reviewing. I'm still building and testing the Python prefida code and will include that plus documentation in a forthcoming PR. Thanks.