source: trunk/MagicSoft/Mars/NEWS@ 9213

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