EMCAL/AliEMCALTriggerBoard.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /*
  17
  18
  19 EMCal trigger board super class
  20 run the sliding window algorithm
  21 Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
  22 */
  23
  24 #include "AliEMCALTriggerBoard.h"
  25 #include "AliEMCALTriggerPatch.h"
  26 #include "AliLog.h"
  27
  28 #include <TClonesArray.h>
  29 #include <iostream>
  30 #include <cstdlib>
  31
  32 using namespace std;
  33
  34 ClassImp(AliEMCALTriggerBoard)
  35
  36 //_______________
  37 AliEMCALTriggerBoard::AliEMCALTriggerBoard() : TObject(),
  38 fRegion(0x0),
  39 fMap(0x0),
  40 fRegionSize(0x0),
  41 fSubRegionSize(0x0),
  42 fPatchSize(0x0),
  43 fPatches(0x0)
  44 {
  45
  46 }
  47
  48 //_______________
  49 AliEMCALTriggerBoard::AliEMCALTriggerBoard(const TVector2& RS) : TObject(),
  50 fRegion(0x0),
  51 fMap(0x0),
  52 fRegionSize(    new TVector2( RS ) ),
  53 fSubRegionSize( new TVector2() ),
  54 fPatchSize(     new TVector2() ),
  55 fPatches( new TClonesArray("AliEMCALTriggerPatch",10) )
  56 {
  57         // Ctor
  58
  59   fRegion = (int**)malloc( (int)fRegionSize->X() * sizeof( int* ) );
  60
  61   if (!fRegion) printf("Error: malloc could not allocate %d bytes for fRegion\n",
  62                        int(fRegionSize->X() * sizeof( int* )));
  63
  64   fMap = (int**)malloc( (int)fRegionSize->X() * sizeof( int* ) );
  65
  66   if (!fMap) printf("Error: malloc could not allocate %d bytes for fMap\n",
  67                     int(fRegionSize->X() * sizeof( int* )));
  68
  69   for (Int_t i=0;i<fRegionSize->X();i++)
  70   {
  71     if(fRegion){
  72       fRegion[i] = (int*)malloc( (int)fRegionSize->Y() * sizeof( int ) );
  73
  74       if (!fRegion[i]) printf("Error: malloc could not allocate %d bytes for fRegion[%d]\n",
  75                             i,int(fRegionSize->Y() * sizeof( int )));
  76     }
  77     if(fMap){
  78       fMap[i] = (int*)malloc( (int)fRegionSize->Y() * sizeof( int ) );
  79
  80       if (!fMap[i]) printf("Error: malloc could not allocate %d bytes for fMap[%d]\n",
  81                            i,int(fRegionSize->Y() * sizeof( int )));
  82     }
  83   }
  84
  85         // Initialize region matrix
  86         ZeroRegion();
  87         if(fMap){
  88         for (int i=0; i<fRegionSize->X(); ++i)
  89                 for (int j=0; j<fRegionSize->Y(); ++j) fMap[i][j] = 0;
  90   }
  91 }
  92
  93 //_______________
  94 AliEMCALTriggerBoard::~AliEMCALTriggerBoard()
  95 {
  96         // Dtor
  97
  98    for (Int_t i=0;i<fRegionSize->X();i++)
  99    {
 100       if (fRegion[i]) {free(fRegion[i]); fRegion[i] = 0;}
 101       if (   fMap[i]) {free(fMap[i]);    fMap[i] = 0;}
 102    }
 103
 104    free(fRegion); fRegion = 0x0;
 105    free(fMap);    fMap = 0x0;
 106
 107    if(fPatches)fPatches->Delete();
 108
 109    delete fPatches;
 110 }
 111
 112 //_______________
 113 void AliEMCALTriggerBoard::ZeroRegion()
 114 {
 115         // Initilize fRegion
 116
 117   if(fRegion){
 118     for (Int_t i=0;i<int(fRegionSize->X());i++) for (Int_t j=0;j<int(fRegionSize->Y());j++) fRegion[i][j] = 0;
 119   }
 120   else {
 121     AliFatal("fRegion was not previously initialized");
 122   }
 123
 124 }
 125
 126 //_______________
 127 void AliEMCALTriggerBoard::SlidingWindow(Int_t thres)
 128 {
 129         // Sliding window
 130         for (int i = 0; i <= int(fRegionSize->X() - fPatchSize->X() * fSubRegionSize->X()); i += int(fSubRegionSize->X())) {
 131                 for (int j = 0; j <= int(fRegionSize->Y() - fPatchSize->Y() * fSubRegionSize->Y()); j += int(fSubRegionSize->Y())) {
 132                         //
 133                         int sum = 0;
 134
 135                         for (int k = 0; k < int(fPatchSize->X() * fSubRegionSize->X()); k++) {
 136                                 for (int l = 0; l < int(fPatchSize->Y() * fSubRegionSize->Y()); l++) {
 137                                         //
 138                                         sum += fRegion[i + k][j + l];
 139                                 }
 140                         }
 141
 142                         if (sum > thres) {
 143                                 AliDebug(999, Form("Adding new patch at (%2d,%2d)", i, j));
 144                                 new((*fPatches)[fPatches->GetEntriesFast()]) AliEMCALTriggerPatch(i, j, sum);
 145                         }
 146                 }
 147         }
 148 }
 149
 150 //__________
 151 void AliEMCALTriggerBoard::Scan()
 152 {
 153         // Dump
 154
 155         cout << "     ";
 156         for (Int_t i=0; i<int(fRegionSize->X()); i++) printf("%8d ",i);
 157         cout << "\n";
 158         for (Int_t i=0; i<int(fRegionSize->X())-5; i++) printf("-------");
 159         cout << "\n";
 160
 161         for (Int_t i=0; i<int(fRegionSize->Y()); i++)
 162         {
 163                 if (i && !(i%12))
 164                 {
 165                         for (Int_t j=0; j<int(fRegionSize->X())-5; j++) printf("-------");
 166                         cout << endl;
 167                 }
 168
 169                 printf("%3d |",i);
 170                 for (Int_t j=0; j<int(fRegionSize->X()); j++)
 171                 {
 172                         printf("%2d/%5d ", fMap[j][i], fRegion[j][i]);
 173                 }
 174                 cout << endl;
 175         }
 176 }
 177
 178 //__________
 179 void AliEMCALTriggerBoard::Reset()
 180 {
 181         //
 182         fPatches->Delete();
 183         ZeroRegion();
 184 }
 185