GSI Forum: Event Generators » [FIXED] Problem with UseEvtGenDirect

Home » PANDA » PandaRoot » Event Generators » [FIXED] Problem with UseEvtGenDirect

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

[FIXED] Problem with UseEvtGenDirect [message #14294]

Thu, 03 January 2013 16:34

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

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 2.235.190*

Which momentum value are you using?

Report message to a moderator

Re: Problem with UseEvtGenDirect [message #14296 is a reply to message #14295]

Thu, 03 January 2013 16:49

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.

Report message to a moderator

Re: Problem with UseEvtGenDirect [message #14297 is a reply to message #14296]

Thu, 03 January 2013 17:13

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?

Report message to a moderator

Re: Problem with UseEvtGenDirect [message #14298 is a reply to message #14297]

Thu, 03 January 2013 17:22

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

Report message to a moderator

Re: Problem with UseEvtGenDirect [message #14299 is a reply to message #14298]

Thu, 03 January 2013 17:24

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.

Report message to a moderator

Previous Topic:	New pbarpSystem and pbardSystem implementations
Next Topic:	Update to new EvtGen

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Dec 01 01:29:49 CET 2024

Total time taken to generate the page: 0.00852 seconds