Ticket #363: Calorimeter_Dither.patch

.h

@@ -52 +52 @@
   Double_t fTrackECalWeightTime, fTrackHCalWeightTime;
   
   Int_t fTowerTrackHits, fTowerPhotonHits;
+  
+  Bool_t fDitherCenter;
 
   TFractionMap fFractionMap; //!
   TBinMap fBinMap; //!

.cc

@@ -152 +152 @@
   fEFlowTrackOutputArray = ExportArray(GetString("EFlowTrackOutputArray", "eflowTracks"));
   fEFlowPhotonOutputArray = ExportArray(GetString("EFlowPhotonOutputArray", "eflowPhotons"));
   fEFlowNeutralHadronOutputArray = ExportArray(GetString("EFlowNeutralHadronOutputArray", "eflowNeutralHadrons"));
-
-
+  
+  fDitherCenter = GetBool("DitherTowerCenters", true);
 }
 
 //------------------------------------------------------------------------------
@@ -433 +433 @@
   energy = ecalEnergy + hcalEnergy;
   time = (TMath::Sqrt(ecalEnergy)*ecalTime + TMath::Sqrt(hcalEnergy)*hcalTime)/(TMath::Sqrt(ecalEnergy) + TMath::Sqrt(hcalEnergy)); 
 
-//  eta = fTowerEta;
-//  phi = fTowerPhi;
-
-  eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]);
-  phi = gRandom->Uniform(fTowerEdges[2], fTowerEdges[3]);
+  if(fDitherCenter)
+  {
+     eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]);
+     phi = gRandom->Uniform(fTowerEdges[2], fTowerEdges[3]);
+  }
+  else
+  {
+     eta = fTowerEta;
+     phi = fTowerPhi;
+  }
 
   pt = energy / TMath::CosH(eta);