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)

---

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 G4cerr, G4cout, and G4endl.

{
  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

---