source: trunk/MagicSoft/Mars/mjobs/MJCalibrateSignal.cc@ 8361

Last change on this file since 8361 was 8360, checked in by tbretz, 20 years ago
*** empty log message ***
File size: 28.1 KB
Line 
1/* ======================================================================== *\
2!
3! *
4! * This file is part of MARS, the MAGIC Analysis and Reconstruction
5! * Software. It is distributed to you in the hope that it can be a useful
6! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
7! * It is distributed WITHOUT ANY WARRANTY.
8! *
9! * Permission to use, copy, modify and distribute this software and its
10! * documentation for any purpose is hereby granted without fee,
11! * provided that the above copyright notice appear in all copies and
12! * that both that copyright notice and this permission notice appear
13! * in supporting documentation. It is provided "as is" without express
14! * or implied warranty.
15! *
16!
17!
18! Author(s): Thomas Bretz, 1/2004 <mailto:tbretz@astro.uni-wuerzburg.de>
19!
20! Copyright: MAGIC Software Development, 2000-2007
21!
22!
23\* ======================================================================== */
24
25/////////////////////////////////////////////////////////////////////////////
26//
27// MJCalibrateSignal
28//
29// This class is reading the output written by callisto. It calibrates
30// signal and time.
31//
32// The signal and time extractors are read from the callisto-output. In
33// pricipal you could overwrite these default using the resource file,
34// but this is NOT recommended!
35//
36/////////////////////////////////////////////////////////////////////////////
37#include "MJCalibrateSignal.h"
38
39#include <TEnv.h>
40#include <TFile.h>
41
42#include "MLog.h"
43#include "MLogManip.h"
44
45#include "MDirIter.h"
46#include "MRunIter.h"
47#include "MParList.h"
48#include "MTaskList.h"
49#include "MEvtLoop.h"
50
51#include "MStatusDisplay.h"
52
53#include "MGeomCam.h"
54#include "MHCamEvent.h"
55#include "MPedestalCam.h"
56#include "MBadPixelsCam.h"
57#include "MArrivalTimeCam.h"
58
59#include "MCalibrationQECam.h"
60#include "MCalibrationBlindCam.h"
61#include "MCalibrationChargeCam.h"
62#include "MCalibrationRelTimeCam.h"
63#include "MCalibrationChargePINDiode.h"
64#include "MCalibrationPulseTimeCam.h"
65#include "MCalibrationPatternDecode.h"
66
67#include "MCalibrationChargeCalc.h"
68#include "MCalibrationRelTimeCalc.h"
69
70#include "MCalibrationIntensityChargeCam.h"
71#include "MCalibrationIntensityBlindCam.h"
72#include "MCalibrationIntensityRelTimeCam.h"
73#include "MCalibrationIntensityQECam.h"
74#include "MCalibrationIntensityConstCam.h"
75#include "MBadPixelsIntensityCam.h"
76
77#include "MHCalibrationChargeCam.h"
78#include "MHCalibrationChargeBlindCam.h"
79#include "MHCalibrationChargePINDiode.h"
80#include "MHCalibrationRelTimeCam.h"
81#include "MHCalibrationPulseTimeCam.h"
82#include "MHCamera.h"
83
84#include "MCalibCalcFromPast.h"
85
86#include "MReadReports.h"
87#include "MReadMarsFile.h"
88#include "MRawFileRead.h"
89#include "MContinue.h"
90#include "MTriggerPatternDecode.h"
91#include "MFTriggerPattern.h"
92#include "MFilterList.h"
93#include "MGeomApply.h"
94#include "MPedestalSubtract.h"
95//#include "MMcPedestalCopy.h"
96#include "MPointingPosCalc.h"
97#include "MPedCalcFromLoGain.h"
98#include "MPedestalSubtract.h"
99#include "MExtractor.h"
100#include "MExtractTimeAndCharge.h"
101#include "MExtractPINDiode.h"
102#include "MExtractBlindPixel.h"
103#include "MFCosmics.h"
104#include "MTaskEnv.h"
105#include "MCalibrateData.h"
106#include "MCalibrateRelTimes.h"
107#include "MBadPixelsMerge.h"
108#include "MBadPixelsCalc.h"
109#include "MBadPixelsTreat.h"
110#include "MFillH.h"
111#include "MWriteRootFile.h"
112
113ClassImp(MJCalibrateSignal);
114
115using namespace std;
116
117// --------------------------------------------------------------------------
118//
119// Default constructor.
120//
121// Sets:
122// - fIsInterlaced to kTRUE
123// - fIsRelTimesUpdate to kFALSE
124// - fIsHiLoCalibration to kFALSE
125// - fPulsePosCheck to kTRUE
126//
127MJCalibrateSignal::MJCalibrateSignal(const char *name, const char *title)
128 : fExtractor(0), fIsInterlaced(kTRUE), fIsRelTimesUpdate(kFALSE)
129{
130 fName = name ? name : "MJCalibrateSignal";
131 fTitle = title ? title : "Tool to calibrate data";
132
133 SetPulsePosCheck();
134 //fCruns = NULL;
135}
136
137MJCalibrateSignal::~MJCalibrateSignal()
138{
139 if (fExtractor)
140 delete fExtractor;
141}
142
143void MJCalibrateSignal::SetExtractor(const MExtractor *ext)
144{
145 if (fExtractor)
146 delete fExtractor;
147
148 fExtractor = ext ? (MExtractor*)ext->Clone() : NULL;
149}
150
151Bool_t MJCalibrateSignal::WriteResult(TObjArray &cont)
152{
153 if (IsNoStorage())
154 return kTRUE;
155
156 const TString name(Form("signal%08d.root", fSequence.GetSequence()));
157 return WriteContainer(cont, name);
158}
159
160Bool_t MJCalibrateSignal::ReadCalibration(TObjArray &l, MBadPixelsCam &cam, MExtractor* &ext2, MExtractor* &ext3, TString &geom) const
161{
162 TString fname = Form("%s/calib%08d.root", fPathIn.Data(), fSequence.GetSequence());
163
164 *fLog << inf << "Reading from file: " << fname << endl;
165
166 TFile file(fname, "READ");
167 if (!file.IsOpen())
168 {
169 *fLog << err << dbginf << "ERROR - Could not open file " << fname << endl;
170 return kFALSE;
171 }
172
173 TObject *o = file.Get("ExtractSignal");
174 if (o && !o->InheritsFrom(MExtractor::Class()))
175 {
176 *fLog << err << dbginf << "ERROR - ExtractSignal read from " << fname << " doesn't inherit from MExtractor!" << endl;
177 return kFALSE;
178 }
179 ext3 = o ? (MExtractor*)o->Clone() : NULL;
180
181 o = file.Get("ExtractTime");
182 if (o && !o->InheritsFrom(MExtractor::Class()))
183 {
184 *fLog << err << dbginf << "ERROR - ExtractTime read from " << fname << " doesn't inherit from MExtractor!" << endl;
185 return kFALSE;
186 }
187 ext2 = o ? (MExtractor*)o->Clone() : NULL;
188 if (!ext3 && !ext2)
189 {
190 *fLog << err << dbginf << "ERROR - Neither ExtractSignal nor ExrtractTime found in " << fname << "!" << endl;
191 return kFALSE;
192 }
193
194 o = file.Get("MGeomCam");
195 if (o && !o->InheritsFrom(MGeomCam::Class()))
196 {
197 *fLog << err << dbginf << "ERROR - MGeomCam read from " << fname << " doesn't inherit from MGeomCam!" << endl;
198 return kFALSE;
199 }
200 geom = o ? o->ClassName() : "";
201
202 TObjArray cont(l);
203 cont.Add(&cam);
204 return ReadContainer(cont);
205}
206
207// --------------------------------------------------------------------------
208//
209// MJCalibration allows to setup several option by a resource file:
210// MJCalibrateSignal.RawData: yes,no
211//
212// For more details see the class description and the corresponding Getters
213//
214Bool_t MJCalibrateSignal::CheckEnvLocal()
215{
216 SetInterlaced(GetEnv("Interlaced", fIsInterlaced));
217 SetRelTimesUpdate(GetEnv("RelTimesUpdate", fIsRelTimesUpdate));
218
219 return MJCalib::CheckEnvLocal();
220}
221
222Bool_t MJCalibrateSignal::Process(MPedestalCam &pedcamab, MPedestalCam &pedcambias,
223 MPedestalCam &pedcamextr)
224{
225 if (!fSequence.IsValid())
226 {
227 *fLog << err << "ERROR - Sequence invalid..." << endl;
228 return kFALSE;
229 }
230
231 *fLog << inf;
232 fLog->Separator(GetDescriptor());
233 *fLog << "Calculate calibrated data from Sequence #";
234 *fLog << fSequence.GetSequence() << endl << endl;
235
236
237 //if (!CheckEnv())
238 // return kFALSE;
239
240 CheckEnv();
241
242 // --------------------------------------------------------------------------------
243
244 MDirIter iter;
245 if (fSequence.IsValid())
246 {
247 if (fSequence.SetupDatRuns(iter, 0, IsUseRawData())<=0)
248 return kFALSE;
249 }
250
251 // Read File
252 MCalibrationIntensityChargeCam ichcam;
253 MCalibrationIntensityQECam iqecam;
254 MCalibrationIntensityBlindCam iblcam;
255 MCalibrationIntensityRelTimeCam itmcam;
256 MCalibrationIntensityConstCam icncam;
257 MBadPixelsIntensityCam ibdcam;
258
259 MHCalibrationChargeCam hchacam;
260 MHCalibrationChargeBlindCam hbndcam;
261 MHCalibrationChargePINDiode hpndiod;
262 MHCalibrationRelTimeCam hrelcam;
263 //MHCalibrationHiLoCam hilocam;
264 //MHCalibrationPulseTimeCam hpulcam;
265
266 hchacam.SetOscillations(kFALSE);
267 hbndcam.SetOscillations(kFALSE);
268 hrelcam.SetOscillations(kFALSE);
269
270 MCalibrationChargeCam calcam;
271 MCalibrationQECam qecam;
272 MCalibrationBlindCam bndcam;
273 MCalibrationChargePINDiode pind;
274 MCalibrationRelTimeCam tmcam;
275 //MCalibrationHiLoCam hilcam;
276 MCalibrationPulseTimeCam pulcam;
277
278 MBadPixelsCam badpix;
279
280 TObjArray interlacedcont;
281 if (fIsInterlaced)
282 {
283 interlacedcont.Add(&ichcam);
284 interlacedcont.Add(&iqecam);
285 interlacedcont.Add(&iblcam);
286 interlacedcont.Add(&itmcam);
287 interlacedcont.Add(&ibdcam);
288 interlacedcont.Add(&icncam);
289 interlacedcont.Add(&hchacam);
290 if (IsUseBlindPixel())
291 interlacedcont.Add(&hbndcam);
292 if (IsUsePINDiode())
293 interlacedcont.Add(&hpndiod);
294 if (fIsRelTimesUpdate)
295 interlacedcont.Add(&hrelcam);
296 }
297
298 MExtractor *extractor1=fExtractor;
299 MExtractor *extractor2=0;
300 MExtractor *extractor3=0;
301 TString geom;
302
303 TObjArray calibcont;
304 calibcont.Add(&calcam);
305 calibcont.Add(&qecam);
306 calibcont.Add(&bndcam);
307 calibcont.Add(&tmcam);
308 if (IsUseBlindPixel())
309 calibcont.Add(&bndcam);
310
311 if (!ReadCalibration(calibcont, badpix, extractor2, extractor3, geom))
312 return kFALSE;
313
314 *fLog << all;
315 if (!geom.IsNull())
316 *fLog << inf << "Camera geometry found in file: " << geom << endl;
317 else
318 *fLog << inf << "No Camera geometry found using default <MGeomCamMagic>" << endl;
319
320 if (extractor3)
321 {
322 *fLog << underline << "Signal Extractor found in calibration file" << endl;
323 extractor3->Print();
324 *fLog << endl;
325 }
326 else
327 *fLog << inf << "No Signal Extractor: ExtractSignal in file." << endl;
328
329
330 if (extractor1)
331 {
332 *fLog << underline << "Modified Signal Extractor set by user." << endl;
333 extractor1->Print();
334 *fLog << endl;
335 }
336 else
337 {
338 *fLog << inf << "No modified Signal Extractor set by user... using ExtractSignal." << endl;
339 extractor1 = extractor3 ? (MExtractor*)extractor3->Clone() : 0;
340 }
341
342 if (extractor2)
343 {
344 *fLog << underline << "Time Extractor found in calibration file" << endl;
345 extractor2->Print();
346 *fLog << endl;
347 }
348 else
349 *fLog << inf << "No Time Extractor: ExtractTime in file." << endl;
350
351 // This is necessary for the case in which it is not in the files
352 MBadPixelsCam badcam;
353 icncam.SetBadPixels(&badpix);
354
355 // Setup Parlist
356 MParList plist;
357 plist.AddToList(this); // take care of fDisplay!
358 plist.AddToList(&badcam);
359 //plist.AddToList(&hilcam);
360 plist.AddToList(&calibcont); // Using AddToList(TObjArray *)
361 plist.AddToList(&interlacedcont); // Using AddToList(TObjArray *)
362 plist.AddToList(&pulcam);
363
364 // Setup Tasklist
365 MTaskList tlist;
366 plist.AddToList(&tlist);
367
368 MReadReports readreal;
369 readreal.AddTree("Events", "MTime.", MReadReports::kMaster);
370 readreal.AddTree("Trigger");
371 readreal.AddTree("Camera");
372 readreal.AddTree("Drive");
373 readreal.AddTree("CC");
374 readreal.AddTree("Currents");
375
376 MReadMarsFile readmc("Events");
377 readmc.DisableAutoScheme();
378
379 MRawFileRead rawread(NULL);
380
381 MRead *read = 0;
382 switch (GetDataFlag())
383 {
384 case kIsUseRawData: read = &rawread; break;
385 case kIsUseMC: read = &readmc; break;
386 case kIsUseRootData: read = &readreal; break;
387 }
388 read->AddFiles(iter);
389
390 const TString fname(Form("%s{s/_D_/_Y_}{s/\\.raw$/.root}{s/\\.raw\\.gz$/.root}", fPathOut.Data()));
391
392
393 // Skips MC which have no contents. This are precisely the
394 // events which fullfilled the MC Lvl1 trigger and an
395 // arbitrary cut (typically at 50phe) to speed up simulation
396 MContinue contmc("MRawEvtData.GetNumPixels<0.5", "ContEmptyMC");
397
398 //MPointingPosInterpolate pextr;
399 //pextr.AddFiles(&iter);
400
401 MGeomApply apply; // Only necessary to create geometry
402 if (!geom.IsNull())
403 apply.SetGeometry(geom);
404 MBadPixelsMerge merge(&badpix);
405
406 // Make sure that pedcamab has the correct name
407 pedcamab.SetName("MPedestalFundamental");
408 pedcamextr.SetName("MPedestalFromExtractorRndm");
409 pedcambias.SetName("MPedestalFromExtractor");
410 plist.AddToList(&pedcamextr);
411 plist.AddToList(&pedcambias);
412 plist.AddToList(&pedcamab);
413
414 MArrivalTimeCam timecam;
415 plist.AddToList(&timecam);
416
417 // Check for interleaved events
418 MCalibrationPatternDecode caldec;
419 MTriggerPatternDecode decode;
420
421 // This will make that for data with version less than 5, where
422 // trigger patterns were not yet correct, all the events in the real
423 // data file will be processed. In any case there are no interleaved
424 // calibration events in such data, so this is fine.
425 MFTriggerPattern ftp;
426 ftp.SetDefault(kTRUE);
427 // ftp.RequireCalibration();
428 ftp.DenyCalibration();
429 ftp.DenyPedestal();
430 // ftp.DenyPinDiode();
431 ftp.SetInverted();
432
433 // This will skip interleaved calibration events and pedestal events (if any)
434 // --> tlist2
435 MContinue conttp(&ftp, "ContTrigPattern");
436
437 // Create the pedestal subtracted raw-data
438 MPedestalSubtract pedsub;
439 pedsub.SetPedestalCam(&pedcamab);
440
441 // Do signal and pedestal calculation
442 MPedCalcFromLoGain pedlo1("MPedCalcFundamental");
443 pedlo1.SetPedestalUpdate(kTRUE);
444 pedlo1.SetNamePedestalCamOut("MPedestalFundamental");
445
446 MPedCalcFromLoGain pedlo2("MPedCalcWithExtractorRndm");
447 pedlo2.SetPedestalUpdate(kTRUE);
448 pedlo2.SetRandomCalculation(kTRUE);
449 pedlo2.SetNamePedestalCamIn("MPedestalFundamental");
450 pedlo2.SetNamePedestalCamOut("MPedestalFromExtractorRndm");
451
452 MPedCalcFromLoGain pedlo3("MPedCalcWithExtractor");
453 pedlo3.SetPedestalUpdate(kTRUE);
454 pedlo3.SetRandomCalculation(kFALSE);
455 pedlo3.SetNamePedestalCamIn("MPedestalFundamental");
456 pedlo3.SetNamePedestalCamOut("MPedestalFromExtractor");
457
458 if (!extractor1)
459 {
460 *fLog << err << "ERROR - extractor1 == NULL" << endl;
461 return kFALSE;
462 }
463
464 extractor1->SetPedestals(&pedcamab);
465
466 // Setup to use the hi-gain extraction window in the lo-gain
467 // range (the start of the lo-gain range is added automatically
468 // by MPedCalcFromLoGain)
469 //
470 // The window size of the extractor is not yet initialized,
471 // so we have to stick to the extraction range
472 //
473 // Even if we would like to use a range comparable to the
474 // hi-gain extraction we use the lo-gain range to make
475 // sure that exclusions (eg. due to switching noise)
476 // are correctly handled.
477 //
478 pedlo1.SetRangeFromExtractor(*extractor1);
479
480 if (extractor1->InheritsFrom("MExtractTimeAndCharge"))
481 {
482 pedlo2.SetExtractor((MExtractTimeAndCharge*)extractor1);
483 pedlo3.SetExtractor((MExtractTimeAndCharge*)extractor1);
484 }
485 else
486 {
487 pedlo2.SetRangeFromExtractor(*extractor1);
488 pedlo3.SetRangeFromExtractor(*extractor1);
489 }
490
491 if (extractor2)
492 extractor2->SetPedestals(&pedcamab);
493
494 if (extractor3)
495 extractor3->SetPedestals(&pedcamab);
496
497 //------------------------------
498 MFTriggerPattern ftp2;
499 ftp2.SetDefault(kTRUE);
500 ftp2.DenyCalibration();
501 if (!extractor1->HasLoGain())
502 ftp2.RequirePedestal();
503
504 pedlo1.SetFilter(&ftp2);
505 pedlo2.SetFilter(&ftp2);
506 pedlo3.SetFilter(&ftp2);
507
508 MContinue contftp2(&ftp2, "ContPedestal");
509 //------------------------------
510
511 MFCosmics fcosmics;
512 fcosmics.SetNamePedestalCam("MPedestalFundamental");
513 MContinue contcos(&fcosmics, "ContTrigEvts");
514 contcos.SetInverted();
515
516 //MMcPedestalCopy pcopy;
517 MTaskEnv taskenv1("ExtractSignal");
518 MTaskEnv taskenv2("ExtractTime");
519 MTaskEnv taskenv3("ExtractInterlaced");
520 taskenv1.SetDefault(extractor1);
521 taskenv2.SetDefault(extractor2);
522 taskenv3.SetDefault(extractor3);
523
524 MFilterList flistftp2("PedestalFilter");
525 flistftp2.SetInverted();
526 flistftp2.AddToList(&ftp2);
527
528 if (!extractor1->HasLoGain())
529 {
530 taskenv1.SetFilter(&flistftp2);
531 taskenv2.SetFilter(&flistftp2);
532 taskenv3.SetFilter(&flistftp2);
533 }
534
535 //
536 // This is new calibration to photo-electrons, hard-coded
537 // as decided at the Wuerzburg software meeting 26.01.05
538 //
539 MCalibrateData calib;
540 calib.SetSignalType(MCalibrateData::kPhe);
541 //calib.AddPedestal("Fundamental");
542 calib.AddPedestal("FromExtractor");
543 calib.AddPedestal("FromExtractorRndm");
544 calib.SetPedestalFlag(MCalibrateData::kEvent);
545 //----------------------------------------------------------
546
547 MExtractPINDiode pinext;
548 MExtractBlindPixel bldext;
549
550 MFTriggerPattern fcalib("CalibFilter");
551 fcalib.SetDefault(kFALSE);
552 fcalib.RequireCalibration();
553
554 MCalibrationChargeCalc chcalc;
555 chcalc.SetExtractor(extractor3);
556
557 MCalibrationRelTimeCalc recalc;
558 MCalibCalcFromPast pacalc;
559
560 pinext.SetPedestals(&pedcamab);
561 bldext.SetPedestals(&pedcamab);
562 chcalc.SetPedestals(&pedcamextr);
563
564 pacalc.SetChargeCalc(&chcalc);
565 if (fIsRelTimesUpdate)
566 pacalc.SetRelTimeCalc(&recalc);
567 pacalc.SetCalibrate(&calib);
568
569 //
570 // Intensity Calibration histogramming
571 //
572 MFillH filpin(&hpndiod, "MExtractedSignalPINDiode", "FillPINDiode");
573 MFillH filbnd(&hbndcam, "MExtractedSignalBlindPixel", "FillBlindCam");
574 MFillH filcam(&hchacam, "MExtractedSignalCam", "FillChargeCam");
575 MFillH filtme(&hrelcam, "MArrivalTimeCam", "FillRelTime");
576 //MFillH filhil(&hilocam, "MExtractedSignalCam", "FillHiLoRatio");
577 //MFillH filpul(&hpulcam, "MRawEvtData", "FillPulseTime");
578 filpin.SetBit(MFillH::kDoNotDisplay);
579 filbnd.SetBit(MFillH::kDoNotDisplay);
580 filcam.SetBit(MFillH::kDoNotDisplay);
581 filtme.SetBit(MFillH::kDoNotDisplay);
582 //filhil.SetBit(MFillH::kDoNotDisplay);
583 //filpul.SetBit(MFillH::kDoNotDisplay);
584
585 MCalibrateRelTimes caltm;
586 MBadPixelsCalc bpcal;
587 MBadPixelsTreat treat;
588
589 //bpcal.SetNamePedPhotCam("MPedPhotFromExtractor");
590 bpcal.SetNamePedPhotCam("MPedPhotFromExtractorRndm");
591
592 //treat.AddNamePedPhotCam("MPedPhotFundamental");
593 treat.AddNamePedPhotCam("MPedPhotFromExtractor");
594 treat.AddNamePedPhotCam("MPedPhotFromExtractorRndm");
595 if (!extractor2 && !extractor1->InheritsFrom("MExtractTimeAndCharge"))
596 treat.SetProcessTimes(kFALSE);
597
598 MHCamEvent evt0( 0, "PedFLG", "Fundamental Pedestal from Lo Gain;;P [cnts/sl]");
599 MHCamEvent evt1( 2, "PedRmsFLG", "RMS from Extractor applied to ped.;;\\sigma_{p} [cnts/sl]");
600 MHCamEvent evt2( 0, "Extra'd", "Extracted Signal;;S [cnts/sl]");
601 // MHCamEvent evt3(4, "PedPhot", "Calibrated Pedestal;;P [phe]");
602 MHCamEvent evt4( 5, "PedRMS", "Calibrated RMS from Extractor applied to ped.;;\\sigma_{p} [phe]");
603 MHCamEvent evt5( 0, "Interp'd", "Interpolated Signal scaled with A/A_{0};;S [phe]");
604 MHCamEvent evt6(102, "Unsuitable", "Fraction of unsuitable events per Pixel;;[1]");
605 // MHCamEvent evt7( 6, "Times", "Calibrated Arrival Time;;T [fadc sl]");
606 MHCamEvent evt8( 0, "Conv", "Calibration Conv. Factors;;[phe/fadc cnts]");
607 MHCamEvent evt9( 7, "PulsePos", "Pulse Position of cosmics (>50phe);;T");
608 MHCamEvent evtR( 4, "HiLoCal", "Hi-/Lo-Gain ratio;;Ratio");
609 MHCamEvent evtO( 7, "HiLoOff", "Lo-/Hi-Gain Offset;;Offset");
610 evt2.SetErrorSpread(kFALSE);
611 evt5.SetErrorSpread(kFALSE);
612 evt6.SetErrorSpread(kFALSE);
613 evt6.SetThreshold();
614
615 MFillH fill0(&evt0, "MPedestalFundamental", "FillPedFLG");
616 MFillH fill1(&evt1, "MPedestalFromExtractorRndm", "FillPedRmsFLG");
617 MFillH fill2(&evt2, "MExtractedSignalCam", "FillExtracted");
618 // MFillH fill3(&evt3, "MPedPhotFundamental", "FillPedPhot");
619 MFillH fill4(&evt4, "MPedPhotFromExtractorRndm", "FillPedRMS");
620 MFillH fill5(&evt5, "MSignalCam", "FillInterpolated");
621 MFillH fill6(&evt6, "MBadPixelsCam", "FillUnsuitable");
622 // MFillH fill7(&evt7, "MSignalCam", "FillTimes");
623 MFillH fill8(&evt8, "MCalibrationIntensityConstCam", "FillConv");
624 MFillH fill9(&evt9, "MSignalCam", "FillPulse");
625 MFillH fillR(&evtR, "MExtractedSignalCam", "FillHiLoCal");
626 MFillH fillO(&evtO, "MArrivalTimeCam", "FillHiLoOff");
627
628 MTaskEnv fillflorian("FinalFantasy");
629 fillflorian.SetDefault();
630
631 // The second rule is for the case reading raw-files!
632 MWriteRootFile write(2, fname, fOverwrite?"RECREATE":"NEW");
633 // Run Header
634 write.AddContainer("MRawRunHeader", "RunHeaders");
635 write.AddContainer("MBadPixelsCam", "RunHeaders");
636 write.AddContainer("MGeomCam", "RunHeaders");
637 // Monte Carlo Headers
638 write.AddContainer("MMcRunHeader", "RunHeaders", kFALSE);
639 write.AddContainer("MMcFadcHeader", "RunHeaders", kFALSE);
640 write.AddContainer("MMcTrigHeader", "RunHeaders", kFALSE);
641 write.AddContainer("MMcConfigRunHeader", "RunHeaders", kFALSE);
642 write.AddContainer("MMcCorsikaRunHeader", "RunHeaders", kFALSE);
643 // Monte Carlo
644 write.AddContainer("MMcEvt", "Events", kFALSE);
645 write.AddContainer("MMcTrig", "Events", kFALSE);
646 // Data tree
647 write.AddContainer("MSignalCam", "Events");
648 // write.AddContainer("MPedPhotFundamental", "Events");
649 write.AddContainer("MPedPhotFromExtractor", "Events");
650 write.AddContainer("MPedPhotFromExtractorRndm", "Events");
651 write.AddContainer("MTime", "Events", kFALSE);
652 write.AddContainer("MRawEvtHeader", "Events");
653 // Slow-Control: Current-tree
654 write.AddContainer("MTimeCurrents", "Currents", kFALSE);
655 write.AddContainer("MCameraDC", "Currents", kFALSE);
656 write.AddContainer("MReportCurrents", "Currents", kFALSE);
657 // Slow-Control: Camera-tree
658 write.AddContainer("MReportCamera", "Camera", kFALSE);
659 write.AddContainer("MTimeCamera", "Camera", kFALSE);
660 write.AddContainer("MCameraAUX", "Camera", kFALSE);
661 write.AddContainer("MCameraCalibration", "Camera", kFALSE);
662 write.AddContainer("MCameraCooling", "Camera", kFALSE);
663 write.AddContainer("MCameraHV", "Camera", kFALSE);
664 write.AddContainer("MCameraLV", "Camera", kFALSE);
665 write.AddContainer("MCameraLids", "Camera", kFALSE);
666 // Slow-Control: Trigger-tree
667 write.AddContainer("MReportTrigger", "Trigger", kFALSE);
668 write.AddContainer("MTimeTrigger", "Trigger", kFALSE);
669 // Slow-Control: Drive-tree
670 write.AddContainer("MReportDrive", "Drive", kFALSE);
671 write.AddContainer("MTimeDrive", "Drive", kFALSE);
672 // Slow-Control: Central Control-tree
673 write.AddContainer("MReportCC", "CC", kFALSE);
674 write.AddContainer("MTimeCC", "CC", kFALSE);
675
676 // Write the special MC tree
677 MWriteRootFile writemc(2, fname, fOverwrite?"RECREATE":"NEW");
678 writemc.SetName("WriteMC");
679 writemc.AddContainer("MMcEvtBasic", "OriginalMC");
680
681 // Now setup tasklist for events
682 MTaskList tlist2;
683
684 tlist2.AddToList(&caldec);
685 tlist2.AddToList(&decode);
686 tlist2.AddToList(&apply);
687 tlist2.AddToList(&merge);
688 tlist2.AddToList(&pedsub);
689 tlist2.AddToList(&ftp2);
690 tlist2.AddToList(&pedlo1);
691 tlist2.AddToList(&pedlo2);
692 tlist2.AddToList(&pedlo3);
693
694 //-----------------------------------------------------------
695
696 MTaskList tlist3;
697 tlist3.SetFilter(&fcalib);
698
699 if (fIsInterlaced)
700 {
701 tlist2.AddToList(&fcalib);
702 tlist2.AddToList(&tlist3);
703
704 if (IsUsePINDiode())
705 tlist3.AddToList(&pinext);
706 if (IsUseBlindPixel())
707 tlist3.AddToList(&bldext);
708 tlist3.AddToList(&taskenv3);
709 tlist3.AddToList(&pacalc);
710 tlist3.AddToList(&filcam);
711 if (fIsRelTimesUpdate)
712 tlist3.AddToList(&filtme);
713 if (IsUseBlindPixel())
714 tlist3.AddToList(&filbnd);
715 if (IsUsePINDiode())
716 tlist3.AddToList(&filpin);
717 tlist3.AddToList(&chcalc);
718 if (fIsRelTimesUpdate)
719 tlist3.AddToList(&recalc);
720 }
721
722 // Continue for all non-cosmic events
723 tlist2.AddToList(&fill0); // fill pedestal events
724 tlist2.AddToList(&fill1); // fill pedestal events
725 if (!extractor1->HasLoGain())
726 tlist2.AddToList(&contftp2); // remove pedestal events from processing
727 tlist2.AddToList(&conttp); // remove calib events from processing
728 if (extractor1)
729 tlist2.AddToList(&taskenv1);
730 if (extractor2)
731 tlist2.AddToList(&taskenv2);
732 tlist2.AddToList(&contcos);
733 /*
734 if (fIsHiLoCalibration)
735 {
736 plist.AddToList(&hilocam);
737 tlist2.AddToList(&filhil);
738 }
739
740 if (fIsPulsePosCheck)
741 {
742 plist.AddToList(&hpulcam);
743 tlist2.AddToList(&filpul);
744 }
745 */
746 tlist2.AddToList(&fill2);
747 tlist2.AddToList(&calib);
748 if (extractor2 || extractor1->InheritsFrom("MExtractTimeAndCharge"))
749 tlist2.AddToList(&caltm);
750 if (fIsInterlaced)
751 tlist2.AddToList(&fill8);
752 tlist2.AddToList(&bpcal);
753 tlist2.AddToList(&treat);
754 tlist2.AddToList(&fill6);
755 // tlist2.AddToList(&fill3);
756 tlist2.AddToList(&fill4);
757 tlist2.AddToList(&fill5);
758 //if (extractor2 || extractor1->InheritsFrom("MExtractTimeAndCharge"))
759 // tlist2.AddToList(&fill7);
760 tlist2.AddToList(&fill9);
761 if (extractor1->HasLoGain())
762 {
763 tlist2.AddToList(&fillR);
764 tlist2.AddToList(&fillO);
765 }
766 tlist2.AddToList(&fillflorian);
767
768 // Setup List for Drive-tree
769 //MPointingPosCalc pcalc;
770
771 // Now setup main tasklist
772 tlist.AddToList(read);
773
774 if (IsUseMC())
775 {
776 tlist.AddToList(&writemc);
777 tlist.AddToList(&contmc);
778 }
779
780 //if (IsUseRootData())
781 // tlist2.AddToList(&pextr);
782 tlist.AddToList(&tlist2, IsUseRootData() ? "Events" : "All");
783
784 //if (IsUseMC())
785 // tlist.AddToList(&pcalc, "Drive");
786
787 tlist.AddToList(&write);
788
789 // Create and setup the eventloop
790 MEvtLoop evtloop(fName);
791 evtloop.SetParList(&plist);
792 evtloop.SetDisplay(fDisplay);
793 evtloop.SetLogStream(fLog);
794 if (!SetupEnv(evtloop))
795 return kFALSE;
796
797 // Execute first analysis
798 const Bool_t rc = evtloop.Eventloop(fMaxEvents);
799
800 // make sure owned object are deleted
801 if (extractor1 && extractor1!=fExtractor)
802 delete extractor1;
803 if (extractor2)
804 delete extractor2;
805 if (extractor3)
806 delete extractor3;
807
808 // return if job failed
809 if (!rc)
810 {
811 *fLog << err << GetDescriptor() << ": Failed." << endl;
812 return kFALSE;
813 }
814
815 // if everything went ok write and display result
816 DisplayResult(plist);
817
818 /*
819 if (fIsPixelCheck)
820 {
821 if (fIsPulsePosCheck)
822 hpulcam[fCheckedPixId].DrawClone("");
823
824 //if (fIsHiLoCalibration)
825 // hilocam[fCheckedPixId].DrawClone("");
826 }
827 */
828 interlacedcont.Add(&pulcam);
829
830 //if (fIsHiLoCalibration)
831 // interlacedcont.Add(&hilcam);
832
833 //if (fIsPulsePosCheck)
834 // interlacedcont.Add(plist.FindObject("MHCalibrationPulseTimeCam"));
835
836 //if (fIsHiLoCalibration)
837 // interlacedcont.Add(plist.FindObject("MHCalibrationHiLoCam"));
838
839 if (!WriteResult(interlacedcont))
840 return kFALSE;
841
842 // return if job went ok
843 *fLog << all << GetDescriptor() << ": Done." << endl;
844 *fLog << endl << endl;
845
846 return kTRUE;
847}
848
849
850void MJCalibrateSignal::DisplayResult(MParList &plist)
851{
852 if (!fDisplay || !fIsHiLoCalibration)
853 return;
854 /*
855 MCalibrationHiLoCam *hcam = (MCalibrationHiLoCam*)plist.FindObject("MCalibrationHiLoCam");
856 MGeomCam *geom = (MGeomCam*)plist.FindObject("MGeomCam");
857 if (!hcam || !geom)
858 return;
859
860 // Create histograms to display
861 MHCamera disp1(*geom, "HiLoConv", "Ratio Amplification HiGain vs. LoGain (Charges)");
862 MHCamera disp2(*geom, "HiLoDiff", "Arrival Time Diff. HiGain vs. LoGain (Times)");
863
864 disp1.SetCamContent(*hcam, 0);
865 disp1.SetCamError( *hcam, 1);
866 disp2.SetCamContent(*hcam, 5);
867 disp2.SetCamError( *hcam, 6);
868
869 disp1.SetYTitle("R [1]");
870 disp2.SetYTitle("\\Delta T [FADC sl.]");
871
872 TCanvas &c1 = fDisplay->AddTab("HiLoConv");
873 c1.Divide(2,3);
874
875 disp1.CamDraw(c1, 1, 2, 1);
876 disp2.CamDraw(c1, 2, 2, 1);
877*/
878}
879
Note: See TracBrowser for help on using the repository browser.