diff --git a/cmake_modules/Cuts.cmake b/cmake_modules/Cuts.cmake index de8c5aa1..62ec87a1 100644 --- a/cmake_modules/Cuts.cmake +++ b/cmake_modules/Cuts.cmake @@ -7,7 +7,7 @@ set(Cuts_INCLUDE_DIR ${Cuts_INSTALL_DIR}/include) set(Cuts_LIBRARY_DIR ${Cuts_INSTALL_DIR}/lib) ExternalProject_Add(Cuts_Ext - GIT_REPOSITORY "https://git.cbm.gsi.de/pwg-c2f/analysis/cuts.git" + GIT_REPOSITORY "https://github.com/HeavyIonAnalysis/AnalysisTreeCuts" GIT_TAG "master" UPDATE_DISCONNECTED ${UPDATE_DISCONNECTED} SOURCE_DIR "${EXTERNAL_DIR}/Cuts_src" diff --git a/cmake_modules/QnTools.cmake b/cmake_modules/QnTools.cmake index 8f78f526..50360ad4 100644 --- a/cmake_modules/QnTools.cmake +++ b/cmake_modules/QnTools.cmake @@ -5,7 +5,7 @@ find_package(QnTools QUIET) include(FetchContent) FetchContent_Declare(QnTools GIT_REPOSITORY https://github.com/HeavyIonAnalysis/QnTools.git - GIT_TAG "v0.0.1" + GIT_TAG "f39fe1fd34d9b01037e28e369a25c9fe03f0b4e4" # PR: Fix bootstrapping error propagation of scaling and pow operations UPDATE_DISCONNECTED ${UPDATE_DISCONNECTED} ) diff --git a/setups/hades-analysis-config.yml b/setups/hades-analysis-config.yml deleted file mode 100644 index 68d7de21..00000000 --- a/setups/hades-analysis-config.yml +++ /dev/null @@ -1,81 +0,0 @@ -_layouts: - -_axes: - - &axis_pT - name: mdc_vtx_tracks/pT - nb: 16 - lo: 0.0 - hi: 1.6 - - &axis_rapidity - name: mdc_vtx_tracks/rapidity - nb: 15 - lo: -0.1 - hi: 1.49 - -hades_analysis: - event-variables: - - event_header/selected_tof_rpc_hits_centrality - axes: - - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } - q-vectors: - - name: u - type: track - phi: mdc_vtx_tracks/phi - weight: Ones - corrections: - - recentering - - twist-and-rescale - axes: - - *axis_pT - - *axis_rapidity - cuts: - mdc_vtx_tracks/pid: { equals: 2212 } - qa: - - {name: mdc_vtx_tracks/phi, nb: 100, lo: -4., hi: 4.} - - *axis_pT - - *axis_rapidity - - [*axis_rapidity, *axis_pT] - - name: M - type: track - phi: mdc_vtx_tracks/phi - weight: Ones - corrections: - - recentering - - twist-and-rescale - axes: - - *axis_rapidity - cuts: - mdc_vtx_tracks/pid: { equals: 2212 } - mdc_vtx_tracks/pT: { range: [0.0, 1.6] } - qa: - - {name: mdc_vtx_tracks/phi, nb: 100, lo: -4., hi: 4.} - - *axis_pT - - *axis_rapidity - - [*axis_rapidity, *axis_pT] - - name: W1 - type: track - phi: forward_wall_hits/phi - weight: Ones - corrections: - - recentering - - twist-and-rescale - cuts: - forward_wall_hits/ring: { range: [1.0, 5.0] } - - name: W2 - type: track - phi: forward_wall_hits/phi - weight: Ones - corrections: - - recentering - - twist-and-rescale - cuts: - forward_wall_hits/ring: { range: [6.0, 7.0] } - - name: W3 - type: track - phi: forward_wall_hits/phi - weight: Ones - corrections: - - recentering - - twist-and-rescale - cuts: - forward_wall_hits/ring: { range: [8.0, 10.0] } \ No newline at end of file diff --git a/setups/hades/au123-correction-debug.yml b/setups/hades/au123-correction-debug.yml new file mode 100644 index 00000000..588a50a0 --- /dev/null +++ b/setups/hades/au123-correction-debug.yml @@ -0,0 +1,154 @@ +_layouts: + +_axes: + - &axis_rapidity + name: mdc_vtx_tracks_rapidity/rapidity + nb: 30 + lo: -0.01 + hi: 1.49 + - &axis_pT_non_uniform + name: mdc_vtx_tracks_rapidity/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_pT_uniform + name: mdc_vtx_tracks_rapidity/pT + nb: 16 + lo: 0.0 + hi: 1.6 +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_rapidity/pid: { equals: 2212 } + - &good_track_cuts + mdc_vtx_tracks_rapidity/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_rapidity/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_rapidity/chi2: {range: [0.0, 100.0]} + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: non_uniform_none + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: none + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_non_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_non_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_non_uniform] + ########################################################################### + - name: non_uniform_m + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_non_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_non_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_non_uniform] + ########################################################################### + - name: uniform_none + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: none + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_uniform] + ########################################################################### + - name: uniform_m + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_uniform] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag + ############################################################# + - name: Mf + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mf_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mb_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} diff --git a/setups/hades/correction-agag-123.yml b/setups/hades/correction-agag-123.yml new file mode 100644 index 00000000..847720b8 --- /dev/null +++ b/setups/hades/correction-agag-123.yml @@ -0,0 +1,205 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0.2, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 1.5] + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pi_pos_cuts + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pi_neg_cuts + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &r1_cuts + forward_wall_hits/rnd_sub: { equals: 0 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &r2_cuts + forward_wall_hits/rnd_sub: { equals: 1 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/time: { range: [22.681788, 27.002128] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/time: { range: [22.681788, 26.684456] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header_extra/selected_tof_rpc_hits_centrality + axes: + - { name: event_header_extra/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + - name: protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: + *protons_cuts + cuts-quality: + *good_particle_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + - name: pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pions_pT + - *axis_pions_rapidity + cuts-pid: + *pi_pos_cuts + cuts-quality: + *good_particle_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pions_pT + - *axis_pions_rapidity + - [*axis_pions_rapidity, *axis_pions_pT] +# - name: pi_neg +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: mdc_vtx_tracks_extra/efficiency +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_pions_pT +# - *axis_pions_rapidity +# cuts-pid: +# *pi_neg_cuts +# cuts-quality: +# *good_particle_cuts +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_pions_pT +# - *axis_pions_rapidity +# - [*axis_pions_rapidity, *axis_pions_pT] + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mf_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mb_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - name: R1 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts: +# *r1_cuts +# - name: R2 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts: +# *r2_cuts + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w1_cuts + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w2_cuts + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w3_cuts \ No newline at end of file diff --git a/setups/hades/correction-agag-158.yml b/setups/hades/correction-agag-158.yml new file mode 100644 index 00000000..12e600df --- /dev/null +++ b/setups/hades/correction-agag-158.yml @@ -0,0 +1,220 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0.2, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 1.5] + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pi_pos_cuts + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pi_neg_cuts + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.17, 1.37]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.27, 0.47]} + - &r1_cuts + forward_wall_hits/rnd_sub: { equals: 0 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &r2_cuts + forward_wall_hits/rnd_sub: { equals: 1 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/time: { range: [22.681788, 27.002128] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/time: { range: [22.681788, 26.684456] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ############################################################# + - name: protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: + *protons_cuts + cuts-quality: + *good_particle_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ############################################################# +# - name: pi_pos +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: mdc_vtx_tracks_extra/efficiency +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_pions_pT +# - *axis_pions_rapidity +# cuts-pid: +# *pi_pos_cuts +# cuts-quality: +# *good_particle_cuts +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_pions_pT +# - *axis_pions_rapidity +# - [*axis_pions_rapidity, *axis_pions_pT] +# ############################################################# +# - name: pi_neg +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: mdc_vtx_tracks_extra/efficiency +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_pions_pT +# - *axis_pions_rapidity +# cuts-pid: +# *pi_neg_cuts +# cuts-quality: +# *good_particle_cuts +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_pions_pT +# - *axis_pions_rapidity +# - [*axis_pions_rapidity, *axis_pions_pT] + ############################################################# + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mf_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mb_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# +# - name: R1 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts: +# *r1_cuts +# ############################################################# +# - name: R2 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts: +# *r2_cuts + ############################################################# + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w1_cuts + ############################################################# + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w2_cuts + ############################################################# + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w3_cuts +# - name: wall_plain +# type: psi +# phi: event_header/fw_angle_plain +# weight: Ones +# norm: m \ No newline at end of file diff --git a/setups/hades/correction-auau-123-all.yml b/setups/hades/correction-auau-123-all.yml new file mode 100644 index 00000000..46240ef8 --- /dev/null +++ b/setups/hades/correction-auau-123-all.yml @@ -0,0 +1,136 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_protons_theta + name: mdc_vtx_tracks_extra/theta + nb: 12 + lo: 0.3 + hi: 1.5 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/beta: { range: [0.84, 1.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/beta: { range: [0.85, 1.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/beta: { range: [0.80, 1.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ############################################################# + - name: protons + type: track + phi: mdc_vtx_tracks_extra/phi +# weight: mdc_vtx_tracks_extra/efficiency + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: +# - *axis_protons_pT +# - *axis_protons_rapidity + - *axis_protons_theta + cuts: + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ############################################################# + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mf_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mb_cuts + cuts-quality: + *good_particle_cuts + cuts-pid: + *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w1_cuts + ############################################################# + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w2_cuts + ############################################################# + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts: + *w3_cuts \ No newline at end of file diff --git a/setups/hades/correction-auau-123-mult.yml b/setups/hades/correction-auau-123-mult.yml new file mode 100644 index 00000000..802e05f0 --- /dev/null +++ b/setups/hades/correction-auau-123-mult.yml @@ -0,0 +1,118 @@ +_layouts: + +_axes: + - &axis_pT + name: mdc_vtx_tracks/pT + nb: 20 + lo: 0.0 + hi: 2.0 + - &axis_rapidity + name: mdc_vtx_tracks/rapidity + nb: 15 + lo: -0.01 + hi: 1.49 +_cuts: + - &protons_cuts + mdc_vtx_tracks/pid: { equals: 2212 } + mdc_vtx_tracks/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks/pid: { equals: 2212 } + mdc_vtx_tracks/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks/rapidity: {range: [1.09, 1.29]} + mdc_vtx_tracks/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks/chi2: {range: [0.0, 100.0]} + - &mb_cuts + mdc_vtx_tracks/pid: { equals: 2212 } + mdc_vtx_tracks/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks/rapidity: {range: [0.19, 0.39]} + mdc_vtx_tracks/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/beta: { range: [0.84, 1.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/beta: { range: [0.85, 1.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/beta: { range: [0.80, 1.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits + axes: + - { name: event_header/selected_tof_rpc_hits, bin-edges: [0, 10, 17, 23, 29, 36, 45, 55, 68, 82, 99, 117, 136, 157, 180, 312, 499] } + q-vectors: + - name: u + type: track + phi: mdc_vtx_tracks/phi + weight: Ones + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT + - *axis_rapidity + cuts: + *protons_cuts + qa: + - {name: mdc_vtx_tracks/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT + - *axis_rapidity + - [*axis_rapidity, *axis_pT] + - name: Mf + type: track + phi: mdc_vtx_tracks/phi + weight: Ones + corrections: + - recentering + - twist-and-rescale + cuts: + *mf_cuts + qa: + - {name: mdc_vtx_tracks/phi, nb: 100, lo: -4., hi: 4.} + - name: Mb + type: track + phi: mdc_vtx_tracks/phi + weight: Ones + corrections: + - recentering + - twist-and-rescale + cuts: + *mb_cuts + qa: + - {name: mdc_vtx_tracks/phi, nb: 100, lo: -4., hi: 4.} + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + corrections: + - recentering + - twist-and-rescale + cuts: + *w1_cuts + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + corrections: + - recentering + - twist-and-rescale + cuts: + *w2_cuts + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + corrections: + - recentering + - twist-and-rescale + cuts: + *w3_cuts \ No newline at end of file diff --git a/setups/hades/correction-auau-123-occupancy.yml b/setups/hades/correction-auau-123-occupancy.yml new file mode 100644 index 00000000..3c97871d --- /dev/null +++ b/setups/hades/correction-auau-123-occupancy.yml @@ -0,0 +1,171 @@ +_layouts: + +_axes: + ############## + # Protons axes + ############## + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_protons_c_eff + name: mdc_vtx_tracks_extra/c_eff + nb: 10 + lo: -0.00025 + hi: 0.00475 + ############ + # Pions axes + ############ + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_pions_c_eff + name: mdc_vtx_tracks_extra/c_eff + nb: 10 + lo: -0.0175 + hi: 0.0475 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pi_neg_cuts + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &pi_pos_cuts + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/beta: { range: [0.84, 1.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/beta: { range: [0.85, 1.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/beta: { range: [0.80, 1.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ############################################################# + - name: protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/occ_weight + norm: m + axes: + - *axis_protons_pT + - *axis_protons_rapidity + - *axis_protons_c_eff + cuts-quality: *good_particle_cuts + cuts-pid: *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: pi_neg + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/occ_weight + norm: m + axes: + - *axis_pions_pT + - *axis_pions_rapidity + - *axis_pions_c_eff + cuts-quality: *good_particle_cuts + cuts-pid: *pi_neg_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/occ_weight + norm: m + axes: + - *axis_pions_pT + - *axis_pions_rapidity + - *axis_pions_c_eff + cuts-quality: *good_particle_cuts + cuts-pid: *pi_pos_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-subevent: *mf_cuts + cuts-quality: *good_particle_cuts + cuts-pid: *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-subevent: *mb_cuts + cuts-quality: *good_particle_cuts + cuts-pid: *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w1_cuts + ############################################################# + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w2_cuts + ############################################################# + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w3_cuts \ No newline at end of file diff --git a/setups/hades/correction-auau-123-plains.yml b/setups/hades/correction-auau-123-plains.yml new file mode 100644 index 00000000..bfe26f1d --- /dev/null +++ b/setups/hades/correction-auau-123-plains.yml @@ -0,0 +1,117 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 10 + lo: 0 + hi: 2.0 + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pi_pos_cuts + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pi_neg_cuts + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &r1_cuts + forward_wall_hits/rnd_sub: { equals: 0 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &r2_cuts + forward_wall_hits/rnd_sub: { equals: 1 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &wall_cuts + forward_wall_hits/rnd_sub: { equals: 1 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/beta: { range: [0.84, 1.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/beta: { range: [0.85, 1.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/beta: { range: [0.80, 1.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ############################################################# + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w1_cuts + ############################################################# + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w2_cuts + ############################################################# + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w3_cuts + ############################################################# + - name: Wall + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *wall_cuts \ No newline at end of file diff --git a/setups/hades/correction-auau-123.yml b/setups/hades/correction-auau-123.yml new file mode 100644 index 00000000..d8e09594 --- /dev/null +++ b/setups/hades/correction-auau-123.yml @@ -0,0 +1,223 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 10 + lo: 0 + hi: 2.0 + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 +_cuts: + - &protons_cuts + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pi_pos_cuts + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pi_neg_cuts + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &good_particle_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &r1_cuts + forward_wall_hits/rnd_sub: { equals: 0 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &r2_cuts + forward_wall_hits/rnd_sub: { equals: 1 } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/beta: { range: [0.84, 1.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/beta: { range: [0.85, 1.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/beta: { range: [0.80, 1.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ############################################################# + - name: charged + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m +# corrections: +# - recentering +# - twist-and-rescale + axes: + - *axis_theta + - *axis_uniform_pT + cuts-track-quality: { mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} } + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m +# corrections: +# - recentering +# - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: *protons_cuts + cuts-quality: *good_particle_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ############################################################# + - name: pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m +# corrections: +# - recentering +# - twist-and-rescale + axes: + - *axis_pions_pT + - *axis_pions_rapidity + cuts-quality: *good_particle_cuts + cuts-pid: *pi_pos_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: pi_neg + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m +# corrections: +# - recentering +# - twist-and-rescale + axes: + - *axis_pions_pT + - *axis_pions_rapidity + cuts-quality: *good_particle_cuts + cuts-pid: *pi_neg_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: pi2 + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-quality: *good_particle_cuts + cuts-pid: { mdc_vtx_tracks_extra/abs_pid: {equals: 211} } + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: p2 + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-quality: *good_particle_cuts + cuts-pid: { mdc_vtx_tracks_extra/pid: {equals: 2212} } + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-subevent: *mf_cuts + cuts-quality: *good_particle_cuts + cuts-pid: *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts-subevent: *mb_cuts + cuts-quality: *good_particle_cuts + cuts-pid: *protons_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: W1 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w1_cuts + ############################################################# + - name: W2 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w2_cuts + ############################################################# + - name: W3 + type: track + phi: forward_wall_hits/phi + weight: forward_wall_hits/signal + norm: m + corrections: + - recentering +# - twist-and-rescale + cuts: *w3_cuts \ No newline at end of file diff --git a/setups/hades/correction-mcini.yml b/setups/hades/correction-mcini.yml new file mode 100644 index 00000000..0ce82ef3 --- /dev/null +++ b/setups/hades/correction-mcini.yml @@ -0,0 +1,40 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: particles/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_protons_rapidity + name: particles/rapidity + nb: 15 + lo: -0.75 + hi: 0.75 + +hades_analysis: + event-variables: + - event_header/centrality + axes: + - { name: event_header/centrality, nb: 6, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: protons + type: track + phi: particles/phi + weight: Ones + norm: m + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: + particles/pid: { equals: 2212 } + qa: + - {name: particles/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: event_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/correlation-agag-123.yml b/setups/hades/correlation-agag-123.yml new file mode 100644 index 00000000..49c22775 --- /dev/null +++ b/setups/hades/correlation-agag-123.yml @@ -0,0 +1,104 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] +# - name: pi_pos +# tags: [ un_vector ] +# - name: pi_neg +# tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector, mdc ] + - name: Mb + tags: [ qn_vector, mdc ] + +_axes: + - ¢rality + name: event_header_extra_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + - &ep_components + [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: +# Scalar product + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u1Q1" + axes: [ *centrality ] +# - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] +# n-samples: 100 +# weights-type: observable +# folder: "/SP/u2Q1Q1" +# axes: [ *centrality ] + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/Q1Q1" + axes: [ *centrality ] + # Event plane + # - args: [*ta_u, *ta_qn_ep] + # n-samples: 100 + # weights-type: observable + # folder: "/EP/u1Q1" + # axes: [ *centrality ] + # - args: [*ta_u2, *ta_qn_ep, *ta_qn_ep] + # n-samples: 100 + # weights-type: observable + # folder: "/EP/u2Q1Q1" + # axes: [ *centrality ] + # - args: [*ta_qn_ep, *ta_qn_ep] + # n-samples: 100 + # weights-type: reference + # folder: "/EP/Q1Q1" + # axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-agag-158.yml b/setups/hades/correlation-agag-158.yml new file mode 100644 index 00000000..18f900bf --- /dev/null +++ b/setups/hades/correlation-agag-158.yml @@ -0,0 +1,87 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] + - name: pi_pos + tags: [ un_vector ] + - name: pi_neg + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector, mdc ] + - name: Mb + tags: [ qn_vector, mdc ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + - &ep_components + [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u1Q1" + axes: [ *centrality ] +# - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] +# n-samples: 100 +# weights-type: observable +# folder: "/SP/u2Q1Q1" +# axes: [ *centrality ] + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/Q1Q1" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-all.yml b/setups/hades/correlation-all.yml new file mode 100644 index 00000000..0549db80 --- /dev/null +++ b/setups/hades/correlation-all.yml @@ -0,0 +1,62 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/uQ" + axes: [ *centrality ] + - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u2Q1Q1" + axes: [ *centrality ] + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/QQ" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-auau-123-occupancy.yml b/setups/hades/correlation-auau-123-occupancy.yml new file mode 100644 index 00000000..862cf1a7 --- /dev/null +++ b/setups/hades/correlation-auau-123-occupancy.yml @@ -0,0 +1,54 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] + - name: pi_neg + tags: [ un_vector ] + - name: pi_pos + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector, mdc ] + - name: Mb + tags: [ qn_vector, mdc ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components [ x1,y1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + # Two particles correlations + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u1Q1" + axes: [ *centrality ] + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/Q1Q1" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-auau-123-pi-neg.yml b/setups/hades/correlation-auau-123-pi-neg.yml new file mode 100644 index 00000000..bec11cf4 --- /dev/null +++ b/setups/hades/correlation-auau-123-pi-neg.yml @@ -0,0 +1,94 @@ +_detectors: &detectors + - name: pi_neg + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + - &ep_components + [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/uQ" + axes: [ *centrality ] +# - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] +# n-samples: 100 +# weights-type: observable +# folder: "/SP/u2Q1Q1" +# axes: [ *centrality ] + - args: [*ta_u, *ta_qn_ep] + n-samples: 100 + weights-type: observable + folder: "/EP/uQ" + axes: [ *centrality ] + # Reference x Reference + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/QQ" + axes: [ *centrality ] + - args: [*ta_qn_ep, *ta_qn_ep] + n-samples: 100 + weights-type: reference + folder: "/EP/QQ" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-auau-123-pi-pos.yml b/setups/hades/correlation-auau-123-pi-pos.yml new file mode 100644 index 00000000..fb445c2d --- /dev/null +++ b/setups/hades/correlation-auau-123-pi-pos.yml @@ -0,0 +1,94 @@ +_detectors: &detectors + - name: pi_pos + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + - &ep_components + [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/uQ" + axes: [ *centrality ] +# - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] +# n-samples: 100 +# weights-type: observable +# folder: "/SP/u2Q1Q1" +# axes: [ *centrality ] + - args: [*ta_u, *ta_qn_ep] + n-samples: 100 + weights-type: observable + folder: "/EP/uQ" + axes: [ *centrality ] + # Reference x Reference + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/QQ" + axes: [ *centrality ] + - args: [*ta_qn_ep, *ta_qn_ep] + n-samples: 100 + weights-type: reference + folder: "/EP/QQ" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-auau-123-proton.yml b/setups/hades/correlation-auau-123-proton.yml new file mode 100644 index 00000000..e8be46c3 --- /dev/null +++ b/setups/hades/correlation-auau-123-proton.yml @@ -0,0 +1,94 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + - &ep_components + [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/uQ" + axes: [ *centrality ] + - args: [*ta_u2, *ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u2Q1Q1" + axes: [ *centrality ] +# - args: [*ta_u, *ta_qn_ep] +# n-samples: 100 +# weights-type: observable +# folder: "/EP/uQ" +# axes: [ *centrality ] + # Reference x Reference + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/QQ" + axes: [ *centrality ] +# - args: [*ta_qn_ep, *ta_qn_ep] +# n-samples: 100 +# weights-type: reference +# folder: "/EP/QQ" +# axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-auau-123.yml b/setups/hades/correlation-auau-123.yml new file mode 100644 index 00000000..fcc776ca --- /dev/null +++ b/setups/hades/correlation-auau-123.yml @@ -0,0 +1,102 @@ +_detectors: &detectors + - name: charged + tags: [ un_vector ] + - name: protons + tags: [ un_vector ] + - name: pi_pos + tags: [ un_vector ] + - name: pi_neg + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector, mdc ] + - name: Mb + tags: [ qn_vector, mdc ] + - name: pi2 + tags: [ q2_vector ] + - name: p2 + tags: [ q2_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components [ x1,y1 ] + - &sp2_components [ x2,y2 ] + - &ep_components [ cos1,sin1 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_q2 + query: { tags: { any-in: [ q2_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ recentered ] + weight: ones + - &ta_qn_sp + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ recentered ] + weight: ones + - &ta_qn_ep + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_sp + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + - &ta_rnd_ep + query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u1Q1" + axes: [ *centrality ] + - args: [*ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/Q1Q1" + axes: [ *centrality ] + # Three particles correlations + - args: [*ta_u, *ta_q2, *ta_qn_sp] + n-samples: 100 + weights-type: observable + folder: "/SP/u1Q2Q1" + axes: [ *centrality ] + - args: [*ta_q2, *ta_qn_sp, *ta_qn_sp] + n-samples: 100 + weights-type: reference + folder: "/SP/Q2Q1Q1" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/correlation-mcini.yml b/setups/hades/correlation-mcini.yml new file mode 100644 index 00000000..f20ea407 --- /dev/null +++ b/setups/hades/correlation-mcini.yml @@ -0,0 +1,56 @@ +_detectors: &detectors + - name: protons + tags: [ un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_centrality + nb: 6 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + +_tasks_args: + - &ta_u + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_u2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_u, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/protons" + axes: [ *centrality ] + - args: [*ta_u2, *ta_psi2] + n-samples: 100 + weights-type: observable + folder: "/u2Q2/protons" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/correlation/hades.yml b/setups/hades/correlation-mult.yml similarity index 85% rename from setups/correlation/hades.yml rename to setups/hades/correlation-mult.yml index 2672c0ca..74f6ce3f 100644 --- a/setups/correlation/hades.yml +++ b/setups/hades/correlation-mult.yml @@ -11,16 +11,17 @@ _detectors: &detectors - name: W3 tags: [ qn_vector ] correction-step: rescaled - - name: M + - name: Mf + tags: [ qn_vector ] + correction-step: rescaled + - name: Mb tags: [ qn_vector ] correction-step: rescaled _axes: - - ¢rality - name: event_header_selected_tof_rpc_hits_centrality - nb: 8 - lo: 0 - hi: 40 + - &multiplicity + name: event_header_selected_tof_rpc_hits + bin-edges: [0, 10, 17, 23, 29, 36, 45, 55, 68, 82, 99, 117, 136, 157, 180, 312, 499] _components: - &sp_components @@ -42,12 +43,12 @@ _tasks: weight: ones # temporarily ignored # actions: [xx,xy,yx,yy] - n-samples: 50 + n-samples: 100 weights-type: observable # temporarily ignored # weights-function: xx folder: "/uQ/SP" - axes: [ *centrality ] + axes: [ *multiplicity ] - args: - query: { tags: { any-in: [ un_vector ] } } query-list: *detectors @@ -56,17 +57,17 @@ _tasks: weight: sumw - query: { tags: { any-in: [ qn_vector ] } } query-list: *detectors - components: *ep_components + components: *ep_components # correction-steps: [ rescaled ] weight: ones # temporarily ignored # actions: [xx,xy,yx,yy] - n-samples: 50 + n-samples: 100 weights-type: observable # temporarily ignored # weights-function: xx folder: "/uQ/EP" - axes: [ *centrality ] + axes: [ *multiplicity ] # Reference x Reference - args: - query: { tags: { any-in: [ qn_vector ] } } @@ -81,11 +82,11 @@ _tasks: correction-steps: [ rescaled ] # temporarily ignored # actions: [xx,xy,yx,yy] - n-samples: 50 + n-samples: 100 weights-type: reference # weights-function: xx folder: "/QQ/SP" - axes: [ *centrality ] + axes: [ *multiplicity ] - args: - query: { tags: { any-in: [ qn_vector ] } } query-list: *detectors @@ -99,8 +100,8 @@ _tasks: correction-steps: [ rescaled ] # temporarily ignored # actions: [xx,xy,yx,yy] - n-samples: 50 + n-samples: 100 weights-type: reference # weights-function: xx folder: "/QQ/EP" - axes: [ *centrality ] \ No newline at end of file + axes: [ *multiplicity ] \ No newline at end of file diff --git a/setups/hades/correlation.yml b/setups/hades/correlation.yml new file mode 100644 index 00000000..2172a896 --- /dev/null +++ b/setups/hades/correlation.yml @@ -0,0 +1,145 @@ +_detectors: &detectors + - name: u + tags: [ un_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: R1 + tags: [ rnd_vector ] + - name: R2 + tags: [ rnd_vector ] + - name: Mf_protons + tags: [ qn_vector ] + - name: Mb_protons + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &ep_components + [ cos1,sin1 ] + +_tasks: + - args: + - query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + n-samples: 100 + weights-type: observable + folder: "/uQ/SP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + n-samples: 100 + weights-type: observable + folder: "/uQ/SP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components # + correction-steps: [ rescaled ] + weight: ones + n-samples: 100 + weights-type: observable + folder: "/uQ/EP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components # + correction-steps: [ rescaled ] + weight: ones + n-samples: 100 + weights-type: observable + folder: "/uQ/EP" + axes: [ *centrality ] + # Reference x Reference + - args: + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + n-samples: 100 + weights-type: reference + folder: "/QQ/SP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + n-samples: 100 + weights-type: reference + folder: "/QQ/SP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + - query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + n-samples: 100 + weights-type: reference + folder: "/QQ/EP" + axes: [ *centrality ] + - args: + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + - query: { tags: { any-in: [ rnd_vector ] } } + query-list: *detectors + components: *ep_components + correction-steps: [ rescaled ] + n-samples: 100 + weights-type: reference + folder: "/QQ/EP" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correction-agag-158.yml b/setups/hades/mc-correction-agag-158.yml new file mode 100644 index 00000000..89c0b753 --- /dev/null +++ b/setups/hades/mc-correction-agag-158.yml @@ -0,0 +1,172 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_gen_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_gen_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/time: { range: [22.681788, 27.002128] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/time: { range: [22.681788, 26.684456] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/time: { range: [22.681788, 28.352235] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: pid_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ########################################################################### + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mf_cuts + cuts-quality: + *good_track_cuts + cuts-pid: + *pid_protons_cut + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-subevent: + *mb_cuts + cuts-quality: + *good_track_cuts + cuts-pid: + *pid_protons_cut + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: gen_prim_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ########################################################################### + - name: W1 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-w1: + *w1_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: W2 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-w2: + *w2_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: W3 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + - twist-and-rescale + cuts-w3: + *w3_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag +# qa: +# - {name: sim_header/reaction_plane, nb: 320, lo: -3.2, hi: 3.2} \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-all.yml b/setups/hades/mc-correction-auau-123-all.yml new file mode 100644 index 00000000..41306508 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-all.yml @@ -0,0 +1,69 @@ +_layouts: + +_axes: + - &axis_protons_theta + name: mdc_vtx_tracks_extra/theta + nb: 48 + lo: 0.3 + hi: 1.5 + - &axis_gen_protons_theta + name: sim_particles_extra/theta + nb: 48 + lo: 0.3 + hi: 1.5 + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_gen_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_gen_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### +# - name: accepted +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: Ones +# norm: none +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_protons_theta +# cuts: +# mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + ########################################################################### + - name: protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + cuts: + sim_particles_extra/is_charged: {equals: 1} + axes: + - *axis_gen_protons_theta + ########################################################################### + - name: total + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + cuts: + sim_particles_extra/is_charged: {equals: 1} + sim_particles_extra/theta: {range: [0.3, 1.5]} + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-even.yml b/setups/hades/mc-correction-auau-123-even.yml new file mode 100644 index 00000000..6e142b25 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-even.yml @@ -0,0 +1,80 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_gen_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/abs_ycm + nb: 7 + lo: 0.05 + hi: 0.75 + - &axis_gen_protons_rapidity + name: sim_particles_extra/abs_ycm + nb: 7 + lo: 0.05 + hi: 0.75 + +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: pid_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ############################################################# + - name: gen_prim_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-occupancy.yml b/setups/hades/mc-correction-auau-123-occupancy.yml new file mode 100644 index 00000000..239fcd00 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-occupancy.yml @@ -0,0 +1,227 @@ +_layouts: + +_axes: + - &axis_rec_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_rec_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_tru_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_tru_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_rec_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_rec_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_tru_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_tru_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_rec_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_tru_theta + name: sim_particles_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_rec_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 10 + lo: 0.0 + hi: 2.0 + - &axis_tru_uniform_pT + name: sim_particles_extra/pT + nb: 10 + lo: 0.0 + hi: 2.0 + +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ########################################################################### + - name: rec_charged + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + axes: + - *axis_rec_theta + - *axis_rec_uniform_pT + cuts-quality: { mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} } + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_theta + - [*axis_rec_theta, *axis_rec_protons_pT] + ########################################################################### + - name: rec_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + cuts-pid: *pid_protons_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + - [*axis_rec_protons_rapidity, *axis_rec_protons_pT] + ########################################################################### + - name: rec_pi_neg + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + cuts-pid: *pid_pi_neg_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: rec_pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + cuts-pid: *pid_pi_pos_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: tru_charged + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_theta + - *axis_tru_uniform_pT + cuts-charge: { sim_particles_extra/is_charged: {range: [0.0, 100.0]} } + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: tru_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + cuts-pid: *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: tru_pi_neg + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + cuts-pid: *gen_prim_pi_neg_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################### + - name: tru_pi_pos + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + cuts-pid: *gen_prim_pi_pos_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-pT-theta.yml b/setups/hades/mc-correction-auau-123-pT-theta.yml new file mode 100644 index 00000000..13ea5d8f --- /dev/null +++ b/setups/hades/mc-correction-auau-123-pT-theta.yml @@ -0,0 +1,227 @@ +_layouts: + +_axes: + - &axis_rec_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_rec_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_tru_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_tru_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_rec_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_rec_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_tru_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_tru_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_rec_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_tru_theta + name: sim_particles_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_rec_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + - &axis_tru_uniform_pT + name: sim_particles_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ########################################################################### + - name: rec_charged + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + axes: + - *axis_rec_theta + - *axis_rec_uniform_pT + cuts-quality: { mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} } + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_theta + - [*axis_rec_theta, *axis_rec_protons_pT] + ########################################################################### + - name: rec_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_theta + - *axis_rec_uniform_pT + cuts-pid: *pid_protons_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + - [*axis_rec_protons_rapidity, *axis_rec_protons_pT] + ########################################################################### + - name: rec_pi_neg + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_theta + - *axis_rec_uniform_pT + cuts-pid: *pid_pi_neg_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: rec_pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_theta + - *axis_rec_uniform_pT + cuts-pid: *pid_pi_pos_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: tru_charged + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_theta + - *axis_tru_uniform_pT + cuts-charge: { sim_particles_extra/is_charged: {range: [0.0, 100.0]} } + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: tru_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_theta + - *axis_tru_uniform_pT + cuts-pid: *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: tru_pi_neg + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_theta + - *axis_tru_uniform_pT + cuts-pid: *gen_prim_pi_neg_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################### + - name: tru_pi_pos + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_theta + - *axis_tru_uniform_pT + cuts-pid: *gen_prim_pi_pos_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-plains.yml b/setups/hades/mc-correction-auau-123-plains.yml new file mode 100644 index 00000000..2678424a --- /dev/null +++ b/setups/hades/mc-correction-auau-123-plains.yml @@ -0,0 +1,104 @@ +_layouts: + +_axes: + - &axis_rec_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_rec_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_tru_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_tru_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_rec_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_rec_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_tru_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_tru_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_rec_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_tru_theta + name: sim_particles_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_rec_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 10 + lo: 0.0 + hi: 2.0 + - &axis_tru_uniform_pT + name: sim_particles_extra/pT + nb: 10 + lo: 0.0 + hi: 2.0 + +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-true.yml b/setups/hades/mc-correction-auau-123-true.yml new file mode 100644 index 00000000..87554442 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-true.yml @@ -0,0 +1,146 @@ +_layouts: + +_axes: + - &axis_gen_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_gen_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_gen_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_gen_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 +_cuts: + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &theta_cut + sim_particles_extra/theta: { range: [0.3, 1.5] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ############################################################# + - name: proton + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# + - name: proton_theta + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + cuts-pid: + *gen_prim_protons_cut + cuts-theta: + *theta_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# + - name: pi_pos + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_pions_pT + - *axis_gen_pions_rapidity + cuts-pid: + *gen_prim_pi_pos_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# + - name: pi_pos_theta + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_pions_pT + - *axis_gen_pions_rapidity + cuts-pid: + *gen_prim_pi_pos_cut + cuts-theta: + *theta_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# + - name: pi_neg + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_pions_pT + - *axis_gen_pions_rapidity + cuts-pid: + *gen_prim_pi_neg_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# + - name: pi_neg_theta + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_pions_pT + - *axis_gen_pions_rapidity + cuts-pid: + *gen_prim_pi_neg_cut + cuts-theta: + *theta_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-wall.yml b/setups/hades/mc-correction-auau-123-wall.yml new file mode 100644 index 00000000..4083ef82 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-wall.yml @@ -0,0 +1,235 @@ +_layouts: + +_axes: + - &axis_rec_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_rec_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_tru_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_tru_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_rec_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_rec_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_tru_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_tru_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + reconstructed_forward_wall_hits/ring: { range: [1.0, 5.0] } + reconstructed_forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + reconstructed_forward_wall_hits/ring: { range: [6.0, 7.0] } + reconstructed_forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + reconstructed_forward_wall_hits/ring: { range: [8.0, 10.0] } + reconstructed_forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ########################################################################### + - name: rec_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + cuts-pid: *pid_protons_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + - [*axis_rec_protons_rapidity, *axis_rec_protons_pT] + ########################################################################### + - name: rec_pi_neg + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + cuts-pid: *pid_pi_neg_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: rec_pi_pos + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + cuts-pid: *pid_pi_pos_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_pions_pT + - *axis_rec_pions_rapidity + - [*axis_rec_pions_rapidity, *axis_rec_pions_pT] + ########################################################################### + - name: tru_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + cuts-pid: *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: tru_pi_neg + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + cuts-pid: *gen_prim_pi_neg_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################### + - name: tru_pi_pos + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + cuts-pid: *gen_prim_pi_pos_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_pions_pT + - *axis_tru_pions_rapidity + - [*axis_tru_pions_rapidity, *axis_tru_pions_pT] + ########################################################################## + - name: Mf + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + cuts-subevent: *mf_cuts + cuts-quality: *good_track_cuts + cuts-pid: *pid_protons_cut + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: Mb + type: track + phi: mdc_vtx_tracks_extra/phi + weight: Ones + norm: m + corrections: + - recentering + cuts-subevent: *mb_cuts + cuts-quality: *good_track_cuts + cuts-pid: *pid_protons_cut + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: W1 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + cuts-w1: *w1_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: W2 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + cuts-w2: *w2_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: W3 + type: track + phi: reconstructed_forward_wall_hits/phi + weight: reconstructed_forward_wall_hits/signal + norm: m + corrections: + - recentering + cuts-w3: *w3_cuts + qa: + - {name: reconstructed_forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123-y-pT-theta.yml b/setups/hades/mc-correction-auau-123-y-pT-theta.yml new file mode 100644 index 00000000..2cddb692 --- /dev/null +++ b/setups/hades/mc-correction-auau-123-y-pT-theta.yml @@ -0,0 +1,137 @@ +_layouts: + +_axes: + - &axis_rec_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_rec_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_tru_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.3, 0.35, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.95, 1.5] + - &axis_tru_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.10, 0.13, 0.15, 0.18, 0.21, 0.25, 0.32, 0.55, 1.0] + - &axis_rec_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_rec_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_tru_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_tru_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_rec_theta + name: mdc_vtx_tracks_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_tru_theta + name: sim_particles_extra/theta + nb: 14 + lo: 0.2 + hi: 1.6 + - &axis_rec_uniform_pT + name: mdc_vtx_tracks_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + - &axis_tru_uniform_pT + name: sim_particles_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 8, lo: 0.0, hi: 40 } + q-vectors: + ########################################################################### + - name: rec_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/efficiency + norm: m + axes: + - *axis_rec_uniform_pT + - *axis_rec_protons_rapidity + - *axis_rec_theta + cuts-pid: *pid_protons_cut + cuts-quality: *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_rec_protons_pT + - *axis_rec_protons_rapidity + - [*axis_rec_protons_rapidity, *axis_rec_protons_pT] + ########################################################################### + - name: tru_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_tru_uniform_pT + - *axis_tru_protons_rapidity + - *axis_tru_theta + cuts-pid: *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_tru_protons_pT + - *axis_tru_protons_rapidity + - [*axis_tru_protons_rapidity, *axis_tru_protons_pT] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correction-auau-123.yml b/setups/hades/mc-correction-auau-123.yml new file mode 100644 index 00000000..d22794ce --- /dev/null +++ b/setups/hades/mc-correction-auau-123.yml @@ -0,0 +1,306 @@ +_layouts: + +_axes: + - &axis_protons_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_uniform_protons_pT + name: mdc_vtx_tracks_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + - &axis_pions_pT + name: mdc_vtx_tracks_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_uniform_pions_pT + name: mdc_vtx_tracks_extra/pT + nb: 10 + lo: 0.0 + hi: 1.0 + - &axis_gen_protons_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_uniform_gen_protons_pT + name: sim_particles_extra/pT + nb: 15 + lo: 0.0 + hi: 1.5 + - &axis_gen_pions_pT + name: sim_particles_extra/pT + bin-edges: [0, 0.08, 0.105, 0.13, 0.155, 0.18, 0.21, 0.25, 0.315, 0.535, 1.0] + - &axis_uniform_gen_pions_pT + name: sim_particles_extra/pT + nb: 10 + lo: 0.0 + hi: 1.0 + - &axis_protons_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_pions_rapidity + name: mdc_vtx_tracks_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_gen_protons_rapidity + name: sim_particles_extra/ycm + nb: 15 + lo: -0.75 + hi: 0.75 + - &axis_gen_pions_rapidity + name: sim_particles_extra/ycm + nb: 17 + lo: -0.65 + hi: 1.05 + - &axis_theta + name: mdc_vtx_tracks_extra/theta + nb: 12 + lo: 0.3 + hi: 1.5 + - &axis_gen_theta + name: sim_particles_extra/theta + nb: 12 + lo: 0.3 + hi: 1.5 +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_extra/pid: { equals: 2212 } + - &pid_pi_pos_cut + mdc_vtx_tracks_extra/pid: { equals: 211 } + - &pid_pi_neg_cut + mdc_vtx_tracks_extra/pid: { equals: -211 } + - &mf_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [1.09, 1.29]} + - &mb_cuts + mdc_vtx_tracks_extra/pT: {range: [0.0, 2.0]} + mdc_vtx_tracks_extra/rapidity: {range: [0.19, 0.39]} + - &gen_prim_protons_cut + sim_particles_extra/pid: { equals: 2212 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_pos_cut + sim_particles_extra/pid: { equals: 211 } + sim_particles_extra/is_primary: { equals: 1 } + - &gen_prim_pi_neg_cut + sim_particles_extra/pid: { equals: -211 } + sim_particles_extra/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_extra/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_extra/chi2: {range: [0.0, 100.0]} + - &w1_cuts + forward_wall_hits/ring: { range: [1.0, 5.0] } + forward_wall_hits/signal: { range: [80.0, 999.0] } + - &w2_cuts + forward_wall_hits/ring: { range: [6.0, 7.0] } + forward_wall_hits/signal: { range: [85.0, 999.0] } + - &w3_cuts + forward_wall_hits/ring: { range: [8.0, 10.0] } + forward_wall_hits/signal: { range: [88.0, 999.0] } + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: pid_protons + type: track + phi: mdc_vtx_tracks_extra/phi + weight: mdc_vtx_tracks_extra/sector_weight +# weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_protons_pT + - *axis_protons_rapidity +# - *axis_theta + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_protons_pT + - *axis_protons_rapidity + - [*axis_protons_rapidity, *axis_protons_pT] + ########################################################################## +# - name: pid_pi_pos +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: Ones +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_pions_pT +# - *axis_pions_rapidity +# - *axis_theta +# cuts-pid: +# *pid_pi_pos_cut +# cuts-quality: +# *good_track_cuts +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_protons_pT +# - *axis_protons_rapidity +# - [*axis_protons_rapidity, *axis_protons_pT] +# ########################################################################### +# - name: pid_pi_neg +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: Ones +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# axes: +# - *axis_pions_pT +# - *axis_pions_rapidity +# - *axis_theta +# cuts-pid: +# *pid_pi_neg_cut +# cuts-quality: +# *good_track_cuts +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_protons_pT +# - *axis_protons_rapidity +# - [*axis_protons_rapidity, *axis_protons_pT] +# ########################################################################### +# - name: Mf +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: Ones +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts-subevent: +# *mf_cuts +# cuts-quality: +# *good_track_cuts +# cuts-pid: +# *pid_protons_cut +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} +# ############################################################# +# - name: Mb +# type: track +# phi: mdc_vtx_tracks_extra/phi +# weight: Ones +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts-subevent: +# *mb_cuts +# cuts-quality: +# *good_track_cuts +# cuts-pid: +# *pid_protons_cut +# qa: +# - {name: mdc_vtx_tracks_extra/phi, nb: 100, lo: -4., hi: 4.} + ############################################################# + - name: gen_prim_protons + type: track + phi: sim_particles_extra/phi + weight: Ones + norm: m + axes: + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity +# - *axis_gen_theta + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_protons_pT + - *axis_gen_protons_rapidity + - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ############################################################# +# - name: gen_prim_pi_pos +# type: track +# phi: sim_particles_extra/phi +# weight: Ones +# norm: m +# axes: +# - *axis_gen_pions_pT +# - *axis_gen_pions_rapidity +# - *axis_gen_theta +# cuts-pid: +# *gen_prim_pi_pos_cut +# qa: +# - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_gen_protons_pT +# - *axis_gen_protons_rapidity +# - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] +# ############################################################# +# - name: gen_prim_pi_neg +# type: track +# phi: sim_particles_extra/phi +# weight: Ones +# norm: m +# axes: +# - *axis_gen_pions_pT +# - *axis_gen_pions_rapidity +# - *axis_gen_theta +# cuts-pid: +# *gen_prim_pi_neg_cut +# qa: +# - {name: sim_particles_extra/phi, nb: 100, lo: -4., hi: 4.} +# - *axis_gen_protons_pT +# - *axis_gen_protons_rapidity +# - [*axis_gen_protons_rapidity, *axis_gen_protons_pT] + ########################################################################### +# - name: W1 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts-w1: +# *w1_cuts +# qa: +# - {name: forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} +# ########################################################################### +# - name: W2 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts-w2: +# *w2_cuts +# qa: +# - {name: forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} +# ########################################################################### +# - name: W3 +# type: track +# phi: forward_wall_hits/phi +# weight: forward_wall_hits/signal +# norm: m +# corrections: +# - recentering +# - twist-and-rescale +# cuts-w3: +# *w3_cuts +# qa: +# - {name: forward_wall_hits/phi, nb: 100, lo: -4., hi: 4.} + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag +# qa: +# - {name: sim_header/reaction_plane, nb: 320, lo: -3.2, hi: 3.2} \ No newline at end of file diff --git a/setups/hades/mc-correction-debug.yml b/setups/hades/mc-correction-debug.yml new file mode 100644 index 00000000..ff674871 --- /dev/null +++ b/setups/hades/mc-correction-debug.yml @@ -0,0 +1,200 @@ +_layouts: + +_axes: + - &axis_rapidity + name: mdc_vtx_tracks_rapidity/rapidity + nb: 15 + lo: -0.01 + hi: 1.49 + - &axis_gen_rapidity + name: sim_tracks/rapidity + nb: 1 + lo: -0.01 + hi: 1.49 + - &axis_pT_non_uniform + name: mdc_vtx_tracks_rapidity/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_pT_uniform + name: mdc_vtx_tracks_rapidity/pT + nb: 16 + lo: 0.0 + hi: 1.6 + - &axis_gen_pT_non_uniform + name: sim_tracks/pT + bin-edges: [0, 0.29375, 0.35625, 0.41875, 0.48125, 0.54375, 0.61875, 0.70625, 0.81875, 1.01875, 2.0] + - &axis_gen_pT_uniform + name: sim_tracks/pT + nb: 16 + lo: 0.0 + hi: 1.6 +_cuts: + - &pid_protons_cut + mdc_vtx_tracks_rapidity/pid: { equals: 2212 } + - &gen_prim_protons_cut + sim_tracks/pid: { equals: 2212 } + sim_tracks/is_primary: { equals: 1 } + - &good_track_cuts + mdc_vtx_tracks_rapidity/dca_xy: {range: [-10.0, 10.0]} + mdc_vtx_tracks_rapidity/dca_z: {range: [-10.0, 10.0]} + mdc_vtx_tracks_rapidity/chi2: {range: [0.0, 100.0]} + +hades_analysis: + event-variables: + - event_header/selected_tof_rpc_hits_centrality + axes: + - { name: event_header/selected_tof_rpc_hits_centrality, nb: 12, lo: 0.0, hi: 60 } + q-vectors: + ########################################################################### + - name: rec_non_uniform_none + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: none + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_non_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_non_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_non_uniform] + ########################################################################### + - name: rec_non_uniform_m + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_non_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_non_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_non_uniform] + ########################################################################### + - name: rec_uniform_none + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: none + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_uniform] + ########################################################################### + - name: rec_uniform_m + type: track + phi: mdc_vtx_tracks_rapidity/phi + weight: Ones + norm: m + corrections: + - recentering + - twist-and-rescale + axes: + - *axis_pT_uniform + - *axis_rapidity + cuts-pid: + *pid_protons_cut + cuts-quality: + *good_track_cuts + qa: + - {name: mdc_vtx_tracks_rapidity/phi, nb: 100, lo: -4., hi: 4.} + - *axis_pT_uniform + - *axis_rapidity + - [*axis_rapidity, *axis_pT_uniform] + ############################################################# + - name: gen_non_uniform_none + type: track + phi: sim_tracks/phi + weight: Ones + norm: none + axes: + - *axis_gen_pT_non_uniform + - *axis_gen_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_tracks/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_pT_non_uniform + - *axis_gen_rapidity + - [*axis_gen_rapidity, *axis_gen_pT_non_uniform] + ############################################################# + - name: gen_non_uniform_m + type: track + phi: sim_tracks/phi + weight: Ones + norm: m + axes: + - *axis_gen_pT_non_uniform + - *axis_gen_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_tracks/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_pT_non_uniform + - *axis_gen_rapidity + - [*axis_gen_rapidity, *axis_gen_pT_non_uniform] + ############################################################# + - name: gen_uniform_none + type: track + phi: sim_tracks/phi + weight: Ones + norm: none + axes: + - *axis_gen_pT_uniform + - *axis_gen_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_tracks/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_pT_uniform + - *axis_gen_rapidity + - [*axis_gen_rapidity, *axis_gen_pT_uniform] + ############################################################# + - name: gen_uniform_m + type: track + phi: sim_tracks/phi + weight: Ones + norm: m + axes: + - *axis_gen_pT_uniform + - *axis_gen_rapidity + cuts-pid: + *gen_prim_protons_cut + qa: + - {name: sim_tracks/phi, nb: 100, lo: -4., hi: 4.} + - *axis_gen_pT_uniform + - *axis_gen_rapidity + - [*axis_gen_rapidity, *axis_gen_pT_uniform] + ########################################################################### + - name: psi_rp + type: psi + phi: sim_header/reaction_plane + weight: Ones + norm: mag \ No newline at end of file diff --git a/setups/hades/mc-correlation-all.yml b/setups/hades/mc-correlation-all.yml new file mode 100644 index 00000000..7a676d8d --- /dev/null +++ b/setups/hades/mc-correlation-all.yml @@ -0,0 +1,80 @@ +_detectors: &detectors + - name: protons + tags: [ gen_un_vector ] + - name: total + tags: [ gen_un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + +_tasks_args: + - &ta_accepted + query: { tags: { any-in: [ rec_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_accepted2 + query: { tags: { any-in: [ rec_un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_all + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_all2 + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + +_tasks: +# - args: [*ta_accepted, *ta_psi] +# n-samples: 100 +# weights-type: observable +# folder: "/uQ/reco" +# axes: [ *centrality ] +# - args: [*ta_accepted2, *ta_psi2] +# n-samples: 100 +# weights-type: observable +# folder: "/u2Q2/reco" +# axes: [ *centrality ] + - args: [*ta_all, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/sim" + axes: [ *centrality ] + - args: [*ta_all2, *ta_psi2] + n-samples: 100 + weights-type: observable + folder: "/u2Q2/sim" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-auau-123-occupancy.yml b/setups/hades/mc-correlation-auau-123-occupancy.yml new file mode 100644 index 00000000..d406eb3c --- /dev/null +++ b/setups/hades/mc-correlation-auau-123-occupancy.yml @@ -0,0 +1,79 @@ +_detectors: &detectors + - name: rec_charged + tags: [ rec_vector ] + - name: rec_protons + tags: [ rec_vector ] + - name: rec_pi_pos + tags: [ rec_vector ] + - name: rec_pi_neg + tags: [ rec_vector ] + - name: tru_charged + tags: [ tru_vector ] + - name: tru_protons + tags: [ tru_vector ] + - name: tru_pi_pos + tags: [ tru_vector ] + - name: tru_pi_neg + tags: [ tru_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components [ x1,y1 ] + - &sp2_components [ x2,y2 ] + +_tasks_args: + - &ta_rec + query: { tags: { any-in: [ rec_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_rec2 + query: { tags: { any-in: [ rec_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_tru + query: { tags: { any-in: [ tru_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + - &ta_tru2 + query: { tags: { any-in: [ tru_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_rec, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/rec" + axes: [ *centrality ] + - args: [*ta_tru, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/tru" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-auau-123-wall.yml b/setups/hades/mc-correlation-auau-123-wall.yml new file mode 100644 index 00000000..f5b4c2da --- /dev/null +++ b/setups/hades/mc-correlation-auau-123-wall.yml @@ -0,0 +1,99 @@ +_detectors: &detectors +# ************************** +# particles vectors + - name: rec_protons + tags: [ rec_vector ] + - name: rec_pi_pos + tags: [ rec_vector ] + - name: rec_pi_neg + tags: [ rec_vector ] + - name: tru_protons + tags: [ tru_vector ] + - name: tru_pi_pos + tags: [ tru_vector ] + - name: tru_pi_neg + tags: [ tru_vector ] +# ************************** +# reference vectors + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] +# ************************** +# tru reaction plane + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components [ x1,y1 ] + +_tasks_args: + - &ta_rec + query: { tags: { any-in: [ rec_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_tru + query: { tags: { any-in: [ tru_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_qn + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_rec, *ta_qn] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/rec" + axes: [ *centrality ] + - args: [*ta_rec, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/rec" + axes: [ *centrality ] + - args: [*ta_tru, *ta_qn] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/tru" + axes: [ *centrality ] + - args: [*ta_tru, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/tru" + axes: [ *centrality ] + - args: [*ta_qn, *ta_qn] + n-samples: 100 + weights-type: reference + folder: "/Q1Q1" + axes: [ *centrality ] + - args: [*ta_qn, *ta_psi] + n-samples: 100 + weights-type: reference + folder: "/Q1Q1" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-auau-123-y-pT-theta.yml b/setups/hades/mc-correlation-auau-123-y-pT-theta.yml new file mode 100644 index 00000000..c0ad1ed4 --- /dev/null +++ b/setups/hades/mc-correlation-auau-123-y-pT-theta.yml @@ -0,0 +1,67 @@ +_detectors: &detectors + - name: rec_protons + tags: [ rec_vector ] + - name: tru_protons + tags: [ tru_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 8 + lo: 0 + hi: 40 + +_components: + - &sp_components [ x1,y1 ] + - &sp2_components [ x2,y2 ] + +_tasks_args: + - &ta_rec + query: { tags: { any-in: [ rec_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_rec2 + query: { tags: { any-in: [ rec_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_tru + query: { tags: { any-in: [ tru_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + - &ta_tru2 + query: { tags: { any-in: [ tru_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_rec, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/rec" + axes: [ *centrality ] + - args: [*ta_tru, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/u1Q1/tru" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-debug.yml b/setups/hades/mc-correlation-debug.yml new file mode 100644 index 00000000..af5cc5c0 --- /dev/null +++ b/setups/hades/mc-correlation-debug.yml @@ -0,0 +1,62 @@ +_detectors: &detectors + - name: rec_non_uniform_none + tags: [ un_vector ] + - name: rec_non_uniform_m + tags: [ un_vector ] + - name: rec_uniform_none + tags: [ un_vector ] + - name: rec_uniform_m + tags: [ un_vector ] + - name: gen_non_uniform_none + tags: [ gen_un_vector ] + - name: gen_non_uniform_m + tags: [ gen_un_vector ] + - name: gen_uniform_none + tags: [ gen_un_vector ] + - name: gen_uniform_m + tags: [ gen_un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_reco + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain, recentered, rescaled ] + weight: sumw + - &ta_sim + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_reco, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/reco" + axes: [ *centrality ] + - args: [*ta_sim, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/gen" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-even.yml b/setups/hades/mc-correlation-even.yml new file mode 100644 index 00000000..9621f3a3 --- /dev/null +++ b/setups/hades/mc-correlation-even.yml @@ -0,0 +1,50 @@ +_detectors: &detectors + - name: pid_protons + tags: [ un_vector ] + - name: gen_prim_protons + tags: [ gen_un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_reco + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_sim + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_reco, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/reco" + axes: [ *centrality ] + - args: [*ta_sim, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/gen" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-gen-pi-neg.yml b/setups/hades/mc-correlation-gen-pi-neg.yml new file mode 100644 index 00000000..9000ff39 --- /dev/null +++ b/setups/hades/mc-correlation-gen-pi-neg.yml @@ -0,0 +1,37 @@ +_detectors: &detectors + - name: gen_prim_pi_neg + tags: [ un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_gen + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_gen, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/gen" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-gen-pi-pos.yml b/setups/hades/mc-correlation-gen-pi-pos.yml new file mode 100644 index 00000000..6e4ced20 --- /dev/null +++ b/setups/hades/mc-correlation-gen-pi-pos.yml @@ -0,0 +1,37 @@ +_detectors: &detectors + - name: gen_prim_pi_pos + tags: [ un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_gen + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_gen, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/gen" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-rec-pi-neg.yml b/setups/hades/mc-correlation-rec-pi-neg.yml new file mode 100644 index 00000000..b851d853 --- /dev/null +++ b/setups/hades/mc-correlation-rec-pi-neg.yml @@ -0,0 +1,37 @@ +_detectors: &detectors + - name: pid_pi_neg + tags: [ un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_reco + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_reco, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/reco" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-rec-pi-pos.yml b/setups/hades/mc-correlation-rec-pi-pos.yml new file mode 100644 index 00000000..811d79ce --- /dev/null +++ b/setups/hades/mc-correlation-rec-pi-pos.yml @@ -0,0 +1,37 @@ +_detectors: &detectors + - name: pid_pi_pos + tags: [ un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + +_tasks_args: + - &ta_reco + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_reco, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/reco" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation-true.yml b/setups/hades/mc-correlation-true.yml new file mode 100644 index 00000000..4f5c84a7 --- /dev/null +++ b/setups/hades/mc-correlation-true.yml @@ -0,0 +1,65 @@ +_detectors: &detectors + - name: proton + tags: [ gen_un_vector ] + - name: proton_theta + tags: [ gen_un_vector ] + - name: pi_pos + tags: [ gen_un_vector ] + - name: pi_pos_theta + tags: [ gen_un_vector ] + - name: pi_neg + tags: [ gen_un_vector ] + - name: pi_neg_theta + tags: [ gen_un_vector ] + - name: psi_rp + tags: [ psi_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + +_tasks_args: + - &ta_sim + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + - &ta_sim2 + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + +_tasks: + - args: [*ta_sim, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/gen" + axes: [ *centrality ] + - args: [*ta_sim2, *ta_psi2] + n-samples: 100 + weights-type: observable + folder: "/u2Q2/gen" + axes: [ *centrality ] \ No newline at end of file diff --git a/setups/hades/mc-correlation.yml b/setups/hades/mc-correlation.yml new file mode 100644 index 00000000..c3a1cff7 --- /dev/null +++ b/setups/hades/mc-correlation.yml @@ -0,0 +1,128 @@ +_detectors: &detectors + - name: pid_protons + tags: [ un_vector ] +# - name: pid_pi_pos +# tags: [ un_vector ] +# - name: pid_pi_neg +# tags: [ un_vector ] + - name: gen_prim_protons + tags: [ gen_un_vector ] +# - name: gen_prim_pi_pos +# tags: [ gen_un_vector ] +# - name: gen_prim_pi_neg +# tags: [ gen_un_vector ] + - name: psi_rp + tags: [ psi_vector ] + - name: W1 + tags: [ qn_vector ] + - name: W2 + tags: [ qn_vector ] + - name: W3 + tags: [ qn_vector ] + - name: Mf + tags: [ qn_vector ] + - name: Mb + tags: [ qn_vector ] + +_axes: + - ¢rality + name: event_header_selected_tof_rpc_hits_centrality + nb: 12 + lo: 0 + hi: 60 + +_components: + - &sp_components + [ x1,y1 ] + - &sp2_components + [ x2,y2 ] + +_tasks_args: + - &ta_reco + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: sumw + - &ta_reco2 + query: { tags: { any-in: [ un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_sim + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + - &ta_sim2 + query: { tags: { any-in: [ gen_un_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: sumw + - &ta_psi + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ plain ] + weight: ones + - &ta_psi2 + query: { tags: { any-in: [ psi_vector ] } } + query-list: *detectors + components: *sp2_components + correction-steps: [ plain ] + weight: ones + - &ta_qn + query: { tags: { any-in: [ qn_vector ] } } + query-list: *detectors + components: *sp_components + correction-steps: [ rescaled ] + weight: ones + +_tasks: + - args: [*ta_reco, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/reco" + axes: [ *centrality ] +# - args: [*ta_reco2, *ta_psi] +# n-samples: 100 +# weights-type: observable +# folder: "/u2Q1/reco" +# axes: [ *centrality ] + - args: [*ta_reco2, *ta_psi2] + n-samples: 100 + weights-type: observable + folder: "/u2Q2/reco" + axes: [ *centrality ] +# - args: [*ta_reco, *ta_qn] +# n-samples: 100 +# weights-type: observable +# folder: "/uQ/reco" +# axes: [ *centrality ] + - args: [*ta_sim, *ta_psi] + n-samples: 100 + weights-type: observable + folder: "/uQ/gen" + axes: [ *centrality ] +# - args: [*ta_sim2, *ta_psi] +# n-samples: 100 +# weights-type: observable +# folder: "/u2Q1/gen" +# axes: [ *centrality ] + - args: [*ta_sim2, *ta_psi2] + n-samples: 100 + weights-type: observable + folder: "/u2Q2/gen" + axes: [ *centrality ] +# - args: [*ta_qn, *ta_psi] +# n-samples: 100 +# weights-type: reference +# folder: "/QQ" +# axes: [ *centrality ] +# - args: [*ta_qn, *ta_qn] +# n-samples: 100 +# weights-type: reference +# folder: "/QQ" +# axes: [ *centrality ] \ No newline at end of file diff --git a/setups/observables-all.yml b/setups/observables-all.yml new file mode 100644 index 00000000..ab5fed7f --- /dev/null +++ b/setups/observables-all.yml @@ -0,0 +1,35 @@ +_detectors: &detectors + - name: W1 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W2 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W3 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: Mf + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + - name: Mb + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + +_observables: + - method: three-sub-event + QQ-directory: SP/QQ + uQ-directory: SP/uQ + u-vector: {name: protons, components: [x1, y1], correction-step: RESCALED} + EP-vectors: ep-vectors + resolution-vectors: res-vectors +# - method: three-sub-event +# QQ-directory: QQ/SP +# uQ-directory: uQ/SP +# u-vector: {name: pions, components: [x1, y1], correction-step: RESCALED} +# EP-vectors: ep-vectors +# resolution-vectors: res-vectors \ No newline at end of file diff --git a/setups/observables-auau.yml b/setups/observables-auau.yml new file mode 100644 index 00000000..d7b0ee46 --- /dev/null +++ b/setups/observables-auau.yml @@ -0,0 +1,41 @@ +_detectors: &detectors + - name: W1 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W2 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W3 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: Mf + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + - name: Mb + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + +_observables: + - method: three-sub-event + QQ-directory: QQ/SP + uQ-directory: uQ/SP + u-vector: {name: protons, components: [x1, y1], correction-step: RESCALED} + EP-vectors: ep-vectors + resolution-vectors: res-vectors + - method: three-sub-event + QQ-directory: QQ/SP + uQ-directory: uQ/SP + u-vector: {name: pi_plus, components: [x1, y1], correction-step: RESCALED} + EP-vectors: ep-vectors + resolution-vectors: res-vectors + - method: three-sub-event + QQ-directory: QQ/SP + uQ-directory: uQ/SP + u-vector: {name: pi_minus, components: [x1, y1], correction-step: RESCALED} + EP-vectors: ep-vectors + resolution-vectors: res-vectors \ No newline at end of file diff --git a/setups/observables-calculation-config.yml b/setups/observables-calculation-config.yml new file mode 100644 index 00000000..99a09e9f --- /dev/null +++ b/setups/observables-calculation-config.yml @@ -0,0 +1,29 @@ +_detectors: &detectors + - name: W1 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W2 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: W3 + correction-step: RESCALED + components: [x1,y1] + tags: [ep-vectors, res-vectors] + - name: Mf + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + - name: Mb + correction-step: RESCALED + components: [x1,y1] + tags: [res-vectors] + +_observables: + - method: three-sub-event + QQ-directory: QQ/SP + uQ-directory: uQ/SP + u-vector: {name: pi_plus, components: [x1, y1], correction-step: RESCALED} + EP-vectors: ep-vectors + resolution-vectors: res-vectors \ No newline at end of file diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt index e02dd849..0137b8b8 100644 --- a/src/CMakeLists.txt +++ b/src/CMakeLists.txt @@ -2,6 +2,7 @@ add_subdirectory(QnAnalysisBase) add_subdirectory(QnAnalysisConfig) add_subdirectory(QnAnalysisCorrect) add_subdirectory(QnAnalysisCorrelate) +add_subdirectory(QnAnalysisObservables) #add_subdirectory(correlate) #add_subdirectory(tasks) #add_subdirectory(analyze) diff --git a/src/QnAnalysisCorrelate/CorrelationTaskRunner.cpp b/src/QnAnalysisCorrelate/CorrelationTaskRunner.cpp index c5ee9643..30a09969 100644 --- a/src/QnAnalysisCorrelate/CorrelationTaskRunner.cpp +++ b/src/QnAnalysisCorrelate/CorrelationTaskRunner.cpp @@ -63,7 +63,7 @@ void Qn::Analysis::Correlate::CorrelationTaskRunner::Run() { container_meta.String() = GenCorrelationMeta(correlation); dir->WriteObject(&container, correlation.meta_key.c_str()); - dir->WriteObject(&container_meta, (correlation.meta_key + "_meta").c_str()); +// dir->WriteObject(&container_meta, (correlation.meta_key + "_meta").c_str()); } } diff --git a/src/QnAnalysisObservables/CMakeLists.txt b/src/QnAnalysisObservables/CMakeLists.txt new file mode 100644 index 00000000..7f36b032 --- /dev/null +++ b/src/QnAnalysisObservables/CMakeLists.txt @@ -0,0 +1,12 @@ +find_package(Boost REQUIRED COMPONENTS program_options) + +set(SOURCES + FileManager.cpp + Method.cpp + MethodOf3SE.cpp + V1Observables.cpp + MethodOfRS.cpp + ConvertConfigs.cpp) + +add_executable(BuildObservables Calculate.cpp ${SOURCES}) +target_link_libraries(BuildObservables QnTools::Base Boost::program_options ${ROOT_LIBRARIES} yaml-cpp) \ No newline at end of file diff --git a/src/QnAnalysisObservables/Calculate.cpp b/src/QnAnalysisObservables/Calculate.cpp new file mode 100644 index 00000000..56b13c09 --- /dev/null +++ b/src/QnAnalysisObservables/Calculate.cpp @@ -0,0 +1,103 @@ +// +// Created by mikhail on 10/23/20. +// + +#include "ConvertConfigs.hpp" +#include "FileManager.hpp" +#include "V1Observables.hpp" +#include +#include + +int main(int argv, char** argc){ + namespace po = boost::program_options; + std::string input_file; + std::string output_file_name{"output.root"}; + std::string config_file; + po::options_description options("Options"); + options.add_options() + ("help,h", "Help screen") + ("MC,m", "MC data sample") + ("input,i", po::value(&input_file), "Input file list") + ("output,o", po::value(&output_file_name), "Name of output file") + ("config,c", po::value(&config_file),"path to yaml-config file"); + po::variables_map vm; + po::parsed_options parsed = po::command_line_parser(argv, argc).options(options).run(); + po::store(parsed, vm); + po::notify(vm); + if (vm.count("help")){ + std::cout << options << std::endl; + return 0; + } + FileManager::OpenFile( input_file ); + auto config = YAML::LoadFile( config_file ); + auto detectors_configs = config["_detectors"].as>(); + auto observables_configs = config["_observables"].as>(); + std::vector observables; + observables.reserve(observables_configs.size()); + + auto out_file = TFile::Open(output_file_name.c_str(), "recreate" ); + for( const auto &observable_config : observables_configs ){ + observables.emplace_back( observable_config.method ); + observables.back().SetQqCorrelationsDirectory(observable_config.qq_correlations_directory); + observables.back().SetUqCorrelationsDirectory(observable_config.uq_correlations_directory); + observables.back().SetUvector( observable_config.u_vector ); + std::vector ep_vectors; + for( const auto& vector : detectors_configs ) + if( std::count(vector.tags.begin(), vector.tags.end(), observable_config.ep_vectors_tag) > 0 ) + ep_vectors.emplace_back(vector); + observables.back().SetEPvectors(ep_vectors); + std::vector resolution_vectors; + for( const auto& vector : detectors_configs ) + if( std::count(vector.tags.begin(), vector.tags.end(), observable_config.resolution_vectors_tag) > 0 ) + resolution_vectors.emplace_back(vector.name+"_"+vector.correction_step); + observables.back().SetResolutionVectors(resolution_vectors); + observables.back().Calculate(); + out_file->cd(); + observables.back().Write(); + } + out_file->Close(); +// V1Observables obs_sp( V1Observables::METHODS::MethodOf3SE ); +// obs_sp.SetUvector("u_RESCALED", {"x1", "y1"}); +// obs_sp.SetEPvectors({"W1_RESCALED", "W2_RESCALED", "W3_RESCALED"}, {"x1", "y1"}); +// obs_sp.SetResolutionVectors({ +// "W1_RESCALED", +// "W2_RESCALED", +// "W3_RESCALED", +// "Mf_protons_RESCALED", +// "Mb_protons_RESCALED", +// }); +// obs_sp.SetQqCorrelationsDirectory("/QQ/SP"); +// obs_sp.SetUqCorrelationsDirectory("/uQ/SP"); +// obs_sp.Calculate(); +// V1Observables obs_ep( V1Observables::METHODS::MethodOf3SE ); +// obs_ep.SetUvector("u_RESCALED", {"x1", "y1"}); +// obs_ep.SetEPvectors({"W1_RESCALED", "W2_RESCALED", "W3_RESCALED"}, {"cos1", "sin1"}); +// obs_ep.SetResolutionVectors({ +// "W1_RESCALED", +// "W2_RESCALED", +// "W3_RESCALED", +// "Mf_protons_RESCALED", +// "Mb_protons_RESCALED", +// }); +// obs_ep.SetQqCorrelationsDirectory("/QQ/EP"); +// obs_ep.SetUqCorrelationsDirectory("/uQ/EP"); +// obs_ep.Calculate(); +// V1Observables obs_rs( V1Observables::METHODS::MethodOfRS ); +// obs_rs.SetUvector("u_RESCALED", {"x1", "y1"}); +// obs_rs.SetEPvectors({"R1_RESCALED", "R2_RESCALED"}, {"x1", "y1"}); +// obs_rs.SetResolutionVectors({"R1_RESCALED", "R2_RESCALED"}); +// obs_rs.SetQqCorrelationsDirectory("/QQ/SP"); +// obs_rs.SetUqCorrelationsDirectory("/uQ/SP"); +//// obs_rs.Calculate(); +// auto file_out = TFile::Open("agag-158-szymon.root", "recreate"); +// file_out->mkdir("SP"); +// file_out->cd("/SP"); +// obs_sp.Write(); +// obs_rs.Write(); +// file_out->mkdir("EP"); +// file_out->cd("/EP"); +// obs_ep.Write(); +// file_out->Close(); + + return 0; +} \ No newline at end of file diff --git a/src/QnAnalysisObservables/ConvertConfigs.cpp b/src/QnAnalysisObservables/ConvertConfigs.cpp new file mode 100644 index 00000000..c83c45e1 --- /dev/null +++ b/src/QnAnalysisObservables/ConvertConfigs.cpp @@ -0,0 +1,5 @@ +// +// Created by mikhail on 11/23/20. +// + +#include "ConvertConfigs.hpp" \ No newline at end of file diff --git a/src/QnAnalysisObservables/ConvertConfigs.hpp b/src/QnAnalysisObservables/ConvertConfigs.hpp new file mode 100644 index 00000000..3564a9b8 --- /dev/null +++ b/src/QnAnalysisObservables/ConvertConfigs.hpp @@ -0,0 +1,45 @@ +// +// Created by mikhail on 11/23/20. +// + +#ifndef QNANALYSIS_SRC_QNANALYSISOBSERVABLES_CONVERTCONFIGS_H_ +#define QNANALYSIS_SRC_QNANALYSISOBSERVABLES_CONVERTCONFIGS_H_ + +#include "V1Observables.hpp" +#include +namespace YAML { +template<> +struct convert { + static bool decode(const Node& node, VectorConfig& vector_config) { + if( !node.IsMap() ) { + return false; + } + vector_config.name = node["name"].as(); + vector_config.components = node["components"].as>(); + vector_config.correction_step = node["correction-step"].as(); + try { + vector_config.tags = node["tags"].as>(); + } catch (std::exception&) {} + return true; + } +}; + +template<> +struct convert { + static bool decode(const Node& node, VnObservablesConfig& vn_observables_config) { + auto method_type = node["method"].as(); + if( method_type == "three-sub-event" ) + vn_observables_config.method = V1Observables::METHODS::MethodOf3SE; + else if( method_type == "rnd-sub-event" ) + vn_observables_config.method = V1Observables::METHODS::MethodOfRS; + else return false; + vn_observables_config.qq_correlations_directory = node["QQ-directory"].as(); + vn_observables_config.uq_correlations_directory = node["uQ-directory"].as(); + vn_observables_config.u_vector = node["u-vector"].as(); + vn_observables_config.ep_vectors_tag = node["EP-vectors"].as(); + vn_observables_config.resolution_vectors_tag = node["resolution-vectors"].as(); + return true; + } +}; +}// namespace YAML +#endif//QNANALYSIS_SRC_QNANALYSISOBSERVABLES_CONVERTCONFIGS_H_ diff --git a/src/QnAnalysisObservables/FileManager.cpp b/src/QnAnalysisObservables/FileManager.cpp new file mode 100644 index 00000000..d393550c --- /dev/null +++ b/src/QnAnalysisObservables/FileManager.cpp @@ -0,0 +1,7 @@ +// +// Created by mikhail on 10/23/20. +// + +#include "FileManager.hpp" + +FileManager* FileManager::instance_= nullptr; \ No newline at end of file diff --git a/src/QnAnalysisObservables/FileManager.hpp b/src/QnAnalysisObservables/FileManager.hpp new file mode 100644 index 00000000..2d6fd652 --- /dev/null +++ b/src/QnAnalysisObservables/FileManager.hpp @@ -0,0 +1,46 @@ +// +// Created by mikhail on 10/23/20. +// + +#ifndef OBSERVABLESCALCULATOR_SRC_FILE_MANAGER_H_ +#define OBSERVABLESCALCULATOR_SRC_FILE_MANAGER_H_ + +#include + +class FileManager { +public: + static bool OpenFile( const std::string& name ){ + if( Instance()->file_ ) + Instance()->file_->Close(); + Instance()->file_ = TFile::Open(name.c_str()); + return Instance()->file_ != nullptr; + } + template + static T* GetObject(const std::string& name){ + if( !Instance()->file_ ) + throw std::runtime_error("FileManager::GetObject(): file is not specified"); + T* obj; + Instance()->file_->GetObject(name.c_str(), obj); + if( !obj ) + throw std::runtime_error("FileManager::GetObject(): there is no object named "+name+" in file"); + return obj; + } + static void Cd( const std::string& directory ){ + if( !Instance()->file_ ) + throw std::runtime_error("FileManager::Cd(): file is not specified"); + Instance()->file_->cd(directory.c_str()); + } + +private: + TFile* file_{}; + static FileManager* instance_; + static FileManager* Instance(){ + if( !instance_ ) + instance_ = new FileManager(); + return instance_; + } + FileManager() = default; + virtual ~FileManager() = default; +}; + +#endif // OBSERVABLESCALCULATOR_SRC_FILE_MANAGER_H_ diff --git a/src/QnAnalysisObservables/Method.cpp b/src/QnAnalysisObservables/Method.cpp new file mode 100644 index 00000000..2ac0358d --- /dev/null +++ b/src/QnAnalysisObservables/Method.cpp @@ -0,0 +1,70 @@ +// +// Created by mikhail on 10/25/20. +// + +#include "Method.hpp" + +Qn::DataContainer Method::ReadContainerFromFile( const std::string& directory, const std::vector& vectors ){ + Qn::DataContainerStatCollect* obj; + std::string name; + std::string component; + try{ + name = directory+"/"; + component.clear(); + for (const auto& vector : vectors) { + name += vector.name+"."; + component += vector.component_name; + } + name+=component; + obj = FileManager::GetObject(name); + auto result = Qn::DataContainerStatCalculate(*obj); + result.SetErrors(Qn::StatCalculate::ErrorType::BOOTSTRAP); + return result; + } catch (std::exception&) {} + auto combinations = GetAllCombinations(vectors); + int n=0; + while ( n(name); + auto result = Qn::DataContainerStatCalculate(*obj); + result.SetErrors(Qn::StatCalculate::ErrorType::BOOTSTRAP); + return result; + } catch (std::runtime_error&) { + n++; + } + } + throw std::runtime_error( "Method::ReadContainerFromFile: No such a correlation in file "+name ); +} + +std::vector> Method::GetAllCombinations( std::vector elements ){ + std::vector> result; + if( elements.size() == 1 ){ + result = {elements}; + return result; + } + auto el1 = elements.front(); + std::vector> prev_combinations; + elements.erase(elements.begin()); + prev_combinations=GetAllCombinations(elements); + for( size_t i=0; ij ) + result.back().emplace_back(prev_combinations.at(i).at(k-1)); + } + } + } + return result; +} diff --git a/src/QnAnalysisObservables/Method.hpp b/src/QnAnalysisObservables/Method.hpp new file mode 100644 index 00000000..064027aa --- /dev/null +++ b/src/QnAnalysisObservables/Method.hpp @@ -0,0 +1,68 @@ +// +// Created by mikhail on 10/25/20. +// + +#ifndef OBSERVABLESCALCULATOR_SRC_METHOD_H_ +#define OBSERVABLESCALCULATOR_SRC_METHOD_H_ + +#include +#include + +#include "FileManager.hpp" + +struct VectorComponentConfig { + std::string name; + std::string component_name; +}; + +class Method { +public: + virtual ~Method() = default; + virtual void CalculateObservables() {} + void SetUqDirectory(const std::string &uq_deirectory) { + uq_directory_ = uq_deirectory; + } + void SetQqDirectory(const std::string &qq_directory) { + qq_directory_ = qq_directory; + } + void Write(){ + int i=0; + std::string component = u_vector_config_.component_name; + for( const auto& config : ep_vectors_configs_){ + component+=config.component_name; + } + for( auto resolution: resolutions_){ + resolution.Write( std::data( "res."+ observables_names_.at(i)+"."+component ) ); + ++i; + } + i=0; + for(auto flow : observables_){ + flow.Write( std::data( "obs."+ observables_names_.at(i)+"."+component ) ); + ++i; + } + } +protected: + Method() = default; + Method(VectorComponentConfig u_vector_config, + std::vector q_vector_config, + std::vector resolution_q_vectors_configs) + : u_vector_config_(std::move(u_vector_config)), + ep_vectors_configs_(std::move(q_vector_config)), + resolution_q_vectors_configs_(std::move(resolution_q_vectors_configs)) {} + + static Qn::DataContainer ReadContainerFromFile( const std::string&, const std::vector& vectors ); + static std::vector> GetAllCombinations( std::vector elements ); + + std::string uq_directory_; + std::string qq_directory_; + VectorComponentConfig u_vector_config_; + std::vector ep_vectors_configs_; + std::vector resolution_q_vectors_configs_; + + std::vector observables_names_; + + std::vector> resolutions_; + std::vector> observables_; +}; + +#endif // OBSERVABLESCALCULATOR_SRC_METHOD_H_ diff --git a/src/QnAnalysisObservables/MethodOf3SE.cpp b/src/QnAnalysisObservables/MethodOf3SE.cpp new file mode 100644 index 00000000..69818db5 --- /dev/null +++ b/src/QnAnalysisObservables/MethodOf3SE.cpp @@ -0,0 +1,53 @@ +// +// Created by mikhail on 10/24/20. +// + +#include "MethodOf3SE.hpp" + +Qn::DataContainer MethodOf3SE::CalculateResolution3SE(std::vector correlations){ + auto result = Sqrt(correlations.at(0) * correlations.at(1) / correlations.at(2) ) * sqrt(2.0); + return result; +} + +void MethodOf3SE::CalculateObservables(){ + auto combinations = ConstructResolution3SECombinations(ep_vectors_configs_.front(), resolution_q_vectors_configs_); + std::vector> results; + for( auto res_combination_names : combinations ){ + std::vector> set_for_res_calc; + set_for_res_calc.push_back( ReadContainerFromFile( qq_directory_, res_combination_names.at(0) ) ); + set_for_res_calc.push_back( ReadContainerFromFile( qq_directory_, res_combination_names.at(1) ) ); + set_for_res_calc.push_back( ReadContainerFromFile( qq_directory_, res_combination_names.at(2) ) ); + observables_names_.push_back( u_vector_config_.name+"."+ ep_vectors_configs_.front().name+"("+res_combination_names.at(2).front().name+","+res_combination_names.at(2).back().name+")" ); + resolutions_.push_back( CalculateResolution3SE( set_for_res_calc ) ); + } + auto uq_correlation_confg = ep_vectors_configs_; + uq_correlation_confg.insert(uq_correlation_confg.begin(), u_vector_config_); + auto uq_correlation = ReadContainerFromFile( uq_directory_, uq_correlation_confg ); + for( const auto& res : resolutions_ ){ + observables_.push_back(uq_correlation/res*2.0); + } +} + +std::vector>> MethodOf3SE::ConstructResolution3SECombinations( + VectorComponentConfig reff_q, std::vector q_combination){ + auto q1 = std::move(reff_q); + std::vector>> combinations; + for (size_t i=0; i> combination; + combination.push_back({q1, q2}); + combination.push_back({q1, q3}); + combination.push_back({q2, q3}); + combinations.push_back(combination); + } + } + return combinations; +} diff --git a/src/QnAnalysisObservables/MethodOf3SE.hpp b/src/QnAnalysisObservables/MethodOf3SE.hpp new file mode 100644 index 00000000..273a6de5 --- /dev/null +++ b/src/QnAnalysisObservables/MethodOf3SE.hpp @@ -0,0 +1,30 @@ +// +// Created by mikhail on 10/24/20. +// + +#ifndef OBSERVABLESCALCULATOR_SRC_METHOD_OF_3SE_H_ +#define OBSERVABLESCALCULATOR_SRC_METHOD_OF_3SE_H_ + +#include + +#include +#include + +#include "FileManager.hpp" +#include "Method.hpp" + +class MethodOf3SE : public Method { +public: + MethodOf3SE(const VectorComponentConfig&u_vector_config, + const std::vector &q_vector_config, + const std::vector &resolution_q_vectors_configs) + : Method(u_vector_config, q_vector_config, resolution_q_vectors_configs) {} + ~MethodOf3SE() override = default; + void CalculateObservables() override; + +protected: + static Qn::DataContainer CalculateResolution3SE(std::vector correlations); + static std::vector>> ConstructResolution3SECombinations(VectorComponentConfig reff_q, std::vector q_combination); +}; + +#endif // OBSERVABLESCALCULATOR_SRC_METHOD_OF_3SE_H_ diff --git a/src/QnAnalysisObservables/MethodOfRS.cpp b/src/QnAnalysisObservables/MethodOfRS.cpp new file mode 100644 index 00000000..6da9fcbd --- /dev/null +++ b/src/QnAnalysisObservables/MethodOfRS.cpp @@ -0,0 +1,22 @@ +// +// Created by mikhail on 11/21/20. +// + +#include "MethodOfRS.hpp" +MethodOfRS::MethodOfRS(const VectorComponentConfig& u_vector_config, const std::vector& q_vector_config, const std::vector& resolution_q_vectors_configs) : Method(u_vector_config, q_vector_config, resolution_q_vectors_configs) {} +MethodOfRS::~MethodOfRS() = default; +void MethodOfRS::CalculateObservables() { + if( std::size(resolution_q_vectors_configs_) != 2 ) + throw std::runtime_error( "MethodOfRS::CalculateObservables(): " + "only 2 vectors must be in set for resolution correction calculation" ); + for( auto& res_vec : resolution_q_vectors_configs_ ){ + res_vec.component_name = ep_vectors_configs_.front().component_name; + } + auto qq_container = ReadContainerFromFile(qq_directory_, resolution_q_vectors_configs_); + for( const auto& ep_vector : ep_vectors_configs_ ){ + auto uq_correlation = ReadContainerFromFile( uq_directory_, {u_vector_config_, ep_vector} ); + resolutions_.emplace_back(Sqrt( qq_container )); + observables_.emplace_back( uq_correlation*sqrt(2.0) / resolutions_.back() ); + observables_names_.emplace_back( u_vector_config_.name+"."+ ep_vector.name+"("+resolution_q_vectors_configs_.at(0).name+","+resolution_q_vectors_configs_.at(1).name+")" ); + } +} diff --git a/src/QnAnalysisObservables/MethodOfRS.hpp b/src/QnAnalysisObservables/MethodOfRS.hpp new file mode 100644 index 00000000..809b7939 --- /dev/null +++ b/src/QnAnalysisObservables/MethodOfRS.hpp @@ -0,0 +1,17 @@ +// +// Created by mikhail on 11/21/20. +// + +#ifndef QNANALYSIS_SRC_QNANALYSISOBSERVABLES_METHOD_OF_RS_H_ +#define QNANALYSIS_SRC_QNANALYSISOBSERVABLES_METHOD_OF_RS_H_ + +#include "Method.hpp" + +class MethodOfRS : public Method { + public: + MethodOfRS(const VectorComponentConfig& u_vector_config, const std::vector& q_vector_config, const std::vector& resolution_q_vectors_configs); + ~MethodOfRS() override; + void CalculateObservables() override; +}; + +#endif//QNANALYSIS_SRC_QNANALYSISOBSERVABLES_METHOD_OF_RS_H_ diff --git a/src/QnAnalysisObservables/V1Observables.cpp b/src/QnAnalysisObservables/V1Observables.cpp new file mode 100644 index 00000000..4e56b46c --- /dev/null +++ b/src/QnAnalysisObservables/V1Observables.cpp @@ -0,0 +1,5 @@ +// +// Created by mikhail on 10/25/20. +// + +#include "V1Observables.hpp" diff --git a/src/QnAnalysisObservables/V1Observables.hpp b/src/QnAnalysisObservables/V1Observables.hpp new file mode 100644 index 00000000..becfa00e --- /dev/null +++ b/src/QnAnalysisObservables/V1Observables.hpp @@ -0,0 +1,107 @@ +// +// Created by mikhail on 10/25/20. +// + +#ifndef OBSERVABLESCALCULATOR_SRC_CALCULATOR_H_ +#define OBSERVABLESCALCULATOR_SRC_CALCULATOR_H_ + +#include +#include + +#include "Method.hpp" +#include "MethodOf3SE.hpp" +#include "MethodOfRS.hpp" + +struct VectorConfig{ + std::string name; + std::string correction_step; + std::vector components; + std::vector tags; +}; + +class V1Observables { + public: + enum class METHODS{ + MethodOf3SE, + MethodOfRS + }; + explicit V1Observables(METHODS method_type) : method_type_(method_type) {} + virtual ~V1Observables() = default; + void SetUvector( std::string name, std::vector components ){ + u_vector_name_ = std::move(name); + u_vector_components_ = std::move(components); + } + void SetUvector( const VectorConfig& u_vector ){ + u_vector_name_ = u_vector.name+"_"+u_vector.correction_step; + u_vector_components_ = u_vector.components; + }; + void SetEPvectors( std::vector names, std::vector components ){ + ep_vectors_names_ = std::move(names); + ep_vectors_components_ = std::move(components); + } + void SetEPvectors( const std::vector& vectors){ + for( const auto& vector : vectors ) { + ep_vectors_names_.emplace_back(vector.name + "_" + vector.correction_step); + } + ep_vectors_components_ = vectors.front().components; + } + void SetResolutionVectors( std::vector names){ + resolution_q_vectors_ = std::move(names); + } + void SetUqCorrelationsDirectory(const std::string& uq_correlations_directory) { + uq_correlations_directory_ = uq_correlations_directory; + } + void SetQqCorrelationsDirectory(const std::string& qq_correlations_directory) { + qq_correlations_directory_ = qq_correlations_directory; + } + void Calculate(){ + std::vector resolution_qs; + for( const auto& name : resolution_q_vectors_ ){ + resolution_qs.push_back({name}); + } + for( const auto& component : u_vector_components_ ) + for( const auto& ep_vector : ep_vectors_names_ ) + for (const auto& ep_component : ep_vectors_components_) { + switch (method_type_) { + case METHODS::MethodOf3SE: + methods_.emplace_back(new MethodOf3SE( + {u_vector_name_, component}, {{ep_vector, ep_component}}, resolution_qs)); + break; + case METHODS::MethodOfRS: + methods_.emplace_back(new MethodOfRS( + {u_vector_name_, component}, {{ep_vector, ep_component}}, resolution_qs)); + break; + } + methods_.back()->SetQqDirectory(qq_correlations_directory_); + methods_.back()->SetUqDirectory(uq_correlations_directory_); + } + for( auto method : methods_ ) + method->CalculateObservables(); + } + void Write(){ + for( auto method : methods_ ) + method->Write(); + } + + private: + std::vector methods_; + METHODS method_type_; + std::string uq_correlations_directory_; + std::string qq_correlations_directory_; + std::string u_vector_name_; + std::vector u_vector_components_; + std::vector ep_vectors_names_; + std::vector ep_vectors_components_; + std::vector resolution_q_vectors_; +}; + +struct VnObservablesConfig { + V1Observables::METHODS method; + std::string uq_correlations_directory; + std::string qq_correlations_directory; + VectorConfig u_vector; + std::string ep_vectors_tag; + std::string resolution_vectors_tag; +}; + +#endif // OBSERVABLESCALCULATOR_SRC_CALCULATOR_H_