source: trunk/MagicSoft/Mars/NEWS@ 8907

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