cetin 9 months ago
parent
commit
60e2490880
  1. 39
      moore_options/get_resolution_and_eff_data.py
  2. 10
      moore_options/get_tracking_losses.py

39
moore_options/get_resolution_and_eff_data.py

@ -36,11 +36,11 @@ from RecoConf.hlt1_tracking import (
get_global_materiallocator,
)
decay = "test"
decay = "B"
options.evt_max = -1
options.ntuple_file = f"data/resolutions_and_effs_{decay}_endVelo_momEff.root"
options.ntuple_file = f"data/resolutions_and_effs_{decay}_endVelo_idealstateP.root"
options.input_type = "ROOT"
@ -56,6 +56,7 @@ elif decay == "test":
options.dddb_tag = "dddb-20210617"
options.conddb_tag = "sim-20210617-vc-md100"
options.simulation = True
options.output_level = 3
def run_tracking_resolution():
@ -91,19 +92,19 @@ def run_tracking_resolution():
LinksToLHCbIDs=links_to_lhcbids,
)
res_checker_forward = check_track_resolution(tracks["Forward"], suffix="Forward")
res_checker_best_long = check_track_resolution(
tracks["BestLong"],
suffix="BestLong",
)
res_checker_best_forward = check_track_resolution(
dict(v1=fitted_forward_tracks),
suffix="BestForward",
)
res_checker_seed = check_track_resolution(
tracks["Seed"],
suffix="Seed",
)
# res_checker_forward = check_track_resolution(tracks["Forward"], suffix="Forward")
# res_checker_best_long = check_track_resolution(
# tracks["BestLong"],
# suffix="BestLong",
# )
# res_checker_best_forward = check_track_resolution(
# dict(v1=fitted_forward_tracks),
# suffix="BestForward",
# )
# res_checker_seed = check_track_resolution(
# tracks["Seed"],
# suffix="Seed",
# )
eff_checker_forward = check_tracking_efficiency(
"Forward",
@ -139,10 +140,10 @@ def run_tracking_resolution():
)
data = [
res_checker_forward,
res_checker_best_long,
res_checker_best_forward,
res_checker_seed,
# res_checker_forward,
# res_checker_best_long,
# res_checker_best_forward,
# res_checker_seed,
eff_checker_forward,
eff_checker_match,
eff_checker_best_long,

10
moore_options/get_tracking_losses.py

@ -40,6 +40,9 @@ tested by mc_matching_example.py
decay = "B"
options.evt_max = -1
if decay == "B":
options.input_files = glob.glob("/auto/data/guenther/Bd_Kstee/*.xdigi")
elif decay == "BJpsi":
@ -55,12 +58,7 @@ options.dddb_tag = "dddb-20210617"
options.simulation = True
options.input_type = "ROOT"
options.ntuple_file = (
f"data/tracking_losses_ntuple_{decay}_default_radlength_beginVelo.root"
)
options.evt_max = -1
options.ntuple_file = f"data/tracking_losses_ntuple_{decay}_endVelo_idealstateP.root"
def run_tracking_losses():

Loading…
Cancel
Save