Implementation of the ALICE3 DQ data model, table reader and efficiency task

Author: axti98

ALICE3/TableProducer/CMakeLists.txt

@@ -60,3 +60,8 @@ o2physics_add_dpl_workflow(alice3-tracking-translator
                     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/TableProducer/alice3-dq-table-maker.cxx

@@ -88,3 +88,8 @@ o2physics_add_dpl_workflow(alice3-strangeness
                     SOURCES alice3-strangeness.cxx
                     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)

ALICE3/Tasks/CMakeLists.txt

@@ -3703,6 +3703,93 @@ 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;
@@ -6702,6 +6789,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;