Hi! I am trying to simulate the current signal caused by the alpha particle incident on a gas medium sensor, and I have generated gas files for 85% Ar and 15% CO2 through the generate.c.
The TrackSrim file I imported is Alpha_in-Ar.txt, and I encountered the issue shown in the above picture.
Hi,
I can think of two explanations off-hand:
- Did you perhaps forget to associate the gas medium to the ComponentAnalyticField object?
- Or is the starting point of your track outside the active area?
Thank you for your answer. I will check if it is caused by one of these two reasons.
Thanks for your reply! I have changed the starting point to:
The output is:
But I still can
t see any current signal output.  New users cant post links so I put the screenshot of my code here.
Hi,
GetElectron is not really meaningful for TrackTrim (all electrons in a cluster are supposed to be at the same location, xc, yc, zc).
You can change the inner for loop to
for (int k = 0; k < nc; ++k) {
drift.DriftElectron(xc, yc, zc, tc);
}
By the way, the format of the GetCluster function that you are using is still supported, but I would like to phase it out. Instead of the while loop you can do:
for (const auto& cluster : track.GetClusters()) {
for (int k = 0; k < cluster.n; ++k) {
drift.DriftElectron(cluster.x, cluster.y, cluster.z, cluster.t);
}
}
Finally, with DriftLineRKF you will always get the same drift line for a given starting point (there are no fluctuations in the drift path).
So, an even more efficient way would be something like this:
for (const auto& cluster : track.GetClusters()) {
drift.SetElectronSignalScalingFactor(cluster.n);
drift.DriftElectron(cluster.x, cluster.y, cluster.z, cluster.t);
}
See also this example:
Thanks a lot.