/scratch0/jsnuveri/BDSIM/BDSIMgit/bdsim/src/BDSCCDPixelSD.cc

#include "BDSGlobalConstants.hh" 
#include "BDSExecOptions.hh"
#include "BDSDebug.hh"
#include "BDSCCDPixelSD.hh"
#include "BDSCCDPixelHit.hh"
#include "G4VPhysicalVolume.hh"
#include "G4LogicalVolume.hh"
#include "G4Track.hh"
#include "G4Step.hh"
#include "G4StepPoint.hh"
#include "G4ParticleDefinition.hh"
#include "G4VTouchable.hh"
#include "G4TouchableHistory.hh"
#include "G4ios.hh"
#include "G4RotationMatrix.hh"
#include "G4ThreeVector.hh"

#include "G4AffineTransform.hh"

#include <vector>

#include "G4SDManager.hh"



BDSCCDPixelSD::BDSCCDPixelSD(G4String name)
  :G4VSensitiveDetector(name),CCDPixelCollection(NULL)
{
  _collectionName="CCDPixel";
  collectionName.insert(_collectionName);
}

BDSCCDPixelSD::~BDSCCDPixelSD()
{;}

void BDSCCDPixelSD::Initialize(G4HCofThisEvent*)
{
  // Create CCDPixel hits collection
  CCDPixelCollection = new BDSCCDPixelHitsCollection(SensitiveDetectorName,collectionName[0]);
}

G4bool BDSCCDPixelSD::ProcessHits(G4Step*aStep,G4TouchableHistory*)
{
#ifdef BDSDEBUG
  G4cout << __METHOD_NAME__ << "processing hits for sensitive detector name " << SensitiveDetectorName << G4endl;  
#endif
  G4TouchableHistory* hist = (G4TouchableHistory*) 
    (aStep->GetPreStepPoint()->GetTouchable());
  G4Track* theTrack = aStep->GetTrack();
  //only store a hit for optical photons.
  if(theTrack->GetDefinition()->GetPDGEncoding() != 0) return false;
  //Do not store hit if the particle is not on the boundary 
  if(aStep->GetPreStepPoint()->GetStepStatus()!=fGeomBoundary) return false;
  
  G4int nCCDPixel=hist->GetVolume()->GetCopyNo();
  G4double weight=theTrack->GetWeight();
  
#ifdef BDSDEBUG
  G4cout << __METHOD_NAME__ << " Storing hit: nCCDPixel weight" << G4endl;
  G4cout << __METHOD_NAME__ << " " << nCCDPixel <<" "  << weight << G4endl;
#endif

  BDSCCDPixelHit* pixelHit
    = new BDSCCDPixelHit(nCCDPixel, weight);
  CCDPixelCollection->insert(pixelHit);
  
#ifdef BDSDEBUG
  G4cout << __METHOD_NAME__ << " entries in hits collection after inserting hit: " << CCDPixelCollection->entries() << G4endl;
#endif
  //The hit was stored, so the return value is "true".
  return true;
}

void BDSCCDPixelSD::EndOfEvent(G4HCofThisEvent*HCE)
{
  G4SDManager * SDman = G4SDManager::GetSDMpointer();
  G4int HCID = SDman->GetCollectionID(collectionName[0]);
  HCE->AddHitsCollection(HCID, CCDPixelCollection );
}

---