PID.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_track_PID_H_
#define ALICEO2_track_PID_H_
 
#ifndef GPUCA_GPUCODE_DEVICE
#include <cstdint>
#endif
 
#include "GPUCommonDef.h"
#include "GPUCommonRtypes.h"
#include "CommonConstants/PhysicsConstants.h"
 
namespace o2
{
namespace track
{
namespace o2cp = o2::constants::physics;
 
namespace pid_constants // GPUs currently cannot have static constexpr array members
{
typedef uint8_t ID;
static GPUglobalconstexpr() ID NIDsTot = 19;
 
#if !defined(GPUCA_GPUCODE_DEVICE) || defined(GPUCA_GPU_DEBUG_PRINT)
GPUconstexpr() const char* sNames[NIDsTot + 1] = 
  {"Electron", "Muon", "Pion", "Kaon", "Proton", "Deuteron", "Triton", "He3", "Alpha",
   "Pion0", "Photon", "K0", "Lambda", "HyperTriton", "Hyperhydrog4", "XiMinus", "OmegaMinus", "HyperHelium4", "HyperHelium5", nullptr};
#endif
 
GPUconstexpr() const float sMasses[NIDsTot] = 
  {o2cp::MassElectron, o2cp::MassMuon, o2cp::MassPionCharged, o2cp::MassKaonCharged,
   o2cp::MassProton, o2cp::MassDeuteron, o2cp::MassTriton, o2cp::MassHelium3,
   o2cp::MassAlpha, o2cp::MassPionNeutral, o2cp::MassPhoton,
   o2cp::MassKaonNeutral, o2cp::MassLambda, o2cp::MassHyperTriton, o2cp::MassHyperhydrog4, o2cp::MassXiMinus, o2cp::MassOmegaMinus, o2cp::MassHyperHelium4, o2cp::MassHyperHelium5};
 
GPUconstexpr() const float sMasses2[NIDsTot] = 
  {o2cp::MassElectron * o2cp::MassElectron,
   o2cp::MassMuon* o2cp::MassMuon,
   o2cp::MassPionCharged* o2cp::MassPionCharged,
   o2cp::MassKaonCharged* o2cp::MassKaonCharged,
   o2cp::MassProton* o2cp::MassProton,
   o2cp::MassDeuteron* o2cp::MassDeuteron,
   o2cp::MassTriton* o2cp::MassTriton,
   o2cp::MassHelium3* o2cp::MassHelium3,
   o2cp::MassAlpha* o2cp::MassAlpha,
   o2cp::MassPionNeutral* o2cp::MassPionNeutral,
   o2cp::MassPhoton* o2cp::MassPhoton,
   o2cp::MassKaonNeutral* o2cp::MassKaonNeutral,
   o2cp::MassLambda* o2cp::MassLambda,
   o2cp::MassHyperTriton* o2cp::MassHyperTriton,
   o2cp::MassHyperhydrog4* o2cp::MassHyperhydrog4,
   o2cp::MassXiMinus* o2cp::MassXiMinus,
   o2cp::MassOmegaMinus* o2cp::MassOmegaMinus,
   o2cp::MassHyperHelium4* o2cp::MassHyperHelium4,
   o2cp::MassHyperHelium5* o2cp::MassHyperHelium5};
 
GPUconstexpr() const float sMasses2Z[NIDsTot] = 
  {o2cp::MassElectron, o2cp::MassMuon,
   o2cp::MassPionCharged, o2cp::MassKaonCharged,
   o2cp::MassProton, o2cp::MassDeuteron,
   o2cp::MassTriton, o2cp::MassHelium3 / 2.,
   o2cp::MassAlpha / 2.,
   0, 0, 0, 0, o2cp::MassHyperTriton, o2cp::MassHyperhydrog4,
   o2cp::MassXiMinus, o2cp::MassOmegaMinus,
   o2cp::MassHyperHelium4 / 2., o2cp::MassHyperHelium5 / 2.};
 
GPUconstexpr() const int sCharges[NIDsTot] = 
  {1, 1, 1, 1, 1, 1, 1, 2, 2,
   0, 0, 0, 0, 1, 1,
   1, 1,
   2, 2};
} // namespace pid_constants
 
class PID
{
 public:
  // particle identifiers, continuos starting from 0
  typedef pid_constants::ID ID;
 
  static GPUglobalconstexpr() ID Electron = 0;
  static GPUglobalconstexpr() ID Muon = 1;
  static GPUglobalconstexpr() ID Pion = 2;
  static GPUglobalconstexpr() ID Kaon = 3;
  static GPUglobalconstexpr() ID Proton = 4;
  static GPUglobalconstexpr() ID Deuteron = 5;
  static GPUglobalconstexpr() ID Triton = 6;
  static GPUglobalconstexpr() ID Helium3 = 7;
  static GPUglobalconstexpr() ID Alpha = 8;
 
  static GPUglobalconstexpr() ID First = Electron;
  static GPUglobalconstexpr() ID Last = Alpha;    
  static GPUglobalconstexpr() ID NIDs = Last + 1; 
 
  // PID for derived particles
  static GPUglobalconstexpr() ID PI0 = 9;
  static GPUglobalconstexpr() ID Photon = 10;
  static GPUglobalconstexpr() ID K0 = 11;
  static GPUglobalconstexpr() ID Lambda = 12;
  static GPUglobalconstexpr() ID HyperTriton = 13;
  static GPUglobalconstexpr() ID Hyperhydrog4 = 14;
  static GPUglobalconstexpr() ID XiMinus = 15;
  static GPUglobalconstexpr() ID OmegaMinus = 16;
  static GPUglobalconstexpr() ID HyperHelium4 = 17;
  static GPUglobalconstexpr() ID HyperHelium5 = 18;
  static GPUglobalconstexpr() ID FirstExt = PI0;
  static GPUglobalconstexpr() ID LastExt = HyperHelium5;
  static GPUglobalconstexpr() ID NIDsTot = pid_constants::NIDsTot; 
  static_assert(NIDsTot == LastExt + 1, "Incorrect NIDsTot, please update!");
 
  GPUdDefault() PID() = default;
  GPUd() PID(ID id) : mID(id) {}
  GPUd() PID(const char* name);
  GPUdDefault() PID(const PID& src) = default;
  GPUhdDefault() PID& operator=(const PID& src) = default;
 
  GPUd() ID getID() const { return mID; }
  GPUd() operator ID() const { return getID(); }
 
  GPUd() float getMass() const { return getMass(mID); }
  GPUd() float getMass2() const { return getMass2(mID); }
  GPUd() float getMass2Z() const { return getMass2Z(mID); }
  GPUd() int getCharge() const { return getCharge(mID); }
 
  GPUd() static float getMass(ID id) { return pid_constants::sMasses[id]; }
  GPUd() static float getMass2(ID id) { return pid_constants::sMasses2[id]; }
  GPUd() static float getMass2Z(ID id) { return pid_constants::sMasses2Z[id]; }
  GPUd() static int getCharge(ID id) { return pid_constants::sCharges[id]; }
#if !defined(GPUCA_GPUCODE_DEVICE) || defined(GPUCA_GPU_DEBUG_PRINT)
  GPUd() const char* getName() const
  {
    return getName(mID);
  }
  GPUd() static const char* getName(ID id) { return pid_constants::sNames[id]; }
#endif
 
 private:
  ID mID = Pion;
 
#if !defined(GPUCA_GPUCODE_DEVICE) || defined(GPUCA_GPU_DEBUG_PRINT)
  // are 2 strings equal ? (trick from Giulio)
  GPUdi() static constexpr bool sameStr(char const* x, char const* y)
  {
    return !*x && !*y ? true : /* default */ (*x == *y && sameStr(x + 1, y + 1));
  }
 
  GPUdi() static constexpr ID nameToID(char const* name, ID id)
  {
    return id > LastExt ? id : sameStr(name, pid_constants::sNames[id]) ? id
                                                                        : nameToID(name, id + 1);
  }
#endif
 
  ClassDefNV(PID, 2);
};
} // namespace track
} // namespace o2
 
#endif