- Timestamp:
- 10/17/06 18:18:40 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.