303 lines
8.1 KiB
C
303 lines
8.1 KiB
C
#ifndef _H5UNITS_H_
|
|
#define _H5UNITS_H_
|
|
/*!
|
|
Defines Units and conversion factors
|
|
|
|
\note Derived from GEANT 4 header files G4SIunits.hh
|
|
|
|
The basic units are those of the International System:
|
|
|
|
meter
|
|
second
|
|
kilogram
|
|
ampere
|
|
degree kelvin
|
|
the amount of substance (mole)
|
|
luminous intensity (candela)
|
|
radian
|
|
steradian
|
|
|
|
We also define some accelerator specific units
|
|
|
|
MeV
|
|
beta
|
|
|
|
Idea of usage: For each attribute (for example x) we have two additional
|
|
informations, namely the attributes unit and the attributes conversion factor,
|
|
which is the factor to convert the attributes unit into SI units.
|
|
|
|
Example:
|
|
|
|
We have a lenght attribute x and the values in x are stored in mm we then would
|
|
have to add the following information to the attribute x: METER, 1000
|
|
|
|
Using the following prototypes
|
|
|
|
H5PartWriteUnit (H5PartFile *f, char *name, const int unit );
|
|
H5PartWriteUnitConversionFac (H5PartFile *f, char *name, const double fac);
|
|
|
|
one would write
|
|
|
|
H5PartWriteUnit (fp,"x", METER );
|
|
H5PartWriteUnitConversionFac (fp, "x", 1000);
|
|
|
|
Note: in case of BETA and GAMMA one has the know the particles mass
|
|
not sure if that is the best solution
|
|
|
|
Author: Andreas Adelmann (PSI)
|
|
Date: 15.8.2006
|
|
|
|
*/
|
|
|
|
// Define SI units
|
|
// \note enum type would be much better, have to check f77 implications
|
|
|
|
#define METER 1
|
|
#define SECOND 2
|
|
#define KILOGRAM 3
|
|
#define AMPERE 4
|
|
#define DEGREE_KELVIN 5
|
|
#define MOLE 6
|
|
#define CANDELA 7
|
|
#define RADIAN 8
|
|
#define STERADIAN 9
|
|
|
|
// Define some accelerator specific units
|
|
|
|
#define MEV 10
|
|
#define BETA 11
|
|
#define GAMMA 12
|
|
|
|
//
|
|
// Length [L]
|
|
//
|
|
static const double meter = 1.;
|
|
static const double meter2 = meter*meter;
|
|
static const double meter3 = meter*meter*meter;
|
|
|
|
static const double millimeter = 0.001*meter;
|
|
static const double millimeter2 = millimeter*millimeter;
|
|
static const double millimeter3 = millimeter*millimeter*millimeter;
|
|
|
|
static const double centimeter = 10.*millimeter;
|
|
static const double centimeter2 = centimeter*centimeter;
|
|
static const double centimeter3 = centimeter*centimeter*centimeter;
|
|
|
|
static const double kilometer = 1000.*meter;
|
|
static const double kilometer2 = kilometer*kilometer;
|
|
static const double kilometer3 = kilometer*kilometer*kilometer;
|
|
|
|
static const double parsec = 3.0856775807e+16*meter;
|
|
|
|
static const double micrometer = 1.e-6 *meter;
|
|
static const double nanometer = 1.e-9 *meter;
|
|
static const double angstrom = 1.e-10*meter;
|
|
static const double fermi = 1.e-15*meter;
|
|
|
|
static const double barn = 1.e-28*meter2;
|
|
static const double millibarn = 1.e-3 *barn;
|
|
static const double microbarn = 1.e-6 *barn;
|
|
static const double nanobarn = 1.e-9 *barn;
|
|
static const double picobarn = 1.e-12*barn;
|
|
|
|
// symbols
|
|
static const double mm = millimeter;
|
|
static const double mm2 = millimeter2;
|
|
static const double mm3 = millimeter3;
|
|
|
|
static const double cm = centimeter;
|
|
static const double cm2 = centimeter2;
|
|
static const double cm3 = centimeter3;
|
|
|
|
static const double m = meter;
|
|
static const double m2 = meter2;
|
|
static const double m3 = meter3;
|
|
|
|
static const double km = kilometer;
|
|
static const double km2 = kilometer2;
|
|
static const double km3 = kilometer3;
|
|
|
|
static const double pc = parsec;
|
|
|
|
//
|
|
// Angle
|
|
//
|
|
static const double radian = 1.;
|
|
static const double milliradian = 1.e-3*radian;
|
|
static const double degree = (3.14159265358979323846/180.0)*radian;
|
|
|
|
static const double steradian = 1.;
|
|
|
|
// symbols
|
|
static const double rad = radian;
|
|
static const double mrad = milliradian;
|
|
static const double sr = steradian;
|
|
static const double deg = degree;
|
|
|
|
//
|
|
// Time [T]
|
|
//
|
|
static const double second = 1.;
|
|
static const double nanosecond = 1.e-9 *second;
|
|
static const double millisecond = 1.e-3 *second;
|
|
static const double microsecond = 1.e-6 *second;
|
|
static const double picosecond = 1.e-12*second;
|
|
|
|
static const double hertz = 1./second;
|
|
static const double kilohertz = 1.e+3*hertz;
|
|
static const double megahertz = 1.e+6*hertz;
|
|
|
|
// symbols
|
|
static const double ns = nanosecond;
|
|
static const double s = second;
|
|
static const double ms = millisecond;
|
|
|
|
//
|
|
// Mass [E][T^2][L^-2]
|
|
//
|
|
static const double kilogram = 1.;
|
|
static const double gram = 1.e-3*kilogram;
|
|
static const double milligram = 1.e-3*gram;
|
|
|
|
// symbols
|
|
static const double kg = kilogram;
|
|
static const double g = gram;
|
|
static const double mg = milligram;
|
|
|
|
//
|
|
// Electric current [Q][T^-1]
|
|
//
|
|
static const double ampere = 1.;
|
|
static const double milliampere = 1.e-3*ampere;
|
|
static const double microampere = 1.e-6*ampere;
|
|
static const double nanoampere = 1.e-9*ampere;
|
|
|
|
//
|
|
// Electric charge [Q]
|
|
//
|
|
static const double coulomb = ampere*second;
|
|
static const double e_SI = 1.60217733e-19; // positron charge in coulomb
|
|
static const double eplus = e_SI*coulomb ; // positron charge
|
|
|
|
//
|
|
// Energy [E]
|
|
//
|
|
static const double joule = kg*m*m/(s*s);
|
|
|
|
static const double electronvolt = e_SI*joule;
|
|
static const double kiloelectronvolt = 1.e+3*electronvolt;
|
|
static const double megaelectronvolt = 1.e+6*electronvolt;
|
|
static const double gigaelectronvolt = 1.e+9*electronvolt;
|
|
static const double teraelectronvolt = 1.e+12*electronvolt;
|
|
static const double petaelectronvolt = 1.e+15*electronvolt;
|
|
|
|
// symbols
|
|
static const double MeV = megaelectronvolt;
|
|
static const double eV = electronvolt;
|
|
static const double keV = kiloelectronvolt;
|
|
static const double GeV = gigaelectronvolt;
|
|
static const double TeV = teraelectronvolt;
|
|
static const double PeV = petaelectronvolt;
|
|
|
|
//
|
|
// Power [E][T^-1]
|
|
//
|
|
static const double watt = joule/second; // watt = 6.24150 e+3 * MeV/ns
|
|
|
|
//
|
|
// Force [E][L^-1]
|
|
//
|
|
static const double newton = joule/meter; // newton = 6.24150 e+9 * MeV/mm
|
|
|
|
//
|
|
// Pressure [E][L^-3]
|
|
//
|
|
#define pascal hep_pascal // a trick to avoid warnings
|
|
static const double hep_pascal = newton/m2; // pascal = 6.24150 e+3 * MeV/mm3
|
|
static const double bar = 100000*pascal; // bar = 6.24150 e+8 * MeV/mm3
|
|
static const double atmosphere = 101325*pascal; // atm = 6.32420 e+8 * MeV/mm3
|
|
|
|
//
|
|
// Electric potential [E][Q^-1]
|
|
//
|
|
static const double megavolt = megaelectronvolt/eplus;
|
|
static const double kilovolt = 1.e-3*megavolt;
|
|
static const double volt = 1.e-6*megavolt;
|
|
|
|
//
|
|
// Electric resistance [E][T][Q^-2]
|
|
//
|
|
static const double ohm = volt/ampere; // ohm = 1.60217e-16*(MeV/eplus)/(eplus/ns)
|
|
|
|
//
|
|
// Electric capacitance [Q^2][E^-1]
|
|
//
|
|
static const double farad = coulomb/volt; // farad = 6.24150e+24 * eplus/Megavolt
|
|
static const double millifarad = 1.e-3*farad;
|
|
static const double microfarad = 1.e-6*farad;
|
|
static const double nanofarad = 1.e-9*farad;
|
|
static const double picofarad = 1.e-12*farad;
|
|
|
|
//
|
|
// Magnetic Flux [T][E][Q^-1]
|
|
//
|
|
static const double weber = volt*second; // weber = 1000*megavolt*ns
|
|
|
|
//
|
|
// Magnetic Field [T][E][Q^-1][L^-2]
|
|
//
|
|
static const double tesla = volt*second/meter2; // tesla =0.001*megavolt*ns/mm2
|
|
|
|
static const double gauss = 1.e-4*tesla;
|
|
static const double kilogauss = 1.e-1*tesla;
|
|
|
|
//
|
|
// Inductance [T^2][E][Q^-2]
|
|
//
|
|
static const double henry = weber/ampere; // henry = 1.60217e-7*MeV*(ns/eplus)**2
|
|
|
|
//
|
|
// Temperature
|
|
//
|
|
static const double kelvin = 1.;
|
|
|
|
//
|
|
// Amount of substance
|
|
//
|
|
static const double mole = 1.;
|
|
|
|
//
|
|
// Activity [T^-1]
|
|
//
|
|
static const double becquerel = 1./second ;
|
|
static const double curie = 3.7e+10 * becquerel;
|
|
|
|
//
|
|
// Absorbed dose [L^2][T^-2]
|
|
//
|
|
static const double gray = joule/kilogram ;
|
|
|
|
//
|
|
// Luminous intensity [I]
|
|
//
|
|
static const double candela = 1.;
|
|
|
|
//
|
|
// Luminous flux [I]
|
|
//
|
|
static const double lumen = candela*steradian;
|
|
|
|
//
|
|
// Illuminance [I][L^-2]
|
|
//
|
|
static const double lux = lumen/meter2;
|
|
|
|
//
|
|
// Miscellaneous
|
|
//
|
|
static const double perCent = 0.01 ;
|
|
static const double perThousand = 0.001;
|
|
static const double perMillion = 0.000001;
|
|
|
|
|
|
#endif |