Ignore:
Timestamp:
06/03/08 12:29:27 (17 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mcalib
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

    r8147 r8917  
    1818!   Author(s): Markus Gaug   02/2004 <mailto:markus@ifae.es>
    1919!
    20 !   Copyright: MAGIC Software Development, 2000-2004
     20!   Copyright: MAGIC Software Development, 2000-2008
    2121!
    2222\* ======================================================================== */
     
    6565//  FADC slices extraction results, but can be retrieved converted to High-Gain amplifications
    6666//  by calls to: GetConvertedMean() or GetConvertedSigma()
     67//
     68// Inline Functions:
     69// -----------------
     70//
     71// GetPedRms(): Get the pedestals RMS: Test bit kHiGainSaturation:
     72//  If yes, return square root of fLoGainPedRmsSquare (if greater than 0,
     73//  otherwise -1.), If no,  return fPedRms
     74//
     75// GetConvertedMean(): Get the Low Gain Mean Charge converted to High Gain
     76//  amplification: Returns fLoGainMean multiplied with fConversionHiLo if
     77//  IsHiGainSaturation(), else return fHiGainMean
     78//
     79// GetConvertedSigma(): Get the Low Gain Sigma converted to High Gain
     80//  amplification: Returns fLoGainSigma multiplied with fConversionHiLo
     81//  if IsHiGainSaturation() else return fHiGainSigma
     82//
     83// GetConvertedRSigmaSquare(): Get the reduced Sigma Square:
     84//  If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1.
     85//  Test bit kHiGainSaturation:
     86//  If yes, return fRSigmaSquare, multiplied with fConversionHiLo^2,
     87//  If no , return fRSigmaSquare
     88//
     89// GetPheFFactorMethodRelVar(): Get the relative variance on the number of
     90//  photo-electrons (F-Factor Method): If fPheFFactorMethodVar is smaller
     91//  than 0 (i.e. has not yet been set), return -1. If fPheFFactorMethod
     92//  is 0, return -1. Else returns fPheFFactorMethodVar / fPheFFactorMethod^2
     93//
     94// GetMeanConvFADC2PheErr(): Get the error on the mean conversion factor
     95//  (FFactor  Method): If fMeanConvFADC2PheVar is smaller than 0 (i.e. has
     96//  not yet been set), return -1. Else returns the square root of
     97//  fMeanConvFADC2PheVar
     98//
    6799//
    68100// See also: MCalibrationChargeCam, MCalibrationChargeCalc,
     
    88120/////////////////////////////////////////////////////////////////////////////
    89121#include "MCalibrationChargePix.h"
     122
     123#include <TMath.h>
    90124
    91125#include "MLog.h"
     
    264298}
    265299
    266 
    267 // Inline Functions:
    268 // -----------------
    269 //
    270 // GetConversionHiLoErr:
    271 //   Get the conversion Error Hi-Gain to Low-Gain:
    272 //   If fConversionHiLoVar is smaller than 0 (i.e. has not yet
    273 //   been set), return -1.
     300// --------------------------------------------------------------------------
     301//
     302//  Get the conversion Error Hi-Gain to Low-Gain:
     303//  If fConversionHiLoVar is smaller than 0 (i.e. has not yet
     304//  been set), return -1.
    274305// 
    275 // GetPedRms(): Get the pedestals RMS: Test bit kHiGainSaturation:
    276 //  If yes, return square root of fLoGainPedRmsSquare (if greater than 0,
    277 //  otherwise -1.), If no,  return fPedRms
    278 //
    279 // GetConvertedMean(): Get the Low Gain Mean Charge converted to High Gain
    280 //  amplification: Returns fLoGainMean multiplied with fConversionHiLo if
    281 //  IsHiGainSaturation(), else return fHiGainMean
    282 //
    283 // GetConvertedMeanErr(): Get the Error of the converted Low Gain Mean:
    284 //  Returns -1 if the variable fLoGainMean or fLoGainMeanVar are smaller than 0.
    285 //  Returns the square root of the quadratic sum of the relative variances of
    286 //  the fLoGainMean and fConversionHiLo, mulitplied with GetConvertedMean()
    287 //  in case of HiGain Saturation,
    288 //  else return GetMeanErr()
    289 //
    290 // GetConvertedSigma(): Get the Low Gain Sigma converted to High Gain
    291 //  amplification: Returns fLoGainSigma multiplied with fConversionHiLo
    292 //  if IsHiGainSaturation() else return fHiGainSigma
    293 //
    294 // GetConvertedSigmaErr(): Get the Error of the converted Sigma:
    295 //  Returns -1 if the variable fLoGainSigma or fLoGainSigmaVar are smaller than 0.
    296 //  if IsHiGainSaturatio()
    297 //  returns the square root of the quadratic sum of the relative variances of
    298 //  the fLoGainSigma and fConversionHiLo, mulitplied with GetConvertedSigma()
    299 //  else returns GetSigmaErr()
    300 //
    301 // GetConvertedRSigma(): Get the converted reduced Sigma:
    302 //  If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1.
    303 //  Test bit kHiGainSaturation:
    304 //  If yes, return square root of fRSigmaSquare, multiplied with fConversionHiLo,
    305 //  If no , return square root of fRSigmaSquare
    306 //
    307 // GetRSigma(): Get the reduced Sigma:
    308 //  If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1.
    309 //
    310 // GetConvertedRSigmaSquare(): Get the reduced Sigma Square:
    311 //  If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1.
    312 //  Test bit kHiGainSaturation:
    313 //  If yes, return fRSigmaSquare, multiplied with fConversionHiLo^2,
    314 //  If no , return fRSigmaSquare
     306Float_t MCalibrationChargePix::GetConversionHiLoErr() const
     307{
     308    return fConversionHiLoVar<0 ? -1 : TMath::Sqrt(fConversionHiLoVar);
     309}
     310
     311// --------------------------------------------------------------------------
     312//
     313// Get the Error of the converted Low Gain Mean: Returns -1 if the
     314// variable fLoGainMean or fLoGainMeanVar are smaller than 0. Returns the
     315// square root of the quadratic sum of the relative variances of the
     316// fLoGainMean and fConversionHiLo, mulitplied with GetConvertedMean() in
     317// case of HiGain Saturation, else return GetMeanErr()
     318//
     319Float_t MCalibrationChargePix::GetConvertedMeanErr() const
     320{
     321    if (!IsHiGainSaturation())
     322        return GetMeanErr();
     323
     324    const Float_t logainrelvar = GetLoGainMeanRelVar();
     325    return logainrelvar<0 ? -1 : TMath::Sqrt(logainrelvar + GetConversionHiLoRelVar()) * GetConvertedMean();
     326}
     327
     328// --------------------------------------------------------------------------
     329//
     330// Get the Error of the converted Sigma: Returns -1 if the variable
     331// fLoGainSigma or fLoGainSigmaVar are smaller than 0. if
     332// IsHiGainSaturatio() returns the square root of the quadratic sum of the
     333// relative variances of the fLoGainSigma and fConversionHiLo, mulitplied
     334// with GetConvertedSigma() else returns GetSigmaErr()
     335//
     336Float_t MCalibrationChargePix::GetConvertedSigmaErr() const
     337{
     338    if (!IsHiGainSaturation())
     339        return GetSigmaErr();
     340
     341    if (fLoGainSigmaVar<0 || fLoGainSigma<0)
     342        return -1.;
     343
     344    const Float_t sigmaRelVar = fLoGainSigmaVar/(fLoGainSigma*fLoGainSigma);
     345    return TMath::Sqrt(sigmaRelVar+GetConversionHiLoRelVar()) * GetConvertedSigma();
     346}
     347
     348// --------------------------------------------------------------------------
     349//
     350// Get the converted reduced Sigma: If fRSigmaSquare is smaller than 0
     351// (i.e. has not yet been set), return -1. Test bit kHiGainSaturation: If
     352// yes, return square root of fRSigmaSquare, multiplied with
     353// fConversionHiLo, If no , return square root of fRSigmaSquare
     354//
     355Float_t MCalibrationChargePix::GetConvertedRSigma() const
     356{
     357    if (fRSigmaSquare < 0)
     358        return -1;
     359
     360    const Float_t rsigma = TMath::Sqrt(fRSigmaSquare);
     361    return IsHiGainSaturation() ? rsigma*fConversionHiLo : rsigma ;
     362}
     363
     364// --------------------------------------------------------------------------
     365//
     366Float_t MCalibrationChargePix::GetMeanConvFADC2PheErr() const
     367{
     368    return fMeanConvFADC2PheStatVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheStatVar);
     369}
     370
     371// --------------------------------------------------------------------------
     372//
     373Float_t MCalibrationChargePix::GetMeanConvFADC2PheSystErr() const
     374{
     375    return fMeanConvFADC2PheSystVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheSystVar);
     376}
     377
     378// --------------------------------------------------------------------------
     379//
     380Float_t MCalibrationChargePix::GetMeanConvFADC2PheTotErr() const
     381{
     382    if (fMeanConvFADC2PheSystVar<0 || fMeanConvFADC2PheStatVar<0)
     383        return -1.;
     384
     385    return TMath::Sqrt(fMeanConvFADC2PheSystVar+fMeanConvFADC2PheStatVar);
     386}
     387
     388// --------------------------------------------------------------------------
     389//
     390//  Get the error on the mean total F-Factor of the signal readout
     391//  (F-Factor Method): If fMeanFFactorFADC2PhotVar is smaller than 0
     392//  (i.e. has not yet been set), return -1. Else returns the square
     393//  root of fMeanFFactorFADC2PhotVar
     394//
     395Float_t MCalibrationChargePix::GetMeanFFactorFADC2PhotErr() const
     396{
     397    return fMeanFFactorFADC2PhotVar<0 ? -1. : TMath::Sqrt(fMeanFFactorFADC2PhotVar);
     398}
     399
     400// --------------------------------------------------------------------------
     401//
     402// Get the relative variance of the conversion factor between higain and logain:
     403// - If fConversionHiLo    is 0,              return -1.
     404// - If fConversionHiLoVar is smaller than 0, return -1.
     405// - Else returns: fConversionHiLoVar / fConversionHiLo^2
     406//
     407const Float_t MCalibrationChargePix::GetConversionHiLoRelVar() const
     408{
     409
     410  if (fConversionHiLoVar < 0.)
     411    return -1.;
     412
     413  if (fConversionHiLo == 0.)
     414    return -1.;
     415
     416  return fConversionHiLoVar / (fConversionHiLo * fConversionHiLo);
     417}
     418
     419// --------------------------------------------------------------------------
     420//
     421// Get the relative variance of the electronics pedestal RMS
     422// - returns rel. variance of gkElectronicPedRms
     423//
     424const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const
     425{
     426  return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms;
     427}
     428 
     429
     430// --------------------------------------------------------------------------
     431//
     432// Get the relative variance of the conversion factor between higain and logain:
     433// - If gkFFactor    is 0,              return -1.
     434// - If gkFFactorErr is smaller than 0, return -1.
     435// - Else returns: gkFFactorErr^2 / gkFFactor*^2
     436//
     437const Float_t MCalibrationChargePix::GetFFactorRelVar() const
     438{
     439
     440  if (gkFFactorErr < 0.)
     441    return -1.;
     442
     443  if (gkFFactor == 0.)
     444    return -1.;
     445
     446  return gkFFactorErr * gkFFactorErr / (gkFFactor * gkFFactor);
     447}
     448
     449// --------------------------------------------------------------------------
     450//
     451Float_t MCalibrationChargePix::GetPedErr() const
     452{
     453    return TMath::Sqrt(fPedVar);
     454}
     455
     456// --------------------------------------------------------------------------
     457//
     458Float_t MCalibrationChargePix::GetPedRms() const
     459{
     460    if (!IsHiGainSaturation())
     461        return fPedRms;
     462
     463    return fLoGainPedRmsSquare<0 ? -1 : TMath::Sqrt(fLoGainPedRmsSquare);
     464}
     465
     466// --------------------------------------------------------------------------
     467//
     468// Get the Error of the pedestals RMS:
     469// - Test bit kHiGainSaturation:
     470//   If yes, return square root of (0.25*fLoGainPedRmsSquareVar/ fLoGainPedRmsSquare) (if greater than 0, otherwise -1.)
     471//   If no , return square root of (fPedVar) (if greater than 0, otherwise -1.), divided by 2.
     472//
     473Float_t MCalibrationChargePix::GetPedRmsErr()  const
     474{
     475  if (IsHiGainSaturation())
     476      return fLoGainPedRmsSquareVar < 0. ? -1. :TMath::Sqrt(0.25*fLoGainPedRmsSquareVar/fLoGainPedRmsSquare);
     477  else
     478      return fPedVar < 0. ? -1. : TMath::Sqrt(fPedVar)/2.;
     479}
     480
     481// --------------------------------------------------------------------------
    315482//
    316483// GetPheFFactorMethodErr(): Get the error on the number of photo-electrons
     
    319486//  fPheFFactorMethodVar
    320487//
    321 // GetMeanFFactorFADC2PhotErr(): Get the error on the mean total F-Factor
    322 //  of the signal readout (F-Factor Method): If fMeanFFactorFADC2PhotVar
    323 //  is smaller than 0 (i.e. has not yet been set), return -1. Else returns
    324 //  the square root of fMeanFFactorFADC2PhotVar
    325 //
    326 // GetPheFFactorMethodRelVar(): Get the relative variance on the number of
    327 //  photo-electrons (F-Factor Method): If fPheFFactorMethodVar is smaller
    328 //  than 0 (i.e. has not yet been set), return -1. If fPheFFactorMethod
    329 //  is 0, return -1. Else returns fPheFFactorMethodVar / fPheFFactorMethod^2
    330 //
    331 // GetMeanConvFADC2PheErr(): Get the error on the mean conversion factor
    332 //  (FFactor  Method): If fMeanConvFADC2PheVar is smaller than 0 (i.e. has
    333 //  not yet been set), return -1. Else returns the square root of
    334 //  fMeanConvFADC2PheVar
    335 
    336 // --------------------------------------------------------------------------
    337 //
    338 // Get the relative variance of the conversion factor between higain and logain:
    339 // - If fConversionHiLo    is 0,              return -1.
    340 // - If fConversionHiLoVar is smaller than 0, return -1.
    341 // - Else returns: fConversionHiLoVar / fConversionHiLo^2
    342 //
    343 const Float_t MCalibrationChargePix::GetConversionHiLoRelVar() const
    344 {
    345 
    346   if (fConversionHiLoVar < 0.)
    347     return -1.;
    348 
    349   if (fConversionHiLo == 0.)
    350     return -1.;
    351 
    352   return fConversionHiLoVar / (fConversionHiLo * fConversionHiLo);
    353 }
    354 
    355 // --------------------------------------------------------------------------
    356 //
    357 // Get the relative variance of the electronics pedestal RMS
    358 // - returns rel. variance of gkElectronicPedRms
    359 //
    360 const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const
    361 {
    362   return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms;
    363 }
    364  
    365 
    366 // --------------------------------------------------------------------------
    367 //
    368 // Get the relative variance of the conversion factor between higain and logain:
    369 // - If gkFFactor    is 0,              return -1.
    370 // - If gkFFactorErr is smaller than 0, return -1.
    371 // - Else returns: gkFFactorErr^2 / gkFFactor*^2
    372 //
    373 const Float_t MCalibrationChargePix::GetFFactorRelVar() const
    374 {
    375 
    376   if (gkFFactorErr < 0.)
    377     return -1.;
    378 
    379   if (gkFFactor == 0.)
    380     return -1.;
    381 
    382   return gkFFactorErr * gkFFactorErr / (gkFFactor * gkFFactor);
    383 }
    384  
    385 
    386 
    387 // --------------------------------------------------------------------------
    388 //
    389 // Get the Error of the pedestals RMS:
    390 // - Test bit kHiGainSaturation:
    391 //   If yes, return square root of (0.25*fLoGainPedRmsSquareVar/ fLoGainPedRmsSquare) (if greater than 0, otherwise -1.)
    392 //   If no , return square root of (fPedVar) (if greater than 0, otherwise -1.), divided by 2.
    393 //
    394 Float_t  MCalibrationChargePix::GetPedRmsErr()  const
    395 {
    396   if (IsHiGainSaturation())
    397     if (fLoGainPedRmsSquareVar < 0.)
    398       return -1.;
    399     else
    400       return TMath::Sqrt(0.25*fLoGainPedRmsSquareVar/fLoGainPedRmsSquare);
    401   else
    402     if (fPedVar < 0.)
    403       return -1.;
    404     else
    405       return TMath::Sqrt(fPedVar)/2.;
    406 }
    407 
     488Float_t MCalibrationChargePix::GetPheFFactorMethodErr() const
     489{
     490    return fPheFFactorMethodStatVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodStatVar);
     491}
     492
     493// --------------------------------------------------------------------------
     494//
     495Float_t MCalibrationChargePix::GetPheFFactorMethodSystErr() const
     496{
     497    return fPheFFactorMethodSystVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodSystVar);
     498}
     499
     500// --------------------------------------------------------------------------
     501//
     502Float_t MCalibrationChargePix::GetPheFFactorMethodTotErr() const
     503{
     504    if (fPheFFactorMethodStatVar<0 || fPheFFactorMethodSystVar<0)
     505        return -1.;
     506
     507    return TMath::Sqrt(fPheFFactorMethodStatVar+fPheFFactorMethodSystVar);
     508}
    408509
    409510// --------------------------------------------------------------------------
     
    439540// --------------------------------------------------------------------------
    440541//
     542// Get the reduced Sigma:
     543// If fRSigmaSquare is smaller than 0 (i.e. has not yet been set),
     544// return -1.
     545//
     546Float_t MCalibrationChargePix::GetRSigma() const
     547{
     548    return fRSigmaSquare<0 ? -1 : TMath::Sqrt(fRSigmaSquare);
     549}
     550
     551// --------------------------------------------------------------------------
     552//
    441553// Get the error of the reduced Sigma:
    442554// - If fRSigmaSquareVar is smaller than 0 (i.e. has not yet been set), return -1.
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

    r7189 r8917  
    9696  Float_t GetAbsTimeMean        () const { return fAbsTimeMean;    }
    9797  Float_t GetAbsTimeRms         () const { return fAbsTimeRms;     }
     98
    9899  Float_t GetConversionHiLo     () const { return fConversionHiLo; }
    99   Float_t GetConversionHiLoErr  () const { return fConversionHiLoVar<0 ? -1 : TMath::Sqrt(fConversionHiLoVar);  }
     100  Float_t GetConversionHiLoErr  () const;
    100101  Float_t GetConversionHiLoSigma() const { return fConversionHiLoSigma; }
     102
    101103  Float_t GetConvertedMean      () const { return IsHiGainSaturation() ? fLoGainMean * fConversionHiLo : fHiGainMean; }
    102   Float_t GetConvertedMeanErr   () const
    103   {
    104       if (!IsHiGainSaturation())
    105           return GetMeanErr();
    106       const Float_t logainrelvar = GetLoGainMeanRelVar();
    107       return logainrelvar<0 ? -1 : TMath::Sqrt(logainrelvar + GetConversionHiLoRelVar()) * GetConvertedMean();
    108   }
     104  Float_t GetConvertedMeanErr   () const;
     105
    109106  Float_t GetConvertedSigma() const { return IsHiGainSaturation() ? fLoGainSigma * fConversionHiLo : fHiGainSigma; }
    110   Float_t GetConvertedSigmaErr() const
    111   {
    112     if (!IsHiGainSaturation())
    113         return GetSigmaErr();
     107  Float_t GetConvertedSigmaErr() const;
    114108
    115     if (fLoGainSigmaVar<0 || fLoGainSigma<0)
    116         return -1.;
    117 
    118     const Float_t sigmaRelVar = fLoGainSigmaVar/(fLoGainSigma*fLoGainSigma);
    119     return TMath::Sqrt(sigmaRelVar+GetConversionHiLoRelVar()) * GetConvertedSigma();
    120   }
    121   Float_t GetConvertedRSigma() const
    122   {
    123       if (fRSigmaSquare < 0)
    124           return -1;
    125       const Float_t rsigma = TMath::Sqrt(fRSigmaSquare);
    126       return IsHiGainSaturation() ? rsigma*fConversionHiLo : rsigma ;
    127   }
     109  Float_t GetConvertedRSigma() const;
    128110  Float_t GetConvertedRSigmaErr() const;
    129111  Float_t GetConvertedRSigmaSquare() const
     
    133115      return IsHiGainSaturation() ? fRSigmaSquare*fConversionHiLo*fConversionHiLo : fRSigmaSquare ;
    134116  }
     117
    135118  Float_t GetMeanConvFADC2Phe() const { return fMeanConvFADC2Phe; }
    136   Float_t GetMeanConvFADC2PheErr  () const { return fMeanConvFADC2PheStatVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheStatVar); }
    137   Float_t GetMeanConvFADC2PheSystErr() const { return fMeanConvFADC2PheSystVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheSystVar); }
    138   Float_t GetMeanConvFADC2PheTotErr() const
    139   {
    140       if (fMeanConvFADC2PheSystVar<0 || fMeanConvFADC2PheStatVar<0)
    141           return -1.;
    142       return TMath::Sqrt(fMeanConvFADC2PheSystVar+fMeanConvFADC2PheStatVar);
    143   }
     119  Float_t GetMeanConvFADC2PheErr() const;
     120  Float_t GetMeanConvFADC2PheSystErr() const;
     121  Float_t GetMeanConvFADC2PheTotErr() const;
     122
    144123  Float_t GetFFactorFADC2Phe        () const { return gkFFactor;   }
    145124  Float_t GetMeanConvFADC2PheVar    () const { return fMeanConvFADC2PheStatVar; }
     
    147126
    148127  Float_t GetMeanFFactorFADC2Phot   () const { return fMeanFFactorFADC2Phot;    }
    149   Float_t GetMeanFFactorFADC2PhotErr() const { return fMeanFFactorFADC2PhotVar<0 ? -1. : TMath::Sqrt(fMeanFFactorFADC2PhotVar); }
     128  Float_t GetMeanFFactorFADC2PhotErr() const;
    150129  Float_t GetMeanFFactorFADC2PhotVar() const { return fMeanFFactorFADC2PhotVar; }
    151   Int_t   GetNumSaturated           () const { return fNumSaturated;            }   
    152   Float_t GetPed                    () const { return fPed;                     }
    153   Float_t GetPedErr                 () const { return TMath::Sqrt(fPedVar);     }
    154   Float_t GetPedRms                 () const
    155   {
    156     if (!IsHiGainSaturation())
    157         return fPedRms;
    158     return fLoGainPedRmsSquare<0 ? -1 : TMath::Sqrt(fLoGainPedRmsSquare);
    159   }
    160   Float_t GetPedRmsErr              () const;
    161   Float_t GetPheFFactorMethod       () const { return fPheFFactorMethod; }
    162   Float_t GetPheFFactorMethodErr    () const { return fPheFFactorMethodStatVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodStatVar); }
    163   Float_t GetPheFFactorMethodSystErr() const { return fPheFFactorMethodSystVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodSystVar); }
    164   Float_t GetPheFFactorMethodTotErr () const
    165   {
    166       if (fPheFFactorMethodStatVar<0 || fPheFFactorMethodSystVar<0)
    167           return -1.;
    168       return TMath::Sqrt(fPheFFactorMethodStatVar+fPheFFactorMethodSystVar);
    169   }
     130
     131  Int_t   GetNumSaturated           () const { return fNumSaturated;            }
     132
     133  Float_t GetPed() const { return fPed; }
     134  Float_t GetPedErr() const;
     135  Float_t GetPedRms() const;
     136  Float_t GetPedRmsErr() const;
     137
     138  Float_t GetPheFFactorMethod() const { return fPheFFactorMethod; }
     139  Float_t GetPheFFactorMethodErr() const;
     140  Float_t GetPheFFactorMethodSystErr() const;
     141  Float_t GetPheFFactorMethodTotErr() const;
     142
    170143  Float_t GetPheFFactorMethodVar    () const { return fPheFFactorMethodStatVar; }
    171144  Float_t GetPheFFactorMethodSystVar() const { return fPheFFactorMethodSystVar; }
    172145  Float_t GetPheFFactorMethodRelVar () const { return fPheFFactorMethodStatVar<=0 ? -1 : fPheFFactorMethodStatVar / (fPheFFactorMethod * fPheFFactorMethod); }
    173146  Float_t GetPheFFactorMethodRelSystVar() const { return fPheFFactorMethodSystVar<=0 ? -1. : fPheFFactorMethodSystVar / (fPheFFactorMethod * fPheFFactorMethod); }
    174   Float_t GetRSigma                 () const { return fRSigmaSquare<0 ? -1 : TMath::Sqrt(fRSigmaSquare); }
    175147
    176   Float_t GetRSigmaErr         () const;
    177   Float_t GetRSigmaRelVar      () const;
    178   Float_t GetRSigmaPerCharge   () const;
     148  Float_t GetRSigma() const;
     149  Float_t GetRSigmaErr() const;
     150  Float_t GetRSigmaRelVar() const;
     151  Float_t GetRSigmaPerCharge() const;
    179152  Float_t GetRSigmaPerChargeErr() const;
    180153
Note: See TracChangeset for help on using the changeset viewer.