BesTofDigitizerBrV1 Class Reference

#include <BesTofDigitizerBrV1.hh>

Inheritance diagram for BesTofDigitizerBrV1:

Public Member Functions
	BesTofDigitizerBrV1 ()
	~BesTofDigitizerBrV1 ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)
void	TofPmtInit ()
void	TofPmtAccum (BesTofHit *)
void	Smear (G4int)
Public Member Functions inherited from BesTofDigitizerV
	BesTofDigitizerV ()
	~BesTofDigitizerV ()
void	Initialize ()
virtual void	Digitize (ScintSingle *, BesTofDigitsCollection *)

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 BesTofDigitizerBrV1.hh.

Constructor & Destructor Documentation

BesTofDigitizerBrV1()

BesTofDigitizerBrV1::BesTofDigitizerBrV1

(

)

Definition at line 21 of file BesTofDigitizerBrV1.cc.

{
    BesTofGeoParameter* tofPara = BesTofGeoParameter::GetInstance();
    m_scinLength = tofPara->GetBr1L();
}

~BesTofDigitizerBrV1()

BesTofDigitizerBrV1::~BesTofDigitizerBrV1

(

)

Definition at line 27 of file BesTofDigitizerBrV1.cc.

{
}

Member Function Documentation

Digitize()

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

virtual

Reimplemented from BesTofDigitizerV.

Definition at line 31 of file BesTofDigitizerBrV1.cc.

{
    G4cout<<"BesTofDigitizerBrV1::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 BesTofDigitizerBrV1::Smear

(

G4int

scinNb

)

Definition at line 109 of file BesTofDigitizerBrV1.cc.

{
    G4cout<<"m_t1st:"<<m_t1st<<"  m_z:"<<m_z<<G4endl;
 
    /*G4double tofRes = 0.08;
    G4double atten;
    for(G4int i=0;i<2;i++)
    {
      m_TDC[i] += tofRes * G4RandGauss::shoot();
    }*/
 
    double pp1[10];
    double pp2[10];
    for (int i=0;i<10;i++)
    {
        pp1[i]=m_tofCaliSvc->BTof(scinNb)->getP1(i);
        pp2[i]=m_tofCaliSvc->BTof(scinNb)->getP2(i);
        //G4cout<<"pp1["<<i<<"]:"<<pp1[i]<<"  "<<"pp2["<<i<<"]:"<<pp2[i]<<G4endl;
    }
    double pp[10];
    for (int i=0;i<2;i++)
    {
        m_ADC[i] *= 7.;
        for (int j=0;j<10;j++)
        {
            if (i==0)
                pp[j]=pp1[j];
            else
                pp[j]=pp2[j];
        }
        m_TDC[i] += (pp[0]+pp[1]*m_z)/TMath::Sqrt(m_ADC[i])+
                    pp[2]*m_ADC[i]+
                    pp[3]*m_ADC[i]*m_ADC[i]+
                    pp[4]*m_ADC[i]*m_ADC[i]*m_ADC[i]+
                    pp[5]/(84.2*84.2+m_z*m_z)+
                    pp[6]*m_z+
                    pp[7]*m_z*m_z+
                    pp[8]*m_z*m_z*m_z+
                    pp[9];
    }
}

Referenced by Digitize().

TofPmtAccum()

void BesTofDigitizerBrV1::TofPmtAccum

(

BesTofHit *

hit

)

Definition at line 74 of file BesTofDigitizerBrV1.cc.

{
    G4int trackIndex = hit->GetTrackIndex();
    G4int partId = hit->GetPartId();
    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 posz = pos.z();
 
    G4double atten;
    atten = m_tofCaliSvc->BAtten(scinNb);
    G4double pathL[2];
    pathL[0]=m_scinLength/2-posz;
    pathL[1]=posz+m_scinLength/2;
    for (G4int j=0;j<2;j++)
    {
        //G4cout<<"atten:"<<atten<<G4endl;
        m_ADC[j] += edep*exp(-pathL[j]/atten);
    }
    if (time<m_t1st)
    {
        m_t1st = time;
        m_z = posz;
        m_TDC[0] = m_TDC[1] = m_t1st;
    }
}

Referenced by Digitize().

TofPmtInit()

void BesTofDigitizerBrV1::TofPmtInit

(

)

Definition at line 68 of file BesTofDigitizerBrV1.cc.

{
    Initialize();
    m_t1st = 9999.;
}

Referenced by Digitize().

---

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

• BesTofDigitizerBrV1.hh
• BesTofDigitizerBrV1.cc