Skip to content
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

Jump to bottom

Tack track curvature directly from Delphes and don't recompute it from pT #112

Merged
merged 5 commits into from
Nov 2, 2023
Merged

Tack track curvature directly from Delphes and don't recompute it from pT #112

Changes from 1 commit
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
e60aa2d
Update of convertTrack in DelphesEDM4HepConverter.cc: take the track'…
EmanuelPerez Oct 9, 2023
db82421
Remove no longer used code
tmadlener Nov 2, 2023
5eefb39
[format] Fix formatting
tmadlener Nov 2, 2023
5b2b2ee
Remove no longer necessary passing of magnetic field
tmadlener Nov 2, 2023
48eaa70
Remove no longer necessary constant
tmadlener Nov 2, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 8 additions & 0 deletions converter/src/DelphesEDM4HepConverter.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -556,6 +556,11 @@ namespace k4SimDelphes {
trackState.phi = cand->Phi;
// Same thing under different name in Delphes
trackState.tanLambda = cand->CtgTheta;

/*
// Instead of recomputing the track's omega from its pT, better take it directly
// from the curvature of the delphes track - see below.

// Only do omega when there is actually a magnetic field.
double varOmega = 0;
if (magFieldBz) {
Expand All @@ -567,6 +572,7 @@ namespace k4SimDelphes {
// constant B-field -> relative error on pT is relative error on omega
varOmega = cand->ErrorPT * cand->ErrorPT / cand->PT / cand->PT * trackState.omega * trackState.omega;
}
*/
tmadlener marked this conversation as resolved.
Show resolved Hide resolved

// fill the covariance matrix. There is a conversion of units
// because the covariance matrix in delphes is with the original units
Expand All @@ -579,6 +585,8 @@ namespace k4SimDelphes {
// relating to curvature (index 2)
double scale2 = -2.; // CAREFUL: DELPHES USES THE HALF-CURVATURE

trackState.omega = cand->C * scale2 ;
tmadlener marked this conversation as resolved.
Show resolved Hide resolved

covMatrix[0] = covaFB(0, 0);

covMatrix[1] = covaFB(1, 0);
Expand Down
Loading