Commit in projects/lcdd/trunk/src/lcdd/detectors on MAIN | |||
OpticalCalorimeterHitProcessor.cc | +10 | -8 | 3036 -> 3037 |
Remove debug print. Activate standard calorimeter hit processor outside of if statement so full edep appears in output. (Needs confirmation from Alex Conway.)
--- projects/lcdd/trunk/src/lcdd/detectors/OpticalCalorimeterHitProcessor.cc 2014-02-27 02:27:54 UTC (rev 3036) +++ projects/lcdd/trunk/src/lcdd/detectors/OpticalCalorimeterHitProcessor.cc 2014-03-03 19:20:27 UTC (rev 3037) @@ -17,12 +17,16 @@
} OpticalCalorimeterHitProcessor::~OpticalCalorimeterHitProcessor() {
- if (_cerenGenerator != 0) {
+ if (_cerenGenerator != NULL) {
delete _cerenGenerator; } } bool OpticalCalorimeterHitProcessor::processHits(G4Step* step) {
+ + /* Activate hit processing for standard edep collection. */ + LegacyCalorimeterHitProcessor::processHits(step); +
// FIXME: This should not be here. Unfortunately, it appears that instantiating // it in the constructor does not work due to unintialized Geant4 data structures. // The HitProcessor may need an initialize method called in start of run.
@@ -46,14 +50,12 @@
NCerenPhotons = 0; }
- if (NCerenPhotons <= 0) { - return LegacyCalorimeterHitProcessor::processHits(step); - } else {
+ /* Perform optical hit processing if Cerenkov photons were created by parameterization. */ + if (NCerenPhotons > 0) {
G4ThreeVector myPoint = step->GetPreStepPoint()->GetPosition();
- G4StepPoint* apreStepPoint = step->GetPreStepPoint();
G4double theEdep = double(NCerenPhotons); // get global cell pos from seg
- G4ThreeVector globalCellPos = getGlobalHitPosition(apreStepPoint);
+ G4ThreeVector globalCellPos = getGlobalHitPosition(pPreStepPoint);
// set the seg bins _calorimeter->getSegmentation()->setBins(step);
@@ -83,9 +85,9 @@
theEdep*GeV, theTrack->GetParticleDefinition()->GetAntiPDGEncoding(), theTrack->GetGlobalTime()));
- std::cout << "added contrib with energy: " << theEdep << std::endl; - return true;
} // end Cerenkov photon treatment
+ + return true;
} G4ThreeVector OpticalCalorimeterHitProcessor::getGlobalHitPosition(const G4StepPoint* aPreStepPoint) {
Use REPLY-ALL to reply to list
To unsubscribe from the LCDET-SVN list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=LCDET-SVN&A=1