Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4SingleDiffractiveExcitation Class Reference

#include <G4SingleDiffractiveExcitation.hh>

Inheritance diagram for G4SingleDiffractiveExcitation:
Collaboration diagram for G4SingleDiffractiveExcitation:

Public Member Functions

 G4SingleDiffractiveExcitation (G4double sigmaPt=0.6 *CLHEP::GeV, G4double minExtraMass=250 *CLHEP::MeV, G4double x0mass=250 *CLHEP::MeV)
 
 ~G4SingleDiffractiveExcitation ()
 
G4bool ExciteParticipants (G4VSplitableHadron *aPartner, G4VSplitableHadron *bPartner) const
 
- Public Member Functions inherited from G4QGSDiffractiveExcitation
 G4QGSDiffractiveExcitation ()
 
virtual ~G4QGSDiffractiveExcitation ()
 
virtual G4ExcitedStringString (G4VSplitableHadron *aHadron, G4bool isProjectile) const
 

Detailed Description

Definition at line 50 of file G4SingleDiffractiveExcitation.hh.

Constructor & Destructor Documentation

G4SingleDiffractiveExcitation::G4SingleDiffractiveExcitation ( G4double  sigmaPt = 0.6*CLHEP::GeV,
G4double  minExtraMass = 250*CLHEP::MeV,
G4double  x0mass = 250*CLHEP::MeV 
)

Definition at line 52 of file G4SingleDiffractiveExcitation.cc.

53 : widthOfPtSquare(-2*sqr(sigmaPt)) , minExtraMass(minextraMass), minmass(x0mass)
54 {}
T sqr(const T &x)
Definition: templates.hh:145
G4SingleDiffractiveExcitation::~G4SingleDiffractiveExcitation ( )

Definition at line 56 of file G4SingleDiffractiveExcitation.cc.

57 {}

Member Function Documentation

G4bool G4SingleDiffractiveExcitation::ExciteParticipants ( G4VSplitableHadron aPartner,
G4VSplitableHadron bPartner 
) const
virtual

Reimplemented from G4QGSDiffractiveExcitation.

Definition at line 60 of file G4SingleDiffractiveExcitation.cc.

61 {
62  G4LorentzVector Pprojectile=projectile->Get4Momentum();
63  G4double Mprojectile2=sqr(projectile->GetDefinition()->GetPDGMass() + minExtraMass);
64 
65  G4LorentzVector Ptarget=target->Get4Momentum();
66  G4double Mtarget2=sqr(target->GetDefinition()->GetPDGMass() + minExtraMass);
67  //G4cout << "E proj, target :" << Pprojectile.e() << ", " << Ptarget.e() << G4endl;
68 
69  G4bool KeepProjectile= G4UniformRand() > 0.5;
70 
71  // reset the min.mass of the non diffractive particle to its value, ( minus a bit for rounding...)
72  if ( KeepProjectile )
73  {
74  //cout << " Projectile fix" << G4endl;
75  Mprojectile2 = sqr(projectile->GetDefinition()->GetPDGMass() * (1-perCent) );
76  } else {
77  //cout << " Target fix" << G4endl;
78  Mtarget2=sqr(target->GetDefinition()->GetPDGMass() * (1-perCent) );
79  }
80 
81  // Transform momenta to cms and then rotate parallel to z axis;
82 
83  G4LorentzVector Psum;
84  Psum=Pprojectile+Ptarget;
85 
86  G4LorentzRotation toCms(-1*Psum.boostVector());
87 
88  G4LorentzVector Ptmp=toCms*Pprojectile;
89 
90  if ( Ptmp.pz() <= 0. )
91  {
92  // "String" moving backwards in CMS, abort collision !!
93  //G4cout << " abort Collision!! " << G4endl;
94  return false;
95  }
96 
97  toCms.rotateZ(-1*Ptmp.phi());
98  toCms.rotateY(-1*Ptmp.theta());
99 
100  //G4cout << "Pprojectile be4 boost " << Pprojectile << G4endl;
101  //G4cout << "Ptarget be4 boost : " << Ptarget << G4endl;
102 
103  G4LorentzRotation toLab(toCms.inverse());
104 
105  Pprojectile.transform(toCms);
106  Ptarget.transform(toCms);
107 
108  G4LorentzVector Qmomentum;
109  G4int whilecount=0;
110  do {
111  // Generate pt
112 
113  G4double maxPtSquare=sqr(Ptarget.pz());
114  if (whilecount++ >= 500 && (whilecount%100)==0)
115  //G4cout << "G4SingleDiffractiveExcitation::ExciteParticipants possibly looping"
116  // << ", loop count/ maxPtSquare : "
117  // << whilecount << " / " << maxPtSquare << G4endl;
118 
119  if (whilecount > 1000 )
120  {
121  Qmomentum=G4LorentzVector(0.,0.,0.,0.);
122  //G4cout << "G4SingleDiffractiveExcitation::ExciteParticipants: Aborting loop!" << G4endl;
123  return false; // Ignore this interaction
124  }
125  Qmomentum=G4LorentzVector(GaussianPt(widthOfPtSquare,maxPtSquare),0);
126 
127  // Momentum transfer
128  G4double Xmin = minmass / ( Pprojectile.e() + Ptarget.e() );
129  G4double Xmax=1.;
130  G4double Xplus =ChooseX(Xmin,Xmax);
131  G4double Xminus=ChooseX(Xmin,Xmax);
132 
133  G4double pt2=G4ThreeVector(Qmomentum.vect()).mag2();
134  G4double Qplus =-1 * pt2 / Xminus/Ptarget.minus();
135  G4double Qminus= pt2 / Xplus /Pprojectile.plus();
136 
137  if ( KeepProjectile )
138  {
139  Qminus = (sqr(projectile->GetDefinition()->GetPDGMass()) + pt2 )
140  / (Pprojectile.plus() + Qplus ) - Pprojectile.minus();
141  } else {
142  Qplus = Ptarget.plus() - (sqr(target->GetDefinition()->GetPDGMass()) + pt2 )
143  / (Ptarget.minus() - Qminus );
144  }
145 
146  Qmomentum.setPz( (Qplus-Qminus)/2 );
147  Qmomentum.setE( (Qplus+Qminus)/2 );
148 
149  //G4cout << "Qplus / Qminus " << Qplus << " / " << Qminus<<G4endl;
150  //G4cout << "pt2 " << pt2 << G4endl;
151  //G4cout << "Qmomentum " << Qmomentum << G4endl;
152  //G4cout << " Masses (P/T) : " << (Pprojectile+Qmomentum).mag() <<
153  // " / " << (Ptarget-Qmomentum).mag() << G4endl;
154 
155  } while ( (Ptarget-Qmomentum).mag2() <= Mtarget2 /* Loop checking, 26.10.2015, A.Ribon */
156  || (Pprojectile+Qmomentum).mag2() <= Mprojectile2
157  || (Ptarget-Qmomentum).e() < 0.
158  || (Pprojectile+Qmomentum).e() < 0. );
159 
160  //G4double Ecms=Pprojectile.e() + Ptarget.e();
161 
162  Pprojectile += Qmomentum;
163  Ptarget -= Qmomentum;
164 
165  //G4cout << "Pprojectile.e() : " << Pprojectile.e() << G4endl;
166  //G4cout << "Ptarget.e() : " << Ptarget.e() << G4endl;
167  //G4cout << "end event_______________________________________________"<<G4endl;
168  //G4cout << "Pprojectile with Q : " << Pprojectile << G4endl;
169  //G4cout << "Ptarget with Q : " << Ptarget << G4endl;
170  //G4cout << "Projectile back: " << toLab * Pprojectile << G4endl;
171  //G4cout << "Target back: " << toLab * Ptarget << G4endl;
172 
173  // Transform back and update SplitableHadron Participant.
174  Pprojectile.transform(toLab);
175  Ptarget.transform(toLab);
176 
177  //G4cout << "G4SingleDiffractiveExcitation- Target mass " << Ptarget.mag() << G4endl;
178  //G4cout << "G4SingleDiffractiveExcitation- Projectile mass " << Pprojectile.mag() << G4endl;
179 
180  target->Set4Momentum(Ptarget);
181  projectile->Set4Momentum(Pprojectile);
182 
183  return true;
184 }
Hep3Vector boostVector() const
const XML_Char * target
Definition: expat.h:268
CLHEP::Hep3Vector G4ThreeVector
static constexpr double perCent
Definition: G4SIunits.hh:332
HepLorentzVector & rotateZ(double)
int G4int
Definition: G4Types.hh:78
Hep3Vector vect() const
#define G4UniformRand()
Definition: Randomize.hh:97
bool G4bool
Definition: G4Types.hh:79
HepLorentzRotation & transform(const HepBoost &b)
HepLorentzVector & rotateY(double)
HepLorentzVector & transform(const HepRotation &)
T sqr(const T &x)
Definition: templates.hh:145
double G4double
Definition: G4Types.hh:76
CLHEP::HepLorentzVector G4LorentzVector

Here is the call graph for this function:

Here is the caller graph for this function:


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