source: trunk/MagicSoft/Mars/NEWS@ 9038

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