source: trunk/MagicSoft/Mars/NEWS@ 8970

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