Home » PANDA » PandaRoot » Analysis » Fast simulation - MCTrack block
Fast simulation - MCTrack block [message #16133] |
Tue, 01 April 2014 10:38 |
Elisabetta Prencipe (2)
Messages: 214 Registered: February 2013
|
first-grade participant |
From: *ikp.kfa-juelich.de
|
|
Dear Klaus,
when I run the sim macro from the folder /macro/scrut/ :
root [0] .x simfast.C("outputDs2535","Y2535DstarK.dec",9.38,10000,"pbarpSystem")
something funny happens, if I check the plot of MCTrack.fPx and MCTRack.fPy. I got a spike on 0, not a distribution.
It is not so for MCTrack.fPz and MCTrack.E. Do I misinterpret something here, or what is the explanation for this?
Best regards, Elisabetta
|
|
|
|
|
Re: Fast simulation - MCTrack block [message #16142 is a reply to message #16138] |
Tue, 01 April 2014 11:36 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: *dip0.t-ipconnect.de
|
|
You don't have enough energy, and EvtGen is not letting your pbarpSystem decay:
This top particle is pbarpSystem 4.410 (10.365,0.000,-0.000,9.380)
pbarpSystem -> D'_s1+ D_s- 2.525 (2.525,0.000,0.000,0.000) 3; 1.968 (1.968,0.000,0.000,0.000)
With your momentum you have sqrt(s) = 4.410, but the mass sum of D'_s and D_s is 2.525+1.968 = 4.493. This is the reason why you have such error messages from evtgen.
You should increase your momentum.
|
|
|
|
|
|
|
|
|
|
|
Re: Fast simulation - MCTrack block [message #16159 is a reply to message #16157] |
Tue, 01 April 2014 14:41 |
Elisabetta Prencipe (2)
Messages: 214 Registered: February 2013
|
first-grade participant |
From: *ikp.kfa-juelich.de
|
|
Hi Stefano and Klaus,
that was just a typo this morning (my apologize for that).
In the meantime, I updated the trunk to rev-24270.
when I run:
root [0] .x simfast.C("outputDs2535","Y2535DstarK.dec",9.83,1000,"pbarpSystem0")
[pbarpsystem0 now is everywhere]
all fX,Y,Z,E distribution are a spike to 0. Nothing passes.
I see that in the new simfast macro a new entry is added:
void simfast(TString Prefix, TString Decfile, Float_t Mom, Int_t nEvents = 1000, TString Resonance="pbarpSystem0", int pdgcode = 11 )
what shall I write instead of 11? If I write 88888, it does not work. I get only "red" warning of non existing particles, when I try this macro.
If I do not write anything, it does not work in any case. If I leave 11, same situation.
FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
-W FairPrimaryGenerator: PDG code 88880 not found in database. This warning can be savely ignored.
[this for 1000 times..at least yesterday it looked the normal staff]
Any idea what's wrong here? how shall I run my analysis with the fast simulation tool? It is still not clear to me....
thank you, Elisabetta
|
|
|
|
Re: Fast simulation - MCTrack block [message #16161 is a reply to message #16159] |
Tue, 01 April 2014 14:50 |
Klaus Götzen
Messages: 293 Registered: June 2006 Location: GSI
|
first-grade participant |
From: *gsi.de
|
|
Hi Elisabetta,
you can ignore the pdgcode=11 for EvtGen. It's the particle type for running Box generator, in case you want to do that. I introduced it to not modify always the macro for running particle-gun for different species. (It's actually described in the README like this.)
This error from FairPrimaryGenerator is normal because TDatabasePDG does not know about particles defined for EvtGen. However, myself I'm looking for a switch to turn it off. If you want to immediately get rid of it, you could also add
double sqrts = ...
TDatabasePDG::Instance()->AddParticle("pbarpSystem0","ppbar0",sqrts,0,0.0001,0,"",88880);
TDatabasePDG::Instance()->AddParticle("pbarpSystem","ppbar",sqrts,0,0.0001,0,"",88888);
in the beginning of simfast.C.
Best,
Klaus
|
|
|
|
|
|
|
|
|
|
|
Goto Forum:
Current Time: Sun Oct 06 12:37:39 CEST 2024
Total time taken to generate the page: 0.00874 seconds
|