source: trunk/MagicSoft/Mars/NEWS@ 9068

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