Hello everyone.
As I promised a few of you at last weeks meeting, I have converted the
AMPT file Michael gave me into a ROOT tree that is somewhat more readible
to the general BRAHMSer. I split the file into 4 trees of 2500 events
each, which you can find here:
/brahms/data12/scratch/bjornhs/AMPT/AmptTree_file[1-4].root
Also in the same place,
/brahms/data12/scratch/bjornhs/AMPT/evt10K-particle.dat
is the original AMPT file.(Raw text, so you can 'less' it.) The files
all-in-all contain 10k AMPT MinBias events.
To plot stuff from the ntuples, open a root session and say
root [0] TChain *ampt = new TChain("T1","Ampt Tree");
(The "T1" is important!)
root [1] ampt->Add(Form("/brahms/data12/scratch/bjornhs/AMPT/AmptTree_file*.root"));
root [2] ampt->StartViewer();
Use commands like
root [3] ampt->Draw("fParticles.fY","fParticles.fId==2212")
to plot (in this case) protons vs. rapidity. The PID is in the PDG scheme,
found on p. 206++ in the latest big Rev.
Note that the only quantities in the orig. file are the momenta, PID, and
mass of the particles. (+ imp. par, npart...) Stuff like rapidity,
transverse mass etc. are calc. by me as I fill the tree, so just in case
I made a mistake I supply the code I use for your ref:
fP = TMath::Sqrt(fPx*fPx+fPy*fPy+fPz*fPz);
fPt = TMath::Sqrt(fPx*fPx+fPy*fPy);
fE = TMath::Sqrt(fP*fP+fMass*fMass);
fY = -999;
if (fE != fPz) fY = 0.5*TMath::Log((fE + fPz)/(fE - fPz));
fMt = TMath::Sqrt(fMass*fMass+fPt*fPt);
Enjoy, and send me any questions or comments.
PS: Since I know this one is coming - no, there is unfortunalely no easy
way of turning this data into something BRAG can read. Not unless someone
has a DumpZdat-routine for ROOT...
--
Bjorn H. Samset Phone: 22856465/92051998
PhD student, heavy ion physics Adr: Schouterrassen 6
Inst. of Physics, University of Oslo 0573 Oslo
\|/
----------------------------> -*- <-----------------------------
/|\
This archive was generated by hypermail 2b30 : Mon May 27 2002 - 08:38:08 EDT