Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4FTFSettingDefaultHDP Class Reference

Public Member Functions

 G4FTFSettingDefaultHDP ()
 

Detailed Description

Definition at line 867 of file G4FTFParameters.cc.

Constructor & Destructor Documentation

◆ G4FTFSettingDefaultHDP()

G4FTFSettingDefaultHDP::G4FTFSettingDefaultHDP ( )
inline

Definition at line 872 of file G4FTFParameters.cc.

872 {
873
874 //==================================================================
875 //
876 // Cross sections for elementary processes
877 //
878 // these are for Inelastic interactions, i.e. Xinelastic=(Xtotal-Xelastix)>0.
879 // for elastic, all the A's & B's, Atop & Ymin are zeros
880 // general formula: Pp = A1*exp(B1*Y) + A2*exp(B2*Y) + A3
881 // but if Y<Ymin, then Pp=max(0.,Atop)
882 // for details, see also G4FTFParameters::GetProcProb( ProcN, y )
883 //
884 // baryons
885 //
886 /* JVY, Oct. 31, 2017: Per Alberto R. & Vladimir U., keep this group of parameters FIXED */
887 /* JVY, June 11, 2020: try to open up... */
888 // Process=0 --> Qexchg w/o excitation
889 //
890 HDP.SetDefault( "FTF_BARYON_PROC0_A1", 13.71 );
891 HDP.SetDefault( "FTF_BARYON_PROC0_B1", 1.75 );
892 HDP.SetDefault( "FTF_BARYON_PROC0_A2", -30.69 );
893 HDP.SetDefault( "FTF_BARYON_PROC0_B2", 3.0 );
894 HDP.SetDefault( "FTF_BARYON_PROC0_A3", 0.0 );
895 HDP.SetDefault( "FTF_BARYON_PROC0_ATOP", 1.0 );
896 HDP.SetDefault( "FTF_BARYON_PROC0_YMIN", 0.93 );
897 /* */
898 //
899 // Process=1 --> Qexchg w/excitation
900 //
901 HDP.SetDefault( "FTF_BARYON_PROC1_A1", 25. );
902 HDP.SetDefault( "FTF_BARYON_PROC1_B1", 1. );
903 HDP.SetDefault( "FTF_BARYON_PROC1_A2", -50.34 );
904 HDP.SetDefault( "FTF_BARYON_PROC1_B2", 1.5 );
905 HDP.SetDefault( "FTF_BARYON_PROC1_A3", 0. );
906 HDP.SetDefault( "FTF_BARYON_PROC1_ATOP", 0. );
907 HDP.SetDefault( "FTF_BARYON_PROC1_YMIN", 1.4 );
908 /* */
909 //
910 // NOTE: Process #2 & 3 are projectile & target diffraction
911 // they have more complex definition of A1 & A2
912 // (see example below)
913 // SetParams( 2, 6.0/Xinel, 0.0 ,-6.0/Xinel*16.28, 3.0 , 0.0, 0.0 , 0.93);// Projectile diffraction
914 // SetParams( 3, 6.0/Xinel, 0.0 ,-6.0/Xinel*16.28, 3.0 , 0.0, 0.0 , 0.93);// Target diffraction
915 //
916 // Also, for ( AbsProjectileBaryonNumber > 10 || NumberOfTargetNucleons > 10 )
917 // projectile and/or target diffraction (dissociation) may be switched ON/OFF
918 //
919 HDP.SetDefault( "FTF_BARYON_DIFF_DISSO_PROJ", false );
920 HDP.SetDefault( "FTF_BARYON_DIFF_DISSO_TGT", false ); // as in hadr-string-diff-V10-03-07
921 //
922 /* JVY, Oct. 31, 2017: Per Alberto R. & Vladimir U., keep this group of parameters FIXED */
923 /* JVY, June 11, 2020: try to open up... */
924 // Process=4 --> Qexchg w/additional multiplier in excitation
925 //
926 HDP.SetDefault( "FTF_BARYON_PROC4_A1", 0.6 );
927 HDP.SetDefault( "FTF_BARYON_PROC4_B1", 0. );
928 HDP.SetDefault( "FTF_BARYON_PROC4_A2", -1.2 );
929 HDP.SetDefault( "FTF_BARYON_PROC4_B2", 0.5 );
930 HDP.SetDefault( "FTF_BARYON_PROC4_A3", 0. );
931 HDP.SetDefault( "FTF_BARYON_PROC4_ATOP",0. );
932 HDP.SetDefault( "FTF_BARYON_PROC4_YMIN",1.4 );
933 /* */
934 //
935 // Parameters of participating hadron (baryon) excitation
936 //
937 HDP.SetDefault( "FTF_BARYON_DELTA_PROB_QEXCHG", 0. );
938 HDP.SetDefault( "FTF_BARYON_PROB_SAME_QEXCHG", 0. );
939 HDP.SetDefault( "FTF_BARYON_DIFF_M_PROJ", 1.16, 1.16, 3. ); // it's supposed to be in GeV but do NOT do (*CLHEP::GeV)
940 // because it'll be done in the G4FTFParameters::SetProjMinDiffMass
941 HDP.SetDefault( "FTF_BARYON_NONDIFF_M_PROJ", 1.16, 1.16, 3. ); // do NOT (*CLHEP::GeV) - same as above
942 HDP.SetDefault( "FTF_BARYON_DIFF_M_TGT", 1.16, 1.16, 3. ); // do NOT (*CLHEP::GeV) - same as above
943 HDP.SetDefault( "FTF_BARYON_NONDIFF_M_TGT", 1.16, 1.16, 3. ); // do NOT (*CLHEP::GeV) - same as above
944 HDP.SetDefault( "FTF_BARYON_AVRG_PT2", 0.3, 0.08, 1. ); // do NOT (*CLHEP::GeV*CLHEP::GeV)
945 //
946 // JVY, Oct. 6, 2017: Per Alberto R., keep these two settings fixed (for now)
947 //
948 // HDP.SetDefault( "FTF_BARYON_PROB_DISTR_PROJ", 0.3 );
949 // HDP.SetDefault( "FTF_BARYON_PROB_DISTR_TGT", 0.3 );
950
951
952 // pions
953 //
954 // JVY, Aug.8, 2018 --> Feb.14, 2019 --> June 25, 2019:
955 // Parameters of participating hadron (pions) excitation
956 //
957 /* JVY, June 25, 2019: For now, keep this group of parameters FIXED */
958 // Process=0 --> Qexchg w/o excitation
959 //
960 HDP.SetDefault( "FTF_PION_PROC0_A1", 150.0 );
961 HDP.SetDefault( "FTF_PION_PROC0_B1", 1.8 );
962 HDP.SetDefault( "FTF_PION_PROC0_A2",-247.3 );
963 HDP.SetDefault( "FTF_PION_PROC0_B2", 2.3 );
964 HDP.SetDefault( "FTF_PION_PROC0_A3", 0.0 );
965 HDP.SetDefault( "FTF_PION_PROC0_ATOP", 1.0 );
966 HDP.SetDefault( "FTF_PION_PROC0_YMIN", 2.3 );
967 //
968 // Process=1 --> Qexchg w/excitation
969 //
970 HDP.SetDefault( "FTF_PION_PROC1_A1", 5.77 );
971 HDP.SetDefault( "FTF_PION_PROC1_B1", 0.6 );
972 HDP.SetDefault( "FTF_PION_PROC1_A2", -5.77 );
973 HDP.SetDefault( "FTF_PION_PROC1_B2", 0.8 );
974 HDP.SetDefault( "FTF_PION_PROC1_A3", 0. );
975 HDP.SetDefault( "FTF_PION_PROC1_ATOP", 0. );
976 HDP.SetDefault( "FTF_PION_PROC1_YMIN", 0.0 );
977 /*
978 //
979 // NOTE: Process #2 & 3 are projectile & target diffraction
980 //
981 // Process=2 --> Projectile diffraction
982 //
983 // Q: Would it even make sense to make these configurable ?
984 // The following is hadrcoded:
985 // Projectile Baryon Number > 10 (AbsProjectileBaryonNumber > 10)
986 // ... which is "strange" because projectile is a pion !!!... so it's always OFF
987 // (see also lines 1007-1016)
988 //
989 HDP.SetDefault( "FTF_PION_PROC2_A1", 2.27 );
990 HDP.SetDefault( "FTF_PION_PROC2_B1", 0.5 );
991 HDP.SetDefault( "FTF_PION_PROC2_A2", -98052.0);
992 HDP.SetDefault( "FTF_PION_PROC2_B2", 4.0 );
993 HDP.SetDefault( "FTF_PION_PROC2_A3", 0. );
994 HDP.SetDefault( "FTF_PION_PROC2_ATOP", 0. );
995 HDP.SetDefault( "FTF_PION_PROC2_YMIN", 3.0 );
996 */
997 //
998 // Process=3 --> Target diffraction
999 //
1000 HDP.SetDefault( "FTF_PION_PROC3_A1", 7.0 );
1001 HDP.SetDefault( "FTF_PION_PROC3_B1", 0.9 );
1002 HDP.SetDefault( "FTF_PION_PROC3_A2", -85.28);
1003 HDP.SetDefault( "FTF_PION_PROC3_B2", 1.9 );
1004 HDP.SetDefault( "FTF_PION_PROC3_A3", 0.08);
1005 HDP.SetDefault( "FTF_PION_PROC3_ATOP", 0. );
1006 HDP.SetDefault( "FTF_PION_PROC3_YMIN", 2.2 );
1007 //
1008 // projectile and/or target diffraction (dissociation) may be switched ON/OFF
1009 //
1010 // NOTE: Both projectile and target diffraction are turned OFF if
1011 // a) Number of Target Nucleons > 10 (NumberOfTargetNucleons > 10)
1012 // OR
1013 // b) Projectile Baryon Number > 10 (AbsProjectileBaryonNumber > 10)
1014 // ... which is "strange" because projectile is a pion !!!... so it's always OFF
1015 //
1016 HDP.SetDefault( "FTF_PION_DIFF_DISSO_PROJ", false );
1017 HDP.SetDefault( "FTF_PION_DIFF_DISSO_TGT", false );
1018 //
1019 /* JVY, June 25, 2019: For now keep this group of parameters FIXED */
1020 /* JVY, June 11, 2020: try to open up... */
1021 // Process=4 --> Qexchg w/additional multiplier in excitation
1022 //
1023 HDP.SetDefault( "FTF_PION_PROC4_A1", 1.0 );
1024 HDP.SetDefault( "FTF_PION_PROC4_B1", 0. );
1025 HDP.SetDefault( "FTF_PION_PROC4_A2",-11.02);
1026 HDP.SetDefault( "FTF_PION_PROC4_B2", 1.0 );
1027 HDP.SetDefault( "FTF_PION_PROC4_A3", 0. );
1028 HDP.SetDefault( "FTF_PION_PROC4_ATOP",0. );
1029 HDP.SetDefault( "FTF_PION_PROC4_YMIN",2.4 );
1030 /* */
1031 //
1032 // NOTE; As of geant4-10-05, all these settings beloe are correct
1033 // (and are the same as they were in 10.4.ref06)
1034
1035 HDP.SetDefault( "FTF_PION_DELTA_PROB_QEXCHG", 0.56 ); // in the past used to be 0.35
1036 HDP.SetDefault( "FTF_PION_DIFF_M_PROJ", 1.0, 0.5, 3. ); // in the past used to be 0.5...
1037 HDP.SetDefault( "FTF_PION_NONDIFF_M_PROJ", 1.0, 0.5, 3. ); // so why not set 0.5 as a low limit ?...
1038 // ... or perhaps even lower ?
1039 HDP.SetDefault( "FTF_PION_DIFF_M_TGT", 1.16, 1.16, 3. ); // All (NON)DIFF_M's are supposed to be in GeV but do NOT do (*CLHEP::GeV)
1040 // because it'll be done in the G4FTFParameters::SetProjMinDiffMass
1041 HDP.SetDefault( "FTF_PION_NONDIFF_M_TGT", 1.16, 1.16, 3. );
1042 HDP.SetDefault( "FTF_PION_AVRG_PT2", 0.3, 0.08, 1. ); // do NOT (*CLHEP::GeV*CLHEP::GeV)
1043
1044 //==================================================================
1045 //
1046 // nuclear destruction
1047 //
1048 // NOTE: Settings of most of these parameters are the same
1049 // for different types of projectile hadron
1050 // However, we decided to introduce separate variables
1051 // and configuration cards for each type of projectile
1052 //
1053 // baryons
1054 //
1055 // projectile destruction
1056 //
1057 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_PROJ", 1., 0., 1. ); // in principle, it should be 1./NBRN - FIXME later !
1058 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_NBRN_PROJ", false );
1059 //
1060 // for now, keep fixed p2 & p3 for the proj destruction
1061 // they're defined explicitly in G4FTFParamCollection ctor
1062 //
1063 // target destruction
1064 //
1065 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_TGT", 1., 0., 1. );
1066 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_ADEP_TGT", false );
1067 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P2_TGT", 4.0, 2., 16. );
1068 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P3_TGT", 2.1, 0., 4. );
1069 //
1070 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P1", 0.035, 0., 0.25 );
1071 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P2", 0.04, 0., 0.25 );
1072 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P3", 4.0, 2., 16. );
1073 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P4", 2.5, 0., 4. );
1074 //
1075 HDP.SetDefault( "FTF_BARYON_NUCDESTR_R2", 1.5*CLHEP::fermi*CLHEP::fermi, 0.5*CLHEP::fermi*CLHEP::fermi, 2.*CLHEP::fermi*CLHEP::fermi );
1076 HDP.SetDefault( "FTF_BARYON_EXCI_E_PER_WNDNUCLN", 40.*CLHEP::MeV, 0., 100.*CLHEP::MeV );
1077 HDP.SetDefault( "FTF_BARYON_NUCDESTR_DISP", 0.3, 0.1, 0.4 );
1078 //
1079 // JVY, Oct. 6, 2017: Per Alberto R., this is just a technical parameter,
1080 // and it should NOT be changed
1081 //
1082 // HDP.SetDefault( "FTF_BARYON_NUCDESTR_MAXPT2", 1. * CLHEP::GeV*CLHEP::GeV );
1083
1084
1085 // mesons - these parameters are common for pions, kaons, etc. (per original code)
1086 //
1087 // NOTE: *NO* projectile destruction for mesons !!!
1088 //
1089 // target destruction
1090 //
1091 HDP.SetDefault( "FTF_MESON_NUCDESTR_P1_TGT", 0.00481, 0., 1. );
1092 HDP.SetDefault( "FTF_MESON_NUCDESTR_P1_ADEP_TGT", true );
1093 HDP.SetDefault( "FTF_MESON_NUCDESTR_P2_TGT", 4.0, 2., 16. );
1094 HDP.SetDefault( "FTF_MESON_NUCDESTR_P3_TGT", 2.1, 0., 4. );
1095 //
1096 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P1", 0.035, 0., 0.25 );
1097 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P2", 0.04, 0., 0.25 );
1098 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P3", 4.0, 2., 16. );
1099 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P4", 2.5, 0., 4. );
1100 //
1101 HDP.SetDefault( "FTF_MESON_NUCDESTR_R2", 1.5*CLHEP::fermi*CLHEP::fermi,
1102 0.5*CLHEP::fermi*CLHEP::fermi,
1103 2.*CLHEP::fermi*CLHEP::fermi );
1104 HDP.SetDefault( "FTF_MESON_EXCI_E_PER_WNDNUCLN", 40.*CLHEP::MeV, 0., 100.*CLHEP::MeV );
1105 HDP.SetDefault( "FTF_MESON_NUCDESTR_DISP", 0.3, 0.1, 0.4 );
1106
1107 }
HDP
G4HadronicDeveloperParameters & HDP
Definition: G4FTFParameters.cc:864
G4HadronicDeveloperParameters::SetDefault
G4bool SetDefault(const std::string name, const G4bool value)
Definition: G4HadronicDeveloperParameters.cc:39
The documentation for this class was generated from the following file: