|Re: Beam profile or interaction point [message #8746 is a reply to message #8743]
||Tue, 02 June 2009 16:14
Registered: January 2009
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.
// 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???
PndReaction reaction( BeamTarget, VtxRes, "PhiPhi.root" );
FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
PndPythiaGenerator* pygen=new PndPythiaGenerator(inFile.Data());
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
[Updated on: Tue, 02 June 2009 16:32]
Report message to a moderator