G4ExcitedStringDecay Class Reference

#include <G4ExcitedStringDecay.hh>

Inheritance diagram for G4ExcitedStringDecay:

Public Member Functions
	G4ExcitedStringDecay ()
	G4ExcitedStringDecay (G4VLongitudinalStringDecay *aStringDecay)
virtual	~G4ExcitedStringDecay ()
virtual G4KineticTrackVector *	FragmentStrings (const G4ExcitedStringVector *theStrings)
Public Member Functions inherited from G4VStringFragmentation
	G4VStringFragmentation ()
virtual	~G4VStringFragmentation ()

Detailed Description

Definition at line 38 of file G4ExcitedStringDecay.hh.

Constructor & Destructor Documentation

G4ExcitedStringDecay::G4ExcitedStringDecay

(

)

Definition at line 32 of file G4ExcitedStringDecay.cc.

                                           : G4VStringFragmentation(),
    theStringDecay(0)
{}

G4ExcitedStringDecay::G4ExcitedStringDecay

(

G4VLongitudinalStringDecay *

aStringDecay

)

Definition at line 36 of file G4ExcitedStringDecay.cc.

: G4VStringFragmentation(),
  theStringDecay(aStringDecay)
{}

G4ExcitedStringDecay::~G4ExcitedStringDecay ( )

virtual

Definition at line 48 of file G4ExcitedStringDecay.cc.

{
}

Member Function Documentation

G4KineticTrackVector * G4ExcitedStringDecay::FragmentStrings ( const G4ExcitedStringVector *  theStrings )

virtual

Implements G4VStringFragmentation.

Definition at line 79 of file G4ExcitedStringDecay.cc.

References CLHEP::HepLorentzVector::e(), G4cerr, G4cout, G4endl, and python.hepunit::perMillion.

{
  G4LorentzVector KTsum(0.,0.,0.,0.);

//G4cout<<"theStrings->size() "<<theStrings->size()<<G4endl;
  for ( unsigned int astring=0; astring < theStrings->size(); astring++)
  {
    KTsum+= theStrings->operator[](astring)->Get4Momentum();
  }

  G4KineticTrackVector * theResult = new G4KineticTrackVector;
  G4int attempts(0);
  G4bool success=false;
  G4bool NeedEnergyCorrector=false;
  do {
       //G4cout<<"Check of momentum at string fragmentations. New try."<<G4endl;
    std::for_each(theResult->begin() , theResult->end() , DeleteKineticTrack());
    theResult->clear();

    attempts++;
        //G4cout<<G4endl<<"attempts "<<attempts<<G4endl;
    G4LorentzVector KTsecondaries(0.,0.,0.,0.);
    NeedEnergyCorrector=false;

    for ( unsigned int astring=0; astring < theStrings->size(); astring++)
    {
          //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;
          //G4int Uzhi; G4cin >>Uzhi;
          G4KineticTrackVector * generatedKineticTracks = NULL;
      if ( theStrings->operator[](astring)->IsExcited() )
      {
           //G4cout<<"Fragment String"<<G4endl;
         generatedKineticTracks=FragmentString(*theStrings->operator[](astring));
      } else {
         generatedKineticTracks = new G4KineticTrackVector;
         generatedKineticTracks->push_back(theStrings->operator[](astring)->GetKineticTrack());
      }    

      if (generatedKineticTracks == NULL) 
      {
         G4cerr << "G4VPartonStringModel:No KineticTracks produced" << G4endl;
         continue;
      }

          G4LorentzVector KTsum1(0.,0.,0.,0.);
          for ( unsigned int aTrack=0; aTrack<generatedKineticTracks->size();aTrack++)
      {
              //--debug-- G4cout<<"Prod part "<<(*generatedKineticTracks)[aTrack]->GetDefinition()->GetParticleName()<<" "<<(*generatedKineticTracks)[aTrack]->Get4Momentum()<<G4endl;
             theResult->push_back(generatedKineticTracks->operator[](aTrack));
             KTsum1+= (*generatedKineticTracks)[aTrack]->Get4Momentum();
      }
      KTsecondaries+=KTsum1;
    
          //--debug--G4cout << "String secondaries(" <<generatedKineticTracks->size()<< ")  momentum: "
          //--debug--<< theStrings->operator[](astring)->Get4Momentum() << " " << KTsum1 << G4endl;
      if  ( KTsum1.e() > 0 && std::abs((KTsum1.e()-theStrings->operator[](astring)->Get4Momentum().e()) / KTsum1.e()) > perMillion )
      {
          //--debug--  G4cout << "String secondaries(" <<generatedKineticTracks->size()<< ")  momentum: "
          //--debug--    << theStrings->operator[](astring)->Get4Momentum() << " " << KTsum1 << G4endl;
        NeedEnergyCorrector=true;
      }

//        clean up
      delete generatedKineticTracks;
    }
       //--debug  G4cout << "Initial Strings / secondaries total  4 momentum " << KTsum << " " <<KTsecondaries << G4endl;

    success=true;
        //G4cout<<"success "<<success<<G4endl;
    if ( NeedEnergyCorrector ) success=EnergyAndMomentumCorrector(theResult, KTsum);
        //G4cout<<"success after Ecorr "<<success<<G4endl;
  } while(!success && (attempts < 10));   // It was 100 !!! Uzhi
          //G4cout<<"End frag string"<<G4endl;

#ifdef debug_ExcitedStringDecay
  G4LorentzVector  KTsum1=0;
  for ( unsigned int aTrack=0; aTrack<theResult->size();aTrack++)
  {
      G4cout << " corrected tracks .. " << (*theResult)[aTrack]->GetDefinition()->GetParticleName()
      <<"  " << (*theResult)[aTrack]->Get4Momentum() << G4endl;
      KTsum1+= (*theResult)[aTrack]->Get4Momentum();
  }

  G4cout << "Needcorrector/success " << NeedEnergyCorrector << "/" << success << ", Corrected total  4 momentum " << KTsum1  << G4endl;
  if ( ! success ) G4cout << "failed to correct E/p" << G4endl;  
#endif

  return theResult;
}

---

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

• G4ExcitedStringDecay.hh
• G4ExcitedStringDecay.cc