source: tags/Mars-V1.2/NEWS@ 13369

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