GSI Forum: Event Generators » Change particles Stable/Unstable in DPM

Home » PANDA » PandaRoot » Event Generators » Change particles Stable/Unstable in DPM

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

Change particles Stable/Unstable in DPM [message #16268]

Tue, 08 April 2014 17:48

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

Report message to a moderator

Re: Change particles Stable/Unstable in DPM [message #16283 is a reply to message #16268]

Thu, 10 April 2014 08:10

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

Report message to a moderator

Re: Change particles Stable/Unstable in DPM [message #16837 is a reply to message #16268]

Wed, 25 June 2014 16:39

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

Report message to a moderator

Re: Change particles Stable/Unstable in DPM [message #16838 is a reply to message #16837]

Wed, 25 June 2014 22:25

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.

Report message to a moderator

Re: Change particles Stable/Unstable in DPM [message #16842 is a reply to message #16268]

Thu, 26 June 2014 11:21

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

Report message to a moderator

Re: Change particles Stable/Unstable in DPM [message #16843 is a reply to message #16842]

Thu, 26 June 2014 11:45

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Thanks for your explanations.

Report message to a moderator

Previous Topic:	Event filter for (a combination of) arbitrary event generators
Next Topic:	New generator ChiGen

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Mon Nov 25 23:16:16 CET 2024

Total time taken to generate the page: 0.00556 seconds