source: trunk/MagicSoft/Mars/NEWS@ 9160

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