PHOS/AliPHOSSimParam.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved.      *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 // Base class for the PHOS simulation parameters.
  17 // Do not use in the simulation; use derivative classes instead.
  18 // Author: Dmitri Peressounko, RRC KI
  19
  20 // --- AliRoot header files ---
  21 #include "AliPHOSSimParam.h"
  22 #include "AliLog.h"
  23
  24 ClassImp(AliPHOSSimParam)
  25
  26 AliPHOSSimParam  * AliPHOSSimParam::fgSimParam = 0 ;
  27 //-----------------------------------------------------------------------------
  28 AliPHOSSimParam::AliPHOSSimParam() :
  29   TNamed(),
  30   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
  31   fLightFactor(0.),fAPDFactor(0.),
  32   fAPDNoise(0.),fEMCDigitThreshold(0.),
  33   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
  34   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
  35   fEMCSubtractPedestals(kFALSE),
  36   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(0),
  37   fADCpedestalCpv(0.),fADCchanelCpv(0.),
  38   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
  39   fDigitizeE(0),fCellNonLineaityOn(1)
  40 {
  41   //Default constructor.
  42   for(Int_t i=0; i<10; i++) fDStream[i] = 0 ;
  43 }
  44
  45 //-----------------------------------------------------------------------------
  46 AliPHOSSimParam::AliPHOSSimParam(Int_t) :
  47   TNamed(),
  48   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
  49   fLightFactor(0.),fAPDFactor(0.),
  50   fAPDNoise(0.),fEMCDigitThreshold(0.),
  51   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
  52   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
  53   fEMCSubtractPedestals(kFALSE),
  54   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(0),
  55   fADCpedestalCpv(0.),fADCchanelCpv(0.),
  56   fCPVNoise(0.),fCPVDigitThreshold(0.),
  57   fNADCcpv(0),
  58   fDigitizeE(0),fCellNonLineaityOn(1)
  59 {
  60   //Real (private) constructor
  61   //Set default parameters
  62
  63   //Parameters describing energy deposition and light collection by APD, used in AliPHOSv1
  64   //Photoelectron statistics:
  65   // The light yield is a poissonian distribution of the number of
  66   // photons created in the PbWo4 crystal, calculated using following formula
  67   // NumberOfPhotons = EnergyLost * LightYieldMean* APDEfficiency
  68   // LightYieldMean is parameter calculated to be over 47000 photons per GeV
  69   // APDEfficiency is 0.02655
  70   // k_0 is 0.0045 from Valery Antonenko
  71   // The number of electrons created in the APD is
  72   // NumberOfElectrons = APDGain * LightYield
  73   // The APD Gain is 300
  74   fLightYieldMean = 47000;            //Average number of photoelectrons per GeV
  75   fIntrinsicAPDEfficiency = 0.02655 ; //APD efficiency including geometric coverage
  76 //  fLightYieldAttenuation  = 0.0045 ;  //light attenuation in PWO. Last analysis shows no z-position dependence
  77 //                                      //so we removed this dependence from simulations
  78   fLightFactor            = fLightYieldMean * fIntrinsicAPDEfficiency ; //Average number of photons collected by
  79                             //APD per GeV deposited energy
  80   fAPDFactor              = (13.418/fLightYieldMean/100.) * 300. ; //factor relating light yield and APD response
  81                             //evaluated as (13.418/fLightYieldMean/100) * APDGain ;
  82
  83
  84   //Parameters defining electronic noise calculation and Digits noise thresholds
  85   //used in AliPHOSDigitizer
  86   fAPDNoise           = 0.004 ;  // [GeV]
  87   fEMCDigitThreshold  = 2.5   ;  // [ADC counts]
  88   fEMCADCchannel      = 0.005 ;  // [GeV]
  89   fTOFa               = 0.5e-9 ; // [sec] constant term
  90   fTOFb               = 1.e-9 ;  // [sec/sqrt(GeV)]] stohastic term
  91 //  fCellNonLineaityA   = 0.18 ;   //Amp of non-linearity of cell responce
  92 //  fCellNonLineaityB   = 0.109;   //Scale of non-linearity of cell responce
  93 //  fCellNonLineaityC   = 0.976;   //Overall calibration
  94   fCellNonLineaityA   = 0. ;   //Amp of non-linearity of cell responce
  95   fCellNonLineaityB   = 0.109;   //Scale of non-linearity of cell responce
  96   fCellNonLineaityC   = 1.;   //Overall calibration
  97
  98   fADCpedestalCpv     = 0.012 ;  // [aux units]
  99   fADCchanelCpv       = 0.0012;  // [aux units]
 100   fCPVNoise           = 0.01;    // [aux units]
 101   fCPVDigitThreshold  = 0.09 ;   // [aux units]
 102   fNADCcpv  =  (Int_t)TMath::Power(2,12) ;
 103
 104   fGlobalAltroOffset = 10;
 105   fGlobalAltroThreshold = 5;
 106   fEMCSampleQualityCut = 4.;
 107
 108   //Imput streams for merging. If true => this stream contains digits (and thus noise) and not SDigits.
 109   for(Int_t i=0; i<10; i++){
 110     fDStream[i] = 0 ;
 111   }
 112   fgSimParam = this ;
 113 }
 114
 115 //-----------------------------------------------------------------------------
 116 AliPHOSSimParam::AliPHOSSimParam(const AliPHOSSimParam& ):
 117   TNamed(),
 118   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
 119   fLightFactor(0.),fAPDFactor(0.),
 120   fAPDNoise(0.),fEMCDigitThreshold(0.),
 121   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
 122   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
 123   fEMCSubtractPedestals(kFALSE),
 124   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(1.),
 125   fADCpedestalCpv(0.),fADCchanelCpv(0.),
 126   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
 127   fDigitizeE(0),fCellNonLineaityOn(1)
 128 {
 129   //Copy constructor.
 130   AliError("Should not use copy constructor for singleton") ;
 131   for(Int_t  i=0; i<10; i++){
 132     fDStream[i] = 0 ;
 133   }
 134   fgSimParam = this ;
 135 }
 136 //-----------------------------------------------------------------------------
 137 AliPHOSSimParam * AliPHOSSimParam::GetInstance(){
 138
 139   if(!fgSimParam)
 140     new AliPHOSSimParam(0) ;
 141   return fgSimParam ;
 142 }
 143 //-----------------------------------------------------------------------------
 144 AliPHOSSimParam& AliPHOSSimParam::operator = (const AliPHOSSimParam& simParam)
 145 {
 146   //Assignment operator.
 147
 148   if(this != &simParam) {
 149     AliError("Should not use operator= for singleton\n") ;
 150   }
 151
 152   return *this;
 153 }
 154