| 1 | ** Add changes at the beginning! **
|
|---|
| 2 |
|
|---|
| 3 | 16/04/2003 A. Moralejo
|
|---|
| 4 |
|
|---|
| 5 | MagicSoft/Simulation/Detector/Data/ :
|
|---|
| 6 | Added eco1000.def and axisdev1000.dat for eco-1000 studies
|
|---|
| 7 |
|
|---|
| 8 | 4/04/2003 A. Moralejo
|
|---|
| 9 |
|
|---|
| 10 | init.h, geometry.c, ph2cph.c:
|
|---|
| 11 | introduced variable ct_max_radius, the maximum distance in x or y from the
|
|---|
| 12 | center to the edge of the telescope, measured on the parabolic dish. This
|
|---|
| 13 | is used in a fast check of whether the photon hits the disk in ph2cph.c
|
|---|
| 14 |
|
|---|
| 15 | ph2cph.c:
|
|---|
| 16 | in function Lin2Curv(float x), changed the calculation of the curvilinear
|
|---|
| 17 | coordinates (which before was only valid for ~17 m focus paraboloids).
|
|---|
| 18 |
|
|---|
| 19 | 22/01/2003 A. Moralejo
|
|---|
| 20 |
|
|---|
| 21 | Fixed small bug in reflector.c. The end of run flag was not written to the
|
|---|
| 22 | output when the max_events option was used in the input card.
|
|---|
| 23 |
|
|---|
| 24 | 21/01/2003 A. Moralejo
|
|---|
| 25 |
|
|---|
| 26 | Removed from the output a null byte written right after the ascii label
|
|---|
| 27 | containing the program version number which is at the beginning of the rfl
|
|---|
| 28 | file.
|
|---|
| 29 |
|
|---|
| 30 | 19/12/2002 - 17/01/2003 A. Moralejo
|
|---|
| 31 |
|
|---|
| 32 | Lots of changes. Moved simulation of Mie scattering and Ozone absorption from
|
|---|
| 33 | attenu.f to atm.c, and changed the way it is done. Before it was not correct
|
|---|
| 34 | for large zenith angle (a variation like for Rayleigh scattering was assumed
|
|---|
| 35 | for both Mie scattering and Ozone absorption).
|
|---|
| 36 |
|
|---|
| 37 |
|
|---|
| 38 | 13/12/2002 A. Moralejo
|
|---|
| 39 |
|
|---|
| 40 | attenu.f: Found mistake in Mie absorption calculation (height from sea
|
|---|
| 41 | level was taken as height above observation level). Changed optical depths
|
|---|
| 42 | table for Mie attenuation. Now they are no longer referred to 2 km height,
|
|---|
| 43 | but to sea level. Detector level is taken into account later in calculation.
|
|---|
| 44 |
|
|---|
| 45 | 10/12/2002 A. Moralejo
|
|---|
| 46 |
|
|---|
| 47 | attenu.f: Added comments, removed old/unnecessary code, corrected small
|
|---|
| 48 | mistake in Elterman's table for ozone optical depth.
|
|---|
| 49 |
|
|---|
| 50 | atm.c,h : removed atmospheric model "ATM_ISOTHERMAL" which actually was
|
|---|
| 51 | not even implemented in the code, although it could be selected in the
|
|---|
| 52 | input card.
|
|---|
| 53 |
|
|---|
| 54 | 15/11/2002 A. Moralejo
|
|---|
| 55 |
|
|---|
| 56 | ph2cph.c: Fixed BUG in timing calculation!!! The time to be subtracted
|
|---|
| 57 | (mirror till ground) had the wrong sign!!!
|
|---|
| 58 |
|
|---|
| 59 | reflector.c: in the wavelength range check, changed the range limits to
|
|---|
| 60 | the fixed values 290 and 600 nm. The StarFieldAdder cer files could not be
|
|---|
| 61 | processed otherwise, because the w range is not written in their headers.
|
|---|
| 62 |
|
|---|
| 63 | ph2cph.c: changed sign of phi angle to make it positive. This was necessary
|
|---|
| 64 | since now the vector indicating the reflected photon trajectory has also
|
|---|
| 65 | changed sign! (see note below, 11/11).
|
|---|
| 66 |
|
|---|
| 67 | 14/11/2002 A. Moralejo
|
|---|
| 68 | reflector.c, parms.c: Added wobble mode option. Added wobble flag to the
|
|---|
| 69 | output and also an atmospheric_model flag.
|
|---|
| 70 |
|
|---|
| 71 | 13/11/2002 A. Moralejo
|
|---|
| 72 | reflector.c, header.h: Introduced 3 counters to keep track of what proportion
|
|---|
| 73 | of the C-photons in each event have been produced by electrons, muons or
|
|---|
| 74 | other particles.
|
|---|
| 75 |
|
|---|
| 76 | 11/11/2002 A. Moralejo
|
|---|
| 77 | ph2cph.c: BUG found! and fixed: the direction cosines in the cer files are
|
|---|
| 78 | those of the downgoing directions of photons, hence third component is
|
|---|
| 79 | negative. Before it was assumed positive. This, together with a wrong
|
|---|
| 80 | reflection in the individual mirrors (also changed) resulted in a subtle
|
|---|
| 81 | error which produced a quite de-focused reflector!
|
|---|
| 82 |
|
|---|
| 83 | reflector.c: for the same reason as above, the sign of theta in the calls
|
|---|
| 84 | to makeOmega had to be changed (+pi added to the phi argument).
|
|---|
| 85 |
|
|---|
| 86 | geometry.c, ph2cph.c, reflector.c, init.h: eliminated array ct_Focal[]
|
|---|
| 87 | which contained the focal lengths of individual mirrors. It was redundant
|
|---|
| 88 | since the same information is available in ct_data[].f
|
|---|
| 89 |
|
|---|
| 90 | Oct/2002 A. Moralejo
|
|---|
| 91 | added attach.c to attach the files magic.def, axisdev.dat and
|
|---|
| 92 | reflectivity.dat to the end of the reflector output.
|
|---|
| 93 |
|
|---|
| 94 | Sept/2002 A. Moralejo
|
|---|
| 95 | header.h, reflector.c: added run header (=Corsika's), changed event header.
|
|---|
| 96 |
|
|---|
| 97 | 3/7/2002 A. Moralejo
|
|---|
| 98 | reflector.c: Introduced NaN check in the photon loop.
|
|---|
| 99 | If NAns are found in a photon data block, it is not processed.
|
|---|
| 100 |
|
|---|
| 101 | ph2cph.c: Introduced "check of positiveness" before taking the sqrt
|
|---|
| 102 | in the calculation of the photon trajectory intersection with the
|
|---|
| 103 | global paraboloid.
|
|---|
| 104 |
|
|---|
| 105 |
|
|---|
| 106 |
|
|---|
