Improvement of the binary cascade interface to FTF

1. Improvement of the binary cascade interface to FTF V. Uzhinsky, 27.07.09 G4BinaryCascade::G4BinaryCascade() //---------------------------------------------------------------------------- G4ReactionProductVector * G4BinaryCascade::Propagate( //---------------------------------------------------------------------------- G4KineticTrackVector * secondaries, G4V3DNucleus * nucleus) G4double StartingTime=DBL_MAX; // Search for minimal formation time Uzhi for(iter = secondaries->begin(); iter != secondaries->end(); ++iter) // Uzhi { // Uzhi if((*iter)->GetFormationTime() < StartingTime) // Uzhi StartingTime = (*iter)->GetFormationTime(); // Uzhi } // Uzhi for(iter = secondaries->begin(); iter != secondaries->end(); ++iter) { G4double FormTime = (*iter)->GetFormationTime() - StartingTime; // Uzhi (*iter)->SetFormationTime(FormTime); // Uzhi if( (*iter)->GetState() == G4KineticTrack::undefined ) { FindLateParticleCollision(*iter); } else { theSecondaryList.push_back(*iter); } theProjectileList.push_back(new G4KineticTrack(*(*iter))); } Synchronization of FTF and Bic times.

2. Improvement of the binary cascade interface to FTF G4bool G4BinaryCascade::ApplyCollision(G4CollisionInitialState * collision) //---------------------------------------------------------------------------- G4KineticTrack * primary = collision->GetPrimary(); G4KineticTrackVector target_collection=collision->GetTargetCollection(); G4bool haveTarget=target_collection.size()>0; // for primary resonances, subtract neutron ( = proton) field ( ie. add std::abs(field)) G4int PDGcode=std::abs(primary->GetDefinition()->GetPDGEncoding()); mom4Primary=primary->Get4Momentum(); initial_Efermi=RKprop->GetField(primary->GetDefinition()->GetPDGEncoding(),primary->GetPosition()); if(!haveTarget) // Uzhi { // Uzhi if ( PDGcode > 1000 ) { // Uzhi initial_Efermi = RKprop->GetField(G4Neutron::Neutron()->GetPDGEncoding(), // Uzhi primary->GetPosition()); // Uzhi primary->Update4Momentum(mom4Primary.e() - initial_Efermi); // Uzhi } } else // Uzhi { // Uzhi if ( PDGcode > 1000 && PDGcode != 2112 && PDGcode != 2212 ) { initial_Efermi = RKprop->GetField(G4Neutron::Neutron()->GetPDGEncoding(), primary->GetPosition()); primary->Update4Momentum(mom4Primary.e() - initial_Efermi); } } // Uzhi Putting FTF particles in the potential well.

3. Improvement of the binary cascade interface to FTF if ( theCapturedList.size() == 1 ) // Uzhi { // Uzhi i=theCapturedList.begin(); // Uzhi G4ReactionProduct * aNew = new G4ReactionProduct((*i)- >GetDefinition()); // Uzhi aNew->SetTotalEnergy((*i)->GetDefinition()->GetPDGMass()); // Uzhi aNew->SetMomentum(G4ThreeVector(0));// see boost below.. // Uzhi precompoundProducts->push_back(aNew); // Uzhi } // Uzhi Transferring of the captured particles into produced ones.

4. Improvement of the binary cascade interface to FTF FTF + Binary cascade p+Fe 10,000 events per Plab Plab Kin./Tot.energy Number of violeted MaxDeviation (MeV) 2 1.270880e+03 0 4 3.170299e+03 3 3.377268e-02 6 5.134648e+03 1 1.080265e-02 8 7.116562e+03 0 10 9.105649e+03 1 3.237562e-02 12 1.109835e+04 1 1.140682e-02 14 1.309313e+04 0 16 1.508922e+04 1 2.517834e-01 18 1.708617e+04 0 20 1.908372e+04 0 FTF + Binary cascade Pi- Fe 2 2.004864e+03 2 1.855239e+02 4 4.002434e+03 4 4.646540e-02 6 6.001623e+03 6 5.577090e+02 8 8.001217e+03 1 4.412120e-03 10 1.000097e+04 5 1.054503e-01 12 1.200081e+04 9 2.518046e-02 14 1.400070e+04 5 1.288119e-02 16 1.600061e+04 3 2.747157e-02 18 1.800054e+04 7 6.015416e-02 20 2.000049e+04 8 5.448328e+02 Energy conservation.

5. Energy response and hadronic models, p+Fe

6. Energy response and hadronic models , Pi+Fe

7. Where is the Truth? Pi+Xe, Plab <= 8 GeV/c Pprot>220 Mev/c, 20% of protons Must go to pions! 50% of Pi0 are not registered! Proton multiplicity distribution is wrong!

8. Conclusion The Bertini model predicts too many nucleons! Improved FTF and the binary cascade give promising results for energy fractions. Butthey produce too many Pi0! Proton multiplicitydistributions are wrong. It would be well to tune Bert, BIC and FTF to experimental data at Plab > 1 GeV/c As a temporal solution a Reggeon theory inspiredmodel of nuclear destruction can be installedfor FTF and QGS.

9. Reggeon theory inspired model of nuclear destruction Model of nuclear disintegration in high-energy nucleus nucleus interactions. K. Abdel-Waged, V.V. Uzhinsky Phys.Atom.Nucl.60:828-840,1997, Yad.Fiz.60:925-937,1997. Features of p C interactions at a momentum of 4.2-GeV/c versus the degree of centrality of a collisions between protons and carbon nuclei: Multiplicity of secondary particles. A.S. Galoian, E.N. Kladnitskaya, V.V. Uzhinsky, O.V. Rogachevsky, R. Togoo Phys.Atom.Nucl.66:836-846,2003, Yad.Fiz.66:868-878,2003. Features of p C interactions at a momentum of 4.2-GeV/c versus the degree of centrality of a collision between a proton and a carbon nucleus: Kinematical features of secondaries. A.S. Galoian, E.N. Kladnitskaya, O.V. Rogachevsky, R. Togoo, V.V. Uzhinsky Phys.Atom.Nucl.67:256-272,2004, Yad.Fiz.67:273-289,2004. Multiplicity fluctuations of particles produced in interactions of light nuclei with carbon nuclei at a momentum of 4.2-A-GeV/c and their theoretical interpretation. A.S. Galoian, E. Kladnitskaya, V. Uzhinsky JETP Lett.86:630-633,2008, Pisma Zh.Eksp.Teor.Fiz.86:718-721,2008.

10. Reggeon theory inspired model of nuclear destruction

11. Reggeon theory inspired model of nuclear destruction Now it is implemented in FTF. The code is created,last bugs are fixing.