source: trunk/MagicSoft/Mars/NEWS@ 8891

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