afsi.h

#ifndef AFSI_H
#define AFSI_H
 
#include <math.h>
#include "ascot5.h"
#include "consts.h"
#include "simulate.h"
#include "random.h"
#include "boschhale.h"
#include "diag/hist.h"

typedef enum {
    PPARPPERP,
    EKINXI
} mom_space_basis;

typedef struct {
    int type1;          
    int type2;          
    int thermal1;       
    int thermal2;       
    histogram* beam1;   
    histogram* beam2;   
    real* r;            
    real* phi;          
    real* z;            
    real* vol;          
    size_t volshape[3]; 
    Reaction reaction;  
    real mult;          
} afsi_data;
 
void afsi_run(sim_data* sim, afsi_data* data, int n,
              histogram* prod1, histogram* prod2);
#endif