HISTORY of WimpEvent ==================== -------- Revision 337, 2018-06-19 * Added configure system (autoconf) to make it easier to compile the code. -------- Revisopm 316, 2017-04-07 * Added possibility to read event files with weights from both solar_crnu and WimpAnn. -------- Revision 129, 2013-01-26 * Added more mass points at low masses to batch scripts (and corresponging plotting scripts and table scripts) * Added error catching routines for Earth files (WimpEvent reports an error if event file ends abruptly, as for the Sun files). -------- Revision 124, 2013-01-21 --------- * Fixed in bug in WimpEvent earth files. When a nadir angle was more than 90 degrees it was incorrectly output as negative. This is now fixed. The bug only affected low-energy neutrinos and only earth files. -------- Revision 115, 2013-01-09 --------- * Added option in WimpAnn to split event files. WimpEvent is updated to read these new files (which have longer headers). For the normal user, this option is not recommended, but for power users needing to distribute large runs on e.g. a cluster it could be useful. Please make sure to read the documentation on how the event weights are affected if you use this option. -------- Revision 111, 2013-01-07 ---------- * Version 3.02 created from trunk at revision 111. The main new ingredient is a bug fix relating to neutrinos with lower energy than 1 GeV (which were not included in the summary files as they should). * It also uses the new nusigma where a bug relating to the kinematical region caused unphysical configurations at times. -------- Revision 101, 2012-06-19 ---------- * Corrected header for Sun's azimuth and elevation (the numbers were correct, but in the header the azimuth and elevation were reversed). -------- Revision 97, 2012-06-04 ---------- * Version 3.01 created from trunk at revision 97. -------- Revision 93, 2012-05-28 ---------- * Version 3.00 created from trunk at revision 93. -------- Revision 85, 2012-05-28 ---------- * Added option to use not only SLA_EVP for positional astronomy, but also the more accurate SLA_EPV. SLA_EPV is much slower though and we do not need that accurary for WimpSim. Hence, the default is SLA_EPV. This is set with the parameter wesunacc. * Added azimuth, elevation, right ascension and declination to the particle tracks in the event files. Event format is slightly changed. * Added different ways to calculate the equinoxes and solstices (either via requiring that the declination is zero or max/min) or that the ecliptic longitude of 0, 90, 180 or 270 degrees. The latter is the astronomical definition. The two methods differ by up to a few minutes. * Note that WimpEvent requires WimpAnn event files of version 3.0 to work properly (it will stop if you give it an older file). -------- Revision 78, 2012-05-19 ---------- * A major revision adding support for true modified Julian dates (MJD) as time stamps for the events. This is faciliated through the use of SLALIB for the spherical astronomy calculations. In practice this means that one can now simulate events not over a simplified year as before, but instead over real years taking the Earth's motion and rotation properly into account. For averaged spectra, like those compiled for DarkSUSY, this will hardly matter, but it makes it easier for direct linking to Neutrino Telescope Monte Carlos, like for IceCube. For each instant in time during the simulation, the Sun's right ascension and declination are calculated. From these, the azimuth and elevation are also calculated. These are also stored together with the MJD for each event. The old timing treatment is still available, but for the new MJD treatment one can either enter the MJD or calendar date directly, or ask for specific times like solstices and equinoxes. These are then calculated on the fly using SLALIB. SLALIB includes the different versions of algorithms for e.g. the Sun-Earth positions and WimpEvent uses the most accurate one, SLA_EVP which agrees fairly well with other published calculations. The new interface routines (to calculate equnioxes and solstices etc) are included in WimpEvent with the prefix 'as'. SLALIB is also included under the Gnu Public License. -------- Revision 64, 2011-04-28 ---------- * A bug in nupropea caused WimpEvent simulations for the Earth to take a very long time in some cases. The results were most likely correct, but the CPU time was off (due to an unitialized variable). This mainly affected low-mass runs. -------- Revision 61, 2011-04-21 ---------- * A bug in CC only simulations have been fixed. It did not affect the event files and not CC + NC runs. The bug caused WimpEvent to pick the previous event for the summary tables when a tau (anti)neutrino between 1 and 5 GeV of energy was found. It only occurred in CC only simulations. -------- Revision 30, 2008-10-16 (tagged as wimpsim-2.11) * Changed summary file writing so that the KK channel (100) is written correctly (the file names contained '**' instead of '100'). * Changed makefiles so that we now use the latest full version of DarkSUSY instead of a light version (this does not lead to any physics changes, but is only from a practical point of view). -------- Revision 20, 2008-04-03 * Fixed incorrect Jacobian for integrated yield summary tables -------- Version 2.09, 07-09-20 * Fixed a minor bug that caused neutral current events below 1 GeV (5 GeV for nu_tau) not to be written to the summary tables. * Added more description in file headers for summary files -------- Version 2.08, 07-09-04 * Now works with nusigma-1.15 which includes new CC cross sections for (anti)tau neutrinos which include the tau mass. These are slightly lower than before and will mainly affect the weights and relative number of the events with tau neutrinos in the output event files. The weight changes for tau neutrino events as the total macroscopic cross section is now slightly reduced and the relative number of charged and neutral current events change as it is only the charged current cross sections that are reduced. * A bug for neutral currents is now fixed. Too many neutral currents on protons (as compared to neutrons) were performed. The weights were still OK, it was only the actual simulation of the events that were wrong. This bug only affected you if you looked at neutral current events. All the charged current events were OK (except for the improvement for tau neutrinos mentioned above). -------- Version 2.07, 07-05-22 * Added option to specify the seed to the Pythia random number generator. An option is to specify the seed from the date and time. If the current routines don't work for you, try replacing src/rndate.f with src/rndate-alt.f in the wimpevent directory. -------- Version 2.06, 06-12-01 * Fixed a bug with the time of events. If the starting time was not 0, incorrect times and Sun angles were obtained. This is now fixed. -------- Version 2.05, 06-09-27 * Fixed a bug which caused the path length through the Earth to be wrong (for neutrinos that came from the Sun). It was not very important, as oscillations are not that important through the Earth, but now it should be correct anyway. * Added the nadir angle of the Sun to each event (together with the time stamp). This simplifies the use of the event files when used for other neutrino telescopes than Amanda/IceCube -------- Version 2.04, 06-09-21 * The event files from the Earth are updated to include both nadir angles and azimuthal angles of the neutrinos, leptons and hadronic showers in a correct way. -------- Version 2.03, 06-05-10 * Event files for Earth updated to also include azimuthal angle. This is needed to unambigously distribute the events around the detector. Currently, this angle is set to zero, but the format is as it is going to be. In the next version, this angle will be calculated correctly. -------- Version 2.02, 06-04-30 * Bug in tabupdate which caused monochromatic neutrinos to be written in the wrong table (ifind returned -5 for enu/mx=1). -------- Version 2.01, 06-04-28 * Optimized oscillation routines through the Earth. They now treat the oscillations in fewer steps with the step boundaries at the mantle/core transition. * Summary files for flux type 3 (charged leptons at the detector) have been changed so that they only include mu- and mu+ (nui=3 and 4). * The summary files now contain the appropriate factors of dz and dtheta, so that they are really tabulations of dPhi/dz and d^2Phi/dz dtheta. Before these factors were added in DarkSUSY instead. -------- Version 2.00, 06-04-27 * Almost complete rewrite with many things changed: - WimpAnn is split into two programs, WimpAnn that takes care of the neutrinos to 1 AU from the Sun (or at the center of the Earth) and WimpEvent that takes these (event-based) neutrinos and propagates them to the detector - WimpEvent is now event based (no summary tables used to generate neutrinos). - Neutrino-nucleon interactions are now treated with nusigma instead of Pythia. - All neutrino flavours (and neutrinos and anti-neutrinos) are stored separately. - Neutrino oscillations are treated with a very general three-neutrino oscillation code. - The hadronic shower in the neutrino-nucelon interaction is now included. ====================================================================== -------- Version 1.45, 00-10-05 -------- * Updated to work with Pythia 6.154 -------- Version 1.44, 98-09-15 -------- * Option to create output files in a format directly readable by the Amanda Monte Carlo Schiffer about to be included. Name changed from numu to WimpAnn. -------- Version 1.43, 98-05-14 -------- * Bug in dflux fixed. It gave a factor of 100 too low fluxes for muons at creation. *** 1998-02-09 Version 1.42. Bug in savetab fixed. *** 1998-02-05 Version 1.41. Tested and bugfixed. *** 1998-02-03 Version 1.4. Converstion to Pythia 6.115. Also added *** option of no ntuples and writing to ascii tables *** directly. Converted to real*8 everywhere except the *** var.. arrays used for communication with HBOOK. *** 1997-10-30 Version 1.35. Added option to save small, medium or large *** ntuples. The medium ones keep all relevant experimental *** information. *** 1997-04-10 Version 1.34. Updated alpha and beta for muon stopping. *** 1995-11-30 Version 1.33. Small bug in musim which sometimes gave *** a nonreal angle theta fixed. *** 1995-11-29 Version 1.32. Double precision in muang.f to avoid *** numerical problems with heavy neutralinos *** 1995-11-17 Version 1.31. Particle decay table reset efter Pythia *** has been called since Pythia makes changes in it *** (important for Higgs channels e.g.) *** 1995-06-14 Version 1.3. Muon propagation in rock supported. *** m_mu = 0 no longer assumed. Multiple muons in CC *** interactions are now all stored. *** 1995-06-06 Version 1.2. Switch to Pythia 5.7 / Jetset 7.4 *** Switch to GRV LO structure functions included *** in Pythia. PDFLIB no longer needed. This means that *** extrapolation down to low Q^2 now is included. *** Bug in HADINT fixed. *** 1994-12-03 Version 1.1. Multiple Coulomb scattering included. *** 1993-12 Version 1.0 completed.