source: trunk/MagicSoft/Mars/NEWS@ 8938

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