Home » PANDA » PandaRoot » Event Generators » new FTF generator in pandaroot
new FTF generator in pandaroot [message #17243] |
Thu, 04 September 2014 17:06 |
Aida Galoyan
Messages: 79 Registered: May 2007 Location: Dubna
|
continuous participant |
From: *jinr.ru
|
|
Hi all,
I have implemented new generator FTF of Geant4 in PandaRoot
in directory trunk/pgenerators.
FTF can generate Pbar-P, Pbar-Nucleus, light Anti-nucleus-Nucleus,
Hyperon-P, Hyperon-Nucleus interactions at wide energy range
from 100 MeV to 1000 GeV.
FTF is based on ideas of DPM and extended them on Pbar-Nucleus interactions.
It works faster than UrQMD.
For FTF generator, I created sub-directory FtfEvtGen similar to DpmEvtGen.
Some words, about FTF run.
-------------------------------------------------------------
At first, you need to compile FTF.
For this, in the folder pgenerators/FtfEvtGen give a command:
> make -f FTFmake
You will obtain FTFGen executable file.
For generation of Pbar-P events, you have to use file PbarP.mac .
In the file PbarP.mac, you need to set
after line
#events
Number of Events
after line
#Plab(GeV/c)
Momentum of projectile antiproton
Then, you give command
> ./FTFGen PbarP.mac
and will obtain FTF.root file with Root-Tree, where all
information about produced particles is presented.
------------------------------------------------------------------------ ---
For generation of Pbar-A events, you can use PbarA.mac.
You need to insert some values there.
First of all, Geant4 name of material - G4_Name, for example,
for carbon - G4_C, for gold - G4_Au, for cupper -G4_Cu and so on,
after line
#material
G4_Name
after line
#targetA
mass number of target nucleus
after line
#generator
ftfp or ftfb
after line
#events
Number of Events
after line
#Plab(GeV/c)
Momentum of projectile antiproton
--------------------------------------------------------------------
"ftfp" is combination of FTF generator and simple generator Preco
for low particles proceeding in nucleus.
"ftfb" is composition of FTF generator and Binary Cascade model of Geant4.
"ftfp" works faster, but "ftfb" gives more correct information about
evaporated nucleons.
------------------------------------------------------------------------ -------
Next command must be
> ./FTFGen PbarA.mac
After run, you will obtain "FTF.root" file with all information in the Root-Tree.
------------------------------------------------------------------------ ---------
We can apply "FTF.root" for future simulation in PandaRoot,
using class PndFtfGenerator similar to PndDpmGenerator.
------------------------------------------------------------------------ ----------
Information about FTF run is also presented in the file
pgenerators/FtfEvtGen/README.txt .
If you have questions about FTF, please, write me.
Good luck
Aida
|
|
|
|
Re: new FTF generator in pandaroot [message #17246 is a reply to message #17244] |
Fri, 05 September 2014 15:29 |
Aida Galoyan
Messages: 79 Registered: May 2007 Location: Dubna
|
continuous participant |
From: *jinr.ru
|
|
Hello Anastaciya,
> Did somebody else try it with apr13 external packages?
When I implemented FTF,
I used apr13 external packages in linux cluster gsi : ica.hpc.gsi.de .
SIMPATH is
/cvmfs/fairroot.gsi.de/fairsoft/apr13
About errore
------------------------------------------
Linking binary FTFGen ...
/usr/bin/ld: tmp/main.o: undefined reference to symbol 'aTouchableHistoryAllocator'
/usr/bin/ld: note: 'aTouchableHistoryAllocator' is defined in DSO /panda/fairroot/apr13Build/lib/libG4geometry.so so try adding it to the linker command line
/panda/fairroot/apr13Build/lib/libG4geometry.so: could not read symbols: Invalid operation
collect2: error: ld returned 1 exit status
make: *** [FTFGen] Error 1
------------------------------------------------------------------------ ---------
I think, if your linker requires to add
"aTouchableHistoryAllocator" from libG4geometry.so,
add in the file common.mk
-lG4geometry
in the line
LDLIBS += $(PANDALIBS) -L$(GEANT4_LIBRARY_DIR) -lG4clhep -lG4global -lG4physicslists $(GLIBS)
after -lG4clhep .
And, try to compile FTF.
If again you have problems, please, give command
>echo $GEANT4_LIBRARY_DIR
what "path" do you obtain?
If you go to this "path", can you observe files: libG4clhep.so, libG4global.so, libG4physicslists.so ?
good luck,
Aida
|
|
|
|
|
Re: new FTF generator in pandaroot [message #17297 is a reply to message #17291] |
Tue, 16 September 2014 14:33 |
Anastasia Karavdina
Messages: 76 Registered: May 2010 Location: Mainz, Germany
|
continuous participant |
From: *kph.uni-mainz.de
|
|
Dear Aida,
I had a look at main.cc file in pgenerators/FtfEvtGen and came to conclusion that elastic scattering is included, but only hadronic part of it.
Is it right?
It would be nice to have option to switch off elastic part and use only inelastic mode, as it was before in DPM.
Also I notice, that energy of final particles is stored in MeV, while pandaroot deals with GeV. I suspect due to this reason particles are treated as particles with high energy and propagation with GEANT4 (within pandaroot) is extremely slow.
I introduced conversion to GeV by changing one line in main.cc:
from
Mom.SetPxPyPzE(px,py,pz,e);
to
Mom.SetPxPyPzE(px/GeV,py/GeV,pz/GeV,e/GeV);
But I am not sure if it is correct place to do this conversion. What do you think?
Best regards,
Anastasia
|
|
|
Re: new FTF generator in pandaroot [message #17308 is a reply to message #17297] |
Wed, 17 September 2014 18:13 |
Aida Galoyan
Messages: 79 Registered: May 2007 Location: Dubna
|
continuous participant |
From: *jinr.ru
|
|
Hello Anastasia,
>I had a look at main.cc file in pgenerators/FtfEvtGen and came
>to conclusion that elastic >scattering is included, but only hadronic part of it.
>Is it right?
>It would be nice to have option to switch off elastic part and
>use only inelastic mode, as it was before in DPM.
It is not possible, because for FTF generator in PandaRoot, I apply to
Geant4 lib and can not change compilled sourse files of FTF.
But, you can obtain only inelastic events using FtfEvtGen/main.cc .
For this, you need to add in the main.cc after line
----------------------------------------------
if(std::abs(pd->GetBaryonNumber()) < 2)
----------------------------------------------
line
---------------------------------------------
if( (n>2) || ((abs(id)!=2212)&&(n==2)) ) {
------------------------------------------------
and after line
----------------------------------------------
new((*fEvt)[cnt++]) TParticle(fparticle);
----------------------------------------------
to add
--------------------
}
--------------------
After this, it is needed to compile FTF.
About "MeV" or "GeV", thank you, I didn't draw attention to it.
Just now, I have changed and commited main.cc .
Output particles are stored in GeV.
Good luck
Aida
|
|
|
|
|
|
|
|
|
Re: new FTF generator in pandaroot [message #17332 is a reply to message #17328] |
Mon, 22 September 2014 11:48 |
Aida Galoyan
Messages: 79 Registered: May 2007 Location: Dubna
|
continuous participant |
From: *dubna.ru
|
|
Hi Donghee,
The main difficulty of FTF compillation is find right path
to Geant4 library files.
In your trunk directory, please give:
> echo $GEANT4_LIBRARY_DIR
what "path" do you obtain?
If you go to this "path", do you see lib-files:
libG4clhep.so, libG4global.so, libG4physicslists.so ?
Do you observe there lib-files, proposed by Anastasia:
libG4geometry, libG4materials, libG4particles, libG4processes?
Aida
|
|
|
|
|
Re: new FTF generator in pandaroot [message #17337 is a reply to message #17327] |
Mon, 22 September 2014 23:29 |
donghee
Messages: 385 Registered: January 2009 Location: Germnay
|
first-grade participant |
From: *dip0.t-ipconnect.de
|
|
Hi Aida,
I have a question about stable particle in FTF.
In DPM, some generated particles are stable, which is Pi0, K_S0, Lambda, anti-Lambda, eta.
I am now interested in the fundamental kinematics of both generators.
At pure stand-alone generator study or fast simulation, they(stable particles) are normally turning to be unstable for DPM case.
If I want to do same study with FTFGen at fast simulation or stand-alone generator, Do I need to allow decay in order to compare with DPM generator.
What about the situation for this at FTF generator?
At full simulation, all stable particles will be decayed at GEANT level afterwards, therefore we don't need to care about it.
But generator level or fast simulation case, we have to know correctly how they are handled.
Best wishes,
Donghee
|
|
|
|
Re: new FTF generator in pandaroot [message #17339 is a reply to message #17337] |
Tue, 23 September 2014 10:58 |
Aida Galoyan
Messages: 79 Registered: May 2007 Location: Dubna
|
continuous participant |
From: *dubna.ru
|
|
Hi Donghee,
> In DPM, some generated particles are stable,which is Pi0,K_S0,
> Lambda, anti-Lambda, eta.
> I am now interested in the fundamental kinematics of both
> generators. At pure stand-alone generator study or fast
> simulation, they(stable particles) are normally turning
> to be unstable for DPM case.
Mainly, these particles are stable in DPM. It is a user
responsibility to point out what particles are unstable.
We set stable all long-lived particles -
Pi0, K0_s, K0_l, Lambda and others.
The same is true for FTF generator.
> If I want to do same study with FTFGen at fast simulation or
> stand-alone generator, Do I need to allow decay in order to compare
> with DPM generator. What about the situation for this at FTF
> generator?
The lists of "stable" particles in DPM and FTF coincide.
> At full simulation, all stable particles will be decayed at GEANT
> level afterwards, therefore we don't need to care about it. But
> generator level or fast simulation case, we have to know correctly how
> they are handled.
There is a problem in FTF how to point out that we are
going to consider rho-mesons (and others)
as stable particles. It is needed for various studies.
We implemented such possibility in DPM.
It is more complicated job in the case of FTF.
Best regards!
Aida
|
|
|
|
|
|
Goto Forum:
Current Time: Tue Oct 08 20:27:07 CEST 2024
Total time taken to generate the page: 0.00999 seconds
|