Changeset 5194
- Timestamp:
- 10/06/04 11:44:00 (20 years ago)
- Location:
- trunk/MagicSoft/Mars
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/Changelog
r5193 r5194 22 22 2004/10/05: Markus Gaug 23 23 24 * mcalib/MCalibrationChargePix.[h,cc] 25 * mcalib/MCalibrationChargeCalc.cc 26 - after talking with Florian, we reduced the electronic noise for the 27 low-gain pedestal calculation from 1.5 to 0.8, independent from 28 inner or outer pixels. 29 24 30 * mcalib/MCalibrationChargeCalc.cc 25 31 - fixed small bug for the case that multiple colours are found, but -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc
r5192 r5194 910 910 prms * fSqrtLoGainSamples, 911 911 prmserr * fSqrtLoGainSamples); 912 cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples , aidx);912 cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples); 913 913 } 914 914 else -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc
r5130 r5194 81 81 using namespace std; 82 82 83 const Float_t MCalibrationChargePix::gkElectronicPedRmsInner = 1.5; 84 const Float_t MCalibrationChargePix::gkElectronicPedRmsOuter = 1.8; 85 const Float_t MCalibrationChargePix::gkElectronicPedRmsErr = 0.35; 86 const Float_t MCalibrationChargePix::gkFFactor = 1.15; 87 const Float_t MCalibrationChargePix::gkFFactorErr = 0.02; 83 const Float_t MCalibrationChargePix::gkElectronicPedRms = 0.8; 84 const Float_t MCalibrationChargePix::gkElectronicPedRmsErr = 0.1; 85 const Float_t MCalibrationChargePix::gkFFactor = 1.15; 86 const Float_t MCalibrationChargePix::gkFFactorErr = 0.02; 88 87 89 88 const Float_t MCalibrationChargePix::fgConversionHiLo = 10.; … … 244 243 // 245 244 // Get the relative variance of the electronics pedestal RMS 246 // - If aidx is 0, return rel. variance of gkElectronicPedRmsInner 247 // - If aidx is 1, return rel. variance of gkElectronicPedRmsOuter 248 // 249 const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar(const Int_t aidx) const 250 { 251 252 if (aidx == 0) 253 return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRmsInner / gkElectronicPedRmsInner; 254 255 if (aidx == 1) 256 return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRmsOuter / gkElectronicPedRmsOuter; 257 258 return -1.; 245 // - returns rel. variance of gkElectronicPedRms 246 // 247 const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const 248 { 249 return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms; 259 250 } 260 251 … … 1006 997 // Var (Low Gain Ped RMS Square) = Var(LoGain NSB square) + Var(elec.ped. square) 1007 998 // 1008 void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples , const Int_t aidx)999 void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples) 1009 1000 { 1010 1001 … … 1015 1006 return; 1016 1007 1017 const Float_t elecPedRms = (aidx == 0 ? gkElectronicPedRmsInner : gkElectronicPedRmsOuter ) 1018 * TMath::Sqrt(logainsamples) / fConversionHiLo; 1019 const Float_t elecPedRmsVar = ( GetElectronicPedRmsRelVar(aidx) + GetConversionHiLoRelVar() ) 1008 const Float_t elecPedRms = gkElectronicPedRms * TMath::Sqrt(logainsamples) / fConversionHiLo; 1009 const Float_t elecPedRmsVar = ( GetElectronicPedRmsRelVar() + GetConversionHiLoRelVar() ) 1020 1010 * elecPedRms * elecPedRms; 1021 1011 -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h
r4776 r5194 10 10 private: 11 11 12 static const Float_t gkElectronicPedRmsInner; //! Electronic component of ped. RMS (now set to: 1.5) 13 static const Float_t gkElectronicPedRmsOuter; //! Electronic component of ped. RMS (now set to: 1.8) 14 static const Float_t gkElectronicPedRmsErr; //! Error Electr. component ped. RMS (now set to: 0.35) 15 static const Float_t gkFFactor; //! Laboratory F-factor PMTs (now set to: 1.15) 16 static const Float_t gkFFactorErr; //! Laboratory F-factor Error PMTs (now set to: 0.02) 12 static const Float_t gkElectronicPedRms; //! Electronic component of ped. RMS from opt. receiver on (now set to: 0.8) 13 static const Float_t gkElectronicPedRmsErr; //! Error Electr. component ped. RMS (now set to: 0.1) 14 static const Float_t gkFFactor; //! Laboratory F-factor PMTs (now set to: 1.15) 15 static const Float_t gkFFactorErr; //! Laboratory F-factor Error PMTs (now set to: 0.02) 17 16 18 17 static const Float_t fgConversionHiLo; //! Default fConversionHiLo (now set to: 10.) … … 45 44 enum { kFFactorMethodValid }; // Possible bits to be set 46 45 47 const Float_t GetElectronicPedRmsRelVar( const Int_t aidx) const;46 const Float_t GetElectronicPedRmsRelVar() const; 48 47 const Float_t GetConversionHiLoRelVar() const; 49 48 const Float_t GetFFactorRelVar() const; … … 108 107 109 108 // Calculations 110 void CalcLoGainPedestal ( const Float_t logainsamples , const Int_t aidx);109 void CalcLoGainPedestal ( const Float_t logainsamples); 111 110 Bool_t CalcReducedSigma ( ); 112 111 Bool_t CalcFFactor ( );
Note:
See TracChangeset
for help on using the changeset viewer.