source: trunk/Mars/NEWS@ 9878

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