Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
EIC
Project Juggler
Commits
b52349d9
Commit
b52349d9
authored
Feb 09, 2022
by
Sylvester Joosten
Browse files
Update Juggler for EICD changes
parent
916ea690
Changes
1
Hide whitespace changes
Inline
Side-by-side
JugFast/src/components/ParticlesWithTruthPID.cpp
View file @
b52349d9
...
...
@@ -61,7 +61,7 @@ public:
int
ID
=
0
;
for
(
const
auto
&
trk
:
tracks
)
{
const
eic
::
VectorXYZ
mom
=
eic
::
VectorPolar
(
1.0
/
std
::
abs
(
trk
.
qOverP
()),
trk
.
direction
().
theta
,
trk
.
direction
().
phi
);
eic
::
VectorPolar
(
1.0
/
std
::
abs
(
trk
.
qOverP
()),
trk
.
theta
()
,
trk
.
phi
()
);
const
auto
charge_rec
=
std
::
copysign
(
1.
,
trk
.
qOverP
());
// utility variables for matching
int
best_match
=
-
1
;
...
...
@@ -121,7 +121,7 @@ public:
if
(
msgLevel
(
MSG
::
DEBUG
))
{
if
(
best_match
>
0
)
{
const
auto
&
mcpart
=
mc
[
best_match
];
debug
()
<<
fmt
::
format
(
"Matched track {} with MC particle {}
\n
"
,
trk
.
ID
().
value
,
best_match
)
<<
endmsg
;
debug
()
<<
fmt
::
format
(
"Matched track {} with MC particle {}
\n
"
,
ID
,
best_match
)
<<
endmsg
;
debug
()
<<
fmt
::
format
(
" - Track: (mom: {}, theta: {}, phi: {}, charge: {})"
,
mom
.
mag
(),
mom
.
theta
(),
mom
.
phi
(),
charge_rec
)
<<
endmsg
;
...
...
@@ -130,7 +130,7 @@ public:
mcpart
.
pdgID
())
<<
endmsg
;
}
else
{
debug
()
<<
fmt
::
format
(
"Did not find a good match for track {}
\n
"
,
trk
.
ID
().
value
)
<<
endmsg
;
debug
()
<<
fmt
::
format
(
"Did not find a good match for track {}
\n
"
,
ID
)
<<
endmsg
;
debug
()
<<
fmt
::
format
(
" - Track: (mom: {}, theta: {}, phi: {}, charge: {})"
,
mom
.
mag
(),
mom
.
theta
(),
mom
.
phi
(),
charge_rec
)
<<
endmsg
;
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment