Changeset 3128 for trunk/MagicSoft


Ignore:
Timestamp:
02/13/04 10:50:42 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r3126 r3128  
    44
    55                                                 -*-*- END OF LINE -*-*-
     6 2004/02/13: Markus Gaug
     7
     8   * mcalib/MCalibrationPix.[h,cc]
     9     - added variables gkAverageQE, gkAverageQEErr,
     10       and fAverageQE, fAverageQEErr,
     11     - gkAverageQE initializes to 0.25
     12     - gkAverageQEErr initializes to 0.03
     13     - added Function: SetAverageQE
     14     - call to GetMeanConversionFFactorMethod returns now:
     15       fMeanConversionFFactorMethod/fAverageQE.
     16     - call to GetErrorConversionFFactorMethod returns now:
     17       errors of fMeanConversionFFactorMethod and fAverageQEErr added
     18       quadratically     
     19
     20
    621 2004/02/13: Thomas Bretz
    722

  • trunk/MagicSoft/Mars/mcalib/MCalibrationPix.cc

    r3115 r3128  
    6262const Float_t MCalibrationPix::gkConvFFactorRelErrorLimit = 0.1;
    6363
     64const Float_t MCalibrationPix::gkAverageQE                = 0.25;     
     65const Float_t MCalibrationPix::gkAverageQEErr             = 0.03; 
    6466// --------------------------------------------------------------------------
    6567//
    … …  
    8789
    8890  Clear();
     91
     92  SetAverageQE();
     93
    8994}
    9095
    … …  
    421426}
    422427
     428//
     429// FIXME: This is a preliminary solution, the qe shall be
     430// calibrated itself!
     431//
    423432Float_t MCalibrationPix::GetMeanConversionFFactorMethod()
    424433{
    … …  
    426435  if (!fFactorCalculated)
    427436    CalcFFactorMethod();
    428  
    429   return fMeanConversionFFactorMethod;
     437
     438  return fMeanConversionFFactorMethod/fAverageQE;
    430439 
    431440}
    … …  
    437446    CalcFFactorMethod();
    438447 
    439   return fErrorConversionFFactorMethod;
     448  Float_t var = fErrorConversionFFactorMethod*fErrorConversionFFactorMethod
     449             / (fMeanConversionFFactorMethod * fMeanConversionFFactorMethod);
     450
     451  var += fAverageQEErr * fAverageQEErr
     452       / (fAverageQE * fAverageQE);
     453
     454  if (var > 0)
     455    return -1.;
     456
     457  var = TMath::Sqrt(var);
     458
     459  return var*GetMeanConversionFFactorMethod();
    440460 
    441461}

  • trunk/MagicSoft/Mars/mcalib/MCalibrationPix.h

    r3115 r3128  
    1212private:
    1313
     14  static const Float_t gkAverageQE;           // The average quantum efficieny agreed on for the first analysis
     15  static const Float_t gkAverageQEErr;        // The error of average quantum efficieny
     16 
    1417  static const Float_t gkElectronicPedRms;    // The pure electronic component of the RMS
    1518  static const Float_t gkErrElectronicPedRms; // The error of the pure electronic component of the RMS
    … …  
    2225  static const Float_t gkTimeErrLimit;        // The limit (in units of FADC slices) for acceptance of the fitted time sigma
    2326  static const Float_t gkConvFFactorRelErrorLimit; // The limit (in units of [1]) for acceptance of the rel. error of the conversion factor with the FFactor method
    24  
    25  
     27
    2628  Int_t   fPixId;                     // the pixel Id
    2729
    2830  UInt_t  fFlags;               // Flag for the set Bits
     31 
     32  Float_t fAverageQE;
     33  Float_t fAverageQEErr; 
    2934 
    3035  Float_t fCharge;              // The mean reduced charge after the fit
    … …  
    101106  MHCalibrationPixel *GetHist() const   { return fHist;     }
    102107
    103   // Charges
    104   Float_t GetCharge()              const { return fCharge;          }
    105   Float_t GetErrCharge()           const { return fErrCharge;       }
    106   Float_t GetChargeProb()          const { return fChargeProb;      }   
    107   Float_t GetSigmaCharge()         const { return fSigmaCharge;     }
    108   Float_t GetErrSigmaCharge()      const { return fErrSigmaCharge;  }
    109   Float_t GetRSigmaCharge()        const { return fRSigmaCharge;    }
    110   Float_t GetErrRSigmaCharge()     const { return fErrRSigmaCharge; } 
    111 
    112 
    113   Float_t GetAbsTimeMean()         const { return fAbsTimeMean;    }
    114   Float_t GetAbsTimeMeanErr()      const { return fAbsTimeMeanErr; }
    115   Float_t GetAbsTimeRms()          const { return fAbsTimeRms;     }
    116  
    117   // Conversion Factors
    118   Float_t GetConversionHiLo()                 const  { return fConversionHiLo;        }
    119   Float_t GetConversionHiLoError()            const  { return fConversionHiLoError;   }
    120 
    121   Float_t GetMeanConversionBlindPixelMethod()  const { return fMeanConversionBlindPixelMethod  ; }
    122   Float_t GetErrorConversionBlindPixelMethod() const { return fErrorConversionBlindPixelMethod ; }
    123   Float_t GetSigmaConversionBlindPixelMethod() const { return fSigmaConversionBlindPixelMethod ; }
    124 
    125   Float_t GetMeanConversionFFactorMethod();
    126   Float_t GetErrorConversionFFactorMethod();
    127   Float_t GetSigmaConversionFFactorMethod();
    128 
    129   Float_t GetMeanConversionPINDiodeMethod()    const { return fMeanConversionPINDiodeMethod ;  }
    130   Float_t GetErrorConversionPINDiodeMethod()   const { return fErrorConversionPINDiodeMethod ; }
    131   Float_t GetSigmaConversionPINDiodeMethod()   const { return fSigmaConversionPINDiodeMethod ; }
    132 
    133   Float_t GetMeanConversionCombinedMethod()    const { return fMeanConversionCombinedMethod ;  }
    134   Float_t GetErrorConversionCombinedMethod()   const { return fErrorConversionCombinedMethod ; }
    135   Float_t GetSigmaConversionCombinedMethod()   const { return fSigmaConversionCombinedMethod ; }
    136 
    137   Float_t GetPheFFactorMethod();   
    138   Float_t GetPheFFactorMethodError();
    139 
    140   Int_t   GetPixId()                           const  { return fPixId;   }
    141 
    142   Float_t GetPed()                             const { return fPed;    }
    143   Float_t GetPedRms()                          const { return fPedRms; }
    144 
    145   Float_t GetTotalFFactorFFactorMethod();
    146   Float_t GetTotalFFactorErrorFFactorMethod();
    147  
    148   Float_t GetTotalFFactorBlindPixelMethod();
    149   Float_t GetTotalFFactorErrorBlindPixelMethod();
    150  
    151   Float_t GetTotalFFactorPINDiodeMethod();
    152   Float_t GetTotalFFactorErrorPINDiodeMethod();
    153 
    154   Float_t GetTotalFFactorCombinedMethod();
    155   Float_t GetTotalFFactorErrorCombinedMethod();
    156  
    157   Bool_t IsExcluded()              const;
    158   Bool_t IsExcludeQualityCheck()   const;
    159   Bool_t IsHiGainSaturation()      const;
    160   Bool_t IsChargeValid()        const;
    161   Bool_t IsFitted()                const;
    162   Bool_t IsOscillating();
    163   Bool_t IsBlindPixelMethodValid() const;
    164   Bool_t IsPINDiodeMethodValid()   const;
    165   Bool_t IsFFactorMethodValid();
    166   Bool_t IsCombinedMethodValid();
    167 
    168108  // Setter
    169   void SetPedestal(Float_t ped, Float_t pedrms, Float_t higainsamp, Float_t logainsamp);
    170   void SetConversionHiLo(Float_t c)      { fConversionHiLo      = c;    }
    171   void SetConversionHiLoError(Float_t e)  { fConversionHiLoError = e;    }
    172 
     109  void SetPedestal(const Float_t ped, const Float_t pedrms,
     110                   const Float_t higainsamp, const Float_t logainsamp);
     111  void SetConversionHiLo(const Float_t c)      { fConversionHiLo      = c;    }
     112  void SetConversionHiLoError(const Float_t e)  { fConversionHiLoError = e;    }
     113  void SetAverageQE(const Float_t qe=gkAverageQE, const Float_t err=gkAverageQEErr)
     114    { fAverageQE    = qe;
     115      fAverageQEErr = err;  }
     116 
    173117  // Setters for MC
    174118  void SetConversionFFactorMethod(Float_t c, Float_t err, Float_t sig);
    … …  
    190134  void SetAbsTimeBordersLoGain(Byte_t f, Byte_t l);
    191135 
     136  // Charges
     137  Float_t GetCharge()              const { return fCharge;          }
     138  Float_t GetErrCharge()           const { return fErrCharge;       }
     139  Float_t GetChargeProb()          const { return fChargeProb;      }   
     140  Float_t GetSigmaCharge()         const { return fSigmaCharge;     }
     141  Float_t GetErrSigmaCharge()      const { return fErrSigmaCharge;  }
     142  Float_t GetRSigmaCharge()        const { return fRSigmaCharge;    }
     143  Float_t GetErrRSigmaCharge()     const { return fErrRSigmaCharge; } 
     144
     145
     146  Float_t GetAbsTimeMean()         const { return fAbsTimeMean;    }
     147  Float_t GetAbsTimeMeanErr()      const { return fAbsTimeMeanErr; }
     148  Float_t GetAbsTimeRms()          const { return fAbsTimeRms;     }
     149 
     150  // Conversion Factors
     151  Float_t GetConversionHiLo()                 const  { return fConversionHiLo;        }
     152  Float_t GetConversionHiLoError()            const  { return fConversionHiLoError;   }
     153
     154  Float_t GetMeanConversionBlindPixelMethod()  const { return fMeanConversionBlindPixelMethod  ; }
     155  Float_t GetErrorConversionBlindPixelMethod() const { return fErrorConversionBlindPixelMethod ; }
     156  Float_t GetSigmaConversionBlindPixelMethod() const { return fSigmaConversionBlindPixelMethod ; }
     157
     158  Float_t GetMeanConversionFFactorMethod();
     159  Float_t GetErrorConversionFFactorMethod();
     160  Float_t GetSigmaConversionFFactorMethod();
     161
     162  Float_t GetMeanConversionPINDiodeMethod()    const { return fMeanConversionPINDiodeMethod ;  }
     163  Float_t GetErrorConversionPINDiodeMethod()   const { return fErrorConversionPINDiodeMethod ; }
     164  Float_t GetSigmaConversionPINDiodeMethod()   const { return fSigmaConversionPINDiodeMethod ; }
     165
     166  Float_t GetMeanConversionCombinedMethod()    const { return fMeanConversionCombinedMethod ;  }
     167  Float_t GetErrorConversionCombinedMethod()   const { return fErrorConversionCombinedMethod ; }
     168  Float_t GetSigmaConversionCombinedMethod()   const { return fSigmaConversionCombinedMethod ; }
     169
     170  Float_t GetPheFFactorMethod();   
     171  Float_t GetPheFFactorMethodError();
     172
     173  Int_t   GetPixId()                           const  { return fPixId;   }
     174
     175  Float_t GetPed()                             const { return fPed;    }
     176  Float_t GetPedRms()                          const { return fPedRms; }
     177
     178  Float_t GetTotalFFactorFFactorMethod();
     179  Float_t GetTotalFFactorErrorFFactorMethod();
     180 
     181  Float_t GetTotalFFactorBlindPixelMethod();
     182  Float_t GetTotalFFactorErrorBlindPixelMethod();
     183 
     184  Float_t GetTotalFFactorPINDiodeMethod();
     185  Float_t GetTotalFFactorErrorPINDiodeMethod();
     186
     187  Float_t GetTotalFFactorCombinedMethod();
     188  Float_t GetTotalFFactorErrorCombinedMethod();
     189 
     190  Bool_t IsExcluded()              const;
     191  Bool_t IsExcludeQualityCheck()   const;
     192  Bool_t IsHiGainSaturation()      const;
     193  Bool_t IsChargeValid()        const;
     194  Bool_t IsFitted()                const;
     195  Bool_t IsOscillating();
     196  Bool_t IsBlindPixelMethodValid() const;
     197  Bool_t IsPINDiodeMethodValid()   const;
     198  Bool_t IsFFactorMethodValid();
     199  Bool_t IsCombinedMethodValid();
     200
    192201  // Fill histos
    193202  Bool_t FillChargeHiGain(Float_t q)         const { return fHist->FillChargeHiGain(q); }
Note: See TracChangeset for help on using the changeset viewer.