Changeset f382b7e in git

Timestamp:

Sep 23, 2020, 1:44:18 PM (4 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

2ca681a

Parents:

55f821a

Message:

now filling neutral hadrons in DR towers

File:

• modules/DualReadoutCalorimeter.cc (modified) (1 diff)

modules/DualReadoutCalorimeter.cc

@@ -645 +645 @@
     //cout<<"Creating tower with Pt, Eta, Phi, Energy: "<<pt<<","<<eta<<","<<phi<<","<<neutralEnergy<<endl;
     tower->Momentum.SetPtEtaPhiE(pt, eta, phi, neutralEnergy);
-    tower->Eem = neutralEnergy;
-    tower->Ehad = 0.0;
-    tower->PID = 22;
-
+
+    // if no hadronic energy, use ECAL resolution 
+    if (fHCalTowerEnergy <= fHCalEnergyMin)
+    {
+      tower->Eem = neutralEnergy;
+      tower->Ehad = 0.0;
+      tower->PID = 22;
+    }
+
+    // if hadronic fraction > 0, use HCAL resolution 
+    else
+    {
+      tower->Eem = 0;
+      tower->Ehad = neutralEnergy;
+      tower->PID = 130;
+    }
+
+    fEFlowPhotonOutputArray->Add(tower);
     
-
-    fEFlowPhotonOutputArray->Add(tower);
-
 
     //clone tracks