ElectronTransport.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_ElectronTransport_H_
#define ALICEO2_TPC_ElectronTransport_H_
 
#include "TPCBase/ParameterDetector.h"
#include "TPCBase/ParameterGas.h"
 
#include "TPCBase/Mapper.h"
#include "MathUtils/RandomRing.h"
 
namespace o2
{
namespace tpc
{
 
 
class ElectronTransport
{
 public:
  static ElectronTransport& instance()
  {
    static ElectronTransport electronTransport;
    return electronTransport;
  }
 
  ~ElectronTransport() = default;
 
  void updateParameters(float vdrift = 0);
 
  GlobalPosition3D getElectronDrift(GlobalPosition3D posEle, float& driftTime);
 
  bool isCompletelyOutOfSectorCoarseElectronDrift(GlobalPosition3D posEle, const Sector& sector) const;
 
  bool isElectronAttachment(float driftTime);
 
  float getDriftTime(float zPos, float signChange = 1.f) const;
 
 private:
  ElectronTransport();
 
  math_utils::RandomRing<> mRandomGaus;
  math_utils::RandomRing<> mRandomFlat;
  const ParameterDetector* mDetParam; 
  const ParameterGas* mGasParam;      
  float mVDrift = 0;                  
};
 
inline bool ElectronTransport::isElectronAttachment(float driftTime)
{
  if (mRandomFlat.getNextValue() < mGasParam->AttCoeff * mGasParam->OxygenCont * driftTime) {
    return true; 
  } else {
    return false; 
  }
}
 
inline float ElectronTransport::getDriftTime(float zPos, float signChange) const
{
  float time = (mDetParam->TPClength - signChange * std::abs(zPos)) / mVDrift;
  return time;
}
} // namespace tpc
} // namespace o2
 
#endif // ALICEO2_TPC_ElectronTransport_H_