Changeset 962 for trunk/source/processes/hadronic/models/parton_string/hadronization/src/G4FragmentingString.cc

Timestamp:

Apr 6, 2009, 12:30:29 PM (16 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/parton_string/hadronization/src/G4FragmentingString.cc (modified) (5 diffs)

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

@@ -80 +80 @@
         if ( old.decaying == Left )
         {
+//G4cout<<" Left "<<G4endl;
+//G4cout<<"Pt right "<<Ptright<<G4endl;
+//G4cout<<"Pt left  "<<Ptleft <<G4endl;
                 RightParton= old.RightParton;
                 Ptright    = old.Ptright;
@@ -85 +88 @@
                 Ptleft     = old.Ptleft - momentum->vect();
                 Ptleft.setZ(0.);
+//G4cout<<"Pt right "<<Ptright<<G4endl;
+//G4cout<<"Pt left  "<<Ptleft <<G4endl;
         } else if ( old.decaying == Right )
         {
+//G4cout<<" Right "<<G4endl;
+//G4cout<<"Pt right "<<Ptright<<G4endl;
+//G4cout<<"Pt left  "<<Ptleft <<G4endl;
                 RightParton = newdecay;
                 Ptright     = old.Ptright - momentum->vect();
@@ -92 +100 @@
                 LeftParton  = old.LeftParton;
                 Ptleft      = old.Ptleft;
+//G4cout<<"Pt right "<<Ptright<<G4endl;
+//G4cout<<"Pt left  "<<Ptleft <<G4endl;
         } else
         {
@@ -106 +116 @@
 //---------------------------------------------------------------------------------
 
+G4FragmentingString::G4FragmentingString(const G4FragmentingString &old,  // Uzhi
+                                         G4ParticleDefinition * newdecay) // Uzhi
+{                                                                         // Uzhi
+        decaying=None;                                                    // Uzhi
+        if ( old.decaying == Left )                                       // Uzhi
+        {                                                                 // Uzhi
+                RightParton= old.RightParton;                             // Uzhi
+                LeftParton = newdecay;                                    // Uzhi
+        } else if ( old.decaying == Right )                               // Uzhi
+        {                                                                 // Uzhi
+                RightParton = newdecay;                                   // Uzhi
+                LeftParton  = old.LeftParton;                             // Uzhi
+        } else                                                            // Uzhi
+        {
+                throw G4HadronicException(__FILE__, __LINE__, "G4FragmentingString::G4FragmentingString: no decay Direction defined");
+        }
+}
+
+
+//---------------------------------------------------------------------------------
+
 G4FragmentingString::~G4FragmentingString()
 {}
@@ -215 +246 @@
         return std::sqrt(this->Mass2());
 }
+
+G4double G4FragmentingString::MassT2() const
+{
+        return Pplus*Pminus;
+}