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

using namespace std;

DecayChannel::DecayChannel() {
  fMotherPDG = kNonsensePDG;
  fBranchingRatio = 0.0;
  fNDaughters = 0;
  for (int i = 0; i < kMaxDaughters; i++)
    fDaughtersPDG[i] = kNonsensePDG;
}

DecayChannel::DecayChannel(const DecayChannel &copy) {
  fMotherPDG = copy.fMotherPDG;
  fBranchingRatio = copy.fBranchingRatio;
  fNDaughters = copy.fNDaughters;
  for (int i = 0; i < fNDaughters; i++)
    fDaughtersPDG[i] = copy.fDaughtersPDG[i];
}

DecayChannel::DecayChannel(int mother, double branching, int nDaughters, int *daughters) {
  fMotherPDG = mother;
  fBranchingRatio = branching;
  fNDaughters = 0;
  for (int i = 0; i < nDaughters; i++) {
    if (i >= kMaxDaughters) {
      edm::LogError("DecayChannel")
          << "From explicit constructor: Number of daughters bigger than the maximum allowed one (" << kMaxDaughters
          << ") !!";
    }
    fDaughtersPDG[fNDaughters++] = *(daughters + i);
  }
}

void DecayChannel::SetDaughters(int *daughters, int n) {
  for (int i = 0; i < n; i++) {
    if (i >= kMaxDaughters) {
      edm::LogError("DecayChannel") << "From SetDaughters(): Number of daughters bigger than the maximum allowed one ("
                                    << kMaxDaughters << ") !!";
    }
    fDaughtersPDG[fNDaughters++] = *(daughters + i);
  }
}

void DecayChannel::AddDaughter(int pdg) {
  if (fNDaughters >= kMaxDaughters) {
    edm::LogError("DecayChannel")
        << "From AddDaughter(): Number of daughters is already >= than the maximum allowed one (" << kMaxDaughters
        << ") !!";
  }
  fDaughtersPDG[fNDaughters++] = pdg;
}

int DecayChannel::GetDaughterPDG(int i) {
  if ((i >= fNDaughters) || (i < 0)) {
    edm::LogError("DecayChannel")
        << "From GetDaughterPDG(): Daughter index required is too big or less than zero!! There are only "
        << fNDaughters << " secondaries in this channel !!";
    return kNonsensePDG;
  }
  return fDaughtersPDG[i];
}