TPC dE/dx [message #11557] |
Tue, 08 March 2011 15:54 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: *to.infn.it
|
|
Dear all,
I have tried to write a "quick and dirty" algorithm to calculate TPC dE/dx with a truncated mean -> pid/PidCorr/PndPidTpcInfo.cxx
The energy is taken from PndTpcCluster::amp(), the dx is calculated somehow from the distance of consecutive hits, as a first approximation (in reality it is the sum of the semi-distance with the previous hit and the semi-distance with the next hit).
This is the plot I obtain, after 10k events * (p + e + pi + mu + k) with the box generator:
If I compare the plot to what present in the TPR, page 142 of the pdf, I have some questions:
a) What is the unity of PndTpcCluster::amp()? If I assume it is in GeV as the pandaroot standard, I obtain a around MeV/cm, while I suppose it should be something like keV/cm. And the numbers are still different. (I am using Geant3+Alice settings) Maybe the gas mixture was changed?
b) the plot is worse than the TPR one, because it is using reconstruction and not MC, but it looks like electrons are less separated now than the old plot (where they were quite far). Here they are overlapping with a bit higher de/dx, you are not able to find them like in the TPC plot. This is not clear to me.
Comments from experts are welcome.
|
|
|