Added macros/configs for QnAnalysis

qnanalysis_macros/cuts.C

@@ -0,0 +1,14 @@
+#include "CutsRegistry.hpp"
+
+int cuts() {
+    using AnalysisTree::RegisterCuts;
+    using namespace AnalysisTree::Version1;
+
+    SimpleCut b_mc(Variable("Event","B"), 0., 14.);
+
+    RegisterCuts("models", Cuts("McEvent", {
+                b_mc
+    }));
+
+    return 0;
+}

qnanalysis_macros/models-analysis-config.yml

@@ -0,0 +1,310 @@
+# --- Set up axes for differential u-vectors (for example, for vn(pT, y))
+_axes:
+  - &axis_pT
+    name: Particles/pT
+    nb: 50
+    lo: 0.
+    hi: 5.
+  - &axis_rapidity
+    name: Particles/rapidity
+    nb: 100
+    lo: -5.
+    hi: 5.
+  - &axis_eta
+    name: Particles/eta
+    nb: 100
+    lo: -5.
+    hi: 5.
+# --- Main configuration
+qn_all:
+  # --- Set up eventwise varibles (for example, centrality)
+  event-variables:
+    - Event/B
+  axes:
+    - { name: Event/B, nb: 140, lo: 0., hi: 14. }
+  # --- Configure Q-vectors
+  q-vectors:
+  # --- u(centrality, pT, y) for all charged hadrons in TPC L
+  - name: u_hadrons_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+  # --- u(centrality, pT, y) for all charged hadrons in TPC R
+  - name: u_hadrons_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+  # --- u(centrality, pT, y) for protons in TPC L
+  - name: u_pr_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particle/pid: { equals: 2212 }
+  # --- u(centrality, pT, y) for protons in TPC R
+  - name: u_pr_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particle/pid: { equals: 2212 }
+  # --- u(centrality, pT, y) for pi+ in TPC L
+  - name: u_pip_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particle/pid: { equals: 211 }
+  # --- u(centrality, pT, y) for pi+ in TPC R
+  - name: u_pip_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particle/pid: { equals: 211 }
+  # --- u(centrality, pT, y) for pi- in TPC L
+  - name: u_pim_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particle/pid: { equals: -211 }
+  # --- u(centrality, pT, y) for pi- in TPC R
+  - name: u_pim_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particle/pid: { equals: -211 }
+  # --- u(centrality, pT, y) for K+ in TPC L
+  - name: u_kap_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particle/pid: { equals: 321 }
+  # --- u(centrality, pT, y) for K+ in TPC R
+  - name: u_kap_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particle/pid: { equals: 321 }
+  # --- u(centrality, pT, y) for K- in TPC L
+  - name: u_kam_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particle/pid: { equals: -321 }
+  # --- u(centrality, pT, y) for K- in TPC R
+  - name: u_kam_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    axes:
+    - *axis_pT
+    - *axis_rapidity
+    - *axis_eta
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particle/pid: { equals: -321 }
+  # --- Q(centrality) for all charged hadrons in TPC L
+  - name: Q_hadrons_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/is_charged: { equals: 1 }
+  # --- Q(centrality) for all charged hadrons in TPC R
+  - name: Q_hadrons_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/is_charged: { equals: 1 }
+  # --- Q(centrality) for protons in TPC L
+  - name: Q_pr_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/pid: { equals: 2212 }
+  # --- Q(centrality) for protons in TPC R
+  - name: Q_pr_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/pid: { equals: 2212 }
+  # --- Q(centrality) for pi+ in TPC L
+  - name: Q_pip_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/pid: { equals: 211 }
+  # --- Q(centrality) for pi+ in TPC R
+  - name: Q_pip_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/pid: { equals: 211 }
+  # --- Q(centrality) for pi- in TPC L
+  - name: Q_pim_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/pid: { equals: -211 }
+  # --- Q(centrality) for pi- in TPC R
+  - name: Q_pim_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/pid: { equals: -211 }
+  # --- Q(centrality) for K+ in TPC L
+  - name: Q_kap_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/pid: { equals: 321 }
+  # --- Q(centrality) for K+ in TPC R
+  - name: Q_kap_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/pid: { equals: 321 }
+  # --- Q(centrality) for K- in TPC L
+  - name: Q_kam_TPC_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-1.5, -0.05] }
+      Particles/pid: { equals: -321 }
+  # --- Q(centrality) for K- in TPC R
+  - name: Q_kam_TPC_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [0.05, 1.5] }
+      Particles/pid: { equals: -321 }
+  # --- Q(centrality) for all charged hadrons in FHCal L
+  - name: Q_hadrons_FHCal_L
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [-5., -2.] }
+  # --- Q(centrality) for all charged hadrons in FHCal R
+  - name: Q_hadrons_FHCal_R
+    type: track
+    phi:  Particles/phi
+    weight: Ones
+    norm: m
+    cuts:
+      Particles/eta: { range: [2., 5.] }
+  # --- Q(centrality) vector based on reaction plane
+  - name: Q_RP
+    type: psi
+    phi: Event/PhiRp
+    weight: Ones
+    norm: m