MollerDetectorSD.cc

Go to the documentation of this file.

#include "MollerDetectorSD.hh"
#include "MollerDetectorHit.hh"
#include "G4Step.hh"
#include "G4HCofThisEvent.hh"
#include "G4TouchableHistory.hh"
#include "G4ios.hh"

#include "G4Track.hh"
#include "G4SDManager.hh"
#include "G4Navigator.hh"
#include "G4VProcess.hh"

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

MollerDetectorSD::~MollerDetectorSD(){;}

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

}

G4bool MollerDetectorSD::ProcessHits(G4Step* aStep, G4TouchableHistory*)
{


  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;
    }
}

void MollerDetectorSD::EndOfEvent(G4HCofThisEvent*)
{
}