Ignore:
Timestamp:
Dec 22, 2010, 3:52:27 PM (13 years ago)
Author:
garnier
Message:

geant4 tag 9.4

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/source/processes/hadronic/models/parton_string/diffraction/src/G4FTFParameters.cc

    r1340 r1347  
    2525//
    2626//
    27 // $Id: G4FTFParameters.cc,v 1.14 2010/09/20 15:50:46 vuzhinsk Exp $
    28 // GEANT4 tag $Name: geant4-09-03-ref-09 $
     27// $Id: G4FTFParameters.cc,v 1.15 2010/11/15 10:02:38 vuzhinsk Exp $
     28// GEANT4 tag $Name: geant4-09-04-ref-00 $
    2929//
    3030
    … …  
    4242//**********************************************************************************************
    4343
     44//G4FTFParameters::G4FTFParameters(const G4ParticleDefinition * particle,
     45//                                                   G4double   theA,
     46//                                                   G4double   theZ,
     47//                                                   G4double   s)
    4448G4FTFParameters::G4FTFParameters(const G4ParticleDefinition * particle,
    45                                                    G4double   theA,
    46                                                    G4double   theZ,
     49                                                   G4int   theA,
     50                                                   G4int   theZ,
    4751                                                   G4double   s)
    4852{
    … …  
    6266    G4double sqrLogPlab = LogPlab * LogPlab;
    6367
    64     G4int NumberOfTargetProtons  = (G4int) theZ;
    65     G4int NumberOfTargetNeutrons = (G4int) theA- (G4int) theZ;
     68    G4int NumberOfTargetProtons  = theZ;
     69    G4int NumberOfTargetNeutrons = theA-theZ;
     70//    G4int NumberOfTargetProtons  = (G4int) theZ;
     71//    G4int NumberOfTargetNeutrons = (G4int) theA- (G4int) theZ;
    6672    G4int NumberOfTargetNucleons = NumberOfTargetProtons + NumberOfTargetNeutrons;
    6773
    6874    G4double Xtotal, Xelastic;
    6975
    70     if( absPDGcode > 1000 )                        //------Projectile is baryon --------
     76    if( PDGcode > 1000 )                        //------Projectile is baryon --------
    7177      {       
    7278       G4double XtotPP = 48.0 +  0. *std::pow(Plab, 0.  ) + 0.522*sqrLogPlab - 4.51*LogPlab;
    … …  
    7581       G4double XelPP  = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
    7682       G4double XelPN  = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
     83
     84       Xtotal          = ( NumberOfTargetProtons  * XtotPP +
     85                           NumberOfTargetNeutrons * XtotPN  ) / NumberOfTargetNucleons;
     86       Xelastic        = ( NumberOfTargetProtons  * XelPP  +
     87                           NumberOfTargetNeutrons * XelPN   ) / NumberOfTargetNucleons;
     88      }
     89    else if( PDGcode < -1000 )                 //------Projectile is anti_baryon --------
     90      {       
     91       G4double XtotPP = 38.4 +  77.6*std::pow(Plab,-0.64) + 0.26*sqrLogPlab - 1.2*LogPlab;
     92       G4double XtotPN =  0.  + 133.6*std::pow(Plab,-0.70) + 1.22*sqrLogPlab +13.7*LogPlab;
     93
     94       G4double XelPP  = 10.2 + 52.7*std::pow(Plab,-1.16) + 0.125*sqrLogPlab - 1.28*LogPlab;
     95       G4double XelPN  = 36.5 +  0. *std::pow(Plab, 0.  ) + 0.   *sqrLogPlab -11.9 *LogPlab;
    7796
    7897       Xtotal          = ( NumberOfTargetProtons  * XtotPP +
    … …  
    184203//      Xtotal and Xelastic in mb
    185204
     205// For Pi- P interactions only!
     206if(std::abs(Plab-1.4) < 0.05) {Xtotal=3.500599e+01; Xelastic= 1.150032e+01;}
     207if(std::abs(Plab-1.5) < 0.05) {Xtotal=3.450591e+01; Xelastic= 1.050038e+01;}
     208if(std::abs(Plab-1.6) < 0.05) {Xtotal=3.430576e+01; Xelastic= 9.800433e+00;}
     209if(std::abs(Plab-1.7) < 0.05) {Xtotal=3.455560e+01; Xelastic= 9.300436e+00;}
     210if(std::abs(Plab-1.8) < 0.05) {Xtotal=3.480545e+01; Xelastic= 8.800438e+00;}
     211if(std::abs(Plab-2.0) < 0.05) {Xtotal=3.570503e+01; Xelastic= 8.200370e+00;}
     212if(std::abs(Plab-2.2) < 0.05) {Xtotal=3.530495e+01; Xelastic= 7.800362e+00;}
     213if(std::abs(Plab-2.5) < 0.05) {Xtotal=3.410484e+01; Xelastic= 7.350320e+00;}
     214if(std::abs(Plab-2.75) < 0.05){Xtotal=3.280479e+01; Xelastic= 7.050273e+00;}
     215if(std::abs(Plab-3.0) < 0.05) {Xtotal=3.180473e+01; Xelastic= 6.800258e+00;}
     216if(std::abs(Plab-4.0) < 0.05) {Xtotal=2.910441e+01; Xelastic= 6.100229e+00;}
     217if(std::abs(Plab-5.0) < 0.05) {Xtotal=2.820372e+01; Xelastic= 5.700275e+00;}
     218if(std::abs(Plab-6.0) < 0.05) {Xtotal=2.760367e+01; Xelastic= 5.400255e+00;}
     219if(std::abs(Plab-7.0) < 0.05) {Xtotal=2.725366e+01; Xelastic= 5.150256e+00;}
     220if(std::abs(Plab-8.0) < 0.05) {Xtotal=2.690365e+01; Xelastic= 4.900258e+00;}
     221if(std::abs(Plab-10.0) < 0.05){Xtotal=2.660342e+01; Xelastic= 4.600237e+00;}
     222if(std::abs(Plab-12.0) < 0.05){Xtotal=2.632341e+01; Xelastic= 4.480229e+00;}
     223if(std::abs(Plab-14.0) < 0.05){Xtotal=2.604340e+01; Xelastic= 4.360221e+00;}
     224if(std::abs(Plab-20.0) < 0.05){Xtotal=2.520337e+01; Xelastic= 4.000197e+00;}
     225if(std::abs(Plab-30.0) < 0.05){Xtotal=2.505334e+01; Xelastic= 3.912679e+00;}
     226//
    186227//----------- Geometrical parameters ------------------------------------------------
    187228      SetTotalCrossSection(Xtotal);
    … …  
    189230      SetInelasticCrossSection(Xtotal-Xelastic);
    190231
    191 //G4cout<<"Xtotal, Xelastic "<<Xtotal<<" "<<Xelastic<<G4endl;
     232//G4cout<<"Plab Xtotal, Xelastic Xinel "<<Plab<<" "<<Xtotal<<" "<<Xelastic<<Xtotal-Xelastic)<<G4endl;
    192233//  // Interactions with elastic and inelastic collisions
    193234      SetProbabilityOfElasticScatt(Xtotal, Xelastic);
    … …  
    212253
    213254//----------- Parameters of excitations ---------------------------------------------
    214            if( absPDGcode > 1000 )                        //------Projectile is baryon --------
     255           if( PDGcode > 1000 )                        //------Projectile is baryon --------
    215256             {
    216257              SetMagQuarkExchange(1.84);//(3.63);
    217258              SetSlopeQuarkExchange(0.7);//(1.2);
    218259              SetDeltaProbAtQuarkExchange(0.);
     260              if(NumberOfTargetNucleons > 26) {SetProbOfSameQuarkExchange(1.);}
     261              else                            {SetProbOfSameQuarkExchange(0.);}
     262
     263              SetProjMinDiffMass(1.16);                   // GeV
     264              SetProjMinNonDiffMass(1.16);                // GeV
     265              SetProbabilityOfProjDiff(0.805*std::exp(-0.35*Ylab));// 0.5
     266
     267              SetTarMinDiffMass(1.16);                    // GeV
     268              SetTarMinNonDiffMass(1.16);                 // GeV
     269              SetProbabilityOfTarDiff(0.805*std::exp(-0.35*Ylab));// 0.5
     270
     271              SetAveragePt2(0.15);                        // 0.15 GeV^2
     272             }
     273           if( PDGcode < -1000 )                  //------Projectile is anti_baryon --------
     274             {
     275              SetMagQuarkExchange(0.);
     276              SetSlopeQuarkExchange(0.);
     277              SetDeltaProbAtQuarkExchange(0.);
     278              SetProbOfSameQuarkExchange(0.);
    219279
    220280              SetProjMinDiffMass(1.16);                   // GeV
    221281              SetProjMinNonDiffMass(1.16);                // GeV
    222              
    223 SetProbabilityOfProjDiff(0.805*std::exp(-0.35*Ylab));// 0.5
     282              SetProbabilityOfProjDiff(0.805*std::exp(-0.35*Ylab));// 0.5
    224283
    225284              SetTarMinDiffMass(1.16);                    // GeV
    226285              SetTarMinNonDiffMass(1.16);                 // GeV
    227              
    228 SetProbabilityOfTarDiff(0.805*std::exp(-0.35*Ylab));// 0.5
     286              SetProbabilityOfTarDiff(0.805*std::exp(-0.35*Ylab));// 0.5
    229287
    230288              SetAveragePt2(0.15);                        // 0.15 GeV^2
    … …  
    232290           else if( absPDGcode == 211 || PDGcode ==  111) //------Projectile is Pion -----------
    233291             {
    234               SetMagQuarkExchange(120.); // 210.
    235               SetSlopeQuarkExchange(2.0);
    236               SetDeltaProbAtQuarkExchange(0.6);
     292              SetMagQuarkExchange(240.);
     293              SetSlopeQuarkExchange(2.);
     294              SetDeltaProbAtQuarkExchange(0.56); //(0.35);
    237295
    238296              SetProjMinDiffMass(0.5);                    // GeV
    239               SetProjMinNonDiffMass(0.3);                 // GeV
    240               SetProbabilityOfProjDiff(0.*0.62*std::pow(s/GeV/GeV,-0.51)); // 40/32 X-dif/X-inel
    241 
    242               SetTarMinDiffMass(1.1);                     // GeV
    243               SetTarMinNonDiffMass(1.1);                  // GeV
    244 
    245               SetProbabilityOfTarDiff(2.*0.62*std::pow(s/GeV/GeV,-0.51)); // 40/32 X-dif/X-inel
     297              SetProjMinNonDiffMass(0.5);                 // GeV 0.3
     298              SetProbabilityOfProjDiff(0.);//(0.*0.62*std::pow(s/GeV/GeV,-0.51)); // 40/32 X-dif/X-inel
     299
     300              SetTarMinDiffMass(1.16);                     // GeV
     301              SetTarMinNonDiffMass(1.16);                  // GeV
     302//              SetProbabilityOfTarDiff(1.);//(2.*0.62*std::pow(s/GeV/GeV,-0.51));
     303//              SetProbabilityOfTarDiff(2.6*std::exp(-0.46*Ylab));
     304              SetProbabilityOfTarDiff(0.8*std::exp(-0.6*(Ylab-3.)));
    246305
    247306              SetAveragePt2(0.3);                         // GeV^2
    … …  
    295354      SetMaxNumberOfCollisions(1000.,1.); //(Plab,2.); //3.); ##############################
    296355
    297       SetCofNuclearDestruction(0.3); //1.0);           // for meson projectile
    298 
    299       SetDofNuclearDestruction(0.4);
    300       SetPt2ofNuclearDestruction(0.17*GeV*GeV);
    301       SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV);
    302 
    303       SetExcitationEnergyPerWoundedNucleon(100*MeV);
    304     } else                                             // for baryon projectile
    305     {
    306 //      SetMaxNumberOfCollisions(Plab,0.1); //6.); // ##############################
    307       SetMaxNumberOfCollisions(Plab,4.); //6.); // ##############################
    308 
    309       SetCofNuclearDestruction(0.62*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1))));
    310 
     356//      SetCofNuclearDestruction(0.); //1.0);           // for meson projectile
     357//      SetCofNuclearDestruction(1.*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1))));
     358//G4cout<<Ylab<<" "<<0.62*std::exp(4.*(Ylab-4.5))/(1.+std::exp(4.*(Ylab-4.5)))<<G4endl;
     359//G4int Uzhi; G4cin>>Uzhi;
     360
     361//      SetMaxNumberOfCollisions(Plab,2.); //4.); // ##############################
     362
     363      SetCofNuclearDestruction(1.*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1)))); //0.62 1.0
     364//------------------------------------------
     365//      SetDofNuclearDestruction(0.4);
     366//      SetPt2ofNuclearDestruction(0.17*GeV*GeV);
     367//      SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV);
     368
     369//      SetExcitationEnergyPerWoundedNucleon(100*MeV);
    311370      SetDofNuclearDestruction(0.4);
    312371      SetPt2ofNuclearDestruction((0.035+
    … …  
    315374
    316375      SetExcitationEnergyPerWoundedNucleon(75.*MeV);
     376    } else                                             // for baryon projectile
     377    {
     378      SetMaxNumberOfCollisions(Plab,2.); //4.); // ##############################
     379
     380      SetCofNuclearDestruction(1.*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1)))); //0.62 1.0
     381//G4cout<<Ylab<<" "<<0.62*std::exp(4.*(Ylab-2.1))/(1.+std::exp(4.*(Ylab-2.1)))<<G4endl;
     382//G4int Uzhi; G4cin>>Uzhi;
     383
     384      SetDofNuclearDestruction(0.4);
     385      SetPt2ofNuclearDestruction((0.035+
     386          0.04*std::exp(4.*(Ylab-2.5))/(1.+std::exp(4.*(Ylab-2.5))))*GeV*GeV); //0.09
     387      SetMaxPt2ofNuclearDestruction(1.0*GeV*GeV);
     388
     389      SetExcitationEnergyPerWoundedNucleon(75.*MeV);
    317390    }
    318391
    … …  
    322395//SetPt2ofNuclearDestruction((0.035+0.1*std::exp(4.*(Ylab-3.))/(1.+std::exp(4.*(Ylab-3.))))*GeV*GeV);
    323396
     397//SetMagQuarkExchange(120.); // 210. PipP
     398//SetSlopeQuarkExchange(2.0);
     399//SetDeltaProbAtQuarkExchange(0.6);
     400//SetProjMinDiffMass(0.7);                    // GeV 1.1
     401//SetProjMinNonDiffMass(0.7);                 // GeV
     402//SetProbabilityOfProjDiff(0.85*std::pow(s/GeV/GeV,-0.5)); // 40/32 X-dif/X-inel
     403//SetTarMinDiffMass(1.1);                     // GeV
     404//SetTarMinNonDiffMass(1.1);                  // GeV
     405//SetProbabilityOfTarDiff(0.85*std::pow(s/GeV/GeV,-0.5)); // 40/32 X-dif/X-inel
     406//
     407//SetAveragePt2(0.3);                         // GeV^2
     408//------------------------------------
    324409//SetProbabilityOfElasticScatt(1.,1.); //(Xtotal, Xelastic);
    325410//SetProbabilityOfProjDiff(1.*0.62*std::pow(s/GeV/GeV,-0.51)); // 0->1
    … …  
    328413//SetAvaragePt2ofElasticScattering(0.);
    329414
    330 //SetCofNuclearDestruction(0.6); //(0.4);                 
    331 SetExcitationEnergyPerWoundedNucleon(75.*MeV); //(75.*MeV);
    332 //SetDofNuclearDestruction(0.6); //(0.4);                 
    333 //SetPt2ofNuclearDestruction(0.12*GeV*GeV); //(0.168*GeV*GeV);
     415//SetMaxNumberOfCollisions(4.*(Plab+0.01),Plab); //6.); // ##############################
     416//SetCofNuclearDestruction(0.2); //(0.4);                 
     417//SetExcitationEnergyPerWoundedNucleon(0.*MeV); //(75.*MeV);
     418//SetDofNuclearDestruction(0.4); //(0.4);                 
     419//SetPt2ofNuclearDestruction(0.1*GeV*GeV); //(0.168*GeV*GeV);
    334420
    335421}
Note: See TracChangeset for help on using the changeset viewer.