Changeset 5194


Ignore:
Timestamp:
10/06/04 11:44:00 (20 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r5193 r5194  
    2222 2004/10/05: Markus Gaug
    2323   
     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
    2430  * mcalib/MCalibrationChargeCalc.cc
    2531    - fixed small bug for the case that multiple colours are found, but
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

    r5192 r5194  
    910910                      prms    * fSqrtLoGainSamples,
    911911                      prmserr * fSqrtLoGainSamples);
    912       cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples, aidx);
     912      cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples);
    913913    }
    914914  else
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

    r5130 r5194  
    8181using namespace std;
    8282
    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;
     83const Float_t MCalibrationChargePix::gkElectronicPedRms    = 0.8;
     84const Float_t MCalibrationChargePix::gkElectronicPedRmsErr = 0.1;
     85const Float_t MCalibrationChargePix::gkFFactor             = 1.15;
     86const Float_t MCalibrationChargePix::gkFFactorErr          = 0.02;
    8887
    8988const Float_t MCalibrationChargePix::fgConversionHiLo           = 10.;
     
    244243//
    245244// 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//
     247const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const
     248{
     249  return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms;
    259250}
    260251 
     
    1006997//   Var (Low Gain Ped RMS Square) = Var(LoGain NSB square) + Var(elec.ped. square)
    1007998//
    1008 void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples, const Int_t aidx)
     999void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples)
    10091000{
    10101001
     
    10151006    return;
    10161007
    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() )
    10201010                               * elecPedRms * elecPedRms;
    10211011 
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

    r4776 r5194  
    1010private:
    1111
    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)
    1716
    1817  static const Float_t fgConversionHiLo;         //! Default fConversionHiLo          (now set to: 10.)
     
    4544  enum  { kFFactorMethodValid   };          // Possible bits to be set
    4645
    47   const Float_t GetElectronicPedRmsRelVar(const Int_t aidx) const;
     46  const Float_t GetElectronicPedRmsRelVar() const;
    4847  const Float_t GetConversionHiLoRelVar()  const;
    4948  const Float_t GetFFactorRelVar()     const;
     
    108107
    109108  // Calculations
    110   void   CalcLoGainPedestal       ( const Float_t logainsamples, const Int_t aidx );
     109  void   CalcLoGainPedestal       ( const Float_t logainsamples);
    111110  Bool_t CalcReducedSigma  (  );
    112111  Bool_t CalcFFactor       (  );
Note: See TracChangeset for help on using the changeset viewer.