larg4::AuxDetSD Class Reference

#include "AuxDetSD.h"

Inheritance diagram for larg4::AuxDetSD:

Public Member Functions
	AuxDetSD (G4String name)
virtual	~AuxDetSD ()
void	Initialize (G4HCofThisEvent *)
void	EndOfEvent (G4HCofThisEvent *)
G4bool	ProcessHits (G4Step *, G4TouchableHistory *)
const sim::AuxDetHitCollection &	GetHits () const

Private Attributes
TempHitCollection	temphitCollection
sim::AuxDetHitCollection	hitCollection

Detailed Description

Definition at line 26 of file AuxDetSD.h.

Constructor & Destructor Documentation

larg4::AuxDetSD::AuxDetSD

(

G4String

name

)

Definition at line 33 of file AuxDetSD.cc.

References hitCollection.

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

larg4::AuxDetSD::~AuxDetSD ( )

virtual

Definition at line 40 of file AuxDetSD.cc.

{}

Member Function Documentation

void larg4::AuxDetSD::EndOfEvent

(

G4HCofThisEvent *

)

Definition at line 95 of file AuxDetSD.cc.

References util::counter(), hitCollection, and temphitCollection.

  {
    if (temphitCollection.size() == 0) return; // No hits so nothing to do
#if defined _verbose_
    std::cout << " EndOfEvent number of temp hits: " << temphitCollection.size() << std::endl;
    std::cout << " EndOfEvent number of aux hits: " << hitCollection.size() << std::endl;
#endif
    std::sort(temphitCollection.begin(), temphitCollection.end());
    int geoId = -1;
    int trackId = -1;
    unsigned int counter = 0;
    std::unordered_set<unsigned int> setofIDs;

    for (auto it = temphitCollection.begin(); it != temphitCollection.end(); it++) {
#if defined _verbose_
      std::cout << "geoID: " << it->GetID() << " track ID: " << it->GetTrackID()
                << " Edep: " << it->GetEnergyDeposited() << " Parent Id: " << it->GetParentID()
                << " exit Time:  " << it->GetExitT() << " is first: " << it->IsIsfirstinVolume()
                << " is last:  " << it->IsIslastinVolume();
#endif
      if (it->GetID() == geoId && trackId == it->GetTrackID()) // trackid and detector didn't change
      {
#if defined _verbose_
        std::cout << " A" << std::endl;
#endif
        if (it->GetExitT()) // change exit vector and add to total charge
        {
          hitCollection[counter - 1].SetEnergyDeposited(
            hitCollection[counter - 1].GetEnergyDeposited() + it->GetEnergyDeposited());
          hitCollection[counter - 1].SetExitX(it->GetExitX());
          hitCollection[counter - 1].SetExitY(it->GetExitY());
          hitCollection[counter - 1].SetExitZ(it->GetExitZ());
          hitCollection[counter - 1].SetExitT(it->GetExitT());
        }
        else // just add to total charge
        {
          hitCollection[counter - 1].SetEnergyDeposited(
            hitCollection[counter - 1].GetEnergyDeposited() + it->GetEnergyDeposited());
        }
      }
      else if (setofIDs.find(it->GetParentID()) != setofIDs.end()) {
        setofIDs.insert(it->GetTrackID());
#if defined _verbose_
        std::cout << " A" << std::endl;
#endif
        hitCollection[counter - 1].SetEnergyDeposited(
          hitCollection[counter - 1].GetEnergyDeposited() + it->GetEnergyDeposited());
      }
      else if (it->GetID() != geoId) // new Detector
      {
        geoId = it->GetID();
        trackId = it->GetTrackID();
        setofIDs.clear();
        setofIDs.insert(it->GetTrackID());
#if defined _verbose_
        std::cout << " N" << std::endl;
#endif
        counter++;
        hitCollection.push_back(sim::AuxDetHit(it->GetID(),
                                               it->GetTrackID(),
                                               it->GetEnergyDeposited(),
                                               it->GetEntryX(),
                                               it->GetEntryY(),
                                               it->GetEntryZ(),
                                               it->GetEntryT(),
                                               it->GetExitX(),
                                               it->GetExitY(),
                                               it->GetExitZ(),
                                               it->GetExitT(),
                                               it->GetExitMomentumX(),
                                               it->GetExitMomentumY(),
                                               it->GetExitMomentumZ()));
      }
      else {
        trackId = it->GetTrackID();
#if defined _verbose_
        std::cout << " N" << std::endl;
#endif
        counter++;
        setofIDs.clear();
        setofIDs.insert(it->GetTrackID());
        hitCollection.push_back(sim::AuxDetHit(it->GetID(),
                                               it->GetTrackID(),
                                               it->GetEnergyDeposited(),
                                               it->GetEntryX(),
                                               it->GetEntryY(),
                                               it->GetEntryZ(),
                                               it->GetEntryT(),
                                               it->GetExitX(),
                                               it->GetExitY(),
                                               it->GetExitZ(),
                                               it->GetExitT(),
                                               it->GetExitMomentumX(),
                                               it->GetExitMomentumY(),
                                               it->GetExitMomentumZ()));
      }
    }
#if defined _verbose_
    std::cout << "Number of AuxDetHits: " << counter << std::endl;
#endif
  } // EndOfEvent

const sim::AuxDetHitCollection& larg4::AuxDetSD::GetHits ( ) const

inline

Definition at line 33 of file AuxDetSD.h.

References hitCollection.

Referenced by larg4::LArG4DetectorService::doFillEventWithArtHits().

{ return hitCollection; }

void larg4::AuxDetSD::Initialize

(

G4HCofThisEvent *

)

Definition at line 42 of file AuxDetSD.cc.

References hitCollection, and temphitCollection.

  {
    hitCollection.clear();
    temphitCollection.clear();
  }

G4bool larg4::AuxDetSD::ProcessHits	(	G4Step *	step,
		G4TouchableHistory *
	)

Definition at line 48 of file AuxDetSD.cc.

References edep, temphitCollection, and track.

  {
    G4double edep = step->GetTotalEnergyDeposit() / CLHEP::MeV;
    if (edep == 0.) return false;
    G4Track* track = step->GetTrack();
    const unsigned int trackID = track->GetTrackID();
    unsigned int ID = step->GetPreStepPoint()->GetPhysicalVolume()->GetCopyNo();
    TempHit tmpHit = TempHit(ID,
                             trackID,
                             track->GetParentID(),
                             step->IsFirstStepInVolume(),
                             step->IsLastStepInVolume(),
                             edep,
                             step->GetPreStepPoint()->GetPosition().getX() / CLHEP::cm,
                             step->GetPreStepPoint()->GetPosition().getY() / CLHEP::cm,
                             step->GetPreStepPoint()->GetPosition().getZ() / CLHEP::cm,
                             step->GetPreStepPoint()->GetGlobalTime() / CLHEP::ns,
                             step->GetPostStepPoint()->GetPosition().getX() / CLHEP::cm,
                             step->GetPostStepPoint()->GetPosition().getY() / CLHEP::cm,
                             step->GetPostStepPoint()->GetPosition().getZ() / CLHEP::cm,
                             step->GetPostStepPoint()->GetGlobalTime() / CLHEP::ns,
                             step->GetPostStepPoint()->GetMomentum().getX() / CLHEP::GeV,
                             step->GetPostStepPoint()->GetMomentum().getY() / CLHEP::GeV,
                             step->GetPostStepPoint()->GetMomentum().getZ() / CLHEP::GeV);
    temphitCollection.push_back(tmpHit);

    /*
  sim::AuxDetHit newHit = sim::AuxDetHit(ID,
                                    trackID,
                                    edep,
                                    step->GetPreStepPoint()->GetPosition().getX()/CLHEP::cm,
                                    step->GetPreStepPoint()->GetPosition().getY()/CLHEP::cm,
                                    step->GetPreStepPoint()->GetPosition().getZ()/CLHEP::cm,
                                    step->GetPreStepPoint()->GetGlobalTime()/CLHEP::ns,
                                    step->GetPostStepPoint()->GetPosition().getX()/CLHEP::cm,
                                    step->GetPostStepPoint()->GetPosition().getY()/CLHEP::cm,
                                    step->GetPostStepPoint()->GetPosition().getZ()/CLHEP::cm,
                                    step->GetPostStepPoint()->GetGlobalTime()/CLHEP::ns,
                                    step->GetPostStepPoint()->GetMomentum().getX()/CLHEP::GeV,
                                    step->GetPostStepPoint()->GetMomentum().getY()/CLHEP::GeV,
                                    step->GetPostStepPoint()->GetMomentum().getZ()/CLHEP::GeV
                                    );
  hitCollection.push_back(newHit);
  */
    return true;
  }

Member Data Documentation

sim::AuxDetHitCollection larg4::AuxDetSD::hitCollection

private

Definition at line 37 of file AuxDetSD.h.

Referenced by AuxDetSD(), EndOfEvent(), GetHits(), and Initialize().

TempHitCollection larg4::AuxDetSD::temphitCollection

private

Definition at line 36 of file AuxDetSD.h.

Referenced by EndOfEvent(), Initialize(), and ProcessHits().

---

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

• larg4/v09_19_03/source/larg4/Services/AuxDetSD.h
• larg4/v09_19_03/source/larg4/Services/AuxDetSD.cc