GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » R3BRoot » Simulation Issues » PDG ID of light ions
PDG ID of light ions [message #19419] Thu, 11 August 2016 14:43 Go to next message
Oleg is currently offline  Oleg
Messages: 19
Registered: August 2016
Location: GSI
occasional visitor
From: *gsi.de
Are there PDG IDs of the ions from t to let say Be?
If I see the original PDG tables, there are no,
but in R3BRoot one can set a very long number for t or alpha and they seems to be generated.

Is there a list of PDG IDs available in FAIRRoot/R3BRoot?

Oleg.
Re: PDG ID of light ions [message #19421 is a reply to message #19419] Fri, 12 August 2016 09:09 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
I am not aware of such list. FairRoot uses TDatabasePDG class to handle PDGs.

But PDG for ions will not help you. This long generated PDG is for internal Geant usage. Input ions have to be specified by setting Z and A in R3BIonGenerator or in AsciiGenerator, which will properly register a new ion in transport code.

Best regards,
Dima
Re: PDG ID of light ions [message #19428 is a reply to message #19421] Mon, 15 August 2016 14:00 Go to previous messageGo to next message
Oleg is currently offline  Oleg
Messages: 19
Registered: August 2016
Location: GSI
occasional visitor
From: *gsi.de
We have strange problem with Alphas (4He).
Should it work with AsciiGenerator - as the ion with A = 4 and Z = 2?
Up to now does not work.

Oleg.
Re: PDG ID of light ions [message #19429 is a reply to message #19428] Mon, 15 August 2016 14:14 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
Please post the line of your ASCII file where the alpha is defined. Yes, it should work with Z=2 and A=4. Do not forget -1 as first parameter. Example event:

0 1 0. 0.
  -1  2  4  px py pz vx vy vz 3.727379508


Or, alternatively, the Ion generator in r3ball.C:

R3BIonGenerator *gen = new R3BIonGenerator(2, 4, 2, 10, 0., 0., 1.); // 10 alpha's with pz = 1 GeV
primGen->AddGenerator(gen);
Re: PDG ID of light ions [message #19430 is a reply to message #19429] Mon, 15 August 2016 14:21 Go to previous messageGo to next message
Oleg is currently offline  Oleg
Messages: 19
Registered: August 2016
Location: GSI
occasional visitor
From: *gsi.de
0 1 0.0 0.0
-1 2 4 0.00 0.00 60.00 0.00 0.00 -40.00 3.7273

====================================================================

### Run 0 start.
-I- R3BAsciiGenerator: Reading Event: 0, pBeam = 0GeV, b = 0 fm, multiplicity 1
-W- R3BAsciiGenerator::ReadEvent: Cannot find Ion_4_2 in database!

TG4PrimaryGeneratorAction::TransformPrimaries:
No primary particles found on the stack.
*** TG4Exception: Aborting execution ***
Re: PDG ID of light ions [message #19438 is a reply to message #19430] Wed, 17 August 2016 13:10 Go to previous message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
The fix for not recognised Alpha particle is already committed to central repository (dev branch).

Best regards,
Dima
Previous Topic: Geo generation
Next Topic: ROOT geometry and Garfield geometry
Goto Forum:
  


Current Time: Thu Aug 11 02:10:17 CEST 2022

Total time taken to generate the page: 0.00859 seconds