Changeset 438 for trunk/MagicSoft/Simulation

Timestamp:

10/25/00 09:14:23 (24 years ago)

Author:

magicsol

Message:

The routine that produce poisson random numbers to decide how many phe
form NSB are emmited in each pixel has been replaced. Now a ROOT routine
is used.

File:

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

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

@@ -21 +21 @@
 //
 // $RCSfile: camera.cxx,v $
-// $Revision: 1.12 $
-// $Author: harald $ 
-// $Date: 2000-09-22 17:40:18 $
+// $Revision: 1.13 $
+// $Author: magicsol $ 
+// $Date: 2000-10-25 08:14:23 $
 //
 ////////////////////////////////////////////////////////////////////////
@@ -49 +49 @@
 #include "TROOT.h"
 
+#include "TRandom.h"
 #include "TApplication.h"
 
@@ -692 +693 @@
   if (Write_McTrig){
 
-    McTrig = new MMcTrig* [icontrigger];
+    McTrig = new MMcTrig * [icontrigger];
   
     for (i=0;i<icontrigger;i++) {
@@ -698 +699 @@
     }
 
-    HeaderTrig = new MHeaderTrig* [icontrigger];
+    HeaderTrig = new MHeaderTrig * [icontrigger];
   
     for (i=0;i<icontrigger;i++) {
@@ -740 +741 @@
     
     HeaderTree.Branch("MHeaderTrig","MHeaderTrig", 
-                      &HeaderTrig[0], bsize, split);    
-  }
-
+                 &HeaderTrig[0], bsize, split);    
+  }
   if (Trigger_Loop && Write_McTrig){
     for(char branchname[10],i=0;i<icontrigger;i++){
-      //
-      //   build the name of the branch for the different trigger conditions
-      //
+      
       sprintf(help,"%i",i+1);
       strcpy (branchname, "MHeaderTrig");
-      strcat (branchname,  &help[0]);
+      strcat (branchname, & help[0]);
       strcat (branchname, ".");
-      
       HeaderTree.Branch(branchname,"MHeaderTrig", 
                      &HeaderTrig[i], bsize, split);
@@ -770 +767 @@
 
   }
-  
-
   if(Trigger_Loop && Write_McTrig){
+
     for (int iconcount=0,ithrescount=0;ithrescount<=Trigger_loop_uthres-Trigger_loop_lthres;ithrescount++){
       for (imulticount=0;imulticount<=Trigger_loop_umult-Trigger_loop_lmult;imulticount++){
@@ -789 +785 @@
 
   //  Fill the Header Tree with the current leaves of each branch
-  
   HeaderTree.Fill() ;
-  
+            
+
   //      create a Tree for the Event data stream 
   TTree EvtTree("EvtTree","Events of Run");
 
   if (Write_McEvt){
+
     EvtTree.Branch("MMcEvt","MMcEvt", 
                    &McEvt, bsize, split);  
@@ -801 +798 @@
 
   if(!Trigger_Loop){
+    
     if (Write_RawEvt){
       EvtTree.Branch("MRawEvt","MRawEvt", 
@@ -810 +808 @@
     }    
   }
-  else{                     // trigger lopp 
+  else{
     if (Write_McTrig){
       for(char branchname[10],i=0;i<icontrigger;i++){
@@ -880 +878 @@
                           &cam,
                           nsbrate_phepns );
-
     if (k != 0){
       cout << "Error when reading starfield... \nExiting.\n";
@@ -892 +889 @@
         cam.dpixsizefactor[i] * cam.dpixsizefactor[i];
     }
+
   }
 
@@ -897 +895 @@
   // Read the reflector file with the Cherenkov data
   //                    
-  
+
   // select input file 
 
   if ( Data_From_STDIN ) {
+
     inputfile = stdin;
+
   }
   else{
@@ -909 +909 @@
     if ( inputfile == NULL ) 
       error( SIGNATURE, "Cannot open input file: %s\n", inname );
+
   }
   
@@ -984 +985 @@
 
         ++nshow;
-        if ( fmod ( nshow, 1000. ) == 0. ) 
-          log(SIGNATURE, "Event %d(+%d)\n", nshow, ntshow);
+        log(SIGNATURE, "Event %d(+%d)\n", nshow, ntshow);
         
         // get MCEventHeader
@@ -1037 +1037 @@
         
         mcevth.get_deviations ( &thetaCT, &phiCT );
-        
+
         if ( (thetaCT == 0.) && (phiCT == 0.) ) {
           
@@ -1134 +1134 @@
           exit(1);
         }
-        
-        if ( fmod ( nshow, 1000. ) == 0. )   
-          log(SIGNATURE, "End of this event: %d cphs(+%d). . .\n",
-              ncph, ntcph);
+          
+        log(SIGNATURE, "End of this event: %d cphs(+%d). . .\n",
+            ncph, ntcph);
 
         ntcph += ncph;
@@ -1159 +1158 @@
         }
         
-        if ( fmod ( nshow, 1000. ) == 0. ) 
-          log(SIGNATURE, "Total number of phe:  %d \n", inumphe ) ; 
+        cout << "Total number of phes: " << inumphe <<endl;
         
 
@@ -1175 +1173 @@
         // TRIGGER HERE
 
+
         //
         //   now the noise of the electronic 
@@ -1187 +1186 @@
         if(Trigger_Loop){
           //  Loop over trigger threshold
-          for (int iconcount=0,ithrescount=Trigger_loop_lthres;ithrescount<=Trigger_loop_uthres;ithrescount++) {
+          for (int iconcount=0,ithrescount=Trigger_loop_lthres;ithrescount<=Trigger_loop_uthres;ithrescount++){
             for (i=0;i<TRIGGER_PIXELS;i++)
               fpixelthres[i]=(float) ithrescount;
-
             Trigger.SetThreshold(fpixelthres);
+
             Trigger.Diskriminate();
-
             //
             //   look if in all the signals in the trigger signal branch
@@ -1410 +1408 @@
                            (UShort_t) ncph) ; 
             }
-
             //   We don not count photons out of the camera.    
-                    
+            
+            
             //
             //    write it out to the file outfile
             // 
             
-            //EvtTree.Fill() ; 
-            // huschel
+            EvtTree.Fill() ; 
+            
           
           
@@ -2973 +2971 @@
 
   float t;
+  TRandom random;             // Random numbers generator
+
+  random.SetSeed ((UInt_t) (RandomNumber*1000));
 
   ii = *inphe; // avoid dereferencing
@@ -3028 +3029 @@
 
       // randomize
-      
-      inumnsbphe =  ignpoi( inumnsbphe );
-      
+
+      if (inumnsbphe>0.0){
+        inumnsbphe = random.Poisson (inumnsbphe );
+      }
+
       // create the photoelectrons
       
@@ -3069 +3072 @@
 //
 // $Log: not supported by cvs2svn $
+// Revision 1.12  2000/09/22 17:40:18  harald
+// Added a lot of changes done by oscar.
+//
 // Revision 1.11  2000/09/21 11:47:33  harald
 // Oscar found some smaller errors in the calculation of the pixel shape and
 // corrected it.
 //
+// $Log: not supported by cvs2svn $
 // Revision 1.10  2000/07/04 14:10:20  MagicSol
 // Some changes have been done in the root output file. The RawEvt tree is only