ExGflashEventAction Class Reference

#include <ExGflashEventAction.hh>

Inheritance diagram for ExGflashEventAction:

Public Member Functions
	ExGflashEventAction ()
	~ExGflashEventAction ()
void	BeginOfEventAction (const G4Event *)
void	EndOfEventAction (const G4Event *)
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 39 of file ExGflashEventAction.hh.

Constructor & Destructor Documentation

ExGflashEventAction::ExGflashEventAction

(

)

Definition at line 49 of file ExGflashEventAction.cc.

                                        :
fNevent(0),fDtime(0.0),fCalorimeterCollectionId(-1)
{
}

ExGflashEventAction::~ExGflashEventAction

(

)

Definition at line 54 of file ExGflashEventAction.cc.

References G4cout, and G4endl.

{
  if ( fNevent > 0 ) G4cout << "Internal Real Elapsed Time /event is: "<< fDtime /fNevent<< G4endl;
}

Member Function Documentation

void ExGflashEventAction::BeginOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 60 of file ExGflashEventAction.cc.

References G4cout, G4endl, G4Event::GetEventID(), and G4Timer::Start().

{
  fTimerIntern.Start();
  G4cout<<" ------ Start ExGflashEventAction ----- "<<G4endl;
  fNevent=evt->GetEventID();
  G4cout<<" Start generating event Nr "<<fNevent<<G4endl<<G4endl;   
}

void ExGflashEventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Hits in sensitive Detector

Computing (x,y,z) of vertex of initial particles

For all Hits in sensitive detector

Reimplemented from G4UserEventAction.

Definition at line 68 of file ExGflashEventAction.cc.

References CLHEP::Hep3Vector::cross(), G4TrajectoryContainer::entries(), G4THitsCollection< T >::entries(), G4cout, G4endl, G4SDManager::GetCollectionID(), G4HCofThisEvent::GetHC(), G4Event::GetHCofThisEvent(), G4PrimaryParticle::GetMomentum(), G4PrimaryVertex::GetPosition(), G4PrimaryVertex::GetPrimary(), G4Event::GetPrimaryVertex(), G4Timer::GetRealElapsed(), G4SDManager::GetSDMpointer(), G4Timer::GetSystemElapsed(), G4Event::GetTrajectoryContainer(), G4Timer::GetUserElapsed(), python.hepunit::GeV, CLHEP::Hep3Vector::mag(), G4INCL::Math::max(), G4Timer::Stop(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

{  
  fTimerIntern.Stop();
  G4cout << G4endl;
  G4cout << "******************************************";
  G4cout << G4endl;
  G4cout << "Internal Real Elapsed Time is: "<< fTimerIntern.GetRealElapsed();
  G4cout << G4endl;
  G4cout << "Internal System Elapsed Time: " << fTimerIntern.GetSystemElapsed();
  G4cout << G4endl;
  G4cout << "Internal GetUserElapsed Time: " << fTimerIntern.GetUserElapsed();
  G4cout << G4endl;
  G4cout << "******************************************"<< G4endl;
  fDtime+=fTimerIntern.GetRealElapsed();
  G4cout<<" ------ ExGflashEventAction::End of event nr. "<<fNevent<<"  -----"<< G4endl;   
  
  G4SDManager * SDman = G4SDManager::GetSDMpointer();
  G4String colNam;
  fCalorimeterCollectionId=SDman->GetCollectionID(colNam="ExGflashCollection");
  if (fCalorimeterCollectionId<0) return;
  G4HCofThisEvent * HCE = evt->GetHCofThisEvent();
  ExGflashHitsCollection* THC = 0;
  G4double totE = 0;
  // Read out of the crysta ECAL
  THC=(ExGflashHitsCollection *)(HCE->GetHC(fCalorimeterCollectionId));
  if (THC)
    {
      /// Hits in sensitive Detector
      int n_hit = THC->entries();
      G4cout<<"  " << n_hit<< " hits are stored in ExGflashHitsCollection "<<G4endl;
      G4PrimaryVertex* pvertex=evt->GetPrimaryVertex();   
      ///Computing (x,y,z) of vertex of initial particles  
      G4ThreeVector vtx=pvertex->GetPosition();
    G4PrimaryParticle* pparticle=pvertex->GetPrimary();
    // direction of the Shower
    G4ThreeVector mom=pparticle->GetMomentum()/pparticle->GetMomentum().mag();
    
    //@@@ ExGflashEventAction: Magicnumber
    G4double energyincrystal[100];
    for (int i=0;i<100;i++) energyincrystal[i]=0.;
    
    //@@@ ExGflashEventAction: Magicnumber
    /// For all Hits in sensitive detector
    for (int i=0;i<n_hit;i++)
      {
        G4double estep = (*THC)[i]->GetEdep()/GeV;
        if (estep >0.0)
          {
            totE += (*THC)[i]->GetEdep()/GeV;
            G4int num=(*THC)[i]->GetCrystalNum();
            
            energyincrystal[num]+=(*THC)[i]->GetEdep()/GeV;
            //G4cout << num << G4endl;
            //  G4cout << " Crystal Nummer " <<  (*THC)[i]->GetCrystalNum()  << G4endl;
            //  G4cout <<  (*THC)[i]->GetCrystalNum() /10 <<
            // "  "<<(*THC)[i]->GetCrystalNum()%10 << G4endl;
            
            G4ThreeVector hitpos=(*THC)[i]->GetPos();            
            G4ThreeVector l (hitpos.x(), hitpos.y(), hitpos.z());
            // distance from shower start
            l = l - vtx; 
            // projection on shower axis = longitudinal profile
            G4ThreeVector longitudinal  =  l;  
            // shower profiles (Radial)
            G4ThreeVector radial = vtx.cross(l);
          }
      }
    G4double   max = 0;
    G4int    index = 0;
    //Find crystal with maximum energy
    for (int i=0;i<100;i++) 
      {
        //G4cout << i <<"  " << energyincrystal[i] << G4endl;
        if (max <energyincrystal[i])
          {
            max=energyincrystal[i];
            index = i;
          }
      }  
    //G4cout << index <<" NMAX  " << index << G4endl;  
    
    // 3x3 matrix of crystals around the crystal with the maximum energy deposit
    G4double e3x3 =
      energyincrystal[index]+energyincrystal[index+1]+energyincrystal[index-1]+
      energyincrystal[index-10]+energyincrystal[index-9]+energyincrystal[index-11]+
      energyincrystal[index+10]+energyincrystal[index+11]+energyincrystal[index+9];
    
    // 5x5 matrix of crystals around the crystal with the maximum energy deposit  
    G4double e5x5 =
      energyincrystal[index]+energyincrystal[index+1]+energyincrystal[index-1]+
      energyincrystal[index+2]+energyincrystal[index-2]+
      energyincrystal[index-10]+energyincrystal[index-9]+energyincrystal[index-11]+
      energyincrystal[index-8]+energyincrystal[index-12]+
      energyincrystal[index+10]+energyincrystal[index+11]+energyincrystal[index+9]+
      energyincrystal[index+12]+energyincrystal[index+8];
    
    G4cout << "   e1  " << energyincrystal[index]  
           << "   e3x3  " << e3x3<<  "   GeV  e5x5  "   <<e5x5  <<G4endl;  
  }
  
  G4cout << " Total energy deposited in the calorimeter: " << totE << " (GeV)" << G4endl;
  G4TrajectoryContainer * trajectoryContainer = evt->GetTrajectoryContainer();
  G4int n_trajectories = 0;
  if(trajectoryContainer){ n_trajectories = trajectoryContainer->entries(); }
  G4cout << "    " << n_trajectories  << " trajectories stored in this event." << G4endl;
  
}

---

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

• ExGflashEventAction.hh
• ExGflashEventAction.cc