A Class of ECL Trigger clustering
More...

#include <TrgEclBhabha.h>

Public Member Functions
	TrgEclBhabha ()
	Constructor.

virtual	~TrgEclBhabha ()
	Constructor.

	TrgEclBhabha (TrgEclBhabha &)=delete
	Destructor.

TrgEclBhabha &	operator= (TrgEclBhabha &)=delete
	Assignement operator, deleted.

bool	GetBhabha00 (std::vector< double >)
	Belle 2D Bhabha veto method.

bool	GetBhabha01 ()
	Belle II 3D Bhabha method for veto.

bool	GetBhabha02 ()
	Belle II 3D Bhabha method for selection.

bool	Getmumu ()
	MuMu selection for calibration.

bool	GetTaub2b (double)
	Taub2b selection.

bool	GetTaub2b2 (double)
	Taub2b selection (tighter selection than Taub2b)

bool	GetTaub2b3 (double)
	Taub2b3.

int	GetBhabhaAddition (void)
	Additional Bhabha veto.

std::vector< double >	GetBhabhaComb ()
	Output 2D Bhabha combination.

void	save (int)
	Save.

void	set2DBhabhaThreshold (const std::vector< double > &i2DBhabhaThresholdFWD, const std::vector< double > &i2DBhabhaThresholdBWD)
	set 2D Bhabha Energy Threshold

void	set3DBhabhaSelectionThreshold (const std::vector< double > &i3DBhabhaSelectionThreshold)
	set 3D selection Bhabha Energy Threshold

void	set3DBhabhaVetoThreshold (const std::vector< double > &i3DBhabhaVetoThreshold)
	set 3D veto Bhabha Energy Threshold

void	set3DBhabhaSelectionAngle (const std::vector< double > &i3DBhabhaSelectionAngle)
	set 3D selection Bhabha Energy Angle

void	set3DBhabhaVetoAngle (const std::vector< double > &i3DBhabhaVetoAngle)
	set 3D veto Bhabha Energy Angle

void	setmumuThreshold (int mumuThreshold)
	set mumu bit Threshold

void	setmumuAngle (const std::vector< double > &imumuAngle)
	set mumu bit Angle selection

void	set3DBhabhaAddAngleCut (const std::vector< double > &i3DBhabhaAddAngleCut)
	set 3D Bhabha addtion Angle selection

void	set3DBhabhaVetoInTrackThetaRegion (const std::vector< int > &i3DBhabhaVetoInTrackThetaRegion)
	set ThetaID (low and high) for 3DBhabhaVetoInTrack

void	setTaub2bAngleCut (const std::vector< int > &iTaub2bAngleCut)
	set 2 Cluster angle selection for tau 1x1 decay [0], [1] for low and high of dphi, [2], [3] for low and high of Theta Sum

void	setTaub2bEtotCut (const double iTaub2bEtotCut)
	set total energy cut for taub2b

void	setTaub2bClusterECut (const double iTaub2bClusterECut1, const double iTaub2bClusterECut2)
	set cluster energy cut for taub2b

void	setTaub2b2Cut (const std::vector< int > &iTaub2b2AngleCut, const double iTaub2b2EtotCut, const double iTaub2b2CLEEndcapCut, const double iTaub2b2CLECut)
	set taub2b2 cut

void	setTaub2b3Cut (const std::vector< int > &iTaub2b3AngleCut, const double iTaub2b3EtotCut, const double iTaub2b3CLEb2bCut, const double iTaub2b3CLELowCut, const double iTaub2b3CLEHighCut)
	set taub2b3 cut

int	get3DBhabhaVetoInTrackFlag (void)
	get trigger bit of flag(1bit) whether two clusters statisfy 3D Bhabha veto are in CDCTRG region in theta (="InTrack") or not

int	get3DBhabhaVetoClusterTCId (int cl_idx)
	get each TCID(most energetic TC in a cluster) of two clusters of 3D Bhabha veto

int	get3DBhabhaVetoClusterThetaId (int cl_idx)
	get each TC theta ID(most energetic TC in a cluster) of two clusters of 3D Bhabha veto

double	get3DBhabhaVetoClusterEnergy (int cl_idx)
	get each cluster energy of two clusters of 3D Bhabha veto (GeV)

double	get3DBhabhaVetoClusterTiming (int cl_idx)
	get each cluster timing of two clusters of 3D Bhabha veto (ns)

int	get3DBhabhaSelectionThetaFlag (void)
	get trigger bit(2bits) of flag which shows theta position of clusters of 3DBhabha Selection.

int	get3DBhabhaSelectionClusterTCId (int cl_idx)
	get each TCID(most energetic TC in a cluster) of two clusters of 3D Bhabha selection

int	get3DBhabhaSelectionClusterThetaId (int cl_idx)
	get each TC theta ID(most energetic TC in a cluster) of two clusters of 3D Bhabha selection

double	get3DBhabhaSelectionClusterEnergy (int cl_idx)
	get each cluster energy of two clusters of 3D Bhabha selection (GeV)

double	get3DBhabhaSelectionClusterTiming (int cl_idx)
	get each cluster timing of two clusters of 3D Bhabha selection (ns)

int	getTaub2bAngleFlag (void)
	get taub2b 2 cluster angle cut flag

int	getTaub2bEtotFlag (void)
	get total energy(TC theta id=1-17) flag for taub2b

int	getTaub2bClusterEFlag (void)
	taub2b Cluster energy flag

Private Member Functions
void	get2CLETP (int, int, int &, int &, int &, int &)
	2 cluster energies, phi difference and theta sum

Private Attributes
TrgEclMapping *	_TCMap = nullptr
	Object of TC Mapping.

TrgEclDataBase *	_database = nullptr
	Object of Trigger ECL DataBase.

std::vector< double >	BhabhaComb
	Bhabha Combination.

std::vector< double >	MaxTCId
	Max TC Id.

std::vector< double >	MaxTCThetaId
	Cluster ThetaId.

std::vector< double >	ClusterEnergy
	Cluster Energy.

std::vector< double >	ClusterTiming
	Cluster Timing.

std::vector< ROOT::Math::XYZVector >	ClusterPosition
	Cluster Timing.

std::vector< double >	_2DBhabhaThresholdFWD
	2D Bhabha Energy Threshold

std::vector< double >	_2DBhabhaThresholdBWD
	2D Bhabha Energy Threshold

std::vector< double >	_3DBhabhaSelectionThreshold
	3D Selection Bhabha Energy Threshold

std::vector< double >	_3DBhabhaVetoThreshold
	3D Veto Bhabha Energy Threshold

std::vector< double >	_3DBhabhaSelectionAngle
	3D Selection Bhabha Energy Angle

std::vector< double >	_3DBhabhaVetoAngle
	3D Veto Bhabha Energy Angle

double	_mumuThreshold
	mumu bit Energy Threshold

std::vector< double >	_mumuAngle
	mumu bit Angle

std::vector< double >	m_3DBhabhaAddAngleCut
	Angle selection of additional Bhabha veto in CM frame.

std::vector< int >	m_taub2bAngleCut
	taub2b 2 Cluster angle cut (degree)

double	m_taub2bEtotCut
	taub2b total energy(TC theta id=1-17) cut (GeV)

double	m_taub2bClusterECut1
	taub2b Cluster one of energy cut in b2b in lab (GeV)

double	m_taub2bClusterECut2
	taub2b Cluster one of energy cut in b2b in lab (GeV)

int	m_taub2bAngleFlag
	taub2b 2 cluster angle cut flag

int	m_taub2bEtotFlag
	taub2b total energy(TC theta id=1-17) flag

int	m_taub2bClusterEFlag
	taub2b Cluster energy flag

double	m_taub2b2EtotCut
	taub2b2 total energy cut (GeV)

std::vector< int >	m_taub2b2AngleCut
	taub2b2 two Cluster angle cut (degree)

double	m_taub2b2CLEEndcapCut
	taub2b2 cluster energy cut for endcap cluster (GeV)

double	m_taub2b2CLECut
	taub2b2 cluster energy cut (GeV)

std::vector< int >	m_taub2b3AngleCut
	taub2b3 selection cuts (3,2,1,0) = (dphi low, dphi high, theta_sum low, theta_sum high)

double	m_taub2b3EtotCut
	taub2b3 total energy (TC theta ID =1-17) (GeV)

double	m_taub2b3CLEb2bCut
	taub2b3 cluster energy cut (GeV)

double	m_taub2b3CLELowCut
	taub2b3 cluster energy cut2 (GeV)

double	m_taub2b3CLEHighCut
	taub2b3 cluster energy cut3 (GeV)

int	m_3DBhabhaVetoInTrackFlag
	trigger bit of flag(1bit) whether two clusters statisfy 3D Bhabha veto are in CDCTRG region in theta (="InTrack") or not

std::vector< int >	m_3DBhabhaVetoClusterTCIds
	TCIDs of two clusters of 3D Bhabha veto.

std::vector< int >	m_3DBhabhaVetoClusterThetaIds
	ThetaIds of two clusters of 3D Bhabha veto.

std::vector< double >	m_3DBhabhaVetoClusterEnergies
	Energies of two clusters of 3D Bhabha veto (GeV)

std::vector< double >	m_3DBhabhaVetoClusterTimings
	Timings of two clusters of 3D Bhabha veto (ns)

std::vector< int >	m_3DBhabhaVetoInTrackThetaRegion
	theta region(low and high) of 3D Bhbabha veto InTrack

int	m_3DBhabhaSelectionThetaFlag = std::numeric_limits<int>::quiet_NaN()
	flag which shows theta position of clusters of 3DBhabha Selection.

std::vector< int >	m_3DBhabhaSelectionClusterTCIds
	TCIDs of two clusters of 3D Bhabha selection.

std::vector< int >	m_3DBhabhaSelectionClusterThetaIds
	ThetaIDs of two clusters used for 3D Bhabha selection.

std::vector< double >	m_3DBhabhaSelectionClusterEnergies
	Energies of two clusters used for 3D Bhabha selection (GeV)

std::vector< double >	m_3DBhabhaSelectionClusterTimings
	Timings of two clusters used for 3D Bhabha selection (ns)

Detailed Description

A Class of ECL Trigger clustering

Definition at line 32 of file TrgEclBhabha.h.

Constructor & Destructor Documentation

◆ TrgEclBhabha() [1/2]

TrgEclBhabha ( )

Constructor.

Definition at line 67 of file TrgEclBhabha.cc.

                          :
  _mumuThreshold(20),
  m_3DBhabhaVetoInTrackThetaRegion(3, 15)
{
  BhabhaComb.clear();
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
 
  _TCMap = new TrgEclMapping();
 
  _database = new TrgEclDataBase();
 
  _2DBhabhaThresholdFWD.clear();
  _2DBhabhaThresholdBWD.clear();
  _3DBhabhaSelectionThreshold.clear();
  _3DBhabhaVetoThreshold.clear();
  _3DBhabhaSelectionAngle.clear();
  _3DBhabhaVetoAngle.clear();
  _mumuAngle.clear();
  m_3DBhabhaAddAngleCut.clear();
 
  _2DBhabhaThresholdFWD = {40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 30, 35}; //  (100 MeV)
  _2DBhabhaThresholdBWD = {25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 30, 30}; // (100 MeV)
  _3DBhabhaVetoThreshold      = {30, 45}; // (low, high) (100 MeV)
  _3DBhabhaSelectionThreshold = {25, 40}; // (low, high) (100 MeV)
  _3DBhabhaVetoAngle      = {160, 200, 165, 190}; // (phi_low, phi_high, theta_low, theta_high) (degree)
  _3DBhabhaSelectionAngle = {140, 220, 160, 200}; // (phi_low, phi_high, theta_low, theta_high) (degree)
  _mumuAngle              = {160, 200, 165, 190}; // (phi_low, phi_high, theta_low, theta_high) (degree)
  m_3DBhabhaAddAngleCut   = {150, 210, 160, 200}; // (phi_low, phi_high, theta_low, theta_high) (degree)
 
  m_3DBhabhaVetoInTrackFlag = -10;
  m_3DBhabhaVetoClusterTCIds.clear();
  m_3DBhabhaVetoClusterEnergies.clear();
  m_3DBhabhaVetoClusterTimings.clear();
  m_3DBhabhaVetoClusterThetaIds.clear();
 
  m_3DBhabhaVetoInTrackFlag = -10;
 
  m_taub2bEtotCut = 7.0;
  m_taub2bClusterECut1 = 1.9;
  m_taub2bClusterECut2 = 999.0;
  m_taub2bAngleFlag = 0;
  m_taub2bEtotFlag = 0;
  m_taub2bClusterEFlag = 0;
 
  m_taub2b2EtotCut = 7.0;
  m_taub2b2CLEEndcapCut = 3.0;
  m_taub2b2CLECut = 1.62;
 
  //Taub2b3 by S.Ito
  m_taub2b3EtotCut = 7.0;
  m_taub2b3CLEb2bCut  = 0.14;
  m_taub2b3CLELowCut  = 0.12;
  m_taub2b3CLEHighCut = 4.5;
 
}

◆ ~TrgEclBhabha()

~TrgEclBhabha ( )

virtual

Constructor.

Destructor

Definition at line 128 of file TrgEclBhabha.cc.

{
  delete _TCMap;
  delete _database;
}

◆ TrgEclBhabha() [2/2]

TrgEclBhabha ( TrgEclBhabha & )

delete

Destructor.

Copy constructor, deleted.

Member Function Documentation

◆ get2CLETP()

void get2CLETP	(	int	TCId1,
		int	TCId2,
		int &	energy1,
		int &	energy2,
		int &	dphi,
		int &	thetaSum
	)

private

2 cluster energies, phi difference and theta sum

Definition at line 929 of file TrgEclBhabha.cc.

{
  int lut1 = _database->Get3DBhabhaLUT(TCId1);
  int lut2 = _database->Get3DBhabhaLUT(TCId2);
  energy1 = 15 & lut1;
  energy2 = 15 & lut2;
  lut1 >>= 4;
  lut2 >>= 4;
  int phi1 = 511 & lut1;
  int phi2 = 511 & lut2;
  lut1 >>= 9;
  lut2 >>= 9;
  int theta1 = lut1;
  int theta2 = lut2;
  dphi = abs(phi1 - phi2);
  if (dphi > 180) {dphi = 360 - dphi;}
  thetaSum = theta1 + theta2;
 
  return;
}

◆ get3DBhabhaSelectionClusterEnergy()

double get3DBhabhaSelectionClusterEnergy ( int cl_idx )

inline

get each cluster energy of two clusters of 3D Bhabha selection (GeV)

Definition at line 220 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaSelectionClusterEnergies.size() != 2) {
        return -10;
      }
      return m_3DBhabhaSelectionClusterEnergies[cl_idx];
    }

◆ get3DBhabhaSelectionClusterTCId()

int get3DBhabhaSelectionClusterTCId ( int cl_idx )

inline

get each TCID(most energetic TC in a cluster) of two clusters of 3D Bhabha selection

Definition at line 203 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaSelectionClusterTCIds.size() != 2) {
        return -10;
      }
      return m_3DBhabhaSelectionClusterTCIds[cl_idx];
    }

◆ get3DBhabhaSelectionClusterThetaId()

int get3DBhabhaSelectionClusterThetaId ( int cl_idx )

inline

get each TC theta ID(most energetic TC in a cluster) of two clusters of 3D Bhabha selection

Definition at line 212 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaSelectionClusterThetaIds.size() != 2) {
        return -10;
      }
      return m_3DBhabhaSelectionClusterThetaIds[cl_idx];
    }

◆ get3DBhabhaSelectionClusterTiming()

double get3DBhabhaSelectionClusterTiming ( int cl_idx )

inline

get each cluster timing of two clusters of 3D Bhabha selection (ns)

Definition at line 228 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaSelectionClusterTimings.size() != 2) {
        return -1000;
      }
      return m_3DBhabhaSelectionClusterTimings[cl_idx];
    }

◆ get3DBhabhaSelectionThetaFlag()

int get3DBhabhaSelectionThetaFlag ( void )

inline

get trigger bit(2bits) of flag which shows theta position of clusters of 3DBhabha Selection.

flag=0 : one of clusters goes to ThetaID=1 flag=1 : one of clusters goes to ThetaID=2 flag=2 : one of clusters goes to ThetaID=3 flag=3 : none of clusters fly to ThetaID=1-3 Based on this flag, pre-scale is applied on GDL to have flat entry of Bhabha event in theta for calibration purpose

Definition at line 198 of file TrgEclBhabha.h.

    {
      return m_3DBhabhaSelectionThetaFlag;
    }

◆ get3DBhabhaVetoClusterEnergy()

double get3DBhabhaVetoClusterEnergy ( int cl_idx )

inline

get each cluster energy of two clusters of 3D Bhabha veto (GeV)

Definition at line 175 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaVetoClusterEnergies.size() != 2) {
        return -10;
      }
      return m_3DBhabhaVetoClusterEnergies[cl_idx];
    }

◆ get3DBhabhaVetoClusterTCId()

int get3DBhabhaVetoClusterTCId ( int cl_idx )

inline

get each TCID(most energetic TC in a cluster) of two clusters of 3D Bhabha veto

Definition at line 158 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaVetoClusterTCIds.size() != 2) {
        return -10;
      }
      return m_3DBhabhaVetoClusterTCIds[cl_idx];
    }

◆ get3DBhabhaVetoClusterThetaId()

int get3DBhabhaVetoClusterThetaId ( int cl_idx )

inline

get each TC theta ID(most energetic TC in a cluster) of two clusters of 3D Bhabha veto

Definition at line 167 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaVetoClusterThetaIds.size() != 2) {
        return -10;
      }
      return m_3DBhabhaVetoClusterThetaIds[cl_idx];
    }

◆ get3DBhabhaVetoClusterTiming()

double get3DBhabhaVetoClusterTiming ( int cl_idx )

inline

get each cluster timing of two clusters of 3D Bhabha veto (ns)

Definition at line 183 of file TrgEclBhabha.h.

    {
      if (cl_idx < 0 || cl_idx > 1 || m_3DBhabhaVetoClusterTimings.size() != 2) {
        return -1000;
      }
      return m_3DBhabhaVetoClusterTimings[cl_idx];
    }

◆ get3DBhabhaVetoInTrackFlag()

int get3DBhabhaVetoInTrackFlag ( void )

inline

get trigger bit of flag(1bit) whether two clusters statisfy 3D Bhabha veto are in CDCTRG region in theta (="InTrack") or not

Definition at line 153 of file TrgEclBhabha.h.

    {
      return m_3DBhabhaVetoInTrackFlag;
    }

◆ GetBhabha00()

bool GetBhabha00 ( std::vector< double > PhiRingSum )

Belle 2D Bhabha veto method.

Definition at line 136 of file TrgEclBhabha.cc.

{
  bool BtoBflag = false;
 
 
  vector<int> k011 = {3, 1, 2, 3 };   // (1)  F1+F2 + F3 + B1+B2
  vector<int> k012 = {2, 16, 17};     // (1)  F1+F2 + F3 + B1+B2
 
  vector<int> k021 = {1,  3};         // (2)  F3 + C12
  vector<int> k022 = {1, 15};         // (2)  F3 + C12
 
  vector<int> k03 = {2,  2, 3};       // (3)  F2 + F3
 
  vector<int> k04 = {1,  4};          // (4)  C1 + backward gap
 
  vector<int> k051 = {1,  4 };        // (5)  C1+C11+C12
  vector<int> k052 = {2, 14, 15};     // (5)  C1+C11+C12
 
  vector<int> k061 = {1,  5};         // (6)  C2+C11+C12
  vector<int> k062 = {2,  14, 15};    // (6)  C2+C11+C12
 
  vector<int> k071 = {2,  4, 5};      // (7)  C1+C2+C11
  vector<int> k072 = {1,  14};        // (7)  C1+C2+C11
 
  vector<int> k081 = {1,  5};         // (8)  C2+C10+C11
  vector<int> k082 = {2, 13, 14};     // (8)  C2+C10+C11
 
  vector<int> k091 = {1,  5 };        // (9)  C2+C9+C10
  vector<int> k092 = {2, 12, 13};     // (9)  C2+C9+C10
 
  vector<int> k101 = {2,  5, 6};      // (10) C2+C3+C10
  vector<int> k102 = {1,  13};        // (10) C2+C3+C10
 
  vector<int> k111 = {2,  5, 6};      // (11) C2+C3+C9
  vector<int> k112 = {1,  12};        // (11) C2+C3+C9
 
  vector<int> k121 = {2,  6, 7};      // (9)  C3+C4+C9
  vector<int> k122 = {1,  12};        // (9)  C3+C4+C9
 
  vector<int> k131 = {2,  6, 7};      // (10) C3+C4+C8
  vector<int> k132 = {1,  11};        // (10) C3+C4+C8
 
  vector<int> k141 = {2,  7, 8};      // (11) C4+C5+C8
  vector<int> k142 = {1,  11};        // (11) C4+C5+C8
 
 
  vector<int> k151 = {1,  8};         // (9)  C5+C7+C8
  vector<int> k152 = {2, 10, 11};     // (9)  C5+C7+C8
 
  vector<int> k161 = {2,  8, 9};      // (10) C5+C6+C7
  vector<int> k162 = {2,  9, 10};     // (10) C5+C6+C7
 
  //vector<int> k17 = {2, 14, 15};      // (11) C11+C12 +forward gap
 
  //vector<int> k18 = {1, 16};          // (11) B1 + forward gap
 
  vector<int> kLOM1 = {2,  2, 3 };    // (1)  F1+F2 + F3 + B1+B2
  vector<int> kLOM2 = {2, 16, 17};    // (1)  F1+F2 + F3 + B1+B2
 
  BhabhaComb.clear();
  BhabhaComb.resize(32, 0.0);
  for (int iii = 1; iii <= k011[0]; iii++) { BhabhaComb[0]  += PhiRingSum[k011[iii] - 1];}
  for (int iii = 1; iii <= k012[0]; iii++) { BhabhaComb[1]  += PhiRingSum[k012[iii] - 1];}
  for (int iii = 1; iii <= k021[0]; iii++) { BhabhaComb[2]  += PhiRingSum[k021[iii] - 1]; }
  for (int iii = 1; iii <= k022[0]; iii++) { BhabhaComb[3]  += PhiRingSum[k022[iii] - 1]; }
  for (int iii = 1; iii <= k03[0]; iii++) { BhabhaComb[4]  += PhiRingSum[k03[iii] - 1]; }
  for (int iii = 1; iii <= k04[0]; iii++) { BhabhaComb[5]  += PhiRingSum[k04[iii] - 1]; }
  for (int iii = 1; iii <= k051[0]; iii++) { BhabhaComb[6]  += PhiRingSum[k051[iii] - 1]; }
  for (int iii = 1; iii <= k052[0]; iii++) { BhabhaComb[7]  += PhiRingSum[k052[iii] - 1]; }
  for (int iii = 1; iii <= k061[0]; iii++) { BhabhaComb[8]  += PhiRingSum[k061[iii] - 1]; }
  for (int iii = 1; iii <= k062[0]; iii++) { BhabhaComb[9]  += PhiRingSum[k062[iii] - 1]; }
  for (int iii = 1; iii <= k071[0]; iii++) { BhabhaComb[10]  += PhiRingSum[k071[iii] - 1]; }
  for (int iii = 1; iii <= k072[0]; iii++) { BhabhaComb[11]  += PhiRingSum[k072[iii] - 1]; }
  for (int iii = 1; iii <= k081[0]; iii++) { BhabhaComb[12]  += PhiRingSum[k081[iii] - 1]; }
  for (int iii = 1; iii <= k082[0]; iii++) { BhabhaComb[13]  += PhiRingSum[k082[iii] - 1]; }
  for (int iii = 1; iii <= k091[0]; iii++) { BhabhaComb[14]  += PhiRingSum[k091[iii] - 1]; }
  for (int iii = 1; iii <= k092[0]; iii++) { BhabhaComb[15]  += PhiRingSum[k092[iii] - 1]; }
  for (int iii = 1; iii <= k101[0]; iii++) { BhabhaComb[16]  += PhiRingSum[k101[iii] - 1]; }
  for (int iii = 1; iii <= k102[0]; iii++) { BhabhaComb[17]  += PhiRingSum[k102[iii] - 1]; }
  for (int iii = 1; iii <= k111[0]; iii++) { BhabhaComb[18] += PhiRingSum[k111[iii] - 1]; }
  for (int iii = 1; iii <= k112[0]; iii++) { BhabhaComb[19] += PhiRingSum[k112[iii] - 1]; }
  for (int iii = 1; iii <= k121[0]; iii++) { BhabhaComb[20] += PhiRingSum[k121[iii] - 1]; }
  for (int iii = 1; iii <= k122[0]; iii++) { BhabhaComb[21] += PhiRingSum[k122[iii] - 1]; }
  for (int iii = 1; iii <= k131[0]; iii++) { BhabhaComb[22] += PhiRingSum[k131[iii] - 1]; }
  for (int iii = 1; iii <= k132[0]; iii++) { BhabhaComb[23] += PhiRingSum[k132[iii] - 1]; }
  for (int iii = 1; iii <= k141[0]; iii++) { BhabhaComb[24] += PhiRingSum[k141[iii] - 1]; }
  for (int iii = 1; iii <= k142[0]; iii++) { BhabhaComb[25] += PhiRingSum[k142[iii] - 1]; }
  for (int iii = 1; iii <= k151[0]; iii++) { BhabhaComb[26] += PhiRingSum[k151[iii] - 1]; }
  for (int iii = 1; iii <= k152[0]; iii++) { BhabhaComb[27] += PhiRingSum[k152[iii] - 1]; }
  for (int iii = 1; iii <= k161[0]; iii++) { BhabhaComb[28] += PhiRingSum[k161[iii] - 1]; }
  for (int iii = 1; iii <= k162[0]; iii++) { BhabhaComb[29] += PhiRingSum[k162[iii] - 1]; }
  for (int iii = 1; iii <= kLOM1[0]; iii++) { BhabhaComb[30]  += PhiRingSum[kLOM1[iii] - 1];}
  for (int iii = 1; iii <= kLOM2[0]; iii++) { BhabhaComb[31]  += PhiRingSum[kLOM2[iii] - 1];}
 
 
  BtoBflag  =
    ((BhabhaComb[0]  * 10  >= _2DBhabhaThresholdFWD[0] &&
      BhabhaComb[1]  * 10  >= _2DBhabhaThresholdBWD[0]) ||
     (BhabhaComb[2]  * 10  >= _2DBhabhaThresholdFWD[1] &&
      BhabhaComb[3]  * 10  >= _2DBhabhaThresholdBWD[1]) ||
     (BhabhaComb[6]  * 10  >= _2DBhabhaThresholdFWD[2] &&
      BhabhaComb[7]  * 10  >= _2DBhabhaThresholdBWD[2]) ||
     (BhabhaComb[8]  * 10  >= _2DBhabhaThresholdFWD[3] &&
      BhabhaComb[9]  * 10  >= _2DBhabhaThresholdBWD[3]) ||
     (BhabhaComb[10] * 10  >= _2DBhabhaThresholdFWD[4] &&
      BhabhaComb[11] * 10  >= _2DBhabhaThresholdBWD[4]) ||
     (BhabhaComb[12] * 10  >= _2DBhabhaThresholdFWD[5] &&
      BhabhaComb[13] * 10  >= _2DBhabhaThresholdBWD[5]) ||
     (BhabhaComb[14] * 10  >= _2DBhabhaThresholdFWD[6] &&
      BhabhaComb[15] * 10  >= _2DBhabhaThresholdBWD[6]) ||
     (BhabhaComb[16] * 10  >= _2DBhabhaThresholdFWD[7] &&
      BhabhaComb[17] * 10  >= _2DBhabhaThresholdBWD[7]) ||
     (BhabhaComb[18] * 10  >= _2DBhabhaThresholdFWD[8] &&
      BhabhaComb[19] * 10  >= _2DBhabhaThresholdBWD[8]) ||
     (BhabhaComb[20] * 10  >= _2DBhabhaThresholdFWD[9] &&
      BhabhaComb[21] * 10  >= _2DBhabhaThresholdBWD[9]) ||
     (BhabhaComb[22] * 10  >= _2DBhabhaThresholdFWD[10] &&
      BhabhaComb[23] * 10  >= _2DBhabhaThresholdBWD[10]) ||
     (BhabhaComb[24] * 10  >= _2DBhabhaThresholdFWD[11] &&
      BhabhaComb[25] * 10  >= _2DBhabhaThresholdBWD[11]) ||
     (BhabhaComb[26] * 10  >= _2DBhabhaThresholdFWD[12] &&
      BhabhaComb[27] * 10  >= _2DBhabhaThresholdBWD[12]) ||
     (BhabhaComb[28] * 10  >= _2DBhabhaThresholdFWD[13] &&
      BhabhaComb[29] * 10  >= _2DBhabhaThresholdBWD[13]));
 
  int bhabha01 = 0;
  int bhabha02 = 0;
  int bhabha03 = 0;
  int bhabha04 = 0;
  int bhabha05 = 0;
  int bhabha06 = 0;
  int bhabha07 = 0;
  int bhabha08 = 0;
  int bhabha09 = 0;
  int bhabha10 = 0;
  int bhabha11 = 0;
  int bhabha12 = 0;
  int bhabha13 = 0;
  int bhabha14 = 0;
 
  if ((BhabhaComb[0]  * 10  >= _2DBhabhaThresholdFWD[0] &&
       BhabhaComb[1]  * 10  >= _2DBhabhaThresholdBWD[0]))  {bhabha01 = 1;}
  if ((BhabhaComb[2]  * 10  >= _2DBhabhaThresholdFWD[1] &&
       BhabhaComb[3]  * 10  >= _2DBhabhaThresholdBWD[1]))  {bhabha02 = 1;}
  if ((BhabhaComb[6]  * 10  >= _2DBhabhaThresholdFWD[2] &&
       BhabhaComb[7]  * 10  >= _2DBhabhaThresholdBWD[2]))  {bhabha03 = 1;}
  if ((BhabhaComb[8]  * 10  >= _2DBhabhaThresholdFWD[3] &&
       BhabhaComb[9]  * 10  >= _2DBhabhaThresholdBWD[3]))  {bhabha04 = 1;}
  if ((BhabhaComb[10] * 10  >= _2DBhabhaThresholdFWD[4] &&
       BhabhaComb[11] * 10  >= _2DBhabhaThresholdBWD[4]))  {bhabha05 = 1;}
  if ((BhabhaComb[12] * 10  >= _2DBhabhaThresholdFWD[5] &&
       BhabhaComb[13] * 10  >= _2DBhabhaThresholdBWD[5]))  {bhabha06 = 1;}
  if ((BhabhaComb[14] * 10  >= _2DBhabhaThresholdFWD[6] &&
       BhabhaComb[15] * 10  >= _2DBhabhaThresholdBWD[6]))  {bhabha07 = 1;}
  if ((BhabhaComb[16] * 10  >= _2DBhabhaThresholdFWD[7] &&
       BhabhaComb[17] * 10  >= _2DBhabhaThresholdBWD[7]))  {bhabha08 = 1;}
  if ((BhabhaComb[18] * 10  >= _2DBhabhaThresholdFWD[8] &&
       BhabhaComb[19] * 10  >= _2DBhabhaThresholdBWD[8]))  {bhabha09 = 1;}
  if ((BhabhaComb[20] * 10  >= _2DBhabhaThresholdFWD[9] &&
       BhabhaComb[21] * 10  >= _2DBhabhaThresholdBWD[9]))  {bhabha10 = 1;}
  if ((BhabhaComb[22] * 10  >= _2DBhabhaThresholdFWD[10] &&
       BhabhaComb[23] * 10  >= _2DBhabhaThresholdBWD[10])) {bhabha11 = 1;}
  if ((BhabhaComb[24] * 10  >= _2DBhabhaThresholdFWD[11] &&
       BhabhaComb[25] * 10  >= _2DBhabhaThresholdBWD[11])) {bhabha12 = 1;}
  if ((BhabhaComb[26] * 10  >= _2DBhabhaThresholdFWD[12] &&
       BhabhaComb[27] * 10  >= _2DBhabhaThresholdBWD[12])) {bhabha13 = 1;}
  if ((BhabhaComb[28] * 10  >= _2DBhabhaThresholdFWD[13] &&
       BhabhaComb[29] * 10  >= _2DBhabhaThresholdBWD[13])) {bhabha14 = 1;}
 
  BhabhaComb.clear();
  BhabhaComb.push_back(bhabha01);
  BhabhaComb.push_back(bhabha02);
  BhabhaComb.push_back(bhabha03);
  BhabhaComb.push_back(bhabha04);
  BhabhaComb.push_back(bhabha05);
  BhabhaComb.push_back(bhabha06);
  BhabhaComb.push_back(bhabha07);
  BhabhaComb.push_back(bhabha08);
  BhabhaComb.push_back(bhabha09);
  BhabhaComb.push_back(bhabha10);
  BhabhaComb.push_back(bhabha11);
  BhabhaComb.push_back(bhabha12);
  BhabhaComb.push_back(bhabha13);
  BhabhaComb.push_back(bhabha14);
  BhabhaComb.push_back(0);
  BhabhaComb.push_back(0);
  BhabhaComb.push_back(0);
  BhabhaComb.push_back(0);
 
  return  BtoBflag;
}

◆ GetBhabha01()

bool GetBhabha01 ( )

Belle II 3D Bhabha method for veto.

Definition at line 330 of file TrgEclBhabha.cc.

{
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
  m_3DBhabhaVetoInTrackFlag = -1;
  m_3DBhabhaVetoClusterTCIds.clear();
  m_3DBhabhaVetoClusterEnergies.clear();
  m_3DBhabhaVetoClusterTimings.clear();
  m_3DBhabhaVetoClusterThetaIds.clear();
  //  int EventId = 0;
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    //  EventId = aTRGECLCluster->getEventId();
    int maxTCId    = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    double clustertiming  =  aTRGECLCluster->getTimeAve();
    ROOT::Math::XYZVector clusterposition(aTRGECLCluster->getPositionX(),
                                          aTRGECLCluster->getPositionY(),
                                          aTRGECLCluster->getPositionZ());
    ClusterTiming.push_back(clustertiming);
    ClusterEnergy.push_back(clusterenergy);
    ClusterPosition.push_back(clusterposition);
    MaxTCId.push_back(maxTCId);
  }
  const int ncluster = ClusterEnergy.size();
  //
  //
  //
  BhabhaComb.clear();
  BhabhaComb.resize(18, 0);
 
  int cl_idx1 = -1;
  int cl_idx2 = -1;
  bool BhabhaFlag = false;
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; jcluster ++) {
      bool BtoBFlag = false;
 
      if (icluster == jcluster) {continue;}
 
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
 
      if (dphi > _3DBhabhaVetoAngle[0] &&
          thetaSum > _3DBhabhaVetoAngle[2] &&
          thetaSum < _3DBhabhaVetoAngle[3]) {BtoBFlag = true;}
      if ((ClusterEnergy[icluster] * 100. > _3DBhabhaVetoThreshold[0] * energy1 &&
           ClusterEnergy[jcluster] * 100. > _3DBhabhaVetoThreshold[0] * energy2) &&
          (ClusterEnergy[icluster] * 100. > _3DBhabhaVetoThreshold[1] * energy1 ||
           ClusterEnergy[jcluster] * 100. > _3DBhabhaVetoThreshold[1] * energy2)) {
        if (BtoBFlag) {
          BhabhaFlag = true;
          cl_idx1 = icluster;
          cl_idx2 = jcluster;
        }
      }
    }
  }
 
  if (BhabhaFlag) {
    m_3DBhabhaVetoClusterTCIds.push_back(MaxTCId[cl_idx1]);
    m_3DBhabhaVetoClusterTCIds.push_back(MaxTCId[cl_idx2]);
    m_3DBhabhaVetoClusterEnergies.push_back(ClusterEnergy[cl_idx1]);
    m_3DBhabhaVetoClusterEnergies.push_back(ClusterEnergy[cl_idx2]);
    m_3DBhabhaVetoClusterTimings.push_back(ClusterTiming[cl_idx1]);
    m_3DBhabhaVetoClusterTimings.push_back(ClusterTiming[cl_idx2]);
    int cl_thetaid1 = _TCMap->getTCThetaIdFromTCId(MaxTCId[cl_idx1]);
    int cl_thetaid2 = _TCMap->getTCThetaIdFromTCId(MaxTCId[cl_idx2]);
    m_3DBhabhaVetoClusterThetaIds.push_back(cl_thetaid1);
    m_3DBhabhaVetoClusterThetaIds.push_back(cl_thetaid2);
    // set InTrack flag
    if (cl_thetaid1 >= m_3DBhabhaVetoInTrackThetaRegion[0] &&
        cl_thetaid1 <= m_3DBhabhaVetoInTrackThetaRegion[1] &&
        cl_thetaid2 >= m_3DBhabhaVetoInTrackThetaRegion[0] &&
        cl_thetaid2 <= m_3DBhabhaVetoInTrackThetaRegion[1]) {
      m_3DBhabhaVetoInTrackFlag = 1;
    } else {
      m_3DBhabhaVetoInTrackFlag = 0;
    }
  }
 
  return BhabhaFlag;
}

◆ GetBhabha02()

bool GetBhabha02 ( )

Belle II 3D Bhabha method for selection.

Definition at line 430 of file TrgEclBhabha.cc.

{
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
  m_3DBhabhaSelectionThetaFlag = -1;
  //  int EventId = 0;
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    //  EventId = aTRGECLCluster->getEventId();
    int maxTCId    = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    double clustertiming  =  aTRGECLCluster->getTimeAve();
    ROOT::Math::XYZVector clusterposition(aTRGECLCluster->getPositionX(),
                                          aTRGECLCluster->getPositionY(),
                                          aTRGECLCluster->getPositionZ());
    ClusterTiming.push_back(clustertiming);
    ClusterEnergy.push_back(clusterenergy);
    ClusterPosition.push_back(clusterposition);
    MaxTCId.push_back(maxTCId);
  }
  const int ncluster = ClusterEnergy.size();
  //
  //
  //
  BhabhaComb.clear();
  BhabhaComb.resize(18, 0);
 
  int cl_idx1 = -1;
  int cl_idx2 = -1;
  bool BhabhaFlag = false;
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; ++jcluster) {
      bool BtoBFlag = false;
 
      if (icluster == jcluster) {continue;}
 
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
 
      if (dphi > _3DBhabhaSelectionAngle[0] &&
          dphi < _3DBhabhaSelectionAngle[1] &&
          thetaSum > _3DBhabhaSelectionAngle[2] &&
          thetaSum < _3DBhabhaSelectionAngle[3]) {BtoBFlag = true;}
      if ((ClusterEnergy[icluster] * 100. > _3DBhabhaSelectionThreshold[0] * energy1 &&
           ClusterEnergy[jcluster] * 100. > _3DBhabhaSelectionThreshold[0] * energy2) &&
          (ClusterEnergy[icluster] * 100. > _3DBhabhaSelectionThreshold[1] * energy1 ||
           ClusterEnergy[jcluster] * 100. > _3DBhabhaSelectionThreshold[1] * energy2)) {
        if (BtoBFlag) {
          BhabhaFlag = true;
          cl_idx1 = icluster;
          cl_idx2 = jcluster;
        }
      }
    }
  }
  if (BhabhaFlag) {
    m_3DBhabhaSelectionClusterTCIds.push_back(MaxTCId[cl_idx1]);
    m_3DBhabhaSelectionClusterTCIds.push_back(MaxTCId[cl_idx2]);
    m_3DBhabhaSelectionClusterEnergies.push_back(ClusterEnergy[cl_idx1]);
    m_3DBhabhaSelectionClusterEnergies.push_back(ClusterEnergy[cl_idx2]);
    m_3DBhabhaSelectionClusterTimings.push_back(ClusterTiming[cl_idx1]);
    m_3DBhabhaSelectionClusterTimings.push_back(ClusterTiming[cl_idx2]);
    int cl_thetaid1 = _TCMap->getTCThetaIdFromTCId(MaxTCId[cl_idx1]);
    int cl_thetaid2 = _TCMap->getTCThetaIdFromTCId(MaxTCId[cl_idx2]);
    m_3DBhabhaSelectionClusterThetaIds.push_back(cl_thetaid1);
    m_3DBhabhaSelectionClusterThetaIds.push_back(cl_thetaid2);
    // set theta flag(2bits) for prescale in GDL
    int cl_thetaid0 = (cl_thetaid1 < cl_thetaid2) ? cl_thetaid1 : cl_thetaid2;
    if (cl_thetaid0 <= 0 || cl_thetaid0 >= 18) {
      m_3DBhabhaSelectionThetaFlag = -2;
    } else {
      if (cl_thetaid0 == 1) { m_3DBhabhaSelectionThetaFlag = 0; }
      else if (cl_thetaid0 == 2) { m_3DBhabhaSelectionThetaFlag = 1; }
      else if (cl_thetaid0 == 3) { m_3DBhabhaSelectionThetaFlag = 2; }
      else                    { m_3DBhabhaSelectionThetaFlag = 3; }
    }
  }
 
  return BhabhaFlag;
}

◆ GetBhabhaAddition()

int GetBhabhaAddition ( void )

Additional Bhabha veto.

Definition at line 849 of file TrgEclBhabha.cc.

{
 
  std::vector<int> MaxThetaId;
  MaxThetaId.clear();
  ClusterEnergy.clear();
 
  int bit_bhabha_addition = 0;
 
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    int maxTCId = aTRGECLCluster->getMaxTCId();
    MaxTCId.push_back(maxTCId);
    int maxThetaId = aTRGECLCluster->getMaxThetaId();
    MaxThetaId.push_back(maxThetaId);
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    ClusterEnergy.push_back(clusterenergy);
  }
  int NofCluster1to17 = MaxThetaId.size();
 
  if (NofCluster1to17 == 1) {
    if (MaxThetaId[0] <= 3) {
      bit_bhabha_addition |= 0x01;
    }
  } else if (NofCluster1to17 == 2) {
 
    int energy1 = 0;
    int energy2 = 0;
    int dphi = 0;
    int thetaSum = 0;
    get2CLETP(MaxTCId[0],
              MaxTCId[1],
              energy1,
              energy2,
              dphi,
              thetaSum);
 
    // for hie1
    if ((dphi     > m_3DBhabhaAddAngleCut[0] &&
         dphi     < m_3DBhabhaAddAngleCut[1]) &&
        (thetaSum > m_3DBhabhaAddAngleCut[2] &&
         thetaSum < m_3DBhabhaAddAngleCut[3])) {
      bit_bhabha_addition |= 0x02;
    }
    // for hie2
    if ((dphi     > m_3DBhabhaAddAngleCut[0] &&
         dphi     < m_3DBhabhaAddAngleCut[1]) ||
        (thetaSum > m_3DBhabhaAddAngleCut[2] &&
         thetaSum < m_3DBhabhaAddAngleCut[3])) {
      bit_bhabha_addition |= 0x04;
    }
 
    // for hie3 and hie4
    // get thetaID of lower energy cluster
    int lowe_MaxThetaId = 0;
    if (ClusterEnergy[0] < ClusterEnergy[1]) {
      lowe_MaxThetaId = MaxThetaId[0];
    } else {
      lowe_MaxThetaId = MaxThetaId[1];
    }
    // select cluster in endcap region
    if (lowe_MaxThetaId <= 3 ||
        lowe_MaxThetaId >= 16) {
      // for hie3
      bit_bhabha_addition |= 0x08;
      // for hie4
      // require low energy bound for cluster
      if (ClusterEnergy[0] > 0.5 &&
          ClusterEnergy[1] > 0.5) {
        bit_bhabha_addition |= 0x10;
      }
    }
  }
 
  return bit_bhabha_addition;
}

◆ GetBhabhaComb()

std::vector< double > GetBhabhaComb ( )

inline

Output 2D Bhabha combination.

Definition at line 65 of file TrgEclBhabha.h.

65{return BhabhaComb ;}

◆ Getmumu()

bool Getmumu ( )

MuMu selection for calibration.

Definition at line 527 of file TrgEclBhabha.cc.

{
  bool BtoBFlag = false;
  bool BhabhaFlag = false;
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
  //  int EventId = 0;
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    //  EventId = aTRGECLCluster->getEventId();
    int maxTCId    = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    double clustertiming  =  aTRGECLCluster->getTimeAve();
    ROOT::Math::XYZVector clusterposition(aTRGECLCluster->getPositionX(),
                                          aTRGECLCluster->getPositionY(),
                                          aTRGECLCluster->getPositionZ());
    ClusterTiming.push_back(clustertiming);
    ClusterEnergy.push_back(clusterenergy);
    ClusterPosition.push_back(clusterposition);
    MaxTCId.push_back(maxTCId);
  }
  const int ncluster = ClusterEnergy.size();
  //
  //
  //
  BhabhaComb.clear();
  BhabhaComb.resize(18, 0);
 
  BtoBFlag = false;
 
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; jcluster ++) {
 
      if (icluster == jcluster) {continue;}
 
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
 
      if (dphi > _mumuAngle[0] &&
          dphi < _mumuAngle[1] &&
          thetaSum > _mumuAngle[2] &&
          thetaSum < _mumuAngle[3]) {BtoBFlag = true;}
      if (ClusterEnergy[icluster] * 10. < _mumuThreshold &&
          ClusterEnergy[jcluster] * 10. < _mumuThreshold) {
        if (BtoBFlag) {BhabhaFlag = true;}
      }
    }
  }
  return BhabhaFlag;
}

◆ GetTaub2b()

bool GetTaub2b ( double E_total1to17 )

Taub2b selection.

Definition at line 594 of file TrgEclBhabha.cc.

{
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
 
  m_taub2bAngleFlag    = 0;
  m_taub2bEtotFlag     = 0;
  m_taub2bClusterEFlag = 0;
 
  if (E_total1to17 < m_taub2bEtotCut) {
    m_taub2bEtotFlag = 1;
  }
 
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    int maxTCId    = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    ClusterEnergy.push_back(clusterenergy);
    MaxTCId.push_back(maxTCId);
  }
 
  const int ncluster = ClusterEnergy.size();
 
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; jcluster ++) {
 
      if (icluster == jcluster) {continue;}
 
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
 
      if (dphi     > m_taub2bAngleCut[0] &&
          dphi     < m_taub2bAngleCut[1] &&
          thetaSum > m_taub2bAngleCut[2] &&
          thetaSum < m_taub2bAngleCut[3]) {
        m_taub2bAngleFlag++;
        //
        if ((ClusterEnergy[icluster] < m_taub2bClusterECut1 &&
             ClusterEnergy[jcluster] < m_taub2bClusterECut2) ||
            (ClusterEnergy[icluster] < m_taub2bClusterECut2 &&
             ClusterEnergy[jcluster] < m_taub2bClusterECut1)) {
          m_taub2bClusterEFlag++;
        }
      }
    }
  }
 
  bool Taub2bFlag = false;
  if (m_taub2bAngleFlag    > 0 &&
      m_taub2bEtotFlag     > 0 &&
      m_taub2bClusterEFlag > 0) {
    Taub2bFlag = true;
  }
 
  return Taub2bFlag;
}

◆ GetTaub2b2()

bool GetTaub2b2 ( double E_total1to17 )

Taub2b selection (tighter selection than Taub2b)

Definition at line 667 of file TrgEclBhabha.cc.

{
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
 
  int taub2b2EtotFlag  = 0;
  int taub2b2AngleFlag = 0;
  int taub2b2NCLEndcapFlag = 0;
  int taub2b2AngleCLEFlag = 0;
 
  if (E_total1to17 < m_taub2b2EtotCut) {
    taub2b2EtotFlag = 1;
  }
 
  // countor of cluster for E(cluster)>threshold in endcap
  int ncl_clecut_endcap = 0;
  // cluster array loop
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    int maxTCId           = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    double clusterthetaid = aTRGECLCluster->getMaxThetaId();
    ClusterEnergy.push_back(clusterenergy);
    MaxTCId.push_back(maxTCId);
    // count N(cluster) of E(cluster) > threshold in endcap
    if ((clusterthetaid <=  3 ||
         clusterthetaid >= 16) &&
        clusterenergy > m_taub2b2CLEEndcapCut) {
      ncl_clecut_endcap++;
    }
  }
 
  // bool for N(cluster) of E(cluster)>threshold in endcap
  if (ncl_clecut_endcap < 2) {
    taub2b2NCLEndcapFlag = 1;
  }
 
  // total number of cluster
  const int ncluster = ClusterEnergy.size();
 
  // 2 cluster combination
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; jcluster ++) {
      if (icluster == jcluster) {continue;}
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
      // delta phi and theta sum selection
      if (dphi     > m_taub2b2AngleCut[0] &&
          dphi     < m_taub2b2AngleCut[1] &&
          thetaSum > m_taub2b2AngleCut[2] &&
          thetaSum < m_taub2b2AngleCut[3]) {
        taub2b2AngleFlag++;
        if (ClusterEnergy[icluster] > m_taub2b2CLECut &&
            ClusterEnergy[jcluster] > m_taub2b2CLECut) {
          taub2b2AngleCLEFlag++;
        }
      }
    }
  }
  //
  bool taub2b2Flag = false;
  if (taub2b2EtotFlag      > 0 &&
      taub2b2AngleFlag     > 0 &&
      taub2b2NCLEndcapFlag > 0 &&
      taub2b2AngleCLEFlag  > 0) {
    taub2b2Flag = true;
  }
 
  return taub2b2Flag;
}

◆ GetTaub2b3()

bool GetTaub2b3 ( double E_total1to17 )

Taub2b3.

Definition at line 756 of file TrgEclBhabha.cc.

{
  //
  // Read Cluster Table
  //
  MaxTCId.clear();
  MaxTCThetaId.clear();
  ClusterEnergy.clear();
  ClusterTiming.clear();
  ClusterPosition.clear();
 
  int taub2b3EtotFlag = 0;
  int taub2b3AngleCLEThetaIdFlag = 0;
  int taub2b3CLELowCutFlag  = 1;
  int taub2b3CLEHighCutFlag = 1;
 
  // Total Energy Etot < 7 GeV in lab
  if (E_total1to17 < m_taub2b3EtotCut) {
    taub2b3EtotFlag = 1;
  }
 
  // cluster array loop
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    int maxTCId           = aTRGECLCluster->getMaxTCId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    ClusterEnergy.push_back(clusterenergy);
    MaxTCId.push_back(maxTCId);
    MaxTCThetaId.push_back(aTRGECLCluster->getMaxThetaId());
 
    // All clusters in the event shoule be E > 0.12 GeV in lab.
    if (clusterenergy <= m_taub2b3CLELowCut) {
      taub2b3CLELowCutFlag = 0;
    }
    // The number of clusters with E > 4.5 GeV in lab should be 0.
    if (clusterenergy > m_taub2b3CLEHighCut) {
      taub2b3CLEHighCutFlag = 0;
    }
  }
 
  // total number of cluster
  const int ncluster = ClusterEnergy.size();
 
  // 2 cluster combination
  for (int icluster = 0; icluster < ncluster ; icluster++) {
    for (int jcluster = icluster + 1; jcluster < ncluster; jcluster ++) {
      if (icluster == jcluster) {continue;}
      int energy1 = 0;
      int energy2 = 0;
      int dphi = 0;
      int thetaSum = 0;
      get2CLETP(MaxTCId[icluster],
                MaxTCId[jcluster],
                energy1,
                energy2,
                dphi,
                thetaSum);
 
      // delta phi and theta sum selection in cms
      if (dphi     > m_taub2b3AngleCut[0] &&
          dphi     < m_taub2b3AngleCut[1] &&
          thetaSum > m_taub2b3AngleCut[2] &&
          thetaSum < m_taub2b3AngleCut[3]) {
        // Cluster ThetaID selection
        if (MaxTCThetaId[icluster] >=  2 &&
            MaxTCThetaId[icluster] <= 16 &&
            MaxTCThetaId[jcluster] >=  2 &&
            MaxTCThetaId[jcluster] <= 16) {
          // Cluster energy selection in lab
          if (ClusterEnergy[icluster] > m_taub2b3CLEb2bCut ||
              ClusterEnergy[jcluster] > m_taub2b3CLEb2bCut) {
            taub2b3AngleCLEThetaIdFlag++;
          }
        }
      }
    }
  }
  // all selections
  bool taub2b3Flag = false;
  if (taub2b3EtotFlag            > 0 &&
      taub2b3AngleCLEThetaIdFlag > 0 &&
      taub2b3CLELowCutFlag       > 0 &&
      taub2b3CLEHighCutFlag      > 0) {
    taub2b3Flag = true;
  }
 
  return taub2b3Flag;
}

◆ getTaub2bAngleFlag()

int getTaub2bAngleFlag ( void )

inline

get taub2b 2 cluster angle cut flag

Definition at line 236 of file TrgEclBhabha.h.

236{ return m_taub2bAngleFlag; }

◆ getTaub2bClusterEFlag()

int getTaub2bClusterEFlag ( void )

inline

taub2b Cluster energy flag

Definition at line 240 of file TrgEclBhabha.h.

240{ return m_taub2bClusterEFlag; }

◆ getTaub2bEtotFlag()

int getTaub2bEtotFlag ( void )

inline

get total energy(TC theta id=1-17) flag for taub2b

Definition at line 238 of file TrgEclBhabha.h.

238{ return m_taub2bEtotFlag; }

◆ set2DBhabhaThreshold()

void set2DBhabhaThreshold	(	const std::vector< double > &	i2DBhabhaThresholdFWD,
		const std::vector< double > &	i2DBhabhaThresholdBWD
	)

inline

set 2D Bhabha Energy Threshold

Definition at line 69 of file TrgEclBhabha.h.

    {
      _2DBhabhaThresholdFWD = i2DBhabhaThresholdFWD;
      _2DBhabhaThresholdBWD = i2DBhabhaThresholdBWD;
    }

◆ set3DBhabhaAddAngleCut()

void set3DBhabhaAddAngleCut ( const std::vector< double > & i3DBhabhaAddAngleCut )

inline

set 3D Bhabha addtion Angle selection

Definition at line 100 of file TrgEclBhabha.h.

    {
      m_3DBhabhaAddAngleCut = i3DBhabhaAddAngleCut;
    }

◆ set3DBhabhaSelectionAngle()

void set3DBhabhaSelectionAngle ( const std::vector< double > & i3DBhabhaSelectionAngle )

inline

set 3D selection Bhabha Energy Angle

Definition at line 86 of file TrgEclBhabha.h.

    {
      _3DBhabhaSelectionAngle = i3DBhabhaSelectionAngle;
    };

◆ set3DBhabhaSelectionThreshold()

void set3DBhabhaSelectionThreshold ( const std::vector< double > & i3DBhabhaSelectionThreshold )

inline

set 3D selection Bhabha Energy Threshold

Definition at line 76 of file TrgEclBhabha.h.

    {
      _3DBhabhaSelectionThreshold = i3DBhabhaSelectionThreshold;
    };

◆ set3DBhabhaVetoAngle()

void set3DBhabhaVetoAngle ( const std::vector< double > & i3DBhabhaVetoAngle )

inline

set 3D veto Bhabha Energy Angle

Definition at line 91 of file TrgEclBhabha.h.

    {
      _3DBhabhaVetoAngle = i3DBhabhaVetoAngle;
    };

◆ set3DBhabhaVetoInTrackThetaRegion()

void set3DBhabhaVetoInTrackThetaRegion ( const std::vector< int > & i3DBhabhaVetoInTrackThetaRegion )

inline

set ThetaID (low and high) for 3DBhabhaVetoInTrack

Definition at line 105 of file TrgEclBhabha.h.

    {
      m_3DBhabhaVetoInTrackThetaRegion = i3DBhabhaVetoInTrackThetaRegion;
    }

◆ set3DBhabhaVetoThreshold()

void set3DBhabhaVetoThreshold ( const std::vector< double > & i3DBhabhaVetoThreshold )

inline

set 3D veto Bhabha Energy Threshold

Definition at line 81 of file TrgEclBhabha.h.

    {
      _3DBhabhaVetoThreshold = i3DBhabhaVetoThreshold;
    };

◆ setmumuAngle()

void setmumuAngle ( const std::vector< double > & imumuAngle )

inline

set mumu bit Angle selection

Definition at line 98 of file TrgEclBhabha.h.

98{_mumuAngle = imumuAngle; }

◆ setmumuThreshold()

void setmumuThreshold ( int mumuThreshold )

inline

set mumu bit Threshold

Definition at line 96 of file TrgEclBhabha.h.

96{_mumuThreshold = mumuThreshold; }

◆ setTaub2b2Cut()

void setTaub2b2Cut	(	const std::vector< int > &	iTaub2b2AngleCut,
		const double	iTaub2b2EtotCut,
		const double	iTaub2b2CLEEndcapCut,
		const double	iTaub2b2CLECut
	)

inline

set taub2b2 cut

Definition at line 128 of file TrgEclBhabha.h.

    {
      m_taub2b2AngleCut     = iTaub2b2AngleCut;
      m_taub2b2EtotCut      = iTaub2b2EtotCut;
      m_taub2b2CLEEndcapCut = iTaub2b2CLEEndcapCut;
      m_taub2b2CLECut       = iTaub2b2CLECut;
    }

◆ setTaub2b3Cut()

void setTaub2b3Cut	(	const std::vector< int > &	iTaub2b3AngleCut,
		const double	iTaub2b3EtotCut,
		const double	iTaub2b3CLEb2bCut,
		const double	iTaub2b3CLELowCut,
		const double	iTaub2b3CLEHighCut
	)

inline

set taub2b3 cut

Definition at line 139 of file TrgEclBhabha.h.

    {
      m_taub2b3AngleCut     = iTaub2b3AngleCut;
      m_taub2b3EtotCut      = iTaub2b3EtotCut;
      m_taub2b3CLEb2bCut    = iTaub2b3CLEb2bCut;
      m_taub2b3CLELowCut    = iTaub2b3CLELowCut;
      m_taub2b3CLEHighCut   = iTaub2b3CLEHighCut;
    }

◆ setTaub2bAngleCut()

void setTaub2bAngleCut ( const std::vector< int > & iTaub2bAngleCut )

inline

set 2 Cluster angle selection for tau 1x1 decay [0], [1] for low and high of dphi, [2], [3] for low and high of Theta Sum

Definition at line 111 of file TrgEclBhabha.h.

    {
      m_taub2bAngleCut = iTaub2bAngleCut;
    }

◆ setTaub2bClusterECut()

void setTaub2bClusterECut	(	const double	iTaub2bClusterECut1,
		const double	iTaub2bClusterECut2
	)

inline

set cluster energy cut for taub2b

Definition at line 121 of file TrgEclBhabha.h.

    {
      m_taub2bClusterECut1 = iTaub2bClusterECut1;
      m_taub2bClusterECut2 = iTaub2bClusterECut2;
    }

◆ setTaub2bEtotCut()

void setTaub2bEtotCut ( const double iTaub2bEtotCut )

inline

set total energy cut for taub2b

Definition at line 116 of file TrgEclBhabha.h.

    {
      m_taub2bEtotCut = iTaub2bEtotCut;
    }

Member Data Documentation

◆ _2DBhabhaThresholdBWD

std::vector<double> _2DBhabhaThresholdBWD

private

2D Bhabha Energy Threshold

Definition at line 267 of file TrgEclBhabha.h.

◆ _2DBhabhaThresholdFWD

std::vector<double> _2DBhabhaThresholdFWD

private

2D Bhabha Energy Threshold

Definition at line 265 of file TrgEclBhabha.h.

◆ _3DBhabhaSelectionAngle

std::vector<double> _3DBhabhaSelectionAngle

private

3D Selection Bhabha Energy Angle

Definition at line 273 of file TrgEclBhabha.h.

◆ _3DBhabhaSelectionThreshold

std::vector<double> _3DBhabhaSelectionThreshold

private

3D Selection Bhabha Energy Threshold

Definition at line 269 of file TrgEclBhabha.h.

◆ _3DBhabhaVetoAngle

std::vector<double> _3DBhabhaVetoAngle

private

3D Veto Bhabha Energy Angle

Definition at line 275 of file TrgEclBhabha.h.

◆ _3DBhabhaVetoThreshold

std::vector<double> _3DBhabhaVetoThreshold

private

3D Veto Bhabha Energy Threshold

Definition at line 271 of file TrgEclBhabha.h.

◆ _database

TrgEclDataBase* _database = nullptr

private

Object of Trigger ECL DataBase.

cppcheck-suppress unsafeClassCanLeak

Definition at line 247 of file TrgEclBhabha.h.

◆ _mumuAngle

std::vector<double> _mumuAngle

private

mumu bit Angle

Definition at line 279 of file TrgEclBhabha.h.

◆ _mumuThreshold

double _mumuThreshold

private

mumu bit Energy Threshold

Definition at line 277 of file TrgEclBhabha.h.

◆ _TCMap

TrgEclMapping* _TCMap = nullptr

private

Object of TC Mapping.

Definition at line 244 of file TrgEclBhabha.h.

◆ BhabhaComb

std::vector<double> BhabhaComb

private

Bhabha Combination.

Definition at line 253 of file TrgEclBhabha.h.

◆ ClusterEnergy

std::vector<double> ClusterEnergy

private

Cluster Energy.

Definition at line 259 of file TrgEclBhabha.h.

◆ ClusterPosition

std::vector<ROOT::Math::XYZVector> ClusterPosition

private

Cluster Timing.

Definition at line 263 of file TrgEclBhabha.h.

◆ ClusterTiming

std::vector<double> ClusterTiming

private

Cluster Timing.

Definition at line 261 of file TrgEclBhabha.h.

◆ m_3DBhabhaAddAngleCut

std::vector<double> m_3DBhabhaAddAngleCut

private

Angle selection of additional Bhabha veto in CM frame.

Definition at line 281 of file TrgEclBhabha.h.

◆ m_3DBhabhaSelectionClusterEnergies

std::vector<double> m_3DBhabhaSelectionClusterEnergies

private

Energies of two clusters used for 3D Bhabha selection (GeV)

Definition at line 337 of file TrgEclBhabha.h.

◆ m_3DBhabhaSelectionClusterTCIds

std::vector<int> m_3DBhabhaSelectionClusterTCIds

private

TCIDs of two clusters of 3D Bhabha selection.

Definition at line 333 of file TrgEclBhabha.h.

◆ m_3DBhabhaSelectionClusterThetaIds

std::vector<int> m_3DBhabhaSelectionClusterThetaIds

private

ThetaIDs of two clusters used for 3D Bhabha selection.

Definition at line 335 of file TrgEclBhabha.h.

◆ m_3DBhabhaSelectionClusterTimings

std::vector<double> m_3DBhabhaSelectionClusterTimings

private

Timings of two clusters used for 3D Bhabha selection (ns)

Definition at line 339 of file TrgEclBhabha.h.

◆ m_3DBhabhaSelectionThetaFlag

int m_3DBhabhaSelectionThetaFlag = std::numeric_limits<int>::quiet_NaN()

private

flag which shows theta position of clusters of 3DBhabha Selection.

Definition at line 331 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoClusterEnergies

std::vector<double> m_3DBhabhaVetoClusterEnergies

private

Energies of two clusters of 3D Bhabha veto (GeV)

Definition at line 325 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoClusterTCIds

std::vector<int> m_3DBhabhaVetoClusterTCIds

private

TCIDs of two clusters of 3D Bhabha veto.

Definition at line 321 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoClusterThetaIds

std::vector<int> m_3DBhabhaVetoClusterThetaIds

private

ThetaIds of two clusters of 3D Bhabha veto.

Definition at line 323 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoClusterTimings

std::vector<double> m_3DBhabhaVetoClusterTimings

private

Timings of two clusters of 3D Bhabha veto (ns)

Definition at line 327 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoInTrackFlag

int m_3DBhabhaVetoInTrackFlag

private

trigger bit of flag(1bit) whether two clusters statisfy 3D Bhabha veto are in CDCTRG region in theta (="InTrack") or not

Definition at line 319 of file TrgEclBhabha.h.

◆ m_3DBhabhaVetoInTrackThetaRegion

std::vector<int> m_3DBhabhaVetoInTrackThetaRegion

private

theta region(low and high) of 3D Bhbabha veto InTrack

Definition at line 329 of file TrgEclBhabha.h.

◆ m_taub2b2AngleCut

std::vector<int> m_taub2b2AngleCut

private

taub2b2 two Cluster angle cut (degree)

Definition at line 299 of file TrgEclBhabha.h.

◆ m_taub2b2CLECut

double m_taub2b2CLECut

private

taub2b2 cluster energy cut (GeV)

Definition at line 303 of file TrgEclBhabha.h.

◆ m_taub2b2CLEEndcapCut

double m_taub2b2CLEEndcapCut

private

taub2b2 cluster energy cut for endcap cluster (GeV)

Definition at line 301 of file TrgEclBhabha.h.

◆ m_taub2b2EtotCut

double m_taub2b2EtotCut

private

taub2b2 total energy cut (GeV)

Definition at line 297 of file TrgEclBhabha.h.

◆ m_taub2b3AngleCut

std::vector<int> m_taub2b3AngleCut

private

taub2b3 selection cuts (3,2,1,0) = (dphi low, dphi high, theta_sum low, theta_sum high)

Definition at line 307 of file TrgEclBhabha.h.

◆ m_taub2b3CLEb2bCut

double m_taub2b3CLEb2bCut

private

taub2b3 cluster energy cut (GeV)

Definition at line 311 of file TrgEclBhabha.h.

◆ m_taub2b3CLEHighCut

double m_taub2b3CLEHighCut

private

taub2b3 cluster energy cut3 (GeV)

Definition at line 315 of file TrgEclBhabha.h.

◆ m_taub2b3CLELowCut

double m_taub2b3CLELowCut

private

taub2b3 cluster energy cut2 (GeV)

Definition at line 313 of file TrgEclBhabha.h.

◆ m_taub2b3EtotCut

double m_taub2b3EtotCut

private

taub2b3 total energy (TC theta ID =1-17) (GeV)

Definition at line 309 of file TrgEclBhabha.h.

◆ m_taub2bAngleCut

std::vector<int> m_taub2bAngleCut

private

taub2b 2 Cluster angle cut (degree)

Definition at line 283 of file TrgEclBhabha.h.

◆ m_taub2bAngleFlag

int m_taub2bAngleFlag

private

taub2b 2 cluster angle cut flag

Definition at line 291 of file TrgEclBhabha.h.

◆ m_taub2bClusterECut1

double m_taub2bClusterECut1

private

taub2b Cluster one of energy cut in b2b in lab (GeV)

Definition at line 287 of file TrgEclBhabha.h.

◆ m_taub2bClusterECut2

double m_taub2bClusterECut2

private

taub2b Cluster one of energy cut in b2b in lab (GeV)

Definition at line 289 of file TrgEclBhabha.h.

◆ m_taub2bClusterEFlag

int m_taub2bClusterEFlag

private

taub2b Cluster energy flag

Definition at line 295 of file TrgEclBhabha.h.

◆ m_taub2bEtotCut

double m_taub2bEtotCut

private

taub2b total energy(TC theta id=1-17) cut (GeV)

Definition at line 285 of file TrgEclBhabha.h.

◆ m_taub2bEtotFlag

int m_taub2bEtotFlag

private

taub2b total energy(TC theta id=1-17) flag

Definition at line 293 of file TrgEclBhabha.h.

◆ MaxTCId

std::vector<double> MaxTCId

private

Max TC Id.

Definition at line 255 of file TrgEclBhabha.h.

◆ MaxTCThetaId

std::vector<double> MaxTCThetaId

private

Cluster ThetaId.

Definition at line 257 of file TrgEclBhabha.h.

The documentation for this class was generated from the following files:

trg/ecl/include/TrgEclBhabha.h
trg/ecl/src/TrgEclBhabha.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrgEclBhabha() [1/2]

◆ ~TrgEclBhabha()

◆ TrgEclBhabha() [2/2]

Member Function Documentation

◆ get2CLETP()

◆ get3DBhabhaSelectionClusterEnergy()

◆ get3DBhabhaSelectionClusterTCId()

◆ get3DBhabhaSelectionClusterThetaId()

◆ get3DBhabhaSelectionClusterTiming()

◆ get3DBhabhaSelectionThetaFlag()

◆ get3DBhabhaVetoClusterEnergy()

◆ get3DBhabhaVetoClusterTCId()

◆ get3DBhabhaVetoClusterThetaId()

◆ get3DBhabhaVetoClusterTiming()

◆ get3DBhabhaVetoInTrackFlag()

◆ GetBhabha00()

◆ GetBhabha01()

◆ GetBhabha02()

◆ GetBhabhaAddition()

◆ GetBhabhaComb()

◆ Getmumu()

◆ GetTaub2b()

◆ GetTaub2b2()

◆ GetTaub2b3()

◆ getTaub2bAngleFlag()

◆ getTaub2bClusterEFlag()

◆ getTaub2bEtotFlag()

◆ set2DBhabhaThreshold()

◆ set3DBhabhaAddAngleCut()

◆ set3DBhabhaSelectionAngle()

◆ set3DBhabhaSelectionThreshold()

◆ set3DBhabhaVetoAngle()

◆ set3DBhabhaVetoInTrackThetaRegion()

◆ set3DBhabhaVetoThreshold()

◆ setmumuAngle()

◆ setmumuThreshold()

◆ setTaub2b2Cut()

◆ setTaub2b3Cut()

◆ setTaub2bAngleCut()

◆ setTaub2bClusterECut()

◆ setTaub2bEtotCut()

Member Data Documentation

◆ _2DBhabhaThresholdBWD

◆ _2DBhabhaThresholdFWD

◆ _3DBhabhaSelectionAngle

◆ _3DBhabhaSelectionThreshold

◆ _3DBhabhaVetoAngle

◆ _3DBhabhaVetoThreshold

◆ _database

◆ _mumuAngle

◆ _mumuThreshold

◆ _TCMap

◆ BhabhaComb

◆ ClusterEnergy

◆ ClusterPosition

◆ ClusterTiming

◆ m_3DBhabhaAddAngleCut

◆ m_3DBhabhaSelectionClusterEnergies

◆ m_3DBhabhaSelectionClusterTCIds

◆ m_3DBhabhaSelectionClusterThetaIds

◆ m_3DBhabhaSelectionClusterTimings

◆ m_3DBhabhaSelectionThetaFlag

◆ m_3DBhabhaVetoClusterEnergies

◆ m_3DBhabhaVetoClusterTCIds

◆ m_3DBhabhaVetoClusterThetaIds

◆ m_3DBhabhaVetoClusterTimings

◆ m_3DBhabhaVetoInTrackFlag

◆ m_3DBhabhaVetoInTrackThetaRegion

◆ m_taub2b2AngleCut

◆ m_taub2b2CLECut

◆ m_taub2b2CLEEndcapCut

◆ m_taub2b2EtotCut

◆ m_taub2b3AngleCut

◆ m_taub2b3CLEb2bCut

◆ m_taub2b3CLEHighCut