Changeset 8106 for trunk/MagicSoft/Mars/mcalib
- Timestamp:
- 10/17/06 18:18:40 (18 years ago)
- Location:
- trunk/MagicSoft/Mars/mcalib
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/MCalibrationBlindPix.cc
r5701 r8106 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationBlindPix.cc,v 1.15 2006-10-17 17:15:59 tbretz Exp $ 3 ! -------------------------------------------------------------------------- 2 4 ! 3 5 ! * … … 375 377 { 376 378 377 const Float_t err = fAttErr[fColor];378 379 if ( err < 0.)380 return -1.; 381 382 return err*err*2.3;379 const Float_t ferr = fAttErr[fColor]; 380 381 if (ferr < 0.) 382 return -1.; 383 384 return ferr*ferr*2.3; 383 385 } 384 386 -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc
r7881 r8106 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCalc.cc,v 1.166 2006-10-17 17:15:59 tbretz Exp $ 3 ! -------------------------------------------------------------------------- 2 4 ! 3 5 ! * … … 1807 1809 1808 1810 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*chargecam)[i]; 1809 MCalibrationQEPix &q epix= (MCalibrationQEPix&) (*qecam) [i];1811 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (*qecam) [i]; 1810 1812 MBadPixelsPix &bad = (*badcam) [i]; 1811 1813 … … 1818 1820 const Float_t qerelvar = avphotrelvar + pix.GetPheFFactorMethodRelVar(); 1819 1821 1820 q epix.SetQEFFactor ( qe , fPulserColor );1821 q epix.SetQEFFactorVar ( qerelvar*qe*qe, fPulserColor );1822 q epix.SetFFactorMethodValid( kTRUE , fPulserColor );1823 1824 if (!q epix.UpdateFFactorMethod( qecam->GetPlexiglassQE() ))1822 qpix.SetQEFFactor ( qe , fPulserColor ); 1823 qpix.SetQEFFactorVar ( qerelvar*qe*qe, fPulserColor ); 1824 qpix.SetFFactorMethodValid( kTRUE , fPulserColor ); 1825 1826 if (!qpix.UpdateFFactorMethod( qecam->GetPlexiglassQE() )) 1825 1827 *fLog << warn << GetDescriptor() 1826 1828 << ": Cannot update Quantum efficiencies with the F-Factor Method" << endl; … … 1949 1951 { 1950 1952 1951 MCalibrationChargePix &pix 1952 MCalibrationQEPix &q epix = (MCalibrationQEPix&) (*qecam) [i];1953 MBadPixelsPix &bad 1953 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*chargecam)[i]; 1954 MCalibrationQEPix &qpix = (MCalibrationQEPix&) (*qecam) [i]; 1955 MBadPixelsPix &bad = (*badcam) [i]; 1954 1956 1955 1957 if (bad.IsUnsuitable(MBadPixelsPix::kUnsuitableRun)) 1956 1958 { 1957 q epix.SetFFactorMethodValid(kFALSE,fPulserColor);1959 qpix.SetFFactorMethodValid(kFALSE,fPulserColor); 1958 1960 pix.SetFFactorMethodValid(kFALSE); 1959 1961 continue; -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc
r7059 r8106 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCam.cc,v 1.69 2006-10-17 17:18:40 tbretz Exp $ 3 ! -------------------------------------------------------------------------- 2 4 ! 3 5 ! * … … 680 682 681 683 682 Bool_t MCalibrationChargeCam::GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t & err, Float_t &ffactor)684 Bool_t MCalibrationChargeCam::GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t &ferr, Float_t &ffactor) 683 685 { 684 686 … … 691 693 692 694 mean = conv; 693 err= pix.GetMeanConvFADC2PheErr();695 ferr = pix.GetMeanConvFADC2PheErr(); 694 696 ffactor = pix.GetMeanFFactorFADC2Phot(); 695 697 -
trunk/MagicSoft/Mars/mcalib/MCalibrationIntensityChargeCam.cc
r7188 r8106 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MCalibrationIntensityChargeCam.cc,v 1.23 2006-10-17 17:15:59 tbretz Exp $ 3 ! -------------------------------------------------------------------------- 2 4 ! 3 5 ! * … … 266 268 sigerr[cnt] = pheerr; 267 269 268 Double_t sig 270 Double_t sig1 = 0.; 269 271 Double_t sig2 = 0.; 270 272 Int_t num = 0; 271 273 272 for (Int_t i=0; i<cam->GetSize(); i++)273 { 274 const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)[ i];274 for (Int_t j=0; j<cam->GetSize(); j++) 275 { 276 const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)[j]; 275 277 // 276 278 // Don't use bad pixels … … 280 282 // 281 283 // 282 if (aidx != geom[ i].GetAidx())284 if (aidx != geom[j].GetAidx()) 283 285 continue; 284 286 285 sig 287 sig1 += pix.GetConvertedMean(); 286 288 sig2 += pix.GetConvertedMean() * pix.GetConvertedMean(); 287 289 num++; … … 290 292 if (num > 1) 291 293 { 292 sig 293 294 Double_t var = (sig2 - sig *sig*num) / (num-1);295 var /= sig *sig;294 sig1 /= num; 295 296 Double_t var = (sig2 - sig1*sig1*num) / (num-1); 297 var /= sig1*sig1; 296 298 var += pherelvar; 297 299 298 phepersig[cnt] = phe/sig ;300 phepersig[cnt] = phe/sig1; 299 301 if (var > 0.) 300 302 phepersigerr[cnt] = TMath::Sqrt(var) * phepersig[cnt]; … … 1057 1059 for (Int_t i=0; i<GetSize(); i++) 1058 1060 { 1059 MCalibrationChargeCam *c am = (MCalibrationChargeCam*)GetCam(i);1061 MCalibrationChargeCam *ccam = (MCalibrationChargeCam*)GetCam(i); 1060 1062 // 1061 1063 // Get the calibration pix from the calibration cam 1062 1064 // 1063 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*c am)[npix];1065 MCalibrationChargePix &pix = (MCalibrationChargePix&)(*ccam)[npix]; 1064 1066 // 1065 1067 // Don't use bad pixels … … 1096 1098 if (option.Contains("conversionfactor")) 1097 1099 { 1098 const MCalibrationChargePix &apix = (MCalibrationChargePix&)c am->GetAverageArea(0);1100 const MCalibrationChargePix &apix = (MCalibrationChargePix&)ccam->GetAverageArea(0); 1099 1101 pvar = apix.GetPheFFactorMethod()/pix.GetConvertedMean(); 1100 1102 }
Note:
See TracChangeset
for help on using the changeset viewer.