Ignore:
Timestamp:
08/01/06 14:13:34 (18 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

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

    r5412 r7829  
    8787//     by one of the three methods.
    8888//
    89 // 11) A call to GetQECascadesBlindPixel(zenith), GetQECascadesFFactor(zenith) or
    90 //     GetQECascadesPINDiode(zenith) returns then the normalization multiplied with an average QE
     89// 11) A call to GetQECascadesBlindPixel(), GetQECascadesFFactor() or
     90//     GetQECascadesPINDiode() returns then the normalization multiplied with an average QE
    9191//     folded into a cascades spectrum. This number should be dependent on zenith angle, but
    9292//     this feature is not yet implemented, instead a fixed number gkDefaultAverageQE is used.
     
    210210const Float_t MCalibrationQEPix::gkLightGuidesEffCT1      = 0.94 ;
    211211const Float_t MCalibrationQEPix::gkLightGuidesEffCT1Err   = 0.03 ;
     212
    212213// --------------------------------------------------------------------------
    213214//
     
    447448// -----------------------------------------------------------------
    448449//
    449 // Return the average Default QE (depending on zenith angle)
    450 // FIXME: The zenith angle dependency is not yet implemented
    451 //
    452 const Float_t MCalibrationQEPix::GetAverageQE( const Float_t zenith ) const
    453 {
    454   //  return gkDefaultAverageQE ;
     450// Return the average Default QE
     451//
     452const Float_t MCalibrationQEPix::GetAverageQE() const
     453{
    455454  return fAverageQE;
    456455}
     
    458457// -----------------------------------------------------------------
    459458//
    460 // Return the relative variance of the average Default QE (depending on zenith angle)
    461 // FIXME: The zenith angle dependency is not yet implemented
    462 //
    463 const Float_t MCalibrationQEPix::GetAverageQERelVar( const Float_t zenith ) const
     459// Return the relative variance of the average Default QE
     460//
     461const Float_t MCalibrationQEPix::GetAverageQERelVar() const
    464462{
    465463  return gkDefaultAverageQEErr * gkDefaultAverageQEErr / (gkDefaultAverageQE * gkDefaultAverageQE );
     
    469467//
    470468// Return the relative variance of the average normalization (Blind Pixel Method)
    471 // FIXME: The zenith angle dependency is not yet implemented
    472469//
    473470const Float_t MCalibrationQEPix::GetAvNormBlindPixelRelVar( ) const
     
    479476//
    480477// Return the relative variance of the average normalization (Combined Method)
    481 // FIXME: The zenith angle dependency is not yet implemented
    482478//
    483479const Float_t MCalibrationQEPix::GetAvNormCombinedRelVar( ) const
     
    498494//
    499495// Return the relative variance of the average normalization (PIN Diode Method)
    500 // FIXME: The zenith angle dependency is not yet implemented
    501496//
    502497const Float_t MCalibrationQEPix::GetAvNormPINDiodeRelVar( ) const
     
    624619// ------------------------------------------------------------------------------
    625620//
    626 // Get the light guides efficiency for Cherenkov spectra,
    627 // depending on the zenith angle of the telescope
     621// Get the light guides efficiency for Cherenkov spectra.
    628622// FIXME: Lacking detailed measurement, these number are not yet available
    629623//        for the individual colours and therefore, only one same number is
    630624//        returned, namely gkLightGuidesEffBlue
    631625//
    632 Float_t MCalibrationQEPix::GetLightGuidesEff( const Float_t zenith )  const
     626Float_t MCalibrationQEPix::GetLightGuidesEff()  const
    633627{
    634628  return gkLightGuidesEffBlue;
     
    638632// ------------------------------------------------------------------------------
    639633//
    640 // Get the relative variance of the light guides efficiency for Cherenkov spectra
    641 // depending on the zenith angle of the telescope
     634// Get the relative variance of the light guides efficiency for Cherenkov spectra.
    642635// FIXME: Lacking detailed measurement, these number are not yet available
    643636//        for the individual colours and therefore, only one same number is
    644637//        returned, namely gkLightGuidesEffBlueErr^2 / gkLightGuidesBlueEff^2
    645638//
    646 Float_t MCalibrationQEPix::GetLightGuidesEffRelVar( const Float_t zenith )  const
     639Float_t MCalibrationQEPix::GetLightGuidesEffRelVar()  const
    647640{
    648641  return gkLightGuidesEffBlueErr  * gkLightGuidesEffBlueErr / gkLightGuidesEffBlue / gkLightGuidesEffBlue;
     
    846839// with the blind pixel method and averaged over the results from the different colours.
    847840//
    848 Float_t MCalibrationQEPix::GetQECascadesBlindPixel( const Float_t zenith  )  const
    849 {
    850   return fAvNormBlindPixel * GetAverageQE ( zenith );
     841Float_t MCalibrationQEPix::GetQECascadesBlindPixel()  const
     842{
     843  return fAvNormBlindPixel * GetAverageQE();
    851844}
    852845
     
    857850// different colours.
    858851//
    859 Float_t MCalibrationQEPix::GetQECascadesBlindPixelVar( const Float_t zenith  )  const
    860 {
    861   return ( GetAvNormBlindPixelRelVar()  + GetAverageQERelVar(zenith))
    862     * GetQECascadesBlindPixel(zenith) * GetQECascadesBlindPixel(zenith);
     852Float_t MCalibrationQEPix::GetQECascadesBlindPixelVar()  const
     853{
     854  return ( GetAvNormBlindPixelRelVar()  + GetAverageQERelVar())
     855    * GetQECascadesBlindPixel() * GetQECascadesBlindPixel();
    863856}
    864857
     
    869862// different colours.
    870863//
    871 Float_t MCalibrationQEPix::GetQECascadesBlindPixelErr( const Float_t zenith  )  const
    872 {
    873   const Float_t var = GetQECascadesBlindPixelVar(zenith);
     864Float_t MCalibrationQEPix::GetQECascadesBlindPixelErr()  const
     865{
     866  const Float_t var = GetQECascadesBlindPixelVar();
    874867 
    875868  if (var < 0.)
     
    885878// from the different colours.
    886879//
    887 Float_t MCalibrationQEPix::GetQECascadesCombined( const Float_t zenith  ) const
    888 {
    889   return fAvNormCombined * GetAverageQE ( zenith );
     880Float_t MCalibrationQEPix::GetQECascadesCombined() const
     881{
     882  return fAvNormCombined * GetAverageQE();
    890883}
    891884
     
    896889// different colours.
    897890//
    898 Float_t MCalibrationQEPix::GetQECascadesCombinedErr( const Float_t zenith  ) const
    899 {
    900   const Float_t var = GetQECascadesCombinedVar(zenith);
     891Float_t MCalibrationQEPix::GetQECascadesCombinedErr() const
     892{
     893  const Float_t var = GetQECascadesCombinedVar();
    901894 
    902895  if (var < 0.)
     
    912905// different colours.
    913906//
    914 Float_t MCalibrationQEPix::GetQECascadesCombinedVar( const Float_t zenith  ) const
    915 {
    916   return ( GetAvNormCombinedRelVar()  + GetAverageQERelVar(zenith))
    917     * GetQECascadesCombined(zenith) * GetQECascadesCombined(zenith);
     907Float_t MCalibrationQEPix::GetQECascadesCombinedVar() const
     908{
     909  return ( GetAvNormCombinedRelVar()  + GetAverageQERelVar())
     910    * GetQECascadesCombined() * GetQECascadesCombined();
    918911}
    919912
     
    923916// with the F-Factor method and averaged over the results from the different colours.
    924917//
    925 Float_t MCalibrationQEPix::GetQECascadesFFactor( const Float_t zenith  )  const
    926 {
    927   return fAvNormFFactor * GetAverageQE ( zenith );
     918Float_t MCalibrationQEPix::GetQECascadesFFactor()  const
     919{
     920  return fAvNormFFactor * GetAverageQE();
    928921}
    929922
     
    934927// different colours.
    935928//
    936 Float_t MCalibrationQEPix::GetQECascadesFFactorErr( const Float_t zenith  )  const
    937 {
    938   const Float_t var = GetQECascadesFFactorVar(zenith);
     929Float_t MCalibrationQEPix::GetQECascadesFFactorErr()  const
     930{
     931  const Float_t var = GetQECascadesFFactorVar();
    939932 
    940933  if (var < 0.)
     
    950943// different colours.
    951944//
    952 Float_t MCalibrationQEPix::GetQECascadesFFactorVar( const Float_t zenith  )  const
    953 {
    954   return ( GetAvNormFFactorRelVar()  + GetAverageQERelVar(zenith))
    955     * GetQECascadesFFactor(zenith) * GetQECascadesFFactor(zenith);
     945Float_t MCalibrationQEPix::GetQECascadesFFactorVar()  const
     946{
     947  return ( GetAvNormFFactorRelVar()  + GetAverageQERelVar())
     948    * GetQECascadesFFactor() * GetQECascadesFFactor();
    956949}
    957950
     
    961954// with the PIN Diode method and averaged over the results from the different colours.
    962955//
    963 Float_t MCalibrationQEPix::GetQECascadesPINDiode( const Float_t zenith  )  const
    964 {
    965   return fAvNormPINDiode * GetAverageQE ( zenith );
     956Float_t MCalibrationQEPix::GetQECascadesPINDiode()  const
     957{
     958  return fAvNormPINDiode * GetAverageQE();
    966959}
    967960
     
    972965// different colours.
    973966//
    974 Float_t MCalibrationQEPix::GetQECascadesPINDiodeErr( const Float_t zenith  ) const
    975 {
    976   const Float_t var = GetQECascadesPINDiodeVar(zenith);
     967Float_t MCalibrationQEPix::GetQECascadesPINDiodeErr() const
     968{
     969  const Float_t var = GetQECascadesPINDiodeVar();
    977970 
    978971  if (var < 0.)
     
    988981// different colours.
    989982//
    990 Float_t MCalibrationQEPix::GetQECascadesPINDiodeVar( const Float_t zenith  ) const
    991 {
    992   return ( GetAvNormPINDiodeRelVar()  + GetAverageQERelVar(zenith))
    993     * GetQECascadesPINDiode(zenith) * GetQECascadesPINDiode(zenith);
     983Float_t MCalibrationQEPix::GetQECascadesPINDiodeVar() const
     984{
     985  return ( GetAvNormPINDiodeRelVar()  + GetAverageQERelVar())
     986    * GetQECascadesPINDiode() * GetQECascadesPINDiode();
    994987}
    995988
Note: See TracChangeset for help on using the changeset viewer.