From: hristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 13 May 2008 08:58:07 +0000 (+0000)
Subject: From v4-12-Release: updated beam size and PHOS simulation parameters
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=aad1175b5760f44ba8a92a3402b83a7ca8b8eb7e

From v4-12-Release: updated beam size and PHOS simulation parameters
---

diff --git a/test/fpprod/Config.C b/test/fpprod/Config.C
index ec4d2829215..f7e8dbc1b7e 100644
--- a/test/fpprod/Config.C
+++ b/test/fpprod/Config.C
@@ -40,6 +40,7 @@
 #include "FMD/AliFMDv1.h"
 #include "MUON/AliMUONv1.h"
 #include "PHOS/AliPHOSv1.h"
+#include "PHOS/AliPHOSSimParam.h"
 #include "PMD/AliPMDv1.h"
 #include "T0/AliT0v1.h"
 #include "EMCAL/AliEMCALv2.h"
@@ -199,7 +200,9 @@ void Config()
   //
   // Size of the interaction diamond
   // Longitudinal
-  Float_t sigmaz  = 7.55 / TMath::Sqrt(2.); // [cm]
+  Float_t sigmaz  = 5.4 / TMath::Sqrt(2.); // [cm]
+  if (energy == 900)
+    sigmaz  = 10.5 / TMath::Sqrt(2.); // [cm]
   //
   // Transverse
   Float_t betast  = 10;                 // beta* [m]
@@ -393,6 +396,23 @@ void Config()
         //=================== PHOS parameters ===========================
 
         AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
+        //Set simulation parameters different from the default ones.
+        AliPHOSSimParam* simEmc = AliPHOSSimParam::GetInstance() ;
+  
+        // APD noise of warm (+20C) PHOS:
+        // a2 = a1*(Y1/Y2)*(M1/M2), where a1 = 0.012 is APD noise at -25C,
+        // Y1 = 4.3 photo-electrons/MeV, Y2 = 1.7 p.e/MeV - light yields at -25C and +20C,
+        // M1 = 50, M2 = 50 - APD gain factors chosen for t1 = -25C and t2 = +20C,
+        // Y = MeanLightYield*APDEfficiency.
+
+        Float_t apdNoise = 0.012*2.5; 
+        simEmc->SetAPDNoise(apdNoise);
+
+        //Raw Light Yield at +20C
+        simEmc->SetMeanLightYield(18800);
+
+        //ADC channel width at +18C.
+        simEmc->SetADCchannelW(0.0125);
     }
 
 
@@ -495,6 +515,7 @@ void ProcessEnvironmentVars()
     // Energy
     if (gSystem->Getenv("CONFIG_ENERGY")) {
       energy = atoi(gSystem->Getenv("CONFIG_ENERGY"));
+      cout<<"Energy set to "<<energy<<" GeV"<<endl;
     }
 
     // Random Number seed