EventT0DQM.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#include <reconstruction/modules/EventT0DQM/EventT0DQM.h>
#include <framework/core/HistoModule.h>
 
using namespace Belle2;
 
REG_MODULE(EventT0DQM);
 
//---------------------------------
EventT0DQMModule::EventT0DQMModule(): HistoModule()
{
  setPropertyFlags(c_ParallelProcessingCertified); // parallel processing
  setDescription("Make data quality monitoring plots for EventT0 for bhabha, mu mu, and hadron samples for different trigger (time) sources.");
}
 
//---------------------------------
EventT0DQMModule::~EventT0DQMModule() { }
 
 
//---------------------------------
void EventT0DQMModule::defineHisto()
{
 
  TDirectory* oldDir = gDirectory;
  oldDir->mkdir("EventT0")->cd();
 
  int nBins = 400 ;
  double minT0 = -100 ;
  double maxT0 =  100 ;
 
  m_histEventT0_ECL_bhabha_L1_ECLTRG = new TH1F("m_histEventT0_ECL_bhabha_L1_ECLTRG",
                                                "ECL EventT0, L1TRG from ECL, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_bhabha_L1_ECLTRG = new TH1F("m_histEventT0_CDC_bhabha_L1_ECLTRG",
                                                "CDC EventT0, L1TRG from ECL, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_bhabha_L1_ECLTRG = new TH1F("m_histEventT0_TOP_bhabha_L1_ECLTRG",
                                                "TOP EventT0, L1TRG from ECL, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_bhabha_L1_ECLTRG = new TH1F("m_histEventT0_SVD_bhabha_L1_ECLTRG",
                                                "SVD EventT0, L1TRG from ECL, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_hadron_L1_ECLTRG = new TH1F("m_histEventT0_ECL_hadron_L1_ECLTRG",
                                                "ECL EventT0, L1TRG from ECL, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_hadron_L1_ECLTRG = new TH1F("m_histEventT0_CDC_hadron_L1_ECLTRG",
                                                "CDC EventT0, L1TRG from ECL, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_hadron_L1_ECLTRG = new TH1F("m_histEventT0_TOP_hadron_L1_ECLTRG",
                                                "TOP EventT0, L1TRG from ECL, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_hadron_L1_ECLTRG = new TH1F("m_histEventT0_SVD_hadron_L1_ECLTRG",
                                                "SVD EventT0, L1TRG from ECL, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_mumu_L1_ECLTRG   = new TH1F("m_histEventT0_ECL_mumu_L1_ECLTRG",
                                                "ECL EventT0, L1TRG from ECL, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_mumu_L1_ECLTRG   = new TH1F("m_histEventT0_CDC_mumu_L1_ECLTRG",
                                                "CDC EventT0, L1TRG from ECL, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_mumu_L1_ECLTRG   = new TH1F("m_histEventT0_TOP_mumu_L1_ECLTRG",
                                                "TOP EventT0, L1TRG from ECL, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_mumu_L1_ECLTRG   = new TH1F("m_histEventT0_SVD_mumu_L1_ECLTRG",
                                                "SVD EventT0, L1TRG from ECL, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
 
  m_histEventT0_ECL_bhabha_L1_CDCTRG = new TH1F("m_histEventT0_ECL_bhabha_L1_CDCTRG",
                                                "ECL EventT0, L1TRG from CDC, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_bhabha_L1_CDCTRG = new TH1F("m_histEventT0_CDC_bhabha_L1_CDCTRG",
                                                "CDC EventT0, L1TRG from CDC, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_bhabha_L1_CDCTRG = new TH1F("m_histEventT0_TOP_bhabha_L1_CDCTRG",
                                                "TOP EventT0, L1TRG from CDC, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_bhabha_L1_CDCTRG = new TH1F("m_histEventT0_SVD_bhabha_L1_CDCTRG",
                                                "SVD EventT0, L1TRG from CDC, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_hadron_L1_CDCTRG = new TH1F("m_histEventT0_ECL_hadron_L1_CDCTRG",
                                                "ECL EventT0, L1TRG from CDC, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_hadron_L1_CDCTRG = new TH1F("m_histEventT0_CDC_hadron_L1_CDCTRG",
                                                "CDC EventT0, L1TRG from CDC, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_hadron_L1_CDCTRG = new TH1F("m_histEventT0_TOP_hadron_L1_CDCTRG",
                                                "TOP EventT0, L1TRG from CDC, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_hadron_L1_CDCTRG = new TH1F("m_histEventT0_SVD_hadron_L1_CDCTRG",
                                                "SVD EventT0, L1TRG from CDC, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_mumu_L1_CDCTRG   = new TH1F("m_histEventT0_ECL_mumu_L1_CDCTRG",
                                                "ECL EventT0, L1TRG from CDC, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_mumu_L1_CDCTRG   = new TH1F("m_histEventT0_CDC_mumu_L1_CDCTRG",
                                                "CDC EventT0, L1TRG from CDC, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_mumu_L1_CDCTRG   = new TH1F("m_histEventT0_TOP_mumu_L1_CDCTRG",
                                                "TOP EventT0, L1TRG from CDC, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_mumu_L1_CDCTRG   = new TH1F("m_histEventT0_SVD_mumu_L1_CDCTRG",
                                                "SVD EventT0, L1TRG from CDC, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
 
  m_histEventT0_ECL_bhabha_L1_TOPTRG = new TH1F("m_histEventT0_ECL_bhabha_L1_TOPTRG",
                                                "ECL EventT0, L1TRG from TOP, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_bhabha_L1_TOPTRG = new TH1F("m_histEventT0_CDC_bhabha_L1_TOPTRG",
                                                "CDC EventT0, L1TRG from TOP, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_bhabha_L1_TOPTRG = new TH1F("m_histEventT0_TOP_bhabha_L1_TOPTRG",
                                                "TOP EventT0, L1TRG from TOP, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_bhabha_L1_TOPTRG = new TH1F("m_histEventT0_SVD_bhabha_L1_TOPTRG",
                                                "SVD EventT0, L1TRG from TOP, HLT Bhabha;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_hadron_L1_TOPTRG = new TH1F("m_histEventT0_ECL_hadron_L1_TOPTRG",
                                                "ECL EventT0, L1TRG from TOP, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_hadron_L1_TOPTRG = new TH1F("m_histEventT0_CDC_hadron_L1_TOPTRG",
                                                "CDC EventT0, L1TRG from TOP, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_hadron_L1_TOPTRG = new TH1F("m_histEventT0_TOP_hadron_L1_TOPTRG",
                                                "TOP EventT0, L1TRG from TOP, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_hadron_L1_TOPTRG = new TH1F("m_histEventT0_SVD_hadron_L1_TOPTRG",
                                                "SVD EventT0, L1TRG from TOP, HLT Hadron;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histEventT0_ECL_mumu_L1_TOPTRG   = new TH1F("m_histEventT0_ECL_mumu_L1_TOPTRG",
                                                "ECL EventT0, L1TRG from TOP, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_CDC_mumu_L1_TOPTRG   = new TH1F("m_histEventT0_CDC_mumu_L1_TOPTRG",
                                                "CDC EventT0, L1TRG from TOP, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_TOP_mumu_L1_TOPTRG   = new TH1F("m_histEventT0_TOP_mumu_L1_TOPTRG",
                                                "TOP EventT0, L1TRG from TOP, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
  m_histEventT0_SVD_mumu_L1_TOPTRG   = new TH1F("m_histEventT0_SVD_mumu_L1_TOPTRG",
                                                "SVD EventT0, L1TRG from TOP, HLT mumu;EventT0 [ns];events / 0.5 ns",
                                                nBins, minT0, maxT0);
 
  m_histAlgorithmSourceFractionsHadronL1ECLTRG =
    new TH1D("AlgorithmSourceFractionsHadronL1ECLTRG",
             "Fraction of events with EventT0 from each algorithm for hadronic events triggerd by ECL;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsHadronL1CDCTRG =
    new TH1D("AlgorithmSourceFractionsHadronL1CDCTRG",
             "Fraction of events with EventT0 from each algorithm for hadronic events triggerd by CDC;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsHadronL1TOPTRG =
    new TH1D("AlgorithmSourceFractionsHadronL1TOPTRG",
             "Fraction of events with EventT0 from each algorithm for hadronic events triggerd by TOP;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsBhaBhaL1ECLTRG =
    new TH1D("AlgorithmSourceFractionsBhaBhaL1ECLTRG",
             "Fraction of events with EventT0 from each algorithm for Bhabha events triggerd by ECL;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsBhaBhaL1CDCTRG =
    new TH1D("AlgorithmSourceFractionsBhaBhaL1CDCTRG",
             "Fraction of events with EventT0 from each algorithm for Bhabha events triggerd by CDC;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsBhaBhaL1TOPTRG =
    new TH1D("AlgorithmSourceFractionsBhaBhaL1TOPTRG",
             "Fraction of events with EventT0 from each algorithm for Bhabha events triggerd by TOP;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsMuMuL1ECLTRG =
    new TH1D("AlgorithmSourceFractionsMuMuL1ECLTRG",
             "Fraction of events with EventT0 from each algorithm for #mu#mu events triggerd by ECL;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsMuMuL1CDCTRG =
    new TH1D("AlgorithmSourceFractionsMuMuL1CDCTRG",
             "Fraction of events with EventT0 from each algorithm for #mu#mu events triggerd by CDC;Algorithm;Fraction",
             6, 0, 6);
  m_histAlgorithmSourceFractionsMuMuL1TOPTRG =
    new TH1D("AlgorithmSourceFractionsMuMuL1TOPTRG",
             "Fraction of events with EventT0 from each algorithm for #mu#mu events triggerd by TOP;Algorithm;Fraction",
             6, 0, 6);
 
  for (uint i = 0; i < 6; i++) {
    m_histAlgorithmSourceFractionsHadronL1ECLTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsHadronL1CDCTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsHadronL1TOPTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsBhaBhaL1ECLTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsBhaBhaL1CDCTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsBhaBhaL1TOPTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsMuMuL1ECLTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsMuMuL1CDCTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
    m_histAlgorithmSourceFractionsMuMuL1TOPTRG->GetXaxis()->SetBinLabel(i + 1, c_eventT0Algorithms[i]);
  }
 
  oldDir->cd();
 
}
 
 
//---------------------------------
void EventT0DQMModule::initialize()
{
 
  m_TrgResult.isOptional();
  m_eventT0.isOptional();
 
  REG_HISTOGRAM
 
}
 
 
 
//---------------------------------
void EventT0DQMModule::beginRun()
{
  if (!m_eventT0.isValid()) {
    B2WARNING("Missing EventT0, EventT0DQM is skipped.");
    return;
  }
 
  if (!m_TrgResult.isValid()) {
    B2WARNING("Missing SoftwareTriggerResult, EventT0DQM is skipped.");
    return;
  }
 
  m_histEventT0_ECL_bhabha_L1_ECLTRG->Reset();
  m_histEventT0_CDC_bhabha_L1_ECLTRG->Reset();
  m_histEventT0_TOP_bhabha_L1_ECLTRG->Reset();
  m_histEventT0_SVD_bhabha_L1_ECLTRG->Reset();
 
  m_histEventT0_ECL_hadron_L1_ECLTRG->Reset();
  m_histEventT0_CDC_hadron_L1_ECLTRG->Reset();
  m_histEventT0_TOP_hadron_L1_ECLTRG->Reset();
  m_histEventT0_SVD_hadron_L1_ECLTRG->Reset();
 
  m_histEventT0_ECL_mumu_L1_ECLTRG->Reset();
  m_histEventT0_CDC_mumu_L1_ECLTRG->Reset();
  m_histEventT0_TOP_mumu_L1_ECLTRG->Reset();
  m_histEventT0_SVD_mumu_L1_ECLTRG->Reset();
 
  m_histEventT0_ECL_bhabha_L1_CDCTRG->Reset();
  m_histEventT0_CDC_bhabha_L1_CDCTRG->Reset();
  m_histEventT0_TOP_bhabha_L1_CDCTRG->Reset();
  m_histEventT0_SVD_bhabha_L1_CDCTRG->Reset();
 
  m_histEventT0_ECL_hadron_L1_CDCTRG->Reset();
  m_histEventT0_CDC_hadron_L1_CDCTRG->Reset();
  m_histEventT0_TOP_hadron_L1_CDCTRG->Reset();
  m_histEventT0_SVD_hadron_L1_CDCTRG->Reset();
 
  m_histEventT0_ECL_mumu_L1_CDCTRG->Reset();
  m_histEventT0_CDC_mumu_L1_CDCTRG->Reset();
  m_histEventT0_TOP_mumu_L1_CDCTRG->Reset();
  m_histEventT0_SVD_mumu_L1_CDCTRG->Reset();
 
  m_histEventT0_ECL_bhabha_L1_TOPTRG->Reset();
  m_histEventT0_CDC_bhabha_L1_TOPTRG->Reset();
  m_histEventT0_TOP_bhabha_L1_TOPTRG->Reset();
  m_histEventT0_SVD_bhabha_L1_TOPTRG->Reset();
 
  m_histEventT0_ECL_hadron_L1_TOPTRG->Reset();
  m_histEventT0_CDC_hadron_L1_TOPTRG->Reset();
  m_histEventT0_TOP_hadron_L1_TOPTRG->Reset();
  m_histEventT0_SVD_hadron_L1_TOPTRG->Reset();
 
  m_histEventT0_ECL_mumu_L1_TOPTRG->Reset();
  m_histEventT0_CDC_mumu_L1_TOPTRG->Reset();
  m_histEventT0_TOP_mumu_L1_TOPTRG->Reset();
  m_histEventT0_SVD_mumu_L1_TOPTRG->Reset();
 
  m_histAlgorithmSourceFractionsHadronL1ECLTRG->Reset();
  m_histAlgorithmSourceFractionsHadronL1CDCTRG->Reset();
  m_histAlgorithmSourceFractionsHadronL1TOPTRG->Reset();
  m_histAlgorithmSourceFractionsBhaBhaL1ECLTRG->Reset();
  m_histAlgorithmSourceFractionsBhaBhaL1CDCTRG->Reset();
  m_histAlgorithmSourceFractionsBhaBhaL1TOPTRG->Reset();
  m_histAlgorithmSourceFractionsMuMuL1ECLTRG->Reset();
  m_histAlgorithmSourceFractionsMuMuL1CDCTRG->Reset();
  m_histAlgorithmSourceFractionsMuMuL1TOPTRG->Reset();
 
}
 
 
//---------------------------------
void EventT0DQMModule::event()
{
  if (!m_objTrgSummary.isValid()) {
    B2WARNING("TRGSummary object not available but required to indicate which detector provided the L1 trigger time");
    return;
  }
 
  // Skip this event if there is no event t0, to avoid crashing other DQM
  if (!m_eventT0.isOptional()) {
    B2WARNING("Missing EventT0, EventT0DQM is skipped.");
    return;
  }
 
  // Determine if there is a valid event t0 to use and then extract the information about it
  if (!m_eventT0.isValid()) {
    return ;
  }
 
  if (!m_TrgResult.isValid()) {
    B2WARNING("SoftwareTriggerResult object not available but require to select bhabha/mumu/hadron events for this module");
    return;
  }
 
  const std::map<std::string, int>& fresults = m_TrgResult->getResults();
  if ((fresults.find("software_trigger_cut&skim&accept_bhabha") == fresults.end())   ||
      (fresults.find("software_trigger_cut&skim&accept_mumu_2trk") == fresults.end())   ||
      (fresults.find("software_trigger_cut&skim&accept_hadron") == fresults.end())) {
    B2WARNING("EventT0DQMModule: Can't find required bhabha or mumu or hadron trigger identifier");
    return;
  }
 
  // for L1 timing source is "ecl trigger" -> TRGSummary::ETimingType is 0
  const bool IsECLL1TriggerSource = (m_objTrgSummary->getTimType() == TRGSummary::ETimingType::TTYP_ECL);
  // for L1 timing source is "top trigger" -> TRGSummary::ETimingType is 1
  const bool IsTOPL1TriggerSource = (m_objTrgSummary->getTimType() == TRGSummary::ETimingType::TTYP_TOP);
  // for L1 timing source is "cdc trigger" -> TRGSummary::ETimingType is 3
  const bool IsCDCL1TriggerSource = (m_objTrgSummary->getTimType() == TRGSummary::ETimingType::TTYP_CDC);
 
  B2DEBUG(20, "IsECLL1TriggerSource = " << IsECLL1TriggerSource);
  B2DEBUG(20, "IsTOPL1TriggerSource = " << IsTOPL1TriggerSource);
  B2DEBUG(20, "IsCDCL1TriggerSource = " << IsCDCL1TriggerSource);
 
  // determine if the event was part of the hadron skim or bhabha skim or mumu skim
  const bool IsEvtAcceptedBhabha = (m_TrgResult->getResult("software_trigger_cut&skim&accept_bhabha") ==
                                    SoftwareTriggerCutResult::c_accept);
  const bool IsEvtAcceptedHadron = (m_TrgResult->getResult("software_trigger_cut&skim&accept_hadron") ==
                                    SoftwareTriggerCutResult::c_accept);
  const bool IsEvtAcceptedMumu = (m_TrgResult->getResult("software_trigger_cut&skim&accept_mumu_2trk") ==
                                  SoftwareTriggerCutResult::c_accept);
 
 
  B2DEBUG(20, "bhabha trigger result = " << static_cast<std::underlying_type<SoftwareTriggerCutResult>::type>
          (m_TrgResult->getResult("software_trigger_cut&skim&accept_bhabha"))) ;
  B2DEBUG(20, "hadron trigger result = " << static_cast<std::underlying_type<SoftwareTriggerCutResult>::type>
          (m_TrgResult->getResult("software_trigger_cut&skim&accept_hadron"))) ;
  B2DEBUG(20, "mu mu trigger result = "  << static_cast<std::underlying_type<SoftwareTriggerCutResult>::type>
          (m_TrgResult->getResult("software_trigger_cut&skim&accept_mumu_2trk"))) ;
  B2DEBUG(20, "bhabha trigger comparison bool = " << IsEvtAcceptedBhabha) ;
  B2DEBUG(20, "hadron trigger comparison bool = " << IsEvtAcceptedHadron) ;
  B2DEBUG(20, "mumu trigger comparison bool = "   << IsEvtAcceptedMumu) ;
 
 
  // Set the different EventT0 values, default is -1000 in case there are no information based on a given detector
  const double eventT0ECL =
    m_eventT0->hasTemporaryEventT0(Const::EDetector::ECL) ? m_eventT0->getBestECLTemporaryEventT0()->eventT0 : -1000;
  const double eventT0CDC =
    m_eventT0->hasTemporaryEventT0(Const::EDetector::CDC) ? m_eventT0->getBestCDCTemporaryEventT0()->eventT0 : -1000;
  const double eventT0TOP =
    m_eventT0->hasTemporaryEventT0(Const::EDetector::TOP) ? m_eventT0->getBestTOPTemporaryEventT0()->eventT0 : -1000;
  const double eventT0SVD =
    m_eventT0->hasTemporaryEventT0(Const::EDetector::SVD) ? m_eventT0->getBestSVDTemporaryEventT0()->eventT0 : -1000;
 
  const auto checkForCDCAlgorithm = [cdcEventT0s = m_eventT0->getTemporaryEventT0s(Const::EDetector::CDC)](
  const std::string & algorithm) {
    for (const auto& evtt0 : cdcEventT0s) {
      if (evtt0.algorithm == algorithm) {
        return true;
      }
    }
    return false;
  };
 
  const bool hasCDCHitBasedEventT0 = checkForCDCAlgorithm("hit based");
  const bool hasCDCFullGridEventT0 = checkForCDCAlgorithm("chi2");
  // We are interested if an EventT0 is set, not whether temporary EventT0s exist that might not be used
  const bool hasAnyEventT0 = m_eventT0->hasEventT0();
  const bool hasECLEventT0 = m_eventT0->hasTemporaryEventT0(Const::EDetector::ECL);
  const bool hasSVDEventT0 = m_eventT0->hasTemporaryEventT0(Const::EDetector::SVD);
  const bool hasTOPEventT0 = m_eventT0->hasTemporaryEventT0(Const::EDetector::TOP);
 
  // Fill the plots that used the ECL trigger as the L1 timing source
  if (IsECLL1TriggerSource) {
    // Fill the histograms with the event t0 values
    if (IsEvtAcceptedBhabha) {   // fill the bha bha skim event t0s
      m_histEventT0_ECL_bhabha_L1_ECLTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_bhabha_L1_ECLTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_bhabha_L1_ECLTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_bhabha_L1_ECLTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsBhaBhaL1ECLTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedHadron) {    // fill the hadron skim event t0s
      m_histEventT0_ECL_hadron_L1_ECLTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_hadron_L1_ECLTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_hadron_L1_ECLTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_hadron_L1_ECLTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsHadronL1ECLTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedMumu) {    // fill the mumu skim event t0s
      m_histEventT0_ECL_mumu_L1_ECLTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_mumu_L1_ECLTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_mumu_L1_ECLTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_mumu_L1_ECLTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsMuMuL1ECLTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
  }
  // Fill the plots that used the TOP trigger as the L1 timing source
  else if (IsTOPL1TriggerSource) {
    // Fill the histograms with the event t0 values
    if (IsEvtAcceptedBhabha) {   // fill the bha bha skim event t0s
      m_histEventT0_ECL_bhabha_L1_TOPTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_bhabha_L1_TOPTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_bhabha_L1_TOPTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_bhabha_L1_TOPTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsBhaBhaL1TOPTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedHadron) {    // fill the hadron skim event t0s
      m_histEventT0_ECL_hadron_L1_TOPTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_hadron_L1_TOPTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_hadron_L1_TOPTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_hadron_L1_TOPTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsHadronL1TOPTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedMumu) {    // fill the mumu skim event t0s
      m_histEventT0_ECL_mumu_L1_TOPTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_mumu_L1_TOPTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_mumu_L1_TOPTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_mumu_L1_TOPTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsMuMuL1TOPTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
  }
  // Fill the plots that used the CDC trigger as the L1 timing source
  else if (IsCDCL1TriggerSource) {
    // Fill the histograms with the event t0 values
    if (IsEvtAcceptedBhabha) {   // fill the bha bha skim event t0s
      m_histEventT0_ECL_bhabha_L1_CDCTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_bhabha_L1_CDCTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_bhabha_L1_CDCTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_bhabha_L1_CDCTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsBhaBhaL1CDCTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedHadron) {    // fill the hadron skim event t0s
      m_histEventT0_ECL_hadron_L1_CDCTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_hadron_L1_CDCTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_hadron_L1_CDCTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_hadron_L1_CDCTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsHadronL1CDCTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
 
    if (IsEvtAcceptedMumu) {    // fill the mumu skim event t0s
      m_histEventT0_ECL_mumu_L1_CDCTRG->Fill(eventT0ECL);
      m_histEventT0_CDC_mumu_L1_CDCTRG->Fill(eventT0CDC);
      m_histEventT0_TOP_mumu_L1_CDCTRG->Fill(eventT0TOP);
      m_histEventT0_SVD_mumu_L1_CDCTRG->Fill(eventT0SVD);
      fillHistogram(m_histAlgorithmSourceFractionsMuMuL1CDCTRG, hasAnyEventT0, hasECLEventT0, hasSVDEventT0,
                    hasCDCHitBasedEventT0, hasCDCFullGridEventT0, hasTOPEventT0);
    }
  }
 
  B2DEBUG(20, "eventT0ECL = " << eventT0ECL << " ns") ;
  B2DEBUG(20, "eventT0CDC = " << eventT0CDC << " ns") ;
  B2DEBUG(20, "eventT0TOP = " << eventT0TOP << " ns") ;
  B2DEBUG(20, "eventT0SVD = " << eventT0SVD << " ns") ;
}
 
void EventT0DQMModule::fillHistogram(TH1D* hist, const bool hasAnyT0, const bool hasECLT0, const bool hasSVDT0,
                                     const bool hasCDCHitT0, const bool hasCDCGridT0, const bool hasTOPT0)
{
  hist->Fill(-1); // counting events for normalisation
  hist->Fill(c_eventT0Algorithms[0], hasECLT0);
  hist->Fill(c_eventT0Algorithms[1], hasSVDT0);
  hist->Fill(c_eventT0Algorithms[2], hasCDCHitT0);
  hist->Fill(c_eventT0Algorithms[3], hasCDCGridT0);
  hist->Fill(c_eventT0Algorithms[4], hasTOPT0);
  hist->Fill(c_eventT0Algorithms[5], hasAnyT0);
}