Changeset d1678fd in git for modules/SimpleCalorimeter.cc

Timestamp:

Mar 17, 2021, 5:07:53 PM (3 years ago)

Author:

michele <michele.selvaggi@…>

Branches:

master

Children:

5eda6767

Parents:

9cc5aeb

Message:

fix weighted time and tower position in calorimeter modules

File:

• modules/SimpleCalorimeter.cc (modified) (3 diffs)

modules/SimpleCalorimeter.cc

@@ -394 +394 @@
     fTowerEnergy += energy;
 
-    fTowerTime += energy * position.T();
-    fTowerTimeWeight += energy;
+    fTowerTime += energy * energy * position.T(); //sigma_t ~ 1/E
+    fTowerTimeWeight += energy * energy;
 
     fTower->AddCandidate(particle);
+    fTower->Position = position;
+
   }
 
@@ -409 +411 @@
 {
   Candidate *tower, *track, *mother;
-  Double_t energy, neutralEnergy, pt, eta, phi;
+  Double_t energy, neutralEnergy, pt, eta, phi, r;
   Double_t sigma, neutralSigma;
   Double_t time;
@@ -442 +444 @@
 
   pt = energy / TMath::CosH(eta);
-
-  fTower->Position.SetPtEtaPhiE(1.0, eta, phi, time);
+  r = TMath::Sqrt(fTower->Position.X()*fTower->Position.X()+fTower->Position.Y()*fTower->Position.Y());
+  
+  fTower->Position.SetPtEtaPhiE(r, eta, phi, time);
   fTower->Momentum.SetPtEtaPhiE(pt, eta, phi, energy);