Home » PANDA » PandaRoot » General » Geantino?
|
|
|
|
Re: Geantino? [message #5845 is a reply to message #5831] |
Mon, 11 February 2008 11:11 |
Ralf Kliemt
Messages: 507 Registered: May 2007 Location: GSI, Darmstadt
|
first-grade participant |
From: 141.30.85*
|
|
Hi Tobias,
The geantino (actually we use the rootino) can be accessed via the PDG number 0. in the simualtion.
//ROOTino
CbmBoxGenerator *fBox1 = new CbmBoxGenerator(0, 1);
To give you the right feeling try to check inside root:
root [0] TDatabasePDG* pdgBase = TDatabasePDG::Instance()
root [1] TParticlePDG *particle = pdgBase->GetParticle(0)
root [2] particle->GetName()
(const char* 0x8466b18)"Rootino"
root [3] particle->Charm()
(const Int_t)(0)
root [4] particle->Charge()
(const Double_t)0.00000000000000000e+00
root [5] particle->Lifetime()
(const Double_t)0.00000000000000000e+00
root [6] particle->Mass()
(const Double_t)0.00000000000000000e+00
root [7] particle->Width()
(const Double_t)0.00000000000000000e+00
root [8] particle->Stable()
(const Int_t)(1)
root [9] particle->ParticleClass()
(const char* 0x84652c8)"Unknown"
root [10] particle->Spin()
(const Double_t)0.00000000000000000e+00
root [11] particle->Parity()
(const Int_t)(0)
...
Greetings from Dresden,
Ralf.
[Updated on: Mon, 11 February 2008 11:16] Report message to a moderator
|
|
|
Re: Geantino? [message #5846 is a reply to message #5845] |
Mon, 11 February 2008 11:19 |
asanchez
Messages: 350 Registered: March 2006
|
first-grade participant |
From: *kph.uni-mainz.de
|
|
hallo Tobias,
you can add your charge geantino doing as follows
best regrads
ALicia
CbmParticle *chargedGeant = new CbmParticle("He5", Z,A,1, mass,Z,kFALSE,lifetime);
fRun->AddNewParticle(cgargedGeant );
Int_t pdg;
TDatabasePDG *db= TDatabasePDG::Instance();
TParticlePDG *p=0;
p=db->GetParticle("chargedGeant");
if(p) pdg = p->PdgCode();
CbmParticleGenerator* partGen = new CbmParticleGenerator(pdg, 1, -0.3,0.3,0.
3, 0., 0., 0);
primGen->AddGenerator(partGen);
|
|
|
|
|
|
|
Re: Geantino? [message #5851 is a reply to message #5850] |
Mon, 11 February 2008 12:07 |
asanchez
Messages: 350 Registered: March 2006
|
first-grade participant |
From: *kph.uni-mainz.de
|
|
i have found the following information in geant4
for the charged geantino
Arguments for constructor are as follows
00062 // name mass width charge
00063 // 2*spin parity C-conjugation
00064 // 2*Isospin 2*Isospin3 G-parity
00065 // type lepton number baryon number PDG encoding
00066 // stable lifetime decay table
00067 // shortlived subType anti_encoding
00068 anInstance = new G4ParticleDefinition(
00069 name, 0.0*MeV, 0.0*MeV, +1.*eplus,
00070 0, 0, 0,
00071 0, 0, 0,
00072 "geantino", 0, 0, 0,
00073 true, 0.0, NULL,
00074 false, "geantino", 0
00075 );
maybe it works for the following options
CbmParticle *chargedGeant = new CbmParticle("chargedGeant",
1,0, 0,1,kTRUE,0);
|
|
|
|
Re: Geantino? [message #5855 is a reply to message #5831] |
Mon, 11 February 2008 13:47 |
Tobias Stockmanns
Messages: 489 Registered: May 2007
|
first-grade participant |
From: *ikp.kfa-juelich.de
|
|
Dear Ralf and Alicia,
with both of you ways (my_pdg_table.txt / fRun->AddNewParticle) I am able to create a charged rootino in root. But it seems that this particle is not correctly transformed into Geant because I always get the following error message:
TG4PrimaryGeneratorAction::TransformPrimaries:
G4ParticleTable::FindParticle() failed for ChargedGeantino pdgEncoding=1001000.
*** TG4Exception: Aborting execution ***
Tobias
|
|
|
Goto Forum:
Current Time: Fri Nov 08 06:25:00 CET 2024
Total time taken to generate the page: 0.00829 seconds
|