ParticlePDG.cc

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
#include "GeneratorInterface/Hydjet2Interface/interface/ParticlePDG.h"
#endif

#include <iostream>

using namespace std;

ParticlePDG::ParticlePDG() {
  fPDG = kNonsensePDG;
  fMass = -1.0;
  fWidth = 0.0;
  fNDecayChannels = 0;
  for (int i = 0; i < kMaxDecayChannels; i++)
    fDecayChannels[i] = new DecayChannel();
}

ParticlePDG::ParticlePDG(char *name, int pdg, double mass, double width) {
  for (int i = 0; i < 9; i++)
    if (*(name + i) != '\0')
      fName[i] = *(name + i);
    else
      break;
  fPDG = pdg;
  fMass = mass;
  fWidth = width;
  fNDecayChannels = 0;
  for (int i = 0; i < kMaxDecayChannels; i++)
    fDecayChannels[i] = new DecayChannel();
}

ParticlePDG::~ParticlePDG() {
  for (int i = 0; i < kMaxDecayChannels; i++)
    delete fDecayChannels[i];
}

double ParticlePDG::GetFullBranching() {
  double fullBranching = 0.0;
  for (int i = 0; i < fNDecayChannels; i++)
    fullBranching += fDecayChannels[i]->GetBranching();
  return fullBranching;
}

void ParticlePDG::AddChannel(DecayChannel &channel) {
  if (channel.GetMotherPDG() != fPDG) {
    edm::LogError("ParticlePDG") << " AddChannel() : You try to add a channel which has a different mother PDG";
    return;
  }
  fDecayChannels[fNDecayChannels]->SetMotherPDG(channel.GetMotherPDG());
  fDecayChannels[fNDecayChannels]->SetBranching(channel.GetBranching());
  fDecayChannels[fNDecayChannels]->SetDaughters(channel.GetDaughters(), channel.GetNDaughters());
  fNDecayChannels++;
}