Changeset 9 in svn for trunk/src

Timestamp:

Nov 5, 2008, 7:10:54 PM (16 years ago)

Author:

Xavier Rouby

Message:

first commit

File:

• trunk/src/SmearUtil.cc (modified) (1 diff)

trunk/src/SmearUtil.cc

@@ -213 +213 @@
   
   if(fabs(hadron.Eta()) < MAX_CALO_CEN) {
-    energyS = gRandom->Gaus(energy_hcal, sqrt(
+    float energyS1 = gRandom->Gaus(energy_hcal, sqrt(
                                               pow(HAD_Nhcal,2) +
                                               pow(HAD_Chcal*energy_hcal,2) + 
-                                              pow(HAD_Shcal*sqrt(energy_hcal),2) )) 
-      +
-      gRandom->Gaus(energy_ecal, sqrt(
+                                              pow(HAD_Shcal*sqrt(energy_hcal),2) )) ;
+      
+    float energyS2 = gRandom->Gaus(energy_ecal, sqrt(
                                       pow(HAD_Necal,2) +
                                       pow(HAD_Cecal*energy_ecal,2) + 
-                                      pow(HAD_Secal*sqrt(energy_ecal),2) ) );     
+                                      pow(HAD_Secal*sqrt(energy_ecal),2) ) );
+
+        energyS = (energyS1>0)?energyS1:0 + (energyS2>0)?energyS2:0;
+        cout << "avant: " << energyS1 + energyS2 << "\t\t" << energyS << endl;
   } else { 
     energyS = gRandom->Gaus(energy,