source: trunk/MagicSoft/Mars/NEWS@ 9229

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