Changeset 5319 for trunk/MagicSoft/Simulation/Detector/include-MTrigger/MTrigger.cxx

Timestamp:

10/26/04 15:04:42 (20 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Simulation/Detector/include-MTrigger/MTrigger.cxx (modified) (10 diffs)

trunk/MagicSoft/Simulation/Detector/include-MTrigger/MTrigger.cxx

@@ -74 +74 @@
   //
 
-  fwhm_resp = RESPONSE_FWHM       ; 
-  ampl_resp = RESPONSE_AMPLITUDE  ; 
+  fwhm_resp = RESPONSE_FWHM; 
+  ampl_resp = RESPONSE_AMPLITUDE; 
 
   overlaping_time = TRIGGER_OVERLAPING;
 
-  threshold = CHANNEL_THRESHOLD  ; 
-  
-
-  gate_leng        = TRIGGER_GATE       ; 
-  trigger_multi    = TRIGGER_MULTI      ; 
-  trigger_geometry = TRIGGER_GEOM ; 
+  threshold = CHANNEL_THRESHOLD; 
+  
+
+  gate_leng        = TRIGGER_GATE; 
+  trigger_multi    = TRIGGER_MULTI; 
+  trigger_geometry = TRIGGER_GEOM; 
 
   //
@@ -234 +234 @@
  
 
-  //
-  //   the amplitude of one single photo electron is not a constant. 
-  //   There exists a measured distribution from Razmik. This distribution
-  //   is used to simulate the noise of the amplitude. 
-  //   For this a histogramm (histPmt) is created and filled with the
-  //   values. 
-  // 
-
-  histPmt = new TH1F ("histPmt","Noise of PMT", 40, 0., 40.) ;
-  
-  Stat_t ValRazmik[41] = { 0., 2.14, 2.06, 2.05, 2.05, 2.06, 2.07, 2.08,  2.15,
-                           2.27, 2.40, 2.48, 2.55, 2.50, 2.35, 2.20,  2.10,
-                           1.90, 1.65, 1.40, 1.25, 1.00, 0.80, 0.65,  0.50,
-                           0.35, 0.27, 0.20, 0.18, 0.16, 0.14, 0.12,  0.10, 
-                           0.08, 0.06, 0.04, 0.02, 0.01, 0.005,0.003, 0.001} ; 
-
-  histMean =  histPmt->GetMean() ;   
-  
-  for (i=0;i<41;i++){
-    histPmt->SetBinContent(i,ValRazmik[i]);
-  }
-
-  histMean =  histPmt->GetMean() ; 
+  // Initialize histogram histPmt which represents the distribution of amplitudes
+  // of the PMT response to a single photpelectron.
+
+  InitGainFluctuations();
+
 
   //
@@ -481 +463 @@
   peak_time = ( (Float_t) imax )  / ( (Float_t) TRIG_SLICES_PER_NSEC ) ;
  
-  //
-  //   the amplitude of one single photo electron is not a constant. 
-  //   There exists a measured distribution from Razmik. This distribution
-  //   is used to simulate the noise of the amplitude. 
-  //   For this a histogramm (histPmt) is created and filled with the
-  //   values. 
-  // 
-
+
+  // Initialize histogram histPmt which represents the distribution of amplitudes
+  // of the PMT response to a single photpelectron.
+
+  InitGainFluctuations();
+  
+  //
+  // Change name of histogram (tag with CT id) to avoid confusion in the case of
+  // multi-telescope simulations.
+  //
   char histname[32];
   sprintf(histname, "histPmt_%d", ct_id);
-  histPmt = new TH1F (histname,"Noise of PMT", 40, 0., 40.) ;
-  
-  Stat_t ValRazmik[41] = { 0., 2.14, 2.06, 2.05, 2.05, 2.06, 2.07, 2.08,  2.15,
-                           2.27, 2.40, 2.48, 2.55, 2.50, 2.35, 2.20,  2.10,
-                           1.90, 1.65, 1.40, 1.25, 1.00, 0.80, 0.65,  0.50,
-                           0.35, 0.27, 0.20, 0.18, 0.16, 0.14, 0.12,  0.10, 
-                           0.08, 0.06, 0.04, 0.02, 0.01, 0.005,0.003, 0.001} ; 
-
-  histMean =  histPmt->GetMean() ;   
-  
-  for (i=0;i<41;i++){
-    histPmt->SetBinContent(i,ValRazmik[i]);
-  }
-
-  histMean =  histPmt->GetMean() ; 
+  histPmt->SetName(histname);
+
 
   //
@@ -644 +615 @@
   //
 
-  fwhm_resp = fwhm       ; 
-  ampl_resp = ampl  ; 
+  fwhm_resp = fwhm; 
+  ampl_resp = ampl; 
 
   overlaping_time = overt;
 
     
-  threshold = CHANNEL_THRESHOLD  ; 
-  
-
-  gate_leng        = gate       ; 
-  trigger_multi    = TRIGGER_MULTI      ; 
-  trigger_geometry = TRIGGER_GEOM ; 
+  threshold = CHANNEL_THRESHOLD; 
+  
+
+  gate_leng        = gate; 
+  trigger_multi    = TRIGGER_MULTI; 
+  trigger_geometry = TRIGGER_GEOM; 
 
   cout << endl
@@ -710 +681 @@
   peak_time = ( (Float_t) imax )  / ( (Float_t) TRIG_SLICES_PER_NSEC ) ;
  
-  //
-  //   the amplitude of one single photo electron is not a constant. 
-  //   There exists a measured distribution from Razmik. This distribution
-  //   is used to simulate the noise of the amplitude. 
-  //   For this a histogramm (histPmt) is created and filled with the
-  //   values. 
-  // 
-
-  histPmt = new TH1F ("histPmt","Noise of PMT", 40, 0., 40.) ;
-  
-  Stat_t ValRazmik[41] = { 0., 2.14, 2.06, 2.05, 2.05, 2.06, 2.07, 2.08,  2.15,
-                           2.27, 2.40, 2.48, 2.55, 2.50, 2.35, 2.20,  2.10,
-                           1.90, 1.65, 1.40, 1.25, 1.00, 0.80, 0.65,  0.50,
-                           0.35, 0.27, 0.20, 0.18, 0.16, 0.14, 0.12,  0.10, 
-                           0.08, 0.06, 0.04, 0.02, 0.01, 0.005,0.003, 0.001} ; 
-
-  histMean =  histPmt->GetMean() ;   
-  
-  for (i=0;i<41;i++){
-    histPmt->SetBinContent(i,ValRazmik[i]);
-  }
-
-  histMean =  histPmt->GetMean() ; 
+
+  // Initialize histogram histPmt which represents the distribution of amplitudes
+  // of the PMT response to a single photpelectron.
+
+  InitGainFluctuations();
+
 
   //
@@ -931 +885 @@
     //  the FADC simulation
     //
-    PmtAmp =  (histPmt->GetRandom()/histMean) ; 
+    if (fGainFluctuations)
+      PmtAmp =  (histPmt->GetRandom()/histMean) ; 
+    else
+      PmtAmp = 1.0;
 
     // AM April 2004: removed updating of counters nphotshow, nphotnsb,
@@ -960 +917 @@
     
     //
-    // get the randomized amplitude
-    //
-    PmtAmp =  (histPmt->GetRandom()/histMean) ; 
+    // get the randomized amplitude (unless fGainFluctuations == kFALSE)
+    //
+    if (fGainFluctuations)
+      PmtAmp =  (histPmt->GetRandom()/histMean) ; 
+    else
+      PmtAmp = 1.0;
 
     //
@@ -1061 +1021 @@
   Float_t rausch ; 
 
-  rausch = RESPONSE_AMPLITUDE * factor ; 
+  rausch = ampl_resp * factor ; 
  
   cout<<"MTrigger::SetElecNoise ... generating database for electronic noise."
@@ -1080 +1040 @@
   //    to the signal
   //
-  Float_t rausch ; 
-
-  rausch = RESPONSE_AMPLITUDE * factor ; 
 
   UInt_t startslice;
@@ -1900 +1857 @@
 
 }
+
+//=================================================================
+
+void MTrigger::InitGainFluctuations()
+{
+  //
+  //   The amplitude of one single photo electron is not a constant. 
+  //   There exists a measured distribution from Razmik. This distribution
+  //   is used to simulate the noise of the amplitude. 
+  //   For this a histogramm (histPmt) is created and filled with the
+  //   values. 
+  // 
+
+  histPmt = new TH1F ("histPmt","Noise of PMT", 40, 0., 40.) ;
+  
+  Stat_t ValRazmik[40] = { 2.14, 2.06, 2.05, 2.05, 2.06, 2.07, 2.08,  2.15,
+                           2.27, 2.40, 2.48, 2.55, 2.50, 2.35, 2.20,  2.10,
+                           1.90, 1.65, 1.40, 1.25, 1.00, 0.80, 0.65,  0.50,
+                           0.35, 0.27, 0.20, 0.18, 0.16, 0.14, 0.12,  0.10, 
+                           0.08, 0.06, 0.04, 0.02, 0.01, 0.005,0.003, 0.001} ; 
+
+  for (Int_t i = 0; i < 40; i++)
+    histPmt->SetBinContent(i+1 ,ValRazmik[i]);
+
+  histMean =  histPmt->GetMean() ; 
+
+  fGainFluctuations = kTRUE; // Default value
+
+  return;
+}