Changeset 6560 for trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx

Timestamp:

02/17/05 09:15:28 (20 years ago)

Author:

moralejo

Message:

Set as default option that of writing all event headers to output file,
not only those of the triggered events. To disable it, set the input card
flag "no_write_all_event_headers".

Changed such that output images for events below the minimum number of
photoelectrons nphe2NSB required to simulate the noise (NSB & electronic)
will be empty. This will avoid the problem of these events being processed,
without any noise, later in the chain. Although those images are not in the
output, one can still check in the headers (MMcTrig) how many such events
with less than nphe2NSB photoelectrons would have triggered.

File:

• trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx (modified) (9 diffs)

trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx

@@ -141 +141 @@
 //@: flag: TRUE: write the event headers of all events to output; 
 //@:       FALSE: only triggered showers
-static int Write_All_Event_Headers = FALSE;
+static int Write_All_Event_Headers = TRUE;
 
 static int Write_McEvt  = TRUE;
@@ -2011 +2011 @@
                             fadc_jitter[ict]);
 
-          inumphe=(inumphe<inumphe_CT[ict])?inumphe_CT[ict]:inumphe;
+          inumphe = (inumphe < inumphe_CT[ict])? inumphe_CT[ict] : inumphe;
+          // inumphe will be the maximum number of phe in any of the telescopes
+
 
           if( k != 0 ){ // non-zero return value means error
@@ -2021 +2023 @@
         // NSB simulation
 
-        if(simulateNSB && nphe2NSB<=inumphe)
+        if(simulateNSB && inumphe >= nphe2NSB)
           {
 
@@ -2083 +2085 @@
             }
           
-          }// end if(simulateNSB && nphe2NSB<=inumphe_CT[0]) ...
+          }// end if(simulateNSB && inumphe_CT[0] >= nphe2NSB) ...
 
 
@@ -2151 +2153 @@
         for(int ict=0;ict<ct_Number;ict++) 
           {
-            if (addElecNoise && nphe2NSB<=inumphe)
+            if (addElecNoise && inumphe >= nphe2NSB)
               { 
-                Trigger_CT[ict]->ElecNoise(Trigger_noise) ;
-                Fadc_CT[ict]->ElecNoise() ;
+                Trigger_CT[ict]->ElecNoise(Trigger_noise);
+                Fadc_CT[ict]->ElecNoise();
               }
           }
@@ -2187 +2189 @@
                     // Rise the discrimnator threshold to avoid huge rates
 
-                    if(riseDiskThres>0.0 && simulateNSB && nphe2NSB<=inumphe)
+                    if(riseDiskThres>0.0 && simulateNSB && inumphe >= nphe2NSB)
                       for(int ii=0;ii<ct_NPixels;ii++)
                         if( nsb_phepns_rotated[0][ii]>riseDiskThres)
@@ -2471 +2473 @@
 
               // Rise the discrimnator threshold to avoid huge rates
-              if(riseDiskThres>0.0 && simulateNSB && nphe2NSB<=inumphe)
+              if(riseDiskThres>0.0 && simulateNSB && inumphe >= nphe2NSB)
                 for(int ii=0;ii<ct_NPixels;ii++)
                   {
@@ -2564 +2566 @@
                 EvtHeader[ict]->SetTriggerPattern((UInt_t)MTriggerPattern::kTriggerLvl1);
 
-                //   fill pixel information
-              
-                if (Lev1MT[ict] /* || Write_All_Event_Headers */){
-                   if (addElecNoise) Fadc_CT[ict]->DigitalNoise();
-                  for(UInt_t i=0;
-                      i<((MGeomCam*)(camgeom.UncheckedAt(ict)))->GetNumPixels();
-                      i++){
+                // Fill pixel information
+                // AM 17/2/2005: added condition on inumphe. Noise is not generated for 
+                // events with less than phe2NSB photoelectrons, and then it is better not
+                // to write the images of those events to the output, so that we avoid them
+                // to be processed (with no noise) later in the chain.
+                
+                if (Lev1MT[ict] && inumphe >= nphe2NSB)
+                  {
+                    if (addElecNoise) 
+                      Fadc_CT[ict]->DigitalNoise();
+                    for(UInt_t i=0;
+                        i<((MGeomCam*)(camgeom.UncheckedAt(ict)))->GetNumPixels();
+                        i++){
 //
 // AM 15 01 2004 Commented out "continue" statement, so that also pixels
 // containing no C-photons will be written to the output in case of running
-// camera with no noise added to the signal.          
+// camera with no noise added to the signal.
 //                  if(!Fadc_CT[ict]->IsPixelUsed(i)) continue;
 //
-                    for (j = 0; j < FADC_slices_written; j++)
-                      {
-                        fadcValues->AddAt(Fadc_CT[ict]->GetFadcSignal(i,j),j);
-                        fadcValuesLow->AddAt(Fadc_CT[ict]->GetFadcLowGainSignal(i,j),j);
-                      }
-                    EvtData[ict]->AddPixel(i,fadcValues,0);
-                    EvtData[ict]->AddPixel(i,fadcValuesLow,kTRUE);
+                      for (j = 0; j < FADC_slices_written; j++)
+                        {
+                          fadcValues->AddAt(Fadc_CT[ict]->GetFadcSignal(i,j),j);
+                          fadcValuesLow->AddAt(Fadc_CT[ict]->GetFadcLowGainSignal(i,j),j);
+                        }
+                      EvtData[ict]->AddPixel(i,fadcValues,0);
+                      EvtData[ict]->AddPixel(i,fadcValuesLow,kTRUE);
+                    }
                   }
-                }
               } 
             }
@@ -5023 +5031 @@
 //
 // $Log: not supported by cvs2svn $
+// Revision 1.88  2005/02/11 20:00:01  moralejo
+//
+// Added to output container "MMcEvtBasic" with the most important MC
+// parameters, to be kept for all events (triggered or not) through
+// the whole analysis chain in order to allow the calculation of effective
+// areas.
+//
+// Updated version of MRawRunHeader from 4 to 5. This means the trigger and
+// calibration patterns are correctly set and can be decoded.
+//
 // Revision 1.87  2005/02/10 19:28:10  moralejo
 //