EventPreSelectAlg.h

Go to the documentation of this file.

#pragma once
 
#include <Gaudi/Property.h>
#include <GaudiKernel/Algorithm.h>
#include <array>
#include <string>
 
#include "EventPreSelect/IDimuPreSelectTool.h"
#include "EventWriter/IEventWriterTool.h"
 
class EventPreSelectAlg : public Algorithm {
public:
  using Algorithm::Algorithm;
 
  StatusCode initialize() override;
  StatusCode execute() override;
  StatusCode finalize() override;
 
private:
  bool WhetherSector( double ph, double ph1, double ph2 );
 
private:
  static std::array<int, 43> idmax;
 
  /**number of total events */
  long m_events{ 0 };
 
  /**number of events selected */
  long m_barrelBhabhaNumber{ 0 };
  long m_endcapBhabhaNumber{ 0 };
  long m_barrelDimuNumber{ 0 };
  long m_endcapDimuNumber{ 0 };
  long m_hadronNumber{ 0 };
  long m_barrelDiphotonNumber{ 0 };
  long m_endcapDiphotonNumber{ 0 };
 
  IDimuPreSelectTool*                    m_dimuTool;
  std::map<std::string, IEventWriterTool*> m_preSelTools;
 
  /* Properties */
  BooleanProperty m_output{ this, "Output", false };
  DoubleProperty  m_ecm{ this, "Ecm", 3.686 };
  BooleanProperty m_selectBhabha{ this, "SelectBhabha", false };
  BooleanProperty m_selectDimu{ this, "SelectDimu", false };
  BooleanProperty m_selectHadron{ this, "SelectHadron", false };
  BooleanProperty m_selectDiphoton{ this, "SelectDiphoton", false };
  BooleanProperty m_writeDst{ this, "WriteDst", false };
  BooleanProperty m_writeRec{ this, "WriteRec", false };
  DoubleProperty  m_vr0cut{ this, "Vr0cut", 1.0 };
  DoubleProperty  m_vz0cut{ this, "Vz0cut", 5.0 };
  DoubleProperty  m_pt0HighCut{ this, "Pt0HighCut", 5.0 };
  DoubleProperty  m_pt0LowCut{ this, "Pt0LowCut", 0.05 };
  DoubleProperty  m_energyThreshold{ this, "EnergyThreshold", 0.05 };
  DoubleProperty  m_gammaPhiCut{ this, "GammaPhiCut", 20.0 };
  DoubleProperty  m_gammaThetaCut{ this, "GammaThetaCut", 20.0 };
 
  DoubleProperty  m_bhabhaEmcECut{ this, "BhabhaEmcECut", 0.7 * m_ecm };
  DoubleProperty  m_bhabhaMaxECut{ this, "BhabhaMaxECut", 0.3 * m_ecm };
  DoubleProperty  m_bhabhaSecECut{ this, "BhabhaSecECut", 0.1 * m_ecm };
  DoubleProperty  m_bhabhaDTheCut{ this, "BhabhaDTheCut", 3. };
  DoubleProperty  m_bhabhaDPhiCut1{ this, "BhabhaDPhiCut1", -25. };
  DoubleProperty  m_bhabhaDPhiCut2{ this, "BhabhaDPhiCut2", -4. };
  DoubleProperty  m_bhabhaDPhiCut3{ this, "BhabhaDPhiCut3", 2. };
  DoubleProperty  m_bhabhaDPhiCut4{ this, "BhabhaDPhiCut4", 20. };
  IntegerProperty m_bhabhaMdcHitCutB{ this, "BhabhaMdcHitCutB", 15 };
  IntegerProperty m_bhabhaMdcHitCutE{ this, "BhabhaMdcHitCutE", 5 };
 
  DoubleProperty m_dimuEHighCut{ this, "DimuEHighCut", 0.27 * m_ecm };
  DoubleProperty m_dimuELowCut{ this, "DimuELowCut", 0.027 * m_ecm };
  DoubleProperty m_dimuDTheCut{ this, "DimuDTheCut", 3. };
  DoubleProperty m_dimuDPhiCut{ this, "DimuDPhiCut", 23. };
 
  DoubleProperty m_hadronChaECut{ this, "HadronChaECut", 0.3 * m_ecm };
  DoubleProperty m_hadronNeuECut{ this, "HadronNeuECut", 0.3 * m_ecm };
 
  DoubleProperty m_diphotonEmcECut{ this, "DiphotonEmcECut", 0.7 * m_ecm };
  DoubleProperty m_diphotonSecECut{ this, "DiphotonSecECut", 0.3 * m_ecm };
  DoubleProperty m_diphotonDTheCut{ this, "DiphotonDTheCut", 3. };
  DoubleProperty m_diphotonDPhiCut1{ this, "DiphotonDPhiCut1", -4. };
  DoubleProperty m_diphotonDPhiCut2{ this, "DiphotonDPhiCut2", 2. };
 
  /* Tuples */
  NTuple::Tuple*       m_tuple0; // Trk2Sh2
  NTuple::Item<double> m_esum;
  NTuple::Item<double> m_eemc;
  NTuple::Item<double> m_etot;
  NTuple::Item<double> m_nGood;
  NTuple::Item<double> m_nCharge;
  NTuple::Item<double> m_nGam;
  NTuple::Item<double> m_ptot;
  NTuple::Item<double> m_pp;
  NTuple::Item<double> m_pm;
  NTuple::Item<long>   m_runnb;
  NTuple::Item<long>   m_evtnb;
  NTuple::Item<double> m_maxE;
  NTuple::Item<double> m_secE;
  NTuple::Item<double> m_dThe;
  NTuple::Item<double> m_dPhi;
  NTuple::Item<long>   m_mdcHit1;
  NTuple::Item<long>   m_mdcHit2;
 
  NTuple::Tuple*       m_tuple1; // charged track vertex
  NTuple::Item<double> m_vx0;
  NTuple::Item<double> m_vy0;
  NTuple::Item<double> m_vz0;
  NTuple::Item<double> m_vr0;
  NTuple::Item<double> m_theta0;
  NTuple::Item<double> m_p0;
  NTuple::Item<double> m_pt0;
 
  NTuple::Tuple*       m_tuple2; // photon
  NTuple::Item<double> m_dthe;
  NTuple::Item<double> m_dphi;
  NTuple::Item<double> m_dang;
  NTuple::Item<double> m_eraw;
 
  NTuple::Tuple* m_tuple3; // for dimu selection
};