use leading edge time for TOF, CTOF, PS and PSC instead of hit energy…

… weighted time
JeffersonLab · Mar 21, 2024 · a78ca55 · a78ca55
1 parent 2e3e84a
commit a78ca55
Show file tree

Hide file tree

Showing 4 changed files with 20 additions and 25 deletions.
diff --git a/src/GlueXSensitiveDetectorCTOF.cc b/src/GlueXSensitiveDetectorCTOF.cc
@@ -233,9 +233,11 @@ G4bool GlueXSensitiveDetectorCTOF::ProcessHits(G4Step* step,
          }
 
          if (merge_hit) {
-            // sum the charge, do energy weighting of the time
-            hiter->t_ns = (hiter->dE_GeV * hiter->t_ns + dEtop/GeV * ttop/ns) /
-	      (hiter->dE_GeV += dEtop/GeV);
+	    // Use the time from the earlier hit but add the charge
+	    hiter->dE_GeV += dEtop/GeV;
+	    if (hiter->t_ns*ns > ttop) {
+		    hiter->t_ns = ttop/ns;
+	    }
          }
          else {
             // create new hit 
@@ -264,10 +266,11 @@ G4bool GlueXSensitiveDetectorCTOF::ProcessHits(G4Step* step,
             }
          }
          if (merge_hit) {
-            // sum the charge, do energy weighting of the time
-            hiter->t_ns = (hiter->dE_GeV * hiter->t_ns +
-                           dEbottom/GeV * tbottom/ns) /
-            (hiter->dE_GeV += dEbottom/GeV);
+            // Use the time from the earlier hit but add the charge
+	    hiter->dE_GeV += dEbottom/GeV;
+	    if (hiter->t_ns*ns > tbottom) {
+		    hiter->t_ns = tbottom/ns;
+	    }
          }
          else {
             // create new hit 

diff --git a/src/GlueXSensitiveDetectorFTOF.cc b/src/GlueXSensitiveDetectorFTOF.cc
@@ -262,11 +262,6 @@ G4bool GlueXSensitiveDetectorFTOF::ProcessHits(G4Step* step,
          }
 
          if (merge_hit) {
-            // sum the charge, do energy weighting of the time
-            //hiter->t_ns = (hiter->dE_GeV * hiter->t_ns +
-            //               dEnorth/GeV * tnorth/ns) /
-            //(hiter->dE_GeV += dEnorth/GeV);
-
 	    // Use the time from the earlier hit but add the charge
 	    hiter->dE_GeV += dEnorth/GeV;
 	    if (hiter->t_ns*ns > tnorth) {
@@ -333,11 +328,6 @@ G4bool GlueXSensitiveDetectorFTOF::ProcessHits(G4Step* step,
             }
          }
          if (merge_hit) {
-            // sum the charge, do energy weighting of the time
-            //hiter->t_ns = (hiter->dE_GeV * hiter->t_ns +
-            //               dEsouth/GeV * tsouth/ns) /
-            //(hiter->dE_GeV += dEsouth/GeV);
-
 	    // Use the time from the earlier hit but add the charge
 	    hiter->dE_GeV += dEsouth/GeV;
 	    if (hiter->t_ns*ns > tsouth) {

diff --git a/src/GlueXSensitiveDetectorPS.cc b/src/GlueXSensitiveDetectorPS.cc
@@ -204,10 +204,11 @@ G4bool GlueXSensitiveDetectorPS::ProcessHits(G4Step* step,
          }
       }
       if (merge_hit) {
-         // Add the charge, do energy-weighted time averaging
-         hiter->t_ns = (hiter->t_ns * hiter->dE_GeV + t/ns * dEsum/GeV) /
-                       (hiter->dE_GeV + dEsum/GeV);
-         hiter->dE_GeV += dEsum/GeV;
+	 // Use the time from the earlier hit but add the charge
+	 hiter->dE_GeV += dEsum/GeV;
+	 if (hiter->t_ns*ns > t) {
+		 hiter->t_ns = t/ns;
+	 }
       }
       else {
          // create new hit 

diff --git a/src/GlueXSensitiveDetectorPSC.cc b/src/GlueXSensitiveDetectorPSC.cc
@@ -204,10 +204,11 @@ G4bool GlueXSensitiveDetectorPSC::ProcessHits(G4Step* step,
          }
       }
       if (merge_hit) {
-         // Add the charge, do energy-weighted time averaging
-         hiter->t_ns = (hiter->t_ns * hiter->dE_GeV + t/ns * dEsum/GeV) /
-                       (hiter->dE_GeV + dEsum/GeV);
-         hiter->dE_GeV += dEsum/GeV;
+	 // Use the time from the earlier hit but add the charge
+	 hiter->dE_GeV += dEsum/GeV;
+	 if (hiter->t_ns*ns > t) {
+		 hiter->t_ns = t/ns;
+	 }
       }
       else {
          // create new hit