Changeset 1340 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4Quasmon.cc

Timestamp:

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

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4Quasmon.cc (modified) (7 diffs)

trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4Quasmon.cc

@@ -28 +28 @@
 //
 //
-// $Id: G4Quasmon.cc,v 1.125 2010/06/10 08:37:27 mkossov Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4Quasmon.cc,v 1.127 2010/06/25 14:03:44 mkossov Exp $
+// GEANT4 tag $Name: hadr-chips-V09-03-08 $
 //
 //      ---------------- G4Quasmon ----------------
@@ -2382 +2382 @@
           G4cerr<<"G4Q::HQ:EV-6:TotEVAPORATION:s="<<sPDG<<",T="<<kinE<<",RM="<<retN4Mom.m()
                 <<"<"<<tmpTM<<",tQC="<<transQC<<",E="<<excE<<",sM="<<sumM<<">tM="<<totMass
-                <<",sCB="<<sCB<<",nQ="<<nQuasms<<G4endl;
+                <<",nQ="<<nQuasms<<G4endl;
           throw G4QException("G4Quasmon::HadrQuasm: Why Fail?(6)ProductMasses>totalMass");
 #endif
@@ -4448 +4448 @@
 // Randomize the Resonance masses and calculate probabilities of hadronization for them
 void G4Quasmon::CalculateHadronizationProbabilities
-  (G4double E, G4double kVal, G4LorentzVector k4M,G4bool piF, G4bool gaF, G4bool first)
+  (G4double E, G4double kVal, G4LorentzVector k4M,G4bool piF, G4bool gaF, G4bool )
+  //(G4double E, G4double kVal, G4LorentzVector k4M,G4bool piF, G4bool gaF, G4bool first)
 //   =====================================================================E is not used====
 {                                                                  //       ^
@@ -4529 +4530 @@
       //G4bool pos=curCand->GetPossibility();
 #ifdef pdebug
+      G4bool pPrint= (abs(cPDG)%10 <3 && cPDG <80000000) || (cPDG >80000000 && frM <5000.);
       //G4bool pPrint = abs(cPDG)%10<3 && cPDG<80000000 ||cPDG==90001000||cPDG==90000001||
       //  cPDG==90000002||cPDG==90001001||cPDG==90001002||cPDG==90002001||cPDG==90002002;
       //G4bool pPrint = cPDG==2212 || cPDG==2112 ||cPDG==90001000||cPDG==90000001;
-      G4bool pPrint = false;
+      //G4bool pPrint = false;
+      //G4bool pPrint = true;
       if(pPrint) G4cout<<"G4Q::CHP:==****==>>>c="<<cPDG<<",dUD="<<dUD<<",pos="<<pos<<",eA="
                        <<envA<<",tM="<<totMass<<" > tmpTM+frM="<<tmpTM+frM<<G4endl;
@@ -4715 +4718 @@
                     //if ( (!piF && first && baryn < 3) || 
                     //     (!piF && !first) || 
-                    if ( (!piF && first && baryn < 3) || 
-                         (!piF && !first && baryn < 5 ) || 
+                    /////if ( (!piF && first && baryn < 3) || 
+                    /////     (!piF && !first && baryn < 5 ) || 
+                    if ( (!piF && baryn < 5 ) || 
                          ( piF && abs(dS) < 3) ) // Isotope Focusing for AtRest Reactions
                     //if(!qIso&&!dC||qIso>0&&dC<0||qIso<0&&dC>0)//MediumIsoFocusingForAll
@@ -4889 +4893 @@
                         // ***VTN*** CHECK IsNecessety toRecover theColTotRes to MinMassTot
                         if(mintM2>rtQ2) //==> Check of ResidualTotalNucleus ** Always **
-                        //if(2>3)  // Negligable difference difference
+                        //if(2>3)  // Negligable difference
                         {
                           G4double nz=0.;
@@ -5071 +5075 @@
                         G4double lp=1.-(dkLS+sr+sr)/boundM;//qm=k+sqrt((E*(M-m)/M)^2-m^2)
 #ifdef debug
-                        if(pPrint) G4cout<<"G4Q::CHP:qa_k+sr="<<lp<<",sr="<<sr<<",m="
-                                         <<mex<<",M="<<frM<<G4endl;
+                        if(pPrint) G4cout<<"G4Q::CHP:qa_k+sr="<<lp<<",sr="<<sr
+                                         <<",M="<<frM<<G4endl;
 #endif
                         if(lp>0.&&lp<1.&&lp>lz)