Changeset 1340 for trunk/source/processes/hadronic/models/parton_string/diffraction/src/G4FTFParameters.cc
- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/parton_string/diffraction/src/G4FTFParameters.cc
r1337 r1340 25 25 // 26 26 // 27 // $Id: G4FTFParameters.cc,v 1.1 3 2009/12/16 17:51:15 gunterExp $28 // GEANT4 tag $Name: geant4-09-0 4-beta-01$27 // $Id: G4FTFParameters.cc,v 1.14 2010/09/20 15:50:46 vuzhinsk Exp $ 28 // GEANT4 tag $Name: geant4-09-03-ref-09 $ 29 29 // 30 30 … … 56 56 G4double Plab = std::sqrt(Elab * Elab - ProjectileMass*ProjectileMass); 57 57 58 G4double Ylab=0.5*std::log((Elab+Plab)/(Elab-Plab)); 59 60 Plab/=GeV; // Uzhi 8.07.10 58 61 G4double LogPlab = std::log( Plab ); 59 62 G4double sqrLogPlab = LogPlab * LogPlab; … … 186 189 SetInelasticCrossSection(Xtotal-Xelastic); 187 190 191 //G4cout<<"Xtotal, Xelastic "<<Xtotal<<" "<<Xelastic<<G4endl; 188 192 // // Interactions with elastic and inelastic collisions 189 193 SetProbabilityOfElasticScatt(Xtotal, Xelastic); … … 196 200 197 201 //----------------------------------------------------------------------------------- 202 198 203 SetSlope( Xtotal*Xtotal/16./pi/Xelastic/0.3894 ); // Slope parameter of elastic scattering 199 204 // (GeV/c)^(-2)) … … 209 214 if( absPDGcode > 1000 ) //------Projectile is baryon -------- 210 215 { 211 SetMagQuarkExchange(3.4); //3.8); 212 SetSlopeQuarkExchange(1.2); 213 SetDeltaProbAtQuarkExchange(0.1); //(0.1*4.); 214 215 SetProjMinDiffMass(1.1); // GeV 216 SetProjMinNonDiffMass(1.1); // GeV 217 SetProbabilityOfProjDiff(0.76*std::pow(s/GeV/GeV,-0.35)); 218 219 SetTarMinDiffMass(1.1); // GeV 220 SetTarMinNonDiffMass(1.1); // GeV 221 SetProbabilityOfTarDiff(0.76*std::pow(s/GeV/GeV,-0.35)); 222 223 SetAveragePt2(0.3); // GeV^2 216 SetMagQuarkExchange(1.84);//(3.63); 217 SetSlopeQuarkExchange(0.7);//(1.2); 218 SetDeltaProbAtQuarkExchange(0.); 219 220 SetProjMinDiffMass(1.16); // GeV 221 SetProjMinNonDiffMass(1.16); // GeV 222 223 SetProbabilityOfProjDiff(0.805*std::exp(-0.35*Ylab));// 0.5 224 225 SetTarMinDiffMass(1.16); // GeV 226 SetTarMinNonDiffMass(1.16); // GeV 227 228 SetProbabilityOfTarDiff(0.805*std::exp(-0.35*Ylab));// 0.5 229 230 SetAveragePt2(0.15); // 0.15 GeV^2 224 231 } 225 232 else if( absPDGcode == 211 || PDGcode == 111) //------Projectile is Pion ----------- … … 286 293 if( absPDGcode < 1000 ) 287 294 { 288 SetCofNuclearDestruction(1.); //1.0); // for meson projectile 289 } else if( theA > 20. ) 295 SetMaxNumberOfCollisions(1000.,1.); //(Plab,2.); //3.); ############################## 296 297 SetCofNuclearDestruction(0.3); //1.0); // for meson projectile 298 299 SetDofNuclearDestruction(0.4); 300 SetPt2ofNuclearDestruction(0.17*GeV*GeV); 301 SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV); 302 303 SetExcitationEnergyPerWoundedNucleon(100*MeV); 304 } else // for baryon projectile 290 305 { 291 SetCofNuclearDestruction(0.2); //2); // for baryon projectile and heavy target 292 } else 293 { 294 SetCofNuclearDestruction(0.2); //1.0); // for baryon projectile and light target 306 // SetMaxNumberOfCollisions(Plab,0.1); //6.); // ############################## 307 SetMaxNumberOfCollisions(Plab,4.); //6.); // ############################## 308 309 SetCofNuclearDestruction(0.62*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1)))); 310 311 SetDofNuclearDestruction(0.4); 312 SetPt2ofNuclearDestruction((0.035+ 313 0.04*std::exp(4.*(Ylab-2.5))/(1.+std::exp(4.*(Ylab-2.5))))*GeV*GeV); //0.09 314 SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV); 315 316 SetExcitationEnergyPerWoundedNucleon(75.*MeV); 295 317 } 296 318 297 319 SetR2ofNuclearDestruction(1.5*fermi*fermi); 298 320 299 SetExcitationEnergyPerWoundedNucleon(100*MeV); 300 301 SetDofNuclearDestruction(0.4); 302 SetPt2ofNuclearDestruction(0.17*GeV*GeV); 303 SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV); 321 //SetCofNuclearDestruction(0.47*std::exp(2.*(Ylab-2.5))/(1.+std::exp(2.*(Ylab-2.5)))); 322 //SetPt2ofNuclearDestruction((0.035+0.1*std::exp(4.*(Ylab-3.))/(1.+std::exp(4.*(Ylab-3.))))*GeV*GeV); 323 324 //SetProbabilityOfElasticScatt(1.,1.); //(Xtotal, Xelastic); 325 //SetProbabilityOfProjDiff(1.*0.62*std::pow(s/GeV/GeV,-0.51)); // 0->1 326 //SetProbabilityOfTarDiff(4.*0.62*std::pow(s/GeV/GeV,-0.51)); // 2->4 327 //SetAveragePt2(0.3); //(0.15); 328 //SetAvaragePt2ofElasticScattering(0.); 329 330 //SetCofNuclearDestruction(0.6); //(0.4); 331 SetExcitationEnergyPerWoundedNucleon(75.*MeV); //(75.*MeV); 332 //SetDofNuclearDestruction(0.6); //(0.4); 333 //SetPt2ofNuclearDestruction(0.12*GeV*GeV); //(0.168*GeV*GeV); 334 304 335 } 305 336 //**********************************************************************************************
Note: See TracChangeset
for help on using the changeset viewer.