source: trunk/Mars/NEWS@ 9947

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