TrgEclMaster Class Reference

ETM class. More...

#include <TrgEclMaster.h>

Public Member Functions
	TrgEclMaster (void)
	TrgEclMaster Constructor.
virtual	~TrgEclMaster ()
	TrgEclMaster Destructor.
	TrgEclMaster (TrgEclMaster &)=delete
	Copy constructor, deleted.
TrgEclMaster &	operator= (TrgEclMaster &)=delete
	Assignment operator, deleted.
void	initialize (int)
	initialize
void	simulate01 (int)
	simulates ECL trigger for Global Cosmic data
void	simulate02 (int)
	simulates ECL trigger for Data Analysis
std::string	name (void) const
	returns name.
std::string	version (void) const
	returns version.
void	setRS (std::vector< int >, std::vector< double >, std::vector< double > &, std::vector< std::vector< double > > &)
	ECL bit information for GDL.
void	setClusterMethod (int cluster)
	Get Event timing.
void	setClusterLimit (int limit)
	Set the limit # of Cluster.
void	setBhabhaMethod (int bhabha)
	Set Bhabha.
void	setEventTimingMethod (int EventTiming)
	Set Cluster.
void	setTimeWindow (int timewindow)
	Set Trigger Decision window size.
void	setOverlapWindow (int overlapwindow)
	Set Trigger Decision overlap window size.
void	setNofTopTC (int noftoptc)
	set # of considered TC in energy weighted Timing method
void	makeLowMultiTriggerBit (std::vector< int >, std::vector< double >)
	make LowMultiTriggerBit
void	makeTriggerBit (int, int, int, int, double, int, int, std::vector< int >, int, int, int, int, int, int, int, int, int, int, int, int, int, int, int, int)
	make Trigger bit except for Low Multiplicity related bit
double	setTotalEnergy (std::vector< double >)
	Set Total Energy.
int	getTriggerbit (int i)
	Get ECL Trigger bit.
int	getLowmultibit ()
	Get Low Multiplicity Trigger Bit.
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	setTaub2bAngleCut (const std::vector< int > &itaub2bAngleCut)
	set tau b2b 2 cluster angle cut
void	setTaub2bEtotCut (double itaub2bEtotCut)
	set tau b2b total energy cut
void	setTaub2bClusterECut (double itaub2bClusterECut1, double itaub2bClusterECut2)
	set tau b2b 1Cluster energy cut
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
void	setn300MeVClusterThreshold (int n300MeVCluster)
	set the number of cluster exceeding 300 MeV
void	setECLBurstThreshold (int ECLBurstThreshold)
	set mumu bit Threshold
void	setTotalEnergyThreshold (const std::vector< double > &iTotalEnergy)
	set Total Energy Theshold (low, high, lum)
void	setLowMultiplicityThreshold (const std::vector< double > &iLowMultiThreshold)
	set Low Multiplicity Threshold
void	set3DBhabhaVetoInTrackThetaRegion (const std::vector< int > &i3DBhabhaVetoInTrackThetaRegion)
	set theta ID region (low and high) of 3DBhabhaVeto InTrack for gg selection
void	setEventTimingQualityThresholds (const std::vector< double > &iEventTimingQualityThresholds)
	set energy threshold(low and high) of event timing quality flag (GeV)

Static Public Member Functions
static TrgEclMaster *	getTrgEclMaster (void)
	get pointer of TrgEclMaster object

Private Attributes
std::vector< std::vector< double > >	m_TCEnergy
	Hit TC Energy.
std::vector< std::vector< double > >	m_TCTiming
	Hit TC Timing.
std::vector< std::vector< int > >	m_TCBeamBkgTag
	Hit TC Beam Background tag.
std::vector< int >	m_HitTCId
	Hit TC Energy in time window.
std::vector< double >	m_TCHitEnergy
	Hit TC Energy in time window.
std::vector< double >	m_TCHitTiming
	Hit TC Timing in time window.
std::vector< int >	m_TCHitBeamBkgTag
	Hit TC Beam Background tag in time window.
std::vector< std::vector< double > >	m_PhiRingSum
	Phi ring sum.
std::vector< std::vector< std::vector< double > > >	m_ThetaRingSum
	Theta ring sum.
double	m_TimeWindow
	Hit TC Energy in time window.
double	m_OverlapWindow
	TRG Decision overlap window.
int	m_Clustering
	clutering option
int	m_Bhabha
	Bhabha option.
int	m_EventTiming
	EventTiming option.
int	m_NofTopTC
int	m_ClusterLimit
	The limit number of Cluster.
int	m_Triggerbit [4]
	ECL Trigger bit.
int	m_Lowmultibit
	Low Multiplicity bit.
int	m_PrescaleFactor
	Bhabha Prescale Factor.
int	m_PrescaleCounter
	Bhabha Prescale Countor.
std::vector< double >	m_2DBhabhaThresholdFWD
	2D Bhabha Energy Threshold
std::vector< double >	m_2DBhabhaThresholdBWD
	2D Bhabha Energy Threshold
std::vector< double >	m_3DBhabhaSelectionThreshold
	3D Selection Bhabha Energy Threshold
std::vector< double >	m_3DBhabhaVetoThreshold
	3D Veto Bhabha Energy Threshold
std::vector< double >	m_3DBhabhaSelectionAngle
	3D Selection Bhabha Energy Angle
std::vector< double >	m_3DBhabhaVetoAngle
	3D Veto Bhabha Energy Angle
double	m_mumuThreshold
	mumu bit Energy Threshold
std::vector< double >	m_mumuAngle
	mumu bit Angle
std::vector< double >	m_3DBhabhaAddAngleCut
	Angle selection of additional Bhabha addition in CM frame.
std::vector< int >	m_taub2bAngleCut
	tau b2b 2 cluster angle cut (degree) (dphi low, dphi high, theta_sum low, theta_sum high)
double	m_taub2bEtotCut
	tau b2b total energy (TC theta ID =1-17) (GeV)
double	m_taub2bClusterECut1
	taub2b one Cluster energy selection (GeV)
double	m_taub2bClusterECut2
	taub2b one Cluster energy selection (GeV)
std::vector< int >	m_taub2b2AngleCut
	taub2b2 angle selection(degree) (3,2,1,0) = (dphi low, dphi high, theta_sum low, theta_sum high)
double	m_taub2b2EtotCut
	taub2b2 total energy (TC theta ID =1-17) (GeV)
double	m_taub2b2CLEEndcapCut
	taub2b2 cluster energy cut for endcap cluster (GeV)
double	m_taub2b2CLECut
	taub2b2 cluseter 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 in lab for one of b2b clusters (GeV)
double	m_taub2b3CLELowCut
	taub2b3 cluster energy low cut in lab for all clusters (GeV)
double	m_taub2b3CLEHighCut
	taub2b3 cluster energy high cut in lab for all clusters (GeV)
int	m_n300MeVCluster
	The number of Cluster exceeding 300 MeV.
double	m_ECLBurstThreshold
	ECL Burst Bit Threshold.
std::vector< double >	m_TotalEnergy
	Total Energy Theshold (low, high, lum)
std::vector< double >	m_LowMultiThreshold
	Low Multiplicity Threshold.
std::vector< int >	m_3DBhabhaVetoInTrackThetaRegion
	Theta region (low, high) of 3D Bhabha Veto InTrack.
std::vector< double >	m_EventTimingQualityThresholds
	energy threshold(low, high) for quality flag (GeV)
TrgEclMapping *	m_obj_map
	Mapping object.
TrgEclCluster *	m_obj_cluster
	Cluster object.
TrgEclTiming *	m_obj_timing
	EventTiming object.
TrgEclBhabha *	m_obj_bhabha
	Bhabha object.
TrgEclBeamBKG *	m_obj_beambkg
	Beam Backgroud veto object.
TrgEclDataBase *	m_obj_database
	Beam Backgroud veto object.

Detailed Description

ETM class.

Definition at line 33 of file TrgEclMaster.h.

Constructor & Destructor Documentation

TrgEclMaster()

TrgEclMaster

(

void

)

TrgEclMaster Constructor.

Definition at line 83 of file TrgEclMaster.cc.

                          :
  m_TimeWindow(250.0), m_OverlapWindow(0.0), m_Clustering(1), m_Bhabha(0),
  m_EventTiming(1), m_NofTopTC(3), m_ClusterLimit(6), m_Triggerbit{0, 0, 0, 0},
  m_Lowmultibit(0), m_PrescaleFactor(0), m_PrescaleCounter(0), m_mumuThreshold(20),
  m_taub2bEtotCut(7.0),
  m_taub2bClusterECut1(1.9),
  m_taub2bClusterECut2(999),
  m_taub2b2EtotCut(7.0),
  m_taub2b2CLEEndcapCut(3.0),
  m_taub2b2CLECut(0.162),
  m_taub2b3EtotCut(7.0),
  m_taub2b3CLEb2bCut(0.14),
  m_taub2b3CLELowCut(0.12),
  m_taub2b3CLEHighCut(4.5),
  m_n300MeVCluster(1),
  m_ECLBurstThreshold(200)
{
 
  m_obj_cluster  = new TrgEclCluster();
  m_obj_beambkg  = new TrgEclBeamBKG();
  m_obj_bhabha   = new TrgEclBhabha();
  m_obj_timing   = new TrgEclTiming();
  m_obj_map      = new TrgEclMapping();
  m_obj_database = new TrgEclDataBase();
}

~TrgEclMaster()

~TrgEclMaster ( )

virtual

TrgEclMaster Destructor.

Definition at line 111 of file TrgEclMaster.cc.

{
  delete m_obj_cluster;
  delete m_obj_beambkg;
  delete m_obj_bhabha;
  delete m_obj_timing;
  delete m_obj_map;
  delete m_obj_database;
}

Member Function Documentation

getLowmultibit()

int getLowmultibit ( )

inline

Get Low Multiplicity Trigger Bit.

Definition at line 111 of file TrgEclMaster.h.

{return m_Lowmultibit;}

getTriggerbit()

int getTriggerbit ( int  i )

inline

Get ECL Trigger bit.

Definition at line 109 of file TrgEclMaster.h.

{return m_Triggerbit[i];}

initialize()

void initialize

(

int

)

initialize

Definition at line 140 of file TrgEclMaster.cc.

{
  m_TCEnergy.clear();
  m_TCTiming.clear();
  m_TCBeamBkgTag.clear();
  m_HitTCId.clear();
  m_TCHitEnergy.clear();
  m_TCHitTiming.clear();
  m_TCHitBeamBkgTag.clear();
 
  m_TCEnergy.resize(576);
  m_TCTiming.resize(576);
  m_TCBeamBkgTag.resize(576);
 
  m_TotalEnergy = {5, 10, 30}; // /100 MeV
  m_2DBhabhaThresholdFWD  =
  {40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 30, 35}; // /100 MeV
  m_2DBhabhaThresholdBWD  =
  {25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 30, 30}; // /100 MeV
  m_3DBhabhaVetoThreshold = {30, 45}; //  /100 MeV
  m_3DBhabhaSelectionThreshold = {20, 40}; //  /100 MeV
  m_3DBhabhaVetoAngle      = {160, 200, 165, 190}; //  /100 MeV
  m_3DBhabhaSelectionAngle = {140, 220, 160, 200}; //  /100 MeV
  m_mumuAngle              = {160, 200, 165, 190}; //  degree
  m_3DBhabhaAddAngleCut   = {150, 210, 160, 200};  //  degree
  m_LowMultiThreshold      = {10, 20, 25, 30}; // degree
 
  m_3DBhabhaVetoInTrackThetaRegion = {3, 15};
  m_EventTimingQualityThresholds = {5, 6}; // GeV
  m_ECLBurstThreshold = 200; // 10 ADC
 
  // taub2b cut
  m_taub2bAngleCut = {110, 250, 130, 230}; // degree
  m_taub2bEtotCut = 7.0; // GeV
  m_taub2bClusterECut1 = 1.9; // GeV
  m_taub2bClusterECut2 = 999; // GeV
 
  // taub2b2 cut
  m_taub2b2AngleCut = {120, 240, 140, 220}; // degree
  m_taub2b2EtotCut = 7.0; // GeV
  m_taub2b2CLEEndcapCut = 3.0; // GeV
  m_taub2b2CLECut = 0.162; // GeV
 
  //taub2b3 cut
  m_taub2b3AngleCut   = {120, 240, 140, 220}; // degree
  m_taub2b3EtotCut    = 7.0; // GeV
  m_taub2b3CLEb2bCut  = 0.14; // GeV
  m_taub2b3CLELowCut  = 0.12; // GeV
  m_taub2b3CLEHighCut = 4.5;  // GeV
 
}

makeLowMultiTriggerBit()

void makeLowMultiTriggerBit	(	std::vector< int >	CenterTCId,
		std::vector< double >	clusterenergy
	)

make LowMultiTriggerBit

Definition at line 1275 of file TrgEclMaster.cc.

{
 
  //---------------------------------------------------------------------------------
  // ECL trigger
  //---------------------------------------------------------------------------------
  //Variable(in Tsim) | N(bit) |   Address     | Parameter
  //-------------------------------------- ------------------------------------------
  //                  |      1 |         62    |  N Cluster  >= 3, at least one Cluster >300 MeV (LAB), not 3D ECL Bhabha
  //                  |      1 |         63    | one Cluster >= 2GeV(CM) with Theta Id = 4~14
  //  m_Lowmultibit    |      1 |         64    | one Cluster >= 2GeV(CM) with Theta Id = 2,3,15 or 16 and not a 3D ECL Bhabha
  //                  |      1 |         65    | one Cluster >= 2GeV(CM) with Theta Id = 2, 3, 15 or 16 and not a 3D ECL Bhabha
  //                  |      1 |         66    | one Cluster >= 2GeV(CM) with Theta Id = 1 or 17 and not a 3D ECL Bhabha
  //                  |      1 |         67    | one Cluster >= 2GeV(CM) with Theta Id = 1 or 17 and a 3D ECL Bhabha
  //                  |      1 |         68    | exactly one Cluster >= 1GeV(CM) and one Cluster > 300  MeV (LAB ), in Theta Id 4~15(Barrel)
  //                  |      1 |         69    | exactly one Cluster >= 1GeV(CM) and one Cluster > 300 MeV (LAB), in Theta Id 2, 3 or 16
  //                  |      1 |         70    | 170 < delta phi(CM) < 190 degree, both Clusters > 250 MeV (LAB), and no 2GeV (CM) Cluster
  //                  |      1 |         71    | 170 < delta phi(CM) < 190 degree, one Cluster < 250 MeV (LAB), the other Cluster > 250 MeV(LAB), and no 2GeV (CM) Cluster
  //                  |      1 |         72    | 160 < delta phi(CM) < 200 degree, 160 < Sum Theta (CM)< 200 degree, no 2 GeV(CM) cluster
  //                  |      1 |         73    | No 2GeV (CM) Cluster
  //---------------------------------------------------------------------------------
  m_Lowmultibit = 0;
  int _nClust = CenterTCId.size();
  int _n300MeV = 0;
  int _n2GeV = 0;
  int _n2GeV414 = 0;
  int _n2GeV231516 = 0;
  int _n2GeV117 = 0;
  int _n1GeV415 = 0;
  int _n1GeV2316 = 0;
  int _nClust216 = 0;
  int _n500MeV216 = 0;
  int _n500MeV611 = 0;
  for (int ic = 0; ic < _nClust; ic++) {
    if (clusterenergy[ic] > 0.3) {_n300MeV++;}
    int thetaid = m_obj_map->getTCThetaIdFromTCId(CenterTCId[ic]);
    int lut = m_obj_database->Get3DBhabhaLUT(CenterTCId[ic]);
    int thresh = 15 & lut;
    if (thetaid >= 2 && thetaid <= 16) {_nClust216++;}
    if (thetaid >= 6 && thetaid <= 11) {
      if (clusterenergy[ic] * 100 > 5 * thresh) {
        _n500MeV611++;
      }
    }
 
    if (clusterenergy[ic] > 0.5 && thetaid >= 2 && thetaid <= 16) {_n500MeV216++;}
    if (clusterenergy[ic] * 100 > (thresh * m_LowMultiThreshold[1])) { //200 <MeV
      _n2GeV++;
      if (thetaid >= 4 && thetaid <= 14) {_n2GeV414++;}
      if (thetaid == 2 || thetaid == 3 || thetaid == 15 || thetaid == 16) {_n2GeV231516++;}
      if (thetaid == 1 || thetaid == 17) {_n2GeV117++;}
    }
    if (clusterenergy[ic] * 100 > thresh * m_LowMultiThreshold[0]) { // 100 MeV
      if (thetaid >= 4 && thetaid <= 15) {_n1GeV415++;}
      if (thetaid == 2 || thetaid == 3 || thetaid == 16) {_n1GeV2316++;}
    }
  }
  //---------------------------------------------------------------------
  //..Trigger objects using back-to-back ECL clusters, plus Bhabha vetoes
  //  nPhiPairHigh nPhiPairLow n3DPair nECLBhabha nTrkBhabha
 
  int _nPhiPairHigh = 0;
  int _nPhiPairLow = 0;
  int _n3DPair = 0;
  int _nECLBhabha = 0;
  for (int i0 = 0; i0 < _nClust - 1; i0++) {
    for (int i1 = i0 + 1; i1 < _nClust; i1++) {
      int lut1 = m_obj_database->Get3DBhabhaLUT(CenterTCId[i0]);
      int lut2 = m_obj_database->Get3DBhabhaLUT(CenterTCId[i1]);
 
      int energy1 = 15 & lut1;
      int 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;
 
      //..back to back in phi
      int dphi = abs(phi1 - phi2);
      if (dphi > 180) {dphi = 360 - dphi;}
      int thetaSum = theta1 + theta2;
 
      // cout << dphi << " "  << thetaSum << endl;
      // cout << clusterenergy[i0] << " " << clusterenergy[i1] << endl;
      //  if (dphi > 180.) {dphi = 360 - dphi;}
      if (dphi > 170. && clusterenergy[i0] > m_LowMultiThreshold[2] / 100
          && clusterenergy[i1] >  m_LowMultiThreshold[2] / 100) {_nPhiPairHigh++;}
      if (dphi > 170. &&
          ((clusterenergy[i0] <  m_LowMultiThreshold[2] / 100 &&
            clusterenergy[i1] >  m_LowMultiThreshold[2] / 100) ||
           (clusterenergy[i0] >  m_LowMultiThreshold[2] / 100 &&
            clusterenergy[i1] <  m_LowMultiThreshold[2] / 100))) {_nPhiPairLow++;}
      //..3D
      if (dphi > 160. && thetaSum > 160. && thetaSum < 200) {_n3DPair++;}
      //..ecl Bhabha
      if (dphi > 160 &&
          thetaSum > 165 &&
          thetaSum < 190 &&
          clusterenergy[i0] * 100 > m_3DBhabhaVetoThreshold[0] * energy1 &&
          clusterenergy[i1] * 100 > m_3DBhabhaVetoThreshold[0]  * energy2 &&
          (clusterenergy[i0] * 100 > m_3DBhabhaVetoThreshold[1]  * energy1 ||
           clusterenergy[i1] * 100 > m_3DBhabhaVetoThreshold[1]  * energy2)) {
        _nECLBhabha++;
 
      }
    }
  }
  int bit1 = 0;
  int bit2 = 0;
  int bit3 = 0;
  int bit4 = 0;
  int bit5 = 0;
  int bit6 = 0;
  int bit7 = 0;
  int bit8 = 0;
  int bit9 = 0;
  int bit10 = 0;
  int bit11 = 0;
  int bit12 = 0;
  int bit13 = 0;
  int bit14 = 0;
 
 
  if (_nClust >= 3 && _n300MeV >= m_n300MeVCluster && _nECLBhabha == 0) {
    bit1 = 0x01; //6
  }
  if (_n2GeV414 > 0) {
    bit2 = 0x01; //7
  }
  if (_n2GeV231516 && _nECLBhabha == 0) {
    bit3 = 0x01; //9
  }
  if (_n2GeV231516 && _nECLBhabha != 0) {
    bit4 = 0x01; //10
  }
  if (_n2GeV117 && _nECLBhabha == 0) {
    bit5 = 0x01; //11
  }
  if (_n2GeV117 && _nECLBhabha != 0) {
    bit6 = 0x01; //12
  }
  if (_n1GeV415 == 1 && _n300MeV == 1) {
    bit7 = 0x01; //13
  }
  if (_n1GeV2316 == 1 && _n300MeV == 1) {
    bit8 = 0x01; //14
  }
  if (_nPhiPairHigh > 0 && _n2GeV == 0) {
    bit9  = 0x01; //15
  }
  if (_nPhiPairLow > 0 && _n2GeV == 0) {
    bit10 = 0x01; //16
  }
  if (_n3DPair > 0 && _n2GeV == 0) {
    bit11 = 0x01; //17;
  }
  if (_n2GeV == 0) {
    bit12 = 0x01; //4
  }
  if (_nClust216 >= 3 && _n500MeV216 >  0 && _nECLBhabha == 0) {
    bit13 = 0x01; //6
  }
  if (_n500MeV611 == 1 && _n300MeV == 1) {
    bit14 = 0x01; //6
  }
 
  int  total_bit = 0;
  total_bit |= bit14;
  total_bit <<= 1;
  total_bit |= bit13;
  total_bit <<= 1;
  total_bit |= bit12;
  total_bit <<= 1;
  total_bit |= bit11;
  total_bit <<= 1;
  total_bit |= bit10;
  total_bit <<= 1;
  total_bit |= bit9;
  total_bit <<= 1;
  total_bit |= bit8;
  total_bit <<= 1;
  total_bit |= bit7;
  total_bit <<= 1;
  total_bit |= bit6;
  total_bit <<= 1;
  total_bit |= bit5;
  total_bit <<= 1;
  total_bit |= bit4;
  total_bit <<= 1;
  total_bit |= bit3;
  total_bit <<= 1;
  total_bit |= bit2;
  total_bit <<= 1;
  total_bit |= bit1;
 
  m_Lowmultibit =  total_bit ;
}

makeTriggerBit()

void makeTriggerBit	(	int	hit,
		int	Timing,
		int	RevoFAM,
		int	TimingSource,
		double	etot,
		int	bhabha2D,
		int	physics,
		std::vector< int >	bhabhatype,
		int	ICN,
		int	BGVeto,
		int	ClusterOverflow,
		int	bhabha3D,
		int	lowmultibit,
		int	bhabha3D_sel,
		int	mumubit,
		int	prescale,
		int	burst,
		int	EventTimingQualityFlag,
		int	bhabha3DVetoInTrackFlag,
		int	bhabha3DSelectionThetaFlag,
		int	taub2bFlag,
		int	bit_hie_bhav,
		int	taub2b2Flag,
		int	taub2b3Flag
	)

make Trigger bit except for Low Multiplicity related bit

Definition at line 1119 of file TrgEclMaster.cc.

{
 
  m_Triggerbit[0] = 0;
  m_Triggerbit[1] = 0;
  m_Triggerbit[2] = 0;
  m_Triggerbit[3] = 0;
 
  //  int physics = 0;
  int elow = 0;
  int ehigh = 0;
  int elum = 0;
  int bhabhaveto = 0;
  int Bhabhatype = bhabha2D;
 
  if (etot > 0.5) {
    elow = 0x01;
  }
  if (etot > 1.0) {
    ehigh = 0x01;
  }
  if (etot > 3.0) {
    elum = 0x01;
  }
 
  if (bhabhatype.size() > 14) {
    for (int ibhabha = 0; ibhabha < 13; ibhabha++) {
      int type = 0x00;
      if (bhabhatype[ibhabha] == 1) {type = 0x01;}
 
      Bhabhatype |= type;
      Bhabhatype <<= 1;
 
    }
  }
 
  int bit_hit = hit & 0x01;
  int bit_Timing = (Timing & 0x7F) ;
  int bit_RevoFAM = (RevoFAM & 0x7F) ;
  int bit_TimingSource  = (TimingSource & 0x07) ;
  int bit_physics = (physics &  0x01) ;
  int bit_2Dbhabha = (bhabhaveto & 0x01) ;
  int bit_bhabhatype = (Bhabhatype & 0x3FFF);
  int bit_etot = (((int)etot) & 0x1FFF) ;
  int bit_elow = (elow & 0x01);
  int bit_ehigh = (ehigh & 0x01) ;
  int bit_elum = (elum & 0x01) ;
  int bit_ICN = (ICN & 0x7F) ;
  int bit_BGVeto = (BGVeto & 0x07) ;
  int bit_ClusterOverflow = (ClusterOverflow & 0x01);
  int bit_3Dbhabha = (bhabha3D & 0x01);
 
  int bit_lowmulti1 = lowmultibit & 0x0FFF;
  int bit_lowmulti2 = (lowmultibit >>= 12) & 0x3;
  int bit_3DBhabha_sel = bhabha3D_sel & 0x01;
  int bit_mumu = mumubit & 0x01;
  int bit_prescale = prescale & 0x01;
  int bit_burst = burst & 0x01;
  int bit_clkcc = 0;  // 4 bits for revo counter (set to be 0 in tsim)
  int bit_eventtimingqualityflag = EventTimingQualityFlag & 0x03;
  int bit_bhabha3dvetointrackflag = 0;
  if (bhabha3D == 1) {
    bit_bhabha3dvetointrackflag = bhabha3DVetoInTrackFlag & 0x01;
  }
  int bit_bhabha3dselectionthetaflag = 0;
  if (bhabha3D_sel == 1) {
    bit_bhabha3dselectionthetaflag = bhabha3DSelectionThetaFlag & 0x03;
  }
  int bit_taub2bflag = taub2bFlag & 0x01;
  int bit_taub2b2flag = taub2b2Flag & 0x01;
  int bit_taub2b3flag = taub2b3Flag & 0x01;
 
  int bit_hie_bhav_hie123 = bit_hie_bhav & 0x7;      // for hie1, hie2, hie3
  int bit_hie_bhav_hie4   = (bit_hie_bhav >> 3) & 0x1; // for hie4
 
  m_Triggerbit[2] |= bit_hie_bhav_hie4;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_taub2b3flag;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_taub2b2flag;
  m_Triggerbit[2] <<= 3;
  m_Triggerbit[2] |= bit_hie_bhav_hie123;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_taub2bflag;
  m_Triggerbit[2] <<= 2;
  m_Triggerbit[2] |= bit_bhabha3dselectionthetaflag;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_bhabha3dvetointrackflag;
  m_Triggerbit[2] <<= 2;
  m_Triggerbit[2] |= bit_eventtimingqualityflag;
  m_Triggerbit[2] <<= 4;
  m_Triggerbit[2] |= bit_clkcc;
  m_Triggerbit[2] <<= 2;
  m_Triggerbit[2] |= bit_lowmulti2;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_burst;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_prescale;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_mumu;
  m_Triggerbit[2] <<= 1;
  m_Triggerbit[2] |= bit_3DBhabha_sel;
  m_Triggerbit[2] <<= 10;
  m_Triggerbit[2] |= ((bit_lowmulti1) >> 2) & 0x3FF;
 
  m_Triggerbit[1] |= (bit_lowmulti1 & 0x03);
  m_Triggerbit[1] <<= 1;
  m_Triggerbit[1] |= bit_3Dbhabha;
  m_Triggerbit[1] <<= 1;
  m_Triggerbit[1] |= bit_ClusterOverflow;
  m_Triggerbit[1] <<= 3;
  m_Triggerbit[1] |= bit_BGVeto;
  m_Triggerbit[1] <<= 7;
  m_Triggerbit[1] |= bit_ICN;
  m_Triggerbit[1] <<= 1;
  m_Triggerbit[1] |= bit_elum;
  m_Triggerbit[1] <<= 1;
  m_Triggerbit[1] |= bit_ehigh;
  m_Triggerbit[1] <<= 1;
  m_Triggerbit[1] |= bit_elow;
  m_Triggerbit[1] <<= 13;
  m_Triggerbit[1] |= bit_etot;
  m_Triggerbit[1] <<= 2;
  m_Triggerbit[1] |= ((bit_bhabhatype >> 12) & 0x03);
 
  m_Triggerbit[0] |= (bit_bhabhatype & 0x0FFF);
  m_Triggerbit[0] <<= 1;
  m_Triggerbit[0] |= bit_2Dbhabha;
  m_Triggerbit[0] <<= 1;
  m_Triggerbit[0] |= bit_physics;
  m_Triggerbit[0] <<= 3;
  m_Triggerbit[0] |= bit_TimingSource;
  m_Triggerbit[0] <<= 7;
  m_Triggerbit[0] |= bit_RevoFAM;
  m_Triggerbit[0] <<= 7;
  m_Triggerbit[0] |= bit_Timing;
  m_Triggerbit[0] <<= 1;
  m_Triggerbit[0] |= bit_hit;
 
  // int tmp = (m_Triggerbit[2] >> 24) & 0x03;
  // printf("%5i %5i %i\n", bit_bhabha3dselectionthetaflag, tmp, m_Triggerbit[2]);
}

name()

std::string name

(

void

)

const

returns name.

Definition at line 124 of file TrgEclMaster.cc.

{
  return "TrgEclMaster";
}

set2DBhabhaThreshold()

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

inline

set 2D Bhabha Energy Threshold

Definition at line 113 of file TrgEclMaster.h.

    {
      m_2DBhabhaThresholdFWD = i2DBhabhaThresholdFWD;
      m_2DBhabhaThresholdBWD = i2DBhabhaThresholdBWD;
    }

set3DBhabhaAddAngleCut()

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

inline

set 3D Bhabha addtion Angle selection

Definition at line 148 of file TrgEclMaster.h.

    {
      m_3DBhabhaAddAngleCut = i3DBhabhaAddAngleCut;
    }

set3DBhabhaSelectionAngle()

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

inline

set 3D selection Bhabha Energy Angle

Definition at line 131 of file TrgEclMaster.h.

    {
      m_3DBhabhaSelectionAngle = i3DBhabhaSelectionAngle;
    };

set3DBhabhaSelectionThreshold()

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

inline

set 3D selection Bhabha Energy Threshold

Definition at line 120 of file TrgEclMaster.h.

    {
      m_3DBhabhaSelectionThreshold = i3DBhabhaSelectionThreshold;
    };

set3DBhabhaVetoAngle()

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

inline

set 3D veto Bhabha Energy Angle

Definition at line 136 of file TrgEclMaster.h.

    {
      m_3DBhabhaVetoAngle = i3DBhabhaVetoAngle;
    };

set3DBhabhaVetoInTrackThetaRegion()

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

inline

set theta ID region (low and high) of 3DBhabhaVeto InTrack for gg selection

Definition at line 214 of file TrgEclMaster.h.

    {
      m_3DBhabhaVetoInTrackThetaRegion = i3DBhabhaVetoInTrackThetaRegion;
    }

set3DBhabhaVetoThreshold()

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

inline

set 3D veto Bhabha Energy Threshold

Definition at line 125 of file TrgEclMaster.h.

    {
      m_3DBhabhaVetoThreshold = i3DBhabhaVetoThreshold;
    };

setBhabhaMethod()

void setBhabhaMethod ( int  bhabha )

inline

Set Bhabha.

Definition at line 82 of file TrgEclMaster.h.

{m_Bhabha = bhabha;}

setClusterLimit()

void setClusterLimit ( int  limit )

inline

Set the limit # of Cluster.

Definition at line 80 of file TrgEclMaster.h.

{m_ClusterLimit = limit;}

setClusterMethod()

void setClusterMethod ( int  cluster )

inline

Get Event timing.

Set Cluster

Definition at line 78 of file TrgEclMaster.h.

{m_Clustering = cluster;}

setECLBurstThreshold()

void setECLBurstThreshold ( int  ECLBurstThreshold )

inline

set mumu bit Threshold

Definition at line 199 of file TrgEclMaster.h.

    {
      m_ECLBurstThreshold = (double) ECLBurstThreshold;
    }

setEventTimingMethod()

void setEventTimingMethod ( int  EventTiming )

inline

Set Cluster.

Definition at line 84 of file TrgEclMaster.h.

{m_EventTiming = EventTiming;}

setEventTimingQualityThresholds()

void setEventTimingQualityThresholds ( const std::vector< double > &  iEventTimingQualityThresholds )

inline

set energy threshold(low and high) of event timing quality flag (GeV)

Definition at line 219 of file TrgEclMaster.h.

    {
      m_EventTimingQualityThresholds = iEventTimingQualityThresholds;
    }

setLowMultiplicityThreshold()

void setLowMultiplicityThreshold ( const std::vector< double > &  iLowMultiThreshold )

inline

set Low Multiplicity Threshold

Definition at line 209 of file TrgEclMaster.h.

    {
      m_LowMultiThreshold = iLowMultiThreshold;
    }

setmumuAngle()

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

inline

set mumu bit Angle selection

Definition at line 143 of file TrgEclMaster.h.

    {
      m_mumuAngle = imumuAngle;
    }

setmumuThreshold()

void setmumuThreshold ( int  mumuThreshold )

inline

set mumu bit Threshold

Definition at line 141 of file TrgEclMaster.h.

{m_mumuThreshold = mumuThreshold; }

setn300MeVClusterThreshold()

void setn300MeVClusterThreshold ( int  n300MeVCluster )

inline

set the number of cluster exceeding 300 MeV

Definition at line 194 of file TrgEclMaster.h.

    {
      m_n300MeVCluster = n300MeVCluster;
    }

setNofTopTC()

void setNofTopTC ( int  noftoptc )

inline

set # of considered TC in energy weighted Timing method

Definition at line 90 of file TrgEclMaster.h.

{m_NofTopTC = noftoptc;}

setOverlapWindow()

void setOverlapWindow ( int  overlapwindow )

inline

Set Trigger Decision overlap window size.

Definition at line 88 of file TrgEclMaster.h.

{m_OverlapWindow = overlapwindow;}

setRS()

void setRS	(	std::vector< int >	TCId,
		std::vector< double >	TCHit,
		std::vector< double > &	phiringsum,
		std::vector< std::vector< double > > &	thetaringsum
	)

ECL bit information for GDL.

ECL bit information for GDL Set Phi Ring Sum

Definition at line 1069 of file TrgEclMaster.cc.

{
  //
  //
  // TC Phi & Theta ring sum
  //
  //----------------------------------------
  //               FW    BR     BW   total
  //----------------------------------------
  //TC phi ring     3    12     2      17
  //ID             1-3  4-15  16-17
  //TC Theta ring   32    36    32
  //ID             1-32  1-36  1-32
 
  //----------------------------------------
  //
  thetaringsum.resize(3, std::vector<double>(36, 0.));
  phiringsum.resize(17);
  const int size_hit = TCHit.size();
  for (int iHit = 0; iHit < size_hit; iHit++) {
    int iTCId = TCId[iHit] - 1;
    if (TCHit[iHit] > 0) {
      int iTCThetaId = m_obj_map->getTCThetaIdFromTCId(iTCId + 1) - 1 ;
      int iTCPhiId   = m_obj_map->getTCPhiIdFromTCId(iTCId + 1) - 1 ;
      phiringsum[iTCThetaId] += TCHit[iHit];
      if (iTCThetaId < 3) {
        //fwd
        if (iTCThetaId != 0) {
          thetaringsum[0][iTCPhiId] += TCHit[iHit];
        }
      } else if (iTCThetaId < 15) {
        //barrel
        thetaringsum[1][iTCPhiId] += TCHit[iHit];
      } else {
        //bwd
        thetaringsum[2][iTCPhiId] += TCHit[iHit];
      }
 
    }
 
  }
 
}

setTaub2b2Cut()

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

inline

set taub2b2 cut

Definition at line 170 of file TrgEclMaster.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 181 of file TrgEclMaster.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 tau b2b 2 cluster angle cut

Definition at line 153 of file TrgEclMaster.h.

    {
      m_taub2bAngleCut = itaub2bAngleCut;
    }

setTaub2bClusterECut()

void setTaub2bClusterECut ( double  itaub2bClusterECut1, double  itaub2bClusterECut2  )

inline

set tau b2b 1Cluster energy cut

Definition at line 163 of file TrgEclMaster.h.

    {
      m_taub2bClusterECut1 = itaub2bClusterECut1;
      m_taub2bClusterECut2 = itaub2bClusterECut2;
    }

setTaub2bEtotCut()

void setTaub2bEtotCut ( double  itaub2bEtotCut )

inline

set tau b2b total energy cut

Definition at line 158 of file TrgEclMaster.h.

    {
      m_taub2bEtotCut = itaub2bEtotCut;
    }

setTimeWindow()

void setTimeWindow ( int  timewindow )

inline

Set Trigger Decision window size.

Definition at line 86 of file TrgEclMaster.h.

{m_TimeWindow = timewindow;}

setTotalEnergy()

double setTotalEnergy

(

std::vector< double >

phisum

)

Set Total Energy.

Definition at line 1480 of file TrgEclMaster.cc.

{
 
  double E_phys = 0;
  for (int iii = 0; iii <= 16; iii++) {
    if (iii > 0 && iii < 15) {E_phys += phisum[iii];}
  }
  return E_phys;
}

setTotalEnergyThreshold()

void setTotalEnergyThreshold ( const std::vector< double > &  iTotalEnergy )

inline

set Total Energy Theshold (low, high, lum)

Definition at line 204 of file TrgEclMaster.h.

    {
      m_TotalEnergy = iTotalEnergy;
    }

simulate01()

void simulate01

(

int

m_nEvent

)

simulates ECL trigger for Global Cosmic data

Definition at line 195 of file TrgEclMaster.cc.

{
  //  TrgEclFAM* obj_trgeclfam = new TrgEclFAM();
  //   obj_trgeclfam->setup(m_nEvent, 1);
  // setPRS(obj_trgeclfam);
  //
  //----------
  // TC Etot
  //----------
  //
  // Energy sum of forward and barrel except for extreme forward
  // so Etot is sum of "phi ring ID" = 1-14
  // Etot threshold are
  // 1.0 GeV for Etot01
  // 0.5 GeV for Etot02nt
  // 3.0 GeV for Etot03
  //
  // Read FAM Output
  m_TCTiming.clear();
  m_TCEnergy.clear();
  m_TCBeamBkgTag.clear();
  m_TCEnergy.resize(576);
  m_TCTiming.resize(576);
  m_TCBeamBkgTag.resize(576);
 
  StoreArray<TRGECLHit> trgeclHitArray;
  for (int ii = 0; ii < trgeclHitArray.getEntries(); ii++) {
 
    TRGECLHit* aTRGECLHit = trgeclHitArray[ii];
    int iTCID = (aTRGECLHit->getTCId() - 1);
    double HitTiming  = aTRGECLHit->getTimeAve();
    double HitEnergy  = aTRGECLHit->getEnergyDep();
    double HitBeamBkg = aTRGECLHit->getBeamBkgTag();
 
    m_TCTiming[iTCID].push_back(HitTiming);
    m_TCEnergy[iTCID].push_back(HitEnergy);
    m_TCBeamBkgTag[iTCID].push_back(HitBeamBkg);
  }
  //
  //
  int nBin = 8000 / (m_TimeWindow / 2) ; //8000/125
  /* cppcheck-suppress variableScope */
  double WindowStart;
  /* cppcheck-suppress variableScope */
  double WindowEnd;
  double fluctuation = ((gRandom->Uniform(-1, 0))) * 125;
  /* cppcheck-suppress variableScope */
  double check_window_start;
  /* cppcheck-suppress variableScope */
  double check_window_end;
 
  for (int iBin = 0 ; iBin < nBin; iBin ++) {
 
    check_window_start = iBin * (m_TimeWindow / 3) + fluctuation - 4000;
    WindowStart = check_window_start;
    check_window_end  = check_window_start + m_TimeWindow / 3;
    WindowEnd = WindowStart + m_TimeWindow;
    m_HitTCId.clear();
    m_TCHitTiming.clear();
    m_TCHitEnergy.clear();
    m_TCHitBeamBkgTag.clear();
 
    // prepare TC Hit in time window --
    for (int iTCId = 0; iTCId < 576; iTCId++) {
      const int hitsize =  m_TCTiming[iTCId].size();
      for (int ihit = 0; ihit < hitsize; ihit++) {
        if (m_TCTiming[iTCId][ihit] > check_window_start &&
            m_TCTiming[iTCId][ihit] < check_window_end) {
          m_HitTCId.push_back(iTCId + 1);
 
        }
      }
    }
    if (m_HitTCId.size() == 0) {continue;}
    else {
      m_HitTCId.clear();
      for (int iTCId = 0; iTCId < 576; iTCId++) {
        const int hitsize =  m_TCTiming[iTCId].size();
        for (int ihit = 0; ihit < hitsize; ihit++) {
          if (m_TCTiming[iTCId][ihit] > WindowStart &&
              m_TCTiming[iTCId][ihit] < WindowEnd) {
            m_HitTCId.push_back(iTCId + 1);
            m_TCHitTiming.push_back(m_TCTiming[iTCId][ihit]);
            m_TCHitEnergy.push_back(m_TCEnergy[iTCId][ihit]);
          }
        }
      }
      iBin  = iBin + 2;
    }
    int noftchit = m_HitTCId.size();
    if (noftchit == 0) {continue;}
 
    double eventtiming = 0;
    // Get EventTiming
    m_obj_timing->Setup(m_HitTCId, m_TCHitEnergy, m_TCHitTiming);
    m_obj_timing->SetNofTopTC(m_NofTopTC);
    //    m_obj_timing->setEventTimingQualityThresholds(m_EventTimingQualityThresholds);
 
    eventtiming = m_obj_timing->GetEventTiming(m_EventTiming);
    int timingsource = m_obj_timing->GetTimingSource();
    int EventTimingQualityFlag = m_obj_timing->getEventTimingQualityFlag();
    //--------------------------------------------------
    // Ring sum and Total Energy Sum
    //-------------------------------------------------
    std::vector<std::vector<double>>  thetaringsum;
    std::vector<double>  phiringsum;
 
    thetaringsum.clear();
    phiringsum.clear();
    thetaringsum.resize(3, std::vector<double>(36, 0));
    phiringsum.resize(17, 0);
    setRS(m_HitTCId, m_TCHitEnergy, phiringsum, thetaringsum);
 
    double E_phys = 0;
    double E_total = 0;
    int E_burst = 0;
    for (int iii = 0; iii <= 16; iii++) {
      if (iii > 0 && iii < 15) {E_phys += phiringsum[iii];}
      E_total += phiringsum[iii];
    }
    if (E_total == 0) {continue;}
    int ELow = 0, EHigh = 0, ELum = 0;
 
    if (E_total > m_ECLBurstThreshold / 10) {
      E_burst = 0x01;
    }
    if (E_phys > m_TotalEnergy[0] / 10) { // GeV
      ELow = 0x01;
    }
    if (E_phys > m_TotalEnergy[1] / 10) { // GeV
      EHigh = 0x01;
    }
    if (E_phys > m_TotalEnergy[2] / 10) { // GeV
      ELum = 0x01;
    }
    //--------------
    // Clustering
    //--------------
    //
    // TrgEclCluster m_obj_cluster;
    m_obj_cluster->setClusteringMethod(m_Clustering);
    m_obj_cluster->setClusterLimit(m_ClusterLimit);
    m_obj_cluster->setEventId(m_nEvent);
    m_obj_cluster->setICN(m_HitTCId, m_TCHitEnergy, m_TCHitTiming);// Clustering
    m_obj_cluster->setICN(m_HitTCId); // Belle Cluster Counting
 
    int icn    = m_obj_cluster->getICNFwBr();
    int icnfwd = m_obj_cluster->getICNSub(0);
    int icnbr  = m_obj_cluster->getICNSub(1);
    int icnbwd = m_obj_cluster->getICNSub(2);
    //--------------
    // Low Multiplicity bit
    //--------------
    std::vector<double> ClusterTiming;
    std::vector<double> ClusterEnergy;
    std::vector<int> MaxTCId;
    ClusterTiming.clear();
    ClusterEnergy.clear();
    MaxTCId.clear();
    StoreArray<TRGECLCluster> trgeclClusterArray;
    for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
      TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
      int maxTCId    = aTRGECLCluster->getMaxTCId();
      double clusterenergy  = aTRGECLCluster->getEnergyDep();
      double clustertiming  =  aTRGECLCluster->getTimeAve();
      ClusterTiming.push_back(clustertiming);
      ClusterEnergy.push_back(clusterenergy);
      MaxTCId.push_back(maxTCId);
    }
    //    const int ncluster = ClusterEnergy.size();
 
    makeLowMultiTriggerBit(MaxTCId, ClusterEnergy);
    //--------------
    // Bhabha veto
    //--------------
    m_obj_bhabha->set2DBhabhaThreshold(m_2DBhabhaThresholdFWD, m_2DBhabhaThresholdBWD);
    m_obj_bhabha->set3DBhabhaSelectionThreshold(m_3DBhabhaSelectionThreshold);
    m_obj_bhabha->set3DBhabhaVetoThreshold(m_3DBhabhaVetoThreshold);
    m_obj_bhabha->set3DBhabhaSelectionAngle(m_3DBhabhaSelectionAngle);
    m_obj_bhabha->set3DBhabhaVetoAngle(m_3DBhabhaVetoAngle);
    m_obj_bhabha->setmumuThreshold(m_mumuThreshold);
    m_obj_bhabha->setmumuAngle(m_mumuAngle);
    m_obj_bhabha->set3DBhabhaAddAngleCut(m_3DBhabhaAddAngleCut);
    m_obj_bhabha->setTaub2bAngleCut(m_taub2bAngleCut);
    m_obj_bhabha->setTaub2bEtotCut(m_taub2bEtotCut);
    m_obj_bhabha->setTaub2bClusterECut(m_taub2bClusterECut1,
                                       m_taub2bClusterECut2);
    m_obj_bhabha->setTaub2b2Cut(m_taub2b2AngleCut,
                                m_taub2b2EtotCut,
                                m_taub2b2CLEEndcapCut,
                                m_taub2b2CLECut);
    // for taub2b3
    m_obj_bhabha->setTaub2b3Cut(m_taub2b3AngleCut,
                                m_taub2b3EtotCut,
                                m_taub2b3CLEb2bCut,
                                m_taub2b3CLELowCut,
                                m_taub2b3CLEHighCut);
 
    std::vector<double> vct_bhabha;
    vct_bhabha.clear();
    int bhabha2D = 0 ;
    int bhabha3D_veto = 0 ;
    int bhabha3D_sel = 0;
    int mumu = 0;
    int bhabha3DVetoInTrackFlag = 0;
    int bhabha3DSelectionThetaFlag = 0;
    int taub2bFlag = 0;
    int taub2b2Flag = 0;
    int taub2b3Flag = 0;
 
    bool b_2Dbhabha = m_obj_bhabha->GetBhabha00(phiringsum);
    vct_bhabha = m_obj_bhabha->GetBhabhaComb();
    if (b_2Dbhabha && (icn < 4)) {bhabha2D = 1;}
    bool b_3Dbhabha =  m_obj_bhabha->GetBhabha01();
    if (b_3Dbhabha) {bhabha3D_veto = 1;}
    bool b_3Dbhabha_sel =  m_obj_bhabha->GetBhabha02();
    if (b_3Dbhabha_sel) {bhabha3D_sel = 1;}
    bool b_mumu =  m_obj_bhabha->Getmumu();
    if (b_mumu) {mumu = 1;}
    bhabha3DVetoInTrackFlag = m_obj_bhabha->get3DBhabhaVetoInTrackFlag();
    bhabha3DSelectionThetaFlag = m_obj_bhabha->get3DBhabhaSelectionThetaFlag();
    taub2bFlag  = (m_obj_bhabha->GetTaub2b(E_total)) ? 1 : 0;
    taub2b2Flag = (m_obj_bhabha->GetTaub2b2(E_total)) ? 1 : 0;
    taub2b3Flag = (m_obj_bhabha->GetTaub2b3(E_total)) ? 1 : 0;
 
    //------------------------
    // Beam Background veto (Old cosmic veto)
    //------------------------
    int beambkgtag = 0;
    beambkgtag = m_obj_beambkg->GetBeamBkg(thetaringsum);
 
    int bhabhaprescale = 0;
    if (m_PrescaleFactor == m_PrescaleCounter) {
      bhabhaprescale = 1;
      m_PrescaleCounter = 0;
    } else if (m_PrescaleFactor > m_PrescaleCounter) {
      m_PrescaleCounter ++;
    }
    //-------------
    // Make ECL  Trigger Bit
    //-------------
    int hit = 1; // hit or not
    int Timing = (int)(eventtiming + 0.5);
    int RevoFAM = 0;
    int TimingSource = m_obj_timing->GetTimingSource(); //  FWD(0), Barrel(0), Backward(0);
    int etot = (int)(E_phys * 1000 + 0.5); // total Energy in theta ID [2~15]
    int physics = 0;
    if ((etot > 1000 || icn > 3) && !(bhabha2D == 1)) {physics = 1;}
    std::vector<int> bhabhabit;
    bhabhabit.clear();
    int bhabhabitsize = vct_bhabha.size();
    for (int ibhabha = 0; ibhabha < bhabhabitsize; ibhabha++) {
      bhabhabit.push_back((int)vct_bhabha[ibhabha]);
    }
    int ClusterOverflow = m_obj_cluster->getNofExceedCluster();
    int flagoverflow = 0;
    if (ClusterOverflow > 0) {
      flagoverflow = 1;
    }
 
    makeTriggerBit(hit, Timing, 0, timingsource, E_phys,
                   bhabha2D, physics, bhabhabit, icn, beambkgtag,
                   flagoverflow, bhabha3D_veto, m_Lowmultibit,
                   bhabha3D_sel, mumu, bhabhaprescale, E_burst,
                   EventTimingQualityFlag,
                   bhabha3DVetoInTrackFlag,
                   bhabha3DSelectionThetaFlag,
                   taub2bFlag, 0, taub2b2Flag, taub2b3Flag);
 
    int m_hitEneNum = 0;
    StoreArray<TRGECLTrg> trgEcltrgArray;
    trgEcltrgArray.appendNew();
    m_hitEneNum = trgEcltrgArray.getEntries() - 1;
    //-------------
    // Store
    //-------------
    trgEcltrgArray[m_hitEneNum]->setEventId(m_nEvent);
    trgEcltrgArray[m_hitEneNum]->setPRS01(phiringsum[0]);
    trgEcltrgArray[m_hitEneNum]->setPRS02(phiringsum[1]);
    trgEcltrgArray[m_hitEneNum]->setPRS03(phiringsum[2]);
    trgEcltrgArray[m_hitEneNum]->setPRS04(phiringsum[3]);
    trgEcltrgArray[m_hitEneNum]->setPRS05(phiringsum[4]);
    trgEcltrgArray[m_hitEneNum]->setPRS06(phiringsum[5]);
    trgEcltrgArray[m_hitEneNum]->setPRS07(phiringsum[6]);
    trgEcltrgArray[m_hitEneNum]->setPRS08(phiringsum[7]);
    trgEcltrgArray[m_hitEneNum]->setPRS09(phiringsum[8]);
    trgEcltrgArray[m_hitEneNum]->setPRS10(phiringsum[9]);
    trgEcltrgArray[m_hitEneNum]->setPRS11(phiringsum[10]);
    trgEcltrgArray[m_hitEneNum]->setPRS12(phiringsum[11]);
    trgEcltrgArray[m_hitEneNum]->setPRS13(phiringsum[12]);
    trgEcltrgArray[m_hitEneNum]->setPRS14(phiringsum[13]);
    trgEcltrgArray[m_hitEneNum]->setPRS15(phiringsum[14]);
    trgEcltrgArray[m_hitEneNum]->setPRS16(phiringsum[15]);
    trgEcltrgArray[m_hitEneNum]->setPRS17(phiringsum[16]);
    //
    trgEcltrgArray[m_hitEneNum]->setEtot(E_phys);
    trgEcltrgArray[m_hitEneNum]->setNofTCHit(noftchit);
    //
    trgEcltrgArray[m_hitEneNum]->setBhabha01(vct_bhabha[0]);
    trgEcltrgArray[m_hitEneNum]->setBhabha02(vct_bhabha[1]);
    trgEcltrgArray[m_hitEneNum]->setBhabha03(vct_bhabha[2]);
    trgEcltrgArray[m_hitEneNum]->setBhabha04(vct_bhabha[3]);
    trgEcltrgArray[m_hitEneNum]->setBhabha05(vct_bhabha[4]);
    trgEcltrgArray[m_hitEneNum]->setBhabha06(vct_bhabha[5]);
    trgEcltrgArray[m_hitEneNum]->setBhabha07(vct_bhabha[6]);
    trgEcltrgArray[m_hitEneNum]->setBhabha08(vct_bhabha[7]);
    trgEcltrgArray[m_hitEneNum]->setBhabha09(vct_bhabha[8]);
    trgEcltrgArray[m_hitEneNum]->setBhabha10(vct_bhabha[9]);
    trgEcltrgArray[m_hitEneNum]->setBhabha11(vct_bhabha[10]);
    trgEcltrgArray[m_hitEneNum]->setBhabha12(vct_bhabha[11]);
    trgEcltrgArray[m_hitEneNum]->setBhabha13(vct_bhabha[12]);
    trgEcltrgArray[m_hitEneNum]->setBhabha14(vct_bhabha[13]);
    trgEcltrgArray[m_hitEneNum]->setBhabha15(vct_bhabha[14]);
    trgEcltrgArray[m_hitEneNum]->setBhabha16(vct_bhabha[15]);
    trgEcltrgArray[m_hitEneNum]->setBhabha17(vct_bhabha[16]);
    trgEcltrgArray[m_hitEneNum]->setBhabha18(vct_bhabha[17]);
    //
    trgEcltrgArray[m_hitEneNum]->setICN(icn);
    trgEcltrgArray[m_hitEneNum]->setICNFw(icnfwd);
    trgEcltrgArray[m_hitEneNum]->setICNBr(icnbr);
    trgEcltrgArray[m_hitEneNum]->setICNBw(icnbwd);
    //
    trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[0], 0);
    trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[1], 1);
    trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[2], 2);
    trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[3], 3);
 
    trgEcltrgArray[m_hitEneNum]->setBhabhaVeto(bhabha2D);
    trgEcltrgArray[m_hitEneNum]->setBeamBkgVeto(beambkgtag);
    trgEcltrgArray[m_hitEneNum]->setEventTiming(eventtiming);
 
    trgEcltrgArray[m_hitEneNum]->setHit(hit);
    trgEcltrgArray[m_hitEneNum]->setRevoclk(RevoFAM);
    trgEcltrgArray[m_hitEneNum]->setTimingSource(TimingSource);
    trgEcltrgArray[m_hitEneNum]->setPhysics(physics) ;
    trgEcltrgArray[m_hitEneNum]->set2DBhabha(bhabha2D);
    trgEcltrgArray[m_hitEneNum]->set3DBhabha(bhabha3D_veto);
    trgEcltrgArray[m_hitEneNum]->set3DBhabhaSel(bhabha3D_sel);
    trgEcltrgArray[m_hitEneNum]->setmumuBit(mumu);
    trgEcltrgArray[m_hitEneNum]->setBhabhaPrescaleBit(bhabhaprescale);
 
    trgEcltrgArray[m_hitEneNum]->setELow(ELow)  ;
    trgEcltrgArray[m_hitEneNum]->setEHihg(EHigh);
    trgEcltrgArray[m_hitEneNum]->setELum(ELum)  ;
    trgEcltrgArray[m_hitEneNum]->setClusterOverflow(ClusterOverflow) ;
    trgEcltrgArray[m_hitEneNum]->setLowMultiBit(m_Lowmultibit);
  }
 
  return;
}

simulate02()

void simulate02

(

int

m_nEvent

)

simulates ECL trigger for Data Analysis

Definition at line 549 of file TrgEclMaster.cc.

{
 
  //  TrgEclFAM* obj_trgeclfam = new TrgEclFAM();
  //   obj_trgeclfam->setup(m_nEvent, 1);
  // setPRS(obj_trgeclfam);
  //
  //----------
  // TC Etot
  //----------
  //
  // Energy sum of forward and barrel except for extreme forward
  // so Etot is sum of "phi ring ID" = 1-14
  // Etot threshold are
  // 1.0 GeV for Etot01
  // 0.5 GeV for Etot02nt
  // 3.0 GeV for Etot03
  //
  // Read FAM Output
  m_TCTiming.clear();
  m_TCEnergy.clear();
  m_TCBeamBkgTag.clear();
  m_TCEnergy.resize(576);
  m_TCTiming.resize(576);
  m_TCBeamBkgTag.resize(576);
 
  StoreArray<TRGECLHit> trgeclHitArray;
  for (int ii = 0; ii < trgeclHitArray.getEntries(); ii++) {
 
    TRGECLHit* aTRGECLHit = trgeclHitArray[ii];
    int iTCID = (aTRGECLHit->getTCId() - 1);
    double HitTiming  = aTRGECLHit->getTimeAve();
    double HitEnergy  = aTRGECLHit->getEnergyDep();
    double HitBeamBkg = aTRGECLHit->getBeamBkgTag();
 
    m_TCTiming[iTCID].push_back(HitTiming);
    m_TCEnergy[iTCID].push_back(HitEnergy);
    m_TCBeamBkgTag[iTCID].push_back(HitBeamBkg);
  }
  //
  //
  int nBin = 2 * 8000 / m_TimeWindow ;
  double WindowStart = 0;
  double WindowEnd = 0;
  double fluctuation = ((gRandom->Uniform(-1, 0))) * 125;
 
  int startBin = nBin / 2 - 1; //start previous bin near 0s
 
  int endBin = nBin / 2 + 1; //start next bin near 0s
 
  if (m_EventTiming == 0) {
    m_TimeWindow = 500;
    m_OverlapWindow = 0;
  }
 
  double maxE = 0;
  int max_bin = 0;
 
  for (int iBin = startBin ; iBin <= endBin; iBin ++) {
    WindowStart = iBin * (m_TimeWindow - m_OverlapWindow) + fluctuation - 4000;
 
    if (iBin == 0) {WindowStart = - 4000 + fluctuation;}
    WindowEnd = WindowStart + m_TimeWindow;
 
    double totalE = 0;
    // prepare TC Hit in time window --
    for (int iTCId = 0; iTCId < 576; iTCId++) {
      const int hitsize =  m_TCTiming[iTCId].size();
      for (int ihit = 0; ihit < hitsize; ihit++) {
        if (m_TCTiming[iTCId][ihit] > WindowStart &&
            m_TCTiming[iTCId][ihit] < WindowEnd) {
          totalE += m_TCEnergy[iTCId][ihit] ;
        }
      }
    }
    if (totalE == 0) {continue;}
    if (maxE < totalE) { //select the bin having the highest total energy
      maxE = totalE;
      max_bin = iBin;
    }
  }
 
  WindowStart = max_bin * (m_TimeWindow - m_OverlapWindow) + fluctuation - 4000;
  if (max_bin == 0) {WindowStart = - 4000 + fluctuation;}
  WindowEnd = WindowStart + m_TimeWindow;
 
  m_HitTCId.clear();
  m_TCHitTiming.clear();
  m_TCHitEnergy.clear();
  m_TCHitBeamBkgTag.clear();
 
  // prepare TC Hit in time window --
  for (int iTCId = 0; iTCId < 576; iTCId++) {
    const int hitsize =  m_TCTiming[iTCId].size();
    for (int ihit = 0; ihit < hitsize; ihit++) {
      if (m_TCTiming[iTCId][ihit] > WindowStart &&
          m_TCTiming[iTCId][ihit] < WindowEnd) {
        m_HitTCId.push_back(iTCId + 1);
        m_TCHitTiming.push_back(m_TCTiming[iTCId][ihit]);
        m_TCHitEnergy.push_back(m_TCEnergy[iTCId][ihit]);
        m_TCHitBeamBkgTag.push_back(m_TCBeamBkgTag[iTCId][ihit]);
      }
    }
  }
 
  int noftchit = m_HitTCId.size();
  if (noftchit == 0) { return;}
 
  double eventtiming = 0;
  // Get EventTiming
  m_obj_timing->Setup(m_HitTCId, m_TCHitEnergy, m_TCHitTiming);
  m_obj_timing->SetNofTopTC(m_NofTopTC);
  m_obj_timing->setEventTimingQualityThresholds(m_EventTimingQualityThresholds);
 
  eventtiming = m_obj_timing->GetEventTiming(m_EventTiming);
  int timingsource = m_obj_timing->GetTimingSource();
 
  int EventTimingQualityFlag = m_obj_timing->getEventTimingQualityFlag();
  //--------------------------------------------------
  // Ring sum and Total Energy Sum
  //-------------------------------------------------
  std::vector<std::vector<double>>  thetaringsum;
  std::vector<double>  phiringsum;
 
  thetaringsum.clear();
  phiringsum.clear();
  thetaringsum.resize(3, std::vector<double>(36, 0));
  phiringsum.resize(17, 0);
  setRS(m_HitTCId, m_TCHitEnergy, phiringsum, thetaringsum);
 
  //double E_br;  //variable not used
  //double E_fwd; //variable not used
  //double E_bwd; //variable not used
  double E_phys = 0;
  double E_total = 0;
  int E_burst = 0;
 
  for (int iii = 0; iii <= 16; iii++) {
    if (iii > 0 && iii < 15) {E_phys += phiringsum[iii];}
    //if (iii < 3) {E_fwd += phiringsum[iii];} //TODO variable not used, should be?
    //if (iii > 2 && iii < 15) {E_br += phiringsum[iii];} //TODO not used, should be?
    //if (iii > 14) {E_bwd += phiringsum[iii];} //TODO not used, should be?
    E_total += phiringsum[iii];
  }
 
  if (E_total == 0) {return;}
  int ELow = 0, EHigh = 0, ELum = 0;
  if (E_total > m_ECLBurstThreshold) {
    E_burst = 0x01;
  }
 
  if (E_phys > m_TotalEnergy[0] / 10) { // GeV
    ELow = 0x01;
  }
  if (E_phys > m_TotalEnergy[1] / 10) { // GeV
    EHigh = 0x01;
  }
  if (E_phys > m_TotalEnergy[2] / 10) { // GeV
    ELum = 0x01;
  }
  //--------------
  // Clustering
  //--------------
  //
  // TrgEclCluster m_obj_cluster;
  m_obj_cluster->setClusteringMethod(m_Clustering);
  m_obj_cluster->setClusterLimit(m_ClusterLimit);
  m_obj_cluster->setEventId(m_nEvent);
  m_obj_cluster->setICN(m_HitTCId, m_TCHitEnergy, m_TCHitTiming);// Clustering
  m_obj_cluster->setICN(m_HitTCId); // Belle Cluster Counting
 
  int icn    = m_obj_cluster->getICNFwBr();
  int icnfwd = m_obj_cluster->getICNSub(0);
  int icnbr  = m_obj_cluster->getICNSub(1);
  int icnbwd = m_obj_cluster->getICNSub(2);
 
  int NofCluster1to17 =  m_obj_cluster->getNofCluster();
  //--------------
  // Low Multiplicity bit
  //--------------
  std::vector<double> ClusterTiming;
  std::vector<double> ClusterEnergy;
  std::vector<int> MaxTCId;
  std::vector<int> MaxThetaId;
  ClusterTiming.clear();
  ClusterEnergy.clear();
  MaxTCId.clear();
  MaxThetaId.clear();
  StoreArray<TRGECLCluster> trgeclClusterArray;
  for (int ii = 0; ii < trgeclClusterArray.getEntries(); ii++) {
    TRGECLCluster* aTRGECLCluster = trgeclClusterArray[ii];
    int maxTCId    = aTRGECLCluster->getMaxTCId();
    int maxThetaId = aTRGECLCluster->getMaxThetaId();
    double clusterenergy  = aTRGECLCluster->getEnergyDep();
    double clustertiming  =  aTRGECLCluster->getTimeAve();
    ClusterTiming.push_back(clustertiming);
    ClusterEnergy.push_back(clusterenergy);
    MaxTCId.push_back(maxTCId);
    MaxThetaId.push_back(maxThetaId);
  }
  //    int NofCluster =  m_obj_cluster->getNofCluster();
  makeLowMultiTriggerBit(MaxTCId, ClusterEnergy);
  //--------------
  // Bhabha veto (and mumu and tau b2b trigger)
  //--------------
  m_obj_bhabha->set2DBhabhaThreshold(m_2DBhabhaThresholdFWD,
                                     m_2DBhabhaThresholdBWD);
  m_obj_bhabha->set3DBhabhaSelectionThreshold(m_3DBhabhaSelectionThreshold);
  m_obj_bhabha->set3DBhabhaVetoThreshold(m_3DBhabhaVetoThreshold);
  m_obj_bhabha->set3DBhabhaSelectionAngle(m_3DBhabhaSelectionAngle);
  m_obj_bhabha->set3DBhabhaVetoAngle(m_3DBhabhaVetoAngle);
  m_obj_bhabha->setmumuThreshold(m_mumuThreshold);
  m_obj_bhabha->setmumuAngle(m_mumuAngle);
  m_obj_bhabha->set3DBhabhaAddAngleCut(m_3DBhabhaAddAngleCut);
  m_obj_bhabha->set3DBhabhaVetoInTrackThetaRegion(m_3DBhabhaVetoInTrackThetaRegion);
  m_obj_bhabha->setTaub2bAngleCut(m_taub2bAngleCut);
  m_obj_bhabha->setTaub2bEtotCut(m_taub2bEtotCut);
  m_obj_bhabha->setTaub2bClusterECut(m_taub2bClusterECut1,
                                     m_taub2bClusterECut2);
  m_obj_bhabha->setTaub2b2Cut(m_taub2b2AngleCut,
                              m_taub2b2EtotCut,
                              m_taub2b2CLEEndcapCut,
                              m_taub2b2CLECut);
  //
  m_obj_bhabha->setTaub2b3Cut(m_taub2b3AngleCut,
                              m_taub2b3EtotCut,
                              m_taub2b3CLEb2bCut,
                              m_taub2b3CLELowCut,
                              m_taub2b3CLEHighCut);
 
 
  std::vector<double> vct_bhabha;
  vct_bhabha.clear();
  int bhabha2D = 0 ;
  int bhabha3D_veto = 0 ;
  int bhabha3D_sel = 0;
  int mumu = 0;
  int bhabha3DVetoInTrackFlag = -1;
  int bhabha3DSelectionThetaFlag = -1;
  int taub2bFlag = 0;
  int taub2b2Flag = 0;
  int taub2b3Flag = 0;
 
  bool b_2Dbhabha = m_obj_bhabha->GetBhabha00(phiringsum);
  vct_bhabha = m_obj_bhabha->GetBhabhaComb();
  if (b_2Dbhabha && (icn < 4)) {bhabha2D = 1;}
  bool b_3Dbhabha =  m_obj_bhabha->GetBhabha01();
  if (b_3Dbhabha) {bhabha3D_veto = 1;}
  bool b_3Dbhabha_sel =  m_obj_bhabha->GetBhabha02();
  if (b_3Dbhabha_sel) {bhabha3D_sel = 1;}
  bool b_mumu =  m_obj_bhabha->Getmumu();
  if (b_mumu) {mumu = 1;}
  int bhabhaprescale = 0;
  if (m_PrescaleFactor == m_PrescaleCounter) {
    bhabhaprescale = 1;
    m_PrescaleCounter = 0;
  } else if (m_PrescaleFactor > m_PrescaleCounter) {
    m_PrescaleCounter ++;
  }
  bhabha3DVetoInTrackFlag    = m_obj_bhabha->get3DBhabhaVetoInTrackFlag();
  bhabha3DSelectionThetaFlag = m_obj_bhabha->get3DBhabhaSelectionThetaFlag();
  taub2bFlag  = (m_obj_bhabha->GetTaub2b(E_total)) ? 1 : 0;
  taub2b2Flag = (m_obj_bhabha->GetTaub2b2(E_total)) ? 1 : 0;
  taub2b3Flag = (m_obj_bhabha->GetTaub2b3(E_total)) ? 1 : 0;
 
  //------------------------
  // additional Bhabha veto
  //------------------------
  int bhabha_addition = m_obj_bhabha->GetBhabhaAddition();
  //------------------------
  // hie with additional Bhabha veto
  //------------------------
  int beambkgtag = 0;
  beambkgtag = m_obj_beambkg->GetBeamBkg(thetaringsum);
 
  int bit_hie_bhav = 0;
  if (E_phys > 1.0) {
    bit_hie_bhav |= (~bhabha_addition & 0x01) & (~bhabha_addition >> 4 & 0x01); // hie4
    bit_hie_bhav <<= 1;
    bit_hie_bhav |= (~bhabha_addition & 0x01) & (~bhabha_addition >> 3 & 0x01); // hie3
    bit_hie_bhav <<= 1;
    bit_hie_bhav |= (~bhabha_addition & 0x01) & (~bhabha_addition >> 2 & 0x01); // hie2
    bit_hie_bhav <<= 1;
    bit_hie_bhav |= (~bhabha_addition & 0x01) & (~bhabha_addition >> 1 & 0x01); // hie1
  }
  //-------------
  // Make ECL  Trigger Bit
  //-------------
  int hit = 1; // hit or not
  int Timing = (int)(eventtiming + 0.5);
  int RevoFAM = 0;
  int TimingSource = m_obj_timing->GetTimingSource(); //  FWD(0), Barrel(0), Backward(0);
  int etot = (int)(E_phys * 1000 + 0.5); // total Energy in theta ID [2~15]
 
  //int bhabha2D = BtoBTag ;
  int physics = 0;
  if ((etot > 1000 || icn > 3) && !(bhabha2D == 1)) {physics = 1;}
  std::vector<int> bhabhabit;
  bhabhabit.clear();
  int bhabhabitsize = vct_bhabha.size();
  for (int ibhabha = 0; ibhabha < bhabhabitsize; ibhabha++) {
    bhabhabit.push_back((int)vct_bhabha[ibhabha]);
  }
  int ClusterOverflow = m_obj_cluster->getNofExceedCluster();
  int flagoverflow = 0;
  if (ClusterOverflow > 0) {
    flagoverflow = 1;
  }
 
  makeTriggerBit(hit, Timing, 0, timingsource, E_phys, bhabha2D,
                 physics, bhabhabit, icn, beambkgtag, flagoverflow,
                 bhabha3D_veto, m_Lowmultibit, bhabha3D_sel, mumu,
                 bhabhaprescale, E_burst,
                 EventTimingQualityFlag,
                 bhabha3DVetoInTrackFlag,
                 bhabha3DSelectionThetaFlag,
                 taub2bFlag,
                 bit_hie_bhav,
                 taub2b2Flag,
                 taub2b3Flag);
 
  //----------------------------------------------------
  // ECL trigger
  //----------------------------------------------------
  // Integer GDL "Output word to GDL:"
  //  "bit0-2 = Etot1,2,3"
  //  "bit3 = Bhabha,"
  //  "bit4 = prescaled Bhabha,"
  //  "bit5-8 = ICN(3bits)+FORWARD(1bit) OR ICN(3+1 carry),"
  //  "bit9 = cosmic,"
  //  "bit10 = neutral timing trigger"
  //
  //------------------------------
  //          2         10     16
  //------------------------------
  //  1 0000000000001     1     1   Etot1
  //  2 0000000000010     2     2   Etot2
  //  3 0000000000100     4     4   Etot3
  //  4 0000000001000     8     8   Bhabha
  //  5 0000000010000    16    10   preBhabha
  //  6 0000000100000    32    20   ICN
  //  7 0000001000000    64    40   ICN
  //  8 0000010000000   128    80   ICN
  //  9 0000100000000   256   100   ForwardICN
  // 10 0001000000000   512   200   BeamBkgVeto
  // 11 0010000000000  1024   400   Timing
  //------------------------------
  // int bitEtot1       = 0x0001;
  // int bitEtot2       = 0x0002;
  // int bitEtot3       = 0x0004;
  // int bitBhabha      = 0x0008;
  // int bitPreBhabha   = 0x0010;
  // int bitForwardICN  = 0x0100;
  // int bitBeamBkgVeto = 0x0200;
  // int bitTiming      = 0x0400;
 
  // bool boolEtot[3] = {false};
  // if (E_phys > 1.0) boolEtot[1] = true;
  // bool boolBhabha = (boolBtoBTag && icn > 4);
  // bool boolPreBhabha = false;
  // bool boolForwardICN = icnfwd;
  // bool boolBeamBkgVeto = boolBeamBkgTag;
  // int bit = 0;
  // //
  // // bit 5-7
  // bit = (icn >= 7) ? 0x0007 : icn;
  // bit <<= 5;
  // // bit 0
  // bit |= boolEtot[0] ? bitEtot1 : 0;
  // // bit 1
  // bit |= boolEtot[1] ? bitEtot2 : 0;
  // // bit 2
  // bit |= boolEtot[2] ? bitEtot3 : 0;
  // // bit 3
  // bit |= boolBhabha  ? bitBhabha : 0;
  // // bit 4
  // bit |= boolPreBhabha   ? bitPreBhabha : 0;
  // // bit 8
  // bit |= boolForwardICN  ? bitForwardICN : 0;
  // // bit 9
  // bit |= boolBeamBkgVeto ? bitBeamBkgVeto : 0;
  // // bit 10
  // bit |= bitTiming;
  //
  //  printf("bit = %i \n", bit);
  //----------------------
  //   if (0){ // check bit by "binary" output
  //     int xxx = bit;
  //     int yyy = 0;
  //     int iii = 0;
  //     int ans = 0;
  //     while (xxx > 0) {
  //       yyy = xxx % 2;
  //       ans = ans + yyy * pow(10,iii);
  //       xxx = xxx / 2;
  //       iii = iii++;
  //     }
  //     printf("xxx = %i \n", ans);
  //   }
 
 
  int m_hitEneNum = 0;
  StoreArray<TRGECLTrg> trgEcltrgArray;
  trgEcltrgArray.appendNew();
  m_hitEneNum = trgEcltrgArray.getEntries() - 1;
  //----------------------------------------------------
  // Store
  //----------------------------------------------------
  trgEcltrgArray[m_hitEneNum]->setEventId(m_nEvent);
  trgEcltrgArray[m_hitEneNum]->setPRS01(phiringsum[0]);
  trgEcltrgArray[m_hitEneNum]->setPRS02(phiringsum[1]);
  trgEcltrgArray[m_hitEneNum]->setPRS03(phiringsum[2]);
  trgEcltrgArray[m_hitEneNum]->setPRS04(phiringsum[3]);
  trgEcltrgArray[m_hitEneNum]->setPRS05(phiringsum[4]);
  trgEcltrgArray[m_hitEneNum]->setPRS06(phiringsum[5]);
  trgEcltrgArray[m_hitEneNum]->setPRS07(phiringsum[6]);
  trgEcltrgArray[m_hitEneNum]->setPRS08(phiringsum[7]);
  trgEcltrgArray[m_hitEneNum]->setPRS09(phiringsum[8]);
  trgEcltrgArray[m_hitEneNum]->setPRS10(phiringsum[9]);
  trgEcltrgArray[m_hitEneNum]->setPRS11(phiringsum[10]);
  trgEcltrgArray[m_hitEneNum]->setPRS12(phiringsum[11]);
  trgEcltrgArray[m_hitEneNum]->setPRS13(phiringsum[12]);
  trgEcltrgArray[m_hitEneNum]->setPRS14(phiringsum[13]);
  trgEcltrgArray[m_hitEneNum]->setPRS15(phiringsum[14]);
  trgEcltrgArray[m_hitEneNum]->setPRS16(phiringsum[15]);
  trgEcltrgArray[m_hitEneNum]->setPRS17(phiringsum[16]);
  //
  trgEcltrgArray[m_hitEneNum]->setEtot(E_phys);
  trgEcltrgArray[m_hitEneNum]->setNofTCHit(noftchit);
  //
  trgEcltrgArray[m_hitEneNum]->setBhabha01(vct_bhabha[0]);
  trgEcltrgArray[m_hitEneNum]->setBhabha02(vct_bhabha[1]);
  trgEcltrgArray[m_hitEneNum]->setBhabha03(vct_bhabha[2]);
  trgEcltrgArray[m_hitEneNum]->setBhabha04(vct_bhabha[3]);
  trgEcltrgArray[m_hitEneNum]->setBhabha05(vct_bhabha[4]);
  trgEcltrgArray[m_hitEneNum]->setBhabha06(vct_bhabha[5]);
  trgEcltrgArray[m_hitEneNum]->setBhabha07(vct_bhabha[6]);
  trgEcltrgArray[m_hitEneNum]->setBhabha08(vct_bhabha[7]);
  trgEcltrgArray[m_hitEneNum]->setBhabha09(vct_bhabha[8]);
  trgEcltrgArray[m_hitEneNum]->setBhabha10(vct_bhabha[9]);
  trgEcltrgArray[m_hitEneNum]->setBhabha11(vct_bhabha[10]);
  trgEcltrgArray[m_hitEneNum]->setBhabha12(vct_bhabha[11]);
  trgEcltrgArray[m_hitEneNum]->setBhabha13(vct_bhabha[12]);
  trgEcltrgArray[m_hitEneNum]->setBhabha14(vct_bhabha[13]);
  trgEcltrgArray[m_hitEneNum]->setBhabha15(vct_bhabha[14]);
  trgEcltrgArray[m_hitEneNum]->setBhabha16(vct_bhabha[15]);
  trgEcltrgArray[m_hitEneNum]->setBhabha17(vct_bhabha[16]);
  trgEcltrgArray[m_hitEneNum]->setBhabha18(vct_bhabha[17]);
  //
  trgEcltrgArray[m_hitEneNum]->setICN(icn);
  trgEcltrgArray[m_hitEneNum]->setICNFw(icnfwd);
  trgEcltrgArray[m_hitEneNum]->setICNBr(icnbr);
  trgEcltrgArray[m_hitEneNum]->setICNBw(icnbwd);
  //
  trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[0], 0);
  trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[1], 1);
  trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[2], 2);
  trgEcltrgArray[m_hitEneNum]->setECLtoGDL(m_Triggerbit[3], 3);
 
  trgEcltrgArray[m_hitEneNum]->setBhabhaVeto(bhabha2D);
  trgEcltrgArray[m_hitEneNum]->setBeamBkgVeto(beambkgtag);
  trgEcltrgArray[m_hitEneNum]->setEventTiming(eventtiming);
 
  trgEcltrgArray[m_hitEneNum]->setHit(hit);
  trgEcltrgArray[m_hitEneNum]->setRevoclk(RevoFAM);
  trgEcltrgArray[m_hitEneNum]->setTimingSource(TimingSource);
  trgEcltrgArray[m_hitEneNum]->setPhysics(physics) ;
  trgEcltrgArray[m_hitEneNum]->set2DBhabha(bhabha2D);
  trgEcltrgArray[m_hitEneNum]->set3DBhabha(bhabha3D_veto);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSel(bhabha3D_sel);
  trgEcltrgArray[m_hitEneNum]->setmumuBit(mumu);
  trgEcltrgArray[m_hitEneNum]->setBhabhaPrescaleBit(bhabhaprescale);
 
  trgEcltrgArray[m_hitEneNum]->setELow(ELow)  ;
  trgEcltrgArray[m_hitEneNum]->setEHihg(EHigh);
  trgEcltrgArray[m_hitEneNum]->setELum(ELum)  ;
  trgEcltrgArray[m_hitEneNum]->setClusterOverflow(ClusterOverflow) ;
  trgEcltrgArray[m_hitEneNum]->setLowMultiBit(m_Lowmultibit);
 
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoInTrackFlag(m_obj_bhabha->get3DBhabhaVetoInTrackFlag());
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterTCId(m_obj_bhabha->get3DBhabhaVetoClusterTCId(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterTCId(m_obj_bhabha->get3DBhabhaVetoClusterTCId(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterEnergy(m_obj_bhabha->get3DBhabhaVetoClusterEnergy(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterEnergy(m_obj_bhabha->get3DBhabhaVetoClusterEnergy(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterTiming(m_obj_bhabha->get3DBhabhaVetoClusterTiming(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterTiming(m_obj_bhabha->get3DBhabhaVetoClusterTiming(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterThetaId(m_obj_bhabha->get3DBhabhaVetoClusterThetaId(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaVetoClusterThetaId(m_obj_bhabha->get3DBhabhaVetoClusterThetaId(1), 1);
 
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionThetaFlag(m_obj_bhabha->get3DBhabhaSelectionThetaFlag());
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterTCId(m_obj_bhabha->get3DBhabhaSelectionClusterTCId(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterTCId(m_obj_bhabha->get3DBhabhaSelectionClusterTCId(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterEnergy(m_obj_bhabha->get3DBhabhaSelectionClusterEnergy(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterEnergy(m_obj_bhabha->get3DBhabhaSelectionClusterEnergy(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterTiming(m_obj_bhabha->get3DBhabhaSelectionClusterTiming(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterTiming(m_obj_bhabha->get3DBhabhaSelectionClusterTiming(1), 1);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterThetaId(m_obj_bhabha->get3DBhabhaSelectionClusterThetaId(0), 0);
  trgEcltrgArray[m_hitEneNum]->set3DBhabhaSelectionClusterThetaId(m_obj_bhabha->get3DBhabhaSelectionClusterThetaId(1), 1);
 
  trgEcltrgArray[m_hitEneNum]->setEventTimingQualityFlag(m_obj_timing->getEventTimingQualityFlag());
  trgEcltrgArray[m_hitEneNum]->setEventTimingTCId(m_obj_timing->getEventTimingTCId());
  trgEcltrgArray[m_hitEneNum]->setEventTimingTCThetaId(m_obj_timing->getEventTimingTCThetaId());
  trgEcltrgArray[m_hitEneNum]->setEventTimingTCEnergy(m_obj_timing->getEventTimingTCEnergy());
 
  trgEcltrgArray[m_hitEneNum]->setEtot1to17(E_total);
  trgEcltrgArray[m_hitEneNum]->setTaub2bFlag(taub2bFlag);
  trgEcltrgArray[m_hitEneNum]->setTaub2bAngleFlag(m_obj_bhabha->getTaub2bAngleFlag());
  trgEcltrgArray[m_hitEneNum]->setTaub2bEtotFlag(m_obj_bhabha->getTaub2bEtotFlag());
  trgEcltrgArray[m_hitEneNum]->setTaub2bClusterEFlag(m_obj_bhabha->getTaub2bClusterEFlag());
 
  trgEcltrgArray[m_hitEneNum]->setNofCluster1to17(NofCluster1to17);
 
  trgEcltrgArray[m_hitEneNum]->setDataClockWindowStartTime(WindowStart);
 
  return;
}

version()

std::string version

(

void

)

const

returns version.

Definition at line 132 of file TrgEclMaster.cc.

{
  return std::string("TrgEclMaster 2.1");
}

Member Data Documentation

m_2DBhabhaThresholdBWD

std::vector<double> m_2DBhabhaThresholdBWD

private

2D Bhabha Energy Threshold

Definition at line 279 of file TrgEclMaster.h.

m_2DBhabhaThresholdFWD

std::vector<double> m_2DBhabhaThresholdFWD

private

2D Bhabha Energy Threshold

Definition at line 277 of file TrgEclMaster.h.

m_3DBhabhaAddAngleCut

std::vector<double> m_3DBhabhaAddAngleCut

private

Angle selection of additional Bhabha addition in CM frame.

Definition at line 293 of file TrgEclMaster.h.

m_3DBhabhaSelectionAngle

std::vector<double> m_3DBhabhaSelectionAngle

private

3D Selection Bhabha Energy Angle

Definition at line 285 of file TrgEclMaster.h.

m_3DBhabhaSelectionThreshold

std::vector<double> m_3DBhabhaSelectionThreshold

private

3D Selection Bhabha Energy Threshold

Definition at line 281 of file TrgEclMaster.h.

m_3DBhabhaVetoAngle

std::vector<double> m_3DBhabhaVetoAngle

private

3D Veto Bhabha Energy Angle

Definition at line 287 of file TrgEclMaster.h.

m_3DBhabhaVetoInTrackThetaRegion

std::vector<int> m_3DBhabhaVetoInTrackThetaRegion

private

Theta region (low, high) of 3D Bhabha Veto InTrack.

Definition at line 332 of file TrgEclMaster.h.

m_3DBhabhaVetoThreshold

std::vector<double> m_3DBhabhaVetoThreshold

private

3D Veto Bhabha Energy Threshold

Definition at line 283 of file TrgEclMaster.h.

m_Bhabha

int m_Bhabha

private

Bhabha option.

Definition at line 260 of file TrgEclMaster.h.

m_Clustering

int m_Clustering

private

clutering option

Definition at line 258 of file TrgEclMaster.h.

m_ClusterLimit

int m_ClusterLimit

private

The limit number of Cluster.

Definition at line 266 of file TrgEclMaster.h.

m_ECLBurstThreshold

double m_ECLBurstThreshold

private

ECL Burst Bit Threshold.

Definition at line 326 of file TrgEclMaster.h.

m_EventTiming

int m_EventTiming

private

EventTiming option.

Definition at line 262 of file TrgEclMaster.h.

m_EventTimingQualityThresholds

std::vector<double> m_EventTimingQualityThresholds

private

energy threshold(low, high) for quality flag (GeV)

Definition at line 334 of file TrgEclMaster.h.

m_HitTCId

std::vector<int> m_HitTCId

private

Hit TC Energy in time window.

Definition at line 234 of file TrgEclMaster.h.

m_Lowmultibit

int m_Lowmultibit

private

Low Multiplicity bit.

Definition at line 270 of file TrgEclMaster.h.

m_LowMultiThreshold

std::vector<double> m_LowMultiThreshold

private

Low Multiplicity Threshold.

Definition at line 330 of file TrgEclMaster.h.

m_mumuAngle

std::vector<double> m_mumuAngle

private

mumu bit Angle

Definition at line 291 of file TrgEclMaster.h.

m_mumuThreshold

double m_mumuThreshold

private

mumu bit Energy Threshold

Definition at line 289 of file TrgEclMaster.h.

m_n300MeVCluster

int m_n300MeVCluster

private

The number of Cluster exceeding 300 MeV.

Definition at line 324 of file TrgEclMaster.h.

m_NofTopTC

int m_NofTopTC

private

of considered TC in energy weighted Timing method

Definition at line 264 of file TrgEclMaster.h.

m_obj_beambkg

TrgEclBeamBKG* m_obj_beambkg

private

Beam Backgroud veto object.

Definition at line 345 of file TrgEclMaster.h.

m_obj_bhabha

TrgEclBhabha* m_obj_bhabha

private

Bhabha object.

Definition at line 343 of file TrgEclMaster.h.

m_obj_cluster

TrgEclCluster* m_obj_cluster

private

Cluster object.

Definition at line 339 of file TrgEclMaster.h.

m_obj_database

TrgEclDataBase* m_obj_database

private

Beam Backgroud veto object.

Definition at line 347 of file TrgEclMaster.h.

m_obj_map

TrgEclMapping* m_obj_map

private

Mapping object.

Definition at line 337 of file TrgEclMaster.h.

m_obj_timing

TrgEclTiming* m_obj_timing

private

EventTiming object.

Definition at line 341 of file TrgEclMaster.h.

m_OverlapWindow

double m_OverlapWindow

private

TRG Decision overlap window.

Definition at line 255 of file TrgEclMaster.h.

m_PhiRingSum

std::vector< std::vector<double> > m_PhiRingSum

private

Phi ring sum.

Definition at line 243 of file TrgEclMaster.h.

m_PrescaleCounter

int m_PrescaleCounter

private

Bhabha Prescale Countor.

Definition at line 274 of file TrgEclMaster.h.

m_PrescaleFactor

int m_PrescaleFactor

private

Bhabha Prescale Factor.

Definition at line 272 of file TrgEclMaster.h.

m_taub2b2AngleCut

std::vector<int> m_taub2b2AngleCut

private

taub2b2 angle selection(degree) (3,2,1,0) = (dphi low, dphi high, theta_sum low, theta_sum high)

Definition at line 305 of file TrgEclMaster.h.

m_taub2b2CLECut

double m_taub2b2CLECut

private

taub2b2 cluseter energy cut (GeV)

Definition at line 311 of file TrgEclMaster.h.

m_taub2b2CLEEndcapCut

double m_taub2b2CLEEndcapCut

private

taub2b2 cluster energy cut for endcap cluster (GeV)

Definition at line 309 of file TrgEclMaster.h.

m_taub2b2EtotCut

double m_taub2b2EtotCut

private

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

Definition at line 307 of file TrgEclMaster.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 314 of file TrgEclMaster.h.

m_taub2b3CLEb2bCut

double m_taub2b3CLEb2bCut

private

taub2b3 cluster energy cut in lab for one of b2b clusters (GeV)

Definition at line 318 of file TrgEclMaster.h.

m_taub2b3CLEHighCut

double m_taub2b3CLEHighCut

private

taub2b3 cluster energy high cut in lab for all clusters (GeV)

Definition at line 322 of file TrgEclMaster.h.

m_taub2b3CLELowCut

double m_taub2b3CLELowCut

private

taub2b3 cluster energy low cut in lab for all clusters (GeV)

Definition at line 320 of file TrgEclMaster.h.

m_taub2b3EtotCut

double m_taub2b3EtotCut

private

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

Definition at line 316 of file TrgEclMaster.h.

m_taub2bAngleCut

std::vector<int> m_taub2bAngleCut

private

tau b2b 2 cluster angle cut (degree) (dphi low, dphi high, theta_sum low, theta_sum high)

Definition at line 296 of file TrgEclMaster.h.

m_taub2bClusterECut1

double m_taub2bClusterECut1

private

taub2b one Cluster energy selection (GeV)

Definition at line 300 of file TrgEclMaster.h.

m_taub2bClusterECut2

double m_taub2bClusterECut2

private

taub2b one Cluster energy selection (GeV)

Definition at line 302 of file TrgEclMaster.h.

m_taub2bEtotCut

double m_taub2bEtotCut

private

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

Definition at line 298 of file TrgEclMaster.h.

m_TCBeamBkgTag

std::vector<std::vector<int> > m_TCBeamBkgTag

private

Hit TC Beam Background tag.

Definition at line 231 of file TrgEclMaster.h.

m_TCEnergy

std::vector<std::vector<double> > m_TCEnergy

private

Hit TC Energy.

Definition at line 227 of file TrgEclMaster.h.

m_TCHitBeamBkgTag

std::vector<int> m_TCHitBeamBkgTag

private

Hit TC Beam Background tag in time window.

Definition at line 240 of file TrgEclMaster.h.

m_TCHitEnergy

std::vector<double> m_TCHitEnergy

private

Hit TC Energy in time window.

Definition at line 236 of file TrgEclMaster.h.

m_TCHitTiming

std::vector<double> m_TCHitTiming

private

Hit TC Timing in time window.

Definition at line 238 of file TrgEclMaster.h.

m_TCTiming

std::vector<std::vector<double> > m_TCTiming

private

Hit TC Timing.

Definition at line 229 of file TrgEclMaster.h.

m_ThetaRingSum

std::vector<std::vector<std::vector<double> > > m_ThetaRingSum

private

Theta ring sum.

Definition at line 245 of file TrgEclMaster.h.

m_TimeWindow

double m_TimeWindow

private

Hit TC Energy in time window.

Hit TC Timing in time window TRG Decision Time window

Definition at line 253 of file TrgEclMaster.h.

m_TotalEnergy

std::vector<double> m_TotalEnergy

private

Total Energy Theshold (low, high, lum)

Definition at line 328 of file TrgEclMaster.h.

m_Triggerbit

int m_Triggerbit[4]

private

ECL Trigger bit.

Definition at line 268 of file TrgEclMaster.h.

---

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

• trg/ecl/include/TrgEclMaster.h
• trg/ecl/src/TrgEclMaster.cc