source: trunk/MagicSoft/Simulation/Detector/ReflectorII/Changelog@ 5447

Last change on this file since 5447 was 5447, checked in by moralejo, 20 years ago
*** empty log message ***
  • Property svn:executable set to *
File size: 7.5 KB
Line 
1** Add changes at the beginning! **
2
321/11/2004 A. Moralejo
4
5Added subdirectory writemagicdef under ReflectorII/. It contains the
6program to generate the positions and orientations of the mirrors.
7Directory contains Makefile and writemagicdef.c
8
9Moved back the default magic.def to be the one with no chessboarding,
10but with the final number of mirrors (964). It is more versatile to
11simulate the chessboarding at the camera level. For possible tests, a
12definitions file magic_chessboard.def has been created, but it is NOT
13recommended to use it.
14
15Removed file axisdev_no_chessboard.def, no longer necessary because now
16both the magic,def and the magic_chessboard.def files cobtain 964
17mirrors.
18
19geometry.[c,h]: added reading of token chessboarded_panels from magic
20definitions file. Nothing special is done with the token. It will simply
21be written to the ascii tail of the reflector simulation output, which
22will in turn be read by the camera simulation.
23
2403/11/2004 A. Moralejo
25
26reflector.c, header.c: Fixed problem: the correct telescope theta and
27phi were not written to the output of reflector, because the variables
28were filled before they were changed in the case of wobble mode. This has
29been corrected.
30
3127/09/2004 A. Moralejo
32
33A bit of cleaning in Makefile (remove useless CERNLIB)
34
35
3628/08/2004 A. Moralejo
37
38Changes in configuration files. Renamed Data/magic.def and Data/axisdev.dat
39to magic_no_chessboard.def and axisdev_no_chessboard.def. Now default
40axisdev.dat and magic.def contain 964 mirrors as in the present MAGIC
41configuration, and the panels are in a chessboard pattern
42
4330/04/2003 A. Moralejo
44
45ph2cph.c: Found out that the fix from yesterday, although it "patched" the
46loss of photons, did not remove the real cause for it. Actually in the
47determination of the intersection of the photon with the paraboloid there
48was an overflow which resulted in nonsense coordinates beyond 10 m from
49the mirror dish center. It has been fixed fixed by changing the following
50cut:
51
52 if ( fabs(a) < 1.e-3 ) {
53 xcut[2] = -c / b
54 }
55
56which was formerly < 1.e-6.
57
58I also noticed that the problem from yesterday was not so serious: it only
59removed a significant amount of the photons nearly paralel to the telescope
60axis, but globally, over a cer file there is no big difference.
61
6229/04/2003 A. Moralejo
63
64ph2cph.c:
65Changed determination of which mirror element is hit. Before there was a
66problem when going to very large distances from the mirror center, sometimes
67a wrong mirror was chosen, which was actually further away from the impact
68point than the mirror size, and therefore the photon was discarded later as
69"not hitting any mirror" (about 20% loss for eco1000) Now we measure the
70distance on the paraboloid. The new code is:
71
72 /* calculate the mirror to be used */
73
74 distmirr = 1000000.0f;
75
76 for (i=0; i<ct_NMirrors && distmirr>=ct_RMirror; ++i) {
77 distmirr2 = (float) sqrt(SQR(ct_data[i].sx - sx) +
78 SQR(ct_data[i].sy - sy));
79
80 if (distmirr2 < distmirr) {
81 i_mirror = i;
82 distmirr = distmirr2;
83 }
84 }
85
86Now the "distance" we calculate on the paraboloid is not a real distance,
87because sx,sy are measured along the cut between the paraboloid and planes
88containing the paraboloid axis and the x or y axis... but this should not
89matter much with such "flat" paraboloids. What cannot be done, as was before,
90is check the distance in the x,y,z space and then check whether the photon
91is inside the mirror in the sx,sy space, because then we lose photons for
92which a nearest mirror is chosen which in sx,sy is not close enough to
93contain the photon.
94
95geometry.c:
96Small change in calculation of ct_nax_radius (introduced fabs, but nothing
97changes)
98
9916/04/2003 A. Moralejo
100
101MagicSoft/Simulation/Detector/Data/ :
102Added eco1000.def and axisdev1000.dat for eco-1000 studies
103
1044/04/2003 A. Moralejo
105
106init.h, geometry.c, ph2cph.c:
107introduced variable ct_max_radius, the maximum distance in x or y from the
108center to the edge of the telescope, measured on the parabolic dish. This
109is used in a fast check of whether the photon hits the disk in ph2cph.c
110
111ph2cph.c:
112in function Lin2Curv(float x), changed the calculation of the curvilinear
113coordinates (which before was only valid for ~17 m focus paraboloids).
114
11522/01/2003 A. Moralejo
116
117Fixed small bug in reflector.c. The end of run flag was not written to the
118output when the max_events option was used in the input card.
119
12021/01/2003 A. Moralejo
121
122Removed from the output a null byte written right after the ascii label
123containing the program version number which is at the beginning of the rfl
124file.
125
12619/12/2002 - 17/01/2003 A. Moralejo
127
128Lots of changes. Moved simulation of Mie scattering and Ozone absorption from
129attenu.f to atm.c, and changed the way it is done. Before it was not correct
130for large zenith angle (a variation like for Rayleigh scattering was assumed
131for both Mie scattering and Ozone absorption).
132
133
13413/12/2002 A. Moralejo
135
136attenu.f: Found mistake in Mie absorption calculation (height from sea
137level was taken as height above observation level). Changed optical depths
138table for Mie attenuation. Now they are no longer referred to 2 km height,
139but to sea level. Detector level is taken into account later in calculation.
140
14110/12/2002 A. Moralejo
142
143attenu.f: Added comments, removed old/unnecessary code, corrected small
144mistake in Elterman's table for ozone optical depth.
145
146atm.c,h : removed atmospheric model "ATM_ISOTHERMAL" which actually was
147not even implemented in the code, although it could be selected in the
148input card.
149
15015/11/2002 A. Moralejo
151
152ph2cph.c: Fixed BUG in timing calculation!!! The time to be subtracted
153(mirror till ground) had the wrong sign!!!
154
155reflector.c: in the wavelength range check, changed the range limits to
156the fixed values 290 and 600 nm. The StarFieldAdder cer files could not be
157processed otherwise, because the w range is not written in their headers.
158
159ph2cph.c: changed sign of phi angle to make it positive. This was necessary
160since now the vector indicating the reflected photon trajectory has also
161changed sign! (see note below, 11/11).
162
16314/11/2002 A. Moralejo
164reflector.c, parms.c: Added wobble mode option. Added wobble flag to the
165output and also an atmospheric_model flag.
166
16713/11/2002 A. Moralejo
168reflector.c, header.h: Introduced 3 counters to keep track of what proportion
169of the C-photons in each event have been produced by electrons, muons or
170other particles.
171
17211/11/2002 A. Moralejo
173ph2cph.c: BUG found! and fixed: the direction cosines in the cer files are
174those of the downgoing directions of photons, hence third component is
175negative. Before it was assumed positive. This, together with a wrong
176reflection in the individual mirrors (also changed) resulted in a subtle
177error which produced a quite de-focused reflector!
178
179reflector.c: for the same reason as above, the sign of theta in the calls
180to makeOmega had to be changed (+pi added to the phi argument).
181
182geometry.c, ph2cph.c, reflector.c, init.h: eliminated array ct_Focal[]
183which contained the focal lengths of individual mirrors. It was redundant
184since the same information is available in ct_data[].f
185
186Oct/2002 A. Moralejo
187added attach.c to attach the files magic.def, axisdev.dat and
188reflectivity.dat to the end of the reflector output.
189
190Sept/2002 A. Moralejo
191header.h, reflector.c: added run header (=Corsika's), changed event header.
192
1933/7/2002 A. Moralejo
194reflector.c: Introduced NaN check in the photon loop.
195If NAns are found in a photon data block, it is not processed.
196
197ph2cph.c: Introduced "check of positiveness" before taking the sqrt
198in the calculation of the photon trajectory intersection with the
199global paraboloid.
200
201
202
Note: See TracBrowser for help on using the repository browser.