GSI Forum: General » Geantino?

Home » PANDA » PandaRoot » General » Geantino?

Show: Today's Messages :: Polls :: Message Navigator

Geantino? [message #5831]

Fri, 08 February 2008 09:14

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi developers,

does anybody know how to create (standard / charged) geantinos with our box generator?

Thank's

Tobias

Report message to a moderator

Re: Geantino? [message #5832 is a reply to message #5831]

Fri, 08 February 2008 09:22

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *dip.t-dialin.net

Ahem,
what is exactly a "geantino"?

Report message to a moderator

Re: Geantino? [message #5833 is a reply to message #5832]

Fri, 08 February 2008 09:44

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi Stefano,

a Geantino is an artifical particle from Geant which is used as a geometrical probe. It is only transported through your detector and does not interact with it.

Ciao,

Tobias

Report message to a moderator

Re: Geantino? [message #5844 is a reply to message #5831]

Mon, 11 February 2008 11:09

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi developers,

I have solved part of the geantino problem. The geantino has (in VirtualMC) the pdg code 0. I still do no know how to handle charge geantinos.

Ciao,

Tobias

Report message to a moderator

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);

Report message to a moderator

Re: Geantino? [message #5847 is a reply to message #5846]

Mon, 11 February 2008 11:31

asanchez
Messages: 350
Registered: March 2006

first-grade participant

From: *kph.uni-mainz.de

Sorry,
CbmParticle *chargedGeant = new CbmParticle("chargedGeant", Z,A,1, mass,Z,kFALSE,lifetime);

Report message to a moderator

Re: Geantino? [message #5848 is a reply to message #5846]

Mon, 11 February 2008 11:39

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi Alicia,

thank you for your reply.

I have tested your suggestion. Unfortunately it fails with the following error message:

TG4PrimaryGeneratorAction::TransformPrimaries:
G4ParticleTable::FindParticle() failed for crootino pdgEncoding=10000000.
*** TG4Exception: Aborting execution ***

Do you have any idea how to solve this problem?

Tanks,

Tobias

Report message to a moderator

Re: Geantino? [message #5849 is a reply to message #5845]

Mon, 11 February 2008 11:41

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi Ralf,

is there a charged rootino, too?

Thanks,

Tobias

Report message to a moderator

Re: Geantino? [message #5850 is a reply to message #5849]

Mon, 11 February 2008 11:46

asanchez
Messages: 350
Registered: March 2006

first-grade participant

From: *kph.uni-mainz.de

Hi Tobias, try to do that

CbmParticle *chargedGeant = new CbmParticle("chargedGeant",
Z,A, mass,Z,kFALSE,lifetime);

sorry there was one additional field inside, remove the one
in the constructor.

Report message to a moderator

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);

Report message to a moderator

Re: Geantino? [message #5854 is a reply to message #5831]

Mon, 11 February 2008 13:03

Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt

first-grade participant

From: 141.30.85*

Hi again,

I found this description in the root classes:

Quote:

TDatabasePDG
Particle database manager class

This manager creates a list of particles which by default is
initialised from with the constants used by PYTHIA6 (plus some
other particles added). See definition and the format of the default
particle list in $ROOTSYS/etc/pdg_table.txt

there are 2 ways of redefining the name of the file containing the
particle properties

1. one can define the name in .rootrc file:

Root.DatabasePDG: $(HOME)/my_pdg_table.txt

2. one can use TDatabasePDG::ReadPDGTable method explicitly:

- TDatabasePDG *pdg = new TDatabasePDG();
- pdg->ReadPDGtable(filename)

See TParticlePDG for the description of a static particle properties.
See TParticle for the description of a dynamic particle particle.

This means you could just create your own charged rootino. Additionally I found the rootino at the pdg ID 521 in root/etc/pdg_table.txt.

Ralf.

Report message to a moderator

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

Report message to a moderator

Previous Topic:	svn and our policies
Next Topic:	PdgCode > 1000000000 ?

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Nov 29 12:43:42 CET 2024

Total time taken to generate the page: 0.00765 seconds