source: trunk/Mars/NEWS@ 9876

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