Coverity fixes
[u/mrichter/AliRoot.git] / TUHKMgen / UHKM / ParticlePDG.cxx
1 //
2 //  Copyright   : The FASTMC and SPHMC Collaboration
3 //  Author      : Ionut Cristian Arsene 
4 //  Affiliation : Oslo University, Norway & Institute for Space Sciences, Bucharest, Romania
5 //  e-mail      : i.c.arsene@fys.uio.no
6 //  Date        : 2007/05/30
7 //
8 //  This class is using the particle and decay lists provided by the 
9 //  THERMINATOR (Computer Physics Communications 174 669 (2006)) and
10 //  SHARE (Computer Physics Communications 167 229 (2005)) collaborations.
11 //
12
13 #include <iostream>
14 #include <string>
15 using namespace std;
16 #include "ParticlePDG.h"
17
18 //________________________________________________________________________
19 ParticlePDG::ParticlePDG() :
20   fPDG(kNonsensePDG),
21   fMass(-1.0),
22   fWidth(0.0),
23   fSpin(0.0),
24   fIsospin(0.0),
25   fIsospinZ(0.0),
26   fLightQuarkNumber(0.0),
27   fAntiLightQuarkNumber(0.0),
28   fStrangeQuarkNumber(0.0),
29   fAntiStrangeQuarkNumber(0.0),
30   fCharmQuarkNumber(0.0),
31   fAntiCharmQuarkNumber(0.0),
32   fNDecayChannels(0),
33   fStable(0.0)
34 {
35 //
36 // default constructor
37 //
38   memset(fName,'a',9);
39   for(Int_t i=0; i<kMaxDecayChannels; i++)
40     fDecayChannels[i] = new DecayChannel();
41 }
42
43 //________________________________________________________________________
44 ParticlePDG::ParticlePDG(const Char_t * const name, Int_t pdg, Double_t mass, Double_t width) :
45   fPDG(pdg),
46   fMass(mass),
47   fWidth(width),
48   fSpin(0.0),
49   fIsospin(0.0),
50   fIsospinZ(0.0),
51   fLightQuarkNumber(0.0),
52   fAntiLightQuarkNumber(0.0),
53   fStrangeQuarkNumber(0.0),
54   fAntiStrangeQuarkNumber(0.0),
55   fCharmQuarkNumber(0.0),
56   fAntiCharmQuarkNumber(0.0),
57   fNDecayChannels(0),
58   fStable(0.0)
59 {
60   //
61   // constructor
62   //
63   for(Int_t i=0; i<9; i++)
64     if(*(name+i) != '\0') fName[i] = *(name+i);
65     else break;
66   for(Int_t i=0; i<kMaxDecayChannels; i++)
67     fDecayChannels[i] = new DecayChannel();
68 }
69
70 //________________________________________________________________________
71 ParticlePDG::~ParticlePDG() {
72   //
73   // destructor
74   //
75   for(Int_t i=0; i<kMaxDecayChannels; i++)
76     delete fDecayChannels[i];
77 }
78
79 //________________________________________________________________________
80 Double_t ParticlePDG::GetFullBranching() {
81   //
82   // calculate the sum of branching ratios from all decay channels (should add up to 1)
83   //
84   Double_t fullBranching = 0.0;
85   for(Int_t i=0; i<fNDecayChannels; i++)
86     fullBranching += fDecayChannels[i]->GetBranching();
87   return fullBranching;
88 }
89
90 //________________________________________________________________________
91 void ParticlePDG::AddChannel(DecayChannel* channel) {
92   //
93   // add a decay channel
94   //
95   if(channel->GetMotherPDG() != fPDG) {
96     cout << "ERROR in ParticlePDG::AddChannel() : You try to add a channel which has a different mother PDG" << endl;
97     return;
98   }
99   fDecayChannels[fNDecayChannels]->SetMotherPDG(channel->GetMotherPDG());
100   fDecayChannels[fNDecayChannels]->SetBranching(channel->GetBranching());
101   fDecayChannels[fNDecayChannels]->SetDaughters(channel->GetDaughters(), channel->GetNDaughters());
102   fNDecayChannels++;
103 }