| Re: Beam profile or interaction point [message #8746 is a reply to message #8743] |
Tue, 02 June 2009 16:14   |
donghee
Messages: 385
Registered: January 2009
Location: Germnay
|
first-grade participant |
From: *kph.uni-mainz.de
|
|
Thank you for your helpful answer.
The script looks quite useful in any case.
Actually, my output file is saved with primary vertex position (0,0,0) in the ascii file.
I call this file in the macro level during the simulation run.
There are two choice.
I can make the smearing in the fortran code of generator, this is not convenient,
or I can do it in the panda macro level.
I prefer the latter.
To do that, I have to keep momentum information from ascii file and just put smearing vertex position as you recommended.
| Quote: |
// mean location of the target
Vertex.SetXYZT(0,0,0,0); // zero by default
// smearing of the primary vertex (Gauss with the above mean, and this sigma)
TVector3 VtxRes(.1,.1,.2); // what target smearing do you want???
TParticle BeamTarget(999,1,0,0,0,0,s,Vertex);
PndReaction reaction( BeamTarget, VtxRes, "PhiPhi.root" );
|
| Quote: |
FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
fRun->SetGenerator(primGen);
PndPythiaGenerator* pygen=new PndPythiaGenerator(inFile.Data());
primGen->AddGenerator(pygen);
fRun->SetStoreTraj(kTRUE);
|
I think that I have to combine some function of PndReaction with PndPythiaGenerator.
I'm not clear what do I have to as a next step?
Could you give me some hint?
Thank you! Aida
Best regards,
Donghee
[Updated on: Tue, 02 June 2009 16:32] Report message to a moderator |
|
|
|
GSI Forum
Search
Help
Members
Login
Home
