Geant4 11.1.1
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 198 of file G4FTFTunings.cc.

Constructor & Destructor Documentation

◆ G4FTFSettingDefaultHDP()

G4FTFSettingDefaultHDP::G4FTFSettingDefaultHDP ( )
inline

Definition at line 200 of file G4FTFTunings.cc.

200 {
201 // Cross sections for elementary processes
202 //
203 // these are for Inelastic interactions, i.e. Xinelastic=(Xtotal-Xelastix)>0.
204 // for elastic, all the A's & B's, Atop & Ymin are zeros
205 // general formula: Pp = A1*exp(B1*Y) + A2*exp(B2*Y) + A3
206 // but if Y<Ymin, then Pp=max(0.,Atop)
207 // for details, see also G4FTFParameters::GetProcProb( ProcN, y )
208 //
209 // baryons
210 /* JVY, Oct. 31, 2017: Per Alberto R. & Vladimir U., keep this group of parameters FIXED */
211 /* JVY, June 11, 2020: try to open up... */
212 // Process=0 --> Qexchg w/o excitation
213 HDP.SetDefault( "FTF_BARYON_PROC0_A1", 13.71 );
214 HDP.SetDefault( "FTF_BARYON_PROC0_B1", 1.75 );
215 HDP.SetDefault( "FTF_BARYON_PROC0_A2", -30.69 );
216 HDP.SetDefault( "FTF_BARYON_PROC0_B2", 3.0 );
217 HDP.SetDefault( "FTF_BARYON_PROC0_A3", 0.0 );
218 HDP.SetDefault( "FTF_BARYON_PROC0_ATOP", 1.0 );
219 HDP.SetDefault( "FTF_BARYON_PROC0_YMIN", 0.93 );
220 // Process=1 --> Qexchg w/excitation
221 HDP.SetDefault( "FTF_BARYON_PROC1_A1", 25.0 );
222 HDP.SetDefault( "FTF_BARYON_PROC1_B1", 1.0 );
223 HDP.SetDefault( "FTF_BARYON_PROC1_A2", -50.34 );
224 HDP.SetDefault( "FTF_BARYON_PROC1_B2", 1.5 );
225 HDP.SetDefault( "FTF_BARYON_PROC1_A3", 0.0 );
226 HDP.SetDefault( "FTF_BARYON_PROC1_ATOP", 0.0 );
227 HDP.SetDefault( "FTF_BARYON_PROC1_YMIN", 1.4 );
228 // NOTE: Process #2 & 3 are projectile & target diffraction
229 // they have more complex definition of A1 & A2
230 // (see example below)
231 // SetParams( 2, 6.0/Xinel, 0.0 ,-6.0/Xinel*16.28, 3.0 , 0.0, 0.0, 0.93 ); // Projectile diffraction
232 // SetParams( 3, 6.0/Xinel, 0.0 ,-6.0/Xinel*16.28, 3.0 , 0.0, 0.0, 0.93 ); // Target diffraction
233 // Also, for ( AbsProjectileBaryonNumber > 10 || NumberOfTargetNucleons > 10 )
234 // projectile and/or target diffraction (dissociation) may be switched ON/OFF
235 HDP.SetDefault( "FTF_BARYON_DIFF_DISSO_PROJ", false );
236 HDP.SetDefault( "FTF_BARYON_DIFF_DISSO_TGT", false );
237 /* JVY, Oct. 31, 2017: Per Alberto R. & Vladimir U., keep this group of parameters FIXED */
238 /* JVY, June 11, 2020: try to open up... */
239 // Process=4 --> Qexchg w/additional multiplier in excitation
240 HDP.SetDefault( "FTF_BARYON_PROC4_A1", 0.6 );
241 HDP.SetDefault( "FTF_BARYON_PROC4_B1", 0.0 );
242 HDP.SetDefault( "FTF_BARYON_PROC4_A2", -1.2 );
243 HDP.SetDefault( "FTF_BARYON_PROC4_B2", 0.5 );
244 HDP.SetDefault( "FTF_BARYON_PROC4_A3", 0.0 );
245 HDP.SetDefault( "FTF_BARYON_PROC4_ATOP", 0.0 );
246 HDP.SetDefault( "FTF_BARYON_PROC4_YMIN", 1.4 );
247 // Parameters of participating hadron (baryon) excitation
248 HDP.SetDefault( "FTF_BARYON_DELTA_PROB_QEXCHG", 0.0 );
249 HDP.SetDefault( "FTF_BARYON_PROB_SAME_QEXCHG", 0.0 );
250 HDP.SetDefault( "FTF_BARYON_DIFF_M_PROJ", 1.16, 1.16, 3.0 ); // it's supposed to be in GeV but do NOT do (*CLHEP::GeV)
251 // because it'll be done in the G4FTFParameters::SetProjMinDiffMass
252 HDP.SetDefault( "FTF_BARYON_NONDIFF_M_PROJ", 1.16, 1.16, 3.0 ); // do NOT (*CLHEP::GeV) - same as above
253 HDP.SetDefault( "FTF_BARYON_DIFF_M_TGT", 1.16, 1.16, 3.0 ); // do NOT (*CLHEP::GeV) - same as above
254 HDP.SetDefault( "FTF_BARYON_NONDIFF_M_TGT", 1.16, 1.16, 3.0 ); // do NOT (*CLHEP::GeV) - same as above
255 HDP.SetDefault( "FTF_BARYON_AVRG_PT2", 0.3, 0.08, 1.0 ); // do NOT (*CLHEP::GeV*CLHEP::GeV)
256 // JVY, Oct. 6, 2017: Per Alberto R., keep these two settings fixed (for now)
257 // HDP.SetDefault( "FTF_BARYON_PROB_DISTR_PROJ", 0.3 );
258 // HDP.SetDefault( "FTF_BARYON_PROB_DISTR_TGT", 0.3 );
259 // pions
260 // JVY, Aug.8, 2018 --> Feb.14, 2019 --> June 25, 2019:
261 // Parameters of participating hadron (pions) excitation
262 /* JVY, June 25, 2019: For now, keep this group of parameters FIXED */
263 // Process=0 --> Qexchg w/o excitation
264 HDP.SetDefault( "FTF_PION_PROC0_A1", 150.0 );
265 HDP.SetDefault( "FTF_PION_PROC0_B1", 1.8 );
266 HDP.SetDefault( "FTF_PION_PROC0_A2", -247.3 );
267 HDP.SetDefault( "FTF_PION_PROC0_B2", 2.3 );
268 HDP.SetDefault( "FTF_PION_PROC0_A3", 0.0 );
269 HDP.SetDefault( "FTF_PION_PROC0_ATOP", 1.0 );
270 HDP.SetDefault( "FTF_PION_PROC0_YMIN", 2.3 );
271 // Process=1 --> Qexchg w/excitation
272 HDP.SetDefault( "FTF_PION_PROC1_A1", 5.77 );
273 HDP.SetDefault( "FTF_PION_PROC1_B1", 0.6 );
274 HDP.SetDefault( "FTF_PION_PROC1_A2", -5.77 );
275 HDP.SetDefault( "FTF_PION_PROC1_B2", 0.8 );
276 HDP.SetDefault( "FTF_PION_PROC1_A3", 0.0 );
277 HDP.SetDefault( "FTF_PION_PROC1_ATOP", 0.0 );
278 HDP.SetDefault( "FTF_PION_PROC1_YMIN", 0.0 );
279 /*
280 // NOTE: Process #2 & 3 are projectile & target diffraction
281 // Process=2 --> Projectile diffraction
282 // Q: Would it even make sense to make these configurable ?
283 // The following is hadrcoded:
284 // Projectile Baryon Number > 10 (AbsProjectileBaryonNumber > 10)
285 // ... which is "strange" because projectile is a pion !!!... so it's always OFF
286 // (see also lines 1007-1016)
287 HDP.SetDefault( "FTF_PION_PROC2_A1", 2.27 );
288 HDP.SetDefault( "FTF_PION_PROC2_B1", 0.5 );
289 HDP.SetDefault( "FTF_PION_PROC2_A2", -98052.0 );
290 HDP.SetDefault( "FTF_PION_PROC2_B2", 4.0 );
291 HDP.SetDefault( "FTF_PION_PROC2_A3", 0.0 );
292 HDP.SetDefault( "FTF_PION_PROC2_ATOP", 0.0 );
293 HDP.SetDefault( "FTF_PION_PROC2_YMIN", 3.0 );
294 */
295 // Process=3 --> Target diffraction
296 HDP.SetDefault( "FTF_PION_PROC3_A1", 7.0 );
297 HDP.SetDefault( "FTF_PION_PROC3_B1", 0.9 );
298 HDP.SetDefault( "FTF_PION_PROC3_A2", -85.28 );
299 HDP.SetDefault( "FTF_PION_PROC3_B2", 1.9 );
300 HDP.SetDefault( "FTF_PION_PROC3_A3", 0.08 );
301 HDP.SetDefault( "FTF_PION_PROC3_ATOP", 0.0 );
302 HDP.SetDefault( "FTF_PION_PROC3_YMIN", 2.2 );
303 // projectile and/or target diffraction (dissociation) may be switched ON/OFF
304 // NOTE: Both projectile and target diffraction are turned OFF if
305 // a) Number of Target Nucleons > 10 (NumberOfTargetNucleons > 10)
306 // OR
307 // b) Projectile Baryon Number > 10 (AbsProjectileBaryonNumber > 10)
308 // ... which is "strange" because projectile is a pion !!!... so it's always OFF
309 HDP.SetDefault( "FTF_PION_DIFF_DISSO_PROJ", false );
310 HDP.SetDefault( "FTF_PION_DIFF_DISSO_TGT", false );
311 /* JVY, June 25, 2019: For now keep this group of parameters FIXED */
312 /* JVY, June 11, 2020: try to open up... */
313 // Process=4 --> Qexchg w/additional multiplier in excitation
314 HDP.SetDefault( "FTF_PION_PROC4_A1", 1.0 );
315 HDP.SetDefault( "FTF_PION_PROC4_B1", 0.0 );
316 HDP.SetDefault( "FTF_PION_PROC4_A2", -11.02 );
317 HDP.SetDefault( "FTF_PION_PROC4_B2", 1.0 );
318 HDP.SetDefault( "FTF_PION_PROC4_A3", 0.0 );
319 HDP.SetDefault( "FTF_PION_PROC4_ATOP", 0.0 );
320 HDP.SetDefault( "FTF_PION_PROC4_YMIN", 2.4 );
321 //
322 HDP.SetDefault( "FTF_PION_DELTA_PROB_QEXCHG", 0.56 );
323 HDP.SetDefault( "FTF_PION_DIFF_M_PROJ", 1.0, 0.5, 3.0 );
324 HDP.SetDefault( "FTF_PION_NONDIFF_M_PROJ", 1.0, 0.5, 3.0 );
325 HDP.SetDefault( "FTF_PION_DIFF_M_TGT", 1.16, 1.16, 3.0 ); // All (NON)DIFF_M's are supposed to be in GeV but do NOT do (*CLHEP::GeV)
326 // because it'll be done in the G4FTFParameters::SetProjMinDiffMass
327 HDP.SetDefault( "FTF_PION_NONDIFF_M_TGT", 1.16, 1.16, 3.0 );
328 HDP.SetDefault( "FTF_PION_AVRG_PT2", 0.3, 0.08, 1.0 ); // do NOT (*CLHEP::GeV*CLHEP::GeV)
329 // nuclear destruction
330 // NOTE: Settings of most of these parameters are the same
331 // for different types of projectile hadron
332 // However, we decided to introduce separate variables
333 // and configuration cards for each type of projectile
334 // baryons
335 // projectile destruction
336 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_PROJ", 1.0, 0.0, 1.0 ); // in principle, it should be 1./NBRN - FIXME later !
337 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_NBRN_PROJ", false );
338 // for now, keep fixed p2 & p3 for the proj destruction
339 // they're defined explicitly in G4FTFParamCollection ctor
340 // target destruction
341 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_TGT", 1.0, 0.0, 1.0 );
342 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P1_ADEP_TGT", false );
343 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P2_TGT", 4.0, 2.0, 16.0 );
344 HDP.SetDefault( "FTF_BARYON_NUCDESTR_P3_TGT", 2.1, 0.0, 4.0 );
345 //
346 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P1", 0.035, 0.0, 0.25 );
347 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P2", 0.04, 0.0, 0.25 );
348 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P3", 4.0, 2.0, 16.0 );
349 HDP.SetDefault( "FTF_BARYON_PT2_NUCDESTR_P4", 2.5, 0.0, 4.0 );
350 //
351 HDP.SetDefault( "FTF_BARYON_NUCDESTR_R2", 1.5*CLHEP::fermi*CLHEP::fermi, 0.5*CLHEP::fermi*CLHEP::fermi, 2.0*CLHEP::fermi*CLHEP::fermi );
352 HDP.SetDefault( "FTF_BARYON_EXCI_E_PER_WNDNUCLN", 40.0*CLHEP::MeV, 0.0, 100.0*CLHEP::MeV );
353 HDP.SetDefault( "FTF_BARYON_NUCDESTR_DISP", 0.3, 0.1, 0.4 );
354 // JVY, Oct. 6, 2017: Per Alberto R., this is just a technical parameter,
355 // and it should NOT be changed
356 // HDP.SetDefault( "FTF_BARYON_NUCDESTR_MAXPT2", 1. * CLHEP::GeV*CLHEP::GeV );
357 // mesons - these parameters are common for pions, kaons, etc. (per original code)
358 // NOTE: *NO* projectile destruction for mesons !!!
359 // target destruction
360 HDP.SetDefault( "FTF_MESON_NUCDESTR_P1_TGT", 0.00481, 0.0, 1.0 );
361 HDP.SetDefault( "FTF_MESON_NUCDESTR_P1_ADEP_TGT", true );
362 HDP.SetDefault( "FTF_MESON_NUCDESTR_P2_TGT", 4.0, 2.0, 16.0 );
363 HDP.SetDefault( "FTF_MESON_NUCDESTR_P3_TGT", 2.1, 0.0, 4.0 );
364 //
365 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P1", 0.035, 0.0, 0.25 );
366 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P2", 0.04, 0.0, 0.25 );
367 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P3", 4.0, 2.0, 16.0 );
368 HDP.SetDefault( "FTF_MESON_PT2_NUCDESTR_P4", 2.5, 0.0, 4.0 );
369 //
370 HDP.SetDefault( "FTF_MESON_NUCDESTR_R2", 1.5*CLHEP::fermi*CLHEP::fermi,
371 0.5*CLHEP::fermi*CLHEP::fermi,
372 2.0*CLHEP::fermi*CLHEP::fermi );
373 HDP.SetDefault( "FTF_MESON_EXCI_E_PER_WNDNUCLN", 40.0*CLHEP::MeV, 0.0, 100.0*CLHEP::MeV );
374 HDP.SetDefault( "FTF_MESON_NUCDESTR_DISP", 0.3, 0.1, 0.4 );
375 }
HDP
G4HadronicDeveloperParameters & HDP
Definition: G4FTFTunings.cc:195
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: