GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Event Generators » Change particles Stable/Unstable in DPM
Change particles Stable/Unstable in DPM [message #16268] Tue, 08 April 2014 17:48 Go to next message
Aida Galoyan is currently offline  Aida Galoyan
Messages: 79
Registered: May 2007
Location: Dubna
continuous participant
From: *jinr.ru
Hi all,

According to request of Klaus and Stefano, I have implemented the possibility to change
a status of particles: stable or unstable in DPM.

Now, when you run DPMGen, after lines
-------------------------------------------
Enter Elastic : 0., 1. or 2.
0. - No elastic scattering, only inelastic
1. - Elastic and inelastic interactions
2. - Only elastic scattering, no inelastic one
0.
-----------------------------------------
there will be lines
-----------------------------------------
If you like to change particle status (stable - 1 / unstable - 0)
Enter: particle PDGcode and status
(for example -- Pi0 unstable: 111 0)
To go to event generation, Enter: 0 0
---------------------------------------------
If you Enter
--------------------------------------------
111 0
0 0
--------------------------------------------

Pi0 will be unstable, other particles save their statuses.

Aida
Re: Change particles Stable/Unstable in DPM [message #16283 is a reply to message #16268] Thu, 10 April 2014 08:10 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Hi,


this functionallity now also is available for PndDpmDirect via methods SetStable(int pdg) and SetUnstable(int pdg) and can be used e.g. like:

PndDpmDirect *Dpm= new PndDpmDirect(Mom,1);
Dpm->SetUnstable(111);   // pi0
Dpm->SetUnstable(310);   // K_S0
Dpm->SetStable(3122);    // Lambda
Dpm->SetStable(-3122);   // anti-Lambda
primGen->AddGenerator(Dpm);


Note that the defaults for pi0, K_S0 and (anti-)Lambda is 'stable'.


Best,
Klaus
Re: Change particles Stable/Unstable in DPM [message #16837 is a reply to message #16268] Wed, 25 June 2014 16:39 Go to previous messageGo to next message
Aida Galoyan is currently offline  Aida Galoyan
Messages: 79
Registered: May 2007
Location: Dubna
continuous participant
From: *jinr.ru
Hi all,

When I entered a procedure of changing a particle status -Stable/Unstable, I made a small bug in DpmEvtGen/decay.f .
Now, bug is fixed and commited in trunk.

Aida
Re: Change particles Stable/Unstable in DPM [message #16838 is a reply to message #16837] Wed, 25 June 2014 22:25 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 93.48.228*
Dear Aida,
just to understand, which kind of bug was it? A minor bug or something more critical? I am asking this since I have produced recently many DPM events for my calculations, and I'd like to understand if I should redo my simulation or not.
Re: Change particles Stable/Unstable in DPM [message #16842 is a reply to message #16268] Thu, 26 June 2014 11:21 Go to previous messageGo to next message
Aida Galoyan is currently offline  Aida Galoyan
Messages: 79
Registered: May 2007
Location: Dubna
continuous participant
From: *jinr.ru
Dear Stefano,
I try to explain the bug.

When I introduced in DPM changing of particle status(stable / unstable),
transformations of K0(PDG code 311) and antiK0(PDG code -311) to K0S and K0L
were not implemented. Thus, if you do not consider K0S and K0L in your study on the
generating level, all is normal. However, if you consider subsequent simulation by
Geant3/Geant4, then (I think) Geant will transform K0S and K0L to K0S and K0L.

Now, these transformations are introduced in DPM.

Aida
Re: Change particles Stable/Unstable in DPM [message #16843 is a reply to message #16842] Thu, 26 June 2014 11:45 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Thanks for your explanations.
Previous Topic: Event filter for (a combination of) arbitrary event generators
Next Topic: New generator ChiGen
Goto Forum:
  


Current Time: Wed Dec 07 07:48:39 CET 2022

Total time taken to generate the page: 0.01036 seconds