35mm Z-vtz; tracking resolutions incorporated

eic · May 3, 2024 · 60495f1 · 60495f1
1 parent f05347d
commit 60495f1
Show file tree

Hide file tree

Showing 6 changed files with 32 additions and 4 deletions.
diff --git a/database/dphi_seed_param.dat b/database/dphi_seed_param.dat
@@ -0,0 +1,6 @@
+-1.39 -1.47 11.6033 -0.2199
+-1.55 -1.64 11.0970 -0.1170
+-1.96 -2.10 8.3530 -0.1203
+-2.25 -2.44 10.8800 -0.1672
+-2.66 -2.95 15.2172 -0.5944
+-3.13 -3.64 23.8946 -1.4373
diff --git a/database/dphi_seed_param.reformatted.dat b/database/dphi_seed_param.reformatted.dat
@@ -0,0 +1,6 @@
+-1.39 -1.50 11.6033 -0.2199
+-1.50 -1.80 11.0970 -0.1170
+-1.80 -2.17  8.3530 -0.1203
+-2.17 -2.55 10.8800 -0.1672
+-2.55 -3.05 15.2172 -0.5944
+-3.05 -3.64 23.8946 -1.4373
diff --git a/database/dtheta_seed_param.dat b/database/dtheta_seed_param.dat
@@ -0,0 +1,6 @@
+-1.39 -1.47 5.4002 -0.0000
+-1.55 -1.64 4.0950 -0.1458
+-1.96 -2.10 2.2099 -0.1344
+-2.25 -2.44 2.1430 -0.1443
+-2.66 -2.95 1.9840 -0.1409
+-3.13 -3.64 1.8964 -0.1431
diff --git a/database/dtheta_seed_param.reformatted.dat b/database/dtheta_seed_param.reformatted.dat
@@ -0,0 +1,6 @@
+-1.39 -1.50 5.4002 -0.0000
+-1.50 -1.80 4.0950 -0.1458
+-1.80 -2.17 2.2099 -0.1344
+-2.17 -2.55 2.1430 -0.1443
+-2.55 -3.05 1.9840 -0.1409
+-3.05 -3.64 1.8964 -0.1431
diff --git a/epic/include/epic.default.h b/epic/include/epic.default.h
@@ -14,7 +14,7 @@
 //
 // Optional smearing of primary vertices along the beam line; applies to HEPMC3 input as well;
 // 10cm (proton bunch size) is a reasonable value;
-#define _PRIMARY_VERTEX_SIGMA_               (10*cm)
+#define _PRIMARY_VERTEX_SIGMA_               (35*mm)//(10*cm)
 
 //
 // All this can be overriden by '-i' cmd line switch (HEPMC3 input); also see '-s <statistics>' key;

diff --git a/scripts/reco-epic.C b/scripts/reco-epic.C
@@ -35,7 +35,9 @@ void reco_epic(const char *dfname, const char *cfname = 0)
   // Carelessly remove (0x1 << n)x(0x1 << n) square area "around" these hits;
   reco->SetBlackoutBlowupValue(3);
 
-  auto hmatch = new TH1D("hmatch", "PID evaluation correctness",       2,    0,      2);
+  reco->ImportTrackingSmearing("./database/dtheta_seed_param.reformatted.dat", "./database/dphi_seed_param.reformatted.dat");
+
+  auto hmatch = new TH1D("hmatch", "PID evaluation correctness",       3,    0,      3);
   //auto hthtr1 = new TH1D("thtr1",  "Cherenkov angle (track)",        200,  220,    320);
   auto hthtr1 = new TH1D("thtr1",  "Cherenkov angle (track)",        40,  270, 290);
   // For a dual aerogel configuration;
@@ -52,10 +54,12 @@ void reco_epic(const char *dfname, const char *cfname = 0)
       for(auto mcparticle: event->ChargedParticles()) {
 	if (!mcparticle->IsPrimary()) continue;
 
-	if (mcparticle->GetPDG() == mcparticle->GetRecoPdgCode()) {
+	if (mcparticle->GetRecoCherenkovHitCount() <= 3) {
 	  hmatch->Fill(0.5);
-	} else {
+	} else if (mcparticle->GetPDG() == mcparticle->GetRecoPdgCode()) {
 	  hmatch->Fill(1.5);
+	} else {
+	  hmatch->Fill(2.5);
 	} //if	  
 
 	hthtr1->Fill(1000*mcparticle->GetRecoCherenkovAverageTheta(a1));