source: trunk/MagicSoft/Mars/NEWS@ 9075

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