source: trunk/MagicSoft/Mars/NEWS@ 8993

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