diff --git a/pccal/anode_gainmatch.C b/pccal/anode_gainmatch.C deleted file mode 100755 index fa267be..0000000 --- a/pccal/anode_gainmatch.C +++ /dev/null @@ -1,67 +0,0 @@ -void anode_gainmatch(){ - TFile *f = new TFile("../results_run16.root"); - - TH2F *pc_index_h2d = (TH2F*)(f->Get("hRawPC/PC_Index_Vs_Energy")); - std::cout << pc_index_h2d << std::endl; - TCanvas c("c1","c1",0,0,1600,800); - //TCanvas c_g("cg","cg",0,900,400,400); - c.Divide(2,1); - auto c1=c.cd(1); - pc_index_h2d->Draw("COLZ"); - pc_index_h2d->GetYaxis()->SetRangeUser(240,5000); - auto c2=c.cd(2); - c2->SetLogy(); - TH1F *h_1d=NULL; - int bin_index=1; - std::vector> all_peaks; - std::vector found_wire_list; - while(bin_index<=24) { - h_1d=(TH1F*)(pc_index_h2d->ProjectionY("_py",bin_index,bin_index)); - auto c1 = c.cd(1); - TBox box(pc_index_h2d->GetXaxis()->GetBinLowEdge(bin_index),0,pc_index_h2d->GetXaxis()->GetBinUpEdge(bin_index),pc_index_h2d->GetYaxis()->GetXmax()); - box.SetFillColorAlpha(kYellow+3,0.3); - box.Draw("SAME"); - c1->Modified(); c1->Update(); - //while(c1->WaitPrimitive()); - - TSpectrum s; - auto c2 = c.cd(2); - h_1d->Draw(); - c2->Modified(); c2->Update(); - int npeaks = s.Search(h_1d,8,"",0.02); std::cout << npeaks << std::endl; - if(npeaks>=3) { - std::vector xpeaks(s.GetPositionX(),s.GetPositionX()+npeaks); - std::sort(xpeaks.begin(),xpeaks.end(),std::greater()); - found_wire_list.push_back((int)pc_index_h2d->GetXaxis()->GetBinCenter(bin_index)); - all_peaks.push_back(xpeaks); - } - while(c2->WaitPrimitive()); - bin_index++; - } - c.cd(2)->SetLogy(kFALSE); - gStyle->SetOptFit(1111); - - std::ofstream outfile("anode_gm_coeffs.dat"); - outfile << found_wire_list.at(0) << " " - << 1.0 << " " - << 0.0 << std::endl; - - for(int i=0; iGetParameter(1) << " " - << ((TF1*)g.FindObject("pol1"))->GetParameter(0) << std::endl; - c2->Modified(); - c2->Update(); - while(c2->WaitPrimitive()); - } - outfile.close(); - f->Close(); - return; -} diff --git a/pccal/anode_gm_coeffs.dat b/pccal/anode_gm_coeffs.dat deleted file mode 100755 index a241cfe..0000000 --- a/pccal/anode_gm_coeffs.dat +++ /dev/null @@ -1,24 +0,0 @@ -0 1 0 -0 0.937314 -16.871 -2 0.965461 -1.54376 -3 0.926501 -3.27662 -4 0.905634 2.54577 -5 0.905634 -11.0387 -6 0.853919 6.23079 -7 0.945588 -9.54044 -8 0.884454 -11.8262 -9 0.922501 -3.42538 -10 0.903053 9.28069 -11 0.914653 9.87642 -12 0.965332 13.2526 -13 0.923847 -3.41775 -14 0.93845 25.9901 -15 0.955424 12.324 -16 0.95116 4.99595 -17 0.910745 2.86648 -18 0.941376 4.57217 -19 0.871622 932.111 -20 1.00624 7.86358 -21 0.969834 -45.001 -22 0.89304 -31.5635 -23 0.933226 4.02193 diff --git a/pccal/cathode_gainmatch.C b/pccal/cathode_gainmatch.C deleted file mode 100755 index d52e32a..0000000 --- a/pccal/cathode_gainmatch.C +++ /dev/null @@ -1,69 +0,0 @@ -void cathode_gainmatch(){ - TFile *f = new TFile("../results_run17.root"); - TH2F *pc_index_h2d = (TH2F*)(f->Get("hRawPC/PC_Index_Vs_Energy")); - std::cout << pc_index_h2d << std::endl; - TCanvas c("c1","c1",0,0,1600,800); - //TCanvas c_g("cg","cg",0,900,400,400); - c.Divide(2,1); - auto c1=c.cd(1); - pc_index_h2d->Draw("COLZ"); - pc_index_h2d->GetYaxis()->SetRangeUser(600,pc_index_h2d->GetYaxis()->GetXmax()); - auto c2=c.cd(2); - c2->SetLogy(); - TH1F *h_1d=NULL; - int bin_index=25; - std::vector pulser_heights = {0.01,0.05,0.1,0.15,0.2,0.25,0.3,0.5}; - std::vector> all_peaks; - std::vector found_wire_list; - while(bin_index<=48) { - h_1d=(TH1F*)(pc_index_h2d->ProjectionY("_py",bin_index,bin_index)); - auto c1 = c.cd(1); - TBox box(pc_index_h2d->GetXaxis()->GetBinLowEdge(bin_index),0,pc_index_h2d->GetXaxis()->GetBinUpEdge(bin_index),pc_index_h2d->GetYaxis()->GetXmax()); - box.SetFillColorAlpha(kYellow+3,0.3); - box.Draw("SAME"); - c1->Modified(); c1->Update(); - //while(c1->WaitPrimitive()); - - TSpectrum s; - auto c2 = c.cd(2); - h_1d->Draw(); - c2->Modified(); c2->Update(); - int npeaks = s.Search(h_1d,20,"",0.1); std::cout << npeaks << std::endl; - if(npeaks==8) { - std::vector xpeaks(s.GetPositionX(),s.GetPositionX()+npeaks); - for(int i=0; i<8; i++) { - std::cout << pc_index_h2d->GetXaxis()->GetBinCenter(bin_index) << " " << xpeaks.at(i) << " " << xpeaks.at(i)/pulser_heights.at(i) << std::endl; - } - std::sort(xpeaks.begin(),xpeaks.end(),std::greater()); - found_wire_list.push_back((int)pc_index_h2d->GetXaxis()->GetBinCenter(bin_index)); - all_peaks.push_back(xpeaks); - } - while(c2->WaitPrimitive()); - bin_index++; - } - c.cd(2)->SetLogy(kFALSE); - gStyle->SetOptFit(1111); - - std::ofstream outfile("cathode_gm_coeffs.dat"); - outfile << found_wire_list.at(0) << " " - << 1.0 << " " - << 0.0 << std::endl; - - for(int i=1; iGetParameter(1) << " " - << ((TF1*)g.FindObject("pol1"))->GetParameter(0) << std::endl; - c2->Modified(); - c2->Update(); - while(c2->WaitPrimitive()); - } - outfile.close(); - f->Close(); - return; -} diff --git a/pccal/cathode_gm_coeffs.dat b/pccal/cathode_gm_coeffs.dat deleted file mode 100755 index b1aaa77..0000000 --- a/pccal/cathode_gm_coeffs.dat +++ /dev/null @@ -1,21 +0,0 @@ -24 1 0 -25 0.941896 6.16135 -26 0.980284 2.86886 -27 0.983166 -3.82952 -28 0.978704 -2.89713 -29 0.964947 2.25786 -30 0.94514 0.925074 -31 0.977231 1.6493 -32 0.919527 5.82742 -33 0.972243 2.88061 -34 0.928892 7.61384 -35 0.947376 -0.644223 -36 0.875342 6.066 -38 0.970953 6.262 -40 0.918408 -3.27891 -41 0.913619 4.11288 -42 0.954083 2.21261 -43 0.993037 5.48924 -45 0.926406 -19.719 -46 1.00459 5.14574 -47 0.942483 5.54183 diff --git a/pccal/slope_intercept_26Al.dat b/pccal/slope_intercept_26Al.dat deleted file mode 100755 index 875aaec..0000000 --- a/pccal/slope_intercept_26Al.dat +++ /dev/null @@ -1,50 +0,0 @@ -#Histogram Number Slope Intercept -#Histogram Number Slope Intercept -0 0.931015 -1.35431 -1 1 -1.87356e-10 -2 0.964185 1.49989 -3 0.92638 -1.30621 -4 0.905569 1.00834 -5 0.901182 0.470903 -6 0.853932 3.32687 -7 0.942785 1.08887 -8 0.878904 -0.0107433 -9 0.922662 -2.32259 -10 0.903343 8.38332 -11 0.914227 6.56108 -12 0.961008 23.0982 -13 0.920976 5.22104 -14 0.936584 31.5073 -15 0.959044 5.43267 -16 0.95263 -0.404053 -17 0.90953 4.82833 -18 0.940277 10.3629 -19 0.86746 -17.8678 -20 1.00683 4.76371 -21 0.968342 -43.9496 -22 0.892882 -32.0742 -23 0.933615 1.10704 -24 1 -2.89219e-10 -25 0.942098 -0.105169 -26 0.980862 -0.732032 -27 0.982975 -2.22704 -28 0.978815 -1.51477 -29 0.965245 -2.19515 -30 0.945384 -0.892599 -31 0.977408 -0.908592 -32 0.919546 3.25464 -33 0.972194 2.44956 -34 0.92852 5.44745 -35 0.947098 1.40531 -36 0.875491 -1.13145 -37 1 0 -38 0.970862 2.86019 -39 1 0 -40 0.91793 -3.80615 -41 0.913897 -2.12964 -42 0.954014 -0.760604 -43 0.993616 -1.40278 -44 1 0 -45 0.926169 -21.2016 -46 1.00577 -2.14281 -47 0.943312 -1.26464 diff --git a/scratch/sx3z_vs_phiz/scan_offset.C b/scratch/sx3z_vs_phiz/scan_offset.C deleted file mode 100755 index ebf6e22..0000000 --- a/scratch/sx3z_vs_phiz/scan_offset.C +++ /dev/null @@ -1,81 +0,0 @@ -#include "testmodel.h" - -int quit=0; -void handler(int){quit=1;} - -int colors[] = {kSpring+3, kRed, kGreen+3, kBlue+3, kViolet, kOrange, kSpring-7, kAzure-5}; -void scan_offset(){ - signal(SIGINT,handler); - TCanvas c("c1","c1",0,0,1600,800); - c.Divide(2,1); - - TF1 f1("model",model,-200,200,2); - f1.SetNpx(10000); - std::vector pars = {0.0,1.}; - f1.SetParameters(pars.data()); - f1.SetLineColor(kGreen+2); - f1.SetLineStyle(kLine); - - - - - TFile* f=NULL; - std::vector files; - int ctr=0; - for(int i=12; i<=21; i++) { - auto c1=c.cd(1); - c1->SetGrid(1,1); - f = new TFile(Form("../../results_run%d.root",i)); - - if(i==12) { - //TH2F *h2 = (TH2F*)(f->Get("phicut/pczguess_vs_pc_int")); - TH2F *h23 = (TH2F*)(f->Get("phicut/pczguess_vs_pc_int_A1C2")); - h23->SetLineColorAlpha(kOrange,0.75); - h23->Draw("box SAME"); - - } else { - //TH2F *h2 = (TH2F*)(f->Get("phicut/pczguess_vs_pc_int")); - //TH2F *h2 = (TH2F*)(f->Get("pcz_vs_sx3pczguess_A1C2_strip12")); - TH2F *h2 = (TH2F*)(f->Get("pcz_vs_sx3pczguess_A1C2")); - //TH2F *h2 = (TH2F*)(f->Get("hPCQQQ/PC_XY_Projection_QQQ2")); - if(!h2) continue; - h2->SetTitle(Form("case%d",i)); - //h2->Draw("colz same"); - h2->SetLineColorAlpha(colors[ctr],0.75); - h2->Draw("box same"); - f1.Draw("same"); - } - TF1 eqline("x","x",-200,200); - eqline.Draw("SAME"); - c1->Modified(); - c1->Update(); - ctr+=1; - - - auto c2=c.cd(2); - c2->SetGrid(1,1); - - TH2F *h3 = (TH2F*)(f->Get("sx3phi_vs_pcphi1")); -// TH2F *h2 = (TH2F*)(f->Get("hPCQQQ/PC_XY_Projection_QQQ2")); - if(!h3) continue; - h3->SetTitle(Form("case%d",i)); - h3->Draw("colz"); - eqline.Draw("SAME"); - c2->Modified(); - c2->Update(); - - while(gPad->WaitPrimitive()); - - files.emplace_back(f); - if(i==21) { - i=11; - c.Clear(); - c.Divide(2,1); - ctr=0; - } - if(quit) break; - } - for(auto file : files) { - file->Close(); - } -} diff --git a/scratch/sx3z_vs_phiz/scan_offset_fix.C b/scratch/sx3z_vs_phiz/scan_offset_fix.C deleted file mode 100755 index e94faad..0000000 --- a/scratch/sx3z_vs_phiz/scan_offset_fix.C +++ /dev/null @@ -1,81 +0,0 @@ -#include "testmodel.h" - -int quit=0; -void handler(int){quit=1;} - -int colors[] = {kSpring+3, kRed, kGreen+3, kBlue+3, kViolet, kOrange, kSpring-7, kAzure-5}; -void scan_offset_fix(){ - signal(SIGINT,handler); - TCanvas c("c1","c1",0,0,1600,800); - c.Divide(2,1); - - TF1 f1("model",model,-200,200,2); - f1.SetNpx(10000); - std::vector pars = {0.0,1.}; - f1.SetParameters(pars.data()); - f1.SetLineColor(kGreen+2); - f1.SetLineStyle(kLine); - - - - - TFile* f=NULL; - std::vector files; - int ctr=0; - for(int i=12; i<=21; i++) { - auto c1=c.cd(1); - c1->SetGrid(1,1); - f = new TFile(Form("../../results_run%d.root",i)); - - if(i==12) { - //TH2F *h2 = (TH2F*)(f->Get("phicut/pczguess_vs_pc_int")); - TH2F *h23 = (TH2F*)(f->Get("pczfix_vs_qqqpczguess_A1C2")); - h23->SetLineColorAlpha(kOrange,0.75); - h23->Draw("SAME"); - - } else { - //TH2F *h2 = (TH2F*)(f->Get("phicut/pczguess_vs_pc_int")); - //TH2F *h2 = (TH2F*)(f->Get("pcz_vs_sx3pczguess_A1C2_strip12")); - TH2F *h2 = (TH2F*)(f->Get("pczfix_vs_sx3pczguess_A1C2")); - //TH2F *h2 = (TH2F*)(f->Get("hPCQQQ/PC_XY_Projection_QQQ2")); - if(!h2) continue; - h2->SetTitle(Form("case%d",i)); - //h2->Draw("colz same"); - h2->SetLineColorAlpha(colors[ctr],0.75); - h2->Draw("box same"); - //f1.Draw("same"); - } - TF1 eqline("x","x",-200,200); - eqline.Draw("SAME"); - c1->Modified(); - c1->Update(); - ctr+=1; - - - auto c2=c.cd(2); - c2->SetGrid(1,1); - - TH2F *h3 = (TH2F*)(f->Get("sx3phi_vs_pcphi1")); -// TH2F *h2 = (TH2F*)(f->Get("hPCQQQ/PC_XY_Projection_QQQ2")); - if(!h3) continue; - h3->SetTitle(Form("case%d",i)); - h3->Draw("colz"); - eqline.Draw("SAME"); - c2->Modified(); - c2->Update(); - - while(gPad->WaitPrimitive()); - - files.emplace_back(f); - if(i==21) { - i=11; - c.Clear(); - c.Divide(2,1); - ctr=0; - } - if(quit) break; - } - for(auto file : files) { - file->Close(); - } -} diff --git a/scratch/sx3z_vs_phiz/stitch.C b/scratch/sx3z_vs_phiz/stitch.C deleted file mode 100755 index 4fbaa02..0000000 --- a/scratch/sx3z_vs_phiz/stitch.C +++ /dev/null @@ -1,5 +0,0 @@ -{ - - - -} diff --git a/scratch/sx3z_vs_phiz/testmodel.h b/scratch/sx3z_vs_phiz/testmodel.h deleted file mode 100755 index 458d808..0000000 --- a/scratch/sx3z_vs_phiz/testmodel.h +++ /dev/null @@ -1,43 +0,0 @@ -#include -double model(double* x, double* p) { - double result = x[0]; - double factor = 29.0; - double slope = 0.7; - if(TMath::Abs(x[0]) < 16.2) result=x[0]*slope; - else if(TMath::Abs(x[0]) < 49.8 ) result=x[0]*slope+TMath::Sign(1.0,x[0])*factor; - else if(TMath::Abs(x[0]) < 85.2 ) result=x[0]*slope+TMath::Sign(1.0,x[0])*factor*2; - else result=x[0]*slope+TMath::Sign(1.0,x[0])*factor*2; - return result; -} - -double model_invert(double *y, double *q) { - double result=y[0]; - double slope = 0.7; - double factor = 40.0; - if(TMath::Abs(y[0]) < 16.2/slope) result = y[0]/slope; - else if(TMath::Abs(y[0]) < 49.8/slope ) result=y[0]/slope-TMath::Sign(1.0,y[0])*factor; - else if(TMath::Abs(y[0]) < 85.2/slope ) result=y[0]/slope-TMath::Sign(1.0,y[0])*factor*2; - else result=y[0]/slope-TMath::Sign(1.0,y[0])*factor*2; - return result; - -} - -/*void testmodel() { - TF1 eqline("x","x",-200,200); - eqline.Draw(""); - eqline.SetLineStyle(kDashed); - - //TF1 f1("model",model,-200,200,2); - TF1 f1("model_inv",model_invert,-200,200,2); - eqline.SetNpx(10000); - f1.SetNpx(10000); - std::vector pars = {0.0,1.}; - f1.SetParameters(pars.data()); - f1.SetLineColor(kGreen+2); - f1.SetLineStyle(kLine); - f1.Draw("L SAME"); - - gPad->Modified(); gPad->Update(); - while(gPad->WaitPrimitive()); - -}*/ diff --git a/scratch/test_eloss.h b/scratch/test_eloss.h deleted file mode 100644 index e69de29..0000000 diff --git a/pccal/pc_gm_coeffs.dat b/slope_intercept_results_17F.dat old mode 100755 new mode 100644 similarity index 100% rename from pccal/pc_gm_coeffs.dat rename to slope_intercept_results_17F.dat