source: trunk/source/processes/hadronic/models/parton_string/diffraction/src/G4DiffractiveSplitableHadron.cc@ 1315

//
// ********************************************************************
// * 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.          *
// ********************************************************************
//
//
// $Id: G4DiffractiveSplitableHadron.cc,v 1.8 2009/07/31 11:03:00 vuzhinsk Exp $
// GEANT4 tag $Name: geant4-09-03 $
//

// ------------------------------------------------------------
//      GEANT 4 class implementation file
//
//      ---------------- G4DiffractiveSplitableHadron----------------
//             by Gunter Folger, August 1998.
//       class splitting an interacting particle. Used by FTF String Model.
// ------------------------------------------------------------

#include "G4DiffractiveSplitableHadron.hh"

#include "G4ParticleDefinition.hh"
#include "Randomize.hh"

G4DiffractiveSplitableHadron::G4DiffractiveSplitableHadron()

{
}

G4DiffractiveSplitableHadron::G4DiffractiveSplitableHadron(const G4ReactionProduct & aPrimary)
      :  G4VSplitableHadron(aPrimary)
{
        PartonIndex=-2;
        Parton[0]=NULL;
}

G4DiffractiveSplitableHadron::G4DiffractiveSplitableHadron(const G4Nucleon & aNucleon)
      :  G4VSplitableHadron(aNucleon)
{
        PartonIndex=-2;
        Parton[0]=NULL;
}

G4DiffractiveSplitableHadron::G4DiffractiveSplitableHadron(const G4VKineticNucleon * aNucleon)
      :  G4VSplitableHadron(aNucleon)
{
        PartonIndex=-2;
        Parton[0]=NULL;
}

G4DiffractiveSplitableHadron::~G4DiffractiveSplitableHadron()
{}

const G4DiffractiveSplitableHadron & G4DiffractiveSplitableHadron::operator=(const G4DiffractiveSplitableHadron &)
{
  throw G4HadronicException(__FILE__, __LINE__, "G4DiffractiveSplitableHadron::operator= meant to not be accessable");
  return *this;
}

void G4DiffractiveSplitableHadron::SplitUp()
{
  if (IsSplit()) return;
  Splitting();
// Split once only...
        if (Parton[0] != NULL) return;

// flavours of quark ends
        
        G4int PDGcode=GetDefinition()->GetPDGEncoding();

        G4int stringStart, stringEnd;
        ChooseStringEnds(PDGcode, &stringStart,&stringEnd);

        Parton[0] = new G4Parton(stringStart);
        Parton[1] = new G4Parton(stringEnd);
        PartonIndex=-1;
}

G4Parton * G4DiffractiveSplitableHadron::GetNextParton()
{
        ++PartonIndex;
        if ( PartonIndex > 1 || PartonIndex < 0 ) return NULL;
        
        return Parton[PartonIndex];
}

G4Parton * G4DiffractiveSplitableHadron::GetNextAntiParton()
{
  return NULL; // to be looked at @@
}

//
//----------------------- Implementation--------------------------
//

void G4DiffractiveSplitableHadron::ChooseStringEnds(G4int PDGcode,G4int * aEnd, G4int * bEnd) const
{
        const G4double udspin1= 1./6.;
        const G4double uuspin1= 1./3.;
//      const G4double udspin0= 1./2.; //@
        
        G4int absPDGcode=std::abs(PDGcode);
        
        if ( absPDGcode < 1000 )   //--------------------  Meson -------------
        {
           G4int heavy=  absPDGcode/ 100;
           G4int light= (absPDGcode %100)/10;
           
//          G4int anti= std::pow(-1 , std::max( heavy, light));
            G4int anti= 1 -2 * ( std::max( heavy, light ) % 2 );
            if (PDGcode < 0 ) anti *=-1;
            
            heavy *= anti;
            light *= -1 * anti;
            
            if ( G4UniformRand() < 0.5 ) 
            {
                *aEnd=heavy;
                *bEnd=light;
            } 
            else 
            {
                *aEnd=light;
                *bEnd=heavy;
            }
        }
        else                      //-------------------- Barion --------------
        {
            G4int j1000 = PDGcode/ 1000;
            G4int j100  = (PDGcode % 1000) / 100;
            G4int j10   = (PDGcode % 100) / 10;
        
            G4double random= G4UniformRand();

        
             if ( std::abs(j100) >= std::abs(j10) )
             {              
                if ( random < udspin1 )
                {
                    *aEnd=j1000;
                    *bEnd= Diquark( j100, j10, 1);
                } else if ( random < (udspin1 + uuspin1) )
                {
                    *aEnd= j10;
                    *bEnd= Diquark( j1000, j100, 1);
                } else
                {
                    *aEnd=j100;
                       // Careful, there is no diquark q1=q2, (q1 q2)0, 
                       //   possible for Omega-
                    *bEnd= Diquark( j1000, j10, j1000 != j100 ? 0 : 1);
                }
             } else
             {
// Lambda-like hadrons have two lightest quarks in spin 0
                if ( random < udspin1 )
                {
                    *aEnd=j1000;
                                // as above, but with charmed barions
                    *bEnd= Diquark( j100, j10, j100 != j10 ? 0 : 10);
                } else if ( random < (udspin1 + uuspin1) )
                {
                    *aEnd= j10;
                    *bEnd= Diquark( j1000, j100, 1);
                } else
                {
                    *aEnd=j100;
                    *bEnd= Diquark( j1000, j10, 1);
                }
                
             } 
        
        }   
            
             
            
}

G4int G4DiffractiveSplitableHadron::Diquark(G4int aquark,G4int bquark,G4int Spin) const
{
        G4int diquarkPDG;
        
        diquarkPDG = std::max( std::abs(aquark), std::abs(bquark) )*1000 +
                     std::min( std::abs(aquark), std::abs(bquark) )*100  +
                     2*Spin + 1;
        return ( aquark > 0 && bquark > 0 ) ?  diquarkPDG : -1*diquarkPDG;
}