FilterCalculator.h

/**************************************************************************
 * 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.                  *
 **************************************************************************/
#pragma once
 
#include <hlt/softwaretrigger/core/SoftwareTriggerObject.h>
#include <hlt/softwaretrigger/calculations/SoftwareTriggerCalculation.h>
#include <hlt/dbobjects/HLTPrefilterParameters.h>
#include <hlt/utilities/HLTPrefilter.h>
#include <framework/database/DBObjPtr.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <mdst/dataobjects/Track.h>
#include <mdst/dataobjects/ECLCluster.h>
#include <mdst/dataobjects/TRGSummary.h>
#include <trg/cdc/Unpacker.h>
#include <framework/gearbox/Unit.h>
 
namespace Belle2::SoftwareTrigger {
  class FilterCalculator : public SoftwareTriggerCalculation {
  public:
    FilterCalculator();

    void requireStoreArrays() override;

    void doCalculation(SoftwareTriggerObject& calculationResult) override;
 
  private:
    StoreArray<Track> m_tracks;
    StoreArray<ECLCluster> m_eclClusters;
    StoreObjPtr<TRGSummary> m_l1Trigger;
    StoreArray<CDCTriggerUnpacker::NNBitStream> m_bitsNN;

    double m_looseTrkZ0 = 10 * Unit::cm;
    double m_tightTrkZ0 = 2 * Unit::cm;
    double m_E2min = 0.2;
    double m_E0min = 0.3;
    double m_Ehigh = 2;
    double m_EminLab = 0.18;
    double m_EminLab4Cluster = 0.3;
    double m_EminLab3Cluster = 0.5;
    double m_EsinglePhoton = 1;
    double m_reducedEsinglePhoton = 0.5;
    double m_cosmicMinPt = 0.5 * Unit::GeV;
    double m_cosmicMaxClusterEnergy = 1.0 * Unit::GeV;
    double m_goodMagneticRegionZ0 = 57.;
    double m_goodMagneticRegionD0 = 26.5;

    DBObjPtr<HLTPrefilterParameters> m_hltPrefilterParameters; 
    HLTPrefilter::TimingCutState m_timingPrefilter;
    HLTPrefilter::CDCECLCutState m_cdceclPrefilter;
 
  };
}