GammaRayTelDigitizer Class Reference

#include <GammaRayTelDigitizer.hh>

Inheritance diagram for GammaRayTelDigitizer:

Collaboration diagram for GammaRayTelDigitizer:

Public Member Functions
	GammaRayTelDigitizer (G4String name)
	~GammaRayTelDigitizer ()
void	Digitize ()
void	SetThreshold (G4double val)
Public Member Functions inherited from G4VDigitizerModule
	G4VDigitizerModule (G4String modName)
virtual	~G4VDigitizerModule ()
int	operator== (const G4VDigitizerModule &right) const
int	operator!= (const G4VDigitizerModule &right) const
G4int	GetNumberOfCollections () const
G4String	GetCollectionName (G4int i) const
G4String	GetName () const
void	SetVerboseLevel (G4int val)

Additional Inherited Members
Protected Member Functions inherited from G4VDigitizerModule
void	StoreDigiCollection (G4VDigiCollection *aDC)
void	StoreDigiCollection (G4int DCID, G4VDigiCollection *aDC)
Protected Attributes inherited from G4VDigitizerModule
G4DigiManager *	DigiManager
G4String	moduleName
std::vector< G4String >	collectionName
G4int	verboseLevel

Detailed Description

Definition at line 51 of file GammaRayTelDigitizer.hh.

Constructor & Destructor Documentation

GammaRayTelDigitizer::GammaRayTelDigitizer

(

G4String

name

)

Definition at line 57 of file GammaRayTelDigitizer.cc.

  :G4VDigitizerModule(name)
{


  G4String colName = "DigitsCollection";
  collectionName.push_back(colName);
  Energy_Threshold = 120.*keV;
  TotalEnergy = 0.;
  ACDThreshold = 15*keV;

//create a messenger for this class
  digiMessenger = new GammaRayTelDigitizerMessenger(this);

}

GammaRayTelDigitizer::~GammaRayTelDigitizer

(

)

Definition at line 75 of file GammaRayTelDigitizer.cc.

{
  delete digiMessenger;
}

Member Function Documentation

void GammaRayTelDigitizer::Digitize ( void  )

virtual

Implements G4VDigitizerModule.

Definition at line 82 of file GammaRayTelDigitizer.cc.

{
  DigitsCollection = new GammaRayTelDigitsCollection
    ("GammaRayTelDigitizer","DigitsCollection"); // to create the Digi Collection
  
  G4DigiManager* DigiMan = G4DigiManager::GetDMpointer();
  
  G4int THCID; // TrackerCollection
  G4int CHCID; // CalorimeterCollection
  G4int AHCID; // AnticoincidenceCollection

  TotalEnergy = 0.; // at each event
  
  // TKR Hits collection
  
  THCID = DigiMan->GetHitsCollectionID("TrackerCollection");  
  GammaRayTelTrackerHitsCollection* THC = 0; 
  THC = (GammaRayTelTrackerHitsCollection*)
    (DigiMan->GetHitsCollection(THCID));


  // CAL Hits collection

  CHCID = DigiMan->GetHitsCollectionID("CalorimeterCollection");
  GammaRayTelCalorimeterHitsCollection* CHC = 0;
  CHC = (GammaRayTelCalorimeterHitsCollection*)
         (DigiMan->GetHitsCollection(CHCID));

  // ACD Hits collection

  AHCID = DigiMan->GetHitsCollectionID("AnticoincidenceCollection");  
  GammaRayTelAnticoincidenceHitsCollection* AHC = 0;
  AHC = (GammaRayTelAnticoincidenceHitsCollection*)
         (DigiMan->GetHitsCollection(AHCID));

  
  if (THC)
    {
      G4int n_hit = THC->entries();

      for (G4int i=0;i<n_hit;i++)
    {
      G4double ESil = (*THC)[i]->GetEdepSil();
      G4int NStrip = (*THC)[i]->GetNStrip();
      G4int NPlane = (*THC)[i]->GetNSilPlane();
      G4int IsX = (*THC)[i]->GetPlaneType();
      
       // digi generation only if energy deposit greater than threshold

      if (ESil>Energy_Threshold)
        {
          GammaRayTelDigi* Digi = new GammaRayTelDigi();
          Digi->SetPlaneNumber(NPlane);
          Digi->SetPlaneType(IsX);
          Digi->SetStripNumber(NStrip);
          Digi->SetDigiType(0);
          Digi->SetEnergy(0.);
          DigitsCollection->insert(Digi);   
        }  
    }
    }
  if (CHC)
    {
      G4int n_hit = CHC->entries();
      
      for (G4int i=0;i<n_hit;i++)
    {
      TotalEnergy +=(*CHC)[i]->GetEdepCAL();
    }
      // digi generation only if energy deposit greater than 0.
      
      if (TotalEnergy>0.)
    {
      GammaRayTelDigi* Digi = new GammaRayTelDigi();
      Digi->SetPlaneNumber(0);
      Digi->SetPlaneType(0);
      Digi->SetStripNumber(0);
      Digi->SetDigiType(1);
      Digi->SetEnergy(TotalEnergy);
      DigitsCollection->insert(Digi);   
    }  

    }

  if (AHC)
    {
      G4int n_hit = AHC->entries();
      
      for (G4int i=0;i<n_hit;i++)
    {
      G4double energy = (*AHC)[i]->GetEdepACD();
      G4int type = (*AHC)[i]->GetACDTileNumber();
      
      // digi generation only if energy deposit greater than 0.
      
      if (energy>ACDThreshold)
        {
          GammaRayTelDigi* Digi = new GammaRayTelDigi();
          Digi->SetPlaneNumber(0);
          Digi->SetPlaneType(0);
          Digi->SetStripNumber(type);
          Digi->SetDigiType(2);
          Digi->SetEnergy(energy);
          DigitsCollection->insert(Digi);   
        }  
      
    }
    }
  
  if (THC||AHC||CHC){
    G4cout << "Number of digits in this event =  " 
       << DigitsCollection->entries()  
      // << " " << DigitsCollection->GetName()  
      // << " " << DigitsCollection->GetDMname() 
       << G4endl;
  }
  
  StoreDigiCollection(DigitsCollection);
  
  G4int DCID = -1;
  if(DCID<0)
    { 
      //      DigiMan->List();
      DCID = DigiMan->GetDigiCollectionID("GammaRayTelDigitizer/DigitsCollection");
    }
  
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelDigitizer::SetThreshold ( G4double  val )

inline

Definition at line 59 of file GammaRayTelDigitizer.hh.

{ Energy_Threshold = val;}

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/gammaray_telescope/include/GammaRayTelDigitizer.hh
• source/geant4.10.03.p03/examples/advanced/gammaray_telescope/src/GammaRayTelDigitizer.cc