Changeset 5701 for trunk/MagicSoft/Mars/mcalib
- Timestamp:
- 01/05/05 14:29:59 (20 years ago)
- Location:
- trunk/MagicSoft/Mars/mcalib
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/MCalibrationBlindPix.cc
r5384 r5701 62 62 63 63 using namespace std; 64 64 65 const Float_t MCalibrationBlindPix::fgArea = 100; 65 66 const Float_t MCalibrationBlindPix::fgAttGreen = 1.97; … … 81 82 const Float_t MCalibrationBlindPix::fgCollEffCT1 = 0.90; 82 83 const Float_t MCalibrationBlindPix::fgCollEffErr = 0.05; 84 83 85 // -------------------------------------------------------------------------- 84 86 // … … 226 228 pix.fPixId = fPixId; 227 229 pix.fFlags = fFlags; 228 pix.fArea = fArea; 230 pix.fArea = fArea; 229 231 pix.fAreaErr = fAreaErr; 230 232 pix.fColor = fColor; … … 233 235 // Copy arrays 234 236 // 235 pix.fAtt = fAtt; 236 pix.fAttErr = fAttErr; 237 pix.fQE = fQE; 238 pix.fQEErr = fQEErr; 239 pix.fCollEff = fCollEff; 240 pix.fCollEffErr = fCollEffErr; 241 242 pix.fLambda = fLambda; 243 pix.fLambdaCheck = fLambdaCheck; 244 pix.fLambdaCheckErr = fLambdaCheckErr; 245 pix.fLambdaVar = fLambdaVar; 246 pix.fFluxInsidePlexiglass = fFluxInsidePlexiglass; 247 pix.fFluxInsidePlexiglassVar = fFluxInsidePlexiglassVar; 248 pix.fMu0 = fMu0; 249 pix.fMu0Err = fMu0Err; 250 pix.fMu1 = fMu1; 251 pix.fMu1Err = fMu1Err; 252 pix.fSigma0 = fSigma0; 253 pix.fSigma0Err = fSigma0Err; 254 pix.fSigma1 = fSigma1; 255 pix.fSigma1Err = fSigma1Err; 256 237 pix.fAtt = fAtt; 238 pix.fAttErr = fAttErr; 239 pix.fQE = fQE; 240 pix.fQEErr = fQEErr; 241 pix.fCollEff = fCollEff; 242 pix.fCollEffErr = fCollEffErr; 243 244 pix.fLambda = fLambda; 245 pix.fLambdaCheck = fLambdaCheck; 246 pix.fLambdaCheckErr = fLambdaCheckErr; 247 pix.fLambdaVar = fLambdaVar; 248 pix.fFluxInsidePlexiglass = fFluxInsidePlexiglass; 249 pix.fFluxInsidePlexiglassVar = fFluxInsidePlexiglassVar; 250 pix.fMu0 = fMu0; 251 pix.fMu0Err = fMu0Err; 252 pix.fMu1 = fMu1; 253 pix.fMu1Err = fMu1Err; 254 pix.fSigma0 = fSigma0; 255 pix.fSigma0Err = fSigma0Err; 256 pix.fSigma1 = fSigma1; 257 pix.fSigma1Err = fSigma1Err; 258 257 259 } 258 260 … … 523 525 524 526 *fLog << inf << GetDescriptor() 525 << ": Blind Pixel Nr. " << fPixId << ": Photon flux [ph/mm^2] inside Plexiglass: "526 << Form("%5.3f%s%5.3f",fFluxInsidePlexiglass," +- ",GetFluxInsidePlexiglassErr()) << endl;527 << ": BlindPix #" << fPixId << ": Flux [ph/mm^2] inside Plexiglass: " 528 << Form("%5.3f +- %5.3f",fFluxInsidePlexiglass, GetFluxInsidePlexiglassErr()) << endl; 527 529 528 530 return kTRUE; … … 533 535 534 536 *fLog << all << GetDescriptor() 535 << Form(" %s%3i","BlindPixel: ",GetPixId())536 << Form(" %s%4.2f%s%4.2f"," Lambda: ",GetLambda(),"+-",GetLambdaErr())537 << Form(" %s%4.2f%s%4.2f"," Mu0: " ,GetMu0(), "+-",GetMu0Err())538 << Form(" %s%4.2f%s%4.2f"," Mu1: " ,GetMu1(), "+-",GetMu1Err())539 << Form(" %s%4.2f%s%4.2f"," Sigma0: ",GetSigma0(),"+-",GetSigma0Err())540 << Form(" %s%4.2f%s%4.2f"," Sigma1: ",GetSigma1(),"+-",GetSigma1Err())537 << Form("BlindPixel: %3i" ,GetPixId()) 538 << Form(" Lambda: %4.2f +- %4.2f",GetLambda(),GetLambdaErr()) 539 << Form(" Mu0: %4.2f +- %4.2f" ,GetMu0(), GetMu0Err()) 540 << Form(" Mu1: %4.2f +- %4.2f" ,GetMu1(), GetMu1Err()) 541 << Form(" Sigma0: %4.2f +- %4.2f",GetSigma0(),GetSigma0Err()) 542 << Form(" Sigma1: %4.2f +- %4.2f",GetSigma1(),GetSigma1Err()) 541 543 << endl; 542 544 *fLog << all 543 545 << " Pedestal Fit OK? :" << IsPedestalFitOK() 544 << Form(" %s%4.2f%s%4.2f"," Lambda (Check): " ,GetLambdaCheck(),"+-",GetLambdaCheckErr())546 << Form(" Lambda (Check): %4.2f +- %4.2f",GetLambdaCheck(),GetLambdaCheckErr()) 545 547 << endl; 546 548 *fLog << all 547 549 << " Flux available? :" << IsFluxInsidePlexiglassAvailable() 548 << Form(" %s%4.2f%s%4.2f"," Flux: " ,GetFluxInsidePlexiglass(),"+-",GetFluxInsidePlexiglassErr())550 << Form(" Flux: %4.2f +- %4.2f",GetFluxInsidePlexiglass(),GetFluxInsidePlexiglassErr()) 549 551 << endl; 550 552 } 551 552 553 554 555 556 557 558 559 -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc
r5682 r5701 810 810 Form("%s%2.1f%s","Signal Error bigger than ",fChargeRelErrLimit," times Mean Signal: ")); 811 811 PrintUncalibrated(MBadPixelsPix::kLoGainSaturation, 812 " Pixels with Low Gain Saturation:");812 "Low Gain Saturation: "); 813 813 PrintUncalibrated(MBadPixelsPix::kMeanTimeInFirstBin, 814 814 Form("%s%2.1f%s","Mean Abs. Arr. Time in First ",1.," Bin(s): ")); … … 824 824 825 825 PrintUncalibrated(MBadPixelsPix::kChargeSigmaNotValid, 826 "Signal Sigma smaller than Pedestal RMS: 826 "Signal Sigma smaller than Pedestal RMS: "); 827 827 PrintUncalibrated(MBadPixelsPix::kHiGainOscillating, 828 " Pixels with changing Hi Gain signal over time:");828 "Changing Hi Gain signal over time: "); 829 829 PrintUncalibrated(MBadPixelsPix::kLoGainOscillating, 830 " Pixels with changing Lo Gain signal over time:");830 "Changing Lo Gain signal over time: "); 831 831 PrintUncalibrated(MBadPixelsPix::kHiGainNotFitted, 832 " Pixels with unsuccesful Gauss fit to the Hi Gain:");832 "Unsuccesful Gauss fit to the Hi Gain: "); 833 833 PrintUncalibrated(MBadPixelsPix::kLoGainNotFitted, 834 " Pixels with unsuccesful Gauss fit to the Lo Gain:");834 "Unsuccesful Gauss fit to the Lo Gain: "); 835 835 PrintUncalibrated(MBadPixelsPix::kDeviatingNumPhes, 836 " Pixels with deviating number of phes:");836 "Deviating number of phes: "); 837 837 PrintUncalibrated(MBadPixelsPix::kDeviatingFFactor, 838 " Pixels with deviating F-Factor:");838 "Deviating F-Factor: "); 839 839 840 840 if (!fOutputFile.IsNull()) … … 1218 1218 } 1219 1219 1220 *fLog << inf << GetDescriptor() << ": Mean number of photo-electrons " 1221 << "with area idx " << i << ": " 1220 *fLog << inf << GetDescriptor() << ": Mean number of phes with area idx " << i << ": " 1222 1221 << Form("%7.2f+-%6.2f",mean,sigma) << endl; 1223 1222 … … 1265 1264 if ( nphe < lowlim[aidx] || nphe > upplim[aidx] ) 1266 1265 { 1267 *fLog << warn << GetDescriptor() << ": Number of phes: "1266 *fLog << warn << "Number of phes: " 1268 1267 << Form("%7.2f out of %3.1f sigma limit: ",nphe,fPheErrLimit) 1269 1268 << Form("[%7.2f,%7.2f] pixel%4i",lowlim[aidx],upplim[aidx],i) << endl; … … 1321 1320 1322 1321 *fLog << inf << GetDescriptor() 1323 << ": Average total number phes in area idx " << aidx << ": "1322 << ": Average total phes for area idx " << aidx << ": " 1324 1323 << Form("%7.2f +- 6.2f",areaphes[aidx],TMath::Sqrt(areaweights[aidx])) << endl; 1325 1324 … … 1356 1355 { 1357 1356 *fLog << warn << GetDescriptor() 1358 <<": Mean number phes perarea for sector " << sector << " could not be calculated: "1357 <<": Mean number phes/area for sector " << sector << " could not be calculated: " 1359 1358 << " Mean: " << sectorphes[sector] 1360 1359 << " Variance: " << sectorweights[sector] << endl; … … 1364 1363 1365 1364 *fLog << inf << GetDescriptor() 1366 << ": Av erage number phes per area in sector " << sector << ": "1367 << Form("%5.3f+-%4.3f [phe/mm^2]",sectorphes[sector],TMath::Sqrt(sectorweights[sector]))1365 << ": Avg number phes/mm^2 in sector " << sector << ": " 1366 << Form("%5.3f+-%4.3f",sectorphes[sector],TMath::Sqrt(sectorweights[sector])) 1368 1367 << endl; 1369 1368 … … 1758 1757 if ( ffactor < lowlim[aidx] || ffactor > upplim[aidx] ) 1759 1758 { 1760 *fLog << warn << GetDescriptor() << ":Overall F-Factor "1759 *fLog << warn << "Overall F-Factor " 1761 1760 << Form("%5.2f",ffactor) << " out of range [" 1762 << Form("%5.2f,%5.2f",lowlim[aidx],upplim[aidx]) << "] pixel " << i << endl;1761 << Form("%5.2f,%5.2f",lowlim[aidx],upplim[aidx]) << "] Pixel " << i << endl; 1763 1762 1764 1763 bad.SetUncalibrated( MBadPixelsPix::kDeviatingFFactor );
Note:
See TracChangeset
for help on using the changeset viewer.