Home » PANDA » PandaRoot » Event Generators » PYTHIA6 usage
PYTHIA6 usage [message #23646] |
Thu, 21 February 2019 12:25 |
Anna Skachkova
Messages: 37 Registered: September 2018
|
continuous participant |
From: *jinr.ru
|
|
Dear Ralf and Tobias,
Turning back again to PYTHIA. I haven't use Pandaroot for several years. And now I see that there were changed a lot.
The way which I used before (just reading by PandaRoot the ntuple, generated by Pythia6 - see attached run_sim_sttcombi_DUB.C file) doesn't work any more. I tryed to do the same in the modern format - see tut_sim_my.C. I doesn't work as well. I found that even the files PndPythia6Generator.h and PndPythia6Generator.cxx do not present in /pgenerators in the current version of PandaRoot.
So the main questions are:
1. How to launch PandaRoot in the current version with the Pythia6 as prepared external ntuple (in the way I did it before)? It would be very appreciated if you could send already working example. And it would be good to renew the correspondent page in tutorial correspondingly.
2. It would also be very good to have a possibility to launch PYTHIA directly from Pandaroot in order to have a chance to make a preselection of events which is crucial important for background simulation. The problem is that I do not understand how to switch between a lot of internal Pythia parameters (the main are the used subprocesses which are supposed for simulation). The only information that I found is http://www.hep.by/root/tutorials/pythia/pythiaExample.C.html which for me is absolutely not enough for understanding.
Thank you in advance,
with the best regards,
Anna
On 20.02.2019 18:12, RalfKliemt wrote:
> Dear Anna,
>
> I found a few lines which I fixed. Please go on and test the macro in the attachment.
>
> Cheers!
> Ralf
>
|
|
|
Re: PYTHIA6 usage [message #23652 is a reply to message #23646] |
Thu, 21 February 2019 12:27 |
Anna Skachkova
Messages: 37 Registered: September 2018
|
continuous participant |
From: *jinr.ru
|
|
Dear Ralf,
Thank you for a quick replay. Yes, fort.16 is an event file in txt format.
I do not fully understand - what do you mean by "we have to add the generator to the Master Run". How to do it?
For the moment I've got the error :
/mnt/hard/work_pandaroot2019/Simulation/DYee/Sig/10000$ root -l -b -q tut_sim_my_Ralf.C
Processing tut_sim_my_Ralf.C...
In file included from input_line_10:1:
/mnt/hard/work_pandaroot2019/Simulation/DYee/Sig/10000/tut_sim_my_Ralf.C :59:9: error: no matching member function for call to 'SetGenerator'
fRun->SetGenerator(primGen); //[R.K.] we have to add the generator to the Master Run
~~~~~~^~~~~~~~~~~~
G__MasterDict dictionary payload:769:8: note: candidate function not viable: no known conversion from 'FairFilteredPrimaryGenerator *' to 'FairGenerator *' for 1st argument
void SetGenerator(FairGenerator *boxGen);
^
G__MasterDict dictionary payload:778:8: note: candidate function not viable: no known conversion from 'FairFilteredPrimaryGenerator *' to 'PndBoxGenerator *' for 1st argument
void SetGenerator(PndBoxGenerator *boxGen);
^
G__MasterDict dictionary payload:753:8: note: candidate function not viable: requires 0 arguments, but 1 was provided
void SetGenerator();
^
Thank you,
Anna
[Updated on: Thu, 21 February 2019 12:29] Report message to a moderator
|
|
|
|
|
|
|
|
|
Re: PYTHIA6 usage [message #25190 is a reply to message #25147] |
Thu, 08 October 2020 14:23 |
Anna Skachkova
Messages: 37 Registered: September 2018
|
continuous participant |
From: *jinr.ru
|
|
Dear Ralf,
Thank you for a fast replay. I have tried your script with the different changes, but didn't manage to get the same results. With the exactly(!) your script I've got:
mostly
......===== PndPidCorrelator - Event: 960 - Number of tracks for pid 0 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 961 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 962 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 963 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 964 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 965 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 1 FSC: 0
===== PndPidCorrelator - Event: 966 - Number of tracks for pid 0 - Number of Clusters for pid: EMC: 0 FSC: 2
===== PndPidCorrelator - Event: 967 - Number of tracks for pid 1 - Number of Clusters for pid: EMC: 2 FSC: 0
....
that means, as far as I understand, there are mostly the only one track per event, that is confirmed by the picture from root.
From the picture Pdg.png (MCTrack.fPdgCode) it clear seen the presence of the fist positive muon in 90% and in 10% mu-, while a priory
mu+ and mu- are absolutely indentical in the file, having in account that there are also the other charged particles in at least >50% events.
I'm very confused in this situation. I would be very grateful if anybody could give a hint what I'm doing wrong.
Thanks in advance,
Anna
-
Attachment: Pdg.png
(Size: 0.10KB, Downloaded 209 times)
|
|
|
|
|
|
|
Re: PYTHIA6 direct usage [message #25203 is a reply to message #25200] |
Fri, 09 October 2020 09:39 |
AnnaSkachkova
Messages: 9 Registered: April 2013 Location: Dubna
|
occasional visitor |
From: 31.131.194*
|
|
A, yes, here they are :
===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0 0x00007fe8a55688fa in __GI___waitpid (pid=27576, stat_loc=stat_loc
entry=0x7ffe2cd38098, options=options
entry=0) at ../sysdeps/unix/sysv/linux/waitpid.c:29
#1 0x00007fe8a54d79cb in do_system (line=<optimized out>) at ../sysdeps/posix/system.c:148
#2 0x00007fe8a616a2cd in TUnixSystem::Exec (shellcmd=<optimized out>, this=0x5628b46d07c0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:2119
#3 TUnixSystem::StackTrace (this=0x5628b46d07c0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:2413
#4 0x00007fe8a616cc04 in TUnixSystem::DispatchSignals (this=0x5628b46d07c0, sig=kSigSegmentationViolation) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:3644
#5 <signal handler called>
#6 0x00007fe88cf0cee5 in ROOT::delete_TMCParticle (p=0x9999999999999999) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/monte carlo/pythia6/G__EGPythia6.cxx:270
#7 0x00007fe8a610b631 in TClass::Destructor (this=this
entry=0x5628b7c0edc0, obj=obj
entry=0x9999999999999999, dtorOnly=dtorOnly
entry=false) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/meta/src/TClass .cxx:5198
#8 0x00007fe8a3919fb6 in TBufferFile::ReadFastArray (this=0x5628b9ff8570, start=<optimized out>, cl=0x5628b7c0edc0, n=<optimized out>, isPreAlloc=<optimized out>, streamer=<optimized out>, onFileClass=0x0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:1660
#9 0x00007fe8a3b8c9ac in TStreamerInfo::ReadBuffer<char**> (this=0x5628b7a32c70, b=..., arr=
0x7ffe2cd3ab70: 0x7ffe2cd3ab68, compinfo=0x5628bd0e4ae8, first=first
entry=0, last=last
entry=1, narr=1, eoffset=0, arrayMode=0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoReadBuffer.cxx:1080
#10 0x00007fe8a39e1743 in TStreamerInfoActions::GenericReadAction (buf=..., addr=<optimized out>, config=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoActions.cxx:185
#11 0x00007fe8a392143e in TStreamerInfoActions::TConfiguredAction::operator() (this=<optimized out>, this=<optimized out>, object=0x5628bd0e3700, buffer=...) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/inclu de/TStreamerInfoActions.h:124
#12 TBufferFile::ApplySequence (obj=0x5628bd0e3700, sequence=..., this=0x5628b9ff8570) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3564
#13 TBufferFile::ReadClassBuffer (this=0x5628b9ff8570, cl=0x5628b64394e0, pointer=0x5628bd0e3700, onFileClass=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3485
#14 0x00007fe895307863 in PndPythia6Direct::Streamer (this=0x5628bd0e3700/build/gdb-UN1U6s/gdb-8.0.1/gdb/stack.c:627: internal-error: void print_frame_args(symbol*, frame_info*, int, ui_file*): Assertion `nsym != NULL' failed.
===========================================================
The lines below might hint at the cause of the crash.
You may get help by asking at the ROOT forum http://root.cern.ch/forum
Only if you are really convinced it is a bug in ROOT then please submit a
report at http://root.cern.ch/bugs Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#6 0x00007fe88cf0cee5 in ROOT::delete_TMCParticle (p=0x9999999999999999) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/monte carlo/pythia6/G__EGPythia6.cxx:270
#7 0x00007fe8a610b631 in TClass::Destructor (this=this
entry=0x5628b7c0edc0, obj=obj
entry=0x9999999999999999, dtorOnly=dtorOnly
entry=false) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/meta/src/TClass .cxx:5198
#8 0x00007fe8a3919fb6 in TBufferFile::ReadFastArray (this=0x5628b9ff8570, start=<optimized out>, cl=0x5628b7c0edc0, n=<optimized out>, isPreAlloc=<optimized out>, streamer=<optimized out>, onFileClass=0x0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:1660
#9 0x00007fe8a3b8c9ac in TStreamerInfo::ReadBuffer<char**> (this=0x5628b7a32c70, b=..., arr=
0x7ffe2cd3ab70: 0x7ffe2cd3ab68, compinfo=0x5628bd0e4ae8, first=first
entry=0, last=last
entry=1, narr=1, eoffset=0, arrayMode=0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoReadBuffer.cxx:1080
#10 0x00007fe8a39e1743 in TStreamerInfoActions::GenericReadAction (buf=..., addr=<optimized out>, config=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoActions.cxx:185
#11 0x00007fe8a392143e in TStreamerInfoActions::TConfiguredAction::operator() (this=<optimized out>, this=<optimized out>, object=0x5628bd0e3700, buffer=...) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/inclu de/TStreamerInfoActions.h:124
#12 TBufferFile::ApplySequence (obj=0x5628bd0e3700, sequence=..., this=0x5628b9ff8570) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3564
#13 TBufferFile::ReadClassBuffer (this=0x5628b9ff8570, cl=0x5628b64394e0, pointer=0x5628bd0e3700, onFileClass=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3485
#14 0x00007fe895307863 in PndPythia6Direct::Streamer (this=0x5628bd0e3700/build/gdb-UN1U6s/gdb-8.0.1/gdb/stack.c:627: internal-error: void print_frame_args(symbol*, frame_info*, int, ui_file*): Assertion `nsym != NULL' failed.
===========================================================
Root > *** Error in `/home/anna/WORK/PandaRoot/fairsoft_jun19/build/bin/root.exe': double free or corruption (!prev): 0x00005628ba08f990 ***
Thanks,
Anna
|
|
|
|
Re: PYTHIA6 direct usage [message #25231 is a reply to message #25218] |
Mon, 19 October 2020 15:33 |
Anna Skachkova
Messages: 37 Registered: September 2018
|
continuous participant |
From: *jinr.ru
|
|
Dear Ralf!
Thank you for your hint. I have tried different variations of comments, but unfortunately it didn't help. Everything seems the same. Or I did something wrong?
-I- FairRunTimeDB::InitContainer() FairGeoParSet
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CountLevels>: max level = 12, max placements = 4550
Info in <TGeoManager::CloseGeometry>: 2028308 nodes/ 5654 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Container FairGeoParSet initialized from ROOT file.
-I- FairRunTimeDB::InitContainer() FairBaseParSet
PndPythia6Direct::SetPandaDefaults()
*** Break *** segmentation violation
Aborted (core dumped)
===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0 0x00007fa6baf4d8fa in __GI___waitpid (pid=10010, stat_loc=stat_loc
entry=0x7ffd06de7158, options=options
entry=0) at ../sysdeps/unix/sysv/linux/waitpid.c:29
#1 0x00007fa6baebc9cb in do_system (line=<optimized out>) at ../sysdeps/posix/system.c:148
#2 0x00007fa6bbb4f2cd in TUnixSystem::Exec (shellcmd=<optimized out>, this=0x55cc9e4677c0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:2119
#3 TUnixSystem::StackTrace (this=0x55cc9e4677c0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:2413
#4 0x00007fa6bbb51c04 in TUnixSystem::DispatchSignals (this=0x55cc9e4677c0, sig=kSigSegmentationViolation) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/unix/src/TUnixS ystem.cxx:3644
#5 <signal handler called>
#6 0x00007fa6a28f1ee5 in ROOT::delete_TMCParticle (p=0x9999999999999999) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/monte carlo/pythia6/G__EGPythia6.cxx:270
#7 0x00007fa6bbaf0631 in TClass::Destructor (this=this
entry=0x55cca19a64e0, obj=obj
entry=0x9999999999999999, dtorOnly=dtorOnly
entry=false) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/meta/src/TClass .cxx:5198
#8 0x00007fa6b92fefb6 in TBufferFile::ReadFastArray (this=0x55cca3da61d0, start=<optimized out>, cl=0x55cca19a64e0, n=<optimized out>, isPreAlloc=<optimized out>, streamer=<optimized out>, onFileClass=0x0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:1660
#9 0x00007fa6b95719ac in TStreamerInfo::ReadBuffer<char**> (this=0x55cca17ca6a0, b=..., arr=
0x7ffd06de9c30: 0x7ffd06de9c28, compinfo=0x55cca6e911a8, first=first
entry=0, last=last
entry=1, narr=1, eoffset=0, arrayMode=0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoReadBuffer.cxx:1080
#10 0x00007fa6b93c6743 in TStreamerInfoActions::GenericReadAction (buf=..., addr=<optimized out>, config=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoActions.cxx:185
#11 0x00007fa6b930643e in TStreamerInfoActions::TConfiguredAction::operator() (this=<optimized out>, this=<optimized out>, object=0x55cca6e8fdc0, buffer=...) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/inclu de/TStreamerInfoActions.h:124
#12 TBufferFile::ApplySequence (obj=0x55cca6e8fdc0, sequence=..., this=0x55cca3da61d0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3564
#13 TBufferFile::ReadClassBuffer (this=0x55cca3da61d0, cl=0x55cca01cfe30, pointer=0x55cca6e8fdc0, onFileClass=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3485
#14 0x00007fa6aacec863 in PndPythia6Direct::Streamer (this=0x55cca6e8fdc0/build/gdb-UN1U6s/gdb-8.0.1/gdb/stack.c:627: internal-error: void print_frame_args(symbol*, frame_info*, int, ui_file*): Assertion `nsym != NULL' failed.
===========================================================
The lines below might hint at the cause of the crash.
You may get help by asking at the ROOT forum http://root.cern.ch/forum
Only if you are really convinced it is a bug in ROOT then please submit a
report at http://root.cern.ch/bugs Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#6 0x00007fa6a28f1ee5 in ROOT::delete_TMCParticle (p=0x9999999999999999) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/monte carlo/pythia6/G__EGPythia6.cxx:270
#7 0x00007fa6bbaf0631 in TClass::Destructor (this=this
entry=0x55cca19a64e0, obj=obj
entry=0x9999999999999999, dtorOnly=dtorOnly
entry=false) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/core/meta/src/TClass .cxx:5198
#8 0x00007fa6b92fefb6 in TBufferFile::ReadFastArray (this=0x55cca3da61d0, start=<optimized out>, cl=0x55cca19a64e0, n=<optimized out>, isPreAlloc=<optimized out>, streamer=<optimized out>, onFileClass=0x0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:1660
#9 0x00007fa6b95719ac in TStreamerInfo::ReadBuffer<char**> (this=0x55cca17ca6a0, b=..., arr=
0x7ffd06de9c30: 0x7ffd06de9c28, compinfo=0x55cca6e911a8, first=first
entry=0, last=last
entry=1, narr=1, eoffset=0, arrayMode=0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoReadBuffer.cxx:1080
#10 0x00007fa6b93c6743 in TStreamerInfoActions::GenericReadAction (buf=..., addr=<optimized out>, config=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TStreamerI nfoActions.cxx:185
#11 0x00007fa6b930643e in TStreamerInfoActions::TConfiguredAction::operator() (this=<optimized out>, this=<optimized out>, object=0x55cca6e8fdc0, buffer=...) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/build_for_fair/inclu de/TStreamerInfoActions.h:124
#12 TBufferFile::ApplySequence (obj=0x55cca6e8fdc0, sequence=..., this=0x55cca3da61d0) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3564
#13 TBufferFile::ReadClassBuffer (this=0x55cca3da61d0, cl=0x55cca01cfe30, pointer=0x55cca6e8fdc0, onFileClass=<optimized out>) at /home/anna/WORK/PandaRoot/fairsoft_jun19/tools/root/io/io/src/TBufferFil e.cxx:3485
#14 0x00007fa6aacec863 in PndPythia6Direct::Streamer (this=0x55cca6e8fdc0/build/gdb-UN1U6s/gdb-8.0.1/gdb/stack.c:627: internal-error: void print_frame_args(symbol*, frame_info*, int, ui_file*): Assertion `nsym != NULL' failed.
===========================================================
Root > *** Error in `/home/anna/WORK/PandaRoot/fairsoft_jun19/build/bin/root.exe': double free or corruption (!prev): 0x000055cca3e3d090 ***
Best regards,
Anna
|
|
|
Goto Forum:
Current Time: Sun Oct 13 13:05:03 CEST 2024
Total time taken to generate the page: 0.00861 seconds
|