Changeset 5418 for trunk/MagicSoft/Simulation/Detector/include-MLons

Timestamp:

11/17/04 11:38:28 (20 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Simulation/Detector/include-MLons

Files:

• MLons.cxx (modified) (6 diffs)
• MLons.hxx (modified) (3 diffs)

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

@@ -19 +19 @@
   fFwhmFadc = 2.0;
 
+  fGainFluctuations = 1;
+
   RandomNumber = new TRandom();
   RandomNumber -> SetSeed(0);
@@ -26 +28 @@
 MLons::MLons(Int_t in_shapeT, Float_t in_amplT, Float_t in_FwhmT, 
              Int_t in_shapeF, Float_t in_integF, Float_t in_FwhmF,
-             Float_t in_Fadc_Slices_per_ns)
+             Float_t in_Fadc_Slices_per_ns, Int_t in_Gain_Fluctuations)
 {
   //--------------------------------------------------------------------
@@ -43 +45 @@
   fFadcSlicesPerNanosec = in_Fadc_Slices_per_ns;
 
+  fGainFluctuations = in_Gain_Fluctuations;
+
   MSLStored = new MStarLight;
 
@@ -63 +67 @@
   fIntegFadc   = 0.; 
   fFwhmFadc    = 0.; 
+
+  fGainFluctuations = 1;
 
   MSLStored->Reset();
@@ -221 +227 @@
       strcat(filename_slt, ".slt");
  
-      // If brightness is different it check if the new file has the
+      // If brightness is different it checks if the new file has the
       // required parameters.
       // Note: I could be faster to store the whole trigger and fadc
@@ -236 +242 @@
           return 0;
         }
+
+      if ( !MSLStored->GetGainFluctuations() && fGainFluctuations)
+        {
+          cout << endl << "ERROR! The read NSB database in " << filename_slt << endl;
+          cout << "has been produced without PMT gain fluctuations. This is incompatible" << endl;
+          cout << "with the camera input card you are using. Please change the camera" << endl;
+          cout << "input card (adding the  noise_gain_fluctuations_off  command) or"  << endl;
+          cout << "produce a new NSB database (using program StarResponse) without" << endl;
+          cout << "gain fluctuations. Exiting..." << endl << endl;
+          return 0;
+        }
+
+      else if ( MSLStored->GetGainFluctuations() && !fGainFluctuations)
+        {
+          cout << endl << "ERROR! The read NSB database in " << filename_slt << endl;
+          cout << "has been produced with PMT gain fluctuations. This is incompatible" << endl;
+          cout << "with the camera input card command  noise_gain_fluctuations_off" << endl;
+          cout << "that you are using. Please change the camera input card or produce a" << endl;
+          cout << "new NSB database (using program StarResponse) without gain fluctuations." << endl;
+          cout << "Exiting..." << endl << endl;
+          return 0;
+        }
+
       MSLStored->SetBrightness(in_br);
 

trunk/MagicSoft/Simulation/Detector/include-MLons/MLons.hxx

@@ -57 +57 @@
   Float_t  fFadcSlicesPerNanosec; // The sampling frequency (GHz) of the FADC
 
+  Int_t    fGainFluctuations; // Indicates whether or not  gain fluctuations have
+                              // been simulated for producing the NSB database
+
+
   TRandom  *RandomNumber; // RandomGenerator 
 
@@ -65 +69 @@
   MLons(Int_t in_shapeT, Float_t in_amplT, Float_t in_FwhmT,
         Int_t in_shapeF, Float_t in_integF, Float_t in_FwhmF,
-        Float_t in_Fadc_Slices_per_ns);
+        Float_t in_Fadc_Slices_per_ns, Int_t in_Gain_Fluctuations);
 
   void Reset() ;
@@ -83 +87 @@
   Float_t GetFwhmFadc() ; 
   Float_t GetFadcSlicesPerNanosec() {return fFadcSlicesPerNanosec;}
+  Int_t   GetGainFluctuations() const { return fGainFluctuations; }
+
 
   void GetPath(Char_t *out);