source: trunk/MagicSoft/Mars/NEWS@ 9620

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