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