ParticlePDG.h

Go to the documentation of this file.

/*
  Copyright   : The FASTMC and SPHMC Collaboration
  Author      : Ionut Cristian Arsene 
  Affiliation : Oslo University, Norway & Institute for Space Sciences, Bucharest, Romania
  e-mail      : i.c.arsene@fys.uio.no
  Date        : 2007/05/30

  This class is using the particle and decay lists provided by the 
  THERMINATOR (Computer Physics Communications 174 669 (2006)) and
  SHARE (Computer Physics Communications 167 229 (2005)) collaborations.
*/

#ifndef PARTICLE_PDG
#define PARTICLE_PDG

#include "Rtypes.h"

#ifndef DECAY_CHANNEL
#include "DecayChannel.h"
#endif

const int kMaxDecayChannels = 100;

class ParticlePDG {
private:
  char fName[9];
  int fPDG;
  double fMass;
  double fWidth;
  double fSpin;                    // J
  double fIsospin;                 // I
  double fIsospinZ;                // I3
  double fLightQuarkNumber;        // u, d quark number
  double fAntiLightQuarkNumber;    // u-, d- quark number
  double fStrangeQuarkNumber;      // s quark number
  double fAntiStrangeQuarkNumber;  // s- quark number
  double fCharmQuarkNumber;        // c quark number
  double fAntiCharmQuarkNumber;    // c- quark number
  int fNDecayChannels;
  DecayChannel* fDecayChannels[kMaxDecayChannels];

public:
  ParticlePDG();
  ParticlePDG(char* name, int pdg, double mass, double width);
  ~ParticlePDG();

  void AddChannel(DecayChannel& channel);
  void SetName(char* name) {
    for (int i = 0; i < 9; i++)
      if (*(name + i) != '\0')
        fName[i] = *(name + i);
      else
        break;
  }
  void SetPDG(int value) { fPDG = value; }
  void SetMass(double value) { fMass = value; }
  void SetWidth(double value) { fWidth = value; }
  void SetSpin(double value) { fSpin = value; }
  void SetIsospin(double value) { fIsospin = value; }
  void SetIsospinZ(double value) { fIsospinZ = value; }
  void SetLightQNumber(double value) { fLightQuarkNumber = value; }
  void SetLightAQNumber(double value) { fAntiLightQuarkNumber = value; }
  void SetStrangeQNumber(double value) { fStrangeQuarkNumber = value; }
  void SetStrangeAQNumber(double value) { fAntiStrangeQuarkNumber = value; }
  void SetCharmQNumber(double value) { fCharmQuarkNumber = value; }
  void SetCharmAQNumber(double value) { fAntiCharmQuarkNumber = value; }

  char* GetName() { return fName; }
  int GetPDG() { return fPDG; }
  double GetMass() { return fMass; }
  double GetWidth() { return fWidth; }
  int GetNDecayChannels() { return fNDecayChannels; }
  double GetSpin() { return fSpin; }
  double GetIsospin() { return fIsospin; }
  double GetIsospinZ() { return fIsospinZ; }
  double GetLightQNumber() { return fLightQuarkNumber; }
  double GetLightAQNumber() { return fAntiLightQuarkNumber; }
  double GetStrangeQNumber() { return fStrangeQuarkNumber; }
  double GetStrangeAQNumber() { return fAntiStrangeQuarkNumber; }
  double GetCharmQNumber() { return fCharmQuarkNumber; }
  double GetCharmAQNumber() { return fAntiCharmQuarkNumber; }
  double GetBaryonNumber() {
    return (fLightQuarkNumber + fStrangeQuarkNumber + fCharmQuarkNumber - fAntiLightQuarkNumber -
            fAntiStrangeQuarkNumber - fAntiCharmQuarkNumber) /
           3.;
  }
  double GetStrangeness() { return (fAntiStrangeQuarkNumber - fStrangeQuarkNumber); }
  double GetCharmness() { return (fCharmQuarkNumber - fAntiCharmQuarkNumber); }
  double GetElectricCharge() { return fIsospinZ + (GetBaryonNumber() + GetStrangeness() + GetCharmness()) / 2.; }

  double GetFullBranching();
  DecayChannel* GetDecayChannel(int i) {
    if (0 <= i && i < fNDecayChannels)
      return fDecayChannels[i];
    else
      return nullptr;
  }
};

#endif