About pion and muon productions

From: Hantao Jing <jinght_at_mail.ihep.ac.cn>
Date: Mon, 04 May 2009 22:32:08 +0800

Dear Fluka experts,
      I want to simulate the muon production by the reaction of
proton-Carbon at 1.6GeV. At the same time, I also want to know about
pion+- productions in this reaction. For different-thickness Carbon
targets, I calculated the pion+- production as showed in figure(in
attachments).
 From the figure, we find that there is a very mild depth at pion energy
of about 100MeV for pion production spectra. For thin target, the
contribution from first peak dominates. With increasing the thickness of
Carbon target, namely, 5cm, 10cm and 20cm, the second peak value become
more important and even dominate at last. But I find there isn't this
depth structure for the other calculated results from 'Treating High
Energy Showers'(Alfredo Ferrari and Paola R. Sala) and experimental
data(PHYSICAL REVIEW C 75, 024611 (2007)). Why does the depth structure
appear? Is some wrong with my calculation by FLUKA?

    Secondly, the pion productions mainly are through the following the
three mechanisms, namely,
1. single pion production
$NN \rightarrow \pi NN$;

2. two-body process
$NN \rightarrow \pi d$;

3. pion-pair production
$NN \rightarrow NN \pi \pi$.

In fluka, the above three mechanisms are all included ?

For my input file(see attachments), are there any suggestions to me in
order to obtain the better simulation results for pion and muon production?

*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
muon fluence inside and around a proton-irradiated c target
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
BEAM -1.6 PROTON
BEAMPOS 0.0 0.0 -50.0
GEOBEGIN COMBNAME
    0 0 A cylinder C target inside vacuum
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* BX1
  RPP BX1 -50000.00 50000.00 -50000.00 50000.00 -50000.00 50000.00
* BX2
  RPP BX2 -10000.00 10000.00 -10000.00 10000.00 -10000.00 10000.00
* CY1
  RCC CY1 0.00 0.00 0.00 0.00 0.00 100.00
               30.00
* CY2
  RCC CY2 0.00 0.00 0.00 0.00 0.00 100.00
               60.00
* PL1
  XYP PL1 10.00
  END
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Reg # 1
* BH1; assigned material: Blackhole; mat # (1)
BH1 5 +BX1 -BX2
* Reg # 2
* VA1; assigned material: Vacuum; mat # (2)
VA1 5 +BX2 -CY2
* Reg # 3
* CR1; assigned material: Carbon; mat # (6)
CR1 5 +CY1 +PL1
* Reg # 4
* VA2; assigned material: Vacuum; mat # (2)
VA2 5 +CY2 -PL1
* Reg # 5
* CR2; assigned material: Vacuum; mat # (2)
* CR2 5 +CY2 +PL1
* Reg # 5
* VA3; assigned material: Vacuum; mat # (2)
VA3 5 +CY2 -CY1 +PL1
  END
GEOEND
ASSIGNMAT 1.0 1
ASSIGNMAT 2.0 2
ASSIGNMAT 6.0 3
ASSIGNMAT 2.0 4
ASSIGNMAT 2.0 5

PART-THR -1.0E-6 10.0 14.0 1.0 0.0
PART-THR -1.0E-6 1.0 2.0 1.0 0.0
PART-THR -1.0E-6 23.0 23.0 0.0 0.0

* Boundary crossing uz current in the edge of the target (muon log, one-way)
USRBDX -1.0 212.0 -41.0 3.0 4.0 2827.43uzCurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing uz+ current in the edge of the target (log, one-way)
USRBDX -1.0 10.0 -41.0 3.0 4.0 2827.43uz+FluenUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing uz-current in the edge of the target (log, one-way)
USRBDX -1.0 11.0 -41.0 3.0 4.0 2827.43uz-CurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*
* Boundary crossing piz current in the edge of the target (pion log, one-way)
USRBDX -1.0 209.0 -41.0 3.0 4.0 2827.43pzCurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing piz+ current in the edge of the target (log, one-way)
USRBDX -1.0 13.0 -41.0 3.0 4.0 2827.43pz+FluenUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing piz- current in the edge of the target (log, one-way)
USRBDX -1.0 14.0 -41.0 3.0 4.0 2827.43pz-CurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing piz0 current in the edge of the target (log, one-way)
USRBDX -1.0 23.0 -41.0 3.0 4.0 2827.43pz0CurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
* Boundary crossing proton&neutrons current in the edge of the target (log, one-way)
USRBDX -1.0 207.0 -41.0 3.0 4.0 2827.43pnzCurrUD
USRBDX +1.6 1.0E-8 +100.0 0.0 1.0 &
*
*
RANDOMIZE 1.0
START 50000000.0
STOP

Received on Tue May 05 2009 - 00:12:11 CEST

This archive was generated by hypermail 2.2.0 : Tue May 05 2009 - 00:12:14 CEST