#include <BesMucDigitizer.hh>

Inheritance diagram for BesMucDigitizer:

Public Member Functions
	BesMucDigitizer (G4String modName)

	~BesMucDigitizer ()

virtual void	Digitize ()

Detailed Description

Definition at line 23 of file BesMucDigitizer.hh.

Constructor & Destructor Documentation

◆ BesMucDigitizer()

BesMucDigitizer::BesMucDigitizer ( G4String modName )

Definition at line 31 of file BesMucDigitizer.cc.

:G4VDigitizerModule(modName)
{
  collectionName.push_back("BesMucDigisCollection");
  m_besMucDigisCollection=0;
 
  //Standard unit: length in mm; time in ns; energy in MeV; att in mm;
  //Here: att in mm; vel in m/s; threshold in MeV; mucRes in ps; 
  
 
  //retrieve G4Svc
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  IG4Svc* tmpSvc;
  StatusCode sc=svcLocator->service("G4Svc", tmpSvc);
  m_G4Svc=dynamic_cast<G4Svc *>(tmpSvc);
 
  //get Muc Ntuple from G4Svc
  if(m_G4Svc->MucRootFlag())
  {
    m_tupleMuc = m_G4Svc->GetTupleMuc();
    sc = m_tupleMuc->addItem("part",m_part);
    sc = m_tupleMuc->addItem("seg",m_seg);
    sc = m_tupleMuc->addItem("gap",m_gap);
    sc = m_tupleMuc->addItem("strip",m_strip);
 
  }
 
 
 
 
  Initialize();  
}

◆ ~BesMucDigitizer()

BesMucDigitizer::~BesMucDigitizer ( )

Definition at line 64 of file BesMucDigitizer.cc.

65{

66}

Member Function Documentation

◆ Digitize()

void BesMucDigitizer::Digitize ( )

virtual

Definition at line 82 of file BesMucDigitizer.cc.

{
  Initialize();
  
  m_besMucDigisCollection = new BesMucDigisCollection(moduleName, collectionName[0]);
  G4DigiManager* DigiMan = G4DigiManager::GetDMpointer();
  
  //hits collection ID
  G4int THCID;
  THCID = DigiMan->GetHitsCollectionID("BesMucHitsCollection");
  
  //hits collection
  BesMucHitsCollection* THC = 0;
  THC = (BesMucHitsCollection*) (DigiMan->GetHitsCollection(THCID));
 
  if(THC) {
    G4int trackIndex, part, seg, gap, strip;
    BesMucDigit *aDigit = new BesMucDigit();
 
    G4int n_hit = THC->entries();
    for(G4int i = 0; i < n_hit; i++) {
      aDigit->SetHit((*THC)[i]);
      trackIndex = aDigit->GetTrackIndex();
      part  = aDigit->GetPart();
      seg   = aDigit->GetSeg();
      gap   = aDigit->GetGap();
      strip = aDigit->GetNearestStripNo();
     
      //fill muc Ntuple
      if(m_G4Svc->MucRootFlag())
    {
      m_part = part;
      m_seg = seg;
      m_gap = gap;
      m_strip = strip;
      
      m_tupleMuc->write();
    }
 
      (*THC)[i]->SetPart(part);
      (*THC)[i]->SetSeg(seg);
      (*THC)[i]->SetGap(gap);
      (*THC)[i]->SetStrip(strip);
      if(m_stripID[part][seg][gap][strip] == 0) {
    NewDigi(trackIndex, part, seg, gap, strip);
      }
      else { 
    OldDigi();
      }
    } 
    //Smear();  
    StoreDigiCollection(m_besMucDigisCollection);   
    delete aDigit;
  }
}