EMCAL/AliEMCALTriggerSTU.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
  20
  21 Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
  22 */
  23
  24 #include "AliEMCALTriggerSTU.h"
  25 #include "AliCDBManager.h"
  26 #include "AliCDBEntry.h"
  27 #include "AliEMCALTriggerSTUDCSConfig.h"
  28 #include "AliVZEROCalibData.h"
  29 #include "AliVZEROdigit.h"
  30 #include "AliEMCALTriggerPatch.h"
  31 #include "AliESDVZERO.h"
  32 #include "AliLog.h"
  33
  34 #include <TClonesArray.h>
  35 #include <TSystem.h>
  36 #include <TH2F.h>
  37 #include <TFile.h>
  38 #include <TTree.h>
  39
  40 #include <fstream>
  41 #include <Riostream.h>
  42 #include <cstdlib>
  43
  44 ClassImp(AliEMCALTriggerSTU)
  45
  46 //_______________
  47 AliEMCALTriggerSTU::AliEMCALTriggerSTU() : AliEMCALTriggerBoard()
  48 ,fDCSConfig(0x0)
  49 {
  50         // Ctor
  51         fGammaTh[0] = fGammaTh[1] = 0;
  52         fJetTh[0] = fJetTh[1] = 0;
  53 }
  54
  55 //_______________
  56 AliEMCALTriggerSTU::AliEMCALTriggerSTU(AliEMCALTriggerSTUDCSConfig *dcsConf, const TVector2& RS) : AliEMCALTriggerBoard(RS)
  57 ,fDCSConfig(dcsConf)
  58 {
  59         // Ctor
  60         fGammaTh[0] = fGammaTh[1] = 0;
  61         fJetTh[0] = fJetTh[1] = 0;
  62 }
  63
  64 //_______________
  65 AliEMCALTriggerSTU::~AliEMCALTriggerSTU()
  66 {
  67         // Dtor
  68 }
  69
  70 //_______________
  71 void AliEMCALTriggerSTU::Build( TString& str, Int_t iTRU, Int_t** M, const TVector2* rSize )
  72 {
  73         // Build
  74
  75         str.ToLower();
  76
  77         Int_t ix = (iTRU % 2) ? 24 : 0;
  78
  79         Int_t iy =  iTRU / 2;
  80
  81         Int_t** v = 0x0;
  82
  83         if (str.Contains("map"))
  84         {
  85                 v = fMap;
  86         }
  87         else if (str.Contains("region"))
  88         {
  89                 v = fRegion;
  90         }
  91         else
  92         {
  93                 AliError("Operation not allowed: STU won't be configured properly!");
  94         }
  95
  96   if(v){
  97     for (Int_t i=0; i<rSize->X(); i++)
  98       for (Int_t j=0; j<rSize->Y(); j++) v[i + ix][j + iy * 4] = M[i][j];
  99   }
 100 }
 101
 102 //_______________
 103 void AliEMCALTriggerSTU::L1(int type)
 104 {
 105         // L1
 106
 107         TVector2 s1, s2, s3, s4;
 108         fDCSConfig->GetSegmentation(s1, s2, s3, s4);
 109
 110         switch (type)
 111         {
 112                 case kL1GammaHigh:
 113                 case kL1GammaLow:
 114                         SetSubRegionSize(s1);
 115                         SetPatchSize(s2);
 116                         break;
 117                 case kL1JetHigh:
 118                 case kL1JetLow:
 119                         SetSubRegionSize(s3);
 120                         SetPatchSize(s4);
 121                         break;
 122                 default:
 123                         AliError("Not supported L1 trigger type");
 124                         return;
 125                         break;
 126         }
 127
 128         SlidingWindow(GetThreshold(type));
 129         AliDebug(999, Form("STU type %d sliding window w/ thr %d found %d patches", type, GetThreshold(type), fPatches->GetEntriesFast()));
 130 }
 131
 132 //___________
 133 void AliEMCALTriggerSTU::ComputeThFromV0(int type, const Int_t M[])
 134 {
 135         // Compute threshold from V0
 136
 137         Short_t P[3] = {0};
 138
 139         switch (type)
 140         {
 141                 case kL1GammaHigh:
 142                         P[0] = fDCSConfig->GetG(0, 0);
 143                         P[1] = fDCSConfig->GetG(1, 0);
 144                         P[2] = fDCSConfig->GetG(2, 0);
 145                         break;
 146                 case kL1GammaLow:
 147                         P[0] = fDCSConfig->GetG(0, 1);
 148                         P[1] = fDCSConfig->GetG(1, 1);
 149                         P[2] = fDCSConfig->GetG(2, 1);
 150                         break;
 151                 case kL1JetHigh:
 152                         P[0] = fDCSConfig->GetJ(0, 0);
 153                         P[1] = fDCSConfig->GetJ(1, 0);
 154                         P[2] = fDCSConfig->GetJ(2, 0);
 155                         break;
 156                 case kL1JetLow:
 157                         P[0] = fDCSConfig->GetJ(0, 1);
 158                         P[1] = fDCSConfig->GetJ(1, 1);
 159                         P[2] = fDCSConfig->GetJ(2, 1);
 160                         break;
 161                 default:
 162                         AliError("AliEMCALTriggerSTU::ComputeThFromV0(): Undefined trigger type, pls check!");
 163                         return;
 164         }
 165
 166         ULong64_t v0sum = M[0] + M[1];
 167
 168         ULong64_t sqrV0 = v0sum * v0sum;
 169
 170         sqrV0 *= P[0];
 171
 172         sqrV0 >>= 32;
 173
 174         v0sum *= P[1];
 175
 176         v0sum >>= 16;
 177
 178         SetThreshold(type, (UShort_t)(sqrV0 + v0sum + P[2]));
 179 }
 180
 181 //___________
 182 void AliEMCALTriggerSTU::SetThreshold(int type, Int_t v)
 183 {
 184         // Set threshold
 185
 186         switch (type)
 187         {
 188                 case kL1GammaHigh:
 189                         fGammaTh[0] = v;
 190                         break;
 191                 case kL1GammaLow:
 192                         fGammaTh[1] = v;
 193                         break;
 194                 case kL1JetHigh:
 195                         fJetTh[0] = v;
 196                         break;
 197                 case kL1JetLow:
 198                         fJetTh[1] = v;
 199                         break;
 200                 default:
 201                         AliError("AliEMCALTriggerSTU::SetThreshold(): Undefined trigger type, pls check!");
 202         }
 203 }
 204
 205 //___________
 206 Int_t AliEMCALTriggerSTU::GetThreshold(int type)
 207 {
 208         // Compute threshold FIXME: need an access to the OCDB
 209         // to get f(V0) parameters depending on trigger type
 210
 211         switch (type)
 212         {
 213                 case kL1GammaHigh:
 214                         return fGammaTh[0];
 215                         break;
 216                 case kL1GammaLow:
 217                         return fGammaTh[1];
 218                         break;
 219                 case kL1JetHigh:
 220                         return fJetTh[0];
 221                         break;
 222                 case kL1JetLow:
 223                         return fJetTh[1];
 224                         break;
 225                 default:
 226                         AliError("AliEMCALTriggerSTU::GetThreshold(): Undefined trigger type, pls check!");
 227         }
 228
 229         return 0;
 230 }
 231
 232 //__________
 233 void AliEMCALTriggerSTU::Reset()
 234 {
 235         // Reset
 236
 237         fPatches->Delete();
 238 }