Changeset fc6300d in git for modules/SimpleCalorimeter.cc

Timestamp:

Jun 8, 2015, 2:22:11 PM (9 years ago)

Author:

Michele Selvaggi <selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

9c491e1

Parents:

e8070b6

Message:

added isEcal flag to simple calorimeter and fill Eem and Ehad variables accordingly (answer to #639)

File:

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

modules/SimpleCalorimeter.cc

@@ -149 +149 @@
 
   fEnergySignificanceMin = GetDouble("EnergySignificanceMin", 0.0);
+
+  // flag that says if current calo is Ecal of Hcal (will then fill correct values of Eem and Ehad)
+  fIsEcal = GetBool("IsEcal", false);
 
   // switch on or off the dithering of the center of calorimeter towers
@@ -425 +428 @@
   fTower->Momentum.SetPtEtaPhiE(pt, eta, phi, energy);
 
+  fTower->Eem = (!fIsEcal) ? 0 : energy;
+  fTower->Ehad = (fIsEcal) ? 0 : energy; 
+
   fTower->Edges[0] = fTowerEdges[0];
   fTower->Edges[1] = fTowerEdges[1];
@@ -447 +453 @@
     pt = energy / TMath::CosH(eta);
 
+    tower->Eem = (!fIsEcal) ? 0 : energy;
+    tower->Ehad = (fIsEcal) ? 0 : energy; 
+
     tower->Momentum.SetPtEtaPhiE(pt, eta, phi, energy);
     fEFlowTowerOutputArray->Add(tower);