source: trunk/MagicSoft/Mars/mcalib/MCalibrationQEPix.cc@ 3663

/* ======================================================================== *\
!
! *
! * This file is part of MARS, the MAGIC Analysis and Reconstruction
! * Software. It is distributed to you in the hope that it can be a useful
! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
! * It is distributed WITHOUT ANY WARRANTY.
! *
! * Permission to use, copy, modify and distribute this software and its
! * documentation for any purpose is hereby granted without fee,
! * provided that the above copyright notice appear in all copies and
! * that both that copyright notice and this permission notice appear
! * in supporting documentation. It is provided "as is" without express
! * or implied warranty.
! *
!
!
!   Author(s): Markus Gaug   02/2004 <mailto:markus@ifae.es>
!
!   Copyright: MAGIC Software Development, 2000-2004
!
!
\* ======================================================================== */

/////////////////////////////////////////////////////////////////////////////
//                                                                         //
// MCalibrationQEPix                                                       //
//                                                                         //
// Storage container of the calibrated Quantrum Efficiency of one pixel 
// For the moment, only a fixed average QE is stored:
//
// - Average QE: (email David Paneque, 14.2.04):
//
//  The conversion factor that comes purely from QE folded to a Cherenkov
//  spectrum has to be multiplied by:
//  * Plexiglass window -->> 0.96 X 0.96
//  * PMT photoelectron collection efficiency -->> 0.9
//  * Light guides efficiency -->> 0.94
//
//  Concerning the light guides effiency estimation... Daniel Ferenc 
//  is preparing some work (simulations) to estimate it. Yet so far, he has 
//  been busy with other stuff, and this work is still UNfinished.
//
//  The estimation I did comes from:
//  1) Reflectivity of light guide walls is 85 % (aluminum)
//  2) At ZERO degree light incidence, 37% of the light hits such walls 
//    (0.15X37%= 5.6% of light lost)
//  3) When increasing the light incidence angle, more and more light hits 
//     the walls.
//
//  However, the loses due to larger amount of photons hitting the walls is more 
//  or less counteracted by the fact that more and more photon trajectories cross 
//  the PMT photocathode twice, increasing the effective sensitivity of the PMT.
//
//  Jurgen Gebauer did some quick measurements about this issue. I attach a 
//  plot. You can see that the angular dependence is (more or less) in agreement 
//  with a CosTheta function (below 20-25 degrees),
//  which is the variation of teh entrance window cross section. So, in 
//  first approximation, no loses when increasing light incidence angle; 
//  and therefore, the factor 0.94.
//
//  So, summarizing... I would propose the following conversion factors 
//  (while working with CT1 cal box) in order to get the final number of photons 
//  from the detected measured size in ADC counts.
// 
//  Nph = ADC * FmethodConversionFactor / ConvPhe-PhFactor
// 
//  FmethodConversionFactor ; measured for individual pmts
// 
//  ConvPhe-PhFactor = 0.98 * 0.23 * 0.90 * 0.94 * 0.96 * 0.96 = 0.18
// 
//  I would not apply any smearing of this factor (which we have in nature), 
//  since we might be applying it to PMTs in the totally wrong direction.
// 
// 
/////////////////////////////////////////////////////////////////////////////
#include "MCalibrationQEPix.h"

#include "MLog.h"
#include "MLogManip.h"

ClassImp(MCalibrationQEPix);

using namespace std;

// --------------------------------------------------------------------------
//
// Default Constructor: 
//
MCalibrationQEPix::MCalibrationQEPix(const char *name, const char *title)
{

  fName  = name  ? name  : "MCalibrationQEPix";
  fTitle = title ? title : "Container of the calibrated quantum efficiency ";

  Clear();

}

// ------------------------------------------------------------------------
//
// Invalidate values
//
void MCalibrationQEPix::Clear(Option_t *o)
{

  fQEGreen      =  -1.;
  fQEBlue       =  -1.;
  fQEUV         =  -1.;
  fQECT1        =  -1.;
  
  fQEGreenErr   =  -1.;
  fQEBlueErr    =  -1.;
  fQEUVErr      =  -1.;
  fQECT1Err     =  -1.;
 
  MCalibrationPix::Clear();

}


void MCalibrationQEPix::SetQE( const Float_t qe, const PulserColor_t col )
{

  switch (col)
  {
      case kGREEN:
          fQEGreen = qe;
          break;
      case kBLUE:
          fQEBlue = qe;
          break;
      case kUV:
          fQEUV = qe;
          break;
      case kCT1:
          fQECT1 = qe;
          break;
      default:
          fQECT1 = qe;
          break;
  }
}

void MCalibrationQEPix::SetQEErr( const Float_t qeerr, const PulserColor_t col )
{

  switch (col)
  {
      case kGREEN:
          fQEGreenErr = qeerr;
          break;
      case kBLUE:
          fQEBlueErr  = qeerr;
          break;
      case kUV:
          fQEUVErr    = qeerr;
          break;
      case kCT1:
          fQECT1Err   = qeerr;
          break;
      default:
          fQECT1Err  = qeerr;
          break;
  }
}


Float_t MCalibrationQEPix::GetQE(const PulserColor_t col )  const
{

  switch (col)
  {
      case kGREEN:
          return fQEGreen;
          break;
      case kBLUE:
          return fQEBlue;
          break;
      case kUV:
          return fQEUV;
          break;
      case kCT1:
          return fQECT1;
          break;
      default:
          return fQECT1;
          break;
  }
}

Float_t MCalibrationQEPix::GetQEErr(const PulserColor_t col )  const
{

  switch (col)
  {
      case kGREEN:
          return fQEGreenErr;
          break;
      case kBLUE:
          return fQEBlueErr;
          break;
      case kUV:
          return fQEUVErr;
          break;
      case kCT1:
          return fQECT1Err;
          break;
      default:
          return fQECT1Err;
          break;
  }
}

// --------------------------------------------------------------
//
// The check return kTRUE if:
//
// Not yet implemented! 
//
Bool_t MCalibrationQEPix::CheckQEValidity()
{
 
  SetValid();
  return kTRUE;
}