Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
K
km3io
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Package registry
Container Registry
Model registry
Operate
Environments
Terraform modules
Monitor
Incidents
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
km3py
km3io
Commits
30394b8b
Commit
30394b8b
authored
5 years ago
by
Tamas Gal
Browse files
Options
Downloads
Patches
Plain Diff
Simplify caching
parent
9921ac8f
No related branches found
Branches containing commit
No related tags found
Tags containing commit
1 merge request
!21
Resolve "Reduce the amount of uproot.open (to one)"
Pipeline
#9167
passed with warnings
5 years ago
Stage: test
Stage: coverage
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
km3io/offline.py
+57
-86
57 additions, 86 deletions
km3io/offline.py
with
57 additions
and
86 deletions
km3io/offline.py
+
57
−
86
View file @
30394b8b
import
functools
import
uproot
import
numpy
as
np
import
warnings
...
...
@@ -303,7 +302,7 @@ class Reader:
def
__repr__
(
self
):
return
"
<{}: {} entries>
"
.
format
(
self
.
__class__
.
__name__
,
len
(
self
))
@property
@
cached_
property
def
keys
(
self
):
"""
wrapper for all keys in an offline file.
...
...
@@ -312,9 +311,7 @@ class Reader:
Class
OfflineKeys.
"""
if
self
.
_keys
is
None
:
self
.
_keys
=
OfflineKeys
(
self
.
_tree
)
return
self
.
_keys
return
OfflineKeys
(
self
.
_tree
)
class
OfflineReader
:
...
...
@@ -329,15 +326,6 @@ class OfflineReader:
path-like object that points to the file of ineterst.
"""
self
.
_file_path
=
file_path
self
.
_events
=
None
self
.
_hits
=
None
self
.
_tracks
=
None
self
.
_mc_hits
=
None
self
.
_mc_tracks
=
None
self
.
_keys
=
None
self
.
_best_reco
=
None
self
.
_header
=
None
self
.
_usr
=
None
if
file_path
is
not
None
:
self
.
_tree
=
uproot
.
open
(
self
.
_file_path
)[
MAIN_TREE_NAME
]
...
...
@@ -358,19 +346,18 @@ class OfflineReader:
def
__len__
(
self
):
return
len
(
self
.
_data
)
@property
@
cached_
property
def
header
(
self
):
if
self
.
_header
is
None
:
fobj
=
uproot
.
open
(
self
.
_file_path
)
if
'
Head
'
in
fobj
:
self
.
_header
=
{}
for
n
,
x
in
fobj
[
'
Head
'
].
_map_3c_string_2c_string_3e_
.
items
():
self
.
_header
[
n
.
decode
(
"
utf-8
"
)]
=
x
.
decode
(
"
utf-8
"
).
strip
()
else
:
warnings
.
warn
(
"
Your file header has an unsupported format
"
)
return
self
.
_header
@property
fobj
=
uproot
.
open
(
self
.
_file_path
)
if
'
Head
'
in
fobj
:
header
=
{}
for
n
,
x
in
fobj
[
'
Head
'
].
_map_3c_string_2c_string_3e_
.
items
():
header
[
n
.
decode
(
"
utf-8
"
)]
=
x
.
decode
(
"
utf-8
"
).
strip
()
return
header
else
:
warnings
.
warn
(
"
Your file header has an unsupported format
"
)
@cached_property
def
keys
(
self
):
"""
wrapper for all keys in an offline file.
...
...
@@ -379,11 +366,9 @@ class OfflineReader:
Class
OfflineKeys.
"""
if
self
.
_keys
is
None
:
self
.
_keys
=
OfflineKeys
(
self
.
_tree
)
return
self
.
_keys
return
OfflineKeys
(
self
.
_tree
)
@property
@
cached_
property
def
events
(
self
):
"""
wrapper for offline events.
...
...
@@ -392,13 +377,11 @@ class OfflineReader:
Class
OfflineEvents.
"""
if
self
.
_events
is
None
:
self
.
_events
=
OfflineEvents
(
self
.
keys
.
cut_events_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
events_keys
])
return
self
.
_events
return
OfflineEvents
(
self
.
keys
.
cut_events_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
events_keys
])
@property
@
cached_
property
def
hits
(
self
):
"""
wrapper for offline hits.
...
...
@@ -407,13 +390,11 @@ class OfflineReader:
Class
OfflineHits.
"""
if
self
.
_hits
is
None
:
self
.
_hits
=
OfflineHits
(
self
.
keys
.
cut_hits_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
hits_keys
])
return
self
.
_hits
return
OfflineHits
(
self
.
keys
.
cut_hits_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
hits_keys
])
@property
@
cached_
property
def
tracks
(
self
):
"""
wrapper for offline tracks.
...
...
@@ -422,14 +403,12 @@ class OfflineReader:
Class
OfflineTracks.
"""
if
self
.
_tracks
is
None
:
self
.
_tracks
=
OfflineTracks
(
self
.
keys
.
cut_tracks_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
tracks_keys
],
fitparameters
=
self
.
keys
.
fitparameters
)
return
self
.
_tracks
@property
return
OfflineTracks
(
self
.
keys
.
cut_tracks_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
tracks_keys
],
fitparameters
=
self
.
keys
.
fitparameters
)
@cached_property
def
mc_hits
(
self
):
"""
wrapper for offline mc hits.
...
...
@@ -438,13 +417,11 @@ class OfflineReader:
Class
OfflineHits.
"""
if
self
.
_mc_hits
is
None
:
self
.
_mc_hits
=
OfflineHits
(
self
.
keys
.
cut_hits_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
mc_hits_keys
])
return
self
.
_mc_hits
return
OfflineHits
(
self
.
keys
.
cut_hits_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
mc_hits_keys
])
@property
@
cached_
property
def
mc_tracks
(
self
):
"""
wrapper for offline mc tracks.
...
...
@@ -453,20 +430,16 @@ class OfflineReader:
Class
OfflineTracks.
"""
if
self
.
_mc_tracks
is
None
:
self
.
_mc_tracks
=
OfflineTracks
(
self
.
keys
.
cut_tracks_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
mc_tracks_keys
],
fitparameters
=
self
.
keys
.
fitparameters
)
return
self
.
_mc_tracks
@property
return
OfflineTracks
(
self
.
keys
.
cut_tracks_keys
,
[
self
.
_data
[
key
]
for
key
in
self
.
keys
.
mc_tracks_keys
],
fitparameters
=
self
.
keys
.
fitparameters
)
@cached_property
def
usr
(
self
):
if
self
.
_usr
is
None
:
self
.
_usr
=
Usr
(
self
.
_file_path
)
return
self
.
_usr
return
Usr
(
self
.
_file_path
)
@property
@
cached_
property
def
best_reco
(
self
):
"""
returns the best reconstructed track fit data. The best fit is defined
as the track fit with the maximum reconstruction stages. When
"
nan
"
is
...
...
@@ -481,24 +454,22 @@ class OfflineReader:
numpy recarray
a recarray of the best track fit data (reconstruction data).
"""
if
self
.
_best_reco
is
None
:
keys
=
"
,
"
.
join
(
self
.
keys
.
fit_keys
[:
-
1
])
empty_fit_info
=
np
.
array
(
[
match
for
match
in
self
.
_find_empty
(
self
.
tracks
.
fitinf
)])
fit_info
=
[
i
for
i
,
j
in
zip
(
self
.
tracks
.
fitinf
,
empty_fit_info
[:,
1
])
if
j
is
not
None
]
stages
=
self
.
_get_max_reco_stages
(
self
.
tracks
.
rec_stages
)
fit_data
=
np
.
array
([
i
[
j
]
for
i
,
j
in
zip
(
fit_info
,
stages
[:,
2
])])
rows_size
=
len
(
max
(
fit_data
,
key
=
len
))
equal_size_data
=
np
.
vstack
([
np
.
hstack
([
i
,
np
.
zeros
(
rows_size
-
len
(
i
))
+
np
.
nan
])
for
i
in
fit_data
])
self
.
_best_reco
=
np
.
core
.
records
.
fromarrays
(
equal_size_data
.
transpose
(),
names
=
keys
)
return
self
.
_best_reco
keys
=
"
,
"
.
join
(
self
.
keys
.
fit_keys
[:
-
1
])
empty_fit_info
=
np
.
array
(
[
match
for
match
in
self
.
_find_empty
(
self
.
tracks
.
fitinf
)])
fit_info
=
[
i
for
i
,
j
in
zip
(
self
.
tracks
.
fitinf
,
empty_fit_info
[:,
1
])
if
j
is
not
None
]
stages
=
self
.
_get_max_reco_stages
(
self
.
tracks
.
rec_stages
)
fit_data
=
np
.
array
([
i
[
j
]
for
i
,
j
in
zip
(
fit_info
,
stages
[:,
2
])])
rows_size
=
len
(
max
(
fit_data
,
key
=
len
))
equal_size_data
=
np
.
vstack
([
np
.
hstack
([
i
,
np
.
zeros
(
rows_size
-
len
(
i
))
+
np
.
nan
])
for
i
in
fit_data
])
return
np
.
core
.
records
.
fromarrays
(
equal_size_data
.
transpose
(),
names
=
keys
)
def
_get_max_reco_stages
(
self
,
reco_stages
):
"""
find the longest reconstructed track based on the maximum size of
...
...
This diff is collapsed.
Click to expand it.
Preview
0%
Loading
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Save comment
Cancel
Please
register
or
sign in
to comment