GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » General » pbarpSystem
pbarpSystem [message #15094] Mon, 05 August 2013 12:47 Go to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
Hi all,

I was wondering why the mass of the pbarpSystem in evt.pdl is set to 2.98 GeV. How is this calculated and is this used in the simulation macro?
And why does it has a width of 0.1 GeV?

Additionally, please change the description in evt.pdl for the spin. It's J (total angular momentum quantum number) not S (spin)! It is just confusing..

Best regards,
Simon
Re: pbarpSystem [message #15096 is a reply to message #15094] Mon, 05 August 2013 14:14 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 Simon,


the values for the pbarpSystem are arbitrary and are setup with reasonable values when EvtGen starts - so the default values don't mean anything as far as I know.


Best,
Klaus
Re: pbarpSystem [message #15100 is a reply to message #15096] Mon, 05 August 2013 15:15 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*
Exactly.
And the "spin" is the EvtGen style, I will not touch it. In theory the user should not play with evt.pdl.
Re: pbarpSystem [message #15102 is a reply to message #15094] Mon, 05 August 2013 15:33 Go to previous messageGo to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
I just used it to modify evt.pdl for a dummy particle for my simulation to set mass, width and spin and it took me some time till I recognized, that is not the spin S.
Re: pbarpSystem [message #15106 is a reply to message #15100] Mon, 05 August 2013 16:48 Go to previous messageGo to next message
Jens Sören Lange is currently offline  Jens Sören Lange
Messages: 193
Registered: June 2005
first-grade participant
From: *wlan.uni-giessen.de
Hi Stefano, I disagree.

I think the user should actually not be discouraged to change evt.pdl.

1.) Our evt.pdl is hopelessly out of date, many masses and widths are wrong (it means different from PDG).

2.) The default width of the pbarsystem 0.1 GeV is definitely wrong (should be zero) and if we are not mistaken, it is probably not initialized by EvtGen with a "reasonable value". At least I see different results if you keep the default or you put in zero width "by hand" in the evt.pdl.

The default pbarpsystem mass of 2.98 GeV is maybe the eta_c mass. As Klaus says, it is overwritten when EvtGen starts, but you see what happened here: somebody put in whatever value into evt.pdl and then checked it into svn.

3.) modifying evt.pdl is still the best way to introduce a new particle.

So, I think, if you want to show any PandaRoot result on a conference (we want to show on the CHARM'13 in 3 weeks from now), we have to check and maybe modify the evt.pdl.

cheers, Soeren
Re: pbarpSystem [message #15107 is a reply to message #15106] Mon, 05 August 2013 19:09 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*
Hi Soeren,
I was not clear, ok, I admit it.

I think it makes no sense to change from our side "spin" into "J", this is just semantics. The evt.pdl style, definition and header" come from EvtGen, maybe one could raise the question to the developers (I mean the Warwick team).

With the new EvtGen we have an updated evt.pdl (check EvtGenNew/EvtGen/Private). If we start to modify the definitions of the particles "randomly", this could be risky. Of course one can use the evt.pdl that he wants.
The basic idea should be that particle properties are modified inside the dec file, and not in the common evt.pdl, which is the standard for all the users and analyses.

About pbarpSystem definition, this is what was used for the physics book. I believe it was hardly checked, but it would be good if somebody could check again if the results are changed by modifying its properties.
Re: pbarpSystem [message #15109 is a reply to message #15094] Tue, 06 August 2013 11:19 Go to previous messageGo to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
How to define a new particle in the dec file? I though the dummies in evt.pdl should be used for that. And is this already possible in the apr13 rel?

[Updated on: Tue, 06 August 2013 11:20]

Report message to a moderator

Re: pbarpSystem [message #15114 is a reply to message #15109] Wed, 07 August 2013 00:21 Go to previous messageGo to next message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *web.vodafone.de
Hi Simon,

to simulate a new particle with EvtGen you can use the particles listed as 'dummiesXX' in the file evt.pdl. For instance, let's assume that you want to simulate a particle called Y, with mass = A, width = B, spin and charge set up to 0 (remember to pass the values of mass a width in [GeV]). Then you need to do the following in your file.dec:

###########################
Particle dummy00_1 A B
ChangeMassMin dummy00_1 [min-mass-value]
ChangeMassMax dummy00_1 [max-mass-value]
#
Alias MyY dummy00_1
Decay MyY
1.0 [daughter-one] [daughter-two] [daughte-whatever] [model]
Enddecay

Indeed, if you need to simulate a new particle of spin 1 and charge =0, then you need to use 'dummy01_1' in your file.dec; you can find the list of all dummy particles available in the file evt.pdl. By default, the dummy00_1 gives you a LundID = 51.

I hope it helps, Elisabetta
Re: pbarpSystem [message #15117 is a reply to message #15094] Wed, 07 August 2013 11:24 Go to previous messageGo to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
Thank you very much for this introduction.
But why do i have to use these?
ChangeMassMin dummy00_1 [min-mass-value]
ChangeMassMax dummy00_1 [min-mass-value]
I'd set the max deviation from mass in evt.pdl to 0. Is this possible, too?
Re: pbarpSystem [message #15121 is a reply to message #15117] Wed, 07 August 2013 11:51 Go to previous messageGo to next message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ext.kfa-juelich.de
Hi Simon,

it is your choice. You can also set up 0 as DeltaMass. I always prefer to set up a range. But basically it should work. You can try.

Elisabetta
Re: pbarpSystem [message #15123 is a reply to message #15094] Wed, 07 August 2013 14:04 Go to previous messageGo to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
When setting max deviation from mass in evt.pdl to 0, no mass cut is used. This does not seems possible here, except for setting ChangeMassMax to 0 and ChangeMassMin to something high like 10GeV.
Re: pbarpSystem [message #15124 is a reply to message #15123] Wed, 07 August 2013 14:36 Go to previous messageGo to next message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ext.kfa-juelich.de
And what happens, then, if you comment the lines related to the Delta Mass cut in you file.dec? I see that in PandaRoot, in the file /EvtGen/evt.pdl, there is the possibility to set up a 'max_Dm/GeV'. The command which I suggested you, it allows to set up a min and max mass value. In my case it is taken into account: I just tried and on the prompt I see that the Delta Mass is reset according the the mass values which I just passed by my file.dec, while the sim-macro is running.
If you do not pass any value there, the default value 'max_Dm/GeV' of the dummyXX from the evt.pdl file is supposed to be used.

cheers, Elisabetta

Re: pbarpSystem [message #15126 is a reply to message #15094] Wed, 07 August 2013 14:50 Go to previous messageGo to next message
Simon Reiter is currently offline  Simon Reiter
Messages: 44
Registered: May 2013
continuous participant
From: *zeu.uni-giessen.de
If I delete the two lines i get the following error:
root.exe: /home/panda/pandaroot/apr13/pgenerators/EvtGen/EvtGenBase/EvtIntervalFlatPdf.cc:19: EvtIntervalFlatPdf::EvtInterv$

And if set to both to 0, i get the same error. If commented, it seems to work without errors.
Re: pbarpSystem [message #15127 is a reply to message #15126] Wed, 07 August 2013 16:48 Go to previous message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ikp.kfa-juelich.de
Hi Simon,
to set up both value to 0 is not correct, as EvtGen is written in a way that a mass interval should be passed somewhere. Now, you can choose those limits your own, or you can accept the default value of Delta Mass of the dummyXX in evt.pdl; but it must be written somewhere. So, if you comment those 2 lines you accept the default setting in evt.pdl, and this should be fine.
When I generated my new particle, I set up the value depending on the limit of the PHSP distribution of the daughters of my new particle, and it worked out. I never tried without any Delta mass value.

Just one comment: you can decide in this way to set up the particle properties (not only the ones related to the dummy particles, but also for instance for the D0, Ks, Lambda..) from your file.dec, without modifying evt.pdl. You can redefine for the particle of the evt.pdl list the mass, width, Lund-ID,....this is still fine.

Dummy00_1, for instance, by default is just a particle with mass = 1 GeV/c^2, width = 0.1, LundID= 51, spin =0 and charge =0. No additional properties. From your file.dec you are able to redefine the mass and the width, as from the example I had sent you, and you accept by default the other properties of this dummy00_1 particle (spin, charge). Make sure in the simulation macro to set up properly the momentum value, depending on the mass of your new particle, and the name of the Resonant state, wherever it is called in the macro. You are now using a particle calle "dummyXX_Y", whatever it is in your case.

Elisabetta
Previous Topic: How to apply transformations on the geometry? (Alignment)
Next Topic: Why is there a hardcoded $PATH in config.sh?
Goto Forum:
  


Current Time: Tue Mar 19 05:41:39 CET 2024

Total time taken to generate the page: 0.01008 seconds