source: trunk/Mars/NEWS@ 9863

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