CCalEndOfEventAction Class Reference

#include <CCalEndOfEventAction.hh>

Inheritance diagram for CCalEndOfEventAction:

Collaboration diagram for CCalEndOfEventAction:

Public Member Functions
	CCalEndOfEventAction (CCalPrimaryGeneratorAction *)
	~CCalEndOfEventAction ()
virtual void	StartOfEventAction (const G4Event *evt)
virtual void	EndOfEventAction (const G4Event *evt)
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
virtual void	SetEventManager (G4EventManager *value)
virtual void	BeginOfEventAction (const G4Event *anEvent)

Private Member Functions
void	instanciateSteppingAction ()
void	initialize ()

Private Attributes
G4bool	isInitialized
CCalPrimaryGeneratorAction *	primaryGenerator
CCalSteppingAction *	theSteppingAction
nameType *	SDnames
G4int	numberOfSD
CCaloOrganization *	theOrg

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Definition at line 44 of file CCalEndOfEventAction.hh.

Constructor & Destructor Documentation

CCalEndOfEventAction()

CCalEndOfEventAction::CCalEndOfEventAction

(

CCalPrimaryGeneratorAction *

pg

)

Definition at line 57 of file CCalEndOfEventAction.cc.

                                                                         : 
  isInitialized(false),SDnames(0),numberOfSD(0) {

  primaryGenerator = pg;
#ifdef debug
  G4cout << "Instantiate CCalEndOfEventAction" << G4endl;
#endif

  G4cout << "Now Instantiate stepping action" << G4endl;
  instanciateSteppingAction();
  
  G4cout << "Get Calorimter organisation" << G4endl;
  theOrg = new CCaloOrganization;
  G4cout << "end of instantiation of EndofEventAction" << G4endl;
}

Here is the call graph for this function:

~CCalEndOfEventAction()

CCalEndOfEventAction::~CCalEndOfEventAction

(

)

Definition at line 74 of file CCalEndOfEventAction.cc.

                                            {

  if (theOrg)
    delete theOrg;
  if (SDnames)
    delete[] SDnames;
  G4cout << "CCalEndOfEventAction deleted" << G4endl;
}

Member Function Documentation

EndOfEventAction()

void CCalEndOfEventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 112 of file CCalEndOfEventAction.cc.

                                                             {

#ifdef debug
  G4cout << G4endl << "=== Begin of EndOfEventAction === " << G4endl;
#endif

  if (!isInitialized) initialize();

  theSteppingAction->endOfEvent();
  
  //
  // Look for the Hit Collection 
  //  
  G4HCofThisEvent* allHC = evt->GetHCofThisEvent();
  if (allHC == 0) {
#ifdef debug
    G4cout << "CCalEndOfEventAction: No Hit Collection in this event" 
     << G4endl;
#endif
    return;
  }
    
  //
  // hits info
  //
  
  //Now make summary
  float hcalE[28], ecalE[49], fullE=0., edec=0, edhc=0;
  int i = 0;
  for (i = 0; i < 28; i++) {hcalE[i]=0.;}
  for (i = 0; i < 49; i++) {ecalE[i]=0.;}

  float* edep = new float[numberOfSD];
  int nhit=0;
  for (i = 0; i < numberOfSD; i++){

    //
    // Look for the Hit Collection
    //
    edep[i] = 0;
    int caloHCid = G4SDManager::GetSDMpointer()->GetCollectionID(SDnames[i]);

    if (caloHCid >= 0) {
      CCalG4HitCollection* theHC = 
    (CCalG4HitCollection*) allHC->GetHC(caloHCid);
    
      if (theHC != 0) {

    G4int nentries = theHC->entries();
#ifdef debug
    G4cout << " There are " << nentries << " hits in " << SDnames[i] 
           << " :" << G4endl;
#endif

    if (nentries > 0) {
  
      int j;
      for (j=0; j<nentries; j++){
#ifdef ddebug
        G4cout << "Hit " << j;
#endif
        CCalG4Hit* aHit =  (*theHC)[j];
        float En = aHit->getEnergyDeposit();
        int unitID = aHit->getUnitID();
        int id=-1;
        if (unitID > 0 && unitID < 29) {
          id = unitID - 1; // HCal
          hcalE[id] += En/GeV;
        } else {
          int i0 = unitID/4096;
          int i1 = (unitID/64)%64;
          int i2 = unitID%64;
          if (i0 == 1 && i1 < 8 && i2 < 8) {
        id = i1*7 + i2; // ECal
        ecalE[id] += En/GeV;
          }
        }
#ifdef ddebug
        G4cout << " with Energy = " << En/MeV << " MeV in Unit " << unitID 
           << " " << id << G4endl;
#endif
        fullE   += En/GeV;
        edep[i] += En/GeV;
        nhit++;
      }
#ifdef ddebug
      G4cout << " ===> Total Energy Deposit in this Calorimeter = " 
         << edep[i]*1000.0 << "  MeV " << G4endl; 
#endif
    }
      }
    }
    if (SDnames[i] == "HadronCalorimeter") {
      edhc = edep[i];
    } else if (SDnames[i] == "CrystalMatrix") {
      edec = edep[i];
    }
  }

  delete[] edep;

  G4ThreeVector pos = primaryGenerator->GetParticlePosition();
  float ener = primaryGenerator->GetParticleEnergy()/GeV;
  float x    = pos.x()/mm;
  float y    = pos.y()/mm;
  float z    = pos.z()/mm;

  CCalAnalysis* analysis = CCalAnalysis::getInstance();
  analysis->InsertEnergy(fullE);
  analysis->InsertEnergyHcal(hcalE);
  analysis->InsertEnergyEcal(ecalE);
  analysis->setNtuple(hcalE, ecalE, ener, x, y, z, fullE, edec, edhc);
  analysis->EndOfEvent(nhit);
  for (i = 0; i < numberOfSD; i++){
    int caloHCid = G4SDManager::GetSDMpointer()->GetCollectionID(SDnames[i]);
    if (caloHCid >= 0) {
      CCalG4HitCollection* theHC = 
    (CCalG4HitCollection*) allHC->GetHC(caloHCid);
      if (theHC != 0) {
    G4int nentries = theHC->entries();
    if (nentries > 0) {
      for (G4int k=0; k<nentries; k++) {
        CCalG4Hit* aHit =  (*theHC)[k];
        analysis->InsertTimeProfile(k,aHit->getTimeSlice(),
                    aHit->getEnergyDeposit()/GeV);
      }
    }
      }
    }
  }
}

Here is the call graph for this function:

initialize()

void CCalEndOfEventAction::initialize ( )

private

Definition at line 84 of file CCalEndOfEventAction.cc.

                                      {

  isInitialized = true;
  numberOfSD = CCalSDList::getInstance()->getNumberOfCaloSD();
#ifdef debug
  G4cout << "CCalEndOfEventAction look for " << numberOfSD 
       << " calorimeter-like SD" << G4endl;
#endif
  if (numberOfSD > 0) {
    G4int n = numberOfSD;
    if (n < 2) n = 2;
    SDnames = new nameType[n];
  }
  for (int i=0; i<numberOfSD; i++) {
    SDnames[i] = G4String(CCalSDList::getInstance()->getCaloSDName(i));
#ifdef debug
    G4cout << "CCalEndOfEventAction: found SD " << i << " name "
     << SDnames[i] << G4endl;
#endif
  }       
}

Here is the call graph for this function:

Here is the caller graph for this function:

instanciateSteppingAction()

void CCalEndOfEventAction::instanciateSteppingAction ( )

private

Definition at line 245 of file CCalEndOfEventAction.cc.

                                                    {
    
  G4UserSteppingAction* theUA = const_cast<G4UserSteppingAction*>(G4RunManager::GetRunManager()->GetUserSteppingAction());
        
  if (theUA == 0) {
#ifdef debug
    G4cout << " CCalEndOfEventAction::instanciateSteppingAction creates"
     << " CCalSteppingAction" << G4endl;
#endif
    theSteppingAction = new CCalSteppingAction;  
    G4RunManager::GetRunManager()->SetUserAction(theSteppingAction);
  }   
    
}

Here is the call graph for this function:

Here is the caller graph for this function:

StartOfEventAction()

void CCalEndOfEventAction::StartOfEventAction ( const G4Event *  evt )

virtual

Definition at line 107 of file CCalEndOfEventAction.cc.

                                                                { 
  G4cout << "\n---> Begin of event: " << evt->GetEventID() << G4endl;
}

Here is the call graph for this function:

Member Data Documentation

isInitialized

G4bool CCalEndOfEventAction::isInitialized

private

Definition at line 57 of file CCalEndOfEventAction.hh.

numberOfSD

G4int CCalEndOfEventAction::numberOfSD

private

Definition at line 62 of file CCalEndOfEventAction.hh.

primaryGenerator

CCalPrimaryGeneratorAction* CCalEndOfEventAction::primaryGenerator

private

Definition at line 59 of file CCalEndOfEventAction.hh.

SDnames

nameType* CCalEndOfEventAction::SDnames

private

Definition at line 61 of file CCalEndOfEventAction.hh.

theOrg

CCaloOrganization* CCalEndOfEventAction::theOrg

private

Definition at line 63 of file CCalEndOfEventAction.hh.

theSteppingAction

CCalSteppingAction* CCalEndOfEventAction::theSteppingAction

private

Definition at line 60 of file CCalEndOfEventAction.hh.

---

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

• CCalEndOfEventAction.hh
• CCalEndOfEventAction.cc