BesTofDigitizerEcV4 Class Reference

#include <BesTofDigitizerEcV4.hh>

Inheritance diagram for BesTofDigitizerEcV4:

Classes
struct	Param
struct	HitStruct
struct	StripStruct

Public Member Functions
	BesTofDigitizerEcV4 ()
	~BesTofDigitizerEcV4 ()
void	initial ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)
int	calStrip (double locZ)
bool	underStrip (double locX, double locZ)
double	calTdcRes_charge (double charge_fC)
double	charge2Time (double charge_fC)
double	calAdcRes_charge (double charge_fC)
double	calTdcRes_charge1 (double charge_fC)
double	charge2Time1 (double charge_fC)
double	calAdcRes_charge1 (double charge_fC)
	BesTofDigitizerEcV4 ()
	~BesTofDigitizerEcV4 ()
void	initial ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)
int	calStrip (double locZ)
bool	underStrip (double locX, double locZ)
double	calTdcRes_charge (double charge_fC)
double	charge2Time (double charge_fC)
double	calAdcRes_charge (double charge_fC)
double	calTdcRes_charge1 (double charge_fC)
double	charge2Time1 (double charge_fC)
double	calAdcRes_charge1 (double charge_fC)
	BesTofDigitizerEcV4 ()
	~BesTofDigitizerEcV4 ()
void	initial ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)
int	calStrip (double locZ)
bool	underStrip (double locX, double locZ)
double	calTdcRes_charge (double charge_fC)
double	charge2Time (double charge_fC)
double	calAdcRes_charge (double charge_fC)
double	calTdcRes_charge1 (double charge_fC)
double	charge2Time1 (double charge_fC)
double	calAdcRes_charge1 (double charge_fC)
Public Member Functions inherited from BesTofDigitizerV
	BesTofDigitizerV ()
	~BesTofDigitizerV ()
void	Initialize ()
	BesTofDigitizerV ()
	~BesTofDigitizerV ()
void	Initialize ()
	BesTofDigitizerV ()
	~BesTofDigitizerV ()
void	Initialize ()

Additional Inherited Members
Protected Attributes inherited from BesTofDigitizerV
BesTofDigitsCollection *	m_besTofDigitsCollection
BesTofHitsCollection *	m_THC
ITofCaliSvc *	m_tofCaliSvc
ITofSimSvc *	m_tofSimSvc
ITofQElecSvc *	m_tofQElecSvc
G4double	m_ADC [2]
G4double	m_TDC [2]
G4int	m_trackIndex
G4double	m_globalTime
Static Protected Attributes inherited from BesTofDigitizerV
static bool	m_booked = false
static NTuple::Tuple *	m_tupleTof1 = 0
static NTuple::Item< double >	m_partId
static NTuple::Item< double >	m_scinNb
static NTuple::Item< double >	m_edep
static NTuple::Item< double >	m_nHits
static NTuple::Item< double >	m_time1st0
static NTuple::Item< double >	m_time1st1
static NTuple::Item< double >	m_timelast0
static NTuple::Item< double >	m_timelast1
static NTuple::Item< double >	m_totalPhot0
static NTuple::Item< double >	m_totalPhot1
static NTuple::Item< double >	m_NphAllSteps
static NTuple::Item< double >	m_max0
static NTuple::Item< double >	m_max1
static NTuple::Item< double >	m_tdc0
static NTuple::Item< double >	m_adc0
static NTuple::Item< double >	m_tdc1
static NTuple::Item< double >	m_adc1
static NTuple::Tuple *	m_tupleTof2 = 0
static NTuple::Item< double >	m_eTotal
static NTuple::Item< double >	m_nDigi
static NTuple::Item< double >	m_partIdMPV
static NTuple::Item< double >	m_scinNbMPV
static NTuple::Item< double >	m_edepMPV
static NTuple::Item< double >	m_nDigiOut
static NTuple::Tuple *	m_tupleTof3 = 0
static NTuple::Item< double >	m_forb
static NTuple::Item< double >	m_timeFlight
static NTuple::Item< double >	m_ddT
static NTuple::Item< double >	m_scinSwim
static NTuple::Item< double >	m_scinTime
static NTuple::Item< double >	m_transitTime
static NTuple::Item< double >	m_endTime
static NTuple::Item< double >	m_edepHit

Detailed Description

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-dbg/include/TofSim/BesTofDigitizerEcV4.hh.

Constructor & Destructor Documentation

BesTofDigitizerEcV4() [1/3]

BesTofDigitizerEcV4::BesTofDigitizerEcV4

(

)

Definition at line 36 of file BesTofDigitizerEcV4.cc.

    : tdcRes_const( "tdcRes_const", -1 )
    , adcRes_const( "adcRes_const", -1 )
    , m_rootFlag( "RootFlag", false )
    , m_fileName( "FileName", "mrpc.root" )
    , m_V( "E", 7000 )
    , m_threshold( "Threshold", 5.5e+08 )
    , m_nstep( "nstep", -1 )
    , m_E_weight( "E_weight", -1 )
    , m_saturationFlag( "saturationFlag", true )
    , m_calTdcRes_charge_flag( "calTdcRes_charge_flag", 0 )
    , m_charge2Time_flag( "charge2Time_flag", 0 )
    , m_calAdcRes_charge_flag( "calAdcRes_charge_flag", 0 ) {
 
  // declare properties
  auto opt_svc = Gaudi::svcLocator()->service<IOptionsSvc>( "JobOptionsSvc" );
  opt_svc->bind( "BesTofDigitizerEcV4", &tdcRes_const );
  opt_svc->bind( "BesTofDigitizerEcV4", &adcRes_const );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_rootFlag );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_fileName );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_V );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_threshold );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_nstep );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_E_weight );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_saturationFlag );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_calTdcRes_charge_flag );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_charge2Time_flag );
  opt_svc->bind( "BesTofDigitizerEcV4", &m_calAdcRes_charge_flag );
 
  initial();
 
  if ( m_rootFlag )
  {
    m_file = new TFile( m_fileName.value().c_str(), "RECREATE" );
    m_tree = new TTree( "mrpc", "mrpc" );
 
    m_tree->Branch( "event", &m_event, "event/D" );
    m_tree->Branch( "partId", &m_partId, "partId/D" );
    m_tree->Branch( "module", &m_module, "module/D" );
    m_tree->Branch( "time_leading_sphi", &m_time_leading_sphi, "time_leading_sphi/D" );
    m_tree->Branch( "time_leading_xphi", &m_time_leading_xphi, "time_leading_xphi/D" );
    m_tree->Branch( "time_trailing_sphi", &m_time_trailing_sphi, "time_trailing_sphi/D" );
    m_tree->Branch( "time_trailing_xphi", &m_time_trailing_xphi, "time_trailing_xphi/D" );
    m_tree->Branch( "tdcRes", &m_tdcRes, "tdcRes/D" );
    m_tree->Branch( "tdcRes_charge", &m_tdcRes_charge, "tdcRes_charge/D" );
    m_tree->Branch( "adc", &m_adc, "adc/D" );
    m_tree->Branch( "adcRes", &m_adcRes, "adcRes/D" );
    m_tree->Branch( "adcRes_charge", &m_adcRes_charge, "adcRes_charge/D" );
    m_tree->Branch( "strip", &m_strip, "strip/D" );
    m_tree->Branch( "trkIndex", &m_trkIndex, "trkIndex/D" );
    m_tree->Branch( "tStart", &m_tStart, "tStart/D" );
    m_tree->Branch( "tPropagate_sphi", &m_tPropagate_sphi, "tPropagate_sphi/D" );
    m_tree->Branch( "tPropagate_xphi", &m_tPropagate_xphi, "tPropagate_xphi/D" );
    m_tree->Branch( "tThreshold", &m_tThreshold, "tThreshold/D" );
    m_tree->Branch( "charge", &m_charge, "charge/D" );
    m_tree->Branch( "nhit", &m_nhit, "nhit/I" );
    m_tree->Branch( "ions_hit", m_ions_hit, "ions_hit[nhit]/D" );
    m_tree->Branch( "trkIndex_hit", m_trkIndex_hit, "trkIndex_hit[nhit]/D" );
    m_tree->Branch( "pdgCode_hit", m_pdgCode_hit, "pdgCode_hit[nhit]/D" );
    m_tree->Branch( "gap_hit", m_gap_hit, "gap_hit[nhit]/D" );
    m_tree->Branch( "underStrip_hit", m_underStrip_hit, "underStrip_hit[nhit]/D" );
    m_tree->Branch( "locx_hit", m_locx_hit, "locx_hit[nhit]/D" );
    m_tree->Branch( "locy_hit", m_locy_hit, "locy_hit[nhit]/D" );
    m_tree->Branch( "locz_hit", m_locz_hit, "locz_hit[nhit]/D" );
    m_tree->Branch( "x_hit", m_x_hit, "x_hit[nhit]/D" );
    m_tree->Branch( "y_hit", m_y_hit, "y_hit[nhit]/D" );
    m_tree->Branch( "z_hit", m_z_hit, "z_hit[nhit]/D" );
    m_tree->Branch( "px_hit", m_px_hit, "px_hit[nhit]/D" );
    m_tree->Branch( "py_hit", m_py_hit, "py_hit[nhit]/D" );
    m_tree->Branch( "pz_hit", m_pz_hit, "pz_hit[nhit]/D" );
  }
}

~BesTofDigitizerEcV4() [1/3]

BesTofDigitizerEcV4::~BesTofDigitizerEcV4

(

)

Definition at line 109 of file BesTofDigitizerEcV4.cc.

                                          {
  if ( m_rootFlag )
  {
    m_file->Write();
    m_file->Close();
  }
}

BesTofDigitizerEcV4() [2/3]

BesTofDigitizerEcV4::BesTofDigitizerEcV4

(

)

~BesTofDigitizerEcV4() [2/3]

BesTofDigitizerEcV4::~BesTofDigitizerEcV4

(

)

BesTofDigitizerEcV4() [3/3]

BesTofDigitizerEcV4::BesTofDigitizerEcV4

(

)

~BesTofDigitizerEcV4() [3/3]

BesTofDigitizerEcV4::~BesTofDigitizerEcV4

(

)

Member Function Documentation

calAdcRes_charge() [1/3]

double BesTofDigitizerEcV4::calAdcRes_charge

(

double

charge_fC

)

Definition at line 547 of file BesTofDigitizerEcV4.cc.

                                                               {
  double time = 0;
  if ( charge_fC < 250 )
  { time = 72.6005 * exp( -charge_fC * 0.0302626 + 1.49059 ) + 40.8627; }
  else { time = 32.6233 + 0.00404149 * charge_fC; }
  if ( time < 0 ) time = 0;
  return time; // ps
}

Referenced by Digitize().

calAdcRes_charge() [2/3]

double BesTofDigitizerEcV4::calAdcRes_charge

(

double

charge_fC

)

calAdcRes_charge() [3/3]

double BesTofDigitizerEcV4::calAdcRes_charge

(

double

charge_fC

)

calAdcRes_charge1() [1/3]

double BesTofDigitizerEcV4::calAdcRes_charge1

(

double

charge_fC

)

Definition at line 572 of file BesTofDigitizerEcV4.cc.

                                                                {
  double time = 64.3326 * exp( -charge_fC / 25.4638 + 0.944184 ) + 19.4846;
  if ( time < 0 ) time = 0;
  return time; // ps
}

Referenced by Digitize().

calAdcRes_charge1() [2/3]

double BesTofDigitizerEcV4::calAdcRes_charge1

(

double

charge_fC

)

calAdcRes_charge1() [3/3]

double BesTofDigitizerEcV4::calAdcRes_charge1

(

double

charge_fC

)

calStrip() [1/3]

int BesTofDigitizerEcV4::calStrip

(

double

locZ

)

Definition at line 343 of file BesTofDigitizerEcV4.cc.

                                               {
  int    strip      = -1;
  double stripWidth = m_param.strip_z + m_param.strip_gap; // Strip spread: (24+3)mm
  int    nstrip     = m_param.nstrip;
  // the offset between strip coordinate and gas SD: 0.5
  double length = locZ + stripWidth * nstrip / 2 - 0.5;
  if ( length <= 0 ) { strip = 0; }
  else if ( length < stripWidth * nstrip )
  {
    for ( int i = 0; i < nstrip; i++ )
    {
      if ( length > i * stripWidth && length < ( i + 1 ) * stripWidth )
      {
        strip = i;
        break;
      }
    }
  }
  else { strip = nstrip - 1; }
  if ( strip < 0 ) strip = 0;
  if ( strip > nstrip - 1 ) strip = nstrip - 1;
 
  return strip;
}

Referenced by Digitize().

calStrip() [2/3]

int BesTofDigitizerEcV4::calStrip

(

double

locZ

)

calStrip() [3/3]

int BesTofDigitizerEcV4::calStrip

(

double

locZ

)

calTdcRes_charge() [1/3]

double BesTofDigitizerEcV4::calTdcRes_charge

(

double

charge_fC

)

Definition at line 531 of file BesTofDigitizerEcV4.cc.

                                                               {
  double time = 0;
  if ( charge_fC < 250 )
  { time = 100.764 * exp( -charge_fC * 0.0413966 + 0.377154 ) + 13.814; }
  else { time = 12.8562 + 0.000507299 * charge_fC; }
  if ( time < 0 ) time = 0;
  return time; // ps
}

Referenced by Digitize().

calTdcRes_charge() [2/3]

double BesTofDigitizerEcV4::calTdcRes_charge

(

double

charge_fC

)

calTdcRes_charge() [3/3]

double BesTofDigitizerEcV4::calTdcRes_charge

(

double

charge_fC

)

calTdcRes_charge1() [1/3]

double BesTofDigitizerEcV4::calTdcRes_charge1

(

double

charge_fC

)

Definition at line 556 of file BesTofDigitizerEcV4.cc.

                                                                {
  double result = 0;
  if ( charge_fC > 50. )
  { result = 67.6737 * exp( -charge_fC / 50.9995 - 0.27755 ) + 9.06223; }
  else { result = 176.322 - 2.98345 * charge_fC; }
  if ( result < 0 ) result = 0;
  return result; // ps
}

Referenced by Digitize().

calTdcRes_charge1() [2/3]

double BesTofDigitizerEcV4::calTdcRes_charge1

(

double

charge_fC

)

calTdcRes_charge1() [3/3]

double BesTofDigitizerEcV4::calTdcRes_charge1

(

double

charge_fC

)

charge2Time() [1/3]

double BesTofDigitizerEcV4::charge2Time

(

double

charge_fC

)

Definition at line 540 of file BesTofDigitizerEcV4.cc.

                                                          {
  double time = 0;
  time        = -120.808 / log( charge_fC * 30.1306 ) + 26.6024; // ns
  if ( time < 0 ) time = 0;
  return time;
}

Referenced by Digitize().

charge2Time() [2/3]

double BesTofDigitizerEcV4::charge2Time

(

double

charge_fC

)

charge2Time() [3/3]

double BesTofDigitizerEcV4::charge2Time

(

double

charge_fC

)

charge2Time1() [1/3]

double BesTofDigitizerEcV4::charge2Time1

(

double

charge_fC

)

Definition at line 565 of file BesTofDigitizerEcV4.cc.

                                                           {
  double time = 0;
  time        = -4.50565 / log( charge_fC * 0.0812208 ) + 16.6653; // ns
  if ( time < 0 ) time = 0;
  return time;
}

Referenced by Digitize().

charge2Time1() [2/3]

double BesTofDigitizerEcV4::charge2Time1

(

double

charge_fC

)

charge2Time1() [3/3]

double BesTofDigitizerEcV4::charge2Time1

(

double

charge_fC

)

Digitize() [1/3]

void BesTofDigitizerEcV4::Digitize ( ScintSingle * single_module, BesTofDigitsCollection * DC )

virtual

Reimplemented from BesTofDigitizerV.

Definition at line 149 of file BesTofDigitizerEcV4.cc.

                                                                                           {
  m_besTofDigitsCollection   = DC;
  G4DigiManager* digiManager = G4DigiManager::GetDMpointer();
  G4int          THCID       = digiManager->GetHitsCollectionID( "BesTofHitsCollection" );
  m_THC = (BesTofHitsCollection*)( digiManager->GetHitsCollection( THCID ) );
  if ( !m_THC ) return;
 
  partId = single_module->GetPartId();
  module = single_module->GetModule_mrpc();
  // cout<<"partId= "<<partId<<"  module= "<<module<<endl;
 
  // Process the hits
  int         nstrip = m_param.nstrip;
  StripStruct stripStruct[12];
  for ( int i = 0; i < nstrip; i++ )
  {
    stripStruct[i].m_param = m_param;
    stripStruct[i].setPar( m_param.alpha, m_param.eta, m_param.v_drift, m_threshold,
                           m_saturationFlag );
  }
 
  BesTofHit* hit;
  for ( unsigned int i = 0; i < single_module->GetHitIndexes_mrpc()->size(); i++ )
  {
    hit     = ( *m_THC )[( *( single_module->GetHitIndexes_mrpc() ) )[i]];
    m_event = hit->GetEvent();
 
    HitStruct hitStruct;
    hitStruct.m_param    = m_param;
    hitStruct.trkIndex   = hit->GetG4Index();
    hitStruct.pdgCode    = hit->GetPDGcode();
    hitStruct.ions       = hit->GetIons();
    hitStruct.strip      = calStrip( hit->GetLocPos().z() / mm );
    hitStruct.underStrip = underStrip( hit->GetLocPos().x() / mm, hit->GetLocPos().z() / mm );
    hitStruct.gap        = hit->GetGapNb();
    hitStruct.glbTime    = hit->GetTime() / ns;
    hitStruct.locx       = hit->GetLocPos().x() / mm;
    hitStruct.locy       = hit->GetLocPos().y() / mm;
    hitStruct.locz       = hit->GetLocPos().z() / mm;
    hitStruct.x          = hit->GetPos().x() / mm;
    hitStruct.y          = hit->GetPos().y() / mm;
    hitStruct.z          = hit->GetPos().z() / mm;
    hitStruct.px         = hit->GetMomentum().x() / ( GeV / ( 3e+08 * m / s ) );
    hitStruct.py         = hit->GetMomentum().y() / ( GeV / ( 3e+08 * m / s ) );
    hitStruct.pz         = hit->GetMomentum().z() / ( GeV / ( 3e+08 * m / s ) );
    // hitStruct.print();
 
    if ( hitStruct.ions > 0 && hitStruct.glbTime > 0 )
    { stripStruct[hitStruct.strip].hitStructCol.push_back( hitStruct ); }
  }
 
  // test multihit, study the lower peak in charge
  int hitSize = 0;
 
  for ( int i = 0; i < nstrip; i++ )
  {
    if ( stripStruct[i].hitStructCol.size() == 0 ) continue;
    stripStruct[i].strip = i;
    stripStruct[i].calFirstHit();
    stripStruct[i].avalanche();
    // stripStruct[i].print();
 
    if ( stripStruct[i].tThreshold <= 0 ) continue;
 
    tdc_sphi =
        stripStruct[i].tStart + stripStruct[i].tThreshold + stripStruct[i].tPropagate_sphi;
    tdc_xphi =
        stripStruct[i].tStart + stripStruct[i].tThreshold + stripStruct[i].tPropagate_xphi;
 
    double tdcRes_charge = 0;
    if ( m_calTdcRes_charge_flag == 0 )
    {
      tdcRes_charge = calTdcRes_charge( stripStruct[i].charge * 1000 ); // ps, charge in fC
    }
    else if ( m_calTdcRes_charge_flag == 1 )
    { tdcRes_charge = calTdcRes_charge1( stripStruct[i].charge * 1000 ); }
    else if ( m_calTdcRes_charge_flag == 2 ) { tdcRes_charge = 0; }
 
    tdcRes = sqrt( tdcRes_charge * tdcRes_charge + tdcRes_const * tdcRes_const ) / 1000; // ns
 
    tdc_sphi = G4RandGauss::shoot( tdc_sphi, tdcRes );
    tdc_xphi = G4RandGauss::shoot( tdc_xphi, tdcRes );
 
    if ( m_charge2Time_flag == 0 )
    {
      adc = charge2Time( stripStruct[i].charge * 1000 ); // ns, charge in fC
    }
    else if ( m_charge2Time_flag == 1 ) { adc = charge2Time1( stripStruct[i].charge * 1000 ); }
 
    double adcRes_charge = 0;
    if ( m_calAdcRes_charge_flag == 0 )
    {
      adcRes_charge = calAdcRes_charge( stripStruct[i].charge * 1000 ); // ps, charge in fC
    }
    else if ( m_calAdcRes_charge_flag == 1 )
    { adcRes_charge = calAdcRes_charge1( stripStruct[i].charge * 1000 ); }
    else if ( m_calAdcRes_charge_flag == 2 ) { adcRes_charge = 0; }
 
    adcRes = sqrt( adcRes_charge * adcRes_charge + adcRes_const * adcRes_const ) / 1000;
    adc    = G4RandGauss::shoot( adc, adcRes );
    if ( adc < 0 ) adc = 0;
 
    time_leading_sphi  = tdc_sphi;
    time_leading_xphi  = tdc_xphi;
    time_trailing_sphi = tdc_sphi + adc;
    time_trailing_xphi = tdc_xphi + adc;
 
    // save digi information
    BesTofDigi* digi = new BesTofDigi;
    digi->SetTrackIndex( stripStruct[i].trkIndex );
    digi->SetPartId( partId );
    digi->SetModule( module );
    digi->SetStrip( stripStruct[i].strip );
    int mo = ( partId - 3 ) * 36 + module;
    int st = stripStruct[i].strip;
    if ( m_param.deadChannel[mo][st] == 0 || m_param.deadChannel[mo][st] == 2 )
    {
      // Set dead channel
      digi->SetForwT1( -999 );
      digi->SetForwT2( -999 );
    }
    else
    {
      digi->SetForwT1( time_leading_sphi );
      digi->SetForwT2( time_trailing_sphi );
    }
    if ( m_param.deadChannel[mo][st] == 1 || m_param.deadChannel[mo][st] == 2 )
    {
      digi->SetBackT1( -999 );
      digi->SetBackT2( -999 );
    }
    else
    {
      digi->SetBackT1( time_leading_xphi );
      digi->SetBackT2( time_trailing_xphi );
    }
    m_besTofDigitsCollection->insert( digi );
    // cout<<"Print digi info: "
    //     <<"  partId= "<<partId
    //     <<"  module= "<<module
    //     <<"  strip= "<<stripStruct[i].strip
    //     <<"  time_leading_sphi= "<<time_leading_sphi
    //     <<"  time_leading_xphi= "<<time_leading_xphi
    //     <<"  time_trailing_sphi= "<<time_trailing_sphi
    //     <<"  time_trailing_xphi= "<<time_trailing_xphi
    //     <<endl;
 
    // save digi information
    if ( m_rootFlag )
    {
      m_partId             = partId;
      m_module             = module;
      m_time_leading_sphi  = time_leading_sphi;
      m_time_leading_xphi  = time_leading_xphi;
      m_time_trailing_sphi = time_trailing_sphi;
      m_time_trailing_xphi = time_trailing_xphi;
      m_tdcRes             = tdcRes;
      m_tdcRes_charge      = tdcRes_charge;
      m_adc                = adc;
      m_adcRes             = adcRes;
      m_adcRes_charge      = adcRes_charge;
 
      m_strip           = stripStruct[i].strip;
      m_trkIndex        = stripStruct[i].trkIndex;
      m_tStart          = stripStruct[i].tStart;
      m_tPropagate_sphi = stripStruct[i].tPropagate_sphi;
      m_tPropagate_xphi = stripStruct[i].tPropagate_xphi;
      m_tThreshold      = stripStruct[i].tThreshold;
      m_charge          = stripStruct[i].charge;
 
      m_nhit = stripStruct[i].hitStructCol.size();
      // cout<<"m_nhit= "<<m_nhit<<endl;
      for ( int j = 0; j < m_nhit; j++ )
      {
        m_ions_hit[j]       = stripStruct[i].hitStructCol[j].ions;
        m_trkIndex_hit[j]   = stripStruct[i].hitStructCol[j].trkIndex;
        m_pdgCode_hit[j]    = stripStruct[i].hitStructCol[j].pdgCode;
        m_gap_hit[j]        = stripStruct[i].hitStructCol[j].gap;
        m_underStrip_hit[j] = stripStruct[i].hitStructCol[j].underStrip;
        m_locx_hit[j]       = stripStruct[i].hitStructCol[j].locx;
        m_locy_hit[j]       = stripStruct[i].hitStructCol[j].locy;
        m_locz_hit[j]       = stripStruct[i].hitStructCol[j].locz;
        m_x_hit[j]          = stripStruct[i].hitStructCol[j].x;
        m_y_hit[j]          = stripStruct[i].hitStructCol[j].y;
        m_z_hit[j]          = stripStruct[i].hitStructCol[j].z;
        m_px_hit[j]         = stripStruct[i].hitStructCol[j].px;
        m_py_hit[j]         = stripStruct[i].hitStructCol[j].py;
        m_pz_hit[j]         = stripStruct[i].hitStructCol[j].pz;
      }
      m_tree->Fill();
    }
  }
}

Digitize() [2/3]

virtual void BesTofDigitizerEcV4::Digitize ( ScintSingle * , BesTofDigitsCollection *  )

virtual

Reimplemented from BesTofDigitizerV.

Digitize() [3/3]

virtual void BesTofDigitizerEcV4::Digitize ( ScintSingle * , BesTofDigitsCollection *  )

virtual

Reimplemented from BesTofDigitizerV.

initial() [1/3]

void BesTofDigitizerEcV4::initial

(

)

Definition at line 117 of file BesTofDigitizerEcV4.cc.

                                  {
  m_param = Param();
  m_param.setPar( m_nstep, m_E_weight, m_V );
  m_param.print();
  m_event  = -1;
  partId   = -999;
  module   = -999;
  tdc_sphi = -999;
  tdc_xphi = -999;
  if ( tdcRes_const < 0 ) tdcRes_const = 38;
  // 45; //sqrt(27*27.+30.*30+20.*20); //ps, 27:TDC; 30:gapNb; 20:cables..
  tdcRes = -999;
 
  adc = -999;
  if ( adcRes_const < 0 ) adcRes_const = 27; // ps TDC
  adcRes = -999;
 
  time_leading_sphi  = -999;
  time_leading_xphi  = -999;
  time_trailing_sphi = -999;
  time_trailing_xphi = -999;
 
  cout << "Property:" << endl
       << "  FileName= " << m_fileName << "  E= " << m_V << "  Threshold= " << m_threshold
       << "  nstep= " << m_nstep << "  E_weight= " << m_E_weight
       << "  saturationFlag= " << m_saturationFlag << "  tdcRes_const= " << tdcRes_const
       << "  adcRes_const= " << adcRes_const
       << "  calTdcRes_charge_flag= " << m_calTdcRes_charge_flag
       << "  charge2Time_flag= " << m_charge2Time_flag
       << "  calAdcRes_charge_flag= " << m_calAdcRes_charge_flag << endl;
}

Referenced by BesTofDigitizerEcV4(), BesTofDigitizerEcV4::HitStruct::HitStruct(), and BesTofDigitizerEcV4::StripStruct::StripStruct().

initial() [2/3]

void BesTofDigitizerEcV4::initial

(

)

initial() [3/3]

void BesTofDigitizerEcV4::initial

(

)

underStrip() [1/3]

bool BesTofDigitizerEcV4::underStrip	(	double	locX,
		double	locZ )

Definition at line 368 of file BesTofDigitizerEcV4.cc.

                                                               {
  bool   flag       = 0;
  int    strip      = -1;
  double stripWidth = m_param.strip_z + m_param.strip_gap; // Strip spread: (24+3)mm
  int    nstrip     = m_param.nstrip;
  // the offset between strip coordinate and gas SD: 0.5
  double length = locZ + stripWidth * nstrip / 2 - 0.5;
  if ( length < stripWidth * nstrip )
  {
    for ( int i = 0; i < nstrip; i++ )
    {
      if ( length > i * stripWidth && length < ( i + 1 ) * stripWidth )
      {
        strip = i;
        if ( length > i * stripWidth + m_param.strip_gap / 2 &&
             length < ( i + 1 ) * stripWidth - m_param.strip_gap / 2 &&
             locX > -m_param.strip_x[strip] / 2 && locX < m_param.strip_x[strip] / 2 )
          flag = 1;
        break;
      }
    }
  }
 
  return flag;
}

Referenced by Digitize().

underStrip() [2/3]

bool BesTofDigitizerEcV4::underStrip	(	double	locX,
		double	locZ )

underStrip() [3/3]

bool BesTofDigitizerEcV4::underStrip	(	double	locX,
		double	locZ )

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TofSim/BesTofDigitizerEcV4.hh
• InstallArea/x86_64-el9-gcc13-opt/include/TofSim/BesTofDigitizerEcV4.hh
• Simulation/BOOST/TofSim/include/TofSim/BesTofDigitizerEcV4.hh
• BesTofDigitizerEcV4.cc