[ALICE3,PWGDQ] Implementation of the ALICE3 DQ data model and analysis tasks

ALICE3/TableProducer/CMakeLists.txt

@@ -1,4 +1,4 @@
 # Copyright 2019-2020 CERN and copyright holders of ALICE O2.
 # See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
 # All rights not expressly granted are reserved.
 #
@@ -11,37 +11,37 @@
 
 add_subdirectory(OTF)
 
 o2physics_add_dpl_workflow(alice3-trackselection
                     SOURCES alice3-trackselection.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-trackextension
                     SOURCES alice3-trackextension.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2::ReconstructionDataFormats
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-centrality
                     SOURCES alice3-centrality.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-decaypreselector
                     SOURCES alice3-decaypreselector.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DCAFitter
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-decayfinder
                     SOURCES alice3-decayfinder.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DCAFitter
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-multicharm-table
                     SOURCES alice3-multicharmTable.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DCAFitter
                     COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(alice3-correlatorddbar
                     SOURCES alice3-correlatorDDbar.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DCAFitter
                     COMPONENT_NAME Analysis)
@@ -60,3 +60,8 @@
                     SOURCES alice3TrackingTranslator.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
+
+o2physics_add_dpl_workflow(alice3-dq-table-maker
+                    SOURCES alice3-dq-table-maker.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2::DetectorsBase O2Physics::AnalysisCCDB O2Physics::PWGDQCore
+                    COMPONENT_NAME Analysis)

ALICE3/Tasks/CMakeLists.txt

@@ -89,6 +89,11 @@ o2physics_add_dpl_workflow(alice3-strangeness
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
 
+o2physics_add_dpl_workflow(alice3-dq-efficiency
+                    SOURCES alice3-dq-efficiency.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGDQCore
+                    COMPONENT_NAME Analysis)
+
 o2physics_add_dpl_workflow(alice3-decayer-qa
                     SOURCES alice3DecayerQA.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore

PWGDQ/Core/CutsLibrary.cxx

@@ -3709,6 +3709,92 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
     return cut;
   }
 
+  //---------------------------------------------------------------
+  // ALICE 3 studies composite cuts
+
+  if (!nameStr.compare("alice3StandardKine")) {
+    cut->AddCut(GetAnalysisCut("alice3StandardKine"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3KineSkim")) {
+    cut->AddCut(GetAnalysisCut("alice3KineSkim"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3TrackQuality")) {
+    cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDEl")) {
+    cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDPi")) {
+    cut->AddCut(GetAnalysisCut("alice3iTOFPIDPi"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDKa")) {
+    cut->AddCut(GetAnalysisCut("alice3iTOFPIDKa"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDPr")) {
+    cut->AddCut(GetAnalysisCut("alice3iTOFPIDPr"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDEl")) {
+    cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDPi")) {
+    cut->AddCut(GetAnalysisCut("alice3oTOFPIDPi"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDKa")) {
+    cut->AddCut(GetAnalysisCut("alice3oTOFPIDKa"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDPr")) {
+    cut->AddCut(GetAnalysisCut("alice3oTOFPIDPr"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDEl")) {
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDEl"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDPi")) {
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDPi"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDKa")) {
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDKa"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDPr")) {
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDPr"));
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3DielectronPID")) {
+    cut->AddCut(GetAnalysisCut("alice3StandardKine"));
+    cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDEl"));
+    cut->AddCut(GetAnalysisCut("alice3RICHPIDPiPrRejection"));
+    return cut;
+  }
+
   delete cut;
   LOGF(fatal, Form("Did not find cut %s. Returning nullptr", cutName));
   return nullptr;
@@ -6732,6 +6818,122 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
     return cut;
   }
 
+  //---------------------------------------------------------------
+  // ALICE 3 Analysis Cuts
+
+  if (!nameStr.compare("alice3KineSkim")) {
+    cut->AddCut(VarManager::kPt, 0.05, 1000.0);
+    cut->AddCut(VarManager::kEta, -4.0, 4.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3StandardKine")) {
+    cut->AddCut(VarManager::kPt, 0.1, 1000.0);
+    cut->AddCut(VarManager::kEta, -1.75, 1.75); // Value taken from LoI
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3TrackQuality")) {
+    cut->AddCut(VarManager::kTrackDCAxy, -1.0, 1.0);
+    cut->AddCut(VarManager::kTrackDCAz, -1.0, 1.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3OTPIDEl")) {
+    cut->AddCut(VarManager::kOTnSigmaEl, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3OTPIDPi")) {
+    cut->AddCut(VarManager::kOTnSigmaPi, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3OTPIDKa")) {
+    cut->AddCut(VarManager::kOTnSigmaKa, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3OTPIDPr")) {
+    cut->AddCut(VarManager::kOTnSigmaPr, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDEl")) {
+    cut->AddCut(VarManager::kInnerTOFnSigmaEl, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDPi")) {
+    cut->AddCut(VarManager::kInnerTOFnSigmaPi, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDKa")) {
+    cut->AddCut(VarManager::kInnerTOFnSigmaKa, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3iTOFPIDPr")) {
+    cut->AddCut(VarManager::kInnerTOFnSigmaPr, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDEl")) {
+    cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDEl")) {
+    cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDPi")) {
+    cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDKa")) {
+    cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3oTOFPIDPr")) {
+    cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0, 3.0);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDEl")) {
+    cut->AddCut(VarManager::kRICHnSigmaEl, -3.0, 3.0);
+    cut->AddCut(VarManager::kHasRICHSigEl, 0.5, 1.5);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDPi")) {
+    cut->AddCut(VarManager::kRICHnSigmaPi, -3.0, 3.0);
+    cut->AddCut(VarManager::kHasRICHSigPi, 0.5, 1.5);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDKa")) {
+    cut->AddCut(VarManager::kRICHnSigmaKa, -3.0, 3.0);
+    cut->AddCut(VarManager::kHasRICHSigKa, 0.5, 1.5);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDPr")) {
+    cut->AddCut(VarManager::kRICHnSigmaPr, -3.0, 3.0);
+    cut->AddCut(VarManager::kHasRICHSigPr, 0.5, 1.5);
+    return cut;
+  }
+
+  if (!nameStr.compare("alice3RICHPIDPiPrRejection")) {
+    cut->AddCut(VarManager::kRICHnSigmaPi, -3.0, 3.0, true);
+    cut->AddCut(VarManager::kRICHnSigmaPr, -3.0, 3.0, true);
+    return cut;
+  }
+
   delete cut;
   LOGF(fatal, Form("Did not find cut %s", cutName));
   return nullptr;