| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // Historic fragment from M.Komogorov; clean-u 26 // Historic fragment from M.Komogorov; clean-up still necessary @@@
27 <<
28 #include "G4ExcitedStringDecay.hh" 27 #include "G4ExcitedStringDecay.hh"
29 #include "G4SystemOfUnits.hh" <<
30 #include "G4KineticTrack.hh" 28 #include "G4KineticTrack.hh"
31 #include "G4LundStringFragmentation.hh" <<
32 #include "G4HadronicInteractionRegistry.hh" <<
33 29
34 #include "G4SampleResonance.hh" << 30 G4ExcitedStringDecay::G4ExcitedStringDecay() : G4VStringFragmentation(),
>> 31 theStringDecay(0)
>> 32 {}
35 33
36 //#define debug_G4ExcitedStringDecay << 34 G4ExcitedStringDecay::G4ExcitedStringDecay(G4VLongitudinalStringDecay * aStringDecay)
37 //#define debug_G4ExcitedStringCorr << 35 : G4VStringFragmentation(),
>> 36 theStringDecay(aStringDecay)
>> 37 {}
38 38
39 G4ExcitedStringDecay::G4ExcitedStringDecay(G4V << 39 G4ExcitedStringDecay::G4ExcitedStringDecay(const G4ExcitedStringDecay &)
40 : G4VStringFragmentation(), theStringDecay(p << 40 : G4VStringFragmentation(),
>> 41 theStringDecay(0)
41 { 42 {
42 if(!ptr) { << 43 throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay::copy ctor not accessible");
43 G4HadronicInteraction* p = << 44 }
44 G4HadronicInteractionRegistry::Instance( <<
45 theStringDecay = static_cast<G4VLongitudin <<
46 if(!theStringDecay) { theStringDecay = new <<
47 } <<
48 SetModelName(theStringDecay->GetModelName()) <<
49 } <<
50 45
51 G4ExcitedStringDecay::~G4ExcitedStringDecay() 46 G4ExcitedStringDecay::~G4ExcitedStringDecay()
52 {} << 47 {
>> 48 }
53 49
54 G4KineticTrackVector *G4ExcitedStringDecay::Fr << 50 const G4ExcitedStringDecay & G4ExcitedStringDecay::operator=(const G4ExcitedStringDecay &)
55 { 51 {
56 return theStringDecay->FragmentString(theStr << 52 throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay::operator= meant to not be accessable");
>> 53 return *this;
57 } 54 }
58 55
59 G4KineticTrackVector *G4ExcitedStringDecay::Fr << 56 int G4ExcitedStringDecay::operator==(const G4ExcitedStringDecay &) const
60 { 57 {
61 G4LorentzVector KTsum(0.,0.,0.,0.); << 58 return 0;
>> 59 }
62 60
63 #ifdef debug_G4ExcitedStringDecay << 61 int G4ExcitedStringDecay::operator!=(const G4ExcitedStringDecay &) const
64 G4cout<<G4endl; << 62 {
65 G4cout<<"--------------------------- G4Excit << 63 return 1;
66 G4cout<<"Hadronization of Excited Strings: t << 64 }
67 #endif <<
68 65
69 for ( unsigned int astring=0; astring < theS << 66 G4KineticTrackVector *G4ExcitedStringDecay::FragmentString
70 // for ( unsigned int astring=0; astring < 1 << 67 (const G4ExcitedString &theString)
71 { << 68 {
72 // G4cout<<"theStrings->operator[](astring << 69 if ( theStringDecay == NULL )
73 if ( theStrings->operator[](astring)->IsEx <<
74 {KTsum+= theStrings->operator[](astri <<
75 else {KTsum+=theStrings->operator[](astrin <<
76 } <<
77 70
78 G4LorentzRotation toCms( -1 * KTsum.boostVec << 71 theStringDecay=new G4LundStringFragmentation();
79 G4LorentzRotation toLab(toCms.inverse()); << 72
80 G4LorentzVector Ptmp; << 73 return theStringDecay->FragmentString(theString);
81 KTsum=G4LorentzVector(0.,0.,0.,0.); << 74 }
>> 75
>> 76 G4KineticTrackVector *G4ExcitedStringDecay::FragmentStrings
>> 77 (const G4ExcitedStringVector * theStrings)
>> 78 {
>> 79 G4LorentzVector KTsum(0.,0.,0.,0.);
82 80
>> 81 //G4cout<<"theStrings->size() "<<theStrings->size()<<G4endl;
83 for ( unsigned int astring=0; astring < theS 82 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
84 // for ( unsigned int astring=0; astring < 1 <<
85 { 83 {
86 if ( theStrings->operator[](astring)->IsEx << 84 KTsum+= theStrings->operator[](astring)->Get4Momentum();
87 { <<
88 Ptmp=toCms * theStrings->operator[](astr <<
89 theStrings->operator[](astring)->GetLeft <<
90 85
91 Ptmp=toCms * theStrings->operator[](astr << 86 if( !(KTsum.e()<1) && !(KTsum.e()>-1) )
92 theStrings->operator[](astring)->GetRigh << 87 {
93 << 88 throw G4HadronicException(__FILE__, __LINE__,
94 KTsum+= theStrings->operator[](astring)- << 89 "G4ExcitedStringDecay::FragmentStrings received nan string...");
95 } << 90 }
96 else <<
97 { <<
98 Ptmp=toCms * theStrings->operator[](astr <<
99 theStrings->operator[](astring)->GetKine <<
100 KTsum+= theStrings->operator[](astring)- <<
101 } <<
102 } 91 }
103 92
104 G4SampleResonance BrW; <<
105 const G4ParticleDefinition* TrackDefinition= <<
106 <<
107 G4KineticTrackVector * theResult = new G4Kin 93 G4KineticTrackVector * theResult = new G4KineticTrackVector;
108 G4int attempts(0); 94 G4int attempts(0);
109 G4bool success=false; 95 G4bool success=false;
110 G4bool NeedEnergyCorrector=false; <<
111 do { 96 do {
112 #ifdef debug_G4ExcitedStringDecay << 97 //G4cout<<"Check of momentum at string fragmentations. New try."<<G4endl;
113 G4cout<<"New try No "<<attempts<<" to <<
114 #endif <<
115 <<
116 std::for_each(theResult->begin() , theResult 98 std::for_each(theResult->begin() , theResult->end() , DeleteKineticTrack());
117 theResult->clear(); 99 theResult->clear();
118 100
119 attempts++; 101 attempts++;
120 << 102 //G4cout<<G4endl<<"attempts "<<attempts<<G4endl;
121 G4LorentzVector KTsecondaries(0.,0.,0.,0.); 103 G4LorentzVector KTsecondaries(0.,0.,0.,0.);
122 NeedEnergyCorrector=false; << 104 G4bool NeedEnergyCorrector=false;
123 105
124 for ( unsigned int astring=0; astring < theS 106 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
125 // for ( unsigned int astring=0; astri <<
126 // for ( unsigned int astring=1; astri <<
127 // for ( unsigned int astring=0; astri <<
128 { 107 {
129 #ifdef debug_G4ExcitedStringDecay << 108 //G4cout<<"String No "<<astring+1<<" "<<theStrings->operator[](astring)->Get4Momentum().mag2()<<" "<<theStrings->operator[](astring)->GetRightParton()->GetPDGcode()<<" "<<theStrings->operator[](astring)->GetLeftParton()->GetPDGcode()<<" "<<theStrings->operator[](astring)->Get4Momentum()<<G4endl;
130 G4cout<<"String No "<<astring+1<<" E << 109 //G4int Uzhi; G4cin >>Uzhi;
131 <<
132 G4cout<<"String No "<<astring+1<<" 4 <<
133 <<" "<<theStrings->operator[](astrin <<
134 #endif <<
135 <<
136 G4KineticTrackVector * generatedKine 110 G4KineticTrackVector * generatedKineticTracks = NULL;
>> 111
137 if ( theStrings->operator[](astring)->IsEx 112 if ( theStrings->operator[](astring)->IsExcited() )
138 { 113 {
139 #ifdef debug_G4ExcitedStringDecay << 114 //G4cout<<"Fragment String"<<G4endl;
140 G4cout<<"Fragment String with par <<
141 <<theStrings->operator[](as <<
142 <<theStrings->operator[](as <<
143 <<"Direction "<<theStrings- <<
144 #endif <<
145 generatedKineticTracks=FragmentString 115 generatedKineticTracks=FragmentString(*theStrings->operator[](astring));
146 #ifdef debug_G4ExcitedStringDecay <<
147 G4cout<<"(G4ExcitedStringDecay) N <<
148 <<generatedKineticTracks->s <<
149 #endif <<
150 } else { 116 } else {
151 #ifdef debug_G4ExcitedStringDecay <<
152 G4cout<<" GetTrack from the Str <<
153 #endif <<
154 G4LorentzVector Mom=theStrings->o <<
155 G4KineticTrack * aTrack= new G4Ki <<
156 theStrings->operat <<
157 theStrings->operat <<
158 G4ThreeVector(0), <<
159 <<
160 aTrack->SetPosition(theStrings->o <<
161 <<
162 #ifdef debug_G4ExcitedStringDecay <<
163 G4cout<<" A particle stored in <<
164 #endif <<
165 <<
166 generatedKineticTracks = new G4KineticT 117 generatedKineticTracks = new G4KineticTrackVector;
167 generatedKineticTracks->push_back(aTrac << 118 generatedKineticTracks->push_back(theStrings->operator[](astring)->GetKineticTrack());
168 } 119 }
169 120
170 if (generatedKineticTracks == nullptr || g << 121 if (generatedKineticTracks == NULL)
171 { 122 {
172 // G4cerr << "G4VPartonStringMode << 123 G4cerr << "G4VPartonStringModel:No KineticTracks produced" << G4endl;
173 delete generatedKineticTracks; << 124 continue;
174 generatedKineticTracks = nullptr; <<
175 continue; <<
176 } 125 }
177 126
178 G4LorentzVector KTsum1(0.,0.,0.,0.); 127 G4LorentzVector KTsum1(0.,0.,0.,0.);
179 for ( unsigned int aTrack=0; aTrack< 128 for ( unsigned int aTrack=0; aTrack<generatedKineticTracks->size();aTrack++)
180 { 129 {
181 #ifdef debug_G4ExcitedStringDecay << 130 //G4cout<<"Prod part "<<(*generatedKineticTracks)[aTrack]->GetDefinition()->GetParticleName()<<" "<<(*generatedKineticTracks)[aTrack]->Get4Momentum()<<G4endl;
182 G4cout<<"Prod part No. "<<aTrack+ <<
183 <<(*generatedKineticTracks) <<
184 <<(*generatedKineticTracks) <<
185 <<(*generatedKineticTracks) <<
186 #endif <<
187 // --------------- Sampling mass <<
188 TrackDefinition = (*generatedKine <<
189 <<
190 if (TrackDefinition->IsShortLived()) <<
191 { <<
192 G4double NewTrackMass = <<
193 BrW.SampleMass( TrackDefiniti <<
194 BrW.GetMinimu <<
195 TrackDefiniti <<
196 G4LorentzVector Tmp=G4LorentzVe <<
197 Tmp.setE(std::sqrt(sqr(NewTrack <<
198 <<
199 (*generatedKineticTracks)[aTrac <<
200 <<
201 #ifdef debug_G4ExcitedStringDec <<
202 G4cout<<"Resonance *** "<<aTrac <<
203 <<(*generatedKineticTrack <<
204 <<(*generatedKineticTrack <<
205 <<(*generatedKineticTrack <<
206 #endif <<
207 } <<
208 //------------------------------- <<
209 <<
210 theResult->push_back(generatedKin 131 theResult->push_back(generatedKineticTracks->operator[](aTrack));
211 KTsum1+= (*generatedKineticTracks 132 KTsum1+= (*generatedKineticTracks)[aTrack]->Get4Momentum();
212 } 133 }
213 KTsecondaries+=KTsum1; 134 KTsecondaries+=KTsum1;
214 135
215 #ifdef debug_G4ExcitedStringDecay << 136 if ( KTsum1.e() > 0 && std::abs((KTsum1.e()-theStrings->operator[](astring)->Get4Momentum().e()) / KTsum1.e()) > perMillion )
216 G4cout << "String secondaries(" <<ge <<
217 <<"Init string momentum: "< <<
218 <<"Final hadrons momentum: "< <<
219 #endif <<
220 <<
221 if ( KTsum1.e() > 0 && <<
222 std::abs((KTsum1.e()-theString <<
223 { 137 {
>> 138 //--debug-- G4cout << "String secondaries(" <<generatedKineticTracks->size()<< ") momentum: "
>> 139 //--debug-- << theStrings->operator[](astring)->Get4Momentum() << " " << KTsum1 << G4endl;
224 NeedEnergyCorrector=true; 140 NeedEnergyCorrector=true;
225 } 141 }
226 142
227 #ifdef debug_G4ExcitedStringDecay << 143 // clean up
228 G4cout<<"NeedEnergyCorrection yes/no <<
229 #endif <<
230 <<
231 // clean up <<
232 delete generatedKineticTracks; 144 delete generatedKineticTracks;
233 success=true; <<
234 } 145 }
>> 146 //--DEBUG G4cout << "Strings/secs total 4 momentum " << KTsum << " " <<KTsecondaries << G4endl;
235 147
>> 148 success=true;
>> 149 //G4cout<<"success "<<success<<G4endl;
236 if ( NeedEnergyCorrector ) success=EnergyAnd 150 if ( NeedEnergyCorrector ) success=EnergyAndMomentumCorrector(theResult, KTsum);
237 } while (!success && (attempts < 100)); /* << 151 //G4cout<<"success after Ecorr "<<success<<G4endl;
238 << 152 } while(!success && (attempts < 10)); // It was 100 !!! Uzhi
239 for ( unsigned int aTrack=0; aTrack<theResul << 153 //G4cout<<"End frag string"<<G4endl;
240 { <<
241 Ptmp=(*theResult)[aTrack]->Get4Momentum(); <<
242 Ptmp.transform( toLab); <<
243 (*theResult)[aTrack]->Set4Momentum(Ptmp); <<
244 } <<
245 <<
246 #ifdef debug_G4ExcitedStringDecay <<
247 G4cout<<"End of the strings fragmentation (G <<
248 <<
249 G4LorentzVector KTsum1(0.,0.,0.,0.); <<
250 154
>> 155 #ifdef debug_ExcitedStringDecay
>> 156 G4LorentzVector KTsum1=0;
251 for ( unsigned int aTrack=0; aTrack<theResul 157 for ( unsigned int aTrack=0; aTrack<theResult->size();aTrack++)
252 { 158 {
253 G4cout << " corrected tracks .. " << (*t 159 G4cout << " corrected tracks .. " << (*theResult)[aTrack]->GetDefinition()->GetParticleName()
254 <<" " << (*theResult)[aTrack]->Get4Mome << 160 <<" " << (*theResult)[aTrack]->Get4Momentum() << G4endl;
255 <<" " << (*theResult)[aTrack]->Get4Mome <<
256 KTsum1+= (*theResult)[aTrack]->Get4Momen 161 KTsum1+= (*theResult)[aTrack]->Get4Momentum();
257 } 162 }
258 << 163 G4cout << "Needcorrector/success " << NeedEnergyCorrector << "/" << success << ", Corrected total 4 momentum " << KTsum1 << G4endl;
259 G4cout << "Needcorrector/success " << NeedEn <<
260 << ", Corrected total 4 momentum " < <<
261 if ( ! success ) G4cout << "failed to correc 164 if ( ! success ) G4cout << "failed to correct E/p" << G4endl;
>> 165 #endif
262 166
263 G4cout<<"End of the Hadronization (G4Excited <<
264 #endif <<
265 <<
266 if (!success) <<
267 { <<
268 if (theResult->size() != 0) <<
269 { <<
270 std::for_each(theResult->begin() , theRe <<
271 theResult->clear(); <<
272 delete theResult; theResult=0; <<
273 } <<
274 for ( unsigned int astring=0; astring < th <<
275 // for ( unsigned int astring=0; astring < <<
276 { <<
277 if ( theStrings->operator[](astring)->Is <<
278 { <<
279 Ptmp=theStrings->operator[](astring)-> <<
280 Ptmp.transform( toLab); <<
281 theStrings->operator[](astring)->GetLe <<
282 <<
283 Ptmp=theStrings->operator[](astring)-> <<
284 Ptmp.transform( toLab); <<
285 theStrings->operator[](astring)->GetRi <<
286 } <<
287 else <<
288 { <<
289 Ptmp=theStrings->operator[](astring)-> <<
290 Ptmp.transform( toLab); <<
291 theStrings->operator[](astring)->GetKi <<
292 } <<
293 } <<
294 } <<
295 return theResult; 167 return theResult;
296 } 168 }
297 169
298 << 170 G4bool G4ExcitedStringDecay::EnergyAndMomentumCorrector
299 G4bool G4ExcitedStringDecay:: << 171 (G4KineticTrackVector* Output, G4LorentzVector& TotalCollisionMom)
300 EnergyAndMomentumCorrector(G4KineticTrackVecto << 172 {
301 { <<
302 const int nAttemptScale = 500; 173 const int nAttemptScale = 500;
303 const double ErrLimit = 1.E-5; 174 const double ErrLimit = 1.E-5;
304 if (Output->empty()) return TRUE; << 175 if (Output->empty())
>> 176 return TRUE;
305 G4LorentzVector SumMom; 177 G4LorentzVector SumMom;
306 G4double SumMass = 0; << 178 G4double SumMass = 0;
307 G4double TotalCollisionMass = Total 179 G4double TotalCollisionMass = TotalCollisionMom.m();
308 180
309 std::vector<G4double> HadronMass; G4double << 181 if( !(TotalCollisionMass<1) && !(TotalCollisionMass>-1) )
>> 182 {
>> 183 std::cout << "TotalCollisionMomentum = "<<TotalCollisionMom<<G4endl;
>> 184 throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay received nan mass...");
>> 185 }
310 186
311 #ifdef debug_G4ExcitedStringCorr << 187 //G4cout<<G4endl<<"EnergyAndMomentumCorrector "<<Output->size()<<G4endl;
312 G4cout<<G4endl<<"EnergyAndMomentumCorrecto <<
313 #endif <<
314 // Calculate sum hadron 4-momenta and summ 188 // Calculate sum hadron 4-momenta and summing hadron mass
315 unsigned int cHadron; 189 unsigned int cHadron;
316 for (cHadron = 0; cHadron < Output->size() << 190 for(cHadron = 0; cHadron < Output->size(); cHadron++)
317 { 191 {
318 SumMom += Output->operator[](cHadron) 192 SumMom += Output->operator[](cHadron)->Get4Momentum();
319 HadronM=Output->operator[](cHadron)->G << 193 if( !(SumMom<1) && !(SumMom>-1) )
320 SumMass += Output->operator[](cHadron) << 194 {
>> 195 throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay::EnergyAndMomentumCorrector() received nan momentum...");
>> 196 }
>> 197 SumMass += Output->operator[](cHadron)->GetDefinition()->GetPDGMass();
>> 198 if( !(SumMass<1) && !(SumMass>-1) )
>> 199 {
>> 200 throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay::EnergyAndMomentumCorrector() received nan mass...");
>> 201 }
321 } 202 }
322 203
323 #ifdef debug_G4ExcitedStringCorr << 204 //G4cout<<"SumMass TotalCollisionMass "<<SumMass<<" "<<TotalCollisionMass<<G4endl;
324 G4cout<<"Sum part mom "<<SumMom<<" "<<SumM <<
325 <<"Sum str mom "<<TotalCollisionMom <<
326 G4cout<<"SumMass TotalCollisionMass "<<Sum <<
327 #endif <<
328 205
329 // Cannot correct a single particle 206 // Cannot correct a single particle
330 if (Output->size() < 2) return FALSE; 207 if (Output->size() < 2) return FALSE;
331 208
332 if (SumMass > TotalCollisionMass) return F 209 if (SumMass > TotalCollisionMass) return FALSE;
333 SumMass = SumMom.m2(); 210 SumMass = SumMom.m2();
334 if (SumMass < 0) return FALSE; 211 if (SumMass < 0) return FALSE;
335 SumMass = std::sqrt(SumMass); 212 SumMass = std::sqrt(SumMass);
336 213
337 // Compute c.m.s. hadron velocity and boos << 214 // Compute c.m.s. hadron velocity and boost KTV to hadron c.m.s.
338 // G4ThreeVector Beta = -SumMom.boostVecto << 215 G4ThreeVector Beta = -SumMom.boostVector();
339 G4ThreeVector Beta = -TotalCollisionMom.bo <<
340 Output->Boost(Beta); 216 Output->Boost(Beta);
341 217
342 // Scale total c.m.s. hadron energy (hadro 218 // Scale total c.m.s. hadron energy (hadron system mass).
343 // It should be equal interaction mass 219 // It should be equal interaction mass
344 G4double Scale = 1; 220 G4double Scale = 1;
345 G4int cAttempt = 0; 221 G4int cAttempt = 0;
346 G4double Sum = 0; 222 G4double Sum = 0;
347 G4bool success = false; 223 G4bool success = false;
348 for (cAttempt = 0; cAttempt < nAttemptScal << 224 for(cAttempt = 0; cAttempt < nAttemptScale; cAttempt++)
349 { 225 {
350 Sum = 0; 226 Sum = 0;
351 for (cHadron = 0; cHadron < Output->size << 227 for(cHadron = 0; cHadron < Output->size(); cHadron++)
352 { 228 {
353 HadronM = HadronMass.at(cHadron); <<
354 G4LorentzVector HadronMom = Output->op 229 G4LorentzVector HadronMom = Output->operator[](cHadron)->Get4Momentum();
355 HadronMom.setVect(Scale*HadronMom.vect 230 HadronMom.setVect(Scale*HadronMom.vect());
356 G4double E = std::sqrt(HadronMom.vect( << 231 G4double E = std::sqrt(HadronMom.vect().mag2() + sqr(Output->operator[](cHadron)->GetDefinition()->GetPDGMass()));
357 <<
358 HadronMom.setE(E); 232 HadronMom.setE(E);
359 Output->operator[](cHadron)->Set4Momen 233 Output->operator[](cHadron)->Set4Momentum(HadronMom);
360 Sum += E; 234 Sum += E;
361 } 235 }
362 Scale = TotalCollisionMass/Sum; << 236 Scale = TotalCollisionMass/Sum;
363 #ifdef debug_G4ExcitedStringCorr << 237 #ifdef debug_G4ExcitedStringDecay
364 G4cout << "Scale-1=" << Scale -1 238 G4cout << "Scale-1=" << Scale -1
365 << ", TotalCollisionMass=" << To << 239 << ", TotalCollisionMass=" << TotalCollisionMass
366 << ", Sum=" << Sum << 240 << ", Sum=" << Sum
367 << G4endl; << 241 << G4endl;
368 #endif << 242 #endif
369 if (std::fabs(Scale - 1) <= ErrLimit) 243 if (std::fabs(Scale - 1) <= ErrLimit)
370 { 244 {
371 success = true; 245 success = true;
372 break; 246 break;
373 } 247 }
374 } 248 }
375 << 249 #ifdef debug_G4ExcitedStringDecay
376 #ifdef debug_G4ExcitedStringCorr << 250 if(!success)
377 if (!success) <<
378 { 251 {
379 G4cout << "G4ExcitedStringDecay::EnergyA 252 G4cout << "G4ExcitedStringDecay::EnergyAndMomentumCorrector - Warning"<<G4endl;
380 G4cout << " Scale not unity at end of 253 G4cout << " Scale not unity at end of iteration loop: "<<TotalCollisionMass<<" "<<Sum<<" "<<Scale<<G4endl;
381 G4cout << " Number of secondaries: " < 254 G4cout << " Number of secondaries: " << Output->size() << G4endl;
382 G4cout << " Wanted total energy: " << 255 G4cout << " Wanted total energy: " << TotalCollisionMom.e() << G4endl;
383 G4cout << " Increase number of attempt 256 G4cout << " Increase number of attempts or increase ERRLIMIT"<<G4endl;
384 // throw G4HadronicException(__FILE__, _ << 257 // throw G4HadronicException(__FILE__, __LINE__, "G4ExcitedStringDecay failed to correct...");
385 } 258 }
386 #endif << 259 #endif
387 <<
388 // Compute c.m.s. interaction velocity and 260 // Compute c.m.s. interaction velocity and KTV back boost
389 Beta = TotalCollisionMom.boostVector(); 261 Beta = TotalCollisionMom.boostVector();
390 Output->Boost(Beta); 262 Output->Boost(Beta);
391 263
392 return success; 264 return success;
393 } << 265 }
394 <<
395 266