BesTofDigitizerEcV1 Class Reference

#include <BesTofDigitizerEcV1.hh>

Inheritance diagram for BesTofDigitizerEcV1:

Public Member Functions
	BesTofDigitizerEcV1 ()
	~BesTofDigitizerEcV1 ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)
void	TofPmtInit ()
void	TofPmtAccum (BesTofHit *)
void	Smear (G4int)
Public Member Functions inherited from BesTofDigitizerV
	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 21 of file BesTofDigitizerEcV1.hh.

Constructor & Destructor Documentation

BesTofDigitizerEcV1()

BesTofDigitizerEcV1::BesTofDigitizerEcV1

(

)

Definition at line 21 of file BesTofDigitizerEcV1.cc.

{
    BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
    m_bucketPosR = tofPara->GetBucketPosR();  // 445 ???
}

~BesTofDigitizerEcV1()

BesTofDigitizerEcV1::~BesTofDigitizerEcV1

(

)

Definition at line 27 of file BesTofDigitizerEcV1.cc.

{
}

Member Function Documentation

Digitize()

void BesTofDigitizerEcV1::Digitize ( ScintSingle * scint, BesTofDigitsCollection * DC )

virtual

Reimplemented from BesTofDigitizerV.

Definition at line 31 of file BesTofDigitizerEcV1.cc.

{
    G4cout<<"BesTofDigitizerEcV1::Digitize"<<G4endl;
    m_besTofDigitsCollection = DC;
    G4DigiManager* digiManager = G4DigiManager::GetDMpointer();
    G4int THCID = digiManager->GetHitsCollectionID("BesTofHitsCollection");
    m_THC = (BesTofHitsCollection*) (digiManager->GetHitsCollection(THCID));
    if (m_THC)
    {
        G4int partId, scinNb, nHits;
        BesTofHit* hit;
        partId=scint->GetPartId();
        scinNb=scint->GetScinNb();
        nHits=scint->GetHitIndexes()->size();
        TofPmtInit();
        for (G4int j=0;j<nHits;j++)
        {
            hit= (*m_THC)[( *(scint->GetHitIndexes()) )[j]];
            TofPmtAccum(hit);
        }
 
        Smear(scinNb);
        if ( m_TDC[0]>0 )
        {
            BesTofDigi* digi = new BesTofDigi;
            digi->SetTrackIndex(m_trackIndex);
            digi->SetPartId(partId);
            digi->SetScinNb(scinNb);
            digi->SetForwADC( m_ADC[0]) ;
            digi->SetForwTDC( m_TDC[0]) ;
            digi->SetBackADC( m_ADC[1]) ;
            digi->SetBackTDC( m_TDC[1]) ;
            m_besTofDigitsCollection->insert(digi);
        }
    }
}

Referenced by BesTofDigitizer::Digitize().

Smear()

void BesTofDigitizerEcV1::Smear

(

G4int

scinNb

)

Definition at line 101 of file BesTofDigitizerEcV1.cc.

{
    /*G4double tofRes = 0.08;
    for(G4int i=0;i<2;i++)
    {
      m_TDC[i] += tofRes * G4RandGauss::shoot();
    }*/
 
    double pp[8];
    for (int i=0;i<8;i++)
    {
        pp[i]=m_tofCaliSvc->ETof(scinNb)->getP(i);
    }
    m_ADC[0] *= 7.;
    m_TDC[0] += (pp[0]+pp[1]*m_r)/TMath::Sqrt(m_ADC[0])+
                pp[2]/m_ADC[0]+
                pp[3]*m_r/m_ADC[0]+
                pp[4]*m_r+
                pp[5]*m_r*m_r+
                pp[6]*m_r*m_r*m_r+
                pp[7];
}

Referenced by Digitize().

TofPmtAccum()

void BesTofDigitizerEcV1::TofPmtAccum

(

BesTofHit *

hit

)

Definition at line 74 of file BesTofDigitizerEcV1.cc.

{
    G4int trackIndex = hit->GetTrackIndex();
    G4int scinNb = hit->GetScinNb();
    G4double time = hit->GetTime();
    if (time<m_globalTime)
    {
        m_globalTime = time;
        m_trackIndex = trackIndex;
    }
    G4double edep = hit->GetEdep();
    G4ThreeVector pos = hit->GetPos();
    G4double posx = pos.x();
    G4double posy = pos.y();
    G4double pathL=abs(m_bucketPosR-sqrt(posx*posx+posy*posy));
    G4double atten;
    atten = m_tofCaliSvc->EAtten(scinNb);
    m_ADC[0] += edep*exp(-pathL/atten);
 
    if (time<m_t1st)
    {
        m_t1st = time;
        m_r = sqrt(posx*posx+posy*posy);
        m_TDC[0] = m_t1st;
    }
}

Referenced by Digitize().

TofPmtInit()

void BesTofDigitizerEcV1::TofPmtInit

(

)

Definition at line 68 of file BesTofDigitizerEcV1.cc.

{
    Initialize();
    m_t1st = 9999.;
}

Referenced by Digitize().

---

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

• BesTofDigitizerEcV1.hh
• BesTofDigitizerEcV1.cc