rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Truncate humongous eloss files, clean up Makevertex a little

Browse Source
This commit is contained in:
Sudarsan Balakrishnan 2026-07-07 15:28:32 -04:00
parent 1a16ef41f3
commit 3fdf5e9c41
1 changed files with 2 additions and 41 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
MakeVertex.C
View File

@ -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");