source: trunk/Mars/NEWS@ 9938

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