source: trunk/MagicSoft/Mars/NEWS@ 8939

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