MollerDetectorSD Class Reference

#include <MollerDetectorSD.hh>

Public Member Functions
	MollerDetectorSD (G4String name)
virtual	~MollerDetectorSD ()
virtual void	Initialize (G4HCofThisEvent *HCE)
virtual G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist)
virtual void	EndOfEvent (G4HCofThisEvent *HCE)
Private Attributes
MollerDetectorHitsCollection *	hitsCollection
G4int	HCID

---

Detailed Description

Definition at line 14 of file MollerDetectorSD.hh.

---

Constructor & Destructor Documentation

MollerDetectorSD::MollerDetectorSD

(

G4String

name

)

Definition at line 16 of file MollerDetectorSD.cc.

References HCID.

:G4VSensitiveDetector(name)
{
  G4String HCname;
  collectionName.insert(HCname="hitsColl");
  HCID = -1;
}

MollerDetectorSD::~MollerDetectorSD

(

)

[virtual]

Definition at line 24 of file MollerDetectorSD.cc.

{;}

---

Member Function Documentation

void MollerDetectorSD::EndOfEvent

(

G4HCofThisEvent *

HCE

)

[virtual]

Definition at line 153 of file MollerDetectorSD.cc.

{
}

void MollerDetectorSD::Initialize

(

G4HCofThisEvent *

HCE

)

[virtual]

Definition at line 26 of file MollerDetectorSD.cc.

References HCID, and hitsCollection.

{
  hitsCollection = new MollerDetectorHitsCollection
                      (SensitiveDetectorName,collectionName[0]); 
  if(HCID<0)
  { HCID = G4SDManager::GetSDMpointer()->GetCollectionID(hitsCollection); }
  HCE->AddHitsCollection(HCID,hitsCollection);

}

G4bool MollerDetectorSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)			[virtual]

Definition at line 36 of file MollerDetectorSD.cc.

References hitsCollection, MollerDetectorHit::SetCreatorProcess(), MollerDetectorHit::SetIon(), MollerDetectorHit::SetKineticEnergy(), MollerDetectorHit::SetKineticEnergy2(), MollerDetectorHit::SetMomentum(), MollerDetectorHit::SetMomentum2(), MollerDetectorHit::SetParticleName(), MollerDetectorHit::SetScatAngle(), MollerDetectorHit::SetTrackID(), MollerDetectorHit::SetType(), MollerDetectorHit::SetVertexPos(), MollerDetectorHit::SetVolume(), and MollerDetectorHit::SetWorldPos().

{


  const G4Track* fTrack = aStep->GetTrack();

  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  G4TouchableHistory* theTouchable
    = (G4TouchableHistory*)(preStepPoint->GetTouchable());
  G4VPhysicalVolume* theMotherPhysical = theTouchable->GetVolume(1); // mother
  G4int copyNo = theMotherPhysical->GetCopyNo();


  G4ThreeVector worldPos = preStepPoint->GetPosition();
  G4ThreeVector vertexPos = fTrack->GetVertexPosition();
  G4ThreeVector momentum = fTrack->GetMomentumDirection();
  G4ThreeVector vertexMomentum = fTrack->GetVertexMomentumDirection();


  const G4VProcess* creatorProcess = fTrack->GetCreatorProcess();
  G4String process;
  if (creatorProcess!=0) process = creatorProcess->GetProcessName();

  G4bool ion = (process=="eIoni");

  G4String particle = fTrack->GetDefinition()->GetParticleName();

  // Need to get a number from the volume name somehow...
  //G4String volume = fTrack->GetVolume()->GetName();
  char volumeName[200];  char volumeNameAlt0[200];  char volumeNameAlt1[200];  char volumeNameAlt2[200];  char volumeNameAlt3[200];
  G4int offset=0;        G4int alt0offset=1000;     G4int alt1offset=1000;     G4int alt2offset=2000;     G4int alt3offset=2000;
  G4int volume=-1000;

  for (G4int vol = 0; vol<100; vol++){

    snprintf(volumeName,200,"planeDet_%i",vol); offset = 0;
    snprintf(volumeNameAlt0,200,"UScoll_%i",vol); alt0offset = 1050;
    snprintf(volumeNameAlt1,200,"DScoll_%i",vol); alt1offset = 1100;
    snprintf(volumeNameAlt2,200,"UScoil_%i",vol); alt2offset = 2000;
    snprintf(volumeNameAlt3,200,"DScoil_%i",vol); alt3offset = 2050;

    //Actual detectors will have volumes starting with 0-1000 (hopefully we won't need more than that?)
    //Collimators will have volumes starting with 1000
    //Coils will have volumes starting with 2000
    if (fTrack->GetVolume()->GetName().compare(volumeName)==0){
  volume = vol+offset;
    } else if (fTrack->GetVolume()->GetName().compare(volumeNameAlt0)==0){
  volume = vol+alt0offset;
    } else if (fTrack->GetVolume()->GetName().compare(volumeNameAlt1)==0){
  volume = vol+alt1offset;
    } else if (fTrack->GetVolume()->GetName().compare(volumeNameAlt2)==0){
  volume = vol+alt2offset;
    } else if (fTrack->GetVolume()->GetName().compare(volumeNameAlt3)==0){
  volume = vol+alt3offset;
    }
  }
  //  G4cout<<"The hit volume is : "<<fTrack->GetVolume()->GetName()<<" : Volume = "<<volume<<G4endl;

  G4int trackID = fTrack->GetTrackID();

  G4int partType;

  if (particle=="e-") partType=0;
  else if (particle=="e+") partType=1;
  else if (particle=="proton") partType=2;
  else if (particle=="antiproton") partType=3;
  else if (particle=="gamma") partType=4;
  else partType=5;

  G4double kineE = fTrack->GetKineticEnergy();
  G4double kineE0 = fTrack->GetVertexKineticEnergy();

  
  //G4double rad_dist = sqrt(worldPos.x()*worldPos.x() + 
  //worldPos.y()*worldPos.y());

  G4double scat_ang = sqrt(vertexMomentum.x()*vertexMomentum.x() + 
         vertexMomentum.y()*vertexMomentum.y());

  /* Insert cut here for limiting hit creation */
  /* This default cut ignores gamma particles, which have partType=4*/

  if (kineE0<100||!(partType==4||partType==0)) {
    return true;
  }

  /* Create and define Hit attributes */
  else {
    MollerDetectorHit* newHit = new MollerDetectorHit(copyNo);
    
    
    newHit->SetWorldPos(worldPos);
    newHit->SetVertexPos(vertexPos);
    newHit->SetMomentum(momentum);
    newHit->SetKineticEnergy(kineE);
    
    newHit->SetMomentum2(vertexMomentum);
    newHit->SetKineticEnergy2(kineE0);
    newHit->SetScatAngle(scat_ang);
    
    newHit->SetCreatorProcess(process);
    newHit->SetParticleName(particle);
    
    newHit->SetIon(ion);
    newHit->SetType(partType);
   
    //    newHit->SetVolume(atof(volume));
    newHit->SetVolume(volume);
    newHit->SetTrackID(trackID);
    
    
    hitsCollection->insert( newHit );
    
    return true;
    }
}

Here is the call graph for this function:

---

Field Documentation

G4int MollerDetectorSD::HCID [private]

Definition at line 27 of file MollerDetectorSD.hh.

Referenced by Initialize(), and MollerDetectorSD().

MollerDetectorHitsCollection* MollerDetectorSD::hitsCollection [private]

Definition at line 26 of file MollerDetectorSD.hh.

Referenced by Initialize(), and ProcessHits().

---

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

• MollerDetectorSD.hh
• MollerDetectorSD.cc