source: trunk/MagicSoft/Mars/NEWS@ 9564

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