source: trunk/MagicSoft/Mars/NEWS@ 9098

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