From 3fdf5e9c41b2f2c88ba8e6753c567192adf0d747 Mon Sep 17 00:00:00 2001 From: Sudarsan B Date: Tue, 7 Jul 2026 15:28:32 -0400 Subject: [PATCH] Truncate humongous eloss files, clean up Makevertex a little --- MakeVertex.C | 43 ++----------------------------------------- 1 file changed, 2 insertions(+), 41 deletions(-) diff --git a/MakeVertex.C b/MakeVertex.C index e2a20b0..9678c87 100644 --- a/MakeVertex.C +++ b/MakeVertex.C @@ -9,44 +9,6 @@ Int_t colors[40] = { kRed-7, kGreen-7, kBlue-7, kYellow-7, kMagenta-7, kCyan-7, kOrange-7, kSpring-7, kTeal-7, kAzure-7, kViolet-7, kPink-7, kGray+2 }; -double triangle(double* y, double* p) { - double result = 0; - double offset = 0; - double lowlim=6500; - double hilim = 9500; - double slope = 160.0;//(hilim-lowlim)/20.0; - double z_grid[8] = {-147.998,-101.946,-59.7634,-19.6965,19.6965,59.7634,101.946,147.998}; //crossover z values - for(int i=0;i<7;i++) { - double zavg = (z_grid[i] + z_grid[i+1])*0.5; //midpoint about which we pivot - if(y[0] <= zavg+offset && y[0] > z_grid[i]+offset) { - result = slope*(y[0]-z_grid[i])+lowlim; - } else if(y[0]>zavg+offset && y[0] <= z_grid[i+1]+offset){ - result = -slope*(y[0]-z_grid[i+1])+lowlim; - } - } - return result; -} - -/*double triangle(double* y, double* p) { - double result = 0; - double offset = 0; - double lowlim=7500; - double hilim = 9500; - double slope = 160.0;//(hilim-lowlim)/20.0; - //double z_grid[8] = {-147.998,-101.946,-59.7634,-19.6965,19.6965,59.7634,101.946,147.998}; //crossover z values - double z_grid[8] = {-147.998,-101.946,-59.7634,-19.6965,15.6965,55.7634,101.946,147.998}; //crossover z values - for(int i=0;i<7;i++) { - double zavg = (z_grid[i] + z_grid[i+1])*0.5; //midpoint about which we pivot - if(y[0] <= zavg+offset && y[0] > z_grid[i]+offset) { - result = slope*(y[0]-z_grid[i])+lowlim; - } else if(y[0]>zavg+offset && y[0] <= z_grid[i+1]+offset){ - result = -slope*(y[0]-z_grid[i+1])+lowlim; - } - } - if(y[0]>69) result*=(4./5.); - return result; -}*/ - #include "MakeVertex.h" #include "Armory/ClassPW.h" #include "Armory/HistPlotter.h" @@ -114,7 +76,6 @@ double z_to_crossover_rho(double z) { return 0.000165896*z*z + 4.61626e-08*z + 32.067; } -TF1 f_triangle("f1",triangle,-146,146); // Global instances PW pwinstance; @@ -1238,7 +1199,7 @@ Bool_t MakeVertex::Process(Long64_t entry) //plotter->Fill2D("pcz_vs_sx3pczguess_A2C1_strip"+std::to_string(sx3event.ch2),300,-200,200,600,-200,200,pczguess,pcevent.pos.Z()); } if(pcevent.multi1==1 && pcevent.multi2==2) { - double dEfix = pcevent.Energy1*8500./f_triangle.Eval(pczguess_int); + double dEfix = pcevent.Energy1;//8500./f_triangle.Eval(pczguess_int); double pcz_fix = pcfix_func.Eval(pcevent.pos.Z()); plotter->Fill2D("pcdEA_vs_pczfix_s",300,-200,200,800,0,20000,pcz_fix,pcevent.Energy1,"pcE_vs_sx3guess"); plotter->Fill2D("pcdEA_vs_sx3zguess",300,-200,200,800,0,20000,pczguess_int,pcevent.Energy1,"pcE_vs_sx3guess"); @@ -1460,7 +1421,7 @@ Bool_t MakeVertex::Process(Long64_t entry) //plotter->Fill2D("pcz_vs_sx3pczguess_A2C1_strip"+std::to_string(sx3event.ch2),300,-200,200,600,-200,200,pczguess,pcevent.pos.Z()); } if(pcevent.multi1==1) { - double dEfix = pcevent.Energy1*8500./f_triangle.Eval(pcz_guess_int); + double dEfix = pcevent.Energy1;//8500./f_triangle.Eval(pcz_guess_int); plotter->Fill2D("pcdEA_vs_qqqzguess",300,-200,200,800,0,20000,pcz_guess_int,pcevent.Energy1,"pcE_vs_qqqguess"); plotter->Fill2D("pcdE2A_vs_qqqzguess",300,-200,200,800,0,20000,pcz_guess_int,pcevent.Energy1*sinTheta2,"pcE_vs_qqqguess"); plotter->Fill2D("pcdEAfix_vs_qqqzguess",300,-200,200,800,0,20000,pcz_guess_int,dEfix,"pcE_vs_qqqguess");