source: trunk/MagicSoft/Mars/NEWS@ 9061

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