source: trunk/MagicSoft/Mars/NEWS@ 9467

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