GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Event Generators » [FIXED] Problem with UseEvtGenDirect
[FIXED] Problem with UseEvtGenDirect [message #14294] Thu, 03 January 2013 16:34 Go to next message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
Hi there and best wishes for 2013!

I try to run macro/run/sim_complete.C using the UseEvtGenDirect.

My .dec file looks like this:
noPhotos

Decay pbarpSystem
  1.0 J/psi pi+ pi- PHSP;
Enddecay

End


Running this results in the error message:
root.exe: /vol0/panda/pandasep12/pandaroot/pgenerators/EvtGen/EvtGenBase/EvtIntervalFlatPdf.cc:19: 
EvtIntervalFlatPdf::EvtIntervalFlatPdf(double, double): Assertion `max >= min' failed.


If I replace the "J/psi" by a "phi" everything works.

What is the problem here?

[Updated on: Wed, 10 September 2014 14:52] by Moderator

Report message to a moderator

Re: Problem with UseEvtGenDirect [message #14295 is a reply to message #14294] Thu, 03 January 2013 16:38 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*
Which momentum value are you using?
Re: Problem with UseEvtGenDirect [message #14296 is a reply to message #14295] Thu, 03 January 2013 16:49 Go to previous messageGo to next message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
I used BeamMomentum =10.0;

I tried 15.0 GeV as well. Then if I add a decay like
Decay J/psi
  1.0 e+ e- VLL;
Enddecay
I get the same error message.
Without decay chain the error message is different:
EvtGen:End of decay chain.
EvtGen:Sorry, no luck finding a valid set of masses.  This may be a pathological combo
root.exe: /vol0/panda/pandasep12/pandaroot/pgenerators/EvtGen/EvtGenBase/EvtParticle.cc:466: void EvtParticle::generateMassTree(): Assertion `0' failed.





Re: Problem with UseEvtGenDirect [message #14297 is a reply to message #14296] Thu, 03 January 2013 17:13 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*
I have just tried with the latest trunk w/o any problem.

Could you please check in the log that you have the following text:

<I> PndEvtGenDirect
<I> Particle: pbarpSystem
<I> decfile: /home/spataro/sep12/pandaroot/macro/run/ronald.dec
<I> pbar-Momentum: 10
<I> Rnd Seed: -1
<I> Rnd Seed changed to 3319
EvtGen:Initializing EvtGen
EvtGen:Storing known decay models
Will register EvtPythia
EvtGen:Main decay file name  :/home/spataro/sep12/pandaroot/pgenerators/EvtGen/DECAY.DEC
EvtGen:PDT table file name   :/home/spataro/sep12/pandaroot/pgenerators/EvtGen/evt.pdl
EvtGen:Initializing RadCorr=PHOTOS
EvtGen:No RadCorr engine given in EvtGen::EvtGen constructor, will use default EvtPHOTOS.
EvtGen:In readDecayFile, reading:/home/spataro/sep12/pandaroot/pgenerators/EvtGen/DECAY.DEC
EvtGen:As requested, PHOTOS will be turned on.


and

EvtGen:In readDecayFile, reading:/home/spataro/sep12/pandaroot/macro/run/ronald.dec
EvtGen:As requested, PHOTOS will be turned off.
EvtGen:Given allowed decays, resetting minMass pbarpSystem 1.48 to 3.37581
EvtGen:Reading /home/spataro/sep12/pandaroot/macro/run/ronald.dec to override decay table.

############# Generating with following conditions:

incident 4-mom : (10.9822, 0, 0, 10), m = 4.53966

######################


and copy it here?
Re: Problem with UseEvtGenDirect [message #14298 is a reply to message #14297] Thu, 03 January 2013 17:22 Go to previous messageGo to next message
Ronald Kunne is currently offline  Ronald Kunne
Messages: 32
Registered: October 2009
continuous participant
From: *in2p3.fr
Good call, Stefano!

When I checked this, I saw that it uses "mom" and not "Beam momentum" to determine the incident beam.

So it's solved. Thanks very much.

Greetings,
Ronald
Re: Problem with UseEvtGenDirect [message #14299 is a reply to message #14298] Thu, 03 January 2013 17:24 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
"mom" is the momentum for the generator, i.e. pbar momentum for evtgen, particle gun momentum for box, pbar momemtum for dpm.
BeamMomentum is the momentum for the scaling of the dipole field.
Previous Topic: New pbarpSystem and pbardSystem implementations
Next Topic: Update to new EvtGen
Goto Forum:
  


Current Time: Thu Dec 05 23:47:11 CET 2024

Total time taken to generate the page: 0.00652 seconds