source: trunk/MagicSoft/Mars/NEWS@ 8889

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