FazBrowse GitHub Viewer | Trending |
URL:
| Home
Tools: [Download Repo ZIP]   [Original HTTPS Page]

Provisionary fix for ITS fit with B=0 by shahor02 · Pull Request #10261 · AliceO2Group/AliceO2 · GitHub

/ AliceO2 Public

Provisionary fix for ITS fit with B=0 #10261

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
shahor02 wants to merge 1 commit into from
Conversation 3 Commits 1 Checks 0 Files changed
Closed

Provisionary fix for ITS fit with B=0 #10261

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension .cxx  (1) All 1 file type selected
Viewed files
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Unified
Split
Hide whitespace
Diff view
Unified
Split
Hide whitespace
38 changes: 23 additions & 15 deletions Detectors/ITSMFT/ITS/tracking/src/TrackerTraits.cxx
Show comments
View file
Open in desktop
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters. Learn more about bidirectional Unicode characters
Show hidden characters
Original file line number Diff line number Diff line change
Expand Up @@ -912,7 +912,8 @@ bool TrackerTraits::trackFollowing(TrackITSExt* track, int rof, bool outward, co
/// a TrackParCov
track::TrackParCov TrackerTraits::buildTrackSeed(const Cluster& cluster1, const Cluster& cluster2, const Cluster& cluster3, const TrackingFrameInfo& tf3, float resolution)
{
const float ca = std::cos(tf3.alphaTrackingFrame), sa = std::sin(tf3.alphaTrackingFrame);
float ca, sa;
o2::math_utils::sincos(tf3.alphaTrackingFrame, sa, ca);
const float x1 = cluster1.xCoordinate * ca + cluster1.yCoordinate * sa;
const float y1 = -cluster1.xCoordinate * sa + cluster1.yCoordinate * ca;
const float z1 = cluster1.zCoordinate;
Expand All @@ -922,28 +923,35 @@ track::TrackParCov TrackerTraits::buildTrackSeed(const Cluster& cluster1, const
const float x3 = tf3.xTrackingFrame;
const float y3 = tf3.positionTrackingFrame[0];
const float z3 = tf3.positionTrackingFrame[1];

const float crv = math_utils::computeCurvature(x1, y1, x2, y2, x3, y3);
const float x0 = math_utils::computeCurvatureCentreX(x1, y1, x2, y2, x3, y3);
const float tgl12 = math_utils::computeTanDipAngle(x1, y1, x2, y2, z1, z2);
const float tgl23 = math_utils::computeTanDipAngle(x2, y2, x3, y3, z2, z3);

const float s2 = resolution;
const float fy = 1. / (cluster2.radius - cluster3.radius);
const float& tz = fy;
float cy = 1.e15f;
if (std::abs(getBz()) > o2::constants::math::Almost0) {
cy = (math_utils::computeCurvature(x1, y1, x2, y2 + resolution, x3, y3) - crv) /
(resolution * getBz() * o2::constants::math::B2C) *
20.f; // FIXME: MS contribution to the cov[14] (*20 added)

bool withField = std::abs(getBz()) > o2::constants::math::Almost0;
if (!withField) {
return track::TrackParCov(tf3.xTrackingFrame, tf3.alphaTrackingFrame,
{y3, z3, std::sin(o2::math_utils::atan2(y3 - y1, x3 - x1)), 0.5f * (tgl12 + tgl23), 1.f / o2::track::kMostProbablePt},
{s2,
0.f, s2,
0.f, 0.f, s2 * fy,
0.f, 0.f, 0.f, s2 * tz * tz,
0.f, 0.f, 0.f, 0.f, o2::track::kMostProbablePt * o2::track::kMostProbablePt});
}
const float s2 = resolution;

const float crv = math_utils::computeCurvature(x1, y1, x2, y2, x3, y3);
const float x0 = math_utils::computeCurvatureCentreX(x1, y1, x2, y2, x3, y3);

const float cy = (math_utils::computeCurvature(x1, y1, x2, y2 + resolution, x3, y3) - crv) / (resolution * getBz() * o2::constants::math::B2C) * 20.f; // FIXME: MS contribution to the cov[14] (*20 added)
return track::TrackParCov(tf3.xTrackingFrame, tf3.alphaTrackingFrame,
{y3, z3, crv * (x3 - x0), 0.5f * (tgl12 + tgl23),
std::abs(getBz()) < o2::constants::math::Almost0 ? 1.f / o2::track::kMostProbablePt
: crv / (getBz() * o2::constants::math::B2C)},
{s2, 0.f, s2, s2 * fy, 0.f, s2 * fy * fy, 0.f, s2 * tz, 0.f, s2 * tz * tz, s2 * cy, 0.f,
s2 * fy * cy, 0.f, s2 * cy * cy});
withField ? 1.f / o2::track::kMostProbablePt : crv / (getBz() * o2::constants::math::B2C)},
{s2,
0.f, s2,
s2 * fy, 0.f, s2 * fy * fy,
0.f, s2 * tz, 0.f, s2 * tz * tz, s2 * cy,
0.f, s2 * fy * cy, 0.f, s2 * cy * cy});
}

void TrackerTraits::traverseCellsTree(const int currentCellId, const int currentLayerId)
Expand Down
Back | FazBrowse Home | New Git URL