GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Event Generators » benchmark generation
benchmark generation [message #15667] Tue, 05 November 2013 12:32 Go to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Hi,
I would like to generate a benchmark reaction and have parameterization from the experiment. Can I use EvtGen for that purpose?
Re: benchmark generation [message #15668 is a reply to message #15667] Tue, 05 November 2013 13:00 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
Yes, you need to write your own model, like this.
Re: benchmark generation [message #15681 is a reply to message #15668] Wed, 06 November 2013 17:14 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Thank you Stefano. Can I produce file for that like:
/input/psi2s_jpsi2pi_1k.evt (the same format) to use as benchmark?
Re: benchmark generation [message #15682 is a reply to message #15681] Wed, 06 November 2013 17:38 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Just to understand, exactly which channel you would like to simulate?
Re: benchmark generation [message #15684 is a reply to message #15682] Wed, 06 November 2013 20:39 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
First of all inclusive Lambda or antiLambda and then Lambda_c or antiLambda_c at our energy.
Re: benchmark generation [message #15685 is a reply to message #15684] Wed, 06 November 2013 21:17 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
You need to take the macros in macro/run, and write your own .dec file, no need to write your model since (at least for lambda) the model already exits.
In some presentations at the older collaboration meetings lambda analysis were shown with also the models which could be used (Lia or Simone talks).
Re: benchmark generation [message #15688 is a reply to message #15685] Thu, 07 November 2013 12:11 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Thank you again Stefano, may I ask you hoe I can simulate lambda-Lambdabar channel with the code EvtLambdaLambdaBar.cc (you've send me) in the frame of EvtGen, e.g. any example for benchmark as in sim_complete.C example there is only example with dec or evt files are presented? We want to do this kind of simulation (and then some other channels) for our forward TOF wall.
Re: benchmark generation [message #15689 is a reply to message #15688] Thu, 07 November 2013 12:24 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Your dec file should be like:

noPhotos
Decay pbarpSystem
1.0 anti-Lambda0 Lambda0 LambdaLambdaBar 1.643;
Enddecay

Decay Lambda0
1.0 p+ pi- PHSP;
Enddecay

Decay anti-Lambda0
1.0 anti-p-pi+ PHSP;
Enddecay

End


where 1.643 is the momentum of your antiproton (in this case it is 1.643, but you can oput your value).

Remember to use EvtGenDirect and set the same momentum value in the sim_complete.
Re: benchmark generation [message #15690 is a reply to message #15689] Thu, 07 November 2013 13:05 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Thank you! Is the value of anti-proton momentum 1.643 in cms system? Is it for the reaction p + pbar -> Lambda + Lambda-bar?

In sim_complete(Int_t nEvents = 10, TString SimEngine ="TGeant3", Float_t mom = 7.24) mom (momentum) is also defined in cms?

Re: benchmark generation [message #15691 is a reply to message #15690] Thu, 07 November 2013 13:06 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
No, it is in the lab frame.
Re: benchmark generation [message #15693 is a reply to message #15691] Fri, 08 November 2013 13:59 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Thank you Stefano, I noticed, that Lambda and lambda-bar are decaying in the interaction point, butn they should have a decay path, is it possible to change it? Is there a simple tool or an example (root macro) how to analyze file with simulated events without reconstruction and digitization (sim_complete.root file)?
Re: benchmark generation [message #15694 is a reply to message #15693] Fri, 08 November 2013 14:00 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
They do not decay at interaction point. Why do you say so?
Re: benchmark generation [message #15695 is a reply to message #15694] Fri, 08 November 2013 14:07 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
I just check the track number, for p-bar, pi+ ... it starts from 0,1 ..., but no Lambda MCtrack.
Re: benchmark generation [message #15696 is a reply to message #15694] Fri, 08 November 2013 15:15 Go to previous messageGo to next message
Yuri Naryshkin is currently offline  Yuri Naryshkin
Messages: 42
Registered: August 2008
continuous participant
From: *desy.de
Sorry Stefano, I've check again, the origin vertex for secondary part is not an interaction point.
Re: benchmark generation [message #15697 is a reply to message #15696] Fri, 08 November 2013 17:13 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
Ok, good.
Previous Topic: [CLOSED] DPM-Direct doesn't have Coulomb el.scattering part
Next Topic: EvtGen workshop
Goto Forum:
  


Current Time: Wed Nov 13 13:44:14 CET 2024

Total time taken to generate the page: 0.00831 seconds