- Timestamp:
- Jan 8, 2010, 11:56:51 AM (15 years ago)
- Location:
- trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src
- Files:
-
- 20 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QBesIKJY.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QBesIKJY.cc,v 1.4 2009/11/10 17:13:46 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QBesIKJY ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QCHIPSWorld.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QCHIPSWorld.cc,v 1.33 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QCHIPSWorld ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QCandidate.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QCandidate.cc,v 1.35 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QCandidate ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QChipolino.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QChipolino.cc,v 1.36 2009/08/10 16:36:53 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QChipolino ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QContent.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QContent.cc,v 1.49 2009/08/07 14:20:57 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QContent ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QDecayChan.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QDecayChan.cc,v 1.28 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QDecayChan ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QEnvironment.cc
r1196 r1228 28 28 // 29 29 // 30 // $Id: G4QEnvironment.cc,v 1.1 57 2009/11/16 18:15:01mkossov Exp $31 // GEANT4 tag $Name: geant4-09-03 -cand-01$30 // $Id: G4QEnvironment.cc,v 1.160 2009/12/03 18:09:07 mkossov Exp $ 31 // GEANT4 tag $Name: geant4-09-03 $ 32 32 // 33 33 // ---------------- G4QEnvironment ---------------- … … 1737 1737 //G4int c3Max = 9; // Max#of "no hadrons" steps (reduced below?) 1738 1738 //G4int c3Max = 3; 1739 G4int c3Max = 1; 1739 //G4int c3Max = 1; 1740 G4int c3Max = 0; 1740 1741 //G4int premC = 27; 1741 1742 //G4int premC = 3; … … 1745 1746 //if(envA>1&&envA<19) premC = 36/envA; 1746 1747 //if(envA>1&&envA<25) premC = 48/envA; 1747 if(envA>1&&envA<41) premC = 80/envA;1748 //if(envA>1&&envA<41) premC = 80/envA; 1748 1749 G4int sumstat= 2; // Sum of statuses of all Quasmons 1749 1750 G4bool force = false; // Prototype of the Force Major Flag 1750 1751 G4int cbR =0; // Counter of the "Stoped by Coulomb Barrier" 1751 1752 // 1752 G4int cbRM =1; // MaxCounter of the "StopedByCoulombBarrier" 1753 G4int cbRM =0; //** MaxCounter of "StopedByCoulombBarrier" ** 1754 //G4int cbRM =1; // MaxCounter of the "StopedByCoulombBarrier" 1753 1755 //G4int cbRM =3; // MaxCounter of the "StopedByCoulombBarrier" 1754 1756 //G4int cbRM =9; // MaxCounter of the "Stoped byCoulombBarrier" … … 1762 1764 G4bool first=true; // Flag of the first interaction (only NucMedia) 1763 1765 G4int cAN=0; // Counter of the nucleon absorptions 1764 G4int mcAN=3; // Max for the counter of nucleon absorptions 1766 //G4int mcAN=27; // Max for the counter of nucleon absorptions 1767 G4int mcAN=9; // Max for the counter of nucleon absorptions 1768 //G4int mcAN=3; // Max for the counter of nucleon absorptions 1769 //G4int mcAN=1; // Max for the counter of nucleon absorptions 1765 1770 while (sumstat||totC<totCM) // ===***=== The MAIN "FOREVER" LOOP ===***=== 1766 1771 { -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QException.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QException.cc,v 1.10 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QException ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QHadron.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QHadron.cc,v 1.64 2009/09/02 15:45:19 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QHadron ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QInteraction.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QInteraction.cc,v 1.6 2009/08/28 14:49:10 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ------------------------------------------------------------------ -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QIsotope.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QIsotope.cc,v 1.13 2009/08/28 14:49:10 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QIsotope class ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QNucleus.cc
r1196 r1228 25 25 // 26 26 // 27 // $Id: G4QNucleus.cc,v 1.11 5 2009/11/04 10:52:34mkossov Exp $28 // GEANT4 tag $Name: geant4-09-03 -cand-01$27 // $Id: G4QNucleus.cc,v 1.116 2009/12/14 16:41:52 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QNucleus ---------------- … … 5990 5990 G4cerr<<"***G4QN::DecDiBar: badPDG="<<qPDG<<" or smallM="<<qM<<",2mP="<<dProt 5991 5991 <<",2mN="<<dNeut<<G4endl; 5992 throw G4QException("G4QNucleus::DecayDibar: Unknown PDG code or small Mass of DB"); 5992 // @@ Nothing to do. Just 2 GeV disappears... Very rare! Just to avoid the exception. 5993 //throw G4QException("G4QNucleus::DecayDibar: Unknown PDG code or small Mass of DB"); 5993 5994 } 5994 5995 } -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QPDGCode.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QPDGCode.cc,v 1.63 2009/11/03 16:13:37 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QPDGCode ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QParentCluster.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QParentCluster.cc,v 1.26 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QParentCluster ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QParticle.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QParticle.cc,v 1.34 2009/02/23 09:49:24 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ---------------- G4QParticle ---------------- -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QParton.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QParton.cc,v 1.9 2009/07/17 16:54:57 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ------------------------------------------------------------ -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QPartonPair.cc
r1196 r1228 27 27 // 28 28 // $Id: G4QPartonPair.cc,v 1.2 2009/02/23 09:49:24 mkossov Exp $ 29 // GEANT4 tag $Name: geant4-09-03 -cand-01$29 // GEANT4 tag $Name: geant4-09-03 $ 30 30 // 31 31 // ------------------------------------------------------------ -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QProbability.cc
r1197 r1228 26 26 // 27 27 // $Id: G4QProbability.cc,v 1.3 2009/09/04 16:13:19 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ------------------------------------------------------------ -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4QString.cc
r1196 r1228 26 26 // 27 27 // $Id: G4QString.cc,v 1.17 2009/09/04 14:38:00 mkossov Exp $ 28 // GEANT4 tag $Name: geant4-09-03 -cand-01$28 // GEANT4 tag $Name: geant4-09-03 $ 29 29 // 30 30 // ------------------------------------------------------------ -
trunk/source/processes/hadronic/models/chiral_inv_phase_space/body/src/G4Quasmon.cc
r1196 r1228 28 28 // 29 29 // 30 // $Id: G4Quasmon.cc,v 1.1 19 2009/11/16 18:15:01mkossov Exp $31 // GEANT4 tag $Name: geant4-09-03 -cand-01$30 // $Id: G4Quasmon.cc,v 1.121 2009/12/03 18:09:07 mkossov Exp $ 31 // GEANT4 tag $Name: geant4-09-03 $ 32 32 // 33 33 // ---------------- G4Quasmon ---------------- … … 537 537 } 538 538 G4ThreeVector totBoost = tot4M.boostVector(); // BoostVector for TotalSystem (backward) 539 G4ThreeVector totRBoost= -totBoost; 539 G4ThreeVector totRBoost= -totBoost; // Boost vector for Total System (forward) 540 540 G4int iniPDG =valQ.GetSPDGCode(); 541 541 G4int iniBN =valQ.GetBaryonNumber(); … … 548 548 #endif 549 549 G4QNucleus iniRN(iniP,iniN-1,iniS); 550 G4double iniRM = iniRN.GetMZNS(); 550 G4double iniRM = iniRN.GetMZNS(); // Mass of Residual Quasmon when neutron is rad 551 551 if(iniBN<2||envA>0) iniRM=0.; 552 552 G4double iniQM =G4QPDGCode(iniPDG).GetMass(); // Minimum mass of Quasmon … … 565 565 <<envQC<<",Q+E="<<quen<<",tM="<<totPDG<<totQC<<totM<<"<"<<totMass<<G4endl; 566 566 #endif 567 G4int tQ = valQ.GetTot(); 567 G4int tQ = valQ.GetTot(); // Total number of quarks for current Quasmon 568 568 G4int bQ = abs(valQ.GetBaryonNumber()); // Baryon number of the current Quasmon 569 G4QContent cQ = valQ; 570 G4int s = 4; 571 if (bQ) s = 3*bQ + 2; 572 if (tQ> s) cQ.DecQAQ((tQ-s)/2); 573 #ifdef debug 574 G4int rsPDG = cQ.GetSPDGCode(); 569 G4QContent cQ = valQ; // Temporary copy of Quasmon QC 570 G4int s = 4; // Mesonic 571 if (bQ) s = 3*bQ + 2; // Barionic 572 if (tQ> s) cQ.DecQAQ((tQ-s)/2); // Reduce QC to minimum QC 573 #ifdef debug 574 G4int rsPDG = cQ.GetSPDGCode(); // PDG for the lowest residual Quasmon state 575 575 G4cout<<"G4Q::HQ:eN="<<envN<<",eZ="<<envZ<<",Q="<<rsPDG<<cQ<<",piF="<<piF<<",gaF="<<gaF 576 576 <<G4endl; 577 577 #endif 578 theEnvironment.UpdateClusters(false); // New A-clusters are calculated578 theEnvironment.UpdateClusters(false); // New A-clusters are calculated 579 579 //theEnvironment.PrepareCandidates(theQCandidates,false,false);//Calc.PrePreprob's of C 580 580 //G4bool fF=piF||gaF; 581 581 // piF,piF or gaF,gaF is correct. Made to avoid theSpecificInteractionTimeClusterChoice 582 582 theEnvironment.PrepareCandidates(theQCandidates,piF,piF);// Calc.PrePreprob of Candid's 583 ModifyInMatterCandidates(); 584 G4double kMom = 0.; 585 G4double minK = 0.; 586 G4double maxK = quasM/2.; 587 G4double kLS = 0; 588 G4double cost = 0.; 589 G4bool kCond = true; 590 G4bool qCond = true; 591 G4bool pCond = true; 592 G4bool fskip=false; 593 G4bool fred =false; 594 //G4bool ffdc =true; 595 G4LorentzVector k4Mom=zeroLV; 596 G4LorentzVector cr4Mom=zeroLV; 597 G4int kCount =0; 583 ModifyInMatterCandidates(); // Calculate InMediaMasses of Cand. & Possibil. 584 G4double kMom = 0.; // Energy of primary qParton in Q-CMS 585 G4double minK = 0.; // k_min for randomization 586 G4double maxK = quasM/2.; // k_max for randomization 587 G4double kLS = 0; // Energy of primary qParton in LS 588 G4double cost = 0.; // Cos(theta) of k in QS InRespecTo Q direction 589 G4bool kCond = true; // k choice condition 590 G4bool qCond = true; // q choice condition 591 G4bool pCond = true; // Not adoptable parent cluster choice condition 592 G4bool fskip=false; // Flag to skip when sucked 593 G4bool fred =false; // Flag of Environment reduction 594 //G4bool ffdc =true; // Flag of not successful decay in fragment 595 G4LorentzVector k4Mom=zeroLV; // 4-momentum prototype for k 596 G4LorentzVector cr4Mom=zeroLV; // 4-momentum prototype for the ColResQuasmon 597 G4int kCount =0; // Counter of attempts of k for hadronization 598 598 //G4int kCountMax=27; 599 599 //G4int kCountMax=9; 600 //G4int kCountMax=3; 601 G4int kCountMax=1; 602 //G4int kCountMax=0; //@@ *** Close search for the minimum k ***600 //G4int kCountMax=3; // Try different k if they are below minK 601 G4int kCountMax=1; // "No reson to increase it" 602 //G4int kCountMax=0; //@@ *** Close search for the minimum k *** 603 603 // 604 G4int qCountMax=27; // Try different q to come over CoulBar or SepE 604 //G4int qCountMax=27; // Try different q to come over CoulBar or SepE 605 G4int qCountMax=12; // Try different q to come over CoulBar or SepE 605 606 //G4int qCountMax=9; // Try different q to come over CoulBar or SepE 606 607 //G4int qCountMax=3; // Try different q to come over CoulBar or SepE … … 612 613 G4int pCountMax=1; //Try differentHadrons(Parents) forBetterRecoil 613 614 //if(envA>0) pCountMax=3; 614 if(envA>0&&envA<31) 615 //if(envA>0&&envA<31) 616 if(envA>0&&envA<61) 615 617 { 616 pCountMax=60/envA; 617 kCountMax=pCountMax; 618 //pCountMax=60/envA; 619 pCountMax=120/envA; 620 //kCountMax=pCountMax; // See above ! 618 621 } 619 622 G4bool gintFlag=false; // Flag of gamma interaction with one quark
Note: See TracChangeset
for help on using the changeset viewer.