source: trunk/MagicSoft/Mars/NEWS@ 8976

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