source: trunk/MagicSoft/Mars/NEWS@ 9138

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