source: trunk/MagicSoft/Mars/NEWS@ 8553

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