G4Quasmon Class Reference

#include <G4Quasmon.hh>

Public Member Functions
	G4Quasmon (G4QContent qQCont=G4QContent(0, 0, 0, 0, 0, 0), G4LorentzVector q4M=G4LorentzVector(0., 0., 0., 0.), G4LorentzVector ph4M=G4LorentzVector(0., 0., 0., 0.))
	G4Quasmon (const G4Quasmon &right)
	G4Quasmon (G4Quasmon *right)
	~G4Quasmon ()
const G4Quasmon &	operator= (const G4Quasmon &right)
G4bool	operator== (const G4Quasmon &right) const
G4bool	operator!= (const G4Quasmon &right) const
G4double	GetTemper () const
G4double	GetSOverU () const
G4double	GetEtaSup () const
G4LorentzVector	Get4Momentum () const
G4QContent	GetQC () const
G4QPDGCode	GetQPDG () const
G4int	GetStatus () const
G4int	GetCharge () const
G4int	GetBaryonNumber () const
G4int	GetStrangeness () const
void	Set4Momentum (G4LorentzVector Q4M)
void	SetQC (G4QContent QQC)
void	Boost (const G4LorentzVector &theBoost)
void	Boost (const G4ThreeVector &B)
G4QHadronVector *	Fragment (G4QNucleus &nucEnviron, G4int nQ=1)
G4QHadronVector *	DecayQuasmon ()
G4QHadronVector *	DecayQHadron (G4QHadron *hadron)
void	ClearOutput ()
void	InitQuasmon (const G4QContent &qQCont, const G4LorentzVector &q4M)
void	IncreaseBy (const G4Quasmon *pQuasm)
void	IncreaseBy (G4QContent &qQCont, const G4LorentzVector &q4M)
void	ClearQuasmon ()
void	KillQuasmon ()
G4int	CalculateNumberOfQPartons (G4double qMass)

Static Public Member Functions
static void	SetParameters (G4double temper=180., G4double ssin2g=.3, G4double etaetap=.3)
static void	SetTemper (G4double temperature)
static void	SetSOverU (G4double ssin2g)
static void	SetEtaSup (G4double etaetap)
static void	OpenElectromagneticDecays ()
static void	CloseElectromagneticDecays ()

Detailed Description

Definition at line 59 of file G4Quasmon.hh.

Constructor & Destructor Documentation

G4Quasmon() [1/3]

G4Quasmon::G4Quasmon	(	G4QContent	qQCont = G4QContent(0,0,0,0,0,0),
		G4LorentzVector	q4M = G4LorentzVector(0.,0.,0.,0.),
		G4LorentzVector	ph4M = G4LorentzVector(0.,0.,0.,0.)
	)

Definition at line 74 of file G4Quasmon.cc.

 : q4Mom(q4M), valQ(qQCont), theWorld(0), phot4M(ph4M), f2all(0), rEP(0.), rMo(0.)
{
#ifdef debug
  G4cout<<"G4Quasmon:Constructor:QC="<<qQCont<<",Q4M="<<q4M<<",photonE="<<ph4M.e()<<G4endl;
#endif
#ifdef pardeb
  G4cout<<"**>G4Q:Con:(1),T="<<Temperature<<",S="<<SSin2Gluons<<",E="<<EtaEtaprime<<G4endl;
#endif
  if(phot4M.e()>0.) q4Mom+=phot4M; // InCaseOf CaptureByQuark it will be subtracted back
  valQ.DecQAQ(-1);
  status=4;                                   
}

G4Quasmon() [2/3]

G4Quasmon::G4Quasmon

(

const G4Quasmon &

right

)

Definition at line 88 of file G4Quasmon.cc.

                                          :  totMass(0), bEn(0), mbEn(0)
{
  q4Mom                 = right.q4Mom;
  valQ                  = right.valQ;
  //theEnvironment        = right.theEnvironment;
  status                = right.status;
  //theQHadrons (Vector)
  G4int nQH             = right.theQHadrons.size();
  if(nQH) for(G4int ih=0; ih<nQH; ih++)
  {
    G4QHadron* curQH    = new G4QHadron(right.theQHadrons[ih]);
    theQHadrons.push_back(curQH);
  }
  theWorld              = right.theWorld;
  phot4M                = right.phot4M;
  nBarClust             = right.nBarClust;
  nOfQ                  = right.nOfQ;
  //theQCandidates (Vector)
  G4int nQC             = right.theQCandidates.size();
  if(nQC) for(G4int iq=0; iq<nQC; iq++)
  {
    G4QCandidate* curQC = new G4QCandidate(right.theQCandidates[iq]);
    theQCandidates.push_back(curQC);
  }
  f2all                 = right.f2all;
  rEP                   = right.rEP;
  rMo                   = right.rMo;
}

G4Quasmon() [3/3]

G4Quasmon::G4Quasmon

(

G4Quasmon *

right

)

Definition at line 117 of file G4Quasmon.cc.

                                    :  totMass(0), bEn(0), mbEn(0)
{
#ifdef sdebug
  G4cout<<"G4Quasmon::Copy-Constructor: ***CALLED*** E="<<right->theEnvironment<<G4endl;
#endif
  q4Mom                 = right->q4Mom;
  valQ                  = right->valQ;
  //theEnvironment        = right->theEnvironment;
  status                = right->status;
  //theQHadrons (Vector)
  G4int nQH             = right->theQHadrons.size();
#ifdef sdebug
  G4cout<<"G4Quasmon::Copy-Constructor:nQH="<<nQH<<G4endl;
#endif
  if(nQH) for(G4int ih=0; ih<nQH; ih++)
  {
#ifdef debug
    G4cout<<"G4Quasmon:Copy-Constructor:H#"<<ih<<",QH="<<right->theQHadrons[ih]<<G4endl;
#endif
    G4QHadron* curQH    = new G4QHadron(right->theQHadrons[ih]);
    theQHadrons.push_back(curQH);
  }
  theWorld              = right->theWorld;
  phot4M                = right->phot4M;
  nBarClust             = right->nBarClust;
  nOfQ                  = right->nOfQ;
  //theQCandidates (Vector)
  G4int nQC             = right->theQCandidates.size();
#ifdef sdebug
  G4cout<<"G4Quasmon:Copy-Constructor: nCand="<<nQC<<G4endl;
#endif
  if(nQC) for(G4int iq=0; iq<nQC; iq++)
  {
#ifdef sdebug
    G4cout<<"G4Quasmon:Copy-Constructor:C#"<<iq<<",QC="<<right->theQCandidates[iq]<<G4endl;
#endif
    G4QCandidate* curQC = new G4QCandidate(right->theQCandidates[iq]);
    theQCandidates.push_back(curQC);
  }
  f2all                 = right->f2all;
  rEP                   = right->rEP;
  rMo                   = right->rMo;
#ifdef sdebug
  G4cout<<"G4Quasmon:Copy-Constructor: >->-> DONE <-<-<"<<G4endl;
#endif
}

~G4Quasmon()

G4Quasmon::~G4Quasmon

(

)

Definition at line 164 of file G4Quasmon.cc.

{
#ifdef sdebug
  G4cout<<"G4Quasmon::Destructor before theQCandidates delete"<<G4endl;
#endif
  for_each(theQCandidates.begin(), theQCandidates.end(), DeleteQCandidate());
#ifdef sdebug
  G4cout<<"G4Quasmon::Destructor before theQHadrons"<<G4endl;
#endif
  for_each(theQHadrons.begin(), theQHadrons.end(), DeleteQHadron());
#ifdef sdebug
  G4cout<<"G4Quasmon::Destructor === DONE ==="<<G4endl;
#endif
}

Member Function Documentation

Boost() [1/2]

void G4Quasmon::Boost

(

const G4LorentzVector &

theBoost

)

Definition at line 6202 of file G4Quasmon.cc.

{  
  // see CERNLIB short writeup U101 for the algorithm
  G4double bm=boost4M.mag();
  G4double factor = (q4Mom.vect()*boost4M.vect()/(boost4M.e()+bm) - q4Mom.e())/bm;
  q4Mom.setE(q4Mom.dot(boost4M)/bm);
  q4Mom.setVect(factor*boost4M.vect() + q4Mom.vect());
} // End of Boost

Referenced by G4QFragmentation::Fragment(), and G4QIonIonCollision::Fragment().

Boost() [2/2]

void G4Quasmon::Boost ( const G4ThreeVector &  B )

inline

Definition at line 99 of file G4Quasmon.hh.

{q4Mom.boost(B);} // Boosts 4-Momentum using v/c

CalculateNumberOfQPartons()

G4int G4Quasmon::CalculateNumberOfQPartons

(

G4double

qMass

)

G4int b = valQ.GetBaryonNumber(); G4int mq= 3*b; if (!b) mq=2; G4double mean = ((1.+sqrt(1.+qMOverT*qMOverT))/2. - mq)/2.; if(mean<0.) nOfQ=mq;

else nOfQ=mq+2*RandomPoisson(mean);

Definition at line 4423 of file G4Quasmon.cc.

{
  static const G4double mK0  = G4QPDGCode(311).GetMass();
  // @@ Temporary here. To have 3 quarks in Nucleon Temperature should be < M_N/4 (234 MeV)
  // M^2=4*n*(n-1)*T^2 => n = M/2T + 1/2 + T/4M + o(T^3/16M^3)
  // @@ Genius (better than 10**(-3) even for n=2!) but useless
  //G4double qMOver2T = qMass/(Temperature+Temperature);
  //G4double est = qMOver2T+1.+0.125/qMOver2T;
  // @@ Longer but exact
  G4double qMOverT = qMass/Temperature;
  G4int valc = valQ.GetTot();
  // .................................
  // --- Exponent, Double Split, Poisson 1 ----------------
  ///G4int b = valQ.GetBaryonNumber();
  ///G4int mq= 3*b;
  ///if (!b) mq=2;
  ///G4double mean = ((1.+sqrt(1.+qMOverT*qMOverT))/2. - mq)/2.;
  ///if(mean<0.) nOfQ=mq;
  // --- Uncomment up to here =-----------=^^^^^^^^^
  // Exponent ------
  //else nOfQ=mq-2*mean*log(G4UniformRand());
  // Poisson 1 ------
  ///else nOfQ=mq+2*RandomPoisson(mean);
  // Double Split ------
  //else
  //{
  //  G4int imean = static_cast<int>(mean);
  //  G4double dm = mean - imean;
  //  if(G4UniformRand()>dm) nOfQ=mq+imean+imean;
  //  else nOfQ=mq+imean+imean+2;
  //}
  // .........
  // Poisson 2 =-----------=
  //if(valc%2==0)nOfQ = 2*RandomPoisson((1.+sqrt(1.+qMOverT*qMOverT))/4.);// abs(b) is even
  //else   nOfQ = 1+2*RandomPoisson((1.+sqrt(1.+qMOverT*qMOverT))/4.-0.5);// abs(b) is odd
  // Poisson 3 =-----------=
  nOfQ = RandomPoisson((1.+sqrt(1.+qMOverT*qMOverT))/2.);
  G4int     ev = valc%2;
  if      (!ev && nOfQ<2) nOfQ=2; // #of valence quarks is even
  else if ( ev && nOfQ<3) nOfQ=3; // #of valence quarks is odd
  //
#ifdef debug
  G4cout<<"G4Q::Calc#ofQP:QM="<<q4Mom<<qMass<<",T="<<Temperature<<",QC="<<valQ<<",n="<<nOfQ
         <<G4endl;
#endif
  G4int absb = abs(valQ.GetBaryonNumber());
  G4int tabn = 0;
  if(absb)tabn=3*absb;      // Minimal QC for baryonic system fragmentation
  else    tabn=4;           // Minimal QC for mesonic system fragmentation (@@ ?)
  if (nOfQ<tabn) nOfQ=tabn;
  G4int nSeaPairs = (nOfQ-valc)/2;
  G4int stran = abs(valQ.GetS());
  G4int astra = abs(valQ.GetAS());
  if(astra>stran) stran=astra;
  G4int nMaxStrangeSea=static_cast<int>((qMass-stran*mK0)/(mK0+mK0));//KK is min for s-sea
  if (absb) nMaxStrangeSea=static_cast<int>((qMass-absb)/672.); //LambdaK is min for s-sea
#ifdef debug
  G4cout<<"G4Q::Calc#ofQP:"<<valQ<<",INtot="<<valc<<",nOfQ="<<nOfQ<<",SeaPairs="<<nSeaPairs
        <<G4endl;
#endif
  if (nSeaPairs)            // Add/subtract sea pairs to/from initial quark content
  {
#ifdef debug
    G4int morDec=0;
#endif
    if(nSeaPairs>0)valQ.IncQAQ(nSeaPairs,SSin2Gluons);
#ifdef debug
    else    morDec=valQ.DecQAQ(-nSeaPairs);
    if(morDec) G4cout<<"G4Q::Calc#ofQP: "<<morDec<<" pairs can be reduced more"<<G4endl;
#endif
    G4int sSea=valQ.GetS(); // Content of strange quarks
    G4int asSea=valQ.GetAS();
    if(asSea<sSea) sSea=asSea;
    if(sSea>nMaxStrangeSea) // @@@@@@@ Too many strange sea ??????
    {
#ifdef debug
      G4cout<<"G4Q::Calc#ofQP:**Reduce** S="<<sSea<<",aS="<<asSea<<",maxS="<<nMaxStrangeSea
            <<G4endl;
#endif
      sSea-=nMaxStrangeSea; // Strange sea excess
      valQ.DecS(sSea);      // Reduce strange sea to adoptable limit
      valQ.DecAS(sSea);
      valQ.IncQAQ(sSea,0.); // Add notstrange sea ????????
    }
  }
  // @@ Chocolate rule --- Temporary (?)
  //G4int nmin = valc+valc-2; // Chocolate
  //G4int nmin = valc+absb;   // String Junction
  //if(nOfQ<nmin) nOfQ=nmin;
  // --- End of Temporary
#ifdef debug
  G4cout<<"G4Quasmon::Calc#ofQP: *** RESULT IN*** nQ="<<nOfQ<<", FinalQC="<<valQ<<G4endl;
#endif
  return nOfQ;
} // End of "CalculateNumberOfQPartons"

ClearOutput()

void G4Quasmon::ClearOutput ( )

inline

Definition at line 168 of file G4Quasmon.hh.

      {std::for_each(theQHadrons.begin(), theQHadrons.end(), DeleteQHadron());
       theQHadrons.clear();
      }

Referenced by KillQuasmon().

ClearQuasmon()

void G4Quasmon::ClearQuasmon ( )

inline

Definition at line 200 of file G4Quasmon.hh.

{
  static const G4QContent zeroQC(0,0,0,0,0,0);
  static const G4LorentzVector nothing(0.,0.,0.,0.);
  phot4M= nothing;
  valQ  = zeroQC;
  q4Mom = nothing;
  status= 0;
  std::for_each(theQCandidates.begin(), theQCandidates.end(), DeleteQCandidate());
  theQCandidates.clear();
 
}

Referenced by KillQuasmon().

CloseElectromagneticDecays()

void G4Quasmon::CloseElectromagneticDecays ( )

static

Definition at line 189 of file G4Quasmon.cc.

{ElMaDecays=false;}

DecayQHadron()

G4QHadronVector * G4Quasmon::DecayQHadron

(

G4QHadron *

hadron

)

Definition at line 5779 of file G4Quasmon.cc.

{
  G4QHadronVector* theFragments = new G4QHadronVector;  // user is responsible to delete!
  G4QPDGCode theQPDG  = qH->GetQPDG();
  G4int      thePDG   = theQPDG.GetPDGCode();    // Get the PDG code of decaying hadron
  G4int        pap = 0;                          // --- particle
  if(thePDG<0) pap = 1;                          // --- anti-particle
  G4LorentzVector t = qH->Get4Momentum();        // Get 4-momentum of decaying hadron
  G4double m_value = t.m();                      // Get the mass value of decaying Hadron
  // --- Randomize a channel of decay
  G4QDecayChanVector decV = theWorld->GetQParticle(theQPDG)->GetDecayVector();
  G4int nChan = decV.size();
#ifdef debug
  G4cout<<"G4Quasm::DecQHadron: PDG="<<thePDG<<",m="<<m_value<<",("<<nChan<<" channels)"<<G4endl;
#endif
  if(nChan)
  {
    G4int i=0;
    if(nChan>1)
    {
      G4double rnd = G4UniformRand();            // Random value to select a Decay Channel
      for(i=0; i<nChan; i++)
      {
        G4QDecayChan* dC = decV[i];              // The i-th Decay Channel
#ifdef debug
        G4cout<<"G4Quasmon::DecaQHadr:i="<<i<<",r="<<rnd<<"<dl="<<dC->GetDecayChanLimit()
              <<", mm="<<dC->GetMinMass()<<G4endl;
#endif
        if(rnd<dC->GetDecayChanLimit() && m_value>dC->GetMinMass()) break;
      }
      if(i>nChan-1) i=nChan-1;
    }
    G4QPDGCodeVector cV=decV[i]->GetVecOfSecHadrons();// PDGVector of theSelectedDecChannel
    G4int nPart=cV.size();                         // A#of particles to decay in
#ifdef debug
    G4cout<<"G4Quasmon::DecayQHadron: resi="<<i<<",nP="<<nPart<<":"<<cV[0]->GetPDGCode()
          <<","<<cV[1]->GetPDGCode();
    if(nPart>2) G4cout<<","<<cV[2]->GetPDGCode();
    G4cout<<G4endl;
#endif
    if(nPart<2||nPart>3)
    {
      G4cerr<<"---Warning---G4Q::DecayQHadr:n="<<nPart<<",ch#"<<i<<",PDG="<<thePDG<<G4endl;
      theFragments->push_back(qH);                    // Fill as it is (del.equiv.)
      return theFragments;
    }
#ifdef debug
    G4cout<<"G4Q::DecQH:Decay("<<ElMaDecays<<") PDG="<<thePDG<<t<<m_value<<",nP="<<nPart<<G4endl;
#endif
    if(nPart==2)
    {
      G4QHadron* fHadr;
      G4QHadron* sHadr;
      G4int fPDG=cV[0]->GetPDGCode();
      G4int sPDG=cV[1]->GetPDGCode();
      // Radiative decays In2 (eta, eta', Sigma0) are closed if the ElMaDecays=false
      if ( (fPDG != 22 && sPDG != 22) || ElMaDecays) {
#ifdef debug
        G4cout<<"G4Q::DecQH:Yes2,fPDG="<<fPDG<<",sPDG="<<sPDG<<",EMF="<<ElMaDecays<<G4endl;
#endif
        if(cV[0]->GetWidth()==0.)
        { // Randomize only the second Hardon or none
          fHadr = new G4QHadron(cV[0]->GetPDGCode());   // the First Hadron is created *1*
          if(cV[1]->GetWidth()==0.)sHadr = new G4QHadron(sPDG);//theSecondHadron is created
          else
          {
            G4QParticle* sPart=theWorld->GetQParticle(cV[1]);// Pt for theSecondHadron
            G4double sdm = m_value - fHadr->GetMass();  // MaxMassLimit for the 2-nd Hadron
            sHadr = new G4QHadron(sPart,sdm);           // the Second Hadron is created *2*
            if(sPDG<0) sHadr->MakeAntiHadron();
          }
        }
        else                                              // Randomize masses ofBothHadrons
        {
          G4QParticle* sPart=theWorld->GetQParticle(cV[1]);// Pt for theSecondHadron
          G4double mim = sPart->MinMassOfFragm();         // MinMassLimit for theSecondHadr
          G4double fdm = m_value - mim;                   // MaxMassLimit for theFirstHadr
          G4QParticle* fPart=theWorld->GetQParticle(cV[0]);// Pt for the First Hadron
          fHadr = new G4QHadron(fPart,fdm);               // the 1-st Hadron is initialized
          if(fPDG<0) fHadr->MakeAntiHadron();
          G4double fm=fHadr->GetMass();                   // Mass of the first hadron
          G4double sdm = m_value - fm;                    // MaxMassLimit for theSecondHadr
          sHadr = new G4QHadron(sPart,sdm);               // the 2-nd Hadron is initialized
          if(sPDG<0) sHadr->MakeAntiHadron();
#ifdef debug
          G4cout<<"G4Q::DQH:M="<<m_value<<",mi="<<mim<<",fd="<<fdm<<",fm="<<fm<<",sd="<<sdm
                <<",sm="<<sHadr->GetMass()<<G4endl;
#endif
        }
#ifdef debug
        G4cout<<"G4Q::DQH:(DecayIn2)1="<<fHadr->GetMass()<<",2="<<sHadr->GetMass()<<G4endl;
#endif
        if(pap)
        {
          fHadr->MakeAntiHadron();
          sHadr->MakeAntiHadron();
        }
        G4LorentzVector f4Mom = fHadr->Get4Momentum();    // Get First Hadron 4Mom (mass) 
        G4LorentzVector s4Mom = sHadr->Get4Momentum();    // Get Second Hadron 4Mom (mass) 
        if(!qH->DecayIn2(f4Mom,s4Mom))                    // Error in DecayIn2
        {
          delete fHadr;                                   // Delete "new fHadr"
          delete sHadr;                                   // Delete "new sHadr"
#ifdef debug
          G4cerr<<"---Warning---G4Q::DecayQHadron:in2,PDGC="<<thePDG<<", ch#"<<i<<": 4M="
                <<qH->Get4Momentum()<<"("<<qH->GetMass()<<")->"<<f4Mom<<"+"<<s4Mom<<G4endl;
          //throw G4QException("***Exception***G4Q::DecayQHadron: Failed to decay in 2");
#endif
          theFragments->push_back(qH);                    // Fill as it is (del.equiv.)
          return theFragments;
        }
        else
        {
          //qH->SetNFragments(2);
          //theFragments.push_back(qH);                   // Fill with NFr=2 (del.equiv.)
          // Instead
          delete qH;                                      // Delete it (without History)
          //
          fHadr->Set4Momentum(f4Mom);             // Put the randomized 4Mom to 1-st Hadron
          G4QHadronVector* theTmpQHV=DecayQHadron(fHadr); // Try to decay
          G4int nProd=theTmpQHV->size();
#ifdef debug
          G4cout<<"G4Q::DecayQHadr:(DecayIn2) nOfProdForQH1="<<nProd<<G4endl;
#endif
          if(nProd==1) theFragments->push_back((*theTmpQHV)[0]);// Final = no Further Decay
          else for(G4int ip1=0; ip1<nProd; ip1++)
          {
            G4QHadronVector* intTmpQHV = DecayQHadron((*theTmpQHV)[ip1]);
            G4int tmpS=intTmpQHV->size();
            if(tmpS==1)theFragments->push_back((*intTmpQHV)[0]);// Final = no Further Decay
            else
            {
              theFragments->resize(tmpS+theFragments->size());// Resize theFragments length
              copy(intTmpQHV->begin(), intTmpQHV->end(), theFragments->end()-tmpS);
            }
#ifdef debug
            G4cout<<"G4Q::DecayQHadr:(DecayIn2) Copy Sec11 nProd="<<tmpS<<G4endl;
#endif
            intTmpQHV->clear();
            delete intTmpQHV;
          }
          theTmpQHV->clear();
          delete theTmpQHV;
          sHadr->Set4Momentum(s4Mom);             // Put the randomized 4Mom to 2-nd Hadron
          theTmpQHV=DecayQHadron(sHadr);          // Try to decay
          nProd=theTmpQHV->size();
#ifdef debug
          G4cout<<"G4Q::DecayQHadr:(DecayIn2) nOfProdForQH2="<<nProd<<G4endl;
#endif
          if(nProd==1) theFragments->push_back((*theTmpQHV)[0]);// Final = no Further Decay
          else for(G4int ip1=0; ip1<nProd; ip1++)
          {
            G4QHadronVector* intTmpQHV = DecayQHadron((*theTmpQHV)[ip1]);
            G4int tmpS=intTmpQHV->size();
            if(tmpS==1)theFragments->push_back((*intTmpQHV)[0]);// Final = no Further Decay
            else
            {
              theFragments->resize(tmpS+theFragments->size());// Resize theFragments length
              copy(intTmpQHV->begin(), intTmpQHV->end(), theFragments->end()-tmpS);
            }
#ifdef debug
            G4cout<<"G4Q::DecayQHadr:(DecayIn2) Copy Sec12 nProd="<<tmpS<<G4endl;
#endif
            intTmpQHV->clear();
            delete intTmpQHV;
          }
          theTmpQHV->clear();
          delete theTmpQHV;
        }
#ifdef debug
        G4cout<<"G4Q::DecQHadr: DecayIn2 is made with nH="<<theFragments->size()<<G4endl;
#endif
      }
      else
      {
#ifdef debug
        if(thePDG==89999003||thePDG==90002999)G4cerr<<"*G4Q::DQH:8999003/90002999"<<G4endl;
#endif
        theFragments->push_back(qH);               // Fill hadron as it is (del.equivalent)
      }
    }
    else if(nPart==3)
    {
      G4QHadron* fHadr;
      G4QHadron* sHadr;
      G4QHadron* tHadr;
      G4int fPDG=cV[0]->GetPDGCode();
      G4int sPDG=cV[1]->GetPDGCode();
      G4int tPDG=cV[2]->GetPDGCode();
      //The radiative decays of the GS hadrons In3 are closed if ElMaDecays=false
      if ( (fPDG != 22 && sPDG != 22 && tPDG != 22) || ElMaDecays)
      {
#ifdef debug
        G4cout<<"G4Q::DQH:Y,f="<<fPDG<<",s="<<sPDG<<",t="<<tPDG<<",F="<<ElMaDecays<<G4endl;
#endif
        if(cV[0]->GetWidth()==0.)                        // Don't randomize theFirstHardon
        {
          fHadr = new G4QHadron(fPDG);                   // theFirst Hadron is created  *1*
          if(cV[1]->GetWidth()==0.)
          {
            sHadr = new G4QHadron(sPDG);                 // theSecond Hadron is created *2*
            if(cV[2]->GetWidth()==0.)tHadr = new G4QHadron(tPDG);//theThirdHadron isCreated
            else
            {
              G4QParticle* tPart=theWorld->GetQParticle(cV[2]);// Pt for the3-rdH
              G4double tdm = m_value-fHadr->GetMass()-sHadr->GetMass();// MaxMass for the 2d Hadr
              tHadr = new G4QHadron(tPart,tdm);                  //the3rdHadron is created
              if(tPDG<0) tHadr->MakeAntiHadron();
            }
          }
          else                                              // Randomize 2nd & 3rd Hadrons
          {
            m_value-=fHadr->GetMass();                       // Reduce the residual MaxMass
            G4QParticle* tPart=theWorld->GetQParticle(cV[2]);// Pt for the 3-rd Hadron
            G4double mim = tPart->MinMassOfFragm();         // MinMassLimit for the 3rd Hd
            G4double sdm = m_value - mim;                   // MaxMassLimit for the 2nd Hd
            G4QParticle* sPart=theWorld->GetQParticle(cV[1]);// Pt for the 2-nd Hadron
            sHadr = new G4QHadron(sPart,sdm);               // theSecondHadron is created
            if(sPDG<0) sHadr->MakeAntiHadron();
            G4double tdm = m_value - sHadr->GetMass();      // MaxMassLimit for the 3-rd H
            tHadr = new G4QHadron(tPart,tdm);               // the Third Hadron is created
            if(tPDG<0) tHadr->MakeAntiHadron();
          }
        }
        else  // Randomize masses of all three Hadrons
        {
          G4QParticle* sPart=theWorld->GetQParticle(cV[1]); // Pt for theSecondHadr
          G4double smim = sPart->MinMassOfFragm();          // MinMassLim for SecondHadron
          G4QParticle* tPart=theWorld->GetQParticle(cV[2]); // Pt for the Third Hadron
          G4double tmim = tPart->MinMassOfFragm();          // MinMassLimit for theThirdHd
          G4double fdm = m_value - smim - tmim;             // MaxMassLimit for theFirstHd
          G4QParticle* fPart=theWorld->GetQParticle(cV[0]); // Pt for the First Hadron
          fHadr = new G4QHadron(fPart,fdm);                 // the First Hadron is created
          if(fPDG<0) fHadr->MakeAntiHadron();
          m_value-=fHadr->GetMass();                        // Reduce the residual MaxMass
          G4double  sdm = m_value - tmim;                   // MaxMassLimit for theSecondH
          sHadr = new G4QHadron(sPart,sdm);                 // theSecondHadron is created
          if(sPDG<0) sHadr->MakeAntiHadron();
          G4double  tdm = m_value - sHadr->GetMass();       // MaxMassLimit for theThird H
          tHadr = new G4QHadron(tPart,tdm);                 // the Third Hadron is created
          if(tPDG<0) tHadr->MakeAntiHadron();
        }     
#ifdef debug
        G4cout<<"G4Quasmon::DecayQHadron:3Dec. m1="<<fHadr->GetMass()
              <<",m2="<<sHadr->GetMass()<<",m3="<<tHadr->GetMass()<<G4endl;
#endif
        if(pap)
        {
          fHadr->MakeAntiHadron();
          sHadr->MakeAntiHadron();
          tHadr->MakeAntiHadron();
        }
        G4LorentzVector f4Mom = fHadr->Get4Momentum(); // Get 4M of the First Hadron (mass)
        G4LorentzVector s4Mom = sHadr->Get4Momentum(); // Get 4M of the SecondHadron (mass)
        G4LorentzVector t4Mom = tHadr->Get4Momentum(); // Get 4M of the Third Hadron (mass)
        if(!qH->DecayIn3(f4Mom,s4Mom,t4Mom))
        {
          delete fHadr;                                // Delete "new fHadr"
          delete sHadr;                                // Delete "new sHadr"
          delete tHadr;                                // Delete "new tHadr"
          G4cerr<<"---Warning---G4Q::DecayQHadron:in3,PDGC="<<thePDG<<", ch#"<<i<<G4endl;
          theFragments->push_back(qH);             // Fill as it is (delete equivalent)
          return theFragments;
        }
        else
        {
          //qH->SetNFragments(3);
          //theFragments.push_back(q);              // Fill with NFr=3 (del.equiv.)
          // Instead
          delete qH;
          //
          fHadr->Set4Momentum(f4Mom);             // Put the randomized 4Mom to 1-st Hadron
          G4QHadronVector* theTmpQHV=DecayQHadron(fHadr); // Try to decay
          G4int nProd=theTmpQHV->size();
#ifdef debug
          G4cout<<"G4Q::DecayQHadr:(DecayIn3) nOfProdForQH1="<<nProd<<G4endl;
#endif
          if(nProd==1) theFragments->push_back((*theTmpQHV)[0]);// Final = no Further Decay
          else for(G4int ip1=0; ip1<nProd; ip1++)
          {
            G4QHadronVector* intTmpQHV = DecayQHadron((*theTmpQHV)[ip1]);
            G4int tmpS=intTmpQHV->size();
            if(tmpS==1)theFragments->push_back((*intTmpQHV)[0]);// Final = no Further Decay
            else
            {
              theFragments->resize(tmpS+theFragments->size());// Resize theFragments length
              copy(intTmpQHV->begin(), intTmpQHV->end(), theFragments->end()-tmpS);
            }
#ifdef debug
            G4cout<<"G4Q::DecayQHadr:(DecayIn3) Copy Sec11 nProd="<<tmpS<<G4endl;
#endif
            intTmpQHV->clear();
            delete intTmpQHV;
          }
          theTmpQHV->clear();
          delete theTmpQHV;
 
          sHadr->Set4Momentum(s4Mom);             // Put the randomized 4Mom to 2-nd Hadron
          theTmpQHV=DecayQHadron(sHadr);          // Try to decay
          nProd=theTmpQHV->size();
#ifdef debug
          G4cout<<"G4Q::DecayQHadr:(DecayIn3) nOfProdForQH2="<<nProd<<G4endl;
#endif
          if(nProd==1) theFragments->push_back((*theTmpQHV)[0]);// Final = no Further Decay
          else for(G4int ip1=0; ip1<nProd; ip1++)
          {
            G4QHadronVector* intTmpQHV = DecayQHadron((*theTmpQHV)[ip1]);
            G4int tmpS=intTmpQHV->size();
            if(tmpS==1)theFragments->push_back((*intTmpQHV)[0]);// Final = no Further Decay
            else
            {
              theFragments->resize(tmpS+theFragments->size());// Resize theFragments length
              copy(intTmpQHV->begin(), intTmpQHV->end(), theFragments->end()-tmpS);
            }
#ifdef debug
            G4cout<<"G4Q::DecayQHadr:(DecayIn3) Copy Sec12 nProd="<<tmpS<<G4endl;
#endif
            intTmpQHV->clear();
            delete intTmpQHV;
          }
          theTmpQHV->clear();
          delete theTmpQHV;
 
          tHadr->Set4Momentum(t4Mom);             // Put the randomized 4Mom to 3-rd Hadron
          theTmpQHV=DecayQHadron(tHadr);          // Try to decay
          nProd=theTmpQHV->size();
#ifdef debug
          G4cout<<"G4Q::DecayQHadr:(DecayIn3) nOfProdForQH3="<<nProd<<G4endl;
#endif
          if(nProd==1) theFragments->push_back((*theTmpQHV)[0]);// Final = no Further Decay
          else for(G4int ip1=0; ip1<nProd; ip1++)
          {
            G4QHadronVector* intTmpQHV = DecayQHadron((*theTmpQHV)[ip1]);
            G4int tmpS=intTmpQHV->size();
            if(tmpS==1)theFragments->push_back((*intTmpQHV)[0]);// Final = no Further Decay
            else
            {
              theFragments->resize(tmpS+theFragments->size());// Resize theFragments length
              copy(intTmpQHV->begin(), intTmpQHV->end(), theFragments->end()-tmpS);
            }
#ifdef debug
            G4cout<<"G4Q::DecayQHadr:(DecayIn3) Copy Sec13 nProd="<<tmpS<<G4endl;
#endif
            intTmpQHV->clear();
            delete intTmpQHV;
          }
          theTmpQHV->clear();
          delete theTmpQHV;
 
        }
#ifdef debug
        G4cout<<"G4Q::DecQHadr: DecayIn3 is made with nH="<<theFragments->size()<<G4endl;
#endif
      }
      else theFragments->push_back(qH);            // Fill hadron as it is (del.equivalent)
    }
  }
  else
  {
#ifdef debug
    G4cout<<"G4Quas::DecQHadr:Fill PDG= "<<thePDG<<t<<m_value<<" as it is ***0***>>"<<G4endl;
#endif
    if(thePDG==89999003||thePDG==90002999)G4cerr<<"-War-G4Q::DQH:8999003/90002999"<<G4endl;
    theFragments->push_back(qH);                   // Fill as it is (delete equivalent)
  }
#ifdef debug
  G4cout<<"G4Q::DecQHadr:=-= HADRON IS DECAYED =-= with nH="<<theFragments->size()<<G4endl;
#endif
  return theFragments;
} // End of "DecayOutHadron"

Referenced by G4QFragmentation::Breeder(), G4QIonIonCollision::Breeder(), and DecayQHadron().

DecayQuasmon()

G4QHadronVector * G4Quasmon::DecayQuasmon

(

)

Definition at line 3769 of file G4Quasmon.cc.

{
  G4QHadron* thisQuasmon = new G4QHadron(valQ,q4Mom);  // create a Hadron for this Quasmon
  FillHadronVector(thisQuasmon);                       // Fill it as a hadron
  G4QHadronVector* theFragments = new G4QHadronVector; // user is responsible to delete!
  G4int nHadrs=theQHadrons.size();
#ifdef debug
  G4cout<<"G4Q::DecayQuasmon:After decay (FillHadronVector byItself) nH="<<nHadrs<<G4endl;
#endif
  if(nHadrs) for (int hadron=0; hadron<nHadrs; hadron++)
  {
    G4QHadron* curHadr = new G4QHadron(theQHadrons[hadron]);
    theFragments->push_back(curHadr);                  // (user must delete)
  }
#ifdef pdebug
  else G4cerr<<"*******G4Quasmon::DecayQuasmon: *** Nothing is in the output ***"<<G4endl;
#endif
  valQ=G4QContent(0,0,0,0,0,0);                        // Wipe the Quasmon out
  q4Mom=G4LorentzVector(0.,0.,0.,0.);                  // ... with its 4-momentum
  return theFragments;
} // End of "DecayQuasmon"

Fragment()

G4QHadronVector * G4Quasmon::Fragment	(	G4QNucleus &	nucEnviron,
		G4int	nQ = 1
	)

Definition at line 6178 of file G4Quasmon.cc.

{
#ifdef debug
  G4cout<<"G4Quasmon::Fragment called E="<<nucEnviron<<nucEnviron.GetProbability()<<G4endl;
#endif
  G4int nQs=nQ;
  HadronizeQuasmon(nucEnviron,nQs);
  G4int nHadrs=theQHadrons.size();
#ifdef debug
  G4cout<<"G4Quasm::Fragm:after HadronizeQuasmon nH="<<nHadrs<<",Env="<<nucEnviron<<G4endl;
#endif
  G4QHadronVector* theFragments = new G4QHadronVector;// user is responsible for delition !
  if(nHadrs) for (int hadron=0; hadron<nHadrs; hadron++)
  {
    G4QHadron* curHadr = new G4QHadron(theQHadrons[hadron]);
    theFragments->push_back(curHadr);         // (delete equivalent - user)
  }
#ifdef pdebug
  else G4cerr<<"*******G4Quasmon::Fragment *** Nothing is in the output ***"<<G4endl;
#endif
  return theFragments;
} // End of "Fragment"

Referenced by G4QInelastic::PostStepDoIt(), and G4QDiffractionRatio::ProjFragment().

Get4Momentum()

G4LorentzVector G4Quasmon::Get4Momentum ( ) const

inline

Definition at line 161 of file G4Quasmon.hh.

{return q4Mom;}

Referenced by G4QEnvironment::AddQuasmon(), G4QEnvironment::Fragment(), G4QFragmentation::Fragment(), G4QIonIonCollision::Fragment(), and IncreaseBy().

GetBaryonNumber()

G4int G4Quasmon::GetBaryonNumber ( ) const

inline

Definition at line 164 of file G4Quasmon.hh.

{return valQ.GetBaryonNumber();}

Referenced by G4QEnvironment::AddQuasmon(), and G4QFragmentation::Fragment().

GetCharge()

G4int G4Quasmon::GetCharge ( ) const

inline

Definition at line 163 of file G4Quasmon.hh.

{return valQ.GetCharge();}

Referenced by G4QEnvironment::AddQuasmon(), G4QEnvironment::Fragment(), G4QFragmentation::Fragment(), and G4QEnvironment::G4QEnvironment().

GetEtaSup()

G4double G4Quasmon::GetEtaSup ( ) const

inline

Definition at line 160 of file G4Quasmon.hh.

{return EtaEtaprime;}

GetQC()

G4QContent G4Quasmon::GetQC ( ) const

inline

Definition at line 162 of file G4Quasmon.hh.

{return valQ;}

Referenced by G4QEnvironment::Fragment(), G4QFragmentation::Fragment(), G4QIonIonCollision::Fragment(), G4QEnvironment::G4QEnvironment(), and IncreaseBy().

GetQPDG()

G4QPDGCode G4Quasmon::GetQPDG ( ) const

inline

Definition at line 166 of file G4Quasmon.hh.

{return G4QPDGCode(valQ);}

GetSOverU()

G4double G4Quasmon::GetSOverU ( ) const

inline

Definition at line 159 of file G4Quasmon.hh.

{return SSin2Gluons;}

GetStatus()

G4int G4Quasmon::GetStatus ( ) const

inline

Definition at line 167 of file G4Quasmon.hh.

{return status;}

GetStrangeness()

G4int G4Quasmon::GetStrangeness ( ) const

inline

Definition at line 165 of file G4Quasmon.hh.

{return valQ.GetStrangeness();}

GetTemper()

G4double G4Quasmon::GetTemper ( ) const

inline

Definition at line 158 of file G4Quasmon.hh.

{return Temperature;}

IncreaseBy() [1/2]

void G4Quasmon::IncreaseBy ( const G4Quasmon *  pQuasm )

inline

Definition at line 180 of file G4Quasmon.hh.

{
  valQ  += pQuasm->GetQC();
  q4Mom += pQuasm->Get4Momentum();
  status= 3;
}

IncreaseBy() [2/2]

void G4Quasmon::IncreaseBy ( G4QContent &  qQCont, const G4LorentzVector &  q4M  )

inline

Definition at line 186 of file G4Quasmon.hh.

{
  valQ  += qQCont;
  q4Mom += q4M;
  status= 3;
}

InitQuasmon()

void G4Quasmon::InitQuasmon ( const G4QContent &  qQCont, const G4LorentzVector &  q4M  )

inline

Definition at line 193 of file G4Quasmon.hh.

{
  valQ  = qQCont;
  q4Mom = q4M;
  status= 3;
}

KillQuasmon()

void G4Quasmon::KillQuasmon ( )

inline

Definition at line 213 of file G4Quasmon.hh.

{
  ClearQuasmon();
  ClearOutput();
}

OpenElectromagneticDecays()

void G4Quasmon::OpenElectromagneticDecays ( )

static

Definition at line 186 of file G4Quasmon.cc.

{ElMaDecays=true;}

operator!=()

G4bool G4Quasmon::operator!= ( const G4Quasmon &  right ) const

inline

Definition at line 157 of file G4Quasmon.hh.

{return this != &rhs;}

operator=()

const G4Quasmon & G4Quasmon::operator=

(

const G4Quasmon &

right

)

Definition at line 202 of file G4Quasmon.cc.

{
  if(this != &right)                          // Beware of self assignment
  {
    q4Mom                 = right.q4Mom;
    valQ                  = right.valQ;
    //theEnvironment        = right.theEnvironment;
    status                = right.status;
    //theQHadrons (Vector)
    G4int iQH             = theQHadrons.size();
    if(iQH) for(G4int jh=0; jh<iQH; jh++) delete theQHadrons[jh];
    theQHadrons.clear();
    G4int nQH             = right.theQHadrons.size();
    if(nQH) for(G4int ih=0; ih<nQH; ih++)
    {
      G4QHadron* curQH    = new G4QHadron(right.theQHadrons[ih]);
      theQHadrons.push_back(curQH);
    }
    theWorld              = right.theWorld;
    phot4M                = right.phot4M;
    nBarClust             = right.nBarClust;
    nOfQ                  = right.nOfQ;
    //theQCandidates (Vector)
    G4int iQC             = theQCandidates.size();
    if(iQC) for(G4int jq=0; jq<iQC; jq++) delete theQCandidates[jq];
    theQCandidates.clear();
    G4int nQC             = right.theQCandidates.size();
    if(nQC) for(G4int iq=0; iq<nQC; iq++)
    {
      G4QCandidate* curQC = new G4QCandidate(right.theQCandidates[iq]);
      theQCandidates.push_back(curQC);
    }
    f2all                 = right.f2all;
    rEP                   = right.rEP;
    rMo                   = right.rMo;
  }
  return *this;
} // End of "="

operator==()

G4bool G4Quasmon::operator== ( const G4Quasmon &  right ) const

inline

Definition at line 156 of file G4Quasmon.hh.

{return this == &rhs;}

Set4Momentum()

void G4Quasmon::Set4Momentum ( G4LorentzVector  Q4M )

inline

Definition at line 96 of file G4Quasmon.hh.

{q4Mom=Q4M;} // Set new value for the Quasmon 4mom

Referenced by G4QFragmentation::Fragment(), and G4QIonIonCollision::Fragment().

SetEtaSup()

void G4Quasmon::SetEtaSup ( G4double  etaetap )

static

Definition at line 200 of file G4Quasmon.cc.

{EtaEtaprime=etaetap;}

SetParameters()

void G4Quasmon::SetParameters ( G4double  temper = 180., G4double  ssin2g = .3, G4double  etaetap = .3  )

static

Definition at line 192 of file G4Quasmon.cc.

{
  Temperature=temperature;
  SSin2Gluons=ssin2g;
  EtaEtaprime=etaetap;
}

Referenced by G4ChiralInvariantPhaseSpace::ApplyYourself(), G4QAtomicElectronScattering::G4QAtomicElectronScattering(), G4QCaptureAtRest::G4QCaptureAtRest(), G4QInelastic::G4QInelastic(), G4QStringChipsParticleLevelInterface::Propagate(), G4StringChipsInterface::Propagate(), G4StringChipsParticleLevelInterface::Propagate(), G4QAtomicElectronScattering::SetParameters(), G4QCaptureAtRest::SetParameters(), and G4QInelastic::SetParameters().

SetQC()

void G4Quasmon::SetQC ( G4QContent  QQC )

inline

Definition at line 97 of file G4Quasmon.hh.

{valQ=QQC;}  // Set new Quark Cont for the Quasmon

Referenced by G4QFragmentation::Fragment(), and G4QIonIonCollision::Fragment().

SetSOverU()

void G4Quasmon::SetSOverU ( G4double  ssin2g )

static

Definition at line 199 of file G4Quasmon.cc.

{SSin2Gluons=ssin2g;}

SetTemper()

void G4Quasmon::SetTemper ( G4double  temperature )

static

Definition at line 198 of file G4Quasmon.cc.

{Temperature=temperature;}

---

The documentation for this class was generated from the following files:

• G4Quasmon.hh
• G4Quasmon.cc