v7.4.1 Feb, 7 2023
  1) Fixed bug in Cross Sec for PEarth, now running for VLE tune event generation lists
  2) Changes in container to fixed problems genie: verbosity reduced, IBD list for VLE tune added, push release docker image in CI 
  3) Linked to km3net-dataformat v3.0.0

v7.4.0 Oct, 5 2022
  1) New coordinate systems with origin fixed on the seabed
  2) Possibility to move the detector can in water along the Z axis
  3) Removed link to old versions of PROPOSAL
  4) Activated step-by-step propagation in water also for PROPOSAL and set vcut_infront = 0.05
  5) Fixed bug in definition of PROPOSAL geometry and sector affecting simulations with origin not on the seabed
  6) GENIE monitor: file name with same prefix of other outputs and refresh every 1000 events
  7) Added polarization in the track class
  8) Added option to choose the output format and removed old options
  9) Added random initialization in CreaBin

v7.3.0 Apr, 19 2022
  1) The status signs have been changed to positive, in accordance with km3net-dataformat
  2) Treatment of the Earth's curvature for CORSIKA has been revised, now all showers are rotated around the Earth instead of shifted horizontally in xy. This means that both (pos_x,pos_y,pos_z) and (dir_x,dir_y,dir_z) are modified accordingly for each track (previously: only pos_x,pos_y).
  3) The position of the (pos_x,pos_y,pos_z) is now consistently at the 1st interaction point
  4) CORSIKA CURVOUT option is now correctly accounted for
  5) Some of the code structure has been factorized out into separate functions to improve readablility and testability
  6) Bug in treatment of the MUPROD option from CORSIKA has been corrected
  7) Separate directory holding helper function libraries `utils` has been added
  8) The code is now linkable to GENIE v3.2.0
  9) Minor bug fixes and aesthetic improvements
  10) Addded doxygen documentation of the classes to https://simulation.pages.km3net.de/gseagen/gseagen.pdf
  11) Updated to km3net-dataformat v2.2.0 (now tracks have counter attribute to store CORSIKA information)

v7.2.1 Jan, 31 2022
  1) Pscale is now multiplied by genie weight
  2) CORSIKA track ordering when the muon track has sufficient and insufficient range has been unified (relevant to CORSIKA only)
  3) Corrected the usage of PropagationTolerance
  4) Fixed the bug with negative z positions for parent tracks (relevant to CORSIKA only)
  5) Fixed the bug with shifting of OtherTracks (relevant to CORSIKA only)
  5) Fixed the bug with duplicate track IDs (relevant to CORSIKA only)
  6) Status 200 is now used (instead of 0) for the primary when using corant input (unified with direct CORSIKA readout)
  7) Removed the unnecessary assert, which was causing problems for tau propagation 

v7.2.0 Dec, 20 2021
  1) Added the option to reduce CORSIKA verbosity
  2) Changed proposal geometry to spherical, fixed bug in propagation in rock
  3) Propagation in water is now done in two steps: first till the can top height, then until the can is reached/missed. 
  4) ROOT functions have been dropped in favour of the std library
  5) Code has been optimized for speed and efficiency: aside from the previously available option to just retry showers from start, they can now also be retried by reshifting them at the can top height (saving time on unnecessary propagation)
  6) -rt "proj" option now samples muons on the curved (taking into account Earth curvature) sea surface (before: volume above and below this surface)
  7) fDistaMax was not saved in the event and thus included in Agen calculation for CORSIKA processing, now it is
  8) Corrected the values passed to the fDistaMax formula (primary vertex is now properly translated back up from the sea surface to the first interaction point)
  9) dtypes were reduced for PropaTrack and gSeaTrack
  10) The main functionality of KM3NeTWrite is now in the .cxx, not in the header file
  11) The sea-level primary position is used for shifting when processing CORSIKA
  12) Only every 100000-th event is displayed for SilentCORSIKA (useful for TeV prod)
  13) Fixed the MUSIC bug by catching the muons "going back in time" (going upwards and gaining energy)
  14) The size of the cross-sections file has been reduced to a minimum to decrease the size of the singularity image
  15) Updated to newer km3net-dataformat
  16) More information about CORSIKA settings has been added to the header
  17) The documentation in the git wiki has been extended by the description of the propagation
  18) For PROPOSAL a problem with the final muons has been fixed (sometimes not only muons were picked as final particles)
  19) Unit conversion for PROPOSAL tracks has been corrected and now the track time values are correct
  20) A bug present since v7.0.0 has been fixed: there were duplicate muons because of incorrect CORSIKA file readout
  21) Fixed a bug that allowed muons at the can to have status -1001 or -2 (intermediate value introduced at the can top height, which should never appear in the output files)
  22) Lateral deflection of the muons has been re-evaluated. Before simply a fixed value of 100m was used, regardless of muon energy. Now a fitted function is used.
  23) Fixed bug in neutrino costheta
  24) Track lenght calculation has been corrected
  25) In "proj" option now projected values of DistaMax are used to increase the can (previously: full value for each dimension)
  26) The propagation tolerance was set to 5m. The 5m tolerance applies in 1D and for slant distances a projected tolerance is used.
  27) Native output for CORSIKA bug is now fixed
  28) Corrected unnecessary unit conversion bug in PropagationMu()
  29) If option -write 0 is used, the mother id's of the secondaries will be set to the primary (as no other parent will be stored) for CORSIKA processing. For GENIE, the mother id's will be set to grandmother id's instead.
  30) added --muon-range-tolerance option that allows to try and propagate muons, which have low chance of reaching the can (may still be useful for low-E productions)
  31) For PROPOSAL each propagation stage is done in one step. For other propagators it starts with a big initial step, then proceeds with smaller steps.
  32) Changed the mother_id to have a consistent file when not all particles are saved (-write 0)
  33) Changed default for -wgt option to 1


v7.0.3 Aug, 12 2021
  1) Changed the default maximal number of records read from CORSIKA files to 0
  2) Corrected the nuclear mass formula in GSeaCORSIKAFileFlux::CorsikaRestMass to include nuclei with odd mass numbers, e.g. Helium-3

v7.0.2 May, 4 2021
  1) Using track event type interaction volume if Glashow resonace is simulated

v7.0.1 Mar, 14 2021
  1) Include the more correct fDistaMax computation in GSeaAntaresFileFlux driver 

v7.0.0 Mar, 10 2021
  1) XSecMean in the third entry of w2list replaced with XSec
  2) Linkable with PROPOSAL 5.0.0 and later
  3) Linkable with JPP/master
  4) start_run and mc_run_id inlcuded in km3net-dataformat output
  5) tauola libraries added to setenv.sh
  6) Fix bug filling fourth field of 'genvol' in km3net-dataformat output
  7) Put back Bx,By,InterId,Scattid in "neutrino" tag (remove info from "w2list") and "target" and "vertex" tags for .evt files.
  8) Reduce verbosity of of PDG codes when filing .evt files.
  9) Add lenght in "track_*" tag.
 10) Introduce the "track_aux" tag to store particles that dont produce light in .evt files.
 11) Use w2list_gseagen definitions.
 12) Remove InCan variable. Now everything is handled with Status.
 13) Neutrinos are not loaded when using corant files.
 14) Mother and grandmother added as "fake" particles when using corant files.
 15) Remove 'can' and 'coord_origin' from headers.
 16) Added -time running opt to define generation time interval (no longer defined in the astro file)
 17) GAstro lib: implemented precession and nutation and improved Gal-Eq transformations
 18) Clean up configure. Now boost, tbb and cern libraries are also detected.
 19) Add event time to km3net-dataformat (valid from v1.3.1) output.
 20) Added tool gCalcAstro to check astronomical coordinates
 21) Generation weights are always per seconds
 22) Added GSeaCORSIKAFileFlux flux drive to read directly CORSIKA files
 23) Improved calculation of the convergence angle and fixed bug in the orientation of the UTM system axes

v6r1 Nov, 29 2019
 1) Added HONDA fluxes
 2) Flux Drivers: code re-organization
 3) Added configure script
 4) setenv.csh and setenv.sh not in the rep but created by configure
 5) Linkable with genie version 3.XX.XX
 6) files included in the MAkefile are in ser/make directory
 7) dir dat: no symbolic links to xsec
 8) GGenerateEventGenie: CalcXSecWater and CalcXsecTot compute xsec in the global valid energy range (and not in the using one)
 9) gCalcCross: check if using Emin and Emax are in the valid energy range; comp xsec are now per nucleon
10) Added the possibility to throw neutrino using the projected can
11) Removed limit in the flavor generation
12) Removed solid angle in the extended source weights
13) Changed default name for xsec file
14) Fixed bug in detector coordinates input
15) Default flux changed from 1E-9*pow(x,-2) into 1E-5*pow(x,-2)

v5r2 Jul, 5 2018
 1) PropaMuon: new model for multiple scattering, improved check to stop propagation at the required 
    depth taking into account large deviations due to multiple scattering and bug fixed for Emu<Ecut 
 2) New DistToCan and Propagation methods in GGenerateEvent: prevent infinite loop
 3) New option -rt to set RT radius
 4) Removed particles at rest from the output
 5) Fixed bug in option -l: worked only if input before option -c
 6) New option -l: now it is possible to change NAbsLenght
 7) Fixed bud in class GSeaEvent: TargetZ and TargetA defined but not filled 
 8) Linked to PROPOSAL
 9) Propagators initialized only if used
 10) Linked to NuWro, added interface class for NuWro (not yet implemented)
 11) Code reorganization: xsec entering in PEarth are non calulated by PEarth class but the generator interface
 12) New option -all to write all interacting neutrinos
 13) New option -origin the origin position coordinates (in meters) in the detector rest frame
 14) Origin position set on seabed when can set from a genometry file
 15) Changed TopVolume to fix bug when can1=0 
 16) Generation area is centered on the can (-rt = "can" and -rt = value) or on the interaction vol (-rt = "vol").
 17) RT is the semidiagonal of the can (or the interaction vol) and does not depend on the origin position
 18) GPEarth changes the rock radius and not the ocean one in order to have a water layer equal to SiteDepth as suggested by Alfonso Garcia
 19) PEarth: calculated using only CC and not Default EventGeneratorList

v5r1 Jan 16, 2018
 1) Fixed bug in calculation of track_in direction cosines
 2) Fixed bug in writing output root file
 3) Added new output root file format (simple tree)
 4) Help not printed in case of fatal error in input parameters. Users had difficulties in finding error msg
 5) New RT value for nu_mu in CC
 6) Part of the code depending on the antares libs are compiled if ANTARES env does not exist but KM3NET_THRONG_DIR exists
 7) .detx geometry files are read if ANTARES or KM3NET_THRONG_DIR env exists  
 8) Added GGenerateEvent class which run the GENIE standard evt generation driver. In the main added check to prevent infinite loop
 9) Propagation through the Earth working but activated only for checks 
 10) Linkable  with genie 2.12.X 
 11) Removed dependence on genie version (no config and xsec files)
 12) Added tools (gCalcCross)
 13) Fixed bug in interaction volume calculation in evt file

v4r2 Apr 24, 2017 
 1) Added the possibility to save particles generated in the bedrock immediately below the seabed
 2) Fixed bug reading declination with DEC: option
 3) Fixed bug reading -l option
 4) Fixed bug in the computation of xsec entering in the PEarth
 5) Added cross section and interaction length in pure water in output files
 6) Library to simulate muon propagation (now called PropaMuon) is stand-alone
 7) Linked with genie 2.10.10 and added new xsec splines and config files
 8) Implemented simulation of neutrino propagation through the Earth (in progress)

v4r1 May 11, 2016
 1) Part of the code depending on the antares libs is not compiled is ANTARES env does not exist
 2) Part of the code using MUISC compiled only if MUSICSRCPATH env exists
 3) Output Root file implemented
 4) Implemented library to simulate muon propagation
 5) Possibility to simulate point (and extended) sources

v3r3 January 10, 2016
Improvements over v3r2:
 1) Linked with genie 2.10.2 and added new xsec splines 
 2) Added the possibility to change the density of the target media. Input ignored for Mantle and Core. 

v3r2 June 3, 2015
Improvements over v3r1:

 1) Event weights computed only if the neutrino produces a detectable event (before calculated for each generated neutrino) 
 2) Added default for MUSICPATH = $GSEAGEN/dat and  GMURNGFL = $GSEAGEN/dat/muon_rmax_music.dat
 3) Added default for GALGCONF = $GSEAGEN/config
 4) Added symbolic links to cross section files (main Makefile)
 5) Added default for XSecName = $GSEAGEN/dat/gxpl-seawater.xml 
 6) Added elements in Rock composition, added Mantle and Core composition 
 7) The user can re-define SeaWater, Rock, Mantle and Core composition from a xml file (env parameter MEDIACOMP)
 8) Column depth calculated for the true neutrino track (before for a track crossing the origin, assumed as vertex) 
 9) PEarth calculated taking into accaunt the composition of the Earth layers (Core, Mantle, Rock and SeaWater)
10) Cross section in w2list is the mean cross section along in the neutrino path: PEarth=exp(-Na*XSecMean*ColumnDepth)
11) Implemented evaluation of systemayics (neutrino interaction)
12) Possibility to set ctmin=ctmax
13) Possibility to set Emin=Emax, but in this case the code set Emax=Emax+0.001 

v3r1 Febrary 6, 2015 
Improvements over v2r1:

1) Functions PropMusic and DistToCan have been written. Removed link to genhen and geasim.
2) Possibility to simulate tau (anti)neutrinos. Generation volume: can if NC interaction; can + muon range if CC or CC+NC interaction.
3) New run option: -C. The generation volume is forced to be the can (also in case of muon and tau nu CC interaction).
4) News in run option -f: possibility to define the real neutrino flux as a math formula. Default flux added (power law with alpha=-2).
5) Muons not propagated if DistToCan greater than the maximum range in water

Linked with genie2.8.6, cross section file updated.