source: trunk/MagicSoft/Mars/NEWS@ 8955

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