source: trunk/MagicSoft/Mars/NEWS@ 8526

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