diff --git a/PWGUD/Tasks/SG_FIT_Analyzer.cxx b/PWGUD/Tasks/SG_FIT_Analyzer.cxx index 9b5830865df..e83e61d8fb0 100644 --- a/PWGUD/Tasks/SG_FIT_Analyzer.cxx +++ b/PWGUD/Tasks/SG_FIT_Analyzer.cxx @@ -33,7 +33,7 @@ struct SG_FIT_Analyzer { // UDTutorial01 ConfigurableAxis etaAxis{"etaAxis", {300, -1.5, 1.5}, ""}; ConfigurableAxis sigTPCAxis{"sigTPCAxis", {100, -100.0, 100.0}, ""}; ConfigurableAxis sigTOFAxis{"sigTOFAxis", {100, -100.0, 100.0}, ""}; - ConfigurableAxis multAxis{"multAxis", {51,-.5,50.5}, ""}; + ConfigurableAxis multAxis{"multAxis", {51, -.5, 50.5}, ""}; ConfigurableAxis FitAxis{"FitAxis", {5000, -0.5, 9999.5}, ""}; ConfigurableAxis ZDCAxis{"FITAxis", {1000, -2.5, 199.5}, ""}; // initialize histogram registry @@ -181,89 +181,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, {{axiszdc},{axiszdc}}}); - registry.add("ZDC/CZNAC", "ZNA vs ZNC, C Gap", {HistType::kTH2F, {{axiszdc},{axiszdc}}}); - registry.add("ZDC/ACZNAC", "ZNA vs ZNC, AC Gap", {HistType::kTH2F, {{axiszdc},{axiszdc}}}); - registry.add("ZDC/AZNAFT0A", "ZNA vs FT0A, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNAFV0A", "ZNA vs FV0A, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNCFV0A", "ZNC vs FV0A, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNAFT0A", "ZNA vs FT0A, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNAFV0A", "ZNA vs FV0A, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNCFV0A", "ZNC vs FV0A, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNAFT0A", "ZNA vs FT0A, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNAFV0A", "ZNA vs FV0A, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNCFV0A", "ZNC vs FV0A, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNAFT0C", "ZNA vs FT0C, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNAFT0C", "ZNA vs FT0C, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNAFT0C", "ZNA vs FT0C, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNCFT0A", "ZNC vs FT0A, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNCFT0A", "ZNC vs FT0A, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNCFT0A", "ZNC vs FT0A, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNCFT0C", "ZNC vs FT0C, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNCFT0C", "ZNC vs FT0C, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNCFT0C", "ZNC vs FT0C, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNAFDDA", "ZNA vs FDDA, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNAFDDA", "ZNA vs FDDA, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNAFDDA", "ZNA vs FDDA, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNAFDDC", "ZNA vs FDDC, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNAFDDC", "ZNA vs FDDC, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNAFDDC", "ZNA vs FDDC, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNCFDDA", "ZNC vs FDDA, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNCFDDA", "ZNC vs FDDA, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNCFDDA", "ZNC vs FDDA, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/AZNCFDDC", "ZNC vs FDDC, A Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/CZNCFDDC", "ZNC vs FDDC, C Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("ZDC/ACZNCFDDC", "ZNC vs FDDC, AC Gap", {HistType::kTH2F, {{axiszdc},{axisfit}}}); - registry.add("FIT/TAFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TAFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - 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, {{axisfit},{axisfit}}}); - registry.add("FIT/CFT0AC", "Amp FT0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/ACFT0AC", "Amp FT0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/AFTV0A", "Amp FT0A vs Amp FV0A, A Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/CFTV0A", "Amp FT0A vs Amp FV0A, C Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/ACFTV0A", "Amp FT0A vs Amp FV0A, AC Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/AFTV0C", "Amp FV0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/CFTV0C", "Amp FV0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/ACFTV0C", "Amp FV0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/AFDV0A", "Amp FDDA vs Amp FV0A, A Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/CFDV0A", "Amp FDDA vs Amp FV0A, C Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/ACFDV0A", "Amp FDDA vs Amp FV0A, AC Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/AFDT0A", "Amp FDDA vs Amp FT0A, A Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/CFDT0A", "Amp FDDA vs Amp FT0A, C Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/ACFDT0A", "Amp FDDA vs Amp FT0A, AC Gap", {HistType::kTH2F, {{axisfit},{axisfit}}}); - registry.add("FIT/TAFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TAFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TAFDDC", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TCFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TCFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TCFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TCFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/TCFDDC", "Time vs Amp FDDC", {HistType::kTH2F, {{40, -5.5, 34.5},{axisfit}}}); - registry.add("FIT/MFV0A", "Track number vs Amp FV0A", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("FIT/MFT0A", "Track number vs Amp FT0A", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("FIT/MFT0C", "Track number vs Amp FT0C", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("FIT/MFDDA", "Track number vs Amp FDDA", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("FIT/MFDDC", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MAZNA", "Track number vs Amp FV0A", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MAZNC", "Track number vs Amp FT0A", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MCZNC", "Track number vs Amp FT0C", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MCZNA", "Track number vs Amp FDDA", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MACZNA", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult},{axisfit}}}); - registry.add("ZDC/MACZNC", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult},{axisfit}}}); + registry.add("ZDC/AZNAC", "ZNA vs ZNC, A Gap", {HistType::kTH2F, {{axiszdc}, {axiszdc}}}); + registry.add("ZDC/CZNAC", "ZNA vs ZNC, C Gap", {HistType::kTH2F, {{axiszdc}, {axiszdc}}}); + registry.add("ZDC/ACZNAC", "ZNA vs ZNC, AC Gap", {HistType::kTH2F, {{axiszdc}, {axiszdc}}}); + registry.add("ZDC/AZNAFT0A", "ZNA vs FT0A, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNAFV0A", "ZNA vs FV0A, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNCFV0A", "ZNC vs FV0A, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNAFT0A", "ZNA vs FT0A, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNAFV0A", "ZNA vs FV0A, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNCFV0A", "ZNC vs FV0A, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNAFT0A", "ZNA vs FT0A, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNAFV0A", "ZNA vs FV0A, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNCFV0A", "ZNC vs FV0A, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNAFT0C", "ZNA vs FT0C, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNAFT0C", "ZNA vs FT0C, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNAFT0C", "ZNA vs FT0C, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNCFT0A", "ZNC vs FT0A, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNCFT0A", "ZNC vs FT0A, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNCFT0A", "ZNC vs FT0A, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNCFT0C", "ZNC vs FT0C, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNCFT0C", "ZNC vs FT0C, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNCFT0C", "ZNC vs FT0C, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNAFDDA", "ZNA vs FDDA, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNAFDDA", "ZNA vs FDDA, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNAFDDA", "ZNA vs FDDA, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNAFDDC", "ZNA vs FDDC, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNAFDDC", "ZNA vs FDDC, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNAFDDC", "ZNA vs FDDC, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNCFDDA", "ZNC vs FDDA, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNCFDDA", "ZNC vs FDDA, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNCFDDA", "ZNC vs FDDA, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/AZNCFDDC", "ZNC vs FDDC, A Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/CZNCFDDC", "ZNC vs FDDC, C Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("ZDC/ACZNCFDDC", "ZNC vs FDDC, AC Gap", {HistType::kTH2F, {{axiszdc}, {axisfit}}}); + registry.add("FIT/TAFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TAFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + 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, {{axisfit}, {axisfit}}}); + registry.add("FIT/CFT0AC", "Amp FT0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/ACFT0AC", "Amp FT0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/AFTV0A", "Amp FT0A vs Amp FV0A, A Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/CFTV0A", "Amp FT0A vs Amp FV0A, C Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/ACFTV0A", "Amp FT0A vs Amp FV0A, AC Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/AFTV0C", "Amp FV0A vs Amp FT0C, A Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/CFTV0C", "Amp FV0A vs Amp FT0C, C Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/ACFTV0C", "Amp FV0A vs Amp FT0C, AC Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/AFDV0A", "Amp FDDA vs Amp FV0A, A Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/CFDV0A", "Amp FDDA vs Amp FV0A, C Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/ACFDV0A", "Amp FDDA vs Amp FV0A, AC Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/AFDT0A", "Amp FDDA vs Amp FT0A, A Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/CFDT0A", "Amp FDDA vs Amp FT0A, C Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/ACFDT0A", "Amp FDDA vs Amp FT0A, AC Gap", {HistType::kTH2F, {{axisfit}, {axisfit}}}); + registry.add("FIT/TAFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TAFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TAFDDC", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TCFV0A", "Time vs Amp FV0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TCFT0A", "Time vs Amp FT0A", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TCFT0C", "Time vs Amp FT0C", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TCFDDA", "Time vs Amp FDDA", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/TCFDDC", "Time vs Amp FDDC", {HistType::kTH2F, {{40, -5.5, 34.5}, {axisfit}}}); + registry.add("FIT/MFV0A", "Track number vs Amp FV0A", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("FIT/MFT0A", "Track number vs Amp FT0A", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("FIT/MFT0C", "Track number vs Amp FT0C", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("FIT/MFDDA", "Track number vs Amp FDDA", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("FIT/MFDDC", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MAZNA", "Track number vs Amp FV0A", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MAZNC", "Track number vs Amp FT0A", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MCZNC", "Track number vs Amp FT0C", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MCZNA", "Track number vs Amp FDDA", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MACZNA", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult}, {axisfit}}}); + registry.add("ZDC/MACZNC", "Track number vs Amp FDDC", {HistType::kTH2F, {{axismult}, {axisfit}}}); } // 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) { @@ -276,9 +276,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.); @@ -297,227 +298,236 @@ 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.); -// 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.); + // 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("ZDC/tCZNA"))->Fill(dgcand.timeZNA(), 1.); - registry.get(HIST("ZDC/tCZNC"))->Fill(dgcand.timeZNC(), 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("ZDC/tCZNA"))->Fill(dgcand.timeZNA(), 1.); + registry.get(HIST("ZDC/tCZNC"))->Fill(dgcand.timeZNC(), 1.); } if (dgcand.gapSide() == 0) { - registry.get(HIST("ZDC/tAZNA"))->Fill(dgcand.timeZNA(), 1.); - registry.get(HIST("ZDC/tAZNC"))->Fill(dgcand.timeZNC(), 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("ZDC/tAZNA"))->Fill(dgcand.timeZNA(), 1.); + registry.get(HIST("ZDC/tAZNC"))->Fill(dgcand.timeZNC(), 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/tACZNA"))->Fill(dgcand.timeZNA(), 1.); - registry.get(HIST("ZDC/tACZNC"))->Fill(dgcand.timeZNC(), 1.); - 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/tACZNA"))->Fill(dgcand.timeZNA(), 1.); + registry.get(HIST("ZDC/tACZNC"))->Fill(dgcand.timeZNC(), 1.); + 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) { @@ -526,70 +536,100 @@ 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/etaA"))->Fill(vtrk.Eta(), 1.); - if (dgcand.gapSide() == 1 ) registry.get(HIST("tracks/etaC"))->Fill(vtrk.Eta(), 1.); - if (dgcand.gapSide() == 2 ) registry.get(HIST("tracks/etaAC"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 0) + registry.get(HIST("tracks/etaA"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 1) + registry.get(HIST("tracks/etaC"))->Fill(vtrk.Eta(), 1.); + if (dgcand.gapSide() == 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.); auto signalTOF = track.tofSignal() * track.sign() / 1.E3; registry.get(HIST("tracks/TOFSignalvspAll"))->Fill(vtrk.Mag(), signalTOF, 1.); if (track.isPVContributor()) { - 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) + 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.); registry.get(HIST("tracks/QCPVC"))->Fill(0., 1.); registry.get(HIST("tracks/QCPVC"))->Fill(1., track.hasITS() * 1.); registry.get(HIST("tracks/QCPVC"))->Fill(2., track.hasTPC() * 1.); 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));