Changeset 1340 for trunk/source/processes/hadronic/models/parton_string/hadronization/include/G4ExcitedStringDecay.hh

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/parton_string/hadronization/include/G4ExcitedStringDecay.hh (modified) (3 diffs)

trunk/source/processes/hadronic/models/parton_string/hadronization/include/G4ExcitedStringDecay.hh

@@ -25 +25 @@
 //
 //
-// $Id: G4ExcitedStringDecay.hh,v 1.12 2010/06/21 17:50:48 vuzhinsk Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4ExcitedStringDecay.hh,v 1.13 2010/08/05 08:44:37 vuzhinsk Exp $
+// GEANT4 tag $Name: geant4-09-03-ref-09 $
 //
 #ifndef G4ExcitedStringDecay_h
@@ -42 +42 @@
       G4ExcitedStringDecay();
       G4ExcitedStringDecay(G4VLongitudinalStringDecay * aStringDecay);
-      ~G4ExcitedStringDecay();
+      virtual ~G4ExcitedStringDecay();
 
   private:
@@ -62 +62 @@
 };
 
-inline
-G4KineticTrackVector *G4ExcitedStringDecay::
-FragmentString(const G4ExcitedString &theString)
-{
-        if ( theStringDecay == NULL ) 
-            theStringDecay=new G4LundStringFragmentation();
-            
-        return theStringDecay->FragmentString(theString);
-}
-        
-
-inline
-G4KineticTrackVector *G4ExcitedStringDecay::
-FragmentStrings(const G4ExcitedStringVector * theStrings)
-{
-  G4KineticTrackVector * theResult = new G4KineticTrackVector;
-
-  G4LorentzVector KTsum(0.,0.,0.,0.);
-  G4LorentzVector KTsecondaries(0.,0.,0.,0.);
-  G4bool NeedEnergyCorrector=false;
-
-  for ( unsigned int astring=0; astring < theStrings->size(); astring++)
-  {
-        KTsum+= theStrings->operator[](astring)->Get4Momentum();
-
-        if( !(KTsum.e()<1) && !(KTsum.e()>-1) )
-        {
-          throw G4HadronicException(__FILE__, __LINE__, 
-                                   "G4ExcitedStringDecay::FragmentStrings received nan string...");
-        }
-        G4KineticTrackVector * generatedKineticTracks = NULL;
-  
-        if ( theStrings->operator[](astring)->IsExcited() )
-        {
-             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++)
-        {
-                theResult->push_back(generatedKineticTracks->operator[](aTrack));
-                KTsum1+= (*generatedKineticTracks)[aTrack]->Get4Momentum();
-        }
-        KTsecondaries+=KTsum1;
-        
-        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 << "Strings/secs total  4 momentum " << KTsum << " " <<KTsecondaries << G4endl;
-
-  G4bool success=true;
-  if ( NeedEnergyCorrector ) success=EnergyAndMomentumCorrector(theResult, KTsum);
-
-
-#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;
-}
-
 #endif