Message: Re: neutron beam hit Zinc target Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: neutron beam hit Zinc target 

Forum: Hadronic Processes
Re: Question neutron beam hit Zinc target (Tuoyu Cao)
Re: None Re: neutron beam hit Zinc target (Michael H. Kelsey)
Re: Question Re: neutron beam hit Zinc target (Tuoyu Cao)
Date: 01 Mar, 2013
From: Michael H. Kelsey <Michael H. Kelsey>

Tuoyu Cao writes:
> Here is the code I use to define Zinc.
> 
>   G4Isotope* iso_Zn64 = new G4Isotope("Zn64", iz=30, ia=64, density= 64*g/mole);
>   G4Isotope* iso_Zn66 = new G4Isotope("Zn66", iz=30, ia=66, density= 66*g/mole);
>   G4Isotope* iso_Zn67 = new G4Isotope("Zn67", iz=30, ia=67, density= 67*g/mole);
>   G4Isotope* iso_Zn68 = new G4Isotope("Zn68", iz=30, ia=68, density= 68*g/mole);
>   G4Isotope* iso_Zn70 = new G4Isotope("Zn70", iz=30, ia=70, density= 70*g/mole);
> 
>   G4Element* ele_Zn = new G4Element("ele Zn", symbol="eZn",ncomponents=5);
> 
>   ele_Zn->AddIsotope(iso_Zn64,abundance=48.60*perCent);
>   ele_Zn->AddIsotope(iso_Zn66,abundance=27.90*perCent);
>   ele_Zn->AddIsotope(iso_Zn67,abundance=04.10*perCent);
>   ele_Zn->AddIsotope(iso_Zn68,abundance=18.80*perCent);
>   ele_Zn->AddIsotope(iso_Zn70,abundance=00.60*perCent);
> 
>   G4Material* mat_ele_Zn = new G4Material("mix_Zn", density= 7.14*g/cm3, ncomponents=1);
>   mat_ele_Zn->AddElement(ele_Zn, fractionmass =1);
> 
> But G4 uses nat_zinc instead. What should I do if I want to simulate one
> isotope individually?

If you want to simulate a pure isotope, just create the G4Element object
with one single entry, the isotop you want with an abundance of 1.00.

> *********************************************************************************************************
> * G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
> *********************************************************************************************************
> Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
>     0     -5 mm   2.72 mm  -4.68 mm    2.5 MeV     0 eV      0 fm      0 fm   Absorber   initStep
>     1   3.57 mm   2.72 mm  -4.68 mm   2.48 MeV     0 eV   8.57 mm   8.57 mm   Absorber   hadElastic
>     :----- List of secondaries ----------------
>        Zn65[0.0]:  energy =  16.3 keV  time = 392.5 ps 
>     :------------------------------------------
>     2   4.01 mm   2.58 mm     -5 mm   2.48 MeV     0 eV    565 um   9.13 mm   Absorber   Transportation
>     3   5.42 mm   2.14 mm     -6 mm   2.48 MeV     0 eV   1.78 mm   1.09 cm      World   OutOfWorld
> 
> shows there is only one secondary, which is Zn65, and it seems like it
> comes from hadElastic.... Which is very confusing. Is this an error in
> G4. Could you please explain this.

I don't think it's an error, but I agree that it is quite confusing.  One of
the difficulties is that the tracking table doesn't tell you exactly what
nucleus from the material was chosen for that interaction.  

The fact that the neutron continues on its way (steps 2 and 3) means that
there really was an elastic or quasielastic interaction.  But in that case,
the Zn65 must have been a recoil, not production.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Question: Re: neutron beam hit Zinc target   (Tuoyu Cao - 01 Mar, 2013)
1 None: Re: neutron beam hit Zinc target   (Mike Kelsey - 01 Mar, 2013)
2 None: Re: neutron beam hit Zinc target   (Koi, Tatsumi - 01 Mar, 2013)
 Add Message Add Message
to: "Re: neutron beam hit Zinc target"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]