Print

Print


Commit in projects/lcdd/trunk/src/lcdd/detectors on MAIN
OpticalCalorimeterHitProcessor.cc+10-83036 -> 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 3036 -> 3037
--- 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) {
SVNspam 0.1


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