GammaRayTelEventAction Class Reference

#include <GammaRayTelEventAction.hh>

Inheritance diagram for GammaRayTelEventAction:

Public Member Functions
	GammaRayTelEventAction ()
virtual	~GammaRayTelEventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
void	SetDrawFlag (G4String val)
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
void	SetEventManager (G4EventManager *value)

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Definition at line 52 of file GammaRayTelEventAction.hh.

Constructor & Destructor Documentation

GammaRayTelEventAction::GammaRayTelEventAction

(

)

Definition at line 75 of file GammaRayTelEventAction.cc.

References G4DigiManager::AddNewModule(), and G4DigiManager::GetDMpointer().

  :trackerCollID(-1),calorimeterCollID(-1),                
  anticoincidenceCollID(-1), drawFlag("all")
{ 
  G4DigiManager * fDM = G4DigiManager::GetDMpointer();
  GammaRayTelDigitizer * myDM = new GammaRayTelDigitizer( "GammaRayTelDigitizer" );
  fDM->AddNewModule(myDM);
}

GammaRayTelEventAction::~GammaRayTelEventAction ( )

virtual

Definition at line 86 of file GammaRayTelEventAction.cc.

{
}

Member Function Documentation

void GammaRayTelEventAction::BeginOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 92 of file GammaRayTelEventAction.cc.

References G4cout, G4endl, G4SDManager::GetCollectionID(), G4Event::GetEventID(), and G4SDManager::GetSDMpointer().

{

  G4int evtNb = evt->GetEventID();
  G4cout << "Event: " << evtNb << G4endl;
  G4SDManager * SDman = G4SDManager::GetSDMpointer();  

  if (trackerCollID==-1) {
    trackerCollID = SDman->GetCollectionID("TrackerCollection");
  }
  if(anticoincidenceCollID==-1) {
    anticoincidenceCollID =
      SDman->GetCollectionID("AnticoincidenceCollection");
  }
  if(calorimeterCollID==-1) {
    calorimeterCollID =
      SDman->GetCollectionID("CalorimeterCollection");
  }
}

void GammaRayTelEventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 114 of file GammaRayTelEventAction.cc.

References GammaRayTelDigitizer::Digitize(), G4TDigiCollection< T >::entries(), G4THitsCollection< T >::entries(), G4DigiManager::FindDigitizerModule(), G4cout, G4endl, G4DigiManager::GetDigiCollection(), G4DigiManager::GetDigiCollectionID(), G4DigiManager::GetDMpointer(), G4Event::GetEventID(), G4HCofThisEvent::GetHC(), G4Event::GetHCofThisEvent(), python.hepunit::keV, python.hepunit::mm, outFile, test::x, and z.

{
  G4int event_id = evt->GetEventID();
  
  
  G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
  GammaRayTelTrackerHitsCollection* THC = 0;
//  GammaRayTelCalorimeterHitsCollection* CHC = 0;
//  GammaRayTelAnticoincidenceHitsCollection* AHC = 0;


  G4DigiManager * fDM = G4DigiManager::GetDMpointer();

  if (HCE)
    {
      THC = (GammaRayTelTrackerHitsCollection*)(HCE->GetHC(trackerCollID));
//      CHC = (GammaRayTelCalorimeterHitsCollection*)
//  (HCE->GetHC(calorimeterCollID));
//      AHC = (GammaRayTelAnticoincidenceHitsCollection*)
//  (HCE->GetHC(anticoincidenceCollID));
      
      if (THC)
    {
      int n_hit = THC->entries();
      G4cout << "Number of tracker hits in this event =  " << n_hit << G4endl;
      G4double ESil=0;
      G4int NStrip, NPlane, IsX;
      
      // This is a cycle on all the tracker hits of this event
      
      for (int i=0;i<n_hit;i++) 
        {
          // Here we put the hit data in a an ASCII file for 
          // later analysis 
          ESil = (*THC)[i]->GetEdepSil();
          NStrip = (*THC)[i]->GetNStrip();
          NPlane = (*THC)[i]->GetNSilPlane();
          IsX = (*THC)[i]->GetPlaneType();
          
#ifdef G4STORE_DATA
          outFile << std::setw(7) << event_id << " " << 
        ESil/keV << " " << NStrip << 
        " " << NPlane << " " << IsX << " " <<
        (*THC)[i]->GetPos().x()/mm <<" "<<
        (*THC)[i]->GetPos().y()/mm <<" "<<
        (*THC)[i]->GetPos().z()/mm <<" "<<
        G4endl;
#else     
          G4cout << std::setw(7) << event_id << " " << 
        ESil/keV << " " << NStrip << 
        " " << NPlane << " " << IsX << " " <<
        (*THC)[i]->GetPos().x()/mm <<" "<<
        (*THC)[i]->GetPos().y()/mm <<" "<<
        (*THC)[i]->GetPos().z()/mm <<" "<<
        G4endl;
#endif
          
#ifdef G4ANALYSIS_USE
          
          // Here we fill the histograms of the Analysis manager
          GammaRayTelAnalysis* analysis = GammaRayTelAnalysis::getInstance();
          
          if(IsX) 
        {
          if (analysis->GetHisto2DMode()=="position")
            analysis->InsertPositionXZ((*THC)[i]->GetPos().x()/mm,(*THC)[i]->GetPos().z()/mm);
          else
            analysis->InsertPositionXZ(NStrip, NPlane);           
          if (NPlane == 0) analysis->InsertEnergy(ESil/keV);
          analysis->InsertHits(NPlane);
        } 
          else 
        {
          if (analysis->GetHisto2DMode()=="position")
            analysis->InsertPositionYZ((*THC)[i]->GetPos().y()/mm,(*THC)[i]->GetPos().z()/mm);  
          else 
            analysis->InsertPositionYZ(NStrip, NPlane);           
          if (NPlane == 0) analysis->InsertEnergy(ESil/keV);
          analysis->InsertHits(NPlane);
        }
          
#ifdef G4ANALYSIS_USE
          analysis->setNtuple( ESil/keV, NPlane, (*THC)[i]->GetPos().x()/mm,
                   (*THC)[i]->GetPos().y()/mm,
                   (*THC)[i]->GetPos().z()/mm);
#endif
          
#endif
      
        }
#ifdef G4ANALYSIS_USE
      GammaRayTelAnalysis* analysis = GammaRayTelAnalysis::getInstance();
      analysis->EndOfEvent(n_hit);
#endif
      
    }
      
      GammaRayTelDigitizer * myDM = 
    (GammaRayTelDigitizer*)fDM->FindDigitizerModule( "GammaRayTelDigitizer" );
      myDM->Digitize();
      
#ifdef write_outfile
      // the whole block is needed only when outfile is active; protect block to avoid
      //     compilations warnings from gcc4.6, Gunter Folger

      G4int myDigiCollID = fDM->GetDigiCollectionID("DigitsCollection");
      
      // G4cout << "digi collecion" << myDigiCollID << G4endl;
      
      GammaRayTelDigitsCollection * DC = (GammaRayTelDigitsCollection*)fDM->GetDigiCollection( myDigiCollID );
      
      if(DC) {
    //    G4cout << "Total Digits " << DC->entries() << G4endl;
    G4int n_digi =  DC->entries();
    G4int NStrip, NPlane, IsX;
    for (G4int i=0;i<n_digi;i++) {
      // Here we put the digi data in a an ASCII file for 
      // later analysis
      NStrip = (*DC)[i]->GetStripNumber();
      NPlane = (*DC)[i]->GetPlaneNumber();
      IsX = (*DC)[i]->GetPlaneType();
      
      outFile << std::setw(7) << event_id << " " << NStrip << 
        " " << NPlane << " " << IsX << " " << G4endl;   
    }
      }
#endif      
    }
}

void GammaRayTelEventAction::SetDrawFlag ( G4String  val )

inline

Definition at line 63 of file GammaRayTelEventAction.hh.

{drawFlag = val;};

---

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

• GammaRayTelEventAction.hh
• GammaRayTelEventAction.cc