SimEnergyDepositSD.cc

Go to the documentation of this file.

//
//               __        __ __  __  __
//   ____ ______/ /_____ _/ // / / /_/ /__
//  / __ `/ ___/ __/ __ `/ // /_/ __/ //_/
// / /_/ / /  / /_/ /_/ /__  __/ /_/ ,<
// \__,_/_/   \__/\__, /  /_/  \__/_/|_|
//               /____/
//
// artg4tk: art based Geant 4 Toolkit
//
//=============================================================================
// SimEnergyDeposit.cc: Class representing a sensitive tracking detector
// Author: Hans Wenzel (Fermilab)
//=============================================================================
#include "larg4/Services/SimEnergyDepositSD.h"
#include "Geant4/G4Cerenkov.hh"
#include "Geant4/G4Event.hh"
#include "Geant4/G4EventManager.hh"
#include "Geant4/G4HCofThisEvent.hh"
#include "Geant4/G4SDManager.hh"
#include "Geant4/G4Scintillation.hh"
#include "Geant4/G4Step.hh"
#include "Geant4/G4SteppingManager.hh"
#include "Geant4/G4ThreeVector.hh"
#include "Geant4/G4VSolid.hh"
#include "Geant4/G4VVisManager.hh"
#include "Geant4/G4ios.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
namespace larg4 {

  SimEnergyDepositSD::SimEnergyDepositSD(G4String name) : G4VSensitiveDetector(name)
  {
    hitCollection.clear();
  }

  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  SimEnergyDepositSD::~SimEnergyDepositSD() {}

  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  void SimEnergyDepositSD::Initialize(G4HCofThisEvent* HCE)
  {
    hitCollection.clear();
  }
  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  G4bool SimEnergyDepositSD::ProcessHits(G4Step* aStep, G4TouchableHistory*)
  {
    G4double edep = aStep->GetTotalEnergyDeposit() / CLHEP::MeV;

    if (edep == 0.) return false;
    //std::cout << "7777777777777777:   "<< aStep->GetTotalEnergyDeposit()/CLHEP::MeV << "   " << aStep->GetTotalEnergyDeposit() <<std::endl;
    const int electronsperMeV = 10000;
    int nrelec = (int)round(edep * electronsperMeV);
    if (aStep->GetTrack()->GetDynamicParticle()->GetCharge() == 0) return false;
    G4int photons = 0;
    G4SteppingManager* fpSteppingManager =
      G4EventManager::GetEventManager()->GetTrackingManager()->GetSteppingManager();
    G4StepStatus stepStatus = fpSteppingManager->GetfStepStatus();
    if (stepStatus != fAtRestDoItProc) {
      G4ProcessVector* procPost = fpSteppingManager->GetfPostStepDoItVector();
      size_t MAXofPostStepLoops = fpSteppingManager->GetMAXofPostStepLoops();
      for (size_t i3 = 0; i3 < MAXofPostStepLoops; i3++) {
        if (!(*procPost)[i3]) continue;
        /*
             if ((*procPost)[i3]->GetProcessName() == "Cerenkov") {
             G4Cerenkov* proc =(G4Cerenkov*) (*procPost)[i3];
             photons+=proc->GetNumPhotons();
             }
           */
        if ((*procPost)[i3]->GetProcessName() == "Scintillation") {
          G4Scintillation* proc1 = (G4Scintillation*)(*procPost)[i3];
          photons += proc1->GetNumPhotons();
        }
      }
    }
    geo::Point_t start = geo::Point_t(aStep->GetPreStepPoint()->GetPosition().x() / CLHEP::cm,
                                      aStep->GetPreStepPoint()->GetPosition().y() / CLHEP::cm,
                                      aStep->GetPreStepPoint()->GetPosition().z() / CLHEP::cm);
    geo::Point_t end = geo::Point_t(aStep->GetPostStepPoint()->GetPosition().x() / CLHEP::cm,
                                    aStep->GetPostStepPoint()->GetPosition().y() / CLHEP::cm,
                                    aStep->GetPostStepPoint()->GetPosition().z() / CLHEP::cm);
    sim::SimEnergyDeposit newHit =
      sim::SimEnergyDeposit(photons,
                            nrelec,
                            1.0,
                            edep,
                            start,
                            end,
                            aStep->GetPreStepPoint()->GetGlobalTime() / CLHEP::ns,
                            aStep->GetPostStepPoint()->GetGlobalTime() / CLHEP::ns,
                            aStep->GetTrack()->GetTrackID(),
                            aStep->GetTrack()->GetParticleDefinition()->GetPDGEncoding(),
                            aStep->GetTrack()->GetTrackID() //original track id
      );
    hitCollection.push_back(newHit);
    return true;
  } // end ProcessHits
} // end namespace  larg4