#include "Garfield/MediumMagboltz.hh"
#include "Garfield/FundamentalConstants.hh"

using namespace Garfield;

int main(int argc, char * argv[]) {

  const double pressure = 0.562; // 0.75 mbar
  const double temperature = 313.;
 
  // Setup the gas.
  MediumMagboltz gas("C3H8");
  gas.SetTemperature(temperature);
  gas.SetPressure(pressure);
  // Set the field range to be covered by the gas table. 
  const size_t nE = 20;
  const double emin =    1.;
  const double emax = 5000.;
  // Flag to request logarithmic spacing.
  constexpr bool useLog = true;
  gas.SetFieldGrid(emin, emax, nE, useLog); 

  const int ncoll = 20;
  // Run Magboltz to generate the gas table.
  gas.GenerateGasTable(ncoll);
  // Save the table. 
  gas.WriteGasFile("C3H8_20_1_5000_313K_0.75mbar.gas");

}