99 if (
Parton[0] != NULL )
return;
103 G4int stringStart, stringEnd;
129 return Parton[ PartonInd ];
140 return Parton[ PartonInd ];
163 G4int* bEnd )
const {
164 G4int absPDGcode = std::abs( PDGcode );
166 if ( absPDGcode < 1000 ) {
167 G4int heavy(0), light(0);
168 if(!((absPDGcode == 111)||(absPDGcode == 221)||(absPDGcode == 331)))
170 heavy = absPDGcode/100;
171 light = (absPDGcode % 100)/10;
174 if (PDGcode < 0 ) anti *= -1;
181 else {heavy = 2; light = -2;}
191 G4int j1000 = PDGcode/1000;
192 G4int j100 = (PDGcode % 1000)/100;
193 G4int j10 = (PDGcode % 100)/10;
196 if((j1000 == j100) && (j1000 == j10)) SuppresUUDDSS=1.;
199 const G4int maxNumberOfLoops = 1000;
200 G4int loopCounter = 0;
207 if(( j100 == j10 ) && (
G4UniformRand() > SuppresUUDDSS ))
continue;
209 if( j100 == j10 ) {*bEnd =
Diquark( j100, j10, 1 );}
212 else {*bEnd =
Diquark( j100, j10, 1 );}
215 else if(random < 0.66667)
217 if(( j1000 == j10 ) && (
G4UniformRand() > SuppresUUDDSS ))
continue;
219 if( j1000 == j10 ) {*bEnd =
Diquark( j1000, j10, 1 );}
222 else {*bEnd =
Diquark( j1000, j10, 1 );}
227 if(( j1000 == j100 ) && (
G4UniformRand() > SuppresUUDDSS ))
continue;
229 if( j1000 == j100 ) {*bEnd =
Diquark( j1000, j100, 1 );}
232 else {*bEnd =
Diquark( j1000, j100, 1 );}
236 ++loopCounter < maxNumberOfLoops );
237 if ( loopCounter >= maxNumberOfLoops ) {
238 *aEnd = j10; *bEnd =
Diquark( j1000, j100, 1 );
277 G4int diquarkPDG =
std::max( std::abs( aquark ), std::abs( bquark ) ) * 1000 +
278 std::min( std::abs( aquark ), std::abs( bquark ) ) * 100 +
280 return ( aquark > 0 && bquark > 0 ) ? diquarkPDG : -1*diquarkPDG;
void ChooseStringEnds(G4int PDGcode, G4int *aEnd, G4int *bEnd) const
G4Parton * GetNextParton()
G4int Diquark(G4int aquark, G4int bquark, G4int Spin) const
G4DiffractiveSplitableHadron()
G4int GetPDGEncoding() const
G4Parton * GetNextAntiParton()
void SetFirstParton(G4int PDGcode)
const G4ParticleDefinition * GetDefinition() const
void SetSecondParton(G4int PDGcode)
~G4DiffractiveSplitableHadron()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
T min(const T t1, const T t2)
brief Return the smallest of the two arguments