source: trunk/MagicSoft/Mars/NEWS@ 9141

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