Changeset d1678fd in git for modules/Calorimeter.cc

Timestamp:

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

Author:

michele <michele.selvaggi@…>

Branches:

master

Children:

5eda6767

Parents:

9cc5aeb

Message:

fix weighted time and tower position in calorimeter modules

File:

• modules/Calorimeter.cc (modified) (4 diffs)

modules/Calorimeter.cc

@@ -450 +450 @@
 
     fTower->AddCandidate(particle);
+    fTower->Position = position;
   }
 
@@ -461 +462 @@
 {
   Candidate *track, *tower, *mother;
-  Double_t energy, pt, eta, phi;
+  Double_t energy, pt, eta, phi, r;
   Double_t ecalEnergy, hcalEnergy;
   Double_t ecalNeutralEnergy, hcalNeutralEnergy;
@@ -511 +512 @@
   for(size_t i = 0; i < fTower->ECalEnergyTimePairs.size(); ++i)
   {
-    weight = TMath::Sqrt(fTower->ECalEnergyTimePairs[i].first);
+    weight = TMath::Power((fTower->ECalEnergyTimePairs[i].first),2);
     sumWeightedTime += weight * fTower->ECalEnergyTimePairs[i].second;
     sumWeight += weight;
@@ -517 +518 @@
   }
 
+  r = TMath::Sqrt(fTower->Position.X()*fTower->Position.X()+fTower->Position.Y()*fTower->Position.Y());
+
   if(sumWeight > 0.0)
   {
-    fTower->Position.SetPtEtaPhiE(1.0, eta, phi, sumWeightedTime / sumWeight);
+    fTower->Position.SetPtEtaPhiE(r, eta, phi, sumWeightedTime / sumWeight);
   }
   else
   {
-    fTower->Position.SetPtEtaPhiE(1.0, eta, phi, 999999.9);
+    fTower->Position.SetPtEtaPhiE(r, eta, phi, 999999.9);
   }