Context Navigation

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

Visit:

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

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

Download in other formats: