source: trunk/MagicSoft/Mars/NEWS@ 8898

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