LauParticlePDG.cc

Go to the documentation of this file.

// Copyright University of Warwick 2008 - 2013.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

// Authors:
// Thomas Latham
// John Back
// Paul Harrison

#include <iostream>
using std::cerr;
using std::endl;

#include "TMath.h"
#include "TParticlePDG.h"

#include "LauConstants.hh"
#include "LauDatabasePDG.hh"
#include "LauParticlePDG.hh"


ClassImp(LauParticlePDG)

LauParticlePDG::LauParticlePDG(const TParticlePDG* particle) :
        theParticle_(particle),
        antiParticle_(0),
        mass_(0.0),
        kaon_(kFALSE),
        pion_(kFALSE),
        ks_(kFALSE),
        pi0_(kFALSE),
        D_(kFALSE),
        B_(kFALSE)
{
        // store the mass and kaon? and pion?
        if ( TMath::Abs(this->code()) == 211 ) { // charged pion
                mass_ = LauConstants::mPi;
                pion_ = kTRUE;
        } else if ( this->code() == 111 ) { // neutral pion
                mass_ = LauConstants::mPi0;
                pion_ = kTRUE;
                pi0_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 321 ) { // charged kaon
                mass_ = LauConstants::mK;
                kaon_ = kTRUE;
        } else if ( this->code() == 130 || this->code() == 310 || TMath::Abs(this->code()) == 311 ) { // neutral kaon
                mass_ = LauConstants::mK0;
                kaon_ = kTRUE;
                if ( this->code() == 310 ) {
                        ks_ = kTRUE;
                }
        } else if ( this->code() == 221 ) { // eta
                mass_ = LauConstants::mEta;
        } else if ( this->code() == 331 ) { // eta'
                mass_ = LauConstants::mEtaPrime;
        } else if ( TMath::Abs(this->code()) == 411 ) { // charged D
                mass_ = LauConstants::mD;
                D_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 421 ) { // neutral D
                mass_ = LauConstants::mD0;
                D_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 431) { // charged D_s
                mass_ = LauConstants::mDs;
                D_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 511 ) { // neutral B
                mass_ = LauConstants::mB0;
                B_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 521 ) { // charged B
                mass_ = LauConstants::mB;
                B_ = kTRUE;
        } else if ( TMath::Abs(this->code()) == 531 ) { // neutral B_s
                mass_ = LauConstants::mBs0;
                B_ = kTRUE;
        }
}

LauParticlePDG::~LauParticlePDG()
{
}

Int_t LauParticlePDG::code() const
{
        if (theParticle_ == 0) {
                cerr<<"ERROR in LauParticlePDG::code : Null pointer."<<endl;
                return -999;
        }
        return theParticle_->PdgCode();
}

Int_t LauParticlePDG::conjugateCode() const
{
        if (antiParticle_ == 0) {
                cerr<<"ERROR in LauParticlePDG::conjugateCode : Null pointer."<<endl;
                return -999;
        }
        return antiParticle_->code();
}

Int_t LauParticlePDG::charge() const
{
        if (theParticle_ == 0) {
                cerr<<"ERROR in LauParticlePDG::charge : Null pointer."<<endl;
                return -999;
        }
        Double_t chargeTimesThree = theParticle_->Charge();
        Int_t theCharge = static_cast<Int_t>(chargeTimesThree/3.0);
        return theCharge;
}

TString LauParticlePDG::string() const
{
        if (theParticle_ == 0) {
                cerr<<"ERROR in LauParticlePDG::string : Null pointer."<<endl;
                return "";
        }
        return theParticle_->GetName();
}

TString LauParticlePDG::stringAlphaNum() const
{
        TString name = this->string();
        name.ReplaceAll("*","st");
        name.ReplaceAll("_","");
        name.ReplaceAll("(","");
        name.ReplaceAll(")","");
        name.ReplaceAll("+","p");
        name.ReplaceAll("-","m");
        return name;
}