source: trunk/MagicSoft/Mars/NEWS@ 8913

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