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

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