Changeset 39022b6 in git for modules/Calorimeter.cc

Timestamp:

Jun 25, 2013, 2:43:27 PM (11 years ago)

Author:

pavel <pavel@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

be42bf4

Parents:

6085382

Message:

use log normal for calorimeter resolution

File:

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

modules/Calorimeter.cc

@@ -155 +155 @@
   fEFlowTrackOutputArray = ExportArray(GetString("EFlowTrackOutputArray", "eflowTracks"));
   fEFlowTowerOutputArray = ExportArray(GetString("EFlowTowerOutputArray", "eflowTowers"));
+
+  
+
 }
 
@@ -373 +376 @@
   Double_t energy, pt, eta, phi;
   Double_t ecalEnergy, hcalEnergy;
-
+  Double_t sE,mE,sH,mH,SE=0,ME=0,SH=0,MH=0;
+  
   if(!fTower) return;
 
-  ecalEnergy = gRandom->Gaus(fTowerECalEnergy, fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy));
+ // cout<<"new tower -----------------------"<<endl;
+
+  
+  sE = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy);
+  mE = fTowerECalEnergy;
+ 
+  sH = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy);
+  mH = fTowerHCalEnergy;
+  
+ // cout<<sE<<","<<mE<<endl;
+ // cout<<sH<<","<<mH<<endl;
+ // cout<<"---"<<endl;
+ 
+  if(mE>0)SE = TMath::Sqrt(TMath::Log((1+(sE*sE)/(mE*mE))));
+  if(mE>0)ME = TMath::Log(mE)-0.5*SE*SE;
+  
+  if(mH>0)SH = TMath::Sqrt(TMath::Log((1+(sH*sH)/(mH*mH))));
+  if(mH>0)MH = TMath::Log(mH)-0.5*SH*SH;
+  
+ // cout<<SE<<","<<ME<<endl;
+ // cout<<SH<<","<<MH<<endl;
+ // cout<<"---"<<endl;
+  
+//  ecalEnergy = gRandom->Gaus(fTowerECalEnergy, fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy));
+  ecalEnergy = (ME != 0.0 ? Draw_lognormal(ME,SE) : 0.0);
+ // cout<<ecalEnergy<<endl;
+  
   if(ecalEnergy < 0.0) ecalEnergy = 0.0;
 
-  hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy));
+  hcalEnergy = (MH != 0.0 ? Draw_lognormal(MH,SH) : 0.0);
+  //hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy));
   if(hcalEnergy < 0.0) hcalEnergy = 0.0;
 
+//  cout<<hcalEnergy<<endl; 
+   
   energy = ecalEnergy + hcalEnergy;
 
-  // eta = fTowerEta;
-  // phi = fTowerPhi;
+ // eta = fTowerEta;
+ // phi = fTowerPhi;
 
   eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]);
@@ -468 +501 @@
 
 //------------------------------------------------------------------------------
+
+Double_t Calorimeter::Draw_lognormal(Double_t mu, Double_t sigma)
+  {
+     Double_t g = gRandom->Gaus(0, 1);
+     return TMath::Exp(mu + sigma * g);
+  }