Context Navigation

source: trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc@ 5863

Visit:

Last change on this file since 5863 was 5843, checked in by gaug, 20 years ago
* empty log message *
File size: 30.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): Markus Gaug 11/2003 <mailto:markus@ifae.es>
19	!
20	! Copyright: MAGIC Software Development, 2000-2004
21	!
22	!
23	\* ======================================================================== */
24	/////////////////////////////////////////////////////////////////////////////
25	//
26	// MCalibrationChargeCam
27	//
28	// Storage container for charge calibration results from the signal distribution
29	// fits (see MHCalibrationChargeCam and MHCalibrationChargePix), the calculation
30	// of reduced sigmas and number of photo-electrons (this class) and conversion
31	// factors sum FADC slices to photo-electrons (see MCalibrationChargeCalc)
32	//
33	// Individual pixels have to be cast when retrieved e.g.:
34	// MCalibrationChargePix &avpix = (MCalibrationChargePix&)(*fChargeCam)[i]
35	//
36	// Averaged values over one whole area index (e.g. inner or outer pixels for
37	// the MAGIC camera), can be retrieved via:
38	// MCalibrationChargePix &avpix = (MCalibrationChargePix&)fChargeCam->GetAverageArea(i)
39	//
40	// Averaged values over one whole camera sector can be retrieved via:
41	// MCalibrationChargePix &avpix = (MCalibrationChargePix&)fChargeCam->GetAverageSector(i)
42	//
43	// Note the averageing has been done on an event-by-event basis. Resulting
44	// Sigma's of the Gauss fit have been multiplied with the square root of the number
45	// of involved pixels in order to make a direct comparison possible with the mean of
46	// sigmas.
47	//
48	// Final numbers of uncalibrated or unreliable pixels can be retrieved via the commands:
49	// GetNumUncalibrated(aidx) and GetNumUnreliable(aidx) where aidx is the area index (0 for
50	// inner and 1 for outer pixels in the MAGIC camera).
51	//
52	// The following "calibration" constants are used for the calibration of each pixel
53	// (see MCalibrate):
54	//
55	// - MCalibrationQEPix::GetMeanConvFADC2Phe(): The mean conversion factor from the
56	// summed FADC slices to the number of photo-electrons (in first order independent
57	// on colour and intensity)
58	// - MCalibrationQEPix::GetMeanFFactorFADC2Phot(): The mean F-Factor of the total
59	// readout chain dividing the signal to noise of the incoming number of photons
60	// (= sqrt(number photons)) by the signal to noise of the outgoing summed FADC slices
61	// signal (= MCalibrationPix::GetMean() / MCalibrationChargePix::GetRSigma() )
62	//
63	// The following calibration constants can be retrieved directly from this class:
64	//
65	// - GetConversionFactorFFactor ( Int_t idx, Float_t &mean, Float_t &err, Float_t &sigma );
66	//
67	// where:
68	// - idx is the pixel software ID
69	// - "mean" is the mean conversion constant, to be multiplied with the retrieved signal
70	// in order to get a calibrated number of PHOTONS.
71	// - "err" is the pure statistical uncertainty about the MEAN
72	// - "sigma", if mulitplied with the square root of signal, gives the approximate sigma of the
73	// retrieved mean number of incident Cherenkov photons.
74	//
75	// Note, Conversion is ALWAYS (included the F-Factor method!) from summed FADC slices to PHOTONS.
76	//
77	// See also: MCalibrationChargePix, MCalibrationChargeCalc, MCalibrationQECam
78	// MHCalibrationChargePix, MHCalibrationChargeCam
79	// MCalibrationBlindPix, MCalibrationChargePINDiode
80	//
81	/////////////////////////////////////////////////////////////////////////////
82	#include "MCalibrationChargeCam.h"
83	#include "MCalibrationChargePix.h"
84
85	#include <TOrdCollection.h>
86	#include <TH1D.h>
87	#include <TF1.h>
88
89	#include "MLog.h"
90	#include "MLogManip.h"
91
92	#include "MGeomCam.h"
93	#include "MGeomPix.h"
94
95	#include "MBadPixelsCam.h"
96	#include "MBadPixelsPix.h"
97
98	#include "MCalibrationQECam.h"
99	#include "MCalibrationQEPix.h"
100
101	#include "MHCamera.h"
102
103	ClassImp(MCalibrationChargeCam);
104
105	using namespace std;
106	// --------------------------------------------------------------------------
107	//
108	// Default constructor.
109	//
110	// Calls:
111	// - Clear()
112	//
113	MCalibrationChargeCam::MCalibrationChargeCam(const char name, const char title)
114	{
115	fName = name ? name : "MCalibrationChargeCam";
116	fTitle = title ? title : "Storage container for the Calibration Information in the camera";
117
118	Clear();
119	}
120
121	/*
122	// --------------------------------------------------------------------------
123	//
124	// Creates new MCalibrationCam only for the averaged areas:
125	// the rest has to be retrieved directly, e.g. via:
126	//
127	TObject MCalibrationChargeCam::Clone(const char ) const
128	{
129
130	//
131	// FIXME, this might be done faster and more elegant, by direct copy.
132	//
133	MCalibrationChargeCam *cam = new MCalibrationChargeCam(fName,fTitle);
134
135	cam->fNumUnsuitable = fNumUnsuitable;
136	cam->fNumUnreliable = fNumUnreliable;
137	cam->fNumHiGainFADCSlices = fNumHiGainFADCSlices;
138	cam->fNumLoGainFADCSlices = fNumLoGainFADCSlices;
139	cam->fPulserColor = fPulserColor;
140
141	cam->fFlags = fFlags;
142
143	cam->fNumPhotonsBlindPixelMethod = fNumPhotonsBlindPixelMethod;
144	cam->fNumPhotonsFFactorMethod = fNumPhotonsFFactorMethod;
145	cam->fNumPhotonsPINDiodeMethod = fNumPhotonsPINDiodeMethod;
146	cam->fNumPhotonsBlindPixelMethodErr = fNumPhotonsBlindPixelMethodErr;
147	cam->fNumPhotonsFFactorMethodErr = fNumPhotonsFFactorMethodErr;
148	cam->fNumPhotonsPINDiodeMethodErr = fNumPhotonsPINDiodeMethodErr;
149
150	for (Int_t i=0; i<GetSize(); i++)
151	cam->fPixels->AddAt((*this)[i].Clone(),i);
152
153	for (Int_t i=0; i<GetAverageAreas(); i++)
154	{
155	cam->fAverageAreas->AddAt(GetAverageArea(i).Clone(),i);
156	cam->fAverageBadAreas->AddAt(GetAverageBadArea(i).Clone(),i);
157	}
158	for (Int_t i=0; i<GetAverageSectors(); i++)
159	{
160	cam->fAverageSectors->AddAt(GetAverageSector(i).Clone(),i);
161	cam->fAverageBadSectors->AddAt(GetAverageBadSector(i).Clone(),i);
162	}
163
164	return cam;
165	}
166	*/
167
168	// -------------------------------------------------------------------
169	//
170	// Add MCalibrationChargePix's in the ranges from - to to fPixels
171	//
172	void MCalibrationChargeCam::Add(const UInt_t a, const UInt_t b)
173	{
174	for (UInt_t i=a; i<b; i++)
175	{
176	fPixels->AddAt(new MCalibrationChargePix,i);
177	(*this)[i].SetPixId(i);
178	}
179	}
180
181	// -------------------------------------------------------------------
182	//
183	// Add MCalibrationChargePix's in the ranges from - to to fAverageAreas
184	//
185	void MCalibrationChargeCam::AddArea(const UInt_t a, const UInt_t b)
186	{
187	for (UInt_t i=a; i<b; i++)
188	{
189	fAverageAreas->AddAt(new MCalibrationChargePix,i);
190	GetAverageArea(i).SetPixId(i);
191	}
192	}
193
194	// -------------------------------------------------------------------
195	//
196	// Add MCalibrationChargePix's in the ranges from - to to fAverageSectors
197	//
198	void MCalibrationChargeCam::AddSector(const UInt_t a, const UInt_t b)
199	{
200	for (UInt_t i=a; i<b; i++)
201	{
202	fAverageSectors->AddAt(new MCalibrationChargePix,i);
203	GetAverageSector(i).SetPixId(i);
204	}
205	}
206
207
208	// --------------------------------------
209	//
210	// Sets all variable to 0.
211	// Sets all flags to kFALSE
212	// Calls MCalibrationCam::Clear()
213	//
214	void MCalibrationChargeCam::Clear(Option_t *o)
215	{
216
217	SetFFactorMethodValid ( kFALSE );
218
219	fNumPhotonsBlindPixelMethod = 0.;
220	fNumPhotonsFFactorMethod = 0.;
221	fNumPhotonsPINDiodeMethod = 0.;
222	fNumPhotonsBlindPixelMethodErr = 0.;
223	fNumPhotonsFFactorMethodErr = 0.;
224	fNumPhotonsPINDiodeMethodErr = 0.;
225
226	MCalibrationCam::Clear();
227
228	return;
229	}
230
231	// -----------------------------------------------
232	//
233	// Sets the kFFactorMethodValid bit from outside
234	//
235	void MCalibrationChargeCam::SetFFactorMethodValid(const Bool_t b)
236	{
237	b ? SETBIT(fFlags, kFFactorMethodValid) : CLRBIT(fFlags, kFFactorMethodValid);
238	}
239
240
241	// --------------------------------------------------------------------------
242	//
243	// Test bit kFFactorMethodValid
244	//
245	Bool_t MCalibrationChargeCam::IsFFactorMethodValid() const
246	{
247	return TESTBIT(fFlags,kFFactorMethodValid);
248	}
249
250	// --------------------------------------------------------------------------
251	//
252	// Print first the well fitted pixels
253	// and then the ones which are not FitValid
254	//
255	void MCalibrationChargeCam::Print(Option_t *o) const
256	{
257
258	*fLog << all << GetDescriptor() << ":" << endl;
259	int id = 0;
260
261	*fLog << all << "Calibrated pixels:" << endl;
262	*fLog << all << endl;
263
264	TIter Next(fPixels);
265	MCalibrationChargePix *pix;
266	while ((pix=(MCalibrationChargePix*)Next()))
267	{
268
269	if (!pix->IsExcluded())
270	{
271
272	*fLog << all
273	<< Form("%s%3i","Pixel: ",pix->GetPixId())
274	<< Form("%s%4.2f%s%4.2f"," Ped.Rms: ",pix->GetPedRms(),"+-",pix->GetPedRmsErr())
275	<< Form("%s%4.2f%s%4.2f"," Charge: " ,pix->GetConvertedMean(),"+-",pix->GetConvertedSigma())
276	<< Form("%s%4.2f%s%4.2f"," Red.Sigma: ",pix->GetConvertedRSigma(),"+-",pix->GetConvertedRSigmaErr())
277	<< Form("%s%4.2f%s%4.2f"," Num.Phes: ",pix->GetPheFFactorMethod(),"+-",pix->GetPheFFactorMethodErr())
278	<< Form("%s%4.2f%s%4.2f"," Conv.FADC2Phe: ",pix->GetMeanConvFADC2Phe(),"+-",pix->GetMeanConvFADC2PheErr())
279	<< " Saturated? :" << pix->IsHiGainSaturation()
280	<< endl;
281	id++;
282	}
283	}
284
285	*fLog << all << id << " pixels" << endl;
286	id = 0;
287
288
289	*fLog << all << endl;
290	*fLog << all << "Excluded pixels:" << endl;
291	*fLog << all << endl;
292
293	id = 0;
294
295	TIter Next4(fPixels);
296	while ((pix=(MCalibrationChargePix*)Next4()))
297	{
298	if (pix->IsExcluded())
299	{
300	*fLog << all << pix->GetPixId() << " ";
301	id++;
302
303	if (!(id % 25))
304	*fLog << endl;
305	}
306	}
307
308	*fLog << endl;
309	*fLog << all << id << " Excluded pixels " << endl;
310	*fLog << endl;
311
312	*fLog << all << endl;
313	*fLog << all << "Averaged Areas:" << endl;
314	*fLog << all << endl;
315
316	TIter Next5(fAverageAreas);
317	while ((pix=(MCalibrationChargePix*)Next5()))
318	{
319	*fLog << all << Form("%s%3i","Area Idx: ",pix->GetPixId())
320	<< " Ped. Rms: " << pix->GetPedRms() << " +- " << pix->GetPedRmsErr()
321	<< " Mean signal: " << pix->GetMean() << " +- " << pix->GetMeanErr()
322	<< " Sigma signal: " << pix->GetSigma() << " +- "<< pix->GetSigmaErr()
323	<< " Reduced Sigma: " << pix->GetRSigma()
324	<< " Nr Phe's: " << pix->GetPheFFactorMethod()
325	<< endl;
326	}
327
328	*fLog << all << endl;
329	*fLog << all << "Averaged Sectors:" << endl;
330	*fLog << all << endl;
331
332	TIter Next6(fAverageSectors);
333	while ((pix=(MCalibrationChargePix*)Next6()))
334	{
335	*fLog << all << Form("%s%3i","Sector: ",pix->GetPixId())
336	<< " Ped. Rms: " << pix->GetPedRms() << " +- " << pix->GetPedRmsErr()
337	<< " Mean signal: " << pix->GetMean() << " +- " << pix->GetMeanErr()
338	<< " Sigma signal: " << pix->GetSigma() << " +- "<< pix->GetSigmaErr()
339	<< " Reduced Sigma: " << pix->GetRSigma()
340	<< " Nr Phe's: " << pix->GetPheFFactorMethod()
341	<< endl;
342	}
343	*fLog << all << endl;
344	}
345
346
347	// --------------------------------------------------------------------------
348	//
349	// The types are as follows:
350	//
351	// Fitted values:
352	// ==============
353	//
354	// 0: Fitted Charge (see MCalibrationPix::GetMean())
355	// 1: Error of fitted Charge (see MCalibrationPix::GetMeanErr())
356	// 2: Sigma of fitted Charge (see MCalibrationPix::GetSigma())
357	// 3: Error of Sigma of fitted Charge (see MCalibrationPix::GetSigmaErr())
358	//
359	// Useful variables derived from the fit results:
360	// =============================================
361	//
362	// 4: Probability Gauss fit Charge distribution (see MCalibrationPix::GetProb())
363	// 5: Reduced Sigma of fitted Charge (see MCalibrationChargePix::GetRSigma())
364	// 6: Error Reduced Sigma of fitted Charge (see MCalibrationChargePix::GetRSigmaErr())
365	// 7: Reduced Sigma per Charge (see MCalibrationChargePix::GetRSigmaPerCharge())
366	// 8: Error of Reduced Sigma per Charge (see MCalibrationChargePix::GetRSigmaPerChargeErr())
367	//
368	// Results of the F-Factor calibration Method:
369	// ===========================================
370	//
371	// 9: Nr. Photo-electrons from F-Factor Method (see MCalibrationChargePix::GetPheFFactorMethod())
372	// 10: Error Nr. Photo-el. from F-Factor Method (see MCalibrationChargePix::GetPheFFactorMethodErr())
373	// 11: Conversion factor from F-Factor Method (see MCalibrationChargePix::GetMeanConvFADC2Phe()
374	// 12: Error conv. factor from F-Factor Method (see MCalibrationChargePix::GetMeanConvFADC2PheErr()
375	// 13: Overall F-Factor from F-Factor Method (see MCalibrationChargePix::GetMeanFFactorFADC2Phot()
376	// 14: Error F-Factor from F-Factor Method (see MCalibrationChargePix::GetMeanFFactorFADC2PhotErr()
377	// 15: Pixels valid calibration F-Factor-Method (see MCalibrationChargePix::IsFFactorMethodValid())
378	//
379	// Results of the Low-Gain vs. High-Gain Conversion:
380	// =================================================
381	//
382	// 16: Mean Signal Hi-Gain / Mean Signal Lo-Gain (see MCalibrationPix::GetHiLoMeansDivided())
383	// 17: Error Signal High-Gain / Signal Low Gain (see MCalibrationPix::GetHiLoMeansDividedErr())
384	// 18: Sigma High-Gain / Sigma Low Gain (see MCalibrationPix::GetHiLoSigmasDivided())
385	// 19: Error Sigma High-Gain / Sigma Low Gain (see MCalibrationPix::GetHiLoSigmasDividedErr())
386	//
387	// Localized defects:
388	// ==================
389	//
390	// 20: Excluded Pixels
391	// 21: Number of pickup events in the Hi Gain (see MCalibrationPix::GetHiGainNumPickup())
392	// 22: Number of pickup events in the Lo Gain (see MCalibrationPix::GetLoGainNumPickup())
393	// 23: Number of blackout events in the Hi Gain (see MCalibrationPix::GetHiGainNumBlackout())
394	// 24: Number of blackout events in the Lo Gain (see MCalibrationPix::GetLoGainNumBlackout())
395	//
396	// Other classifications of pixels:
397	// ================================
398	//
399	// 25: Pixels with saturated High-Gain (see MCalibrationPix::IsHiGainSaturation())
400	//
401	// Calculated absolute arrival times (very low precision!):
402	// ========================================================
403	//
404	// 26: Absolute Arrival time of the signal (see MCalibrationChargePix::GetAbsTimeMean())
405	// 27: RMS Ab. Arrival time of the signal (see MCalibrationChargePix::GetAbsTimeRms())
406	//
407	// Used Pedestals:
408	// ===============
409	//
410	// 28: Pedestal for entire signal extr. range (see MCalibrationChargePix::Ped())
411	// 29: Error Pedestal entire signal extr. range (see MCalibrationChargePix::PedErr())
412	// 30: Ped. RMS entire signal extraction range (see MCalibrationChargePix::PedRms())
413	// 31: Error Ped. RMS entire signal extr. range (see MCalibrationChargePix::PedRmsErr())
414	//
415	// Special variables (for data check):
416	// ====================================
417	//
418	// 32: HiGain RMS divided by Mean for every pixel (with inclusion of the excluded pixels)
419	//
420	Bool_t MCalibrationChargeCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
421	{
422
423	if (idx > GetSize())
424	return kFALSE;
425
426	Float_t area = cam[idx].GetA();
427
428	if (area == 0)
429	return kFALSE;
430
431	MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[idx];
432
433	switch (type)
434	{
435	case 0:
436	if (pix.IsExcluded())
437	return kFALSE;
438	val = pix.GetMean();
439	break;
440	case 1:
441	if (pix.IsExcluded())
442	return kFALSE;
443	val = pix.GetMeanErr();
444	break;
445	case 2:
446	if (pix.IsExcluded())
447	return kFALSE;
448	val = pix.GetSigma();
449	break;
450	case 3:
451	if (pix.IsExcluded())
452	return kFALSE;
453	val = pix.GetSigmaErr();
454	break;
455	case 4:
456	if (pix.IsExcluded())
457	return kFALSE;
458	val = pix.GetProb();
459	break;
460	case 5:
461	if (pix.IsExcluded())
462	return kFALSE;
463	if (pix.GetRSigma() == -1.)
464	return kFALSE;
465	val = pix.GetRSigma();
466	break;
467	case 6:
468	if (pix.IsExcluded())
469	return kFALSE;
470	if (pix.GetRSigma() == -1.)
471	return kFALSE;
472	val = pix.GetRSigmaErr();
473	break;
474	case 7:
475	if (pix.IsExcluded())
476	return kFALSE;
477	val = pix.GetRSigmaPerCharge();
478	break;
479	case 8:
480	if (pix.IsExcluded())
481	return kFALSE;
482	val = pix.GetRSigmaPerChargeErr();
483	break;
484	case 9:
485	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
486	return kFALSE;
487	val = pix.GetPheFFactorMethod();
488	break;
489	case 10:
490	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
491	return kFALSE;
492	val = pix.GetPheFFactorMethodErr();
493	break;
494	case 11:
495	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
496	return kFALSE;
497	val = pix.GetMeanConvFADC2Phe();
498	break;
499	case 12:
500	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
501	return kFALSE;
502	val = pix.GetMeanConvFADC2PheErr();
503	break;
504	case 13:
505	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
506	return kFALSE;
507	val = pix.GetMeanFFactorFADC2Phot();
508	break;
509	case 14:
510	if (pix.IsExcluded() \|\| !pix.IsFFactorMethodValid())
511	return kFALSE;
512	val = pix.GetMeanFFactorFADC2PhotErr();
513	break;
514	case 15:
515	if (pix.IsExcluded())
516	return kFALSE;
517	if (pix.IsFFactorMethodValid())
518	val = 1;
519	else
520	return kFALSE;
521	break;
522	case 16:
523	if (pix.IsExcluded())
524	return kFALSE;
525	val = pix.GetHiLoMeansDivided();
526	break;
527	case 17:
528	if (pix.IsExcluded())
529	return kFALSE;
530	val = pix.GetHiLoMeansDividedErr();
531	break;
532	case 18:
533	if (pix.IsExcluded())
534	return kFALSE;
535	val = pix.GetHiLoSigmasDivided();
536	break;
537	case 19:
538	if (pix.IsExcluded())
539	return kFALSE;
540	val = pix.GetHiLoSigmasDividedErr();
541	break;
542	case 20:
543	if (pix.IsExcluded())
544	val = 1.;
545	else
546	return kFALSE;
547	break;
548	case 21:
549	if (pix.IsExcluded())
550	return kFALSE;
551	val = pix.GetHiGainNumPickup();
552	break;
553	case 22:
554	if (pix.IsExcluded())
555	return kFALSE;
556	val = pix.GetLoGainNumPickup();
557	break;
558	case 23:
559	if (pix.IsExcluded())
560	return kFALSE;
561	val = pix.GetHiGainNumBlackout();
562	break;
563	case 24:
564	if (pix.IsExcluded())
565	return kFALSE;
566	val = pix.GetLoGainNumBlackout();
567	break;
568	case 25:
569	if (pix.IsExcluded())
570	return kFALSE;
571	val = pix.IsHiGainSaturation();
572	break;
573	case 26:
574	if (pix.IsExcluded())
575	return kFALSE;
576	val = pix.GetAbsTimeMean();
577	break;
578	case 27:
579	if (pix.IsExcluded())
580	return kFALSE;
581	val = pix.GetAbsTimeRms();
582	break;
583	case 28:
584	if (pix.IsExcluded())
585	return kFALSE;
586	val = pix.GetPed();
587	break;
588	case 29:
589	if (pix.IsExcluded())
590	return kFALSE;
591	val = pix.GetPedErr();
592	break;
593	case 30:
594	if (pix.IsExcluded())
595	return kFALSE;
596	val = pix.GetPedRms();
597	break;
598	case 31:
599	if (pix.IsExcluded())
600	return kFALSE;
601	val = pix.GetPedErr()/2.;
602	break;
603	case 32:
604	val = pix.GetMean() == 0. ? 0. : pix.GetRms()/pix.GetMean();
605	break;
606	default:
607	return kFALSE;
608	}
609
610	return val!=-1.;
611	}
612
613
614
615	Bool_t MCalibrationChargeCam::GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t &err, Float_t &ffactor)
616	{
617
618	MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[ipx];
619
620	Float_t conv = pix.GetMeanConvFADC2Phe();
621
622	if (conv < 0.)
623	return kFALSE;
624
625	mean = conv;
626	err = pix.GetMeanConvFADC2PheErr();
627	ffactor = pix.GetMeanFFactorFADC2Phot();
628
629	return kTRUE;
630	}
631
632
633	// --------------------------------------------------------------------------
634	//
635	// Calculates the average conversion factor FADC counts to photons for pixel sizes.
636	// The geometry container is used to get the necessary
637	// geometry information (area index). The MCalibrationQECam container is necessary for
638	// the quantum efficiency information.
639	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
640	// in the calculation of the size average.
641	//
642	// Returns a TArrayF of dimension 2:
643	// arr[0]: averaged conversion factors (default: -1.)
644	// arr[1]: Error (rms) of averaged conversion factors (default: 0.)
645	//
646	TArrayF MCalibrationChargeCam::GetAveragedConvFADC2PhotPerArea ( const MGeomCam &geom, const MCalibrationQECam &qecam,
647	const UInt_t ai, MBadPixelsCam *bad)
648	{
649
650	const Int_t np = GetSize();
651
652	Double_t mean = 0.;
653	Double_t mean2 = 0.;
654	Int_t nr = 0;
655
656	MHCamera convcam(geom,"ConvFactors","Conversion Factors;Conv Factor [phot/FADC cnts];channels");
657
658	for (int i=0; i<np; i++)
659	{
660	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
661	continue;
662
663	if (bad && (*bad)[i].IsUncalibrated(MBadPixelsPix::kDeviatingNumPhes))
664	continue;
665
666	const UInt_t aidx = geom[i].GetAidx();
667
668	if (ai != aidx)
669	continue;
670
671	const MCalibrationQEPix &qepix = (MCalibrationQEPix&)qecam[i];
672	if (!qepix.IsAverageQEFFactorAvailable())
673	continue;
674
675	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
676	const Float_t conv = pix.GetMeanConvFADC2Phe()/qepix.GetQECascadesFFactor();
677
678	mean += conv;
679	mean2 += conv*conv;
680	nr ++;
681
682	convcam.Fill(i,conv);
683	convcam.SetUsed(i);
684	}
685
686	Float_t mn = nr ? mean/nr : -1.;
687	Float_t sg = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0.;
688
689	const Int_t aidx = (Int_t)ai;
690
691	TH1D *h = convcam.ProjectionS(TArrayI(),TArrayI(1,&aidx),"_py",750);
692	h->SetDirectory(NULL);
693	h->Fit("gaus","Q");
694	TF1 *fit = h->GetFunction("gaus");
695
696	Float_t ci2 = fit->GetChisquare();
697	Float_t sigma = fit->GetParameter(2);
698
699	if (ci2 > 500. \|\| sigma > sg)
700	{
701	h->Fit("gaus","QM");
702	fit = h->GetFunction("gaus");
703
704	ci2 = fit->GetChisquare();
705	sigma = fit->GetParameter(2);
706	}
707
708	const Int_t ndf = fit->GetNDF();
709
710	if (ci2 < 500. && sigma < sg && ndf > 2)
711	{
712	mn = fit->GetParameter(1);
713	sg = sigma;
714	}
715
716	*fLog << inf << "Conversion Factors to photons area idx: " << aidx << ":" << endl;
717	*fLog << inf << "Mean: " << Form("%4.3f",mn)
718	<< "+-" << Form("%4.3f",fit->GetParError(1))
719	<< " Sigma: " << Form("%4.3f",sg) << "+-" << Form("%4.3f",fit->GetParError(2))
720	<< " Chisquare: " << Form("%4.3f",fit->GetChisquare()) << " NDF : " << ndf << endl;
721
722	delete h;
723	gROOT->GetListOfFunctions()->Remove(fit);
724
725	TArrayF arr(2);
726	arr[0] = mn;
727	arr[1] = sg;
728
729	return arr;
730	}
731
732	// --------------------------------------------------------------------------
733	//
734	// Calculates the average conversion factor FADC counts to photons for camera sectors.
735	// The geometry container is used to get the necessary
736	// geometry information (area index). The MCalibrationQECam container is necessary for
737	// the quantum efficiency information.
738	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
739	// in the calculation of the size average.
740	//
741	// Returns a TArrayF of dimension 2:
742	// arr[0]: averaged conversion factors (default: -1.)
743	// arr[1]: Error (rms) of averaged conversion factors (default: 0.)
744	//
745	TArrayF MCalibrationChargeCam::GetAveragedConvFADC2PhotPerSector( const MGeomCam &geom, const MCalibrationQECam &qecam,
746	const UInt_t sec, MBadPixelsCam *bad)
747	{
748	const Int_t np = GetSize();
749
750	Double_t mean = 0.;
751	Double_t mean2 = 0.;
752	Int_t nr = 0;
753
754	for (int i=0; i<np; i++)
755	{
756	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
757	continue;
758
759	const UInt_t sector = geom[i].GetSector();
760
761	if (sector != sec)
762	continue;
763
764	const MCalibrationQEPix &qepix = (MCalibrationQEPix&)qecam[i];
765	if (!qepix.IsAverageQEFFactorAvailable())
766	continue;
767
768	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
769	const Float_t conv = pix.GetMeanConvFADC2Phe();
770	const Float_t qe = qepix.GetQECascadesFFactor();
771
772	mean += conv/qe;
773	mean2 += conv*conv/qe/qe;
774	nr ++;
775
776	}
777
778	TArrayF arr(2);
779	arr[0] = nr ? mean/nr : -1;
780	arr[1] = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0;
781	return arr;
782	}
783
784	// --------------------------------------------------------------------------
785	//
786	// Calculates the average mean arrival times for pixel sizes.
787	// The geometry container is used to get the necessary
788	// geometry information (area index).
789	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
790	// in the calculation of the size average.
791	//
792	// Returns a TArrayF of dimension 2:
793	// arr[0]: averaged mean arrival times (default: -1.)
794	// arr[1]: Error (rms) of averaged mean arrival times (default: 0.)
795	//
796	TArrayF MCalibrationChargeCam::GetAveragedArrivalTimeMeanPerArea(const MGeomCam &geom,
797	const UInt_t ai, MBadPixelsCam *bad)
798	{
799
800	const Int_t np = GetSize();
801
802	Double_t mean = 0.;
803	Double_t mean2 = 0.;
804	Int_t nr = 0;
805
806	for (int i=0; i<np; i++)
807	{
808	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
809	continue;
810
811	const UInt_t aidx = geom[i].GetAidx();
812
813	if (ai != aidx)
814	continue;
815
816	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
817	const Float_t time = pix.GetAbsTimeMean();
818
819	mean += time ;
820	mean2 += time*time;
821	nr ++;
822
823	}
824
825	TArrayF arr(2);
826	arr[0] = nr ? mean/nr : -1;
827	arr[1] = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0;
828	return arr;
829	}
830
831	// --------------------------------------------------------------------------
832	//
833	// Calculates the average mean arrival times for camera sectors.
834	// The geometry container is used to get the necessary
835	// geometry information (area index).
836	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
837	// in the calculation of the size average.
838	//
839	// Returns a TArrayF of dimension 2:
840	// arr[0]: averaged mean arrival times (default: -1.)
841	// arr[1]: Error (rms) of averaged mean arrival times (default: 0.)
842	//
843	TArrayF MCalibrationChargeCam::GetAveragedArrivalTimeMeanPerSector(const MGeomCam &geom,
844	const UInt_t sec, MBadPixelsCam *bad)
845	{
846	const Int_t np = GetSize();
847
848	Double_t mean = 0.;
849	Double_t mean2 = 0.;
850	Int_t nr = 0;
851
852	for (int i=0; i<np; i++)
853	{
854	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
855	continue;
856
857	const UInt_t sector = geom[i].GetSector();
858
859	if (sector != sec)
860	continue;
861
862	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
863	const Float_t time = pix.GetAbsTimeMean();
864
865	mean += time;
866	mean2 += time*time;
867	nr ++;
868
869	}
870
871	TArrayF arr(2);
872	arr[0] = nr ? mean/nr : -1;
873	arr[1] = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0;
874	return arr;
875	}
876
877	// --------------------------------------------------------------------------
878	//
879	// Calculates the average arrival time RMSs for pixel sizes.
880	// The geometry container is used to get the necessary
881	// geometry information (area index).
882	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
883	// in the calculation of the size average.
884	//
885	// Returns a TArrayF of dimension 2:
886	// arr[0]: averaged arrival time RMSs (default: -1.)
887	// arr[1]: Error (rms) of averaged arrival time RMSs (default: 0.)
888	//
889	TArrayF MCalibrationChargeCam::GetAveragedArrivalTimeRmsPerArea ( const MGeomCam &geom,
890	const UInt_t ai, MBadPixelsCam *bad)
891	{
892
893	const Int_t np = GetSize();
894
895	Double_t mean = 0.;
896	Double_t mean2 = 0.;
897	Int_t nr = 0;
898
899	for (int i=0; i<np; i++)
900	{
901	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
902	continue;
903
904	const UInt_t aidx = geom[i].GetAidx();
905
906	if (ai != aidx)
907	continue;
908
909	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
910	const Float_t rms = pix.GetAbsTimeRms();
911
912	mean += rms;
913	mean2 += rms*rms;
914	nr ++;
915
916	}
917
918	TArrayF arr(2);
919	arr[0] = nr ? mean/nr : -1;
920	arr[1] = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0;
921	return arr;
922	}
923
924	// --------------------------------------------------------------------------
925	//
926	// Calculates the average arrival time RMSs for camera sectors.
927	// The geometry container is used to get the necessary
928	// geometry information (area index).
929	// If the bad pixel container is given all pixels which have the flag 'kUnsuitableRun' are ignored
930	// in the calculation of the size average.
931	//
932	// Returns a TArrayF of dimension 2:
933	// arr[0]: averaged arrival time RMSs (default: -1.)
934	// arr[1]: Error (rms) of averaged arrival time RMSs (default: 0.)
935	//
936	TArrayF MCalibrationChargeCam::GetAveragedArrivalTimeRmsPerSector( const MGeomCam &geom, const UInt_t sec, MBadPixelsCam *bad)
937	{
938	const Int_t np = GetSize();
939
940	Double_t mean = 0.;
941	Double_t mean2 = 0.;
942	Int_t nr = 0;
943
944	for (int i=0; i<np; i++)
945	{
946	if (bad && (*bad)[i].IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
947	continue;
948
949	const UInt_t sector = geom[i].GetSector();
950
951	if (sector != sec)
952	continue;
953
954	const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*this)[i];
955	const Float_t rms = pix.GetAbsTimeRms();
956
957	mean += rms;
958	mean2 += rms*rms;
959	nr ++;
960	}
961
962	TArrayF arr(2);
963	arr[0] = nr ? mean/nr : -1;
964	arr[1] = nr>1 ? TMath::Sqrt((mean2 - mean*mean/nr)/(nr-1)) : 0;
965
966	return arr;
967	}

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

Download in other formats: