source: trunk/MagicSoft/Mars/NEWS@ 8988

Last change on this file since 8988 was 8988, checked in by tbretz, 16 years ago
*** empty log message ***
File size: 118.6 KB
Line 
1 -*-*- END -*-*-
2
3== <cvs> ==
4
5 ;general
6
7 * should work now with root 5.20/00
8
9
10
11
12== Version 2.1 (2007/09/03) ==
13
14 ;Database
15
16 * introduced average temperature and wind speed from star-files
17
18 * introduced average values from Pyrometer (mainly Cloudiness)
19
20 * introduced the ratio of accepted calibration events in the
21 calibration run
22
23 * introduced the average rates of the different trigger pattern,
24 e.g. the rate of pedestal and calibration events, from the signal-file
25
26 ;general
27
28 * fixed GroupBy options in plot*.C macros
29
30 * only display KVA points (not the Tuorla points) in plotoptical.C
31
32 * all programs now allow to add resources to the ones from the resource
33 file or overwrite them from the command line. This can be done like this
34 ganymed --config=ganymed.rc --rc=Cut1.Param1:0.25 --rc=Cut1.Param3:5.5
35 if something seems to be wrong you can debug what the program is doing
36 using the option --debug-env=3
37
38 * increased the diversity of debug levels. That also means that with -v3
39 you get now less output, i.e. only the output which is important for
40 an analysis. To get more output which helps for debugging try 4, 5 and 6.
41
42 * fixed a problem in the task synchronizing the reading of the reports
43 and the data stream. The effect was that in some circumstances the last
44 events of some trees could have been skipped. Since this is only a
45 very tiny part of the data no severe effect to your analysis results
46 is expected (an exception might be very short, i.e. singel runs,
47 datasets).
48
49 * Dataset file now allow to include the sequence file directly within
50 the dataset file and also allow to exclude single runs from the dataset.
51 (for more details see Class reference of MDataSet)
52
53 * Sequence files now allow for excluding runs from the analysis
54 (for more details see Class reference of MSequence)
55
56 * added code to process the sum-trigger flag
57
58 * the histogram classes MH3 and MHn now support profile histograms
59
60 * the histogram classes MH3 and MHn now support labels at the axis
61
62 * the Random Forest Train-classes (MJTRain*) now store the used datasets
63 in the output file
64
65 * the MFMagicCuts have a new option to allow a linear area cut (for
66 special studies)
67
68 * The liniking of the shared object is now done into a file defined
69 by 'mktemp'. This should be a local file system which accelerates
70 liniking a lot.
71
72 * The code has been prepared for compilation with root 5.18/00d
73
74 * The MHEnergyEst histogram now shows the distribution of
75 (Eest-Emc)/Est and the distributions (Eest-Emc)/Eest vs. Eest
76 and (Eest-Emc)/Emc vs Emc.
77
78 * fixed mars (event display) for merpped raw-data and MC data
79
80 * adapted the synatx of sequence files to allow for the new run/file
81 scheme. For more details see class description of MSequence
82
83 ;showplot
84
85 * Added support for tiff files
86
87 * Added support for csv files (tables containing tab number,
88 name and title)
89
90 * Tabs now can have a title which is displayed in the postscript file
91
92 * When multiple file are written at once (e.g. pngs from a display)
93 it is now possible to use the tab name and or tab number in the
94 file name. For details see showplos's help.
95
96 ;merpp
97
98 * Implemented cc file version 2008-02-20/0
99
100 * Implemented cc file version 2008-05-19/0
101
102 * With new options merpping can be restricted (in addition to the
103 run-number) also by telescope and file number.
104 For details see "merpp --help"
105
106 * Implemented raw data format 11 (the most important change, except
107 from plenty of changes in the format definition, is the addition
108 of the file-number in the header) If nobody makes a mistake this
109 format should have the advantage that it is by definition forward
110 compatible, so old programs will be able to read newer formats.
111
112 * Added Pyrometer information to default output (stored in a new Tree
113 Pyrometer, needs a complete new merpp, i.e. new calibration)
114
115 * Merpping (also update) can now be restricted to one single report.
116 For example:
117 If you want to update an old file with Pyrometer information call
118 merpp --only=Pyrometer --update
119
120 ;callisto
121
122 * in the pedestal processing now at least 50 processed events are
123 required for succeeding.
124
125 * the signal and calib files now contain also the sequence information
126
127 * the two MPedestalCam in the calib file now have proper names
128 to distinguish them from each other (MPedestalCam, MPedestalExtracted)
129
130 * Both pedestal cams now contain the number of events processed
131
132 * In Mars V2.0 the extraction of the fundamental pedestal was broken,
133 i.e. for the extraction of the calibration pulses basically the average
134 over all events of the first slice instead of all slices was used.
135 Since pedestal and calibration constants are continously recalculated
136 this only effected the very first events of every sequence.
137
138 * For the extraction of the position of the maximum sample now the
139 extraction range is set automatically to the extraction window,
140 before the whole accessible range was used.
141
142 * If the number of samples in the events and the number of samples
143 in the run-header disagrees now an error is raised (this can
144 happen in not well simulated Monte Carlo files)
145
146 * In the case of MUX-data pedestal events with the lvl1 trigger flag
147 could survive into the pedestal calculation -- fixed.
148
149 * Tab "TrigPat" added showing the distribution of the trigger pattern
150 as found in the run-files (normalized with the run length)
151
152 * The selection of events from the trigger pattern has been improved.
153 If the trigger pattern is corrupted events won't pass anymore.
154
155 * Added a new command line option --dev-null to suppress output of
156 Y-files (this is useful for test cases)
157
158 * added a new tab "CalPos" showing the arrival time of extracted
159 interleaved calibration events
160
161 ;star
162
163 * the star file now contains also the sequence information
164
165 * Added a new command line option --dev-null to suppress output of
166 I-files (this is useful for test cases)
167
168 * Tab "Rate" added showing the distribution of the trigger pattern
169 (should only be Trig (Lvl1/Lvl2) and Sum (Sum only)) as found in
170 the Y-files (normalized with the run length)
171
172 * The effective on-time calculation doesn't use events with only
173 sum-trigger anymore
174
175 * The data in the MHWeather tab has been reorganized. The never
176 working solar radiation has been removed and the data from
177 the pyrometer (cloudiness, air and sky temperature) is
178 displayed in addition.
179
180 ;ganymed/sponde
181
182 * Now data points are even shown correctly if they have error bars larger
183 than one sigma. Points with null or negative significances are omitted.
184
185 * DrawNicePlot now allows to set a projection range
186
187 * Updates pointing models. One is applied after the implementation
188 of the new LUTS in August, the other one is applied after the
189 next major change of the pointing model in October. Note that
190 these models are preliminary and give wrong results below Zd=65deg
191
192 ;sponde
193
194 * fixed resource files ('''Spectrum''' instead of '''MJSpectrum''')
195
196 * a new option "--force-runtime" is available. In case of very short
197 datasets (for light-curves) which are in the order of a few
198 minutes the calculated effective on-time is not very acurate because
199 it is quantizised in th eorder of a minute. Therefore you can switch
200 to using the real run-time instead (remakrt: this ignores any
201 dead-time!)
202
203 * a more analytical way is used to "fill" Monte Carlos into regions
204 of impact parameters which have not been simulated because the
205 trigger efficiency is assumed to be 0.
206
207 * estimated sensitivity curves for hi and lo zenith angles are shown
208 scaled to the correct observation time and collection area
209
210 * the Monte Carlo events after cuts are now written to the output file
211 if an output file given
212
213
214
215 *** Version 2.0 (2007/09/03)
216
217 - database: the position used for the pulse pos check was missing
218 in the callisto output and could thus not be filled into the DB.
219 Also the old values seemed not exactly the PulsePos used for
220 teh check.
221
222 - general: Resource file now allow an Include-Resource, i.e. you can
223 read a resource file with default settings, include it in your
224 resource file and overwrite the settings in your file. More than
225 one include file is allowed. Inclusions can be iterative.
226 Include: mydefaults.rc yourdefaults.rc
227 The resources in the first file have higher priority than the
228 second file.
229
230 - general: Now the output files (calib*.root, etc) also contain the
231 resource file (to check it open it in the TBrowser and choose Print()
232 from the context menu)
233
234 - general: Dataset files now allow to overwrite the default path
235 to search for sequence and data files. Note, that the program
236 option and individual resources will overwrite it. For example:
237 SequencePath: /magic/montecarlo/sequences
238 DataPath: /magic/montecarlo/star
239
240 - general: Fixed a problem with variables in MDataPhrase. It seems
241 that only optimdisp (didn't work) was effected.
242
243 - general: Some improvements to the optim and train processing. For
244 example the ststud window now got a proper title.
245
246 - general: The "Status Display: " was removed from the window title
247 within an Eventloop.
248
249 - general: added a new flexible histogram-class (MHn) which can be used
250 for example in testing after optimizations or training to produce
251 user defined histograms.
252
253 - general: The Random Forest now allows attaching an evaluation
254 function which is applied to its output. This allows to train
255 on whatever you like (eg. train on log10(energy), but return
256 energy)
257
258 - general: The Random Forest training of energy (trainenergy.c) and
259 disp (traindisp.C) has been improved by more graphical output to
260 judge the result.
261
262 - general: The training of energy (trainenergy.C) has been improved
263 a lot by training on a different quantity than energy. It gives
264 a result which is highly independant of the source spectrum
265 and shows lower bias than all other options.
266
267 - general: Dataset files are allowed to contain collections of
268 datasets. For more details see the class reference of MDataSet.
269
270 - general: The default WobbleMode in a datset file now (if not
271 overwritten) is "auto", i.e. wobble mode is set if the dataset
272 doesn't contain off-sequences.
273
274 - merpp: Merpp can read file format version 9 now. Please realize
275 that for file format version 9 the extraction range is different
276 and must be set to 0/49 instead of 15/64 in callisto_mux.rc until
277 the first and last 15 slices have been removed from all MUX data.
278
279 - mars: Similar to writing movies you can now select events which
280 should be displayed in your display. Therefore add the following
281 lines to your mars.rc:
282 MFEvtNumber.FileName: ganymed00223552.root
283 MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5
284 The input file can be a ganymed-file after (ganymed*.root) or
285 before (ganymed*-summary.root) file. For example to select all events
286 from your on sample after cuts (excluding the ThetaSq cut) use
287 MFEvtNumber.Selector: DataType.fVal>0.5
288 To show all events from your on-sample after quality cuts
289 use the summary file instead.
290
291 - mars: new ways to call mars are implemented. Instead
292 of calling
293 mars filename.root
294 you can now also use
295 mars sequence.txt
296 or
297 mars sequence.txt inputpath
298 or replace sequence.txt by the sequence number if the default
299 sequence file should be used.
300
301 - mars: The default size has been increased
302
303 - mars: There are two new heckboxes which allow to switch off
304 the calculated image parameters
305
306 - callisto: was broken for MCs... fixed.
307
308 - callisto: improved calculation of spline coefficients a lot. This
309 leads to a further improvement of the event rate calibrating MUX
310 data of about 15% (175evt/s instead of 150evt/s)
311
312 - callisto: finally the raw data we read is now real 16bit (with
313 10bit precision) for MUX data and not a stripped 8bit number
314 anymore. Due to scaling of the numbers (without precision loss
315 into an [0-256[ range, the final numbers for pedestal, pedestalrms
316 and extacted signal don't change. However, at all places where
317 raw data is directly accessed (eg. saturation and variation
318 limits in the signal- and pedestal-extraction) the numbers
319 refer still to the data's 16bit range.
320
321 - callisto: The maximum arrival time difference used in the bad
322 pixel treatment is now in units of nanoseconds, which leads to 3.0ns
323 for both, old FADC and MUX-FADC, data.
324
325 - callisto: Be aware that callisto is now calibrating the arrival time
326 in nanoseconds rather than time slices. Also all following tasks
327 have been changed to assume ns units now. If you star data
328 calibrated with an older version it might be necessary to adapt
329 some timimng parameters (Muon Analysis and MFSoftwareTrigger)
330 to obtain optimal results. The plot in the PulsePos tab
331 will now be in ns, too. This might also mean that the Muon
332 analysis in star might not work as perfect as expected as long
333 as old files read in.
334
335 - callisto: The calibration constants of earlier updates got lost
336 somehow. All constants have been updated.
337
338 - star: The PSF is now determined from the profile of the ArcWidth
339 instead of arcwidth/radius. The old way gave to much weight
340 to the bins with low statistics. The reference lines have been
341 updated.
342
343 - star: For speed reasons events suitable for the muon analysis
344 are now also preselected by fConcCOG<0.1 which seems to be a very
345 good preselector for muons.
346
347 - star: is now calculating a new timing parameter the Slope along
348 the major and minor axis of the shower. Therefore a numerical
349 solution of a line-fit is used. To use MHillasExt.fSlopeLong
350 and MHillasExt.fSlopeTrans don't forget to multiply it with
351 the sign stored in MHillasSrc.fCosDeltaAlpha as usual.
352
353 - star: The parameters fInnerSize, fInnerLeakage1 and fInnerLeakage2
354 have been removed from MNewImagePar. They have never been good for
355 anything.
356
357 - star: Implemented a new recursive algorithm which is supposed to be
358 faster. Also implemented new cleaning oiptions, which are:
359 + CleanLevel0:
360 The cleaning level above single core pixels can be kept
361 + KeepIsolatedPixels:
362 Define whether single core pixels should be kept
363 + TimeLevel1:
364 The coincidence window in nanoseconds to single used neighbors
365 + TimeLevel2:
366 The coincidence window in nanoseconds to two used neighbors
367 + PostCleanType:
368 MOde in which a time dependent post cleaning should be applied
369 0: No post cleaning
370 1: Require for each used pixels at least one used neighbor
371 within a time-window of CleanLevel1
372 2: Require for each used pixels at least two used neighbor
373 within a time-window of CleanLevel2
374 3: Require first two and second one used neighbor
375 within a time-window of CleanLevel1 and fCleanLevel2
376 respectively
377 + RecoverIsolatedPixels:
378 - recoveres isolated core pixels if they have a used neighbor
379
380 - mars/star: The new default cleaning is:
381 + CleanLevel1: 6.0
382 + CleanLevel2: 3.0
383 + CleanTime1: 1.75
384 + CleanTime2: 1.75
385 + PostCleanType: 3
386
387 - star/ganymed: The old spark-cuts have been replaced by new ones.
388 These new spark cuts have been cross checked with the sequences
389 84720 (before splitter), 101041 (after splitter) and 223539 (new
390 FADCs). In all cases they seperate pretty well, but of course not
391 perfect. Monte Carlos have suggested to tighten the cuts a little
392 bit further at low sizes. This has been done.
393
394 - ganymed: IMPORTANT - The spark cuts have been changed to match
395 all available data more or less well. This might mean that for
396 YOUR data there are still sparks visible.
397 Whenever you do an analysis you have to make sure that NO sparks
398 survive your cuts!
399
400 - ganymed: The old hadronness and size cuts have gotten new indices
401 (10, 11). The old indices 8 and 9 now mean an additional condition
402 for ghostbusting based on the formula (dist-c[9])*c[8]-slope.
403 If you want to switch off this cut just move the line far away enough,
404 e.g. by setting c[9] to -500. Please update your ganymed.rc files
405 accordingly!
406 The unit of c[9] is deg. It is the dist at which your cut will
407 cross the slope==0. The unit of the slope is ns/deg. So if
408 you readin files calibrated with an older Mars-version you
409 have to adapt the default value for c[9] to the FADCs:
410 - Multiply by two for the MUX FADCs (run number > 200.000)
411 - Multiply by 0.3 for the Siegen FADCs (run number < 200.000)
412 - Multiply by 42 for the meaning of life (run number = 200.000)
413
414 - ganymed: The old Cut0 is now called CutQ and is calculated before
415 the source position is caluclated. If you need a quality cut
416 including source dependant data (be very carefull with this!)
417 you have to use Cut0 instead.
418
419 - ganymed: now properly supports three off regions in wobble mode.
420 It will be the default from now on in ganymed_wobble.rc. To change
421 the number of off-regions use
422 NumOffSourcePositions: 1
423 in your ganymed.rc. Also by default no off-theta cut will be done
424 anymore so that you will get a background level higher than the
425 signal level in your ThetaSq plot for ThetaSq>0.2.
426 Due to this be carefull changing the scaling mode away from
427 "None". You will underestimate your signal if your scale interval
428 is not properly set.
429 (Rem: ThetaSqN is no longer needed)
430
431 - ganymed: is now storing the events of all (on- and off-)
432 source positions even in the case of more than one off-source
433 position. This should allow sponde to work properly and the
434 energy estimation is done correctly in all cases (without any
435 "nearest position" trick)
436
437 - ganymed: From 85340 (19.3.2006) on a real starguider calibration
438 with a real pointing model for the starguider is done, calibrated
439 with more than thousand tpoints taken from this date on.
440
441 - ganymed: was broken for Wobble-MCs... fixed.
442
443 - ganymed: The default cuts (MFMagicCuts) and the disp-parametrization
444 has changed. Note that your old ganymed.rc will be incompatible
445 with the new algorithms! Please check the ganymed.rc files and
446 the class description of MFMagicCuts for more information.
447 It is recommended at the moment to use identical Disp-Parametrizations
448 for on/off- and wobble-mode. If cuts are optimized only
449 the Area-Cut (parameters 2,3,4) and the theta-sq cut (parameter 1)
450 should be optimized. All other values should be kept as they are.
451 If you want to do an anylsis without timing a parameters, you can
452 switch off the timing cuts by seeting parameter 7 to a very small
453 value (eg. -99). To remove the slope parameter from the disp
454 paremtrization set parameter 8 to 0.
455
456 - ganymed: The distribution of observation time versus zenith
457 angle is now shown for off- AND on-data.
458
459 - ganymed: The contents of the source position plot are no longer
460 averaged, thus a lot of (fake) events between the two wobble
461 positions have disappeared.
462
463 - ganymed: There is a new tab "CutT" it contains the VsSize plots
464 with the CutT, whihc is defined in ganymed.rc, applied. It
465 is resonable to use your theta-cut for it, to see if there is
466 something strange (sparks!) in your signal-region.
467
468 - ganymed: To choose a dataset from a collection in a dataset file
469 use the new command line option "--dataset=12345"
470
471 - optim, sponde: should now properly support three off-regions.
472 Just produce your ganymed summary files with three off-regions.
473 optim and sponde will automatically get all off events from
474 all off-regions from this file. Please do not forget to
475 set the scaling in optimwobble.C properly:
476 fit.SetScaleUser(1./3);
477 it is not yet read in automatically (will follow soon).
478 In the case of fixed Scale you can also use
479 cuts.SetThetaCut(MFMagicCuts::kOn);
480 no off-cut needed.
481 Be carefull in case of three off-regions when switching on
482 scale mode. Make sure that the scale interval is set correctly.
483
484 - sponde: sponde.rc and sponde_onoff.rc are now idetical
485
486 - sponde: the resouce files have been rewritten with a lot of
487 comments and a well working parametrization for an energy estimator
488
489 - sponde: In the estimated energy versus monte carlo energy plot the
490 contents where exchanged.... fixed.
491
492 - sponde: the so called "accurate"-mode has been removed. It didn't
493 give any improvement in accuracy, only decreased execution speed.
494
495 - sponde: the so called "simple"-mode has been removed. It didn't
496 give any improvement in simple.
497
498 - sponde: the so called "refill"-mode has been removed. It was anyhow
499 not implemented.
500
501 - sponde: now checks whether the theta distribution of your on-data
502 and the theta-distribution of your Monte Carlo sample (after
503 weighting) fits. If it doesn't fit properly (eg. the Monte Carlo
504 sample is incomplete) execution is stopped. Execution can be forced
505 using the new option --force-theta. Use this option with care!
506
507 - sponde: Proper collection areas can now be constructed also from
508 Monte Carlo samples generated with different maximum impact
509 parameters. Note that in previous version you neither got
510 a warning or failure, nor was there any obvious sign that the
511 collection area was overestimated due to usage of files with
512 different maximum impact parameters.
513
514 - sponde: If MC files with different lower energy limits are used
515 the primary MC spectrum is artificially completed down to the
516 lowest energy used at all. WARNING: that this gives correct
517 collection areas ONLY if none of the events in this region would
518 survive your cuts at all.
519
520 - sponde: the output file now contains more information about
521 the spectrum (eg. the full 2D collection area histogram).
522 Note, that this information can only be written to the file
523 if it is stored automatically via command line argument.
524 If you only store the status display from within the display
525 the information is lost.
526
527 - sponde: added a new tab "Disp". It shows the residual of the disp
528 (Disp-Dist) versus several parameters and it is meant to judge the
529 quality of the disp estimator. In the ideal case the residual is zero
530 and doesn't depend on any variable. The tab is displayed after
531 CutQ and Cut0.
532
533 - sponde: added a new tab "Energy". It shows the residual of the
534 energy (lg(estimated energy)-lg(monte carlo energy)) versus
535 several parameters and it is meant to judge the quality of the
536 energy estimator. In the ideal case the residual is zero
537 and doesn't depend on any variable. In reality even a good estimator
538 can show residuals versus Monte Carlo energy. The tab is displayed
539 after all cuts.
540
541 - sponde: added a new tab "EventDist" showing the unweighted real
542 absolute number of events of your sample after cuts. This tells
543 you how many events with this energy you had in your MC files.
544 The same information you get from the error bars of the weighted
545 histograms, but this is less intuitive.
546
547 - sponde: The spectrum plots now show the crab- and 1553-spectrum
548 for comparison. It is not meant to show these curves in
549 publications, they are only for production.
550
551 - sponde: The OriginalMC tree with the events produced by corsika
552 is now processed only once
553
554 - sponde: Finally fixed most annoying bug, which has effected the
555 last bin of the collection area. It was the problem that binnings
556 of root histograms are from 1 to n (included) but there was a C-like
557 loop in sponde from 0 to n (excluded). Now the loop also includes
558 correctly under- and overflow-bins. The worse thing with this is,
559 that for the events in the highest bin not only the bin-content
560 was wrong, but also the weights applied to these events, which could
561 in pricipal effect also other parts of the distribution (eg. due
562 to energy estimation). Fortunately it effected only a few events
563 in most cases.
564
565
566
567 *** Version 1.2 (2007/05/14)
568
569 - database: The database now has two new values Unsuitable50 and
570 Unsuitable01. They express the number of pixels which are
571 unsuitable for more than 50%, respectively more than 1%, of
572 all calibrated event. It is a more accurate number than the
573 previously used numbers, because they only take the first
574 calibration (ignoring the interleaved calibrations) into account.
575
576 - database: The database now has two other new values UnsuitableMax and
577 DeadMax. They express the maximum number of pixels which were
578 unsuitable, respectively dead, during the sequence. Because
579 of high pedestal rms (cars passing) a few events with very high
580 numbers of unsuitable pixels can happen. Not to suffer from this
581 effect we don't take the highest 0.1% of the numbers into account.
582
583 - general: fixed a bug which caused callisto and star to stop working
584 properly because the callisto output was currupted
585
586 - general: a script called scripts/makemovie has been added which
587 simplified the process of producing several movies from a single
588 sequence with different setup.
589
590 - callisto: Some more code cleanup. Started to remove old obsolete
591 code from cvs.
592
593 - callisto: The extraction of the pedestal with the extractor was
594 not random. It seems that extracting at one position is still biased
595 (for example in the first or last slices we could still suffer from
596 switching noise) Now it is completely random.
597
598 - callisto: In MPedCalcPedRun (which is the pedestal extraction from
599 pedestal files) individual pixels could be skipped due to high
600 variations. This was already the case for the pedestal extraction
601 from the lo-gains since the beginning and properly handled
602 there, but was introduced in the extraction from the pedestal files
603 a while ago. To calculate the average value it was still divided
604 by the number of events procesed not by the number of summands
605 really summed. This let to a pedestal which was a tiny amount too
606 small (for a few piels in the order of 0.1%). This led to an slightly
607 positive offset of the randomly extracted pedestal for a few pixels.
608 I doubt that this has a big effect on the result, because the effect
609 on the individual numbers is quite small. Due to this it is suggested
610 that you rerun your calibration to make sure you are without any
611 avoidable bias.
612
613 - callisto: the ArrTimeRmsLimit is now is a check of the deviation from
614 the median rms of the absolute arrival time and expressed in significance
615 levels. This replaces the old meaning of an abolute upper limit. An
616 absolute upper limit doesn't make much sense, because in our data
617 all arrival times are shifted at the same time (in fact it might be
618 a shift of the artificial trigger w.r.t. to the light flash).
619
620 - callisto: Added a new tab "BadPixTm" which shows the time evolution
621 of the number of unsuitable pixels over the whole sequence. This
622 is usefull mainly to judge if an intermediate calibration had
623 problems.
624
625 - callisto: Added a new tab "DeadPixTm" which shows the time evolution
626 of the number of dead pixels over the whole sequence. Dead pixels
627 in this context are unmapped pixels, i.e. pixels which could not
628 be interpolated, and thus are ignored in the further analysis.
629
630 - callisto: It is now possible to use the position of the maximum
631 of the spline as arrival time instead of the leading edge even
632 when using the integral for the signal. Therefor set the extraction
633 type to kIntegralRel ("IntegralRelative") and the relative height
634 ("HeightTm") to -1:
635 MJPedestalC1.ExtractSignal.ExtractionType: IntegralRelative
636 MJPedestalC1.ExtractSignal.HeightTm: -1
637
638
639
640
641 *** Version 1.1 (2007/04/27)
642
643 - general: The new class MGeomCamDwarf contains the camera geometry
644 for our prefered Dwarf camera and can be used to create general
645 geometries of roundish or hexagonal cameras.
646
647 - mars: The default cleaning is now abolute 8.5/4.0.
648
649 - mars: The resource file mars.rc now contains examples for MUX data
650
651 - merpp: Implemented file format version 8 (MUX FADC data). The 16-bit
652 data containing only 10-bit information is still cooked down to
653 8-bit, by cutting away the lowest 8-bit.
654
655 - merpp: When reading raw data the data is now stored in a single
656 array. The lo-gain array is obsolete. The interface stays the same.
657
658 - callisto: To process the new MUX data use the resource file
659 callisto_mux.rc
660
661 - callisto: The -root and -raw options to read merpped raw data have
662 been removed
663
664 - callisto: If the -mc option is given the sequence-name can now be
665 replaced by the MC sequence number
666
667 - callisto: If the sequence file contains "MonteCarlo: Yes" the
668 -mc option can now be omitted
669
670 - callisto: Started further simplification by removing obsolete
671 dependancies on MCalibIntensity*Cams.
672
673 - callisto: removed unsused intensity calibration from MJCalibration
674 for further simplification
675
676 - callisto: for simplification removed references to
677 MBadPixelsIntensityCam, MCalibrationIntensityQECam,
678 MCalibrationIntensityRelTimeCam, MCalibrationIntensityTestCam,
679 MCalibrationIntensityChargeCam, and MCalibrationIntensityBlindCam
680 (this should at minimum change nothing at all, and at maximum
681 improve thing if these containers where somewhere not handled
682 properly)
683
684 - callisto: Added the option to write an event movie. Please use
685 this option with care (huge files) but exessively to understand
686 our events. For further documentation please refer to the
687 IACT wiki
688 (currently at http://www.astro.uni.wuerzburg.de/mediawiki-1.9.0)
689
690 - callisto: Improved binning for the pulse position check histogram
691
692 - callisto: The code calculating the pulse position checknow takes
693 the pedestal into account
694
695 - callisto: The arrival time is now stored with full precision this
696 helps to get rid of artefacts in histograms. For the signal this
697 is not necessray because it is displayed in log-scale anyhow
698
699 - callisto: Fixed a problem with internal exclusions. Pixels which
700 have once been marked as unsuitable didn't get the necessary
701 histograms filled anymore to be rechecked with each interleaved
702 calibration.
703
704 - callisto: The relative arrival time resolution (the rms of the
705 relative arrival time offset) is now checked against its
706 median and the median deviation, instead of the average
707 plus an absolute offset.
708
709 - callisto: The relative arrival time offset is now calculated
710 w.r.t. its median not w.r.t. to an artificial reference pixel.
711
712 - callisto: The relative time calibration now marks pixels with
713 outlaying offset as unsuitable
714
715 - callisto: The interleaved relative time calibration has been switched
716 on. The advantage is that pixels which are marked as outliers
717 now can recover. An improvement over time is not expected, because
718 the relative time calibration seems stable.
719
720 - ganymed: If the dataset file containes "MonteCarlo: Yes" the "-mc"
721 option can now be omitted.
722
723 - ganymed: If either the dataset file containes "MonteCarlo: Yes"
724 or the option "-mc" is specified as default the monte carlo datacenter
725 paths (/magic/montecarlo/...) are used if not overwritten in the
726 dataset file.
727
728 - ganymed: The output files now contain run- and event-number for all
729 events, accessible as RunNumber.fVal and EvtNumber.fVal
730
731 - ganymed: Added a plot showing the source position in the camera also
732 for wobble mode
733
734 - ganymed: In addition to the existing task which can be setup to
735 be executed after Cut1 a final task called "EstimateEnergy"
736 has been added, executed before writing the final output
737 (ganymed00000000.root). The task setup in all three cases
738 can be choosen individually. (For all three(!) cases you
739 can for example choose MRanForestCalc. For examples see all
740 sponde*.rc and all ganymed*rc)
741
742 - ganymed: If a container "MEstimatedEnergy" is found in the
743 parameter list (eg. because it was setup by MEstimateEnergy)
744 it will be written to the output file.
745
746
747
748 *** Version 1.0 (2007/02/16)
749
750 - database: Website to plot values from the database has been added.
751
752 - database: The rate of spark events as found by star is now kept
753 and the rate after cleaning is now with spark events subtracted
754
755 - general: The changes to the interpretation of phrases lead to
756 a problem with the access of MHMatrix objects (eg. "M[0]")... fixed.
757
758 - general: Added example files how you can create your own calendar
759 (eg. as christmas present) with Mars. Just run
760 root macros/tutorials/calendar.C
761
762 - general: Added a possibility to display grouped data in the plotdb.C
763 macro (and similar macros). This allows to plot hourly, nightly,
764 weekly, monthly and yearly averages.
765 For more details see the plotall function of the macro
766
767 - general: fixed the error display of the camera display. It showed
768 wrong results.
769
770 - general: fixed a memory leak when reading many MStatusDisplays
771
772 - general: added a macro (macros/tutorials/mirrordelay.C) plotting the delay
773 between a spherical and a parabolic mirror
774
775 - callisto: We have a new pulse position check which takes the properties
776 of the extractor automatically into account. So no need to set
777 ExtractWinLeft/Right anymore. Again a check for pulses too much to the
778 left is done. If possible the lo-gain extraction start and hi-gain
779 extraction end is changed accordingly. If no proper extraction range
780 can be set (hi-gain too much to the left) an error is raised.
781
782 - callisto: New swapped pixels have been found and corrected
783
784 - callisto: Redone the signal extraction this includes:
785 + Unified the algorithms to extract the pedestal for all pedestal
786 extractions done
787 + Unification of hi- and lo-gain into one array
788 + Unification of hi- and lo-gain extraction for Spline and
789 Digital Filter (means now exactly the same algorithm is used,
790 not a copy-and-pasted one)
791 + Store pedestal subtracted data in a common place
792 + Improved extraction algorithms. If the algorithms search a
793 maximum but don't find one they don't return the maximum anymore
794 + Improved the extraction such that a information for signal
795 and time which is consistent is returned
796 + replaced the numerical spline with a (not slower) analytical
797 solution
798 + Lowered the saturation limit from 250 to 245 to make sure that
799 we are far away in all pixels from problems with the clock-noise
800 + The start of the lo-gain extraction is calculated from the hi-gain
801 arrival time. In case of saturation an estimation is used.
802 + the hi-gain window was extended to slice 16. This allows extraction
803 of very late pulses in single pixels or of huge hadron showers.
804 + Now uses the lo-gain extraction range for pedestal extraction,
805 not the hi-gain range in the lo-gain
806
807 - callisto: New hi-/lo-gain calibration constants have been produced
808 for all periods from the statistic of the whole period. They are
809 now used automatically depending on the periods. The old constants
810 where deviation in avarega about 2%-3%. The constants are very
811 stable from period to period (<<1%) except for exchanged channels.
812 In a few cases the constants used so far were off by 1%-2% (maybe
813 temperature effects?)
814
815 - callisto: To gain obsolete space we do not write the
816 MPedestalFundamental anymore. Also calibration and bad pixel
817 treatment is not done for it anymore because it has never been used.
818
819 - callisto: Fixed a bug in the seuqence handling which caused the
820 light condition to be ignored
821
822 - star: Added two new tabs "Sparkless" and "Sparks" these tabs show
823 basically the same as "Cleaned" but the distribution devided into
824 the events fullfilling a spark cut or being rejected by it.
825 This cut can be changed from star.rc (see example inside) by
826 the "SparkCut" directive
827
828 - star: Star displays now the average individual pixel-rate (IPR) versus
829 time and the average discriminator threshold (DT) of all pixels.
830
831 - optim: fixed a problem with the optim-macros which was due to a bug
832 in the new phrase parsing
833
834 - ganymed: To get a nice FalseSource plot (instead of the all-information
835 debug plot displayed as a standard plot) use the context menu
836 somewhere between the two upper plots, click on "DrawNicePlot"
837 and enjoy the result. (It is the same context menu which contains
838 the other options)
839 PLEASE ALWAYS USE THIS PLOT FOR ANY KIND OF PRESENTATION!
840
841 - ganymed/optim/train: The cut program, trainings- and optimization
842 macros could run with less files found than expected from the
843 sequence-files. No the execution is stopped in this case.
844
845 - ganymed: now stores the dataset in the output file with all information
846 contained (as sequence number, sequence path, data path) so that the
847 used dataset file can be recreated complete from the output
848
849 - sponde: now also allows the uasge of random forest
850
851 - sponde: added a plot showing the cut efficiency
852
853 - sponde: Using root versions < 5.13/04 still resulted in wrong binomial
854 errors (eg. effecting the collection area) because the errors of the
855 two histograms were swapped in the calculation, which makes a difference
856 in the case the histograms contains weighted events. Aditionally the
857 totalarea was incorrectly taken into account, which results in too small
858 errors, which makes a difference at low statistics were the errors
859 or the collection area dominate.
860
861
862
863 *** Version 0.10.3 (2006/10/23)
864
865 - database: Implemented ObservationMode flag in websites.
866
867 - database: For sequences there is now also the stoptime available on the
868 websites.
869
870 - database: Removed ManuallyChanged flag from sequence-info websites,
871 as it is not needed anymore with the new sequence building algorithm.
872
873 - general: MRolke.[cc,h]: Modification of TRolke from root_v5.12.00b.
874 Class for calculation of confidence intervals (upper limits).
875 Use model 3 or 4 for upper limit calculation for MAGIC observations.
876 Please read the class instruction of TRolke for more information.
877
878 - general: The programs now return 0xfe if the requested resource file
879 doesn't exist.
880
881 - general: Errors comming from root itself are now not output to
882 stderr anymore but to our default (colored) log-stream.
883
884 - general: Changed the interface to TFormula such that now everything
885 is parsed by TFormula, which should allow to use phrases even
886 without paranthesis correctly parsed.
887
888 - callisto: replaced the AdHoc number of effective integration slices
889 in the digital filter by an average of the correct number, averaged
890 over all possible set of weights. This also gives correct numbers
891 for the lo-gains which were totally wrong before.
892
893 - callisto: drastically reduced memory consumption by not storing
894 all intermediate calibration results in memory anymore (we might
895 write them into a file in the future)
896
897 - ganymed: To get a nice Theta-Sq plot (instead of the all-information
898 debug plot displayed as a standard plot) use the context menu
899 somewhere between the two upper plots, click on "DrawNicePlot"
900 and enjoy the result. (It is the same context menu which contains
901 the "DrawAll" option)
902 PLEASE ALWAYS USE THIS PLOT FOR ANY KIND OF PRESENTATION!
903
904
905
906 *** Version 0.10.2 (2006/10/15)
907
908 - general: changed the building of the sequence files such, that
909 the source-/project name in case of wobble-sources is not the
910 source/project of the first data run anymore but the new generalized
911 wobble source/project
912
913 - general: Added some functions to MAstro to calculate sun and
914 moon properties
915
916 - general: Added two new macros (starvisday.C and starvisyear.C)
917 which allow to plot the visibility of several stars in one day
918 or of one star along one year
919
920
921
922 *** Version 0.10.1 (2006/10/11)
923
924 - general: HTML logging output was distorted... fixed.
925
926 - general: Changed the writing of the sequencefiles such that only
927 as many pedestal files as giving at least 1000 pedestal events which
928 are nearby the first calibration run are used as pedestal runs
929
930 - sponde: In the calculation of the collection area(s) and the
931 distribution for MOnte Carlo and estimated energy the error
932 calculation was wrong because root didn't take the errors
933 properly into account... fixed.
934
935
936
937 *** Version 0.10 (2006/10/10)
938
939 - general: all executables now check consistency of the root version
940 installed with the root version used for compilation. The return
941 code in this case is 255.
942
943 - general: The algorithm to build seqeunces has been redesigned
944 completly. It is now more flexible in terms of a resource file
945 resources/sequences.rc which can be used to partly overwrite
946 the default algorithm and defined the important part of the
947 algorithm. It also takes correctly into account several ways
948 of taking wobble-data. Currently the new algorithm doesn't take
949 the number of pedestal events anymore into account (so a sequence
950 could have much less than 1000 pedestal events and data runs are
951 no-more used instead)
952
953 - general: a lot of small changes to make the code compatible with
954 root 5.12/00 and gcc 4.1
955
956 - general: Fixed some warnings thrown if more warnings are switched
957 on in the compiler
958
959 - general: Added a function to MMath which returns the Median of
960 the distribution abs(y[i]-median) but at 68:32 division instead
961 of 50:50, which is somehow the median counterpart of the rms.
962
963 - general: MAstroCatalog allows now to read calatog files (as bsc5)
964 even if they are compressed with gzip (recommended compression
965 is with "-1")
966
967 - showplot:
968 + batch mode creation of image files now works with root 5.12/00
969 + added support for writing bmp (though it never produces
970 bmp which can be read by programs like xv or gimp)
971 + added support for xml
972 + added new option to start the display with a different size, by
973 --display-width=nnn or --display-height=nnn
974 + added options to define the canvas size when the display is started
975 --canvas-width=nnn or --canvas-height=nnn
976 the option can also be used for batch-mode conversion
977 + added an option which automatically fits the display size to the
978 desktop-size: --auto-size
979
980 - merpp: didn't recognize files with the extension .raw.gz - fixed.
981
982 - merpp: allow the reports to have 1000ms due to a bug in the camera
983 program.
984
985 - callisto: New swapped pixels have been found and corrected. For
986 details see MRawRunHeader::FixAssignment()
987
988 - callisto: Callisto now raises an error if one of the input files could
989 not be opened.
990
991 - callisto: In some rare cases it could happen that in the signal
992 extraction of the digital filter the extractor returned nonsense
993 results due a wrong initializion of a variable (this happens in
994 cases in which the extracted signal is exactly zero or the
995 digital filter couldn't slide because of a too small extraction
996 window, which could happen in the lo-gain extraction)
997
998 - callisto: The lo-gains have been marked as valid in some cases
999 when they got not extracted at all (deltsumlo=deltatimelo=0) because
1000 the allowed extraction range (fLoGainLast-fLoGainFirst) was
1001 smaller than the number of slices (fLoGainWindowSize) needed
1002 by the lo-gain extractor
1003
1004 - callisto: Updated absolute calibration constants (muon calibration)
1005
1006 - callisto: fixed a bug which caused the hi-gain saturation not
1007 to be handled properly in the spline (the position of the
1008 saturation was calculated incorrectly)
1009
1010 - callisto: fixed a bug which gave weird results if hi- _and_
1011 lo-gain could not be extracted properly. In some cases hi-
1012 and lo-gain have not been extracted at all and the random result
1013 got calibrated. In other cases the hi-gain was saturating and the
1014 lo-gain has not been extracted. For this a completly arbitraty
1015 scale factor was applied to the result from extracting the
1016 charge of the saturating hi-gain.
1017
1018 - callisto: improves handling of extracted signal if lo- and/or
1019 hi-gain could not be extracted properly. If no appropriate
1020 calibration of either the signal or the arrival time is possible
1021 the pixel is now markes as unsuitable instead of estimated.
1022
1023 - callisto: Fixed a problem in the spline extraction which gave
1024 weird results if the signal was extracted from the first (and
1025 highest) slice (as a workaround extraction from the first
1026 slice has been forbidden)
1027
1028 - callisto: Fixed a bug regarding the use of the QE which doesn't
1029 effect the result because the default is used anyhow.
1030 (The bug was there since 14/07/2005)
1031
1032 - callisto: Changed the default for fgLoGainStartShift to -2.4
1033
1034 - callisto: Fixed a bug which caused all hi-/lo-gain calibration
1035 factors to be always 10 for all pixels as soon as the first
1036 calibration update has happend.
1037
1038 - callisto: Fixed the hi-/lo-gain offset for the Digital Filter. It
1039 was off by more than one slice (which gives wrong values for
1040 pixels which have saturated the hi-gain) This was a problem
1041 mainly for timing studies, not so much for the automatic standard
1042 analysis. Also the values for the Spline have been changed slightly:
1043 + Digital Filter from 1.70 to 0.95
1044 + Spline from 1.39 to 1.30
1045 The numbers were sorted out by an extensive test in the
1046 test-database.
1047
1048 - callisto: Removed the tab showing the average calibrated arrival
1049 time. It never gave useful information.
1050
1051 - callisto: Now tab "HiLoCal" which shows the ratio between hi- and
1052 lo-gain pulse for data pulses if hi- and lo-gain have saturated
1053 and both were extracted because the hi-gain was above the
1054 LoGainSwitch.
1055
1056 - callisto: Now tab "HiLoOff" which shows the offset between hi- and
1057 lo-gain pulse position (after application of LoGainOffset, which
1058 means in the ideal case it is 0) for data pulses if hi- and lo-gain
1059 have saturated and both were extracted because the hi-gain was above
1060 the LoGainSwitch.
1061
1062 - callisto: After extensive tests a new threshold level for determination
1063 of the pulse position has been found at 50phe. This is due to the
1064 fact that for lower pulses the arrival time changes its average
1065 behaviour dramatically. It tends to the average of the extraction
1066 range and starts sticking to non-floating point numbers (5, 6, 7, ...)
1067
1068 - callisto: Added two more plots to the calibration part. The average
1069 extracted signal and the average extracted arrival time as
1070 extracted by the signal extractor (This can help debugging problems
1071 because it doesn't involve a fit)
1072
1073 - star: changed the fit for the effective on time such that initial
1074 values are calculated automatically now instead of using build
1075 in values. This makes the fit more indepedant of the underlaying
1076 rates. With a test of ~350 sequences in the test database it could
1077 be shown that the new fit gives the same result +/-1sek.
1078 The highest deviation was +5s the lowest -10s. You can now control
1079 the bahaviour with two resource values:
1080 MHEffectiveOnTime.FistBin: 3
1081 MHEffectiveOnTime.NumEvents: 12000
1082 Use FirstBin=1 and NumEvents=120 to fit the effective on time
1083 from the resulting gamma candidates
1084
1085 - star: The Center histogram in MHHillas is now by default displayed
1086 with the colz draw-option
1087
1088 - ganymed: now produces valid error codes if failed
1089
1090 - ganymed: MFMagicCuts now supports a hadronness cut. Switch it on
1091 with the option "hadronness" instead of "area" or "both" if
1092 you want to use both cuts. The hadronness cut includes a fixed
1093 cut in hadronness (index 8) and a fixed cut in size (index 9).
1094
1095 - ganymed: MFMagicCuts: fixed abug introduced 2005/7/14 which caused
1096 the M3Long of the anti source to be not properly evaluated
1097 (actually is was evaluated as a unique random number) in case
1098 of optimization with anti-theta cut switched on. This might have
1099 lead to improper, but not necessarily wrong, optimization results.
1100
1101 - ganymed: ganymed_onoff.rc, ganymed_wobble.rc - updated with comments
1102 how to use random forest
1103
1104 - ganymed: use the same spark cuts as in ganymed_wobble.rc as in
1105 ganymed_onoff.rc. Seems to be unintentional that they were different.
1106
1107 - sponde: Since the introduction of the possibility to define a
1108 normalization energy instead of using the integral the default
1109 was not using the integral anymore but a normalization energy of
1110 1GeV. It is now the integral again.
1111
1112 - sponde: Fixed a problem in the weighting which did not allow
1113 a MC slope of -1.
1114
1115
1116
1117 *** Version 0.9.6 (2006/05/24)
1118
1119 - general: changed the binning of the azimuthal profile of the camera
1120 such that it is aligned with the 60deg boundaries.
1121
1122 - merpp: better handling of problems with the TH, TD part of the
1123 CC-REPORT for files older than 200507190 and 200412210
1124 respectively
1125
1126 - merpp: improved handling of RECEIVERS-COM-ERROR in CC-REPORT
1127
1128 - merpp: implementd latest changes in starguider reports (sinc 9th May)
1129
1130 - callisto: if an arrival time is extracted outside the fadc range is
1131 was set to the boundary before. Now a random number is assigned
1132 instead to avoid unfortunate coincidences of arrival times just
1133 by this arbitrary value.
1134
1135 - callisto: added a new option to callisto.rc which allows
1136 to set the maximum arrival time difference allowed for the
1137 badpixel algorithm to determin whether a pixels belongs to a shower:
1138 MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 0.9
1139
1140 - callisto: changed the old value (0.5) for MaxArrivalTimeDiff
1141 to 0.9 as suggested by Stefan due to a recent study on the fraction
1142 of pixels surviving the image cleaning
1143
1144 - callisto: switched off error in case of "pulse is too much to the left"
1145
1146 - callisto: The extraction window is not adapted anymore. This means that
1147 the pulse is extracted from the full range. The draw-back is that the
1148 background (pedestal) will become slightly higher. Switching on/off
1149 the PulsePosCheck in the callisto.rc only effects filling the
1150 histograms.
1151 WARNING: If the pulse is too much to the left data is calibrated!
1152 Check the PulsePosCheck-PulsePosition AND the Inhomogeneity parameter
1153 in the database to make sure that your data is valid!
1154
1155 - star: Show also the number of stars correlated correctly by the
1156 starguider
1157
1158 - ganymed: implemented a new class (MHThetaSqN) which allows to use more
1159 than one off-source region in wobble-mode. To use it add the following to
1160 your ganymed_wobble.rc:
1161 + MJCut.NameHist: MHThetaSqN (switch the new feature on)
1162 + MHThetaSqN.NumOffSourcePos: 3 (define number of off-regions)
1163 + MHThetaSqN.DoOffCut: Yes,No (switch on/off the off-cut)
1164 + Cut1.ThetaCut: None
1165 + MHThetaSqN.SignificanceCutLevel: 2.0 (increase off-cut by 2.0/1.7)
1166
1167 - ganymed: in the case the latest report is older than a default (currently
1168 one minute) and the current report will be skip the starguider
1169 correction and calibration will be reset. The maximum age can be
1170 setup from ganymed.rc by (minutes):
1171 MPointingDevCalc.MaxAge: 1.5
1172
1173 - ganymed: A first implementation showing all size-bins for the theta-sq
1174 resp. alpha-plots in a single tab automatically
1175
1176
1177
1178 *** Version 0.9.5.1 (2006/05/05)
1179
1180 - general: Added a new option "zdiff" to the Makefile which runs
1181 "cvs diff" as "make diff" would do, but with compression
1182
1183 - general: accelerated MTFillMatrix (used for eample in all training
1184 and optimization macros) by skipping some obsolete calles in the
1185 eventloop the first and second loop could be accelerated by ~20%
1186
1187 - general: Added a missing feature in the MFilterList class which
1188 prevented MFEnergySlope from working correctly in trainenergy.C
1189
1190 - general: Accelerated the random forest training and usage a bit
1191
1192 - merpp: Adapted to new raw data file format version 6
1193
1194 - merpp: fixed a problem with reading the starguider reports
1195 after 6.3.2006
1196
1197 - callisto: At some times six pixels had to be exchanged manually in the
1198 data. If the pixels to be swaped are both not found, an error was raised
1199 and calibration was stopped. Now it is assumed they are intentially not
1200 mapped and the Error has been changed to a Warning.
1201
1202 - star: Added new image cleaning based on the arrival time of the shower.
1203 The new cleaning can be accessed using the parameter "Time" instead of
1204 the the old one (like "Absolute) in the star.rc-file.
1205 Recommended parameters:
1206 + MImgCleanStd.CleanLevel1: 8.2
1207 + MImgCleanStd.CleanLevel2: 5
1208 + MImgCleanStd.CleanRings: 2
1209 + MImgCleanStd.KeepSinglePixels: No
1210
1211 - star: Simplified the calculation of the effective on time. By letting
1212 the fit start in the third bin the result should become more stable.
1213 The fit quality is now ignored for the theta-plot which is not used
1214 in the analysis later. Still better pre-cut which, for example,
1215 remove events triggered by car-light (rate changing on a short
1216 time-scale) are missing.
1217
1218 - macros: optimonoff.C, optimwobble.C: implemented a new strategy
1219 for optimization (kWeakSource) which is suggested by the Berlin
1220 people for weak sources.
1221
1222 - macros: plotdb.C:
1223 + if no dataset and no time-interval is given only sequences with
1224 an effective on-time above five minutes are plotted anymore
1225 + added a new tab "RelTime" displaying the relative on-time
1226 (effective on-time divided by the observation time)
1227
1228 - macros: optim*.C: by skipping some obsolete calles in the eventloop
1229 the optimization became a lot faster
1230
1231 - ganymed: In the second loop the MHNewImagePar histograms disapeared.
1232 To display them just delete the pads displayed on top.
1233
1234 - ganymed: fixed a bug which could cause strange behaviour in wobble
1235 mode if no starguider information was available due to an incorrect
1236 initialization of the starguider calibration
1237
1238 - ganymed: The scale region in the False Source Plot (MHDisp) can now
1239 be changed from within the plot. It can also be setup from the
1240 resource file using:
1241 + MHDisp.fScaleMin: 0.325
1242 + MHDisp.fScaleMax: 0.475
1243
1244 - ganymed, showplot: Fixed a problem with the fit to the off-data
1245 (HistOff) which could cause a time lag or an inifinite loop.
1246 The bug did not effect the results at all.
1247
1248 - sponde: Added a new tab showing the effective collection area after
1249 trigger (including the software trigger)
1250
1251 - sponde: MMcSpectrumWeight has a new option which allows to define
1252 the scaling factor not by the integral of the events but at
1253 a defined energy:
1254 e.g. MMcSpectrumWeight.NormEnergy: 200
1255
1256 - MTFillMatrix, sponde: by skipping some obsolete calles in the
1257 eventloops the loop could be accelerated by ~20%
1258
1259
1260
1261 *** Version 0.9.5 (2006/03/13)
1262
1263 - macros: fixed a typo optimwoble.C (SetThetaCut instead SetAlphaCut)
1264 and added the standard Disp-parametrization
1265
1266 - macros: added a new macro datacenter/macros/plotstat.C allowing
1267 to plot statistics about the data processing
1268
1269 - macros: added a new macro to plot the muon correction coefficients
1270 versus period (macros/plot/mucal.C)
1271
1272 - showplot: path inflation for ganymed files was broken
1273
1274 - callisto: Fixed an important bug. When the number of saturating pixels
1275 were calculated the bad pixel information already available from the
1276 calibration was ignored. This yields wrong values in case pixels
1277 are broken and saturate all the time or randomly.
1278
1279 - callisto: Updated Muon calibration constants. Spline data now also uses
1280 different values for different periods
1281
1282 - star: The time-evolution shown in the effective on-time plot now displays
1283 the axis range (if there are no holes it should be similar to the
1284 effective on-time)
1285
1286 - star: Added two new image parameters:
1287 + MNewImagePar.fConcCOG: The ratio of the equivalent signals of
1288 the three pixels next to the center of gravity and the size
1289 + MNewImagePar.fConcCore: The ratio of the signals of all
1290 pixels inside or touching the ellipse and the size
1291
1292 - ganymed: now stops in requested source isn't found in catalog file
1293
1294 - ganymed: Many improvements to the phi-plot in wobble mode:
1295 + with a template the ideal background is displayed
1296 + a possible camera inhomogeneity is taken into account
1297 + the anti theta cut is better supported
1298 + the phi cut is automatically calculated from the theta cut
1299 See the examples in ganymed_wobble.rc for more details.
1300
1301 - ganymed: fixed plot for E^2*dN/dE. The spectrum was multiplied with two
1302 times the bin-width instead of the absolute energy
1303
1304 - ganymed: New cuts for on/off and wobble mode
1305
1306 - ganymed: The smearing of the False Source plot has been increased
1307 to 0.06
1308
1309 - ganymed: The starguider calibration can now be changed from the
1310 resource file.
1311 See the examples in ganymed*.rc for more details.
1312
1313 - ganymed: added a new option for on-/off-observations
1314 MJCut.RandomSourcePosition
1315 If switched on the distribution of observation time of the source
1316 position in the camera is calculated from the on-data and randomly
1317 applied with this distrbution to the off-data.
1318
1319 - ganymed: by skipping some obsolete calles in the eventloop the
1320 first and second loop could be accelerated by ~20%
1321
1322 - ganymed: The default in ganymed_onoff.rc is now to determin the source
1323 position in the off-data randomly according to its distribution in the
1324 on-data.
1325
1326
1327
1328 *** Version 0.9.4.3 (2006/02/13)
1329
1330 - general: Fixed a bug in most of the programs which prevented the
1331 compile date in the startup message to be displayed
1332
1333 - general: Simplified program calls. The following abbreviations
1334 are now possible:
1335 showlog callisto 55599
1336 showplot signal 55599
1337 star 55599
1338 callisto 55999
1339 ganymed 2
1340 In this cases the corrsponding root-, log-, sequence- or dataset-
1341 filename is automatically compiled using the default datacenter
1342 paths ("/magic/*")
1343
1344 - general: Fixed a long outstanding bug which made axis tick marks
1345 disappear if MDataFormula was used (eg. abs(x))
1346
1347 - general: Updated MMath with new functions to calculate the results of
1348 a exponential, logarithmic and powerlaw fits analytically.
1349
1350 - general: Updated some macros with comments:
1351 + macros/optim/optimdisp.C
1352 + macros/optim/optimenergy.C
1353 + macros/optim/optimwobble.C:
1354
1355 - general: MTFillMatrix (the class to fill one or two MHMatrix from
1356 files) now allows adding a pre-cut like in the optimization. E.g. this
1357 is useful to perform g/h-separation cuts before training the random
1358 forest.
1359
1360 - general: Between the rund 53300 and 68754 the pixels 553-558 were not
1361 correctly assigned in the raw data assignement table. A fix to correct
1362 the assignment tabel have been implemented.
1363
1364 - general: MAstroCatalog now allows rotating of the sky-grid and stars
1365 by 90, 180 and 270 deg.
1366
1367 - general: the macro plotdb.C now allows to highlight sequences belonging
1368 to a datset. For more details see the descritions in the macro.
1369
1370 - StatusDisplay: Fixed the algorithm to calculate the size ratio
1371 of the display, it gave wrong results espacially for big sizes.
1372 Now the embedded canvas should always have a ratio width/height=3/2.
1373 The "Size" options have been updated with values fitting the 3/2
1374 size ratio.
1375
1376 - RanForest:
1377 + Updated the random forest classes to support also the
1378 regression method implemented by Thomas H.
1379 + added new tutorial macro how to train the random forest
1380 for energy estimation (macros/optim/rfenergyest.C)
1381 + new classes to train the random forest (still in development)
1382 mjtrain/MJTrainEnergy, mjtrain/MJTrainDisp, mjtrain/MJTrainSeparation
1383 + new tutorial macros for random forest training in macros/train
1384 trainenergy.C, traindisp.C, trainseparation.C
1385
1386 - mars: now allows to setup a cut in the resource file mars.rc to choose
1387 the events to be displayed, eg:
1388 Cut.Condition: {0} && {1}
1389 Cut.0: MHillas.fSize>100
1390 Cut.1: MHillas.fSize<900
1391
1392 - merpp: Made aware of the new .rep format 200510250. This includes
1393 interpretation of the Receiver Board Temperatures (MCameraRecTemp),
1394 the discriminator thresholds (MCameraTD) and the disciminator
1395 delays (MCameraTD). All of them come with the CC report and are therefor
1396 written to the CC tree.
1397
1398 - callisto: fixed a bug in the calculation of the arrival times in
1399 the Digital Filter. They were all shifted by the same amount
1400 fHiGainFirst. This caused an offset in the order of the position of
1401 the left border of the extraction window in the times calculated
1402 from high-gain and from low-gain. The bug effected only the absolute
1403 arrival times. The bug was there since 30.1.2005.
1404
1405 - callisto: Changed the ExtractionWindow for the pulse position check
1406 such that the camera inhomogeneity becomes as less as possible and
1407 comparable to the pulso position check switched off. More details
1408 on the study will be given at our result web-page.
1409
1410 - ganymed: ganymed.rc, ganymed_onoff.rc:
1411 changed to new cuts trained with 1ES1218 and Crab-data
1412 for maximum sigma*log10(excess)
1413
1414 - ganymed: In addition to the Hadronness calculator (CalcHadronness)
1415 a new option was implemented to estimate Disp (CalcDisp)
1416
1417 - ganymed: Implemented two new options which allow
1418 to overwrite the default path to search for the sequence files
1419 --ins=/magic/sequences
1420 and where to search for the image data
1421 --ind=/magic/data/star
1422
1423 - ganymed: MSrcPosCorrect now allows overwriting of the misfocussing
1424 dx/dy from the resource file by
1425 MSrcPosCorrect.Dx
1426 MSrcPosCorrect.Dy
1427
1428 - ganymed: fixed reading the source catalog. Only as many chars
1429 as the requested source name had were compared (If a source
1430 "1ES1218" and "1ES1218+304" were in the catalog and "1ES1218+304"
1431 were requested it could happen that the coordinates of
1432 "1ES1218" were read) The source coordinates are now printed
1433 in addition to its name to the logging output.
1434
1435 - ganymed: changed the background fit in the case of a ThetaSq function
1436 from a polynomial to a exponential, this should give a big improvement
1437 for on-only analysis'
1438
1439 - ganymed/sponde: Calculation of the error of single size-/energy-
1440 bins now uses Li/Ma (5) instead of LiMa (17)
1441
1442 - sponde: Added a plot E^2*dN/dE
1443
1444 - sponde: The energy estimator plot should now show values like
1445 they are commonly used.
1446
1447 - sponde: Now MMcSpectrumWeight also excepts formulas with two X
1448 (a powerlaw with cutoff didn't work before)
1449
1450
1451
1452 *** Version 0.9.4.2 (2005/10/01)
1453
1454 - general: added new tutorial macro macros/tutorials/threshold.C:
1455 how to calculate the threshold
1456
1457 - general: added a new game. Start it from the interpreter with
1458 MagicJam j; (other games are: MineSweeper, MagicSnake, MagicDomino
1459 and MagicReversi)
1460
1461 - general: neighbor pixels in MGeoms are now sorted clockwise
1462
1463 - merpp: made aware of new format 20050829-0
1464
1465 - callisto: New scale factors from muon analysis, stored in new
1466 resource file resources/calibration.rc
1467
1468 - callisto: MBadPixelsCalc now raises an error instead of simply stopping
1469 the eventloop if something went wrong. This is necessary for the
1470 automatic processing
1471
1472 - callisto: Implemented a new treatment for the time information
1473 developed by S.Ruegamer, this should give much better response than
1474 the old very primitive algorithm. This is especially important
1475 for bad pixels due to bright stars as Ceta-Tauri (btw: the old
1476 treatment was buggy! As neighbor pixels to interpolate the
1477 arrival times always the pixels 0 to 5 were taken)
1478
1479 - callisto: In addition to the new treatment of arrival times the
1480 decision of the treatment whether a pixel contains shower
1481 information or not is used to improve the signal interpolation
1482 further. Previously interpolated pixels showed a non-zero signal
1483 do the peak searching signal extractors. Having a signal pixel
1484 beside raised the signal above the cleaning level. Thus the
1485 pixel survived cleaning to often but with a too low mean signal.
1486 By replacing empty pixels with a gaus (mean and rms from the
1487 corresponding Area entries in MPedPhotCam) this effect can
1488 be supressed.
1489
1490 - callisto: The status "unmapped" (for pixel which cannot be interpolated)
1491 got lost when the events were stored by the calibration... fixed.
1492 (This was no problem for normal image cleaning because the pixel
1493 content was artificially set to 0 not surviving standard cleaning
1494 algorithm)
1495
1496 - callisto: Fixed a bug which caused the random and peak-search pedestal
1497 extracted with an extractor to be exchanged for the first (roughly)
1498 500 events (*Extractor and *ExtractorRndm).
1499 (They were simply exchanged in callisto.cc) This bug might have been
1500 introduced in Mars 0.9.4 when the order of calculation of the two
1501 types of pedestal in callisto got exchanged. The bug only effects
1502 the first seconds of data of each sequence.
1503
1504 - star: Taking the arrival time of Muons into account the calculation
1505 of the parameters could be made independant of the noise
1506 (pedestal rms). The dependance was due to the bias of the peak search
1507 extractor. Now only a ZA dependamcy is left which is not yet fully
1508 understood.
1509
1510 - ganymed: The sequences were still not sorted correctly in any case,
1511 because they were sorted by the full qualified path name. Now they
1512 are sorted by their sequence number.
1513
1514 - ganymed: new values for the spark-event-cut
1515
1516 - sponde: required missing MEnergyEst from ganymed.root... fixed.
1517
1518
1519
1520 *** Version 0.9.4.1 (2005/08/18)
1521
1522 - callisto: fixed some problems with the calibration in case of
1523 inteleaved events. Therefor the final and some intermediate
1524 fits are skipped which would take place on partly filled #
1525 histograms and gave wrong results or failed completely.
1526
1527 - callisto: If the intermediate finalization of the histograms
1528 calculating the mean charge of the calibration signal fails
1529 it is counted now and printed in PostProcess of MCalibCalcFromPast.
1530
1531 - ganymed: Fixed some bugs which caused problems in On-only mode.
1532 Still the false source plot doesn't give reasonable results.
1533
1534
1535
1536 *** Version 0.9.4 (2005/08/05)
1537
1538 - general: Fixed the ZA binning. It did not correctly fit the
1539 MC binning
1540
1541 - general: Added a class to correct for the missfocussing in
1542 April (MSrcPosCorrect)
1543
1544 - general: Random Forest energy estimator should now work properly
1545
1546 - general: MHCamera now displays the profiles in deg instead of mm
1547
1548 - general: MH::SetPalette offers a lot of new palettes
1549
1550 - general: MHillas - the case of CorrXY==0 is now handled properly
1551
1552 - general: implemented the possibility to change the line and
1553 marker style of a sky-grid drawn by MAstroCatalog
1554
1555 - general: Runs in sequences and datasets are now automatically
1556 ordered (it is assumed that the order of their filenames correctly
1557 represents the order of observations) to ensure correct reading
1558 of the subsystem data
1559
1560 - general: Runs cannot be added to a sequence twice anymore. Also
1561 Sequences cannot be added twice anymore to a dataset.
1562
1563 - general: fixed a bug in MAstro::Deg2Hms reported by Abelardo.
1564 It gave wrong results by 1./2pi
1565
1566 - showplot: got a new option to start a root interpreter, too
1567
1568 - mars: now displays a rough estimate of Disp and the third moment
1569
1570 - mars: now displays the event time MTime if available
1571
1572 - mars: show muon parameters graphically
1573
1574 - mars: now the file to open can be given as commandline
1575 argument
1576
1577 - merpp: Added support for new arehucas version >= 200507190
1578 implementing values for three temperature sensors in the
1579 counting house (MCameraAUX) and the name of the HV settings file
1580 (MCameraHV).
1581
1582 - merrp: Also fixed some problems with wrong versions numbers stated
1583 by arehucas and implemented some fixed for old versions (<200407070)
1584
1585 - callisto: MCalibrationHiLoCam can now be printed from its context
1586 menu, eg in the TBrowser
1587
1588 - callisto: fixed logain offset (fgOffsetLoGain) from 1.7 to
1589 - 1.39 (MExtractTimeAndChargeSpline)
1590 - 1.40 (MExtractTimeAndChargeDigitalFilter)
1591 This is important mainly for timing studies.
1592
1593 - callisto: Changed limits in MHCalibrationChargeCalc from
1594 - -100.125 to -98 (fgChargeHiGainFirst)
1595 - 1899.875 to 1902. (fgChargeHiGainLast)
1596 - -100.25 to -99 (fgChargeLoGainFirst)
1597 - 899.75 to 901. (fgChargeLoGainLast)
1598 Introduced new limits:
1599 - fgNumLoGainBlackoutLimit: 0.05
1600 - fgLoGainBlackoutLimit: 3.5
1601 - fgLoGainPickupLimit: 3.5
1602
1603 - callisto: The handling of automatic weights for the digital filter
1604 has been improved. There are now new calibration weights special
1605 for Monte Carlo calibrations.
1606
1607 - callisto: Updated Muon scale with latest results (1% difference
1608 for most data, 5% for Mar/Apr'05)
1609
1610 - callisto: use a new criterium to exclude bad pixels: If the high-gain
1611 was saturated and the blackout-events in the low-gain exceed the
1612 fNumLoGainBlackoutLimit, the pixel is declared unsuitable.
1613 This excludes those pixels which have a saturating high-gain
1614 channel, but the low-gain switch does not switch often enough
1615 to make the distribution reliable.
1616
1617 - callisto: fix a bug counting the number of saturated events. Up to now,
1618 the number of saturated slices was counted (which is one for a not too
1619 high number), but for some (pathological) pixels, many more slices
1620 saturated and produced wrong limits.
1621
1622 - callisto: New options in in callisto.rc for MCalibrationRelTimeCalc:
1623 + MCalibrationRelTimeCam.CheckFitResults: Yes
1624 + MCalibrationRelTimeCam.CheckDeviatingBehavior: Yes
1625 + MCalibrationRelTimeCam.CheckHistOverflow: Yes
1626 + MCalibrationRelTimeCam.CheckOscillations: Yes
1627
1628 - callisto: introduce max. number of events for intercalibration,
1629 the rest gets skipped. Default: 4096
1630 The calibration causes too many un-reliable pixels if more
1631 than about 5000 events are treated (@500 Hz) because of the
1632 mode hopping of the VCSels. However, in the past, some
1633 calibration runs have been taken (erroneously) with more
1634 than 5000 events, especially the intensity scans where
1635 a good precision is needed.
1636
1637 - callisto: Set UseExtractorRes to kTRUE, will raise the number of
1638 photo-electrons for the spline extractor by 1-2% (correctly).
1639
1640 - callisto: The intercalibration is now done using a full f-factor
1641 method instead of a simple update. This makes things slower but
1642 more reliable in terms of gain-fluctuations and sudden VCSEL gain
1643 changes.
1644
1645 - callisto: MPedestalY3 now uses the same extraction range as the data
1646 as calculated by MPedestalY2
1647
1648 - callisto: A check of the fraction of allowed unsuitable (10%) and
1649 unreliable (30%) pixels has been addded when the conversion is
1650 calculated from interleaved events
1651
1652 - star: fixed a bug which caused MEffectiveOnTime containers not to
1653 be written. In the worst case they were skipped almost all, in the
1654 best case they were all ok. It depends on how often the new value
1655 was assigned just with an event which was skipped later by another
1656 task
1657
1658 - star: new histograms MHPointing and MHRate
1659
1660 - star: New tab (Currents) displaying the time evolution of the mean
1661 DC currents
1662
1663 - star: New tab (MeanRms) displaying the time volution of the mean
1664 Pedestal Rms
1665
1666 - star: Output now contains MReportStarguider
1667
1668 - star: Improved muon analysis. Due to the absense of a constant
1669 for the gaus fit for arc width and the fact that the arcwidth
1670 plot was a histogram instead of a profile the resulting psf
1671 became pedestal rms dependant
1672
1673 - ganymed: added new plot (MHPhi) for wobble mode as described in
1674 http://www.astro.uni-wuerzburg.de/results/ringmethod/
1675
1676 - ganymed: corrects now (correctly?) for the misspointing in
1677 May/April 05
1678
1679 - ganymed: Two new default templates for ganymed.rc are avilable:
1680 ganymed_onoff.rc and ganymed_wobble.rc
1681
1682 - ganymed: The first version of MFMagicCuts has been released
1683
1684 - ganymed: the Conc1 plot was incorrectly scaled in MHVsSize
1685
1686 - ganymed, sponde: There are two new options to determin the number
1687 of bins for the signal region (NumBinsSignal) and the number of
1688 total bins (NumBinsTotal) in the MHThetaSq histogram
1689
1690 - ganymed: optimized palettes for MHDisp
1691
1692 - ganymed: Fixed an important bug when using ThetaSq as fit function:
1693 The result for the sigma of the gauss(sqrt(x)) fit was a factor
1694 sqrt(2) to large
1695
1696 - ganymed: The false source plot (MHDisp) is now based on Disp
1697 and a background model determined in the first loop is
1698 subtracted
1699
1700 - ganymed: MFMagicCuts are now available in the cvs. The correct
1701 solution for the disp ambiguity problem is now calculated using
1702 the third moment. The third moment cut therefor is skipped.
1703
1704 - ganymed: Output now contains "Disp [MParameterD]" which sign is
1705 already determined by the third moment along the major axis and
1706 so its actual sign is aligned along the major axis of the shower.
1707
1708 - ganymed: If starguider data is available it is taken into account
1709 when calculating the source position. This is especially important
1710 for wobble data at high zenith angles and culmination data.
1711
1712 - ganymed: If available the starguider data is taken into account as
1713 good as it is possible without a starguider calibration. To switch
1714 it off add
1715 MPointingDevCalc.MaxAbsDev: -1
1716 to your ganymed.rc
1717
1718 - sponde: the zenith angle distribution is now weighted instead of
1719 choosen by default. This gurantees higher statistics in the highest
1720 and lowest energy bins. To switch this behaviour off use...
1721
1722 - sponde: now automatically displays the energy threshold
1723
1724
1725
1726 *** Version 0.9.3 (2005/06/03)
1727
1728 - general: MMcEvt now derived from MMcEvtBasic which should
1729 have no influence on compatibility with older camera files
1730
1731 - general: Implemented a new option to MReadReports which allows
1732 to force reading of the first event of some trees (eg. Drive-tree)
1733 before any other tree.
1734
1735 - general: Made a fix to MPointingPosCalc which assumes that the
1736 run-type is real-data if Process is called before ReInit, such
1737 that also the very first Drive-events in a run are treated
1738 correctly.
1739
1740 - shoplot now displays the filename in the title bar of the window
1741
1742 - merpp: new options: --auto-time-start --auto-time-stop which
1743 allow to cut updating with slow-control at the start- or
1744 (previousle --auto-time only allowed: and) stop-time of
1745 the run stored in the runheader.
1746
1747 - callisto: changed default of fgFallTimeHiGain in
1748 MExtractTimeAndChargeSpline from 1.5 to 0.5. 0.5 was
1749 already set as default in callisto_Dec04Jan05.rc
1750
1751 - callisto: the digital filter (MExtractTimeAndChargeDigitalFilter)
1752 now changes the filename for new weights automatically for
1753 cosmics, MC and for different calibration colors.
1754
1755 - callisto: a problem with the order of some code cause the
1756 name of the weights file in MPedestalY2 to be ignored.
1757 Wrongly the calibration_UV weights were used instead
1758 of the cosmics weights set in callisto.rc.
1759
1760 - callisto: should allow ct1 pulser data when the color is
1761 explicitly set also for pretty new runs
1762
1763 - callisto: a new option (MCalibrateData.ScaleFactor) allows
1764 to introduce an additional scale factor (should be around 1)
1765 to adapt for slight deviations in different extractors from
1766 the MC calibration. For MC calibrations the scale factor
1767 is fixed to 1
1768
1769 - callisto: MFCosmics should handle very high signals a bit
1770 more accurate now
1771
1772 - callisto: when calculating the mean pedestal rms for different
1773 pixel sizes or sectors also kUnsuitableEvt is now taken into
1774 account - not only kUnsuitableRun
1775
1776 - callisto: '-raw' option is now the default
1777
1778 - callisto: A dead pixel check has been implemented in MJPedestalC2.
1779 This has the advantage that mean values in the calibration (eg
1780 arrival times) are calculated more accuratly.
1781
1782 - callisto: In the resource file callisto_Dec04Jan05.rc
1783 MJPedestalY2.ExtractWinRight has been reduced from 4.0 to 2.0
1784
1785 - callisto: new Hi-/Lo-Gain intercalibration constants
1786 hilocalib_sp1.root and hilocalib_sp1_mc.root
1787
1788 - callisto: changed default for MJPedestalY2.MaxEvents
1789 from 2000 to 5000 like in callisto_Dec04Jan05.txt
1790
1791 - callisto: in MCalibrationChargeCalc the limit fgPheErrLowerLimit
1792 has changed from 9.0 to 6.0. This is necessary to exclude two very
1793 ringing pixels in September which - when kept - destroy the image
1794 cleaning around them (Markus G.)
1795
1796 - callisto: (!) The Muon analysis of all data taken so far has
1797 suggested that different extractors give slightly different
1798 absolute results after calibration, due to there different
1799 handling of the calibration- and cosmics pulse shape. Although
1800 the over all calibration seems to be a bit to pesimistic (or
1801 the MC seems to be a little bit - 3-5% - too much attenuated).
1802 There for correction factors have been introduced which correct
1803 the calibration constants by
1804 0.956937799 for Digital Filter (callisto.rc)
1805 and
1806 0.904159132 for Spline (callisto_Dec04Jan05.rc)
1807 calibration. The main advantage is, that both data types can
1808 be easily analysed together and compared to the same MC sample.
1809 The factors are choosen such, that the reference is now a
1810 Monte Carlo sample calibrated with Digital Filter (callisto.rc)
1811
1812 - star: fixed a bug which caused the first bin of the histograms
1813 vs time to disapear in PostProcess. This gave wrong results
1814 on the screen. The output and all following programs not
1815 directly accessing the histograms haven't been affected
1816
1817 - star, ganymed: Made sure that always the first Drive-tree event
1818 is read and Processed before the first data event was read.
1819
1820 - ganymed: implemented the possibility to replace the standard
1821 false source plot with a different one. For example: If
1822 you want to use a Disp based FS plot use the resource:
1823 NameHistFS: MHDisp
1824 For details about the requirements of these plots look
1825 at the corresponding classes
1826
1827 - MJOptimize: For On/Off and Wobble optimization MJOptimize
1828 has been replaced by MJOptimizeCuts. To use MHThetaSq instead
1829 of MHAlpha call MJOptimizeCuts("MHThetaSq"). Make sure in this
1830 case MAlphaFitter is setup correctly. For On-Only optimization
1831 MJOptimize has been Replaced by MJOptimizeCuts, too. Call
1832 RunOn() instead of Run(). For energy estimation MJOptimize
1833 has been replaced by MJOptimizeEnergy
1834
1835 - sponde: The input MC spectrum can now be weighted to fake a
1836 different spectrum. This is done via MMcSpectrumWeight. For
1837 more details see the class description and sponde.rc
1838
1839 - sponde: The paremeter comparsion plots are not scaled by
1840 their entries anymore. Instead the MC plot is scaled by using
1841 the result spectrum of the analysis. If the input MC spectrum
1842 and the result spectrum has different slopes the absolut
1843 normalization is normally wrong.
1844
1845
1846
1847 *** Version 0.9.2 (2005/05/20)
1848
1849 - general: Changed the default paths for calibrated data and image
1850 files. (The implemented access to these files doesn't yet exist)
1851
1852 - general: support for MRunIter has been removed from the job classes
1853 (use the setter functions of MSeqeunce instead)
1854
1855 - general: ProcessFile has been renamed to Process in all job classes,
1856 because ProcessFile is missleading
1857
1858 - general: The storage of the palette in MHCamera is now possible
1859
1860 - general: new bits describing unsuitable and unreliable pixels have
1861 been introduced in MBadPixelsPix: kPreviouslyExcluded,
1862 kDeadPedestalRms, kFluctuatingArrivalTimes
1863
1864 - general: fixed reading of CC-REPORTS (weather data) for Arehucas file
1865 versions newer than V200404070
1866
1867 - general: support for cc file versions newer than V200504130. Added
1868 support for active loads (MCameraActiveLoad) and central pixel
1869 (MCameraCentralPixel)
1870
1871 - general: added new histograms for image parameters versus size
1872 (MHVsSize)
1873
1874 - general: The eventloop now automatically calls PrintStatistics
1875
1876 - Status Display: Improved support for printing:
1877 + A default can now be set in .rootrc (for more details see
1878 MStatusDisplay::PrintPS)
1879 + A print dialog is opened so that printer and command line
1880 can be canged
1881
1882 - showlog: improved such that it handles different verbosity levels
1883 more accurate now. This mean you can suppress output when
1884 using showlog.
1885
1886 - showplot: improved
1887 + The new graphics formats have been implemented
1888 (pdf, svg, png, jpg and xpm)
1889 + Implemented an interface for printing. Because you can change
1890 the printing command from the command line you can use all
1891 the nice postscript tool. To print a booklet for example call:
1892 showplot -b --print --print-cmd="cat %f" filename.root | lpr
1893 showplot -b --print --print-cmd="psbook %f | psnup -2 | lpr" filename.root
1894 showplot -b --null --print --print-cmd="psbook %f" filename.root > book.ps
1895
1896 - callisto: fixed a bug in the recognition of the calibration bit
1897 pattern. Concerns only intensity calibration. (The strength of the
1898 intensity has been updated in some cases)
1899
1900 - callisto: add a filter against PIN Diode events when running over
1901 calibration runs since some of the recent calibration runs have
1902 been taken with Pin Diode.
1903
1904 - callisto: added support for the runs 39942, 39944, 44834, 39941,
1905 39943, 44833 in the calibration (MCalibColorSet)
1906
1907 - callisto: MJCalibration.MHCalibrationChargeCam.ProbLimit has
1908 been set to 1e-18 in callisto_Dec04_Jan05.rc
1909
1910 - callisto: in callisto_Dec04Jan05.rc: set value for
1911 MJPedestalY2.MaxEvents from 2000 to 5000 to get better
1912 statistics for the calculation of the mean pulse time
1913 of the cosmisc
1914
1915 - callisto: the charge limit in the calibration has been raised
1916 from 2.5 to 4.5
1917
1918 - callisto: new setup read from callisto.rc
1919 + MCalibCalcFromPast.fNumEventsDump
1920 + MHCalibrationChargeCam.NumHiGainSaturationLimit
1921 + MHCalibrationChargeCam.NumLoGainSaturationLimit
1922
1923 - callisto: MExtractTimeAndChargeSpline extractor:
1924 + The integration ranges were shifted by 0.2 FADC slices to the right.
1925 + unreliable handling of floating point values caused different
1926 results on different platform or with different compiler options:
1927 The new algorithm is numerically exact and more accurate because
1928 it doesn't calculate the integration position by multiple times
1929 adding the step-size (which results in numerical uncertanties
1930 exspecially if multiplied with large numbers)
1931 + A lot of fixes have been introduced which effects integrating the
1932 spline at the edges of the valid range. In this case any memory
1933 was randomly accessed. This behaviour can be tested replacing
1934 all MArrayF by TArrayF which perform a range check (but of course
1935 they are a lot slower)
1936 ! No result obtained with the Spline before can be trusted! Due to
1937 random memory access it might by completely random!
1938
1939 - callisto: set new defaults in MExtractTimeAndChargeDigitalFilter:
1940 + extraction range in the high gain to 4.
1941 + weights file to cosmics_weights46.dat
1942
1943 - callisto:
1944 + introduce the check for fluctuating arrival times.
1945 + set charge limit to from 2.5 to 4.5
1946 + introduce asymmetric exclusion due to "deviating number of phes":
1947 Default: -9.0 +4.5
1948 + make deviating number of phes unsuitable (see David's findings).
1949
1950 - callisto: fixed the display of unsuitable pixels. The old change
1951 to a real profile in MHCamera broke it.
1952
1953 - callisto: The cosmics filter is used in MJCalibration now
1954 independant of the pulser color (formaly it was used only for
1955 CT1-pulser data). The threshold of number of max empty pixels
1956 has been raised from 2% to 5% in MJCalibration. The maximum
1957 fraction of rejected event is 50% otherwise the calibration has
1958 failed. The filter now ignores unsuitable pixels.
1959 Setup the filter from callisto.rc using "ContCosmics"
1960 (see MFCosmics::ReadEnv)
1961
1962 - callisto: changed default for fgOffsetFromLogain back
1963 from -2.8 to -1.8 in digital filter.
1964
1965 - callisto: for the common extractors
1966 + MExtractFixedWindow
1967 * MExtractTimeAndChargeSlidingWindow
1968 * MExtractTimeAndChargeSpline
1969 * MExtractTimeAndChargeDigitalFilter
1970 the hi- and lo-gain reolution is stored in fResolutionPerPheHiGain
1971 and fResolutionPerPheLoGain
1972
1973 - callisto: Updated the hi-/lo-gain intercalibration constants
1974 hilocalib_sp1.root, hilocalib_df46_mc.root, hilocalib_df4.root,
1975 hilocalib_df6.root
1976
1977 - callisto: The default lo-gain extraction window for the
1978 MExtractTimeAndChargeSlidingWindow extractor has been increased
1979 from 6 to 8 to take the stretch of the pusle into account
1980
1981 - callisto (!!!): fixed an important bug in the setting of the
1982 incoming pedestal pointer to the signal extractor extracting
1983 pedestals. This bug caused the pedestals not to be subtracted
1984 correctly in the calibration afterwards. The effect was about 10%
1985 too high conversion factors for the digital filter and the spline
1986 and more than 40% for the sliding window.
1987 The bug existed since end of February! (Markus G.)
1988
1989 - star: added muon support to star. A new tab "MHMuonPar" is
1990 displayed. The lower right plot is an estimate of the point spread
1991 function. A reference value is displayed. To get the real PSF
1992 you must compare with MC. The setup for the muon analysis
1993 is done in star.rc. More informations can be found in
1994 MHSingleMuon and MHCalibParCalc
1995
1996 - star: The weather data is now displayed in star
1997
1998 - star: MEventRateCalc handles the calculation of the event rate more
1999 accurate now in case of the start of a new run inside a sequence
2000
2001 - ganymed: a bug in MSrcPosCalc has been fixed (The changes were taken
2002 from a bugfix of Wolfgang to Loc0LocToCam; Wolfgang: the results are
2003 changed only minimally
2004
2005 - ganymed: (MJCut and MJOptimize) now displayes the number of
2006 excess events versus size. The energy estimation is done in
2007 MJSpectrum (sponde)
2008
2009 - ganymed: got support for using other variables than Alpha, eg. Theta.
2010 Therefor you need a class deriving from MHAlpha which supports
2011 this variable (one is already existing: MHTheta, which is now the
2012 default in ganymed). It is setup through ganymed.rc
2013
2014 - added a full featured spectrum program (sponde) which reads
2015 ganymed output and monte carlos and compiles a spectrum
2016
2017
2018
2019 *** Version 0.9.1 (2005/04/14)
2020
2021 - Fix bug on the arrival time reconstruction when using Digital
2022 Filter (in some marginal cases (10%) times were reconstructed shifted
2023 by 0.2 slices)
2024
2025 - all executables now return 2 where in previous versions -1 has
2026 been returned (gave problems with shell-scripts)
2027
2028 - the Hi-/Lo-Gain intercalibration for monte carlo files is now
2029 also done with a intercalibration input file
2030
2031 - a bug in the time fine adjustment of the digital filter has been fixed
2032
2033 - a lot of small fixes to calling TPad::Divide with a margin==0
2034
2035 - MStatusDisplay now shows the memory usage
2036
2037 - added a memory leak statistic after the eventloop
2038
2039 - MHadronness replaced by MParameterD
2040
2041 - MEnergyEst replaced by MParameterD
2042
2043 - wobble mode can now be switched on from the data-set
2044
2045 - switched on bad pixel detection for a pedestal level
2046 which is more than 5 times the variance
2047
2048 - Fixed default paths to be correct for the new data center
2049
2050 - Fixed handling of drive data. The pointing position in now
2051 copied to MPointingPos inside star for real and MC events
2052
2053 - made sure, that the sequences in a dataset are sorted by
2054 their number to ensure the correct synchronization
2055
2056 - added all resource files necessary to handle spline extraction
2057 of signal for the late 2004 and the early 2005 data
2058
2059 - update the reference lines for the december data
2060
2061 - changed default for fgHiGainFirst from 2 to 0 in spline extractor
2062 This is important for the Dec. and Jan. data where the calibration
2063 pulse position changes much.
2064
2065 - changed default for fgOffsetFromLogain from -1.8 to -2.8 in digital
2066 filter. This corrects only an defective extraction of the MC at
2067 low-gain signals lower than 90 phes.
2068
2069 - fixed some slowdown in calibrating data with interleaved events from
2070 Jan and Dec. As a workaround made the pretty palette the default.
2071 The inverse Deep Blue Sea has a bug which slows down long calibration
2072 runs a lot.
2073
2074 - implemented new values from new starguider report as defined
2075 in TDAS00-07 V7.1 (eg SkyBrighness)
2076
2077 - Fixed the formerly wrong unit of azimuth and zenith angle deviation
2078 in starguider report (If you merpp the same reports with this version
2079 you'll get the correct units (arcmin) but a different result
2080 than with old versions (wrong units))
2081
2082 - Now the status display can also write PNG, JPEG and XPM (other file
2083 formats are still unsupported by root)
2084
2085 - added a lot of scripts and macros in directory datacenter used
2086 for automatic processing
2087
2088
2089
2090 *** Version 0.9 (2005/03/23)
2091
2092 - Calibration is now done autmatically into photo-electrons like
2093 decided by the Collaboration Board
2094
2095 - Reference lines in the display can be read from resource file
2096
2097 - Interlaced calibration events can be handled by callisto
2098
2099 - Hi-Gain vs. Lo-Gain calibration constants can be calculated
2100 and applied in callisto.
2101
2102 - class MJCalibTest can be called by callisto
2103
2104 - callisto can now run with three data types: rawdata, rootdata, MC
2105
2106 - added two classes mreport/MReportStarguider.cc and
2107 mreport/MReportStarguider.h for implementing starguider-data
2108
2109 - the output of merpp-update now contains starguider information:
2110 MReportStarguider and MTimeStarguider
2111
2112 - Fixed bug in MHHadronness.cc which made that the displayed Q-factor
2113 values were lower than they should. This was more visible for high
2114 energy (Size) events.
2115
2116 - Made MExtractTimeAndChargeDigitalFilter the default for callisto
2117
2118 - added a new loop to callisto to support the correct extraction of
2119 the pedestal and pedestal rms through the signal extractor
2120
2121 - added class to perform fourier tranformation of a hexagonal structure
2122 (MHexagonalFT)
2123
2124 - added class to perform fourier tranformation of a hexagonal structure
2125 in a tasklist (MHexagonalFTCalc) example: hft.C
2126
2127 - added class to enhance the MCerPhotEvt from the Magic camera geometry
2128 MGeomCamMagic to a enhanced geometry MGeomCamMagicXT having only
2129 small pixels, example: extendcam.C
2130
2131 - added possibility to write data to memory (TTree) using MWriteRootFile
2132
2133 - added possibility to read a TTree stored only in memory by MReadTree
2134
2135 - new classes to deal with the interleaved calibration events:
2136 MTriggerPattern, MTriggerPatternDecode and MFTriggerPattern
2137
2138 - implemented filters to callisto to skip calibration events
2139
2140 - changed callisto to calculate three kinds of pedestals and pedestal rms:
2141 + fundamental pedestal (MPedPhotFundamental): extracted directly from
2142 the slices by statistics
2143 + pedestal from extractor (MPedPhotFromExtractor): extracted with the
2144 signal extractor without randomization
2145 + random pedestal from extractor (MPedPhotFromExtractorRndm): extracted
2146 with the signal extractor with randomization (extraction window is
2147 randomly fixed)
2148 The three kind of pedestal is not yet implemented (for speed reasons) as
2149 long as nobody really needs it:
2150 + pedestal from extractor with randomizatiuon (MPedPhotFromExtractorRndm):
2151 extracted with the signal extractor with randomization
2152
2153 - callisto now finds the pulse position in the data automatically
2154
2155 - implemented data reduction in calibrated data by reducing the precision
2156 of the stored floating point values (pedestal, rms and signal) to
2157 roughly 0.5%
2158
2159 - fixed showplot to make file format conversion (eg. converting a
2160 status display stored in a root file into a postscript file) in
2161 batch mode possible:
2162 showplot -b --save-as-ps[=outfile.ps] infile.root
2163
2164 - Allow class names as part of filter rules, eg
2165 (MFSupercuts) && (MHillas.fSize>100)
2166
2167 - fixed a lot of stuff to allow all programs to process MonteCarlo
2168 camera files.
2169
2170 - enhanced possibility to debug reading of resource file
2171
2172 - MDataChain should handle priorities of operator more
2173 accurate now
2174
2175 - completely revised the container in which the signal is stored
2176 MCerPhotEvt/Pix, MArrivalTime --> MSignalCam/Pix
2177 * the new implementation is faster and
2178 * needs less storage space
2179
2180 - FIRST version which can process Monte Carlos through the
2181 whole chain of programs (callisto, star, etc)
2182
2183 - implemented new image parameters for the border length of
2184 the shower image (MNewImagePar2)
2185
2186 - implemented wobble-mode in ganymed
2187
2188 - new class (MJOptimize) taking care of all optimization processes
2189
2190
2191
2192 *** Version 0.8.6 (2004/10/05)
2193
2194 - New combined signal and time extractor
2195 MExtractTimeAndChargeDigitalFilter
2196
2197 - The Report for the trigger section in the CC run summary has been
2198 updated, and now the IPRs (Individual Pixel Rates) are available.
2199
2200 - new macros for database interaction: filldotrun.C, filldotrbk.C
2201
2202 - added new program which can convert colored output which was redirected
2203 into a file back into colored text (for example to display it with
2204 more): showlog
2205
2206 - added new program to display an MStatusArray of plots (stored with
2207 MStatusDisplay) from a file: showplot
2208
2209 - Added support for Sequences (MSequence). Sequences are blocks of run
2210 files containing the same calibration and belonging to the same
2211 source
2212
2213 - Added a new task (MTaskEnv) which can be used in setup files to switch
2214 between different tasks (for example to choose a signal extractor)
2215
2216 - Implemented automatic file splitting in MWriteRootFile
2217
2218 - After image cleaning an island index is assigned to all used pixels.
2219 The index corresponds to the order of the islands in size.
2220
2221 - added a cleaning which takes signal height _and_ arrival time into
2222 account: probability cleaning (for more details see MImgCleanStd)
2223
2224 - added a cleaning which takes takes the absolute number of photons
2225 times the area ratio as cleaning level
2226 (for more details see MImgCleanStd)
2227
2228 - Changed calculation of image parameters such, that it can be done
2229 for all used pixels or the pixels with a defined island index only
2230
2231 - implemented new image parameters displaying the number of islands,
2232 saturated hi-gain and lo-gain pixels (MImagePar, MHImagePar)
2233
2234 - event display in executable changed to support also calibrated files
2235 (done with MJCalibrateSignal)
2236
2237 - added program doing the calibration and the signal extraction
2238 and signal calibration: callisto Most part of the setup is done
2239 by a resource file: callisto.rc
2240
2241 - implemented support for pixels with negative IDs (Raw file format
2242 version 4)
2243
2244 - implemented an new option in merpp (--auto-time) to get the time
2245 range for cc and/or caco data from the run header of the raw data
2246
2247 - New star-program (star and MJStar). Taking calibrated data from
2248 callisto/MJCalibrateSignal as input. The setup is done through
2249 a resource file: star.rc
2250
2251 - Added possibility to setup the eventloop in the mars-executable
2252 Camera Display by a resource file: mars.rc
2253
2254 - The the mars-executable Camera Display can now also read
2255 calibrated data files written by callisto/MJCalibrateSignal
2256
2257 - Added a new class which can calculate the source position in
2258 the camera from a pointing model known to be correct and
2259 the pointing model used for observations (MSrcPosFromModel)
2260
2261 - added a new version of the class calculating the effective on-time
2262 MHEffectiveOnTime. It unifies the old classes MHEffOnTime,
2263 MHEffOnTimeTime, MHEffOnTimeTheta, MHTimeDiff, MHTimeDiffTime and
2264 MHTimeDiffTheta.
2265
2266 - added the possibility to count how often a signal is above (or
2267 below) a certain threshold in MHCamEvent and MHCamEventRot (this
2268 can replace MHTriggerLvl0)
2269
2270 - added ore convinient ways to set the axis-titles of your histograms
2271 in your macro to some MH-classes like: MH3, MHVsTime, MHSectorVsTime,
2272 MHPixVsTime, MHVsTime, MHCamEvent and MHCamEventRot
2273
2274
2275
2276 *** Version 0.8.5 (2004/08/06)
2277
2278 - new Signal Extractors: MExtractFixedWindowSpline
2279 and MExtractAmplitudeSpline
2280
2281 - implementes interleave in merpp to be able to merpp parts of a file
2282
2283 - implemented force-functionality to merpp
2284
2285 - star display to display star catalogs (MAstroCatalog)
2286
2287 - star display to 'mirror' stars to the camera (MAstroCamera)
2288
2289 - implemented a new cleaning option which allows cleaning scaled with
2290 the real ratio of the rms depending on the pixel size, not only the
2291 square root of the pixelsize (MImgCleanStd, MCameraData)
2292
2293 - removed obsolete MBlind* classes
2294
2295 - Signal extractor taking the AB flag into account implemented (this makes
2296 an odd number of slices for signal extraction possible
2297 (MExtractSignalABcorr)
2298
2299 - implemented a class extracting the pedestal from the raw-data (events)
2300 itself using empty lo-gain (MpedCalcFromLoGain)
2301
2302 - Added mean values to MPedPhotCam
2303
2304 - Added Makefile for Darwin (this OS is NOT officially supported!)
2305
2306
2307
2308 *** Version 0.8.4 (2004/04/19)
2309
2310 - new fast arrival time extractor using cubic splines:
2311 MExtractTimeFastSpline
2312
2313 - implementes multi-argument support in MDataChain
2314 ("pow(a, b)")
2315
2316 - implemented double reference support in MDataChain
2317 ("MCameraLV.fPowerSupplyA.fVoltagePos5V")
2318
2319 - added support in MF for expressiond like
2320 "MHillas.fWidth<2*<MHillas.fLength"
2321
2322 - MDataChain is now able to support variables like [0], [1], ...
2323 which can be used in fit functions as parameters. The interface
2324 is implemented through the new virtual function
2325 MParContainer::SetVariables
2326
2327 - added new class MArrivalTimeCam/MArrivalTimePix:
2328 Storage containers to hold the ArrivalTime of the camera/pixel
2329
2330 - added new classes for the calibration of relative arrival times:
2331 MCalibrationRelTimePix and MCalibrationRelTimeCam
2332
2333 - added new classes to extract the signal of the blind pixel and the
2334 PIN Diode separately: MExtractBlindPixel and MExtractPINDiode
2335
2336 - added new classes for the bad-pixels treatment (MBadPixels*)
2337 which are more powerfull than the old ones (MBlindPixel*)
2338 and will replace them.
2339
2340 - added filter against cosmics: MFCosmics
2341
2342 - added new class MArrivalTimeCalc2:
2343 Calculates the arrival time as the signal weighted mean of the time
2344 slices in the maximum window
2345
2346 - added a new macro call pointspreadfucntion.C that fits with a 2D
2347 gaussian the DC spot for an star. It gives you the RMS of
2348 the PSF and the center of the star, very useful for misspointing
2349 studies.
2350
2351 - added support for multi-telescope camera files. Added classes
2352 MStereoCalc and MStereoPar, classes to calculate and store shower
2353 parameters calculated from the data of two telescopes, for the MC
2354 studies regarding the MAGIC clone.
2355
2356 - added MMcCalibrationCalc, MMcCalibrationUpdate, classes to obtain
2357 MC data calibrated in photons instead of ADC counts.
2358
2359 - added three macros to show the calibration and analysis chain of
2360 MC data. The macro starmc.C does both calibration and analysis, producing
2361 an output containing Hillas parameters for each event. The macro
2362 mccalibrate.C does just the calibration and writes out a file containing
2363 pixel information in photons instead of FADC data. The macro starmc2.C
2364 reads in such a file and produces an output with Hillas parameters.
2365
2366 - added MPedPhotCalc to compute pedestals in units of photons, see
2367 macro pedphotcalc.C
2368
2369 - added all class to perform the calibration, see macro calibration.C
2370
2371 - added class MFFT to perform Fast Fourier Transforms
2372
2373 - added class MSimulatedAnnealing to perform simulated annealing
2374 minimizations
2375
2376 - added new macro bootcampstandardanalysis.C which holds the skeleton
2377 of the analysis chain discussed in the december bootcamp.
2378
2379 - added support for DAQ raw data file format V2
2380
2381 - implemented an arbitrary time stamp for the events from raw-data
2382 files with file version<3. It is calculated from the run number,
2383 the date of start and the FADC clock tick from the first module
2384
2385 - implemented the correct decoding of the event time-stamp
2386 coming from the DAQ for raw-data files with file-version>2
2387
2388 - added a new classes storing the pedestals in units of photons:
2389 MPedPhotPix, MPedPhotCam
2390
2391 - added support for DC current files from the camera control
2392
2393 - added support for start- and stop-time when merpping report files
2394
2395 - added MReportCC storing the weather data from the CC-REPORT
2396
2397 - added class to fill trainings and test matrices (MTFillMatrix)
2398
2399 - added a filter performing a software trigger (MFSoftwareTrigger)
2400
2401 - Added classes to perform parts of the standard analysis (Jobs):
2402 + MJPedestals
2403 + MJCalibration
2404 + MJExtractSignal
2405 Note: All of them are preliminary!
2406
2407 - Added Iterator to iterate over files giving only the run-number
2408
2409 - Added class (MHSectorVsTime) and example (sectorvstime.C) to plot
2410 the evolution of a sector vs time
2411
2412
2413
2414 *** Version 0.8.3 (2003/12/05)
2415
2416 - adapted Random Forest classes to allow the use with MC
2417 multi-telescope files which can be produced with the coming
2418 camera version 0.7.
2419
2420 - added classes to read and interprete central control report files
2421 (Directory: mreport, Base-Class: MReport, Reader: MReportFileRead, ...)
2422
2423 - enhanced merpp (the executable) to support .raw and .rep files.
2424
2425 - added classes for light calibration (MCalibration* and
2426 MHCalibration*) by Markus Gaug. Further work has to be and will
2427 be done soon.
2428
2429 - Added a preliminary 'synchronizer' (MReadReports) which will
2430 replace MReadMarsFile as soon as it is finished.
2431
2432 - Added a tutorial macro demonstrating how report files are merpped
2433 and how they can be read from the root file in the correct order
2434 of time (readrep.C)
2435
2436 - Added a generalized histogram (similar to MH3) called MVsTime
2437 which allows to plot variables (or rules - MDataChain - of
2438 variables) vs time. For an example have a look at tempvstime.C.
2439
2440 - Changed the name of the event trees:
2441 + Events --> Events
2442 + PedEvents --> Pedestals
2443 + CalEvents --> Calibration
2444
2445 - changed name of MTime in the event trees from MRawEvtTime to MTime
2446
2447 - added displaying the sector indices to the mars event display
2448
2449 - Magic Online Analysis (MOnA) implemented in a first version
2450 (see mona.cc, MOnlineDump and MOnlineDisplay)
2451
2452 - added classes to calculate event rate
2453 (preliminary: MEventRateCalc, MEventRate)
2454
2455 - added new executable 'readdaq' to read *.raw files from the DAQ directly
2456
2457 - added classes dealing with the telescope pointing position (MPointing*)
2458
2459 - implemented the 'final' Mars time stamp based on some experiences with
2460 the new 'Synchonizer'
2461
2462
2463
2464 *** Version 0.8.2 (2003/11/10)
2465
2466 - accelerated standard analysis by a factor of 6-8 ...
2467 + doing precalculations in MGeomCam for the pixel area ratio
2468 + creating a look up table for the single pixels in MCerPhotEvt
2469
2470 - added task to calculates weights in case you want to 'change'
2471 the energy spectrum of a monte carlo file
2472 (MMcWeightEnergySpecCalc - Great work from Marcos Lopez Moya)
2473
2474 - added a new executable star which implements the
2475 star.C macro as precompiled executable. Currently a flexible
2476 setup is not possible, but a 'input card' setup will follow.
2477
2478 - Starting to implement Multi-Telescope-Support based on the
2479 upcomming camera version 0.7. The base for a general support
2480 is implemented and the star-functionality is proven to work.
2481
2482 - Replaced Event-Display in and Camera-Display in the Mars GUI
2483 executable by a more convinient display base on MStatusDisplay.
2484
2485 - fixed the algorithm for the blind pixel treatment (MBlindPixelCalc)
2486
2487
2488
2489 *** Version 0.8.1 (2003/10/20)
2490
2491 - added signal subtraction for pure on data by means of fitting
2492 the background in the off region or by performing a combined
2493 signal/background fit. Provides necessary histograms for
2494 obtaining energy spectra and a light curve.
2495
2496 - added classes to perform and study the selection of the
2497 2nd Level Trigger on MC data (example in triglvl2.C macro)
2498
2499 - added pedestal calculation
2500
2501 - implemented "democratic cleaning" (tail cut value proportional
2502 to pixel area) according to Wolfgang's suggestion.
2503
2504 - added Sigmabar treatment (calculation and basic histogramming)
2505
2506 - added basic routines to apply padding
2507
2508 - added first implementation of MARS-based source visibility plotter
2509 (needs to be linked against slalib, preliminary version)
2510
2511 - enhanced functionality of data chains: random numbers, conversion
2512 functions and access to MHMatrix objects
2513
2514 - fixed some bugs in the CT1 Preproc file reader, which gave wrong
2515 numbers under certain circumstances (runs with 0 events included)
2516
2517 - fixed a bug in the general filter MF (didn't work on some systems
2518 at all)
2519
2520 - fixed a bug in chaining files, which caused the analysis to crash or
2521 ignore files
2522
2523 - implemented more geometry classes describing the parameters used
2524 to simulate the PMT camera
2525
2526 - added new Monte Carlo classes containing the setup of the different
2527 programs in the MC chain.
2528
2529 - added a (preliminary) general task to evaluate chi-square of a given
2530 parameter in an eventloop.
2531
2532 - first implementation of a container describing the observatory
2533 location (Long, Lat, etc)
2534
2535 - added Random Forest - method for g/h-separation
2536
2537 - made compatible with the latest PRO version of root (3.04/02 and 3.05/07)
2538 (this means, that it is compiling, but not yet fully tested)
2539
2540 - added a new status display which can show the present status
2541 of histograms while an eventloop is running (MStatusDisplay,
2542 an example can be found at status.C)
2543
2544 - reorganized directories:
2545 classes describing the image --> mimage
2546 classes containing Monte Carlo Histograms --> mhistmc
2547 classes for Random Forest method --> mranforest
2548
2549 - added (event-)weights to fill histograms
2550
2551 - colors in logging output to screen
2552
2553 - changes to the Image Parameter declaration:
2554 + moved fConc, fConc1 from MHillasSrc to MNewImagePar
2555 + moved fNumUsedPix, fNumCorePix from MHillas to MNewImagePar
2556 + MHillasExt now derives from MParContainer instead of MHillas
2557 --> Files written with the old definition are more or less unreadable
2558
2559 - reading of reflector files
2560
2561 - display data from reflector files in the camera
2562
2563 - implemented currents (and display them in the camera)
2564
2565 - logarithmic scale in camera
2566
2567 - added class to iterate through directory contents (MDirIter)
2568
2569
2570
2571 *** Version 0.8 (2002/11/18)
2572
2573 - added all subdirectories as Include-Paths for ACLIC (s. rootlogon.C)
2574
2575 - corrected 48(!) wrong entries in the Next Neighbor table of the
2576 Magic Camera (all in the inner part of the outer ring, mostly
2577 missing a single NN)
2578
2579 - Fixed a bug in the table describing the Next Neighbours of the
2580 CT1 camera (a pixel 127 doesn't exist and the pixel 126 didn't
2581 have any neighbour)
2582
2583 - changed trigrate.C: now it calculates the pure NSB triggers from
2584 a camera file (generated with no showers), so there is no need
2585 to introduce them by hand in the BgR array (for more details
2586 please ask Abelardo)
2587
2588 - added MHHadronness (and MHadronness) as a general facility to evaluate
2589 the quality of a gamma-hadron seperation method
2590
2591 - changed MWriteAsciiFile to accept rules (MDataChain)
2592
2593 - Calls to MWriteAsciiFile::AddContainer and MWriteAsciiFile::AddRule
2594 must be replaced by calles to MWriteAsciiFile::AddColumn and
2595 MWriteAsciiFile::AddColumns. Be carefull, the arguments have slightly
2596 changed.
2597
2598 - Implemented generalized event-matrices (one row per event)
2599 (MHMatrix)
2600
2601 - implemented the one dimensional composite probabilities
2602 (MHCompProb, MCompProbCalc, comprob.C)
2603
2604 - implemented the Multidimensional Distances (Next Neighbours, Kernel)
2605 (MMultiDimDistCalc, multidimdist.C, multidimdist2.C, MHMatrix)
2606
2607 - Added a camera display displaying the pixel numbers to the
2608 event display gui.
2609
2610 - Added Monte Carlo Informations to event display gui
2611
2612 - Changed the camera display to display the pixel numbering
2613
2614 - Added three buttons to the camera display to change the palette
2615
2616 - Added the number of used and core pixels after image cleaning
2617 to MHillas
2618
2619 - Changed the algorithm used in MHillas according to TDAS 02-03
2620
2621 - Changed the Alpha range according to TDAS 02-03 to [-90°,90°]
2622
2623 - Changed the algorithm used in MHillasSrc according to TDAS 02-03
2624
2625 - Added fCosDeltaAlpha to MHillasSrc
2626
2627 - The numbers of photons used to calculate fConc and fConc1 in
2628 MHillasExt are now scaled with the pixel size, so that one get
2629 a four times smaller value for the bigger pixels in the outer ring.
2630
2631 - added new task to smooth the camera contents (MCameraSmooth)
2632
2633 - added possibility to use interpolated pixel values for blind pixels
2634 instead of removing it completely from the analysis
2635
2636 - Implemented the possibility to save the whole eventloop with
2637 its setup to a file. This information can be read and from
2638 it you can (re)create a Macro. See MEvtLoop::Read, MEvtLoop::Write
2639 and MEvtLoop::MakeMacro. If you find something not working, please
2640 report - this stuff is still in a beta phase.
2641
2642 - MBlindPixelCalc set all Ceta-Tauri pixels to unused not only for the
2643 file having crab starfield, but also for all following files - fixed
2644
2645 - MTaskList::PrintStatistics can now be instructud to print also the
2646 title, too: use PrintStatistics(0, kTRUE)
2647
2648 - Changed the image cleaning so that pixels with to many 'used'
2649 neighbors are left used (to get rid of 'holes' in events)
2650
2651 - Introduced more output to the camera displays
2652
2653 - Added an array of histograms (eg one hist per run), MHArray
2654
2655 - Added a _preliminary_ version of the reader for CT1 PreProc files
2656 (MCT1ReadPreProc)
2657
2658 - Fixed the error calculation of MHMCollectionArea and MHMcIntRate
2659 (thanks to Raquel)
2660
2661 - added the first energy estimator (MEnergyEstParam, estimate.C,
2662 estfit.C) using a parametrization
2663
2664 - Added some new Monte Carlo paremeter containers used to store
2665 setup information from the MC chain (MMcConfigRunHeader,
2666 MMcCorsikaRunHeader, MGeomPMT, MGeomMirror)
2667
2668
2669
2670 *** Version 0.7 (2002/04/30)
2671
2672 - added a bugfix to MCerPhotCalc. In older camera versions (<=40)
2673 the pedestal mean value was saved incorrectly. For files from
2674 this versions we substract 0.5 from the pedestal mean.
2675 WARNING: This may effect your results, so don't wonder...
2676
2677 - First implementation of an algorithm using only triggered events
2678 for the collection area calculation
2679
2680 - Ascii Output can now also be used for parameter containers which
2681 doesn't overload MParCointainer::AsciiWrite
2682
2683 - The Ascii Output is now also capable of writing single data members
2684 of one container
2685
2686 - You are now able to change the order of the values written to the
2687 ascii file
2688
2689 - You can now specify a conversion factor for each data member written
2690 to an ascii file. This may be usefull to change the units of the
2691 data member (eg. degrees instead of millimeters in case of the
2692 hillas parameters)
2693
2694 - Replaced old MHillas by a new structure which allows you to extend
2695 the parameters stored in MHillas very easily:
2696 see MHillas and MHillasExt
2697
2698 - Added classes to handle source dependancy of image parameters:
2699 see MHillasSrc
2700
2701 - Added container (MBinning) to have a standard input for the binning
2702 in different histograms (eg. the Energy bins should be the same in
2703 all histograms)
2704
2705 - Changed Hillas histograms from mm to deg
2706
2707 - Added the flexible binning to the hillas histograms
2708
2709 - Added a filter for the alpha parameter (MFAlpha)
2710
2711 - It is now possible to write single data members of a class object to
2712 an output stream instead of the whole container only
2713
2714 - Added a generalized filter for a data member:
2715 MFDataMember filter("MHillas.fWidth", '<', 0.5);
2716
2717 - Added a generalized Filter parser for more difficult filter rules
2718 MF filter("MHillas.fWidth<0.5 && MHillas.fLength>0.5");
2719
2720 - Added new Monte Carlo histograms:
2721 MHMcEnergyImpact and MHMcEfficiency
2722
2723 - Added many new histograms needed for the flux calculation.
2724
2725 - Added a generalized histograms which can fill and show up to three
2726 data members in either a one-, two- or three-dimensional histogram, eg:
2727 MH3 mh3("MHillas.fWidth", "MHillas.fLength");
2728
2729 - Added:
2730 * calculation of <Theta> as a function of Theta (MHThetabarTheta)
2731 * calculation of <Theta> as a function of time (MHThetabarTime)
2732 they are needed in the flux calculation in
2733 order to select for a given bin in Theta (or time) the appropriate
2734 effective collection area, which is assumed to be given as a function
2735 of Theta and the true energy
2736
2737 - Added calculation of the migration matrix in the energy
2738 (MHMcEnergyMigration):
2739 the migration matrix describes the migration from the true (E_true)
2740 to the estimated energy (E_est); it is needed to determine the
2741 distribution of E_true from the measured distribution of E_est
2742 by unfolding
2743
2744 - changed the color palette in the camera display to DeapSea.
2745
2746
2747
2748 *** Version 0.6 (2002/01/15)
2749
2750 - Introduce the option of disable pixels (MBlindPixels and
2751 MBlindPixelCalc) for the analysis (eg MHillasCalc) to be
2752 able to deal with Star Field NSB.
2753
2754 - Fixed a bug in the pedestals substraction
2755
2756 - Introduced verbosity levels
2757
2758 - Introduced a task (MReadMarsFile) which reads the RunHeader tree, too.
2759
2760 - Fixed several bugs in MFTriggerLvl1 which caused to filter not to
2761 work correctly (thanks to Abelardo)
2762
2763 - Introduced correct treatment of the diffuse Night Sky Background
2764 and NSB (MMcPedestalNSB)
2765
2766 - Corrected treatment of the electronic noise (MMcPedestalCopy)
2767
2768 - Introduced MMcRunHeaders which contains monte carlo run informations
2769
2770 - Added a 'Print'-Task (MPrint) which calls the Print function
2771 of a parameter container for each event
2772
2773 - Removed MHtml (was replaced by THtml of root 3.02) in dohtml.C
2774
2775 - changes to be more efficient with root 3.02/06
2776
2777 - root 3.02/06 fixes:
2778 + the memory leak when reading MRawEvtData
2779 + loosing more time than necessary to update the progress bar
2780
2781 - MWriteAsciiFile is now capable of writing more than one container
2782 in a line. For example: To be able to write Monte Carlo and
2783 Hillas values (energy and alpha) in one line.
2784
2785
2786
2787 *** Version 0.5 (2001/11/06)
2788
2789 - Added new GUI functionality (Analysis)
2790
2791 - Added a camera display to be able to display events in the camera
2792 together with the calculated hillas ellipse
2793
2794 - Added a 'Clone'-Task (MClone) which can clone a parameter container
2795 object at any state of the analysis, so that the output are several
2796 states
2797
2798 - Fixed a bug in the collection area error calculation
2799 (Thanks to Ciro and Abelardo)
2800
2801 - Fixed a bug which causes merpp to crash in some environments
2802
2803 - Implemented auto enabling scheme for Branches to speed up reading
2804 a lot (see MReadTree and MTask)
2805
2806 - Fixed a bug in the Histogram classes which caused histograms to
2807 be written to a file by MWriteRootFile without request
2808
2809 - changed the object I/O from the old stylish root I/O to the new
2810 'auto evolution scheme' This was done only for classes which are not
2811 yet used in Monte Carlo files, so that these files are still readable.
2812
2813 - Added support for a progress bar in MEvtLoop and MReadTree
2814 (see the Analysis Window in Mars for an example)
2815
2816 - Fixed a bug in the package creation which disturbed the conversion
2817 between a mars version unpacked from a tar file and the repository.
2818
2819 - Changed the code to writer merpped (converted into root format)
2820 raw files, to get faster (around 10%)
2821
2822
2823
2824 *** Version 0.4 (2001/09/20)
2825
2826 - This release is made to compile on a root version >= 3.00
2827
2828 - Fixed a bug which causes MReadTree to read the first file only
2829
2830 - Added a task to compute the trigger rate (MMcTriggerRateCalc)
2831
2832 - Added a task to write a container to an Ascii file (MWriteAsciiFile)
2833
2834 - Added a task to write several container to a root file (MWriteRootFile)
2835
2836 - Added calculation of the Enegry Threshold (MMcThresholdCalc)
2837
2838 - Added calculation of the collection area (MMcCollectionAreaCalc)
2839
2840 - fixed some bugs in the Hillas calculation
2841
2842 - added filters to be able to control the task execution dependent on
2843 a parameter (for example: the number of level 1 triggers in a MC-file)
2844
2845
2846
2847 *** Version 0.3 (2001/04/23)
2848
2849 - Many demo macros added
2850
2851 - MCT1ReadAscii is now able to handle more than one file (s. ::AddFile)
2852
2853 - Hillas Calculation added
2854
2855 - Standard image cleaning available
2856
2857 - MReadTree is now able to handle more than one file
2858 (Remark: of the same structure)
2859
2860 - Now we are able to display events in a first version of the
2861 Camera event display. To do this we must run the macro
2862 "macros/readCT1.C" to read data from CT1.
2863
2864
2865
2866 *** Version 0.2 (2001/01/24)
2867
2868 - First Implementation of a Gui to read the raw data. Start the
2869 gui with the command mars.
2870
2871
2872
2873 *** Version 0.1 (2000/12/20)
2874
2875 - First official release 0.1
2876
2877 - Changes:
2878
2879 + changed calling style of MRawPixelIterator from do...while to
2880 while-loop
Note: See TracBrowser for help on using the repository browser.