source: trunk/MagicSoft/Mars/NEWS@ 9500

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