source: trunk/MagicSoft/Mars/NEWS@ 9181

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