GEMAmplification.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef ALICEO2_TPC_GEMAmplification_H_
#define ALICEO2_TPC_GEMAmplification_H_
 
#include "MathUtils/RandomRing.h"
#include "TPCBase/ParameterGas.h"
#include "TPCBase/ParameterGEM.h"
#include "TPCBase/CRU.h"
#include "TPCBase/PadPos.h"
#include "TPCBase/CalDet.h"
 
namespace o2
{
namespace tpc
{
 
 
class GEMAmplification
{
 public:
  static GEMAmplification& instance()
  {
    static GEMAmplification gemAmplification;
    return gemAmplification;
  }
 
  ~GEMAmplification() = default;
 
  void updateParameters();
 
  int getStackAmplification(int nElectrons = 1);
 
  int getEffectiveStackAmplification(int nElectrons = 1);
 
  int getStackAmplification(const CRU& cru, const PadPos& pos, const AmplificationMode mode, int nElectrons = 1);
 
  int getSingleGEMAmplification(int nElectrons, int GEM);
 
  int getElectronLosses(int nElectrons, float probability);
 
  int getGEMMultiplication(int nElectrons, int GEM);
 
 private:
  GEMAmplification();
 
  math_utils::RandomRing<> mRandomGaus;
  math_utils::RandomRing<> mRandomFlat;
  std::array<math_utils::RandomRing<>, 4> mGain;
  math_utils::RandomRing<> mGainFullStack;
 
  const ParameterGEM* mGEMParam; 
  const ParameterGas* mGasParam; 
  const CalPad* mGainMap;        
};
 
inline int GEMAmplification::getStackAmplification(const CRU& cru, const PadPos& pos, const AmplificationMode mode, int nElectrons)
{
  switch (mode) {
    case AmplificationMode::FullMode: {
      return static_cast<int>(static_cast<float>(getStackAmplification(nElectrons)) *
                              mGainMap->getValue(cru, pos.getRow(), pos.getPad()));
      break;
    }
    case AmplificationMode::EffectiveMode: {
      return static_cast<int>(static_cast<float>(getEffectiveStackAmplification(nElectrons)) *
                              mGainMap->getValue(cru, pos.getRow(), pos.getPad()));
      break;
    }
  }
  return nElectrons;
}
} // namespace tpc
} // namespace o2
 
#endif // ALICEO2_TPC_GEMAmplification_H_