Truncate humongous eloss files, clean up Makevertex a little

This commit is contained in:
Sudarsan Balakrishnan 2026-07-07 15:28:32 -04:00
parent 3b99dc05d1
commit 6034c5ceeb
25 changed files with 2 additions and 6395715 deletions

View File

@ -9,44 +9,6 @@ Int_t colors[40] = {
kRed-7, kGreen-7, kBlue-7, kYellow-7, kMagenta-7, kCyan-7, kOrange-7, 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 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 "MakeVertex.h"
#include "Armory/ClassPW.h" #include "Armory/ClassPW.h"
#include "Armory/HistPlotter.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; return 0.000165896*z*z + 4.61626e-08*z + 32.067;
} }
TF1 f_triangle("f1",triangle,-146,146);
// Global instances // Global instances
PW pwinstance; 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()); //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) { 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()); 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_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"); 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()); //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) { 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("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("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"); plotter->Fill2D("pcdEAfix_vs_qqqzguess",300,-200,200,800,0,20000,pcz_guess_int,dEfix,"pcE_vs_qqqguess");

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff