From 80a8746b2f752bcfa1489066cddbc46ca8224ea3 Mon Sep 17 00:00:00 2001
From: Dmitry Kalinkin <dmitry.kalinkin@gmail.com>
Date: Sat, 20 Jul 2024 17:50:34 -0400
Subject: [PATCH] tracking_performances: replace 15 GeV with 20 GeV (supported
by campaigns)
---
benchmarks/tracking_performances/Snakefile | 26 +++++++++++++------
.../Tracking_Performances.C | 6 ++---
.../doCompare_truth_real_widebins_mom.C | 4 +--
3 files changed, 23 insertions(+), 13 deletions(-)
diff --git a/benchmarks/tracking_performances/Snakefile b/benchmarks/tracking_performances/Snakefile
index 5cbf9278..8bbdee90 100644
--- a/benchmarks/tracking_performances/Snakefile
+++ b/benchmarks/tracking_performances/Snakefile
@@ -48,12 +48,22 @@ rule tracking_performance_at_momentum:
input:
script="benchmarks/tracking_performances/Tracking_Performances.C",
# TODO pass as a file list?
- sim=lambda wildcards: expand(
- "sim_output/tracking_performance/{DETECTOR_CONFIG}/{{PARTICLE}}/{ENERGY}/{PHASE_SPACE}/{{PARTICLE}}_{ENERGY}_{PHASE_SPACE}.{INDEX:04d}.eicrecon.tree.edm4eic.root",
- DETECTOR_CONFIG="epic_craterlake_tracking_only",
- ENERGY=f"{float(wildcards.MOMENTUM):.0f}GeV" if float(wildcards.MOMENTUM) >= 1 else f"{float(wildcards.MOMENTUM) * 1000:.0f}MeV",
- PHASE_SPACE=["3to50deg", "45to135deg", "130to177deg"],
- INDEX=range(1),
+ sim=lambda wildcards: branch(
+ wildcards.CAMPAIGN == "local",
+ then=expand(
+ "sim_output/tracking_performance/{DETECTOR_CONFIG}/{{PARTICLE}}/{ENERGY}/{PHASE_SPACE}/{{PARTICLE}}_{ENERGY}_{PHASE_SPACE}.{INDEX:04d}.eicrecon.tree.edm4eic.root",
+ DETECTOR_CONFIG="epic_craterlake_tracking_only",
+ ENERGY=f"{float(wildcards.MOMENTUM):.0f}GeV" if float(wildcards.MOMENTUM) >= 1 else f"{float(wildcards.MOMENTUM) * 1000:.0f}MeV",
+ PHASE_SPACE=["3to50deg", "45to135deg", "130to177deg"],
+ INDEX=range(1),
+ ),
+ otherwise=expand(
+ "EPIC/RECO/{{CAMPAIGN}}/epic_craterlake/SINGLE/{{PARTICLE}}/{ENERGY}/{PHASE_SPACE}/{{PARTICLE}}_{ENERGY}_{PHASE_SPACE}.{INDEX:04d}.eicrecon.tree.edm4eic.root",
+ DETECTOR_CONFIG="epic_craterlake",
+ ENERGY=f"{float(wildcards.MOMENTUM):.0f}GeV" if float(wildcards.MOMENTUM) >= 1 else f"{float(wildcards.MOMENTUM) * 1000:.0f}MeV",
+ PHASE_SPACE=["3to50deg", "45to135deg", "130to177deg"],
+ INDEX=range(1),
+ ),
),
output:
"{SEEDING}/pi-/mom/Performances_mom_{MOMENTUM}_mom_resol_{SEEDING_IGNORE}_{PARTICLE}.root",
@@ -80,14 +90,14 @@ rule tracking_performance_summary_at_eta:
"realseed/pi-/mom/Performances_mom_{MOMENTUM:.1f}_mom_resol_realseed_pi-.root",
"realseed/pi-/dca/Performances_dca_{MOMENTUM:.1f}_dca_resol_realseed_pi-.root",
],
- MOMENTUM=[0.5, 1.0, 2.0, 5.0, 10.0, 15.0],
+ MOMENTUM=[0.5, 1.0, 2.0, 5.0, 10.0, 20.0],
),
script="benchmarks/tracking_performances/doCompare_truth_real_widebins_mom.C",
output:
expand(
"Debug_Plots/{SEEDING}/pi-/mom/{SEEDING}_mom_resol_mom{MOMENTUM:.1f}_{{ETA_MIN}}_eta_{{ETA_MAX}}.png",
SEEDING=["real", "truth"],
- MOMENTUM=[0.5, 1.0, 2.0, 5.0, 10.0, 15.0],
+ MOMENTUM=[0.5, 1.0, 2.0, 5.0, 10.0, 20.0],
),
"Final_Results/pi-/mom/mom_resol_{ETA_MIN}_eta_{ETA_MAX}.png",
"Final_Results/pi-/mom/mom_resol_{ETA_MIN}_eta_{ETA_MAX}.root",
diff --git a/benchmarks/tracking_performances/Tracking_Performances.C b/benchmarks/tracking_performances/Tracking_Performances.C
index 6d1b213c..bfb751d6 100644
--- a/benchmarks/tracking_performances/Tracking_Performances.C
+++ b/benchmarks/tracking_performances/Tracking_Performances.C
@@ -34,11 +34,11 @@ void Tracking_Performances(TString filename="tracking_output",TString particle="
int theta_val[nbins_eta+1] ={3,50,45,135,130,177};
int nfiles = 100;
double eta[nbins_eta+1]={-3.5,-2.5,-1.0,1.0,2.5,3.5};
- double pt[nbins_eta+1]={0.1,0.5,1.0,2.0,5.0,10.0};
+ double pt[nbins_eta+1]={0.5,1.0,2.0,5.0,10.0,20.1};
TH1D *histp[nbins_eta];
- TH3D *h_d0xy_3d= new TH3D("h_d0xy_3d","Transverse Pointing Resolution",500,-0.1,0.1,70,-3.5,3.5,200,0.,20.);
- TH3D *h_d0z_3d= new TH3D("h_d0z_3d","Longitudinal Pointing Resolution",500,-0.1,0.1,70,-3.5,3.5,200,0.,20.);
+ TH3D *h_d0xy_3d= new TH3D("h_d0xy_3d","Transverse Pointing Resolution",500,-0.1,0.1,70,-3.5,3.5,201,0.,20.1);
+ TH3D *h_d0z_3d= new TH3D("h_d0z_3d","Longitudinal Pointing Resolution",500,-0.1,0.1,70,-3.5,3.5,201,0.,20.1);
for (int i=0; i<nbins_eta; i++){
histp[i] = new TH1D(Form("hist_etabin%d",i),Form("hist_etabin%d",i),600,-0.3,0.3);
diff --git a/benchmarks/tracking_performances/doCompare_truth_real_widebins_mom.C b/benchmarks/tracking_performances/doCompare_truth_real_widebins_mom.C
index 4d417428..283dffce 100644
--- a/benchmarks/tracking_performances/doCompare_truth_real_widebins_mom.C
+++ b/benchmarks/tracking_performances/doCompare_truth_real_widebins_mom.C
@@ -24,7 +24,7 @@ void doCompare_truth_real_widebins_mom(TString particle = "pi-",double etamin=-1
gStyle->SetOptStat(1);
const Int_t nfiles = 6;
- double mom[nfiles] ={0.5,1.0,2.0,5.0,10.0,15.0};
+ double mom[nfiles] ={0.5,1.0,2.0,5.0,10.0,20.0};
std::vector<double> momV_truth, momV_real, momresolV_truth, err_momresolV_truth, momresolV_real, err_momresolV_real;
momV_truth.clear(); momV_real.clear(); momresolV_truth.clear(); err_momresolV_truth.clear(); momresolV_real.clear(); err_momresolV_real.clear();
@@ -150,7 +150,7 @@ void doCompare_truth_real_widebins_mom(TString particle = "pi-",double etamin=-1
mgMom->Add(gr1);
mgMom->Add(gr2);
c_mom->cd();
- mgMom->GetXaxis()->SetRangeUser(0.40,15.2);
+ mgMom->GetXaxis()->SetRangeUser(0.40,20.2);
mgMom->GetYaxis()->SetRangeUser(0.,10.0);
mgMom->Draw("AP");
lmom->AddEntry(gr1,"Truth Seeding");
--
GitLab