source: trunk/Mars/NEWS@ 12634

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