Changeset 2339 for trunk/MagicSoft/Simulation/Detector

Timestamp:

09/15/03 11:53:33 (21 years ago)

Author:

blanch

Message:

CAMERA_PIXELS changed to member variable numpix.
Less noise slices to avoid program abort.

File:

• trunk/MagicSoft/Simulation/Detector/include-MFadc/MFadc.cxx (modified) (21 diffs)

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

@@ -18 +18 @@
 #include "MGFadcSignal.hxx"
 
-MFadc::MFadc(Float_t ampl, Float_t fwhm, Float_t amplout, Float_t fwhmout) {
+MFadc::MFadc(Int_t pix, Float_t ampl, Float_t fwhm, Float_t amplout, Float_t fwhmout) {
   //
   //  Constructor overloaded II 
@@ -31 +31 @@
   //     this parameters of the trigger may be changed
   //  3. Then the all signals are set to zero
+  numpix=pix;
   
   fwhm_resp = fwhm; 
@@ -136 +137 @@
   //   if this is the first use, reset all signal for that pixels
   // 
-  if ( iPix > CAMERA_PIXELS ) {
+  if ( iPix > numpix ) {
     cout << " WARNING:  MFadc::Fill() :  iPix greater than CAMERA_PIXELS"
          << endl ;
@@ -201 +202 @@
   //   if this is the first use, reset all signal for that pixels
   // 
-  if ( iPix > CAMERA_PIXELS ) {
+  if ( iPix > numpix ) {
     cout << " WARNING:  MFadc::FillOuter() :  iPix greater than CAMERA_PIXELS"
          << endl ;
@@ -264 +265 @@
   //   if this is the first use, reset all signal for that pixels
   // 
-  if ( iPix > CAMERA_PIXELS ) {
+  if ( iPix > numpix ) {
     cout << " WARNING:  MFadc::Fill() :  iPix greater than CAMERA_PIXELS"
          << endl ;
@@ -298 +299 @@
   //   if this is the first use, reset all signal for that pixels
   // 
-  if ( iPix > CAMERA_PIXELS ) {
+  if ( iPix > numpix ) {
     cout << " WARNING:  MFadc::Fill() :  iPix greater than CAMERA_PIXELS"
          << endl ;
@@ -323 +324 @@
   Int_t i;
 
-  for(i=0;i<CAMERA_PIXELS;i++){
+  for(i=0;i<numpix;i++){
     pedestal[i]= (Float_t)(ped* GenElec->Rndm());
   }
 }
 
-void MFadc::SetPedestals(  Float_t ped[CAMERA_PIXELS]) { 
+void MFadc::SetPedestals(  Float_t *ped) { 
   //  It sets pedestal for each pixel from ped array
 
   Int_t i;
 
-  for(i=0;i<CAMERA_PIXELS;i++){
+  for(i=0;i<numpix;i++){
     pedestal[i]= ped[i];
   }
@@ -346 +347 @@
   Float_t baseline;
 
-  for(j=0;j<CAMERA_PIXELS;j++){
+  for(j=0;j<numpix;j++){
     baseline=0.0;
     for(i=0;i<(Int_t) SLICES_MFADC;i++){
@@ -367 +368 @@
   Int_t i, j;
 
-  for(i=0;i<CAMERA_PIXELS;i++)
+  for(i=0;i<numpix;i++)
     for(j=0;j<(Int_t)SLICES_MFADC;j++)
       sig[i][j]+=pedestal[i];
@@ -388 +389 @@
       // It does not exist, so all pixels will have the same offset
 
-      for(i=0;i<CAMERA_PIXELS;i++){
+      for(i=0;i<numpix;i++){
         if (used[i]){
           fdum=(10*GenOff->Rndm());
@@ -411 +412 @@
       // It does not exist, so all pixels will have the same offset
 
-      for(i=0;i<CAMERA_PIXELS;i++){
+      for(i=0;i<numpix;i++){
         if (used[i]){
           for(j=0;j<(Int_t) SLICES_MFADC;j++)
@@ -433 +434 @@
   UInt_t i;
  
-  cout<<"MFadc::SetElecNoise ... generating database for electroni noise."
+  cout<<"MFadc::SetElecNoise ... generating database for electronic noise."
       <<endl;
 
-  for (i=0;i<CAMERA_PIXELS*(Int_t) SLICES_MFADC*101;i++){
+  for (i=0;i<UInt_t(SLICES_MFADC*1001);i++){
     noise[i]=GenElec->Gaus(0., value  );
       }
@@ -452 +453 @@
   UInt_t startslice;
 
-  startslice=GenElec->Integer(CAMERA_PIXELS*(Int_t) SLICES_MFADC*100);
-
-    for ( Int_t i = 0 ; i < CAMERA_PIXELS; i++) {
+    for ( Int_t i = 0 ; i < numpix; i++) {
         //
         //  but at the beginning we must check if this pixel is
         //  hitted the first time
         //
+
+      startslice=GenElec->Integer((Int_t) SLICES_MFADC*1000);
+
         if ( used[i] == FALSE ) {
             used [i] = TRUE ; 
             
-            memcpy( (Float_t*)&sig[0][0],
-                    (Float_t*)&noise[startslice+i*(Int_t) SLICES_MFADC], 
+            memcpy( (Float_t*)&sig[i][0],
+                    (Float_t*)&noise[startslice], 
                     (Int_t) SLICES_MFADC*sizeof(Float_t));
-    
+      
         }
         //      
@@ -473 +475 @@
           for ( Int_t is=0 ; is< (Int_t)SLICES_MFADC ; is++ ) {
             
-            sig[i][is] += noise[startslice+i*(Int_t) SLICES_MFADC+is] ; 
+            sig[i][is] += noise[startslice+is] ; 
             
           }
-        
     }
 }
@@ -485 +486 @@
   
 
-  for ( Int_t ip=0; ip<CAMERA_PIXELS; ip++ ) {
+  for ( Int_t ip=0; ip<numpix; ip++ ) {
     
     if ( used[ip] == kTRUE ) {
@@ -531 +532 @@
   iFirstSlice = (Int_t) ( t /  WIDTH_FADC_TIMESLICE ) ; 
 
-  for ( Int_t ip=0; ip<CAMERA_PIXELS; ip++ ) {
+  for ( Int_t ip=0; ip<numpix; ip++ ) {
     
     if ( used[ip] == kTRUE ) {
@@ -563 +564 @@
   //    puts the pedestal values into the array offset
   
-  for ( Int_t i=0; i< CAMERA_PIXELS; i++ ) {
+  for ( Int_t i=0; i< numpix; i++ ) {
 
     offset[i] = pedestal[i] ; 
@@ -593 +594 @@
   iFirstSlice = (Int_t) ( t /  WIDTH_FADC_TIMESLICE ) ; 
 
-  for ( Int_t ip=0; ip<CAMERA_PIXELS; ip++ ) {
+  for ( Int_t ip=0; ip<numpix; ip++ ) {
     
     if ( used[ip] == kTRUE ) { 
@@ -600 +601 @@
         if (is< (Int_t)SLICES_MFADC && sig[ip][is]>0.0)
           {
-          output[ip][i++]=(sig[ip][is] > 255. ? 255 :(UChar_t) (sig[ip][is]+0.5));
-          output_lowgain[ip][i++]= 
-            (Int_t)(pedestal[ip]+(sig[ip][is]-pedestal[ip])/HIGH2LOWGAIN) > 255. ? 255 :
-            (UChar_t)(pedestal[ip]+(sig[ip][is]-pedestal[ip])/HIGH2LOWGAIN);
-          
+            output[ip][i]=(sig[ip][is] > 255. ? 255 :(UChar_t) (sig[ip][is]+0.5));
+            output_lowgain[ip][i]= 
+              (Int_t)(pedestal[ip]+(sig[ip][is]-pedestal[ip])/HIGH2LOWGAIN) > 255. ? 255 :
+              (UChar_t)(pedestal[ip]+(sig[ip][is]-pedestal[ip])/HIGH2LOWGAIN);
+            i++;
+          }
+        else if(sig[ip][is]>0.0)
+          {
+            output[ip][i]= (UChar_t)(pedestal[ip]+0.5);
+            output_lowgain[ip][i]= (UChar_t)(pedestal[ip]+0.5);
+            i++;
           }
         else 
           {
             output[ip][i]= 0;
-            output_lowgain[ip][i++]= 0;
+            output_lowgain[ip][i]= 0;
+            i++;
           }
       }
@@ -643 +651 @@
   
   Int_t ic = 0 ; 
-  for ( Int_t i=0 ; i < CAMERA_PIXELS; i++  ) {
+  for ( Int_t i=0 ; i < numpix; i++  ) {
     if ( used [i] == TRUE ) {
 
@@ -706 +714 @@
    return (output_lowgain[pixel][slice]);
 }
-
-
-