Roadmap for hadron analysis in heavy ions reactions

Documentation and macros:

documentation on particle ID: http://hades-wiki.gsi.de/cgi-bin/view/SimAna/HadesParticleID

documentation on parameter production for hadron analysis http://hades-wiki.gsi.de/cgi-bin/view/SimAna/ParameterProduction

example of hadron analysis results: http://hades-wiki.gsi.de/cgi-bin/view/SimAna/C2CNOV02Kick

macros for pidtrackcandidate qa: /d/hades/pid/PID_MACROS/HADRONS/pidtrackcand_macros

macros for pid analysis /d/hades/pid/PID_MACROS/HADRONS/pid_macros/

1. quality check of dst's:

see qa plots produced during the DST production

see macros in ./pidtrackcand_macros

check of beta-vs-momnetum 2D histos in ./pidtrackcand_macros/2DHISTOS

2. parameters for pid analysis:

"relative intensities" for pid reconstructor - from 2DHISTOS

beta-vs_momentum pdf parameters - from calculated values and 2DHISTOS for fitting of 2DHISTOS see

http://hades-wiki.gsi.de/cgi-bin/view/SimAna/ParameterProduction

and ./pid_macros/fitnfill_pvsbetaparams

existing sets for C2C and C1C, now C2C used for ArKCl

3. PID analysis - macros adapted from

Tomek's macros in

http://hades-wiki.gsi.de/cgi-bin/view/SimAna/PidMacros

for description of the procedure see

http://hades-wiki.gsi.de/cgi-bin/view/SimAna/PidImplementation

macros for hadron pid: see ./pid_macros/PID

input

dst's containing PidTrackCand

output

post dst containing PidParticle

ntuple with PID decision (dump of PidParticle)

4. calculation of PID efficiency and purity:

see ./pid_macros/count_pid_effpur

input

pid ntuple with UrQMD data (output of step 3)

output

histograms with efficency and purity

(example see in ./pid_macros/count_pid_effpur)

5. calculation of reconstruction efficiency:

step 5a - filling of ttree with needed info:

for all primary UrQMD particles shows if it was reconstructed

by selected trcaking algorithm

see ./pid_macros/acceptance_ttree

input: UrQMD dst

output: ttree

step 5b - filling of efficiency histograms from ttree

see ./pid_macros/acceptance

6. filling weights to pid ntuple

(weight = PID purity/(PID efficiency * reconstruction efficiency)

see ./pid_macros/weights

input

pid ntuple

root file with PID efficiency and purity

root file with reconstruction efficiency

output

ttree with identical content as pid ntuple and correct weights

7. filling of macros for identified particled: theta, momentum, rapidity,

m_t distributions

see ./pid_macros/refilled

input

pid ttree

output

root file with histograms

8. normalization to one event

./pid_macros/refilled/pidplots/fill_spectra_sec1.C

9. plotting of results

-- PavelTlusty - 27 Jun 2007