Message: What is the correct way to create a G4Element consisting of a single G4Isotope? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question What is the correct way to create a G4Element consisting of a single G4Isotope? 

Forum: Geometry
Date: 24 Oct, 2018
From: Sergio Losilla <Sergio Losilla>

Hi,

I am wondering what is the correct way to create an element consisting of a single isotope. Everything looks fine except for the mass number N:

    auto b10{new G4Element{"B10", "B", /* number of components : */ 1}};
    b10->AddIsotope(new G4Isotope{"B10", 5, 10}, /* relative abundance : */ 1.0);

    std::cout << "Z: " << b10->GetZ() << std::endl;
    std::cout << "N: " << b10->GetN() << std::endl;
    std::cout << "A: " << b10->GetA()/(g/mole) << std::endl;

Output:

    Z: 5
    N: 10.0129
    A: 10.0129

I noticed that in G4Element.cc (lines 97 and 196) the effective mass number is just the effective atomic mass in different units:

    fNeff   = fAeff/(g/mole);

Should I expect that Neff for the element differs from N for the isotope? If not, how should I create the G4Element? If yes, what does Neff mean?

Best regards,

Sergio Losilla

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

1 Feedback: Re: What is the correct way to create a G4Element consisting of a single G4Isotope?   (Michael H. Kelsey - 28 Oct, 2018)
1 None: Re: What is the correct way to create a G4Element consisting of a single G4Isotope?   (Sergio Losilla - 29 Oct, 2018)
2 None: Re: What is the correct way to create a G4Element consisting of a single G4Isotope?   (michel maire - 05 Nov, 2018)
1 None: Re: What is the correct way to create a G4Element consisting of a single G4Isotope?   (Sergio Losilla - 06 Nov, 2018)
 Add Message Add Message
to: "What is the correct way to create a G4Element consisting of a single G4Isotope?"

 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 ]