float process(float Beta, std::ifstream & INPUT, float mass, float energy, int Juse = 1, int Jerror = 1, int intype = 1) { float crosscons = 8 * M_PI * mass * ma