diff --git a/MakeVertex.C b/MakeVertex.C index 2e5662d..ace016e 100755 --- a/MakeVertex.C +++ b/MakeVertex.C @@ -1328,21 +1328,21 @@ Bool_t MakeVertex::Process(Long64_t entry) double Ex_from_proton = apkin_p.getExc(sx3Efix, theta_recon * 180. / M_PI); double Ex_from_alpha = apkin_a.getExc(sx3Efixalpha, theta_recon * 180. / M_PI); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3", 400, 0, 30, 800, 0, 40000, sx3Efix, apSumE * sinTheta, nA_label); + plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3", 800, 0, 30, 1600, 0, 20000, sx3Efix, apSumE * sinTheta, nA_label); plotter->Fill1D(nA_label + "_Ex_from_alphas_SX3" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); plotter->Fill1D(nA_label + "_Ex_from_protons_SX3" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); plotter->Fill2D(nA_label + "_sx3_E_vs_theta_raw_SX3", 180, 0, 180, 400, 0, 30, theta_recon * 180. / M_PI, sx3event.Energy1, nA_label); plotter->Fill2D(nA_label + "_sx3_E_vs_theta_corr_SX3", 180, 0, 180, 400, 0, 30, theta_recon * 180. / M_PI, sx3Efix, nA_label); - if (vtx_gate != "") - { - plotter->Fill1D(nA_label + "_twisted_pcz_recon_SX3" + vtx_gate, 600, -300, 300, pcz_intersect.Z(), nA_label); - plotter->Fill1D(nA_label + "_twisted_vertex_recon_SX3" + vtx_gate, 600, -300, 300, vertex_recon, nA_label); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3" + vtx_gate, 400, 0, 30, 800, 0, 40000, sx3Efix, apSumE * sinTheta, nA_label); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3_alpha" + vtx_gate, 400, 0, 30, 800, 0, 40000, sx3Efixalpha, apSumE * sinTheta, nA_label); - plotter->Fill1D(nA_label + "_Ex_from_alphas_SX3" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); - plotter->Fill1D(nA_label + "_Ex_from_protons_SX3" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); - } + // if (vtx_gate != "") + // { + // plotter->Fill1D(nA_label + "_twisted_pcz_recon_SX3" + vtx_gate, 600, -300, 300, pcz_intersect.Z(), nA_label); + // plotter->Fill1D(nA_label + "_twisted_vertex_recon_SX3" + vtx_gate, 600, -300, 300, vertex_recon, nA_label); + // plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3" + vtx_gate, 400, 0, 30, 800, 0, 40000, sx3Efix, apSumE * sinTheta, nA_label); + // plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_SX3_alpha" + vtx_gate, 400, 0, 30, 800, 0, 40000, sx3Efixalpha, apSumE * sinTheta, nA_label); + // plotter->Fill1D(nA_label + "_Ex_from_alphas_SX3" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); + // plotter->Fill1D(nA_label + "_Ex_from_protons_SX3" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); + // } } } #endif // DIAG_nA_SX3 @@ -1391,21 +1391,21 @@ Bool_t MakeVertex::Process(Long64_t entry) double Ex_from_proton = apkin_p.getExc(qqqEfix, theta_recon * 180. / M_PI); double Ex_from_alpha = apkin_a.getExc(qqqEfixalpha, theta_recon * 180. / M_PI); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ", 400, 0, 30, 800, 0, 40000, qqqEfix, apSumE * sinTheta, nA_label); + plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ", 800, 0, 30, 800, 0, 20000, qqqEfix, apSumE * sinTheta, nA_label); plotter->Fill1D(nA_label + "_Ex_from_alphas_QQQ" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); plotter->Fill1D(nA_label + "_Ex_from_protons_QQQ" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); plotter->Fill2D(nA_label + "_qqq_E_vs_theta_raw_QQQ", 180, 0, 180, 400, 0, 30, theta_recon * 180. / M_PI, qqqevent.Energy1, nA_label); plotter->Fill2D(nA_label + "_qqq_E_vs_theta_corr_QQQ", 180, 0, 180, 400, 0, 30, theta_recon * 180. / M_PI, qqqEfix, nA_label); - if (vtx_gate != "") - { - plotter->Fill1D(nA_label + "_twisted_pcz_recon_QQQ" + vtx_gate, 600, -300, 300, pcz_intersect.Z(), nA_label); - plotter->Fill1D(nA_label + "_twisted_vertex_recon_QQQ" + vtx_gate, 600, -300, 300, vertex_recon, nA_label); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ" + vtx_gate, 400, 0, 30, 800, 0, 40000, qqqEfix, apSumE * sinTheta, nA_label); - plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ_alpha" + vtx_gate, 400, 0, 30, 800, 0, 40000, qqqEfixalpha, apSumE * sinTheta, nA_label); - plotter->Fill1D(nA_label + "_Ex_from_alphas_QQQ" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); - plotter->Fill1D(nA_label + "_Ex_from_protons_QQQ" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); - } + // if (vtx_gate != "") + // { + // plotter->Fill1D(nA_label + "_twisted_pcz_recon_QQQ" + vtx_gate, 600, -300, 300, pcz_intersect.Z(), nA_label); + // plotter->Fill1D(nA_label + "_twisted_vertex_recon_QQQ" + vtx_gate, 600, -300, 300, vertex_recon, nA_label); + // plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ" + vtx_gate, 400, 0, 30, 800, 0, 40000, qqqEfix, apSumE * sinTheta, nA_label); + // plotter->Fill2D(nA_label + "_dE_Ecorr_Anode_QQQ_alpha" + vtx_gate, 400, 0, 30, 800, 0, 40000, qqqEfixalpha, apSumE * sinTheta, nA_label); + // plotter->Fill1D(nA_label + "_Ex_from_alphas_QQQ" + vtx_gate, 200, -10, 10, Ex_from_alpha, nA_label); + // plotter->Fill1D(nA_label + "_Ex_from_protons_QQQ" + vtx_gate, 200, -10, 10, Ex_from_proton, nA_label); + // } } } #endif // DIAG_nA_QQQ diff --git a/anasen_fem/run.py b/anasen_fem/run.py index 1e886a6..43e22d7 100755 --- a/anasen_fem/run.py +++ b/anasen_fem/run.py @@ -3,7 +3,7 @@ import os # val=-174.3 val=0 -count=11 +count=10 while val<174.3+0.1: print(val) os.system("python3 wires_gmsh2d_bc.py "+str(val)) diff --git a/anasen_fem/scalars.dat.names b/anasen_fem/scalars.dat.names index a1c6859..0d58be4 100755 --- a/anasen_fem/scalars.dat.names +++ b/anasen_fem/scalars.dat.names @@ -2,7 +2,7 @@ Metadata for SaveScalars file: ./scalars.dat Elmer version: 26.2 Elmer compilation date: 2026-05-14 Solver input file: wires2d.sif -File started at: 2026/05/19 16:24:56 +File started at: 2026/05/19 22:16:54 Variables in columns of matrix: 1: res: potential difference diff --git a/run_sx3.sh b/run_sx3.sh index fe6e156..20680c6 100755 --- a/run_sx3.sh +++ b/run_sx3.sh @@ -21,7 +21,7 @@ fi #alpha+gas 27Al export DATASET="27Al" #root -q -b -x ../ANASEN_analysis/data/27Al_Data/Run_009_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run09.root; -if [[ 1 -eq 1 ]]; then +if [[ 1 -eq 0 ]]; then #export timecut_low=230.0; export timecut_low=400.0; #export timecut_high=400.0; @@ -84,7 +84,7 @@ export source_vertex=-24.96; root -q -l -b -x ../ANASEN_analysis/data/17F_Data/S export source_vertex=-73.96; root -q -l -b -x ../ANASEN_analysis/data/17F_Data/SourceRun_021_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run21.root; fi #17F reaction data -if [[ 1 -eq 1 ]]; then +if [[ 1 -eq 0 ]]; then export source_vertex=-57.28; root -q -l -b -x ../ANASEN_analysis/data/17F_Data/ProtonRun_035_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run35.root; #export source_vertex=-8.28; root -q -l -b -x ../ANASEN_analysis/data/17F_Data/ProtonRun_036_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root resulrs_run36.root; #export source_vertex=-27.88; root -q -l -b -x ../ANASEN_analysis/data/17F_Data/ProtonRun_037_mapped.root -e 'tree->Process("MakeVertex.C+O")'; mv Analyzer_SX3.root results_run37.root; diff --git a/slope_intercept_results.dat b/slope_intercept_results.dat new file mode 100644 index 0000000..a29898b --- /dev/null +++ b/slope_intercept_results.dat @@ -0,0 +1,49 @@ +#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