G4SPBaryon.cc

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
#include "G4SPBaryon.hh"
#include "Randomize.hh"
#include "G4ParticleTable.hh"

// correcting numbers, HPW Dec 1999

G4int G4SPBaryon::FindQuark(G4int diQuark) const
{
  G4double sum = GetProbability(diQuark);
  G4double random = G4UniformRand();
  G4double running = 0;
  G4int Quark(0);
  typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
  iter i;
  for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
  {
    if (std::abs((*i)->GetDiQuark()) == std::abs(diQuark))
    {
      running += (*i)->GetProbability();
      if (running/sum >= random)
      {
        Quark = (*i)->GetQuark();
        break;
      }
    }
  }
  return Quark;
}

G4double G4SPBaryon::GetProbability(G4int diQuark) const
{
  G4double sum = 0;
  typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
  iter i;
  for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
  {
    if (std::abs((*i)->GetDiQuark()) == std::abs(diQuark))
    {
      sum += (*i)->GetProbability();
    }
  }
  return sum;
}

G4int G4SPBaryon::MatchDiQuarkAndGetQuark(const G4SPBaryon & aBaryon, G4int & aDiQuark) const
{
  G4int result=0;
  typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
  iter i;
  G4double running = 0;
  G4double total = 0;
  for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
  {
    total += aBaryon.GetProbability((*i)->GetDiQuark());
  }
  G4double random = G4UniformRand()*total;
  for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
  {
    running += aBaryon.GetProbability((*i)->GetDiQuark());
    if(random/total<running)
    {
      result = (*i)->GetQuark(); // (diquark annihilated)
      aDiQuark = (*i)->GetDiQuark();
      break;
    }
  }
  return result;
}

void G4SPBaryon::SampleQuarkAndDiquark(G4int & quark, G4int & diQuark) const
{
  typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
  G4double random = G4UniformRand();
  G4double sum = 0;
  iter i;
  for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
  {
    sum += (*i)->GetProbability();
    if (sum > random)
    {
      if (theDefinition->GetPDGEncoding() < 0)
      {
        quark = (*i)->GetDiQuark();
        diQuark = (*i)->GetQuark();
      } else {
        quark = (*i)->GetQuark();
        diQuark = (*i)->GetDiQuark();
      }
      break;
    }
  }
}

void G4SPBaryon::FindDiquark(G4int quark, G4int & diQuark) const
{
  typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
  G4double sum = 0;
  iter i;
  for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
  {
    if (std::abs((*i)->GetQuark()) == std::abs(quark))
    {
      sum += (*i)->GetProbability();
    }
  }
  G4double random = G4UniformRand();
  G4double running = 0;
  for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
  {
    if (std::abs((*i)->GetQuark()) == std::abs(quark))
    {
      running += (*i)->GetProbability();
      if (running/sum >= random)
      {
       diQuark = (*i)->GetDiQuark();
       break;
      }
    }
  }
}


G4SPBaryon::G4SPBaryon(G4Proton * aProton)
{
  theDefinition = aProton;
  thePartonInfo.push_back(new G4SPPartonInfo(2203, 1, 1./3.)); // uu_1, d
  thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 1./6.)); // ud_1, u
  thePartonInfo.push_back(new G4SPPartonInfo(2101, 2, 1./2.)); // ud_0, u
}

G4SPBaryon::G4SPBaryon(G4AntiProton * aAntiProton)
{
  theDefinition = aAntiProton;
  thePartonInfo.push_back(new G4SPPartonInfo(-2203, -1, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-2101, -2, 1./2.));
}

G4SPBaryon::G4SPBaryon(G4Neutron * aNeutron)
{
  theDefinition = aNeutron;
  thePartonInfo.push_back(new G4SPPartonInfo(2103, 1, 1./6.)); // ud_1, d
  thePartonInfo.push_back(new G4SPPartonInfo(2101, 1, 1./2.)); // ud_0, d
  thePartonInfo.push_back(new G4SPPartonInfo(1103, 2, 1./3.)); // dd_1, u
}

G4SPBaryon::G4SPBaryon(G4AntiNeutron * aAntiNeutron)
{
  theDefinition = aAntiNeutron;
  thePartonInfo.push_back(new G4SPPartonInfo(-2103, -1, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-2101, -1, 1./2.));
  thePartonInfo.push_back(new G4SPPartonInfo(-1103, -2, 1./3.));
}

G4SPBaryon::G4SPBaryon(G4Lambda * aLambda)
{
  theDefinition = aLambda;
  thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));  // ud_1, s
  thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./4.));  // su_1, d
  thePartonInfo.push_back(new G4SPPartonInfo(3201, 1, 1./12.)); // su_0, d
  thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./4.));  // sd_1, u
  thePartonInfo.push_back(new G4SPPartonInfo(3101, 2, 1./12.)); // sd_0, u
}

G4SPBaryon::G4SPBaryon(G4AntiLambda * aAntiLambda)
{
  theDefinition = aAntiLambda;
  thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./4.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3201, -1, 1./12.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./4.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3101, -2, 1./12.));
}

G4SPBaryon::G4SPBaryon(G4SigmaPlus * aSigmaPlus)
{
  theDefinition = aSigmaPlus;
  thePartonInfo.push_back(new G4SPPartonInfo(2203, 3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(3203, 2, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(3201, 2, 1./2.));
}

G4SPBaryon::G4SPBaryon(G4AntiSigmaPlus * aAntiSigmaPlus)
{
  theDefinition = aAntiSigmaPlus;
  thePartonInfo.push_back(new G4SPPartonInfo(-2203, -3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3203, -2, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3201, -2, 1./2.));
}

G4SPBaryon::G4SPBaryon(G4SigmaZero * aSigmaZero)
{
  theDefinition = aSigmaZero;
  thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./12.));
  thePartonInfo.push_back(new G4SPPartonInfo(3201, 1, 1./4.));
  thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./12.));
  thePartonInfo.push_back(new G4SPPartonInfo(3101, 2, 1./4.));
}

G4SPBaryon::G4SPBaryon(G4AntiSigmaZero * aAntiSigmaZero)
{
  theDefinition = aAntiSigmaZero;
  thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./12.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3201, -1, 1./4.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./12.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3101, -2, 1./4.));
}

G4SPBaryon::G4SPBaryon(G4SigmaMinus * aSigmaMinus)
{
  theDefinition = aSigmaMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(1103, 3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(3103, 1, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(3101, 1, 1./2.));
}

G4SPBaryon::G4SPBaryon(G4AntiSigmaMinus * aAntiSigmaMinus)
{
  theDefinition = aAntiSigmaMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(-1103, -3, 1./3.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3103, -1, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3101, -1, 1./2.));
}

G4SPBaryon::G4SPBaryon(G4XiMinus * aXiMinus)
{
  theDefinition = aXiMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(3103, 3, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(3101, 3, 1./2.));
  thePartonInfo.push_back(new G4SPPartonInfo(3303, 1, 1./3.));
}

G4SPBaryon::G4SPBaryon(G4AntiXiMinus * aAntiXiMinus)
{
  theDefinition = aAntiXiMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(-3103, -3, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3101, -3, 1./2.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3303, -1, 1./3.));
}

G4SPBaryon::G4SPBaryon(G4XiZero * aXiZero)
{
  theDefinition = aXiZero;
  thePartonInfo.push_back(new G4SPPartonInfo(3203, 3, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(3201, 3, 1./2.));
  thePartonInfo.push_back(new G4SPPartonInfo(3303, 2, 1./3.));
}

G4SPBaryon::G4SPBaryon(G4AntiXiZero * aAntiXiZero)
{
  theDefinition = aAntiXiZero;
  thePartonInfo.push_back(new G4SPPartonInfo(-3203, -3, 1./6.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3201, -3, 1./2.));
  thePartonInfo.push_back(new G4SPPartonInfo(-3303, -2, 1./3.));
}

G4SPBaryon::G4SPBaryon(G4OmegaMinus * anOmegaMinus)
{
  theDefinition = anOmegaMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(3303, 3, 1.));
}

G4SPBaryon::G4SPBaryon(G4AntiOmegaMinus * anAntiOmegaMinus)
{
  theDefinition = anAntiOmegaMinus;
  thePartonInfo.push_back(new G4SPPartonInfo(-3303, -3, 1.));
}

// non static particles
G4SPBaryon::G4SPBaryon(G4ParticleDefinition * aDefinition)
{
  theDefinition = aDefinition;
  if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(2224))// D++
  {
    thePartonInfo.push_back(new G4SPPartonInfo(2203, 2, 1.));
  } 
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-2224))// anti D++
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-2203, -2, 1.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(2214))// D+
  {
    thePartonInfo.push_back(new G4SPPartonInfo(2203, 1, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-2214))// anti D+
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-2203, -1, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(2114))// D0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 1, 2./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 1./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-2114))// anti D0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -1, 2./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 1./3.));
  }     
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(1114))// D-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(1103, 1, 1.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-1114))// anti D-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-1103, -1, 1.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(3224))// S*+
  {
    thePartonInfo.push_back(new G4SPPartonInfo(2203, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-3224))// anti S*+
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-2203, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(3214))// S*0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-3214))// anti S*0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(3114))// S*-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(1103, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 1, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-3224))// anti S*-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-1103, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -1, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(3324))// Xi*0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3303, 2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-3324))// anti Xi*0
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3303, -2, 2./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(3314))// Xi*-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 3, 2./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3303, 1, 1./3.));
  }
  else if(theDefinition == G4ParticleTable::GetParticleTable()->FindParticle(-3314))// anti Xi*-
  {
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -3, 2./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3303, -1, 1./3.));
  }
}

G4SPBaryon::~G4SPBaryon()
{
  for(unsigned int i=0;i<thePartonInfo.size(); i++) delete thePartonInfo[i];
}