Context Navigation

source: trunk/MagicSoft/Mars/mjobs/MJSimulation.cc@ 9343

Visit:

Last change on this file since 9343 was 9343, checked in by tbretz, 16 years ago
* empty log message *
File size: 24.2 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/2009 <mailto:tbretz@astro.uni-wuerzburg.de>
19	!
20	! Copyright: MAGIC Software Development, 2000-2009
21	!
22	!
23	\* ======================================================================== */
24
25	/////////////////////////////////////////////////////////////////////////////
26	//
27	// MJSimulation
28	//
29	//
30	// Force reading a corsika file even if the footer (RUNE-section) is missing
31	// by setting fForceMode to kTRUE or from the resource file by
32	//
33	// ForceMode: Yes
34	//
35	//
36	// To switch off the simulation of the camera electronics, use:
37	//
38	// Camera: Off
39	//
40	// Note, that the border of the camera and the propertied of the cones
41	// are still simulated (simply because it is fast). Furthermore, this
42	// switches off the trigger for the output, i.e. all data which deposits
43	// at least one photon in at least one pixel is written to the output file.
44	//
45	//
46	// In case of a pedestal or calibration run the artificial trigger can
47	// be "switched off" and the cosmics trrigger "switched on" by setting
48	// fForceTrigger to kTRUE or from the resource file by
49	//
50	// ForceTrigger: Yes
51	//
52	//
53	/////////////////////////////////////////////////////////////////////////////
54	#include "MJSimulation.h"
55
56	#include <TEnv.h>
57	#include <TCanvas.h>
58
59	// Core
60	#include "MLog.h"
61	#include "MLogManip.h"
62
63	#include "MArgs.h"
64	//#include "MDirIter.h"
65	#include "MParList.h"
66	#include "MTaskList.h"
67	#include "MEvtLoop.h"
68
69	#include "MStatusDisplay.h"
70
71	// Tasks
72	#include "MCorsikaRead.h"
73	#include "MContinue.h"
74	#include "MFillH.h"
75	#include "MGeomApply.h"
76	#include "MHillasCalc.h"
77	#include "MImgCleanStd.h"
78	#include "MWriteRootFile.h"
79
80	#include "MSimMMCS.h"
81	#include "MSimAbsorption.h"
82	#include "MSimReflector.h"
83	#include "MSimPointingPos.h"
84	#include "MSimPSF.h"
85	#include "MSimGeomCam.h"
86	#include "MSimSignalCam.h"
87	#include "MSimAPD.h"
88	#include "MSimExcessNoise.h"
89	#include "MSimCamera.h"
90	#include "MSimTrigger.h"
91	#include "MSimReadout.h"
92	#include "MSimRandomPhotons.h"
93	#include "MSimBundlePhotons.h"
94	#include "MSimCalibrationSignal.h"
95
96	// Histograms
97	#include "MBinning.h"
98
99	#include "MHn.h"
100	#include "MHCamera.h"
101	#include "MHCamEvent.h"
102	#include "MHPhotonEvent.h"
103
104	// Container
105	#include "MRawRunHeader.h"
106	#include "MParameters.h"
107	#include "MReflector.h"
108	#include "MParEnv.h"
109	#include "MPulseShape.h"
110	#include "MGeomCam.h"
111
112	#include "MPedestalCam.h"
113	#include "MPedestalPix.h"
114
115	ClassImp(MJSimulation);
116
117	using namespace std;
118
119	// --------------------------------------------------------------------------
120	//
121	// Default constructor.
122	//
123	// Sets fRuns to 0, fExtractor to NULL, fDataCheck to kFALSE
124	//
125	MJSimulation::MJSimulation(const char name, const char title)
126	: fForceMode(kFALSE), fCamera(kTRUE) ,fForceTrigger(kFALSE)
127	{
128	fName = name ? name : "MJSimulation";
129	fTitle = title ? title : "Standard analysis and reconstruction";
130	}
131
132	Bool_t MJSimulation::CheckEnvLocal()
133	{
134	fForceMode = GetEnv("ForceMode", fForceMode);
135	fForceTrigger = GetEnv("ForceTrigger", fForceTrigger);
136	fCamera = GetEnv("Camera", fCamera);
137
138	return kTRUE;
139	}
140
141	Bool_t MJSimulation::WriteResult(const MParList &plist)
142	{
143	if (fPathOut.IsNull())
144	{
145	*fLog << inf << "No output path specified via SetPathOut - no output written." << endl;
146	return kTRUE;
147	}
148
149	TObjArray cont;
150	cont.Add(const_cast<TEnv*>(GetEnv()));
151	//cont.Add(const_cast<MSequence*>(&fSequence));
152
153	cont.Add(plist.FindObject("MPulseShape"));
154	cont.Add(plist.FindObject("Reflector"));
155	cont.Add(plist.FindObject("MGeomCam"));
156	cont.Add(plist.FindObject("GeomCones"));
157
158	if (fDisplay)
159	{
160	// TString title = "-- Reflector: ";
161	// title += fSequence.GetSequence();
162	// title += " --";
163	// fDisplay->SetTitle(title, kFALSE);
164	fDisplay->SetTitle("Ceres", kFALSE);
165
166	cont.Add(fDisplay);
167	}
168
169	// const TString oname = Form("reflector%08d.root", fSequence.GetSequence());
170	const TString oname = "ceres.root";
171	return WriteContainer(cont, oname, "RECREATE");
172	}
173
174	void MJSimulation::SetupHist(MHn &hist) const
175	{
176	hist.AddHist("MCorsikaEvtHeader.fTotalEnergy");
177	hist.InitName("Energy");
178	hist.InitTitle("Energy;E [GeV]");
179	hist.SetLog(kTRUE, kTRUE, kFALSE);
180
181	hist.AddHist("MPhotonEvent.GetNumExternal");
182	hist.InitName("Size");
183	hist.InitTitle("Size;S [#]");
184	hist.SetLog(kTRUE, kTRUE, kFALSE);
185
186	hist.AddHist("-MCorsikaEvtHeader.fX/100","-MCorsikaEvtHeader.fY/100");
187	hist.SetDrawOption("colz");
188	hist.InitName("Impact;Impact;Impact");
189	hist.InitTitle("Impact;West <--> East [m];South <--> North [m]");
190
191	hist.AddHist("MCorsikaEvtHeader.fFirstInteractionHeight/100000");
192	hist.InitName("Height");
193	hist.InitTitle("FirstInteractionHeight;h [km]");
194
195	hist.AddHist("MCorsikaEvtHeader.fAzTMath::RadToDeg()", "MCorsikaEvtHeader.fZdTMath::RadToDeg()");
196	hist.InitName("SkyOrigin;Az;Zd");
197	hist.InitTitle("Sky Origin;Az [deg];Zd [deg]");
198	hist.SetDrawOption("colz");
199
200	TString sin2 = "sin(MCorsikaEvtHeader.fZd)sin(MCorsikaRunHeader.fZdMinTMath::DegToRad())";
201	TString cos2 = "cos(MCorsikaEvtHeader.fZd)cos(MCorsikaRunHeader.fZdMinTMath::DegToRad())";
202	TString cos = "cos(MCorsikaEvtHeader.fAz-MCorsikaRunHeader.fAzMin*TMath::DegToRad())";
203
204	TString form = "acos("+sin2+""+cos+"+"+cos2+")TMath::RadToDeg()";
205
206	hist.AddHist(form);
207	hist.InitName("ViewCone");
208	hist.InitTitle("Incident Angle;\\alpha [\\deg]");
209	}
210
211	void MJSimulation::SetupCommonFileStructure(MWriteRootFile &write) const
212	{
213	// Common run headers
214	write.AddContainer("MCorsikaRunHeader", "RunHeaders", kFALSE);
215	write.AddContainer("MRawRunHeader", "RunHeaders");
216	write.AddContainer("MGeomCam", "RunHeaders");
217	write.AddContainer("MMcRunHeader", "RunHeaders");
218
219	// Common events
220	write.AddContainer("MCorsikaEvtHeader", "Events", kFALSE);
221	write.AddContainer("MRawEvtHeader", "Events");
222	write.AddContainer("MMcEvt", "Events");
223	}
224
225	Bool_t MJSimulation::Process(const MArgs &args)
226	{
227	/*
228	if (!fSequence.IsValid())
229	{
230	*fLog << err << "ERROR - Sequence invalid!" << endl;
231	return kFALSE;
232	}
233	*/
234
235	//if (!HasWritePermission(GetPathOut()))
236	// return kFALSE;
237
238	*fLog << inf;
239	fLog->Separator(GetDescriptor());
240
241	if (!CheckEnv())
242	return kFALSE;
243
244	*fLog << warn << "FIXME: Monte Carlo simulation: Sequences not supported yet.";
245	//*fLog << fSequence.GetFileName() << endl;
246	*fLog << endl;
247
248	// --------------------------------------------------------------------------------
249
250	//MDirIter iter;
251	//if (fSequence.GetRuns(iter, MSequence::kCalibrated)<=0)
252	// return kFALSE;
253
254	// Setup Parlist
255	MParList plist;
256	plist.AddToList(this); // take care of fDisplay!
257
258	// setup TaskList
259	MTaskList tasks;
260	plist.AddToList(&tasks);
261
262	// --------------------------------------------------------------------------------
263
264	// FIXME: Allow empty GeomCones!
265	MParEnv env0("Reflector");
266	MParEnv env1("GeomCones");
267	MParEnv env2("MGeomCam");
268	env0.SetClassName("MReflector");
269	env1.SetClassName("MGeomCam");
270	env2.SetClassName("MGeomCam");
271	plist.AddToList(&env0);
272	plist.AddToList(&env1);
273	plist.AddToList(&env2);
274
275	plist.FindCreateObj("MPedPhotCam", "MPedPhotFromExtractorRndm");
276
277	MPulseShape shape;
278	plist.AddToList(&shape);
279
280	// * FIXME * FIXME * FIXME *
281	plist.FindCreateObj("MMcRunHeader");
282
283	MRawRunHeader header;
284	header.SetValidMagicNumber();
285	//header.InitFadcType(3);
286
287	header.SetRunType(MRawRunHeader::kRTMonteCarlo\|MRawRunHeader::kRTData);
288	if (args.GetNumArguments()==1)
289	{
290	if (!args.GetArgumentStr(0).CompareTo("pedestal", TString::kIgnoreCase))
291	{
292	header.SetRunType(MRawRunHeader::kRTMonteCarlo\|MRawRunHeader::kRTPedestal);
293	header.SetSourceInfo("Pedestal");
294	}
295	if (!args.GetArgumentStr(0).CompareTo("calibration", TString::kIgnoreCase))
296	{
297	header.SetRunType(MRawRunHeader::kRTMonteCarlo\|MRawRunHeader::kRTCalibration);
298	header.SetSourceInfo("Calibration");
299	}
300	if (!args.GetArgumentStr(0).CompareTo("pointrun", TString::kIgnoreCase))
301	header.SetRunType(MRawRunHeader::kRTMonteCarlo\|MRawRunHeader::kRTPointRun);
302	}
303
304	// FIXME: Move to MSimPointingPos, MSimCalibrationSignal
305	// Can we use this as input for MSimPointingPos?
306	header.SetObservation("On", "MonteCarlo");
307	plist.AddToList(&header);
308	// ++++++++ FIXME FIXME FIXME +++++++++++++
309
310	/*
311	MPedestalCam pedcam;
312	pedcam.Init(geomcam.GetNumPixels());
313	for (UInt_t i=0; i<geomcam.GetNumPixels(); i++)
314	pedcam[i].Set(128./header.GetScale(), 22.5/header.GetScale());
315	plist.AddToList(&pedcam);
316	*/
317
318	// -------------------------------------------------------------------
319
320	MCorsikaRead read;
321	read.SetForceMode(fForceMode);
322
323	for (int i=0; i<args.GetNumArguments();i ++)
324	read.AddFile(args.GetArgumentStr(i));
325
326	MSimMMCS simmmcs;
327
328	MSimAbsorption absapd("PhotonDetectionEfficiency");
329	MSimAbsorption absmir("MirrorReflectivity");
330	MSimAbsorption cones("ConesAngularAcceptance");
331	cones.SetUseTheta();
332
333	MSimPointingPos pointing;
334
335	MSimReflector reflect;
336	reflect.SetNameGeomCam("GeomCones");
337	// MSimStarField stars;
338
339	MContinue cont1("MPhotonEvent.GetNumPhotons<1", "ContEmpty1");
340	MContinue cont2("MPhotonEvent.GetNumPhotons<1", "ContEmpty2");
341
342	// -------------------------------------------------------------------
343
344	MBinning binse( 100, 1, 100000, "BinningEnergy", "log");
345	MBinning binss( 100, 1, 10000000, "BinningSize", "log");
346	MBinning binsi( 50, -250, 250, "BinningImpact");
347	MBinning binsh( 150, 0, 50, "BinningHeight");
348	MBinning binsaz(720, -360, 360, "BinningAz");
349	MBinning binszd( 70, 0, 70, "BinningZd");
350	MBinning binsvc( 45, 0, 9, "BinningViewCone");
351	MBinning binsew(150, 0, 50, "BinningEvtWidth");
352	MBinning binstr("BinningTrigPos");
353
354	plist.AddToList(&binse);
355	plist.AddToList(&binss);
356	plist.AddToList(&binsi);
357	plist.AddToList(&binsh);
358	plist.AddToList(&binszd);
359	plist.AddToList(&binsaz);
360	plist.AddToList(&binsvc);
361	plist.AddToList(&binsew);
362	plist.AddToList(&binstr);
363
364	MHn mhn1, mhn2, mhn3;
365	SetupHist(mhn1);
366	SetupHist(mhn2);
367	SetupHist(mhn3);
368
369	MH3 mhtp("TriggerPos.fVal-IntendedPulsePos.fVal-MPulseShape.GetPulseWidth");
370	mhtp.SetName("TrigPos");
371	mhtp.SetTitle("Trigger position w.r.t. the first photon hitting a detector");
372
373	MH3 mhew("MPhotonStatistics.fLength");
374	mhew.SetName("EvtWidth");
375	mhew.SetTitle("Time between first and last photon hitting a detector");
376
377	MH3 mhed("MPhotonStatistics.fTimeMedDev");
378	mhed.SetName("EvtMedDev;EvtWidth");
379
380	MFillH fillh1(&mhn1, "", "FillCorsika");
381	MFillH fillh2(&mhn2, "", "FillH2");
382	MFillH fillh3(&mhn3, "", "FillH3");
383	MFillH filltp(&mhtp, "", "FillTriggerPos");
384	MFillH fillew(&mhew, "", "FillEvtWidth");
385	MFillH filled(&mhed, "", "FillMedDev");
386	fillh1.SetNameTab("Corsika", "Distribution as simulated by Corsika");
387	fillh2.SetNameTab("Detectable", "Distribution of events hit the detector");
388	fillh3.SetNameTab("Triggered", "Distribution of triggered events");
389	filltp.SetNameTab("TrigPos", "TriggerPosition w.r.t the first photon");
390	fillew.SetNameTab("EvtWidth", "Time between first and last photon hitting a detector");
391	filled.SetNameTab("MedDev", "Time between first and last photon hitting a detector");
392
393	MHPhotonEvent planeG(1); // Get from MaxImpact
394	MHPhotonEvent plane0(2); // Get from MReflector
395	MHPhotonEvent plane1(2);
396	MHPhotonEvent plane2(2);
397	MHPhotonEvent plane3(2);
398	MHPhotonEvent plane4(2);
399	MHPhotonEvent planeF1(3); // Get from MGeomCam
400	MHPhotonEvent planeF2(header.IsPointRun()?4:3); // Get from MGeomCam
401
402	//MHPSF psf;
403
404	MFillH fillG(&planeG, "MPhotonEvent", "FillGround");
405	MFillH fill0(&plane0, "MirrorPlane0", "FillReflector");
406	//MFillH fill1(&plane1, "MirrorPlane1", "FillCamShadow");
407	MFillH fill2(&plane2, "MirrorPlane2", "FillCandidates");
408	MFillH fill3(&plane3, "MirrorPlane3", "FillReflected");
409	MFillH fill4(&plane4, "MirrorPlane4", "FillFocal");
410	MFillH fillF1(&planeF1, "MPhotonEvent", "FillCamera");
411	MFillH fillF2(&planeF2, "MPhotonEvent", "FillCones");
412	//MFillH fillP(&psf, "MPhotonEvent", "FillPSF");
413
414	fillG.SetNameTab("Ground", "Photon distribution at ground");
415	fill0.SetNameTab("Reflector", "Photon distribution at reflector plane w/o camera shadow");
416	//fill1.SetNameTab("CamShadow", "Photon distribution at reflector plane w/ camera shadow");
417	fill2.SetNameTab("Candidates", "'Can hit' photon distribution at reflector plane w/ camera shadow");
418	fill3.SetNameTab("Reflected", "Photon distribution after reflector projected to reflector plane");
419	fill4.SetNameTab("Focal", "Photon distribution at focal plane");
420	fillF1.SetNameTab("Camera", "Photon distribution which hit the detector");
421	fillF2.SetNameTab("Cones", "Photon distribution after cones");
422	//fillP.SetNameTab("PSF", "Photon distribution after cones for all mirrors");
423
424	// -------------------------------------------------------------------
425
426	// FIXME: Pedestal and Calibration runs should get P and C
427	const TString rule1(Form("s/cer([0-9]+)/%s\\/00$1_R_MonteCarlo.root/", Esc(fPathOut).Data()));
428	const TString rule2(Form("s/cer([0-9]+)/%s\\/00$1_I_MonteCarlo.root/", Esc(fPathOut).Data()));
429	const TString rule3(Form("s/cer([0-9]+)/%s\\/00$1_D_MonteCarlo.root/", Esc(fPathOut).Data()));
430
431	MWriteRootFile write3a( 2, rule3, fOverwrite?"RECREATE":"NEW", "Camera file");
432	MWriteRootFile write3b( 2, rule3, fOverwrite?"RECREATE":"NEW", "Camera file");
433	MWriteRootFile write2a( 2, rule2, fOverwrite?"RECREATE":"NEW", "Signal file");
434	MWriteRootFile write2b( 2, rule2, fOverwrite?"RECREATE":"NEW", "Signal file");
435	MWriteRootFile write1a( 2, rule1, fOverwrite?"RECREATE":"NEW", "Reflector file");
436	MWriteRootFile write1b( 2, rule1, fOverwrite?"RECREATE":"NEW", "Reflector file");
437
438	write1a.SetName("WriteRefData");
439	write1b.SetName("WriteRefMC");
440	write2a.SetName("WriteSigData");
441	write2b.SetName("WriteSigMC");
442	write3a.SetName("WriteCamData");
443	write3b.SetName("WriteCamMC");
444
445	SetupCommonFileStructure(write1a);
446	SetupCommonFileStructure(write2a);
447	SetupCommonFileStructure(write3a);
448
449	// R: Dedicated file structureedicated events
450	write1a.AddContainer("MPhotonEvent", "Events");
451
452	// I: Dedicated file structureedicated events
453	write2a.AddContainer("MPedPhotFromExtractorRndm", "RunHeaders"); // FIXME: Needed for the signal files to be display in MARS
454	write2a.AddContainer("MSignalCam", "Events");
455
456	// D: Dedicated file structureedicated events
457	write3a.AddContainer("ElectronicNoise", "RunHeaders");
458	write3a.AddContainer("MRawEvtData", "Events");
459
460	// Basic MC data
461	write1b.AddContainer("MMcEvtBasic", "OriginalMC");
462	write2b.AddContainer("MMcEvtBasic", "OriginalMC");
463	write3b.AddContainer("MMcEvtBasic", "OriginalMC");
464
465	// -------------------------------------------------------------------
466
467	MGeomApply apply;
468
469	MSimPSF simpsf;
470
471	MSimGeomCam simgeom;
472	simgeom.SetNameGeomCam("GeomCones");
473
474	MSimCalibrationSignal simcal;
475	simcal.SetNameGeomCam("GeomCones");
476
477	// Dark Current: 4MHz = 0.004 GHz
478	// Light of night sky at La Palma: ~ 0.2 ph / cm^2 / sr / ns
479	// 5deg camera: 6e-3 sr
480	// NSB = 0.2*6e-3
481
482	MSimRandomPhotons simnsb;
483	// simnsb.SetFreq(0.005, 0.004); // 5MHz/cm^2, 4MHz
484	simnsb.SetFreq(0, 0.04); // 40MHz fixed for each APD
485	simnsb.SetNameGeomCam("GeomCones");
486
487	// FIXME: How do we handle star-light? We need the rate but we also
488	// need to process it through the mirror!
489
490	MSimAPD simapd;
491	simapd.SetFreq(0.04); // Must be identical to MSimRanodmPhotons!!
492	simapd.SetNameGeomCam("GeomCones");
493
494	MSimExcessNoise simexcnoise;
495	MSimBundlePhotons simsum;
496	MSimSignalCam simsignal;
497	MSimCamera simcam;
498	MSimTrigger simtrig;
499	MSimReadout simdaq;
500
501	MContinue conttrig("TriggerPos.fVal<0", "ContTrigger");
502
503	MParameterD pulpos("IntendedPulsePos");
504	// pulpos.SetVal(40); // [ns]
505	plist.AddToList(&pulpos);
506
507	MParameterD trigger("TriggerPos");
508	trigger.SetVal(0);
509	plist.AddToList(&trigger);
510
511	// -------------------------------------------------------------------
512
513	// FIXME: Remove isolated pixels
514	MImgCleanStd clean(7, 4.5);
515	clean.SetMethod(MImgCleanStd::kAbsolute);
516
517	//clean.SetNamePedPhotCam("MPedPhotFromExtractorRndm");
518
519	MHillasCalc hcalc;
520	hcalc.Disable(MHillasCalc::kCalcConc);
521
522	MHCamEvent evt0a(/10/0, "Signal", "Average signal;;S [ph]");
523	MHCamEvent evt0c(/10/0, "MaxSignal", "Maximum signal;;S [ph]");
524	MHCamEvent evt0d(/11/8, "ArrTm", "Time after first photon;;T [ns]");
525	evt0a.SetErrorSpread(kFALSE);
526	evt0c.SetCollectMax();
527
528	MFillH fillx0a(&evt0a, "MSignalCam", "FillAvgSignal");
529	MFillH fillx0c(&evt0c, "MSignalCam", "FillMaxSignal");
530	MFillH fillx0d(&evt0d, "MSignalCam", "FillArrTm");
531	MFillH fillx1("MHHillas", "MHillas", "FillHillas");
532	//MFillH fillx2("MHHillasExt", "", "FillHillasExt");
533	MFillH fillx3("MHHillasSrc", "MHillasSrc", "FillHillasSrc");
534	//MFillH fillx4("MHImagePar", "MImagePar", "FillImagePar");
535	//MFillH fillx5("MHNewImagePar", "MNewImagePar", "FillNewImagePar");
536
537	// -------------------------------------------------------------------
538
539	if (header.IsDataRun())
540	{
541	tasks.AddToList(&read);
542	tasks.AddToList(&pointing);
543	tasks.AddToList(&simmmcs);
544	if (!fPathOut.IsNull() && !HasNullOut())
545	{
546	tasks.AddToList(&write1b);
547	tasks.AddToList(&write2b);
548	if (fCamera)
549	tasks.AddToList(&write3b);
550	}
551	// if (header.IsPointRun())
552	// tasks.AddToList(&stars);
553	//tasks.AddToList(&print);
554	tasks.AddToList(&fillh1);
555	tasks.AddToList(&fillG);
556	if (!header.IsPointRun())
557	{
558	tasks.AddToList(&absapd);
559	tasks.AddToList(&absmir);
560	}
561	tasks.AddToList(&reflect);
562	if (!header.IsPointRun())
563	{
564	tasks.AddToList(&fill0);
565	//tasks.AddToList(&fill1);
566	tasks.AddToList(&fill2);
567	tasks.AddToList(&fill3);
568	tasks.AddToList(&fill4);
569	tasks.AddToList(&fillF1);
570	}
571	tasks.AddToList(&cones);
572	tasks.AddToList(&fillF2);
573	//if (header.IsPointRun())
574	// tasks.AddToList(&fillP);
575	tasks.AddToList(&cont1);
576	if (!header.IsPointRun())
577	tasks.AddToList(&fillh2);
578	}
579	// -------------------------------
580	tasks.AddToList(&apply);
581	if (header.IsDataRun())
582	{
583	tasks.AddToList(&simpsf);
584	tasks.AddToList(&simgeom);
585	tasks.AddToList(&cont2);
586	}
587	if (fCamera)
588	{
589	if (header.IsPedestalRun() \|\| header.IsCalibrationRun())
590	tasks.AddToList(&simcal);
591	tasks.AddToList(&simnsb);
592	tasks.AddToList(&simapd);
593	tasks.AddToList(&simexcnoise);
594	}
595	tasks.AddToList(&simsum);
596	if (fCamera)
597	{
598	tasks.AddToList(&simcam);
599	if (header.IsDataRun() \|\| fForceTrigger)
600	tasks.AddToList(&simtrig);
601	tasks.AddToList(&conttrig);
602	tasks.AddToList(&simdaq);
603	}
604	tasks.AddToList(&simsignal); // What do we do if signal<0?
605	if (!fPathOut.IsNull() && !HasNullOut())
606	{
607	tasks.AddToList(&write1a);
608	if (!header.IsPedestalRun())
609	tasks.AddToList(&write2a);
610	if (fCamera)
611	tasks.AddToList(&write3a);
612	}
613	// -------------------------------
614	if (fCamera)
615	{
616	// FIXME: MHCollectionArea Trigger Area!
617	if (header.IsDataRun())
618	tasks.AddToList(&fillh3);
619	tasks.AddToList(&filltp);
620	}
621	if (header.IsDataRun())
622	{
623	tasks.AddToList(&fillew);
624	tasks.AddToList(&filled);
625	}
626	if (!header.IsPedestalRun())
627	{
628	tasks.AddToList(&fillx0a);
629	tasks.AddToList(&fillx0c);
630	//tasks.AddToList(&clean);
631	tasks.AddToList(&hcalc);
632	tasks.AddToList(&fillx0d);
633	tasks.AddToList(&fillx1);
634	//tasks.AddToList(&fillx2);
635	tasks.AddToList(&fillx3);
636	//tasks.AddToList(&fillx4);
637	//tasks.AddToList(&fillx5);
638	}
639	//-------------------------------------------
640
641	tasks.SetAccelerator(MTask::kAccDontReset\|MTask::kAccDontTime);
642
643	// Create and setup the eventloop
644	MEvtLoop evtloop(fName);
645	evtloop.SetParList(&plist);
646	evtloop.SetDisplay(fDisplay);
647	evtloop.SetLogStream(&gLog);
648	if (!SetupEnv(evtloop))
649	return kFALSE;
650
651	// FIXME: If pedestal file given use the values from this file
652	//-------------------------------------------
653
654	if (binstr.IsDefault())
655	binstr.SetEdgesLin(150, -shape.GetPulseWidth(),
656	header.GetFreqSampling()/1000.*header.GetNumSamples()+shape.GetPulseWidth());
657
658	header.Print();
659
660	// FIXME: Display acceptance curves!
661
662	if (fDisplay)
663	{
664	TCanvas &c = fDisplay->AddTab("Info");
665	c.Divide(2,2);
666
667	c.cd(1);
668	gPad->SetBorderMode(0);
669	gPad->SetFrameBorderMode(0);
670	gPad->SetGridx();
671	gPad->SetGridy();
672	gROOT->SetSelectedPad(0);
673	shape.DrawClone()->SetBit(kCanDelete);
674
675	if (env0.GetCont() && (header.IsDataRun() \|\| header.IsPointRun()))
676	{
677	c.cd(3);
678	gPad->SetBorderMode(0);
679	gPad->SetFrameBorderMode(0);
680	gPad->SetGridx();
681	gPad->SetGridy();
682	gROOT->SetSelectedPad(0);
683	static_cast<MReflector*>(env0.GetCont())->DrawClone("line")->SetBit(kCanDelete);
684
685	}
686
687	if (fCamera)
688	{
689	if (env1.GetCont())
690	{
691	c.cd(2);
692	gPad->SetBorderMode(0);
693	gPad->SetFrameBorderMode(0);
694	gPad->SetGridx();
695	gPad->SetGridy();
696	gROOT->SetSelectedPad(0);
697	MHCamera c = new MHCamera(static_cast<MGeomCam&>(env1.GetCont()));
698	c->SetStats(kFALSE);
699	c->SetBit(MHCamera::kNoLegend);
700	c->SetBit(kCanDelete);
701	c->Draw();
702
703	}
704
705	if (env2.GetCont())
706	{
707	c.cd(4);
708	gPad->SetBorderMode(0);
709	gPad->SetFrameBorderMode(0);
710	gPad->SetGridx();
711	gPad->SetGridy();
712	gROOT->SetSelectedPad(0);
713	MHCamera c = new MHCamera(static_cast<MGeomCam&>(env2.GetCont()));
714	c->SetStats(kFALSE);
715	c->SetBit(MHCamera::kNoLegend);
716	c->SetBit(kCanDelete);
717	c->Draw();
718	}
719	}
720	}
721
722	//-------------------------------------------
723
724	// Execute first analysis
725	if (!evtloop.Eventloop(fMaxEvents))
726	{
727	gLog << err << GetDescriptor() << ": Failed." << endl;
728	return kFALSE;
729	}
730
731	//-------------------------------------------
732	// FIXME: Display Spline in tab
733	// FIXME: Display Reflector in tab
734	// FIXME: Display Routing(?) in tab
735	// FIXME: Display Camera(s) in tab
736	//-------------------------------------------
737
738	if (!WriteResult(plist))
739	return kFALSE;
740
741	*fLog << all << GetDescriptor() << ": Done." << endl << endl << endl;;
742
743	return kTRUE;
744	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: