PoissonSolverHelpers.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_POISSONSOLVERHELPERS_H_
#define ALICEO2_TPC_POISSONSOLVERHELPERS_H_
 
#include "CommonConstants/MathConstants.h"
 
namespace o2
{
namespace tpc
{
 
enum class CycleType {
  VCycle = 0, 
  WCycle = 1, 
  FCycle = 2  
};
 
enum class GridTransferType {
  Half = 0, 
  Full = 1, 
};
 
enum class RelaxType {
  Jacobi = 0,         
  WeightedJacobi = 1, 
  GaussSeidel = 2     
};
 
struct MGParameters {                                             
  inline static bool isFull3D = true;                             
  inline static CycleType cycleType = CycleType::FCycle;          
  inline static GridTransferType gtType = GridTransferType::Full; 
  inline static RelaxType relaxType = RelaxType::GaussSeidel;     
  inline static int nPre = 2;                                     
  inline static int nPost = 2;                                    
  inline static int nMGCycle = 200;                               
  inline static int maxLoop = 7;                                  
  inline static int gamma = 1;                                    
  inline static bool normalizeGridToOneSector = false;            
};
 
template <typename DataT = double>
struct TPCParameters {
  static constexpr DataT TPCZ0{249.525};                        
  static constexpr DataT IFCRADIUS{83.5};                       
  static constexpr DataT OFCRADIUS{254.5};                      
  static constexpr DataT ZOFFSET{0.2};                          
  static constexpr DataT DVDE{0.0024};                          
  static constexpr DataT EM{-1.602176487e-19 / 9.10938215e-31}; 
  static constexpr DataT E0{8.854187817e-12};                   
  inline static DataT cathodev{-103070.0};                      
  inline static DataT vg1t{-3260};                              
};
 
template <typename DataT = double>
struct GEMFrameParameters {
  static constexpr DataT WIDTHFRAME{1};                             
  static constexpr DataT LENGTHFRAMEIROCBOTTOM{29.195};             
  static constexpr DataT LENGTHFRAMEOROC3TOP{87.048};               
  static constexpr DataT POSBOTTOM[]{83.65, 133.5, 169.75, 207.85}; 
  static constexpr DataT POSTOP[]{133.3, 169.75, 207.85, 247.7};    
};
 
template <typename DataT = double>
struct GridProperties {
  static constexpr DataT RMIN{TPCParameters<DataT>::IFCRADIUS};                  
  static constexpr DataT ZMIN{0};                                                
  static constexpr DataT PHIMIN{0};                                              
  static constexpr DataT RMAX{TPCParameters<DataT>::OFCRADIUS};                  
  static constexpr DataT ZMAX{TPCParameters<DataT>::TPCZ0};                      
  static constexpr DataT PHIMAX{static_cast<DataT>(o2::constants::math::TwoPI)}; 
 
  static constexpr DataT getRMin() { return RMIN; }
  static constexpr DataT getZMin() { return ZMIN; }
  static constexpr DataT getPhiMin() { return PHIMIN; }
 
  static constexpr DataT getGridSpacingR(const unsigned int nR) { return (RMAX - RMIN) / (nR - 1); }
 
  static constexpr DataT getGridSpacingZ(const unsigned int nZ) { return (ZMAX - ZMIN) / (nZ - 1); }
 
  static constexpr DataT getGridSpacingPhi(const unsigned int nPhi) { return PHIMAX / nPhi; }
};
 
} // namespace tpc
} // namespace o2
 
#endif