source: trunk/MagicSoft/Mars/NEWS@ 9579

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