Changeset 1340 for trunk/source/processes/hadronic/models/parton_string/hadronization/src/G4VLongitudinalStringDecay.cc

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/src/G4VLongitudinalStringDecay.cc (modified) (6 diffs)

trunk/source/processes/hadronic/models/parton_string/hadronization/src/G4VLongitudinalStringDecay.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4VLongitudinalStringDecay.cc,v 1.20 2010/06/21 17:50:48 vuzhinsk Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4VLongitudinalStringDecay.cc,v 1.22 2010/09/20 12:46:23 vuzhinsk Exp $
+// GEANT4 tag $Name: geant4-09-03-ref-09 $
 //
 // -----------------------------------------------------------------------------
@@ -71 +71 @@
 
 // Changable Parameters below.
-   SigmaQT = 0.5 * GeV;
+   SigmaQT = 0.5 * GeV;  // 0.5
    
    StrangeSuppress  = 0.44;    //  27 % strange quarks produced, ie. u:d:s=1:1:0.27
@@ -250 +250 @@
            Hadron1 = (minMassHadronizer->*build)(string->GetLeftParton(),
                                                  string->GetRightParton());
+//G4cout<<"Hadron1 "<<Hadron1->GetParticleName()<<G4endl;
            mass= (Hadron1)->GetPDGMass();
         } else
@@ -323 +324 @@
                         G4FragmentingString *&newString)
 {
+//G4cout<<"Start SplitUP"<<G4endl;
        //... random choice of string end to use for creating the hadron (decay)   
        SideOfDecay = (G4UniformRand() < 0.5)? 1: -1;
@@ -341 +343 @@
            HadronDefinition= DiQuarkSplitup(string->GetDecayParton(), newStringEnd);
        }      
+
+//G4cout<<"New had "<<HadronDefinition->GetParticleName()<<G4endl;
 // create new String from old, ie. keep Left and Right order, but replace decay
 
        newString=new G4FragmentingString(*string,newStringEnd); // To store possible
                                                                 // quark containt of new string
+//G4cout<<"SplitEandP "<<G4endl;
        G4LorentzVector* HadronMomentum=SplitEandP(HadronDefinition, string, newString);
 
@@ -360 +365 @@
            delete HadronMomentum;
        }      
+//G4cout<<"End SplitUP"<<G4endl;
        return Hadron;
 }