diff --git a/PWGUD/Tasks/SGSpectraAnalyzer.cxx b/PWGUD/Tasks/SGSpectraAnalyzer.cxx index 29d493f9af5..dbdf8596493 100644 --- a/PWGUD/Tasks/SGSpectraAnalyzer.cxx +++ b/PWGUD/Tasks/SGSpectraAnalyzer.cxx @@ -12,254 +12,255 @@ using namespace o2; using namespace o2::aod; using namespace o2::framework; using namespace o2::framework::expressions; -#define mpion 0.1396 +#define mpion 0.1396 #define mkaon 0.4937 #define mproton 0.9383 -struct SGSpectraAnalyzer{ +struct SGSpectraAnalyzer { HistogramRegistry registry{ "registry", - { - // Pion histograms for each eta bin and gapSide - {"hPtPion_gap0", "Pion pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_gap1", "Pion pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_gap2", "Pion pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hEtaPion_gap0", "Pion eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaPion_gap1", "Pion eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaPion_gap2", "Pion eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hPtPion_etaBin1_gap0", "Pion pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin2_gap0", "Pion pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin3_gap0", "Pion pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin4_gap0", "Pion pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin1_gap1", "Pion pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin2_gap1", "Pion pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin3_gap1", "Pion pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin4_gap1", "Pion pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin1_gap2", "Pion pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin2_gap2", "Pion pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin3_gap2", "Pion pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtPion_etaBin4_gap2", "Pion pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {// Pion histograms for each eta bin and gapSide + {"hPtPion_gap0", "Pion pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_gap1", "Pion pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_gap2", "Pion pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hEtaPion_gap0", "Pion eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaPion_gap1", "Pion eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaPion_gap2", "Pion eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hPtPion_etaBin1_gap0", "Pion pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin2_gap0", "Pion pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin3_gap0", "Pion pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin4_gap0", "Pion pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin1_gap1", "Pion pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin2_gap1", "Pion pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin3_gap1", "Pion pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin4_gap1", "Pion pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin1_gap2", "Pion pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin2_gap2", "Pion pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin3_gap2", "Pion pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtPion_etaBin4_gap2", "Pion pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - // Kaon histograms for each eta bin and gapSide - {"hPtKaon_gap0", "Kaon pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_gap1", "Kaon pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_gap2", "Kaon pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hEtaKaon_gap0", "Kaon eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaKaon_gap1", "Kaon eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaKaon_gap2", "Kaon eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hPtKaon_etaBin1_gap0", "Kaon pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin2_gap0", "Kaon pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin3_gap0", "Kaon pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin4_gap0", "Kaon pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin1_gap1", "Kaon pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin2_gap1", "Kaon pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin3_gap1", "Kaon pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin4_gap1", "Kaon pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin1_gap2", "Kaon pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin2_gap2", "Kaon pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin3_gap2", "Kaon pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtKaon_etaBin4_gap2", "Kaon pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + // Kaon histograms for each eta bin and gapSide + {"hPtKaon_gap0", "Kaon pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_gap1", "Kaon pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_gap2", "Kaon pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hEtaKaon_gap0", "Kaon eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaKaon_gap1", "Kaon eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaKaon_gap2", "Kaon eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hPtKaon_etaBin1_gap0", "Kaon pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin2_gap0", "Kaon pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin3_gap0", "Kaon pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin4_gap0", "Kaon pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin1_gap1", "Kaon pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin2_gap1", "Kaon pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin3_gap1", "Kaon pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin4_gap1", "Kaon pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin1_gap2", "Kaon pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin2_gap2", "Kaon pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin3_gap2", "Kaon pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtKaon_etaBin4_gap2", "Kaon pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - // Proton histograms for each eta bin and gapSide - {"hPtProton_gap0", "Proton pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_gap1", "Proton pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_gap2", "Proton pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hEtaProton_gap0", "Proton eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaProton_gap1", "Proton eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hEtaProton_gap2", "Proton eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, - {"hPtProton_etaBin1_gap0", "Proton pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin2_gap0", "Proton pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin3_gap0", "Proton pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin4_gap0", "Proton pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin1_gap1", "Proton pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin2_gap1", "Proton pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin3_gap1", "Proton pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin4_gap1", "Proton pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin1_gap2", "Proton pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin2_gap2", "Proton pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin3_gap2", "Proton pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, - {"hPtProton_etaBin4_gap2", "Proton pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}} - } - }; -using udtracks = soa::Join; -//using udtracks = soa::Join; -using udtracksfull = soa::Join; -//using UDCollisionsFull = soa::Join; -using UDCollisionsFull = soa::Join; // -//using UDCollisionsFull = soa::Join; -using UDCollisionFull = UDCollisionsFull::iterator; + // Proton histograms for each eta bin and gapSide + {"hPtProton_gap0", "Proton pT in eta [-1, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_gap1", "Proton pT in eta [-1, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_gap2", "Proton pT in eta [-1, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hEtaProton_gap0", "Proton eta Gap 0; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaProton_gap1", "Proton eta Gap 1; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hEtaProton_gap2", "Proton eta Gap 2; eta; Entries", {HistType::kTH1F, {{100, -1, 1}}}}, + {"hPtProton_etaBin1_gap0", "Proton pT in eta [-1, -0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin2_gap0", "Proton pT in eta [-0.5, 0] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin3_gap0", "Proton pT in eta [0, 0.5] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin4_gap0", "Proton pT in eta [0.5, 1] Gap 0; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin1_gap1", "Proton pT in eta [-1, -0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin2_gap1", "Proton pT in eta [-0.5, 0] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin3_gap1", "Proton pT in eta [0, 0.5] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin4_gap1", "Proton pT in eta [0.5, 1] Gap 1; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin1_gap2", "Proton pT in eta [-1, -0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin2_gap2", "Proton pT in eta [-0.5, 0] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin3_gap2", "Proton pT in eta [0, 0.5] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}, + {"hPtProton_etaBin4_gap2", "Proton pT in eta [0.5, 1] Gap 2; pT (GeV/c); Entries", {HistType::kTH1F, {{100, 0, 10}}}}}}; + using udtracks = soa::Join; + // using udtracks = soa::Join; + using udtracksfull = soa::Join; + // using UDCollisionsFull = soa::Join; + using UDCollisionsFull = soa::Join; // + // using UDCollisionsFull = soa::Join; + using UDCollisionFull = UDCollisionsFull::iterator; -void process(UDCollisionFull const& collision, udtracksfull const& tracks){ + void process(UDCollisionFull const& collision, udtracksfull const& tracks) + { TLorentzVector a; int gapSide = collision.gapSide(); - if(gapSide < 0 || gapSide > 2) return; + if (gapSide < 0 || gapSide > 2) + return; for (auto& track : tracks) { - if(!track.isPVContributor())break; + if (!track.isPVContributor()) + break; float nSigmaPi = track.tpcNSigmaPi(); float nSigmaK = track.tpcNSigmaKa(); float nSigmaP = track.tpcNSigmaPr(); if (std::abs(nSigmaPi) < 3.0) { - a.SetXYZM(track.px(), track.py(), track.pz(),mpion); - fillHistograms("Pion", a.Pt(), a.Eta(), gapSide); + a.SetXYZM(track.px(), track.py(), track.pz(), mpion); + fillHistograms("Pion", a.Pt(), a.Eta(), gapSide); } if (std::abs(nSigmaK) < 3.0) { - a.SetXYZM(track.px(), track.py(), track.pz(),mkaon); + a.SetXYZM(track.px(), track.py(), track.pz(), mkaon); fillHistograms("Kaon", a.Pt(), a.Eta(), gapSide); } if (std::abs(nSigmaP) < 3.0) { - a.SetXYZM(track.px(), track.py(), track.pz(),mproton); + a.SetXYZM(track.px(), track.py(), track.pz(), mproton); fillHistograms("Proton", a.Pt(), a.Eta(), gapSide); } } } - void fillHistograms(const std::string& particleType, float pt, float eta, int gapSide) { + void fillHistograms(const std::string& particleType, float pt, float eta, int gapSide) + { if (particleType == "Pion") { - if (gapSide == 0) { - registry.get(HIST("hPtPion_gap0"))->Fill(pt); - registry.get(HIST("hEtaPion_gap0"))->Fill(eta); - } else if (gapSide == 1) { - registry.get(HIST("hPtPion_gap1"))->Fill(pt); - registry.get(HIST("hEtaPion_gap1"))->Fill(eta); - } else if (gapSide == 2) { - registry.get(HIST("hPtPion_gap2"))->Fill(pt); - registry.get(HIST("hEtaPion_gap2"))->Fill(eta); - } + if (gapSide == 0) { + registry.get(HIST("hPtPion_gap0"))->Fill(pt); + registry.get(HIST("hEtaPion_gap0"))->Fill(eta); + } else if (gapSide == 1) { + registry.get(HIST("hPtPion_gap1"))->Fill(pt); + registry.get(HIST("hEtaPion_gap1"))->Fill(eta); + } else if (gapSide == 2) { + registry.get(HIST("hPtPion_gap2"))->Fill(pt); + registry.get(HIST("hEtaPion_gap2"))->Fill(eta); + } } else if (particleType == "Kaon") { - if (gapSide == 0) { - registry.get(HIST("hPtKaon_gap0"))->Fill(pt); - registry.get(HIST("hEtaKaon_gap0"))->Fill(eta); - } else if (gapSide == 1) { - registry.get(HIST("hPtKaon_gap1"))->Fill(pt); - registry.get(HIST("hEtaKaon_gap1"))->Fill(eta); - } else if (gapSide == 2) { - registry.get(HIST("hPtKaon_gap2"))->Fill(pt); - registry.get(HIST("hEtaKaon_gap2"))->Fill(eta); - } + if (gapSide == 0) { + registry.get(HIST("hPtKaon_gap0"))->Fill(pt); + registry.get(HIST("hEtaKaon_gap0"))->Fill(eta); + } else if (gapSide == 1) { + registry.get(HIST("hPtKaon_gap1"))->Fill(pt); + registry.get(HIST("hEtaKaon_gap1"))->Fill(eta); + } else if (gapSide == 2) { + registry.get(HIST("hPtKaon_gap2"))->Fill(pt); + registry.get(HIST("hEtaKaon_gap2"))->Fill(eta); + } } else if (particleType == "Proton") { - if (gapSide == 0) { - registry.get(HIST("hPtProton_gap0"))->Fill(pt); - registry.get(HIST("hEtaProton_gap0"))->Fill(eta); - } else if (gapSide == 1) { - registry.get(HIST("hPtProton_gap1"))->Fill(pt); - registry.get(HIST("hEtaProton_gap1"))->Fill(eta); - } else if (gapSide == 2) { - registry.get(HIST("hPtProton_gap2"))->Fill(pt); - registry.get(HIST("hEtaProton_gap2"))->Fill(eta); - } + if (gapSide == 0) { + registry.get(HIST("hPtProton_gap0"))->Fill(pt); + registry.get(HIST("hEtaProton_gap0"))->Fill(eta); + } else if (gapSide == 1) { + registry.get(HIST("hPtProton_gap1"))->Fill(pt); + registry.get(HIST("hEtaProton_gap1"))->Fill(eta); + } else if (gapSide == 2) { + registry.get(HIST("hPtProton_gap2"))->Fill(pt); + registry.get(HIST("hEtaProton_gap2"))->Fill(eta); + } } std::vector> etaBins = {{-1, -0.5}, {-0.5, 0}, {0, 0.5}, {0.5, 1}}; for (int i = 0; i < 4; ++i) { if (eta > etaBins[i].first && eta <= etaBins[i].second) { - if (i==0){ - if (particleType == "Pion") { - if (gapSide == 0) { - registry.get(HIST("hPtPion_etaBin1_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtPion_etaBin1_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtPion_etaBin1_gap2"))->Fill(pt); + if (i == 0) { + if (particleType == "Pion") { + if (gapSide == 0) { + registry.get(HIST("hPtPion_etaBin1_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtPion_etaBin1_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtPion_etaBin1_gap2"))->Fill(pt); + } + } else if (particleType == "Kaon") { + if (gapSide == 0) { + registry.get(HIST("hPtKaon_etaBin1_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtKaon_etaBin1_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtKaon_etaBin1_gap2"))->Fill(pt); + } + } else if (particleType == "Proton") { + if (gapSide == 0) { + registry.get(HIST("hPtProton_etaBin1_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtProton_etaBin1_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtProton_etaBin1_gap2"))->Fill(pt); + } + } } - } else if (particleType == "Kaon") { - if (gapSide == 0) { - registry.get(HIST("hPtKaon_etaBin1_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtKaon_etaBin1_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtKaon_etaBin1_gap2"))->Fill(pt); - } - } else if (particleType == "Proton") { - if (gapSide == 0) { - registry.get(HIST("hPtProton_etaBin1_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtProton_etaBin1_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtProton_etaBin1_gap2"))->Fill(pt); + if (i == 1) { + if (particleType == "Pion") { + if (gapSide == 0) { + registry.get(HIST("hPtPion_etaBin2_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtPion_etaBin2_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtPion_etaBin2_gap2"))->Fill(pt); + } + } else if (particleType == "Kaon") { + if (gapSide == 0) { + registry.get(HIST("hPtKaon_etaBin2_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtKaon_etaBin2_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtKaon_etaBin2_gap2"))->Fill(pt); + } + } else if (particleType == "Proton") { + if (gapSide == 0) { + registry.get(HIST("hPtProton_etaBin2_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtProton_etaBin2_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtProton_etaBin2_gap2"))->Fill(pt); + } + } } - } - } - if (i==1){ - if (particleType == "Pion") { - if (gapSide == 0) { - registry.get(HIST("hPtPion_etaBin2_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtPion_etaBin2_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtPion_etaBin2_gap2"))->Fill(pt); + if (i == 2) { + if (particleType == "Pion") { + if (gapSide == 0) { + registry.get(HIST("hPtPion_etaBin3_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtPion_etaBin3_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtPion_etaBin3_gap2"))->Fill(pt); + } + } else if (particleType == "Kaon") { + if (gapSide == 0) { + registry.get(HIST("hPtKaon_etaBin3_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtKaon_etaBin3_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtKaon_etaBin3_gap2"))->Fill(pt); + } + } else if (particleType == "Proton") { + if (gapSide == 0) { + registry.get(HIST("hPtProton_etaBin3_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtProton_etaBin3_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtProton_etaBin3_gap2"))->Fill(pt); + } + } } - } else if (particleType == "Kaon") { - if (gapSide == 0) { - registry.get(HIST("hPtKaon_etaBin2_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtKaon_etaBin2_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtKaon_etaBin2_gap2"))->Fill(pt); - } - } else if (particleType == "Proton") { - if (gapSide == 0) { - registry.get(HIST("hPtProton_etaBin2_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtProton_etaBin2_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtProton_etaBin2_gap2"))->Fill(pt); + if (i == 3) { + if (particleType == "Pion") { + if (gapSide == 0) { + registry.get(HIST("hPtPion_etaBin4_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtPion_etaBin4_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtPion_etaBin4_gap2"))->Fill(pt); + } + } else if (particleType == "Kaon") { + if (gapSide == 0) { + registry.get(HIST("hPtKaon_etaBin4_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtKaon_etaBin4_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtKaon_etaBin4_gap2"))->Fill(pt); + } + } else if (particleType == "Proton") { + if (gapSide == 0) { + registry.get(HIST("hPtProton_etaBin4_gap0"))->Fill(pt); + } else if (gapSide == 1) { + registry.get(HIST("hPtProton_etaBin4_gap1"))->Fill(pt); + } else if (gapSide == 2) { + registry.get(HIST("hPtProton_etaBin4_gap2"))->Fill(pt); + } + } } - } - } - if (i==2){ - if (particleType == "Pion") { - if (gapSide == 0) { - registry.get(HIST("hPtPion_etaBin3_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtPion_etaBin3_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtPion_etaBin3_gap2"))->Fill(pt); - } - } else if (particleType == "Kaon") { - if (gapSide == 0) { - registry.get(HIST("hPtKaon_etaBin3_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtKaon_etaBin3_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtKaon_etaBin3_gap2"))->Fill(pt); - } - } else if (particleType == "Proton") { - if (gapSide == 0) { - registry.get(HIST("hPtProton_etaBin3_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtProton_etaBin3_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtProton_etaBin3_gap2"))->Fill(pt); - } - } - } - if (i==3){ - if (particleType == "Pion") { - if (gapSide == 0) { - registry.get(HIST("hPtPion_etaBin4_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtPion_etaBin4_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtPion_etaBin4_gap2"))->Fill(pt); - } - } else if (particleType == "Kaon") { - if (gapSide == 0) { - registry.get(HIST("hPtKaon_etaBin4_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtKaon_etaBin4_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtKaon_etaBin4_gap2"))->Fill(pt); - } - } else if (particleType == "Proton") { - if (gapSide == 0) { - registry.get(HIST("hPtProton_etaBin4_gap0"))->Fill(pt); - } else if (gapSide == 1) { - registry.get(HIST("hPtProton_etaBin4_gap1"))->Fill(pt); - } else if (gapSide == 2) { - registry.get(HIST("hPtProton_etaBin4_gap2"))->Fill(pt); - } - } - } } } } diff --git a/PWGUD/Tasks/SG_FIT_Analyzer.cxx b/PWGUD/Tasks/SG_FIT_Analyzer.cxx index d822ffb4d15..520beafc7ae 100644 --- a/PWGUD/Tasks/SG_FIT_Analyzer.cxx +++ b/PWGUD/Tasks/SG_FIT_Analyzer.cxx @@ -172,89 +172,89 @@ struct SG_FIT_Analyzer { // UDTutorial01 registry.add("FIT/TFT0C", "Time FT0C", {HistType::kTH1F, {{100, -49.5, 50.5}}}); registry.add("FIT/TFDDA", "Time FDDA", {HistType::kTH1F, {{100, -49.5, 50.5}}}); registry.add("FIT/TFDDC", "Time FDDC", {HistType::kTH1F, {{100, -49.5, 50.5}}}); - registry.add("ZDC/AZNAC", "ZNA vs ZNC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{500,-5.5,199.5}}}); - registry.add("ZDC/CZNAC", "ZNA vs ZNC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{500,-5.5,199.5}}}); - registry.add("ZDC/ACZNAC", "ZNA vs ZNC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{500,-5.5,199.5}}}); - registry.add("ZDC/AZNAFT0A", "ZNA vs FT0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/AZNAFV0A", "ZNA vs FV0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/AZNCFV0A", "ZNC vs FV0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNAFT0A", "ZNA vs FT0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/CZNAFV0A", "ZNA vs FV0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNCFV0A", "ZNC vs FV0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNAFT0A", "ZNA vs FT0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/ACZNAFV0A", "ZNA vs FV0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNCFV0A", "ZNC vs FV0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/AZNAFT0C", "ZNA vs FT0C, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/CZNAFT0C", "ZNA vs FT0C, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/ACZNAFT0C", "ZNA vs FT0C, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/AZNCFT0A", "ZNC vs FT0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/CZNCFT0A", "ZNC vs FT0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/ACZNCFT0A", "ZNC vs FT0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,399.5}}}); - registry.add("ZDC/AZNCFT0C", "ZNC vs FT0C, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/CZNCFT0C", "ZNC vs FT0C, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/ACZNCFT0C", "ZNC vs FT0C, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,299.5}}}); - registry.add("ZDC/AZNAFDDA", "ZNA vs FDDA, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNAFDDA", "ZNA vs FDDA, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNAFDDA", "ZNA vs FDDA, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/AZNAFDDC", "ZNA vs FDDC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNAFDDC", "ZNA vs FDDC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNAFDDC", "ZNA vs FDDC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/AZNCFDDA", "ZNC vs FDDA, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNCFDDA", "ZNC vs FDDA, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNCFDDA", "ZNC vs FDDA, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/AZNCFDDC", "ZNC vs FDDC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/CZNCFDDC", "ZNC vs FDDC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("ZDC/ACZNCFDDC", "ZNC vs FDDC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5},{100,-.5,199.5}}}); - registry.add("FIT/TAFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TAFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TTFT0", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{40,-5.5,34.5}}}); - registry.add("FIT/TTFT01", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{40,-5.5,34.5}}}); - registry.add("FIT/TTFT02", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{40,-5.5,34.5}}}); - registry.add("FIT/ACFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5},{100,-.5,99.5}}}); - registry.add("FIT/AFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5},{100,-.5,99.5}}}); - registry.add("FIT/CFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5},{100,-.5,99.5}}}); - registry.add("FIT/AFT0AC", "Amp FT0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/CFT0AC", "Amp FT0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/ACFT0AC", "Amp FT0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/AFTV0A", "Amp FT0A vs Amp FV0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/CFTV0A", "Amp FT0A vs Amp FV0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/ACFTV0A", "Amp FT0A vs Amp FV0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/AFTV0C", "Amp FV0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/CFTV0C", "Amp FV0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/ACFTV0C", "Amp FV0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/AFDV0A", "Amp FDDA vs Amp FV0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/CFDV0A", "Amp FDDA vs Amp FV0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/ACFDV0A", "Amp FDDA vs Amp FV0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/AFDT0A", "Amp FDDA vs Amp FT0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/CFDT0A", "Amp FDDA vs Amp FT0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/ACFDT0A", "Amp FDDA vs Amp FT0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5},{100,-.5,499.5}}}); - registry.add("FIT/TAFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TAFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TAFDDC", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TCFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TCFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TCFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TCFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/TCFDDC", "Time vs Amp FDDC", {HistType::kTH2F, {{40, -5.5, 34.5},{50,-.5,149.5}}}); - registry.add("FIT/MFV0A", "Track number vs Amp FV0A", {HistType::kTH2F, {{50, -.5, 49.5},{250,-.5,499.5}}}); - registry.add("FIT/MFT0A", "Track number vs Amp FT0A", {HistType::kTH2F, {{50, -.5, 49.5},{250,-.5,499.5}}}); - registry.add("FIT/MFT0C", "Track number vs Amp FT0C", {HistType::kTH2F, {{50, -.5, 49.5},{250,-.5,399.5}}}); - registry.add("FIT/MFDDA", "Track number vs Amp FDDA", {HistType::kTH2F, {{50, -.5, 49.5},{250,-.5,499.5}}}); - registry.add("FIT/MFDDC", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5},{250,-.5,499.5}}}); - registry.add("ZDC/MAZNA", "Track number vs Amp FV0A", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); - registry.add("ZDC/MAZNC", "Track number vs Amp FT0A", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); - registry.add("ZDC/MCZNC", "Track number vs Amp FT0C", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); - registry.add("ZDC/MCZNA", "Track number vs Amp FDDA", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); - registry.add("ZDC/MACZNA", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); - registry.add("ZDC/MACZNC", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5},{250,-5.5,149.5}}}); + registry.add("ZDC/AZNAC", "ZNA vs ZNC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {500, -5.5, 199.5}}}); + registry.add("ZDC/CZNAC", "ZNA vs ZNC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {500, -5.5, 199.5}}}); + registry.add("ZDC/ACZNAC", "ZNA vs ZNC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {500, -5.5, 199.5}}}); + registry.add("ZDC/AZNAFT0A", "ZNA vs FT0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/AZNAFV0A", "ZNA vs FV0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/AZNCFV0A", "ZNC vs FV0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNAFT0A", "ZNA vs FT0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/CZNAFV0A", "ZNA vs FV0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNCFV0A", "ZNC vs FV0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNAFT0A", "ZNA vs FT0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/ACZNAFV0A", "ZNA vs FV0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNCFV0A", "ZNC vs FV0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/AZNAFT0C", "ZNA vs FT0C, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/CZNAFT0C", "ZNA vs FT0C, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/ACZNAFT0C", "ZNA vs FT0C, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/AZNCFT0A", "ZNC vs FT0A, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/CZNCFT0A", "ZNC vs FT0A, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/ACZNCFT0A", "ZNC vs FT0A, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 399.5}}}); + registry.add("ZDC/AZNCFT0C", "ZNC vs FT0C, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/CZNCFT0C", "ZNC vs FT0C, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/ACZNCFT0C", "ZNC vs FT0C, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 299.5}}}); + registry.add("ZDC/AZNAFDDA", "ZNA vs FDDA, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNAFDDA", "ZNA vs FDDA, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNAFDDA", "ZNA vs FDDA, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/AZNAFDDC", "ZNA vs FDDC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNAFDDC", "ZNA vs FDDC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNAFDDC", "ZNA vs FDDC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/AZNCFDDA", "ZNC vs FDDA, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNCFDDA", "ZNC vs FDDA, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNCFDDA", "ZNC vs FDDA, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/AZNCFDDC", "ZNC vs FDDC, A Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/CZNCFDDC", "ZNC vs FDDC, C Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("ZDC/ACZNCFDDC", "ZNC vs FDDC, AC Gap", {HistType::kTH2F, {{500, -5.5, 199.5}, {100, -.5, 199.5}}}); + registry.add("FIT/TAFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TAFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TTFT0", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {40, -5.5, 34.5}}}); + registry.add("FIT/TTFT01", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {40, -5.5, 34.5}}}); + registry.add("FIT/TTFT02", "Time FT0A vs Time FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {40, -5.5, 34.5}}}); + registry.add("FIT/ACFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5}, {100, -.5, 99.5}}}); + registry.add("FIT/AFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5}, {100, -.5, 99.5}}}); + registry.add("FIT/CFT0", "Amp FITA vs Amp FITC", {HistType::kTH2F, {{100, -.5, 99.5}, {100, -.5, 99.5}}}); + registry.add("FIT/AFT0AC", "Amp FT0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/CFT0AC", "Amp FT0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/ACFT0AC", "Amp FT0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/AFTV0A", "Amp FT0A vs Amp FV0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/CFTV0A", "Amp FT0A vs Amp FV0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/ACFTV0A", "Amp FT0A vs Amp FV0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/AFTV0C", "Amp FV0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/CFTV0C", "Amp FV0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/ACFTV0C", "Amp FV0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/AFDV0A", "Amp FDDA vs Amp FV0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/CFDV0A", "Amp FDDA vs Amp FV0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/ACFDV0A", "Amp FDDA vs Amp FV0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/AFDT0A", "Amp FDDA vs Amp FT0A, A Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/CFDT0A", "Amp FDDA vs Amp FT0A, C Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/ACFDT0A", "Amp FDDA vs Amp FT0A, AC Gap", {HistType::kTH2F, {{100, -.5, 499.5}, {100, -.5, 499.5}}}); + registry.add("FIT/TAFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TAFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TAFDDC", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TCFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TCFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TCFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TCFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/TCFDDC", "Time vs Amp FDDC", {HistType::kTH2F, {{40, -5.5, 34.5}, {50, -.5, 149.5}}}); + registry.add("FIT/MFV0A", "Track number vs Amp FV0A", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -.5, 499.5}}}); + registry.add("FIT/MFT0A", "Track number vs Amp FT0A", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -.5, 499.5}}}); + registry.add("FIT/MFT0C", "Track number vs Amp FT0C", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -.5, 399.5}}}); + registry.add("FIT/MFDDA", "Track number vs Amp FDDA", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -.5, 499.5}}}); + registry.add("FIT/MFDDC", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -.5, 499.5}}}); + registry.add("ZDC/MAZNA", "Track number vs Amp FV0A", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); + registry.add("ZDC/MAZNC", "Track number vs Amp FT0A", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); + registry.add("ZDC/MCZNC", "Track number vs Amp FT0C", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); + registry.add("ZDC/MCZNA", "Track number vs Amp FDDA", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); + registry.add("ZDC/MACZNA", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); + registry.add("ZDC/MACZNC", "Track number vs Amp FDDC", {HistType::kTH2F, {{50, -.5, 49.5}, {250, -5.5, 149.5}}}); } // define data types using UDCollisionsFull = soa::Join; // UDCollisions - //using UDCollisionsFull = soa::Join; // UDCollisions + // using UDCollisionsFull = soa::Join; // UDCollisions using UDCollisionFull = UDCollisionsFull::iterator; using UDTracksFull = soa::Join; -// using UDTracksFull = soa::Join; + // using UDTracksFull = soa::Join; void process(UDCollisionFull const& dgcand, UDTracksFull const& dgtracks) { @@ -267,9 +267,10 @@ struct SG_FIT_Analyzer { // UDTutorial01 // select PV contributors Partition PVContributors = aod::udtrack::isPVContributor == true; PVContributors.bindTable(dgtracks); - if(PVContributors.size() > 50)return; + if (PVContributors.size() > 50) + return; registry.get(HIST("collisions/multiplicityPVC"))->Fill(PVContributors.size(), 1.); - bool tof = false; + bool tof = false; // relative BC number auto bcnum = dgcand.globalBC() % o2::constants::lhc::LHCMaxBunches; registry.get(HIST("collisions/BC"))->Fill(bcnum, 1.); @@ -288,223 +289,232 @@ struct SG_FIT_Analyzer { // UDTutorial01 int z1pva = 0; int z1pvc = 0; int z1pvac = 0; - registry.get(HIST("FIT/TFT0A"))->Fill(dgcand.timeFT0A(), 1.); - registry.get(HIST("FIT/TFT0C"))->Fill(dgcand.timeFT0C(), 1.); - registry.get(HIST("FIT/TFV0A"))->Fill(dgcand.timeFV0A(), 1.); - registry.get(HIST("FIT/TFDDA"))->Fill(dgcand.timeFDDA(), 1.); - registry.get(HIST("FIT/TFDDC"))->Fill(dgcand.timeFDDC(), 1.); - registry.get(HIST("ZDC/tZNA"))->Fill(dgcand.timeZNA(), 1.); - registry.get(HIST("ZDC/tZNC"))->Fill(dgcand.timeZNC(), 1.); -// if (dgcand.gapSide() == 0) { - registry.get(HIST("FIT/TAFT0A"))->Fill(dgcand.timeFT0A(), dgcand.totalFT0AmplitudeA()); - if (dgcand.gapSide() == 0) registry.get(HIST("FIT/TTFT0"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); - if (dgcand.gapSide() == 1) registry.get(HIST("FIT/TTFT01"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); - if (dgcand.gapSide() == 2) registry.get(HIST("FIT/TTFT02"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); - if (dgcand.gapSide() == 0) - registry.get(HIST("FIT/AFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 1) - registry.get(HIST("FIT/CFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 2) - registry.get(HIST("FIT/ACFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 0) - registry.get(HIST("FIT/AFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 1) - registry.get(HIST("FIT/CFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 2) - registry.get(HIST("FIT/ACFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 0) - registry.get(HIST("FIT/AFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 1) - registry.get(HIST("FIT/CFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 2) - registry.get(HIST("FIT/ACFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); - if (dgcand.gapSide() == 0) - registry.get(HIST("FIT/AFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 1) - registry.get(HIST("FIT/CFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 2) - registry.get(HIST("FIT/ACFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); - if (dgcand.gapSide() == 0) - registry.get(HIST("FIT/AFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); - if (dgcand.gapSide() == 1) - registry.get(HIST("FIT/CFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); - if (dgcand.gapSide() == 2) - registry.get(HIST("FIT/ACFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); + registry.get(HIST("FIT/TFT0A"))->Fill(dgcand.timeFT0A(), 1.); + registry.get(HIST("FIT/TFT0C"))->Fill(dgcand.timeFT0C(), 1.); + registry.get(HIST("FIT/TFV0A"))->Fill(dgcand.timeFV0A(), 1.); + registry.get(HIST("FIT/TFDDA"))->Fill(dgcand.timeFDDA(), 1.); + registry.get(HIST("FIT/TFDDC"))->Fill(dgcand.timeFDDC(), 1.); + registry.get(HIST("ZDC/tZNA"))->Fill(dgcand.timeZNA(), 1.); + registry.get(HIST("ZDC/tZNC"))->Fill(dgcand.timeZNC(), 1.); + // if (dgcand.gapSide() == 0) { + registry.get(HIST("FIT/TAFT0A"))->Fill(dgcand.timeFT0A(), dgcand.totalFT0AmplitudeA()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/TTFT0"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/TTFT01"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/TTFT02"))->Fill(dgcand.timeFT0A(), dgcand.timeFT0C()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/AFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/CFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/ACFT0AC"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/AFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/CFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/ACFTV0A"))->Fill(dgcand.totalFT0AmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/AFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/CFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/ACFTV0C"))->Fill(dgcand.totalFV0AmplitudeA(), dgcand.totalFT0AmplitudeC()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/AFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/CFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/ACFDV0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFV0AmplitudeA()); + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/AFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/CFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); + if (dgcand.gapSide() == 2) + registry.get(HIST("FIT/ACFDT0A"))->Fill(dgcand.totalFDDAmplitudeA(), dgcand.totalFT0AmplitudeA()); float totalA = dgcand.totalFDDAmplitudeA() + dgcand.totalFT0AmplitudeA() + dgcand.totalFV0AmplitudeA(); float totalC = dgcand.totalFDDAmplitudeC() + dgcand.totalFT0AmplitudeC(); - float zna = -1; - float znc = -1; - int an = 0; - int cn = 0; - if( dgcand.energyCommonZNC() > 0) znc = dgcand.energyCommonZNC() ; - if( dgcand.energyCommonZNA() > 0) zna = dgcand.energyCommonZNA() ; - if (zna > 0 && zna < 4) an = 1; - else if (zna > 4 && zna < 10) an = 2; - else if (zna > 10 && zna < 30) an = 3; - else if (zna > 30) an = 4; - if (znc > 0 && znc < 4) cn = 1; - else if (znc > 4 && znc < 10) cn = 2; - else if (znc > 10 && znc < 30) cn = 3; - else if (znc > 30) cn = 4; -//if (zna >0 && znc> 0) LOGF(info, "ZNA %f, ZNC %f", zna, znc); - registry.get(HIST("FIT/TAFT0C"))->Fill(dgcand.timeFT0C(), dgcand.totalFT0AmplitudeC()); - registry.get(HIST("FIT/TAFV0A"))->Fill(dgcand.timeFV0A(), dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/TAFDDA"))->Fill(dgcand.timeFDDA(), dgcand.totalFDDAmplitudeA()); - registry.get(HIST("FIT/TAFDDC"))->Fill(dgcand.timeFDDC(), dgcand.totalFDDAmplitudeC()); -// } + float zna = -1; + float znc = -1; + int an = 0; + int cn = 0; + if (dgcand.energyCommonZNC() > 0) + znc = dgcand.energyCommonZNC(); + if (dgcand.energyCommonZNA() > 0) + zna = dgcand.energyCommonZNA(); + if (zna > 0 && zna < 4) + an = 1; + else if (zna > 4 && zna < 10) + an = 2; + else if (zna > 10 && zna < 30) + an = 3; + else if (zna > 30) + an = 4; + if (znc > 0 && znc < 4) + cn = 1; + else if (znc > 4 && znc < 10) + cn = 2; + else if (znc > 10 && znc < 30) + cn = 3; + else if (znc > 30) + cn = 4; + // if (zna >0 && znc> 0) LOGF(info, "ZNA %f, ZNC %f", zna, znc); + registry.get(HIST("FIT/TAFT0C"))->Fill(dgcand.timeFT0C(), dgcand.totalFT0AmplitudeC()); + registry.get(HIST("FIT/TAFV0A"))->Fill(dgcand.timeFV0A(), dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/TAFDDA"))->Fill(dgcand.timeFDDA(), dgcand.totalFDDAmplitudeA()); + registry.get(HIST("FIT/TAFDDC"))->Fill(dgcand.timeFDDC(), dgcand.totalFDDAmplitudeC()); + // } if (dgcand.gapSide() == 1) { - registry.get(HIST("FIT/TCFT0A"))->Fill(dgcand.timeFT0A(), dgcand.totalFT0AmplitudeA()); - registry.get(HIST("FIT/TCFT0C"))->Fill(dgcand.timeFT0C(), dgcand.totalFT0AmplitudeC()); - registry.get(HIST("FIT/TCFV0A"))->Fill(dgcand.timeFV0A(), dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/TCFDDA"))->Fill(dgcand.timeFDDA(), dgcand.totalFDDAmplitudeA()); - registry.get(HIST("FIT/TCFDDC"))->Fill(dgcand.timeFDDC(), dgcand.totalFDDAmplitudeC()); + registry.get(HIST("FIT/TCFT0A"))->Fill(dgcand.timeFT0A(), dgcand.totalFT0AmplitudeA()); + registry.get(HIST("FIT/TCFT0C"))->Fill(dgcand.timeFT0C(), dgcand.totalFT0AmplitudeC()); + registry.get(HIST("FIT/TCFV0A"))->Fill(dgcand.timeFV0A(), dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/TCFDDA"))->Fill(dgcand.timeFDDA(), dgcand.totalFDDAmplitudeA()); + registry.get(HIST("FIT/TCFDDC"))->Fill(dgcand.timeFDDC(), dgcand.totalFDDAmplitudeC()); } if (dgcand.gapSide() == 0) { - registry.get(HIST("ZDC/AZNA"))->Fill(zna, 1.); - registry.get(HIST("ZDC/AZNC"))->Fill(znc, 1.); - registry.get(HIST("ZDC/AZNAC"))->Fill(zna, znc); - registry.get(HIST("ZDC/AZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/AZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/AZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/AZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/AZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/AZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/AZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/AZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/AZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("ZDC/AZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/AFT0"))->Fill(totalA, totalC); - registry.get(HIST("FIT/AFITA"))->Fill(totalA, 1.); - registry.get(HIST("FIT/AFITC"))->Fill(totalC, 1.); - registry.get(HIST("FIT/AFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - if (an&&cn)registry.get(HIST("FIT/zAFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - if (cn==0){ - registry.get(HIST("FIT/AFT0A0"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFV0A0"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C0"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - if (an==0)registry.get(HIST("FIT/AFT0C00"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/AFDDC0"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - } - else if (cn==1){ - registry.get(HIST("FIT/AFT0A1"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFV0A1"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C1"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/AFDDC1"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - } - else if (cn==2){ - registry.get(HIST("FIT/AFT0A2"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFV0A2"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C2"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/AFDDC2"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - } - else if (cn==3){ - registry.get(HIST("FIT/AFT0A3"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFV0A3"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C3"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/AFDDC3"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - } - else if (cn==4){ - registry.get(HIST("FIT/AFT0A4"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFV0A4"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFT0C4"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/AFDDC4"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - } - if (an&&cn)registry.get(HIST("FIT/zAFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + registry.get(HIST("ZDC/AZNA"))->Fill(zna, 1.); + registry.get(HIST("ZDC/AZNC"))->Fill(znc, 1.); + registry.get(HIST("ZDC/AZNAC"))->Fill(zna, znc); + registry.get(HIST("ZDC/AZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/AZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/AZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/AZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/AZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/AZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/AZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/AZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/AZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("ZDC/AZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/AFT0"))->Fill(totalA, totalC); + registry.get(HIST("FIT/AFITA"))->Fill(totalA, 1.); + registry.get(HIST("FIT/AFITC"))->Fill(totalC, 1.); + registry.get(HIST("FIT/AFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + if (an && cn) + registry.get(HIST("FIT/zAFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + if (cn == 0) { + registry.get(HIST("FIT/AFT0A0"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFV0A0"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C0"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + if (an == 0) + registry.get(HIST("FIT/AFT0C00"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/AFDDC0"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + } else if (cn == 1) { + registry.get(HIST("FIT/AFT0A1"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFV0A1"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C1"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/AFDDC1"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + } else if (cn == 2) { + registry.get(HIST("FIT/AFT0A2"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFV0A2"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C2"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/AFDDC2"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + } else if (cn == 3) { + registry.get(HIST("FIT/AFT0A3"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFV0A3"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C3"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/AFDDC3"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + } else if (cn == 4) { + registry.get(HIST("FIT/AFT0A4"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFV0A4"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFT0C4"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/AFDDC4"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + } + if (an && cn) + registry.get(HIST("FIT/zAFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - registry.get(HIST("FIT/AFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/AFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - registry.get(HIST("FIT/AFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - registry.get(HIST("FIT/MFT0C"))->Fill(PVContributors.size(), dgcand.totalFT0AmplitudeC()); - registry.get(HIST("FIT/MFDDC"))->Fill(PVContributors.size(), dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/MAZNA"))->Fill(PVContributors.size(), zna); - registry.get(HIST("ZDC/MAZNC"))->Fill(PVContributors.size(), znc); + registry.get(HIST("FIT/AFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/AFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + registry.get(HIST("FIT/AFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + registry.get(HIST("FIT/MFT0C"))->Fill(PVContributors.size(), dgcand.totalFT0AmplitudeC()); + registry.get(HIST("FIT/MFDDC"))->Fill(PVContributors.size(), dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/MAZNA"))->Fill(PVContributors.size(), zna); + registry.get(HIST("ZDC/MAZNC"))->Fill(PVContributors.size(), znc); } if (dgcand.gapSide() == 1) { - registry.get(HIST("ZDC/CZNA"))->Fill(zna, 1.); - registry.get(HIST("ZDC/CZNC"))->Fill(znc, 1.); - registry.get(HIST("ZDC/CZNAC"))->Fill(zna, znc); - registry.get(HIST("ZDC/CZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/CZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/CZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/CZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/CZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/CZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/CZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/CZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/CZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("ZDC/CZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/CFT0"))->Fill(totalA, totalC); - registry.get(HIST("FIT/CFITA"))->Fill(totalA, 1.); - registry.get(HIST("FIT/CFITC"))->Fill(totalC, 1.); - registry.get(HIST("FIT/MFT0A"))->Fill(PVContributors.size(), dgcand.totalFT0AmplitudeA()); - registry.get(HIST("FIT/MFV0A"))->Fill(PVContributors.size(), dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/MFDDA"))->Fill(PVContributors.size(), dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/MCZNA"))->Fill(PVContributors.size(), zna); - registry.get(HIST("ZDC/MCZNC"))->Fill(PVContributors.size(), znc); - registry.get(HIST("FIT/CFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/CFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); -if(an&&cn){ - registry.get(HIST("FIT/zCFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/zCFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/zCFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); -} - if (an==0){ - registry.get(HIST("FIT/CFT0A0"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - if (cn==0)registry.get(HIST("FIT/CFT0A00"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFV0A0"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA0"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - } - else if (an==1){ - registry.get(HIST("FIT/CFT0A1"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFV0A1"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA1"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - } - else if (an==2){ - registry.get(HIST("FIT/CFT0A2"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFV0A2"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA2"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - } - else if (an==3){ - registry.get(HIST("FIT/CFT0A3"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFV0A3"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA3"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - } - else if (an==4){ - registry.get(HIST("FIT/CFT0A4"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFV0A4"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDA4"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - } + registry.get(HIST("ZDC/CZNA"))->Fill(zna, 1.); + registry.get(HIST("ZDC/CZNC"))->Fill(znc, 1.); + registry.get(HIST("ZDC/CZNAC"))->Fill(zna, znc); + registry.get(HIST("ZDC/CZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/CZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/CZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/CZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/CZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/CZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/CZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/CZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/CZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("ZDC/CZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/CFT0"))->Fill(totalA, totalC); + registry.get(HIST("FIT/CFITA"))->Fill(totalA, 1.); + registry.get(HIST("FIT/CFITC"))->Fill(totalC, 1.); + registry.get(HIST("FIT/MFT0A"))->Fill(PVContributors.size(), dgcand.totalFT0AmplitudeA()); + registry.get(HIST("FIT/MFV0A"))->Fill(PVContributors.size(), dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/MFDDA"))->Fill(PVContributors.size(), dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/MCZNA"))->Fill(PVContributors.size(), zna); + registry.get(HIST("ZDC/MCZNC"))->Fill(PVContributors.size(), znc); + registry.get(HIST("FIT/CFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/CFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + if (an && cn) { + registry.get(HIST("FIT/zCFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/zCFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/zCFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } + if (an == 0) { + registry.get(HIST("FIT/CFT0A0"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + if (cn == 0) + registry.get(HIST("FIT/CFT0A00"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFV0A0"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA0"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } else if (an == 1) { + registry.get(HIST("FIT/CFT0A1"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFV0A1"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA1"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } else if (an == 2) { + registry.get(HIST("FIT/CFT0A2"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFV0A2"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA2"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } else if (an == 3) { + registry.get(HIST("FIT/CFT0A3"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFV0A3"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA3"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } else if (an == 4) { + registry.get(HIST("FIT/CFT0A4"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFV0A4"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDA4"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + } - registry.get(HIST("FIT/CFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - registry.get(HIST("FIT/CFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + registry.get(HIST("FIT/CFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + registry.get(HIST("FIT/CFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); } if (dgcand.gapSide() == 2) { - registry.get(HIST("ZDC/ACZNA"))->Fill(zna, 1.); - registry.get(HIST("ZDC/ACZNC"))->Fill(znc, 1.); - registry.get(HIST("ZDC/ACZNAC"))->Fill(zna, znc); - registry.get(HIST("ZDC/ACZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/ACZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/ACZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); - registry.get(HIST("ZDC/ACZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); - registry.get(HIST("ZDC/ACZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/ACZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/ACZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); - registry.get(HIST("ZDC/ACZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); - registry.get(HIST("ZDC/ACZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("ZDC/ACZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); - registry.get(HIST("FIT/ACFT0"))->Fill(totalA, totalC); - registry.get(HIST("FIT/ACFITA"))->Fill(totalA, 1.); - registry.get(HIST("FIT/ACFITC"))->Fill(totalC, 1.); - registry.get(HIST("FIT/ACFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("FIT/ACFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - registry.get(HIST("FIT/ACFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); - registry.get(HIST("FIT/ACFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); - registry.get(HIST("FIT/ACFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); - registry.get(HIST("ZDC/MACZNA"))->Fill(PVContributors.size(), zna); - registry.get(HIST("ZDC/MACZNC"))->Fill(PVContributors.size(), znc); + registry.get(HIST("ZDC/ACZNA"))->Fill(zna, 1.); + registry.get(HIST("ZDC/ACZNC"))->Fill(znc, 1.); + registry.get(HIST("ZDC/ACZNAC"))->Fill(zna, znc); + registry.get(HIST("ZDC/ACZNAFT0A"))->Fill(zna, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/ACZNAFT0C"))->Fill(zna, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/ACZNCFT0A"))->Fill(znc, dgcand.totalFT0AmplitudeA()); + registry.get(HIST("ZDC/ACZNCFT0C"))->Fill(znc, dgcand.totalFT0AmplitudeC()); + registry.get(HIST("ZDC/ACZNAFDDA"))->Fill(zna, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/ACZNAFDDC"))->Fill(zna, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/ACZNCFDDA"))->Fill(znc, dgcand.totalFDDAmplitudeA()); + registry.get(HIST("ZDC/ACZNCFDDC"))->Fill(znc, dgcand.totalFDDAmplitudeC()); + registry.get(HIST("ZDC/ACZNAFV0A"))->Fill(zna, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("ZDC/ACZNCFV0A"))->Fill(znc, dgcand.totalFV0AmplitudeA()); + registry.get(HIST("FIT/ACFT0"))->Fill(totalA, totalC); + registry.get(HIST("FIT/ACFITA"))->Fill(totalA, 1.); + registry.get(HIST("FIT/ACFITC"))->Fill(totalC, 1.); + registry.get(HIST("FIT/ACFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("FIT/ACFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + registry.get(HIST("FIT/ACFV0A"))->Fill(dgcand.totalFV0AmplitudeA(), 1.); + registry.get(HIST("FIT/ACFDDA"))->Fill(dgcand.totalFDDAmplitudeA(), 1.); + registry.get(HIST("FIT/ACFDDC"))->Fill(dgcand.totalFDDAmplitudeC(), 1.); + registry.get(HIST("ZDC/MACZNA"))->Fill(PVContributors.size(), zna); + registry.get(HIST("ZDC/MACZNC"))->Fill(PVContributors.size(), znc); } for (auto track : dgtracks) { @@ -513,15 +523,22 @@ if(an&&cn){ registry.get(HIST("tracks/QCAll"))->Fill(2., track.hasTPC() * 1.); registry.get(HIST("tracks/QCAll"))->Fill(3., track.hasTRD() * 1.); registry.get(HIST("tracks/QCAll"))->Fill(4., track.hasTOF() * 1.); - if (track.tofChi2() > -10) tof = true; + if (track.tofChi2() > -10) + tof = true; auto vtrk = TVector3(track.px(), track.py(), track.pz()); registry.get(HIST("tracks/ptAll"))->Fill(track.pt(), 1.); - if (dgcand.gapSide() == 0) registry.get(HIST("tracks/etavsptA"))->Fill(vtrk.Eta(), track.pt(), 1.); - if (dgcand.gapSide() == 1) registry.get(HIST("tracks/etavsptC"))->Fill(vtrk.Eta(), track.pt(), 1.); - if (dgcand.gapSide() == 2) registry.get(HIST("tracks/etavsptAC"))->Fill(vtrk.Eta(), track.pt(), 1.); - if (dgcand.gapSide() == 0 && cn ==4) registry.get(HIST("tracks/etaA"))->Fill(vtrk.Eta(), 1.); - if (dgcand.gapSide() == 1 && an ==4) registry.get(HIST("tracks/etaC"))->Fill(vtrk.Eta(), 1.); - if (dgcand.gapSide() == 2 && an > 2 && cn > 2) registry.get(HIST("tracks/etaAC"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 0) + registry.get(HIST("tracks/etavsptA"))->Fill(vtrk.Eta(), track.pt(), 1.); + if (dgcand.gapSide() == 1) + registry.get(HIST("tracks/etavsptC"))->Fill(vtrk.Eta(), track.pt(), 1.); + if (dgcand.gapSide() == 2) + registry.get(HIST("tracks/etavsptAC"))->Fill(vtrk.Eta(), track.pt(), 1.); + if (dgcand.gapSide() == 0 && cn == 4) + registry.get(HIST("tracks/etaA"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 1 && an == 4) + registry.get(HIST("tracks/etaC"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 2 && an > 2 && cn > 2) + registry.get(HIST("tracks/etaAC"))->Fill(vtrk.Eta(), 1.); auto signalTPC = track.tpcSignal() * track.sign(); registry.get(HIST("tracks/TPCSignalvspAll"))->Fill(vtrk.Mag(), signalTPC, 1.); @@ -534,49 +551,72 @@ if(an&&cn){ registry.get(HIST("tracks/QCPVC"))->Fill(3., track.hasTRD() * 1.); registry.get(HIST("tracks/QCPVC"))->Fill(4., track.hasTOF() * 1.); registry.get(HIST("tracks/ptPVC"))->Fill(track.pt(), 1.); -// registry.get(HIST("tracks/etavsptPVC"))->Fill(vtrk.Eta(), track.pt(), 1.); + // registry.get(HIST("tracks/etavsptPVC"))->Fill(vtrk.Eta(), track.pt(), 1.); if (dgcand.gapSide() == 0) { - registry.get(HIST("tracks/etaApv"))->Fill(vtrk.Eta(), 1.); - if (!an) registry.get(HIST("tracks/eta2Apv"))->Fill(vtrk.Eta(), 1.); - if (!an) z0pva++; - else if (an==4) z1pva++; - pva++; - } + registry.get(HIST("tracks/etaApv"))->Fill(vtrk.Eta(), 1.); + if (!an) + registry.get(HIST("tracks/eta2Apv"))->Fill(vtrk.Eta(), 1.); + if (!an) + z0pva++; + else if (an == 4) + z1pva++; + pva++; + } if (dgcand.gapSide() == 1) { - registry.get(HIST("tracks/etaCpv"))->Fill(vtrk.Eta(), 1.); - if (!cn) registry.get(HIST("tracks/eta2Cpv"))->Fill(vtrk.Eta(), 1.); - if (!cn) z0pvc++; - else if (cn==4) z1pvc++; - pvc++; - } + registry.get(HIST("tracks/etaCpv"))->Fill(vtrk.Eta(), 1.); + if (!cn) + registry.get(HIST("tracks/eta2Cpv"))->Fill(vtrk.Eta(), 1.); + if (!cn) + z0pvc++; + else if (cn == 4) + z1pvc++; + pvc++; + } if (dgcand.gapSide() == 2) { - registry.get(HIST("tracks/etaACpv"))->Fill(vtrk.Eta(), 1.); - if (!an&&!cn) registry.get(HIST("tracks/eta2ACpv"))->Fill(vtrk.Eta(), 1.); - if (!an&&!cn) z0pvac++; - else if (an>=3 && cn>=3) z1pvac++; - pvac++; - } + registry.get(HIST("tracks/etaACpv"))->Fill(vtrk.Eta(), 1.); + if (!an && !cn) + registry.get(HIST("tracks/eta2ACpv"))->Fill(vtrk.Eta(), 1.); + if (!an && !cn) + z0pvac++; + else if (an >= 3 && cn >= 3) + z1pvac++; + pvac++; + } registry.get(HIST("tracks/TPCSignalvspPVC"))->Fill(vtrk.Mag(), signalTPC, 1.); registry.get(HIST("tracks/TOFSignalvspPVC"))->Fill(vtrk.Mag(), signalTOF, 1.); } } - if (pva) registry.get(HIST("collisions/multiplicityPVCA"))->Fill(pva, 1.); - if (pvc) registry.get(HIST("collisions/multiplicityPVCC"))->Fill(pvc, 1.); - if (pvac) registry.get(HIST("collisions/multiplicityPVCAC"))->Fill(pvac, 1.); - if (pva) registry.get(HIST("collisions/multiplicityZ0PVCA"))->Fill(z0pva, 1.); - if (pvc) registry.get(HIST("collisions/multiplicityZ0PVCC"))->Fill(z0pvc, 1.); - if (pvac) registry.get(HIST("collisions/multiplicityZ0PVCAC"))->Fill(z0pvac, 1.); - if (pva) registry.get(HIST("collisions/multiplicityZ1PVCA"))->Fill(z1pva, 1.); - if (pvc) registry.get(HIST("collisions/multiplicityZ1PVCC"))->Fill(z1pvc, 1.); - if (pvac) registry.get(HIST("collisions/multiplicityZ1PVCAC"))->Fill(z1pvac, 1.); -if (tof){ - if (dgcand.gapSide() == 0) registry.get(HIST("FIT/tAFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); - if (dgcand.gapSide() == 1) registry.get(HIST("FIT/tCFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); - registry.get(HIST("collisions/multiplicityTPVC"))->Fill(PVContributors.size(), 1.); - if (pva) registry.get(HIST("collisions/multiplicityTPVCA"))->Fill(pva, 1.); - if (pvc) registry.get(HIST("collisions/multiplicityTPVCC"))->Fill(pvc, 1.); - if (pvac) registry.get(HIST("collisions/multiplicityTPVCAC"))->Fill(pvac, 1.); -} + if (pva) + registry.get(HIST("collisions/multiplicityPVCA"))->Fill(pva, 1.); + if (pvc) + registry.get(HIST("collisions/multiplicityPVCC"))->Fill(pvc, 1.); + if (pvac) + registry.get(HIST("collisions/multiplicityPVCAC"))->Fill(pvac, 1.); + if (pva) + registry.get(HIST("collisions/multiplicityZ0PVCA"))->Fill(z0pva, 1.); + if (pvc) + registry.get(HIST("collisions/multiplicityZ0PVCC"))->Fill(z0pvc, 1.); + if (pvac) + registry.get(HIST("collisions/multiplicityZ0PVCAC"))->Fill(z0pvac, 1.); + if (pva) + registry.get(HIST("collisions/multiplicityZ1PVCA"))->Fill(z1pva, 1.); + if (pvc) + registry.get(HIST("collisions/multiplicityZ1PVCC"))->Fill(z1pvc, 1.); + if (pvac) + registry.get(HIST("collisions/multiplicityZ1PVCAC"))->Fill(z1pvac, 1.); + if (tof) { + if (dgcand.gapSide() == 0) + registry.get(HIST("FIT/tAFT0C"))->Fill(dgcand.totalFT0AmplitudeC(), 1.); + if (dgcand.gapSide() == 1) + registry.get(HIST("FIT/tCFT0A"))->Fill(dgcand.totalFT0AmplitudeA(), 1.); + registry.get(HIST("collisions/multiplicityTPVC"))->Fill(PVContributors.size(), 1.); + if (pva) + registry.get(HIST("collisions/multiplicityTPVCA"))->Fill(pva, 1.); + if (pvc) + registry.get(HIST("collisions/multiplicityTPVCC"))->Fill(pvc, 1.); + if (pvac) + registry.get(HIST("collisions/multiplicityTPVCAC"))->Fill(pvac, 1.); + } // fill FIT histograms for (auto bit = 0; bit < 33; bit++) { registry.get(HIST("FIT/BBFV0A"))->Fill(bit - 16, TESTBIT(dgcand.bbFV0Apf(), bit));