source: trunk/MagicSoft/Mars/NEWS@ 9069

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