1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
21 Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
24 #include "AliEMCALTriggerSTU.h"
25 #include "AliCDBManager.h"
26 #include "AliCDBEntry.h"
27 #include "AliEMCALCalibData.h"
28 #include "AliVZEROCalibData.h"
29 #include "AliVZEROdigit.h"
30 #include "AliEMCALTriggerPatch.h"
31 #include "AliESDVZERO.h"
34 #include <TClonesArray.h>
41 #include <Riostream.h>
44 ClassImp(AliEMCALTriggerSTU)
47 AliEMCALTriggerSTU::AliEMCALTriggerSTU() : AliEMCALTriggerBoard()
50 fV0M[0] = fV0M[1] = 0;
54 AliEMCALTriggerSTU::AliEMCALTriggerSTU(AliEMCALCalibData *calibData, const TVector2& RS) :
55 AliEMCALTriggerBoard(calibData, RS)
58 fV0M[0] = fV0M[1] = 0;
62 AliEMCALTriggerSTU::~AliEMCALTriggerSTU()
68 void AliEMCALTriggerSTU::BuildMap( Int_t iTRU, Int_t** M, const TVector2* rSize )
71 if ( iTRU == 31 ) iTRU = 35;
73 Int_t i2y = iTRU / 3 / 2;
75 if ( ( iTRU / 3 ) % 2 ) // odd (z<0) C side
77 Int_t i1y = 2 - ( iTRU - int( iTRU / 3 ) * 3 ); // 0 1 2 w/ increasing phi
79 for (Int_t i=0; i<rSize->X(); i++)
80 for (Int_t j=0; j<rSize->Y(); j++) fMap[24+i][j + i1y * 4 + i2y * 12] = M[i][j];
84 Int_t i1y = iTRU - int( iTRU / 3 ) * 3;
86 for (Int_t i=0; i<rSize->X(); i++)
87 for (Int_t j=0; j<rSize->Y(); j++) fMap[ i][j + i1y * 4 + i2y * 12] = M[i][j];
92 void AliEMCALTriggerSTU::L1( L1TriggerType_t type )
95 SlidingWindow( type, int(ThresholdFromV0( type )) );
99 void AliEMCALTriggerSTU::PrintADC( L1TriggerType_t type, TVector2& pos, TVector2& idx )
102 Int_t ix = (Int_t) (( pos.X() + fPatchSize->X() ) * fSubRegionSize->X());
103 Int_t iy = (Int_t) (( pos.Y() + fPatchSize->Y() ) * fSubRegionSize->Y());
105 TString subRegionADC[] = {"0->15", "16->31", "32->47", "48->63", "64->79", "80->95"};
111 Int_t iTRU = ( (ix-1) < 24 ) ? 31 - int(pos.Y() / 4) : 15 - int(pos.Y() / 4);
113 printf("TRU #%d row #%d col #%d fastor: ",iTRU,int(idx.Y()),int(idx.X()));
115 for (Int_t i=(Int_t)(pos.X() * fSubRegionSize->X());i<ix;i++)
117 for (Int_t j=(Int_t) (pos.Y() * fSubRegionSize->Y());j<iy;j++)
119 Int_t jtru = ( i < 24 ) ? 31 - j / 4 : 15 - j / 4;
120 printf("TRU#%d_%d ",jtru,fMap[i][j]);
129 //Int_t jTRU = ( (ix-1) < 24 ) ? 31 - (iy-1) / 4 : 15 - (iy-1) / 4;
131 printf("jet found at row : %d and col : %d",int(idx.X()),int(idx.Y()));
134 Int_t nSubRegion = 0;
136 for (Int_t i=(Int_t)(pos.X() * fSubRegionSize->X());i<ix;i++)
138 for (Int_t j=(Int_t)(pos.Y() * fSubRegionSize->Y());j<iy;j++)
140 Int_t itru = ( i < 24 ) ? 31 - j / 4 : 15 - j / 4;
142 Int_t idSR = fMap[i][j]/16;
144 Char_t value = ((itru << 3) & 0xF8) | (idSR & 0x7);
146 Bool_t isFound = kFALSE;
148 for (Int_t k=0;k<nSubRegion && nSubRegion;k++)
150 if (vPair[k] == value) isFound = kTRUE;
156 vPair = (Char_t*)realloc(vPair, nSubRegion * sizeof(Char_t));
157 if (vPair == NULL) {AliError("Error (re)allocating PrintADC() memory");}
159 vPair[nSubRegion-1] = value;
165 for (Int_t i=0;i<nSubRegion;i++)
167 cout << " TRU#" << ((vPair[i] & 0xF8) >> 3) << "_" << subRegionADC[(vPair[i] & 0x7)];
171 if (vPair) delete[] vPair;
175 AliError("AliEMCALTriggerSTU::PrintADC(): Undefined trigger type, pls check!");
180 void AliEMCALTriggerSTU::V0Multiplicity( TTree& treeV0 )
183 AliCDBManager *man = AliCDBManager::Instance();
184 AliCDBEntry *entry = man->Get("VZERO/Calib/Data");
185 AliVZEROCalibData *calibdata = (AliVZEROCalibData*)entry->GetObject();
187 TClonesArray* digitsArray = 0x0;
188 treeV0.SetBranchAddress("VZERODigit",&digitsArray);
193 for (Int_t i=0; i<64; i++)
199 Int_t nEntries = (Int_t)treeV0.GetEntries();
201 for (Int_t e=0; e<nEntries; e++)
205 Int_t nDigits = digitsArray->GetEntriesFast();
207 for (Int_t d=0; d<nDigits; d++)
209 AliVZEROdigit* digit = (AliVZEROdigit*)digitsArray->At(d);
210 Int_t pmNumber = digit->PMNumber();
212 if (adc[pmNumber] > (int(1.0/calibdata->GetMIPperADC(pmNumber)) /2) )
213 mult[pmNumber] += float(adc[pmNumber])*calibdata->GetMIPperADC(pmNumber);
219 for (Int_t j=0; j<32; j++)
221 fV0M[0] += short(mult[j ]+0.5);
222 fV0M[1] += short(mult[j+32]+0.5);
227 void AliEMCALTriggerSTU::FetchFOR( Int_t iTRU, Int_t **R, const TVector2* rSize )
229 // L1 triggers run over the whole EMCal surface
230 // STU builds its own fRegion aggregating TRU fRegion into one
232 // STU I from TRUs O in Olivier's coordinate system
234 if ( iTRU == 31 ) iTRU = 35;
236 Int_t i2y = iTRU / 3 / 2;
238 if ( ( iTRU / 3 ) % 2 ) // C side odd (z<0)
240 Int_t i1y = 2 - ( iTRU - int( iTRU / 3 ) * 3 ); // 0 1 2 w/ increasing phi
242 for (Int_t i=0; i<rSize->X(); i++) // 0:23 0:4
243 for (Int_t j=0; j<rSize->Y(); j++) fRegion[24+i][j + i1y * 4 + i2y * 12] = R[i][j];
247 Int_t i1y = iTRU - int( iTRU / 3 ) * 3;
249 for (Int_t i=0; i<rSize->X(); i++)
250 for (Int_t j=0; j<rSize->Y(); j++) fRegion[ i][j + i1y * 4 + i2y * 12] = R[i][j];
255 void AliEMCALTriggerSTU::PatchGenerator(const TClonesArray* lpos, Int_t val)
259 Int_t vTRU[32][24][4] ;//= {0};
261 for (Int_t i = 0; i < 32; i++)
262 for (Int_t j = 0; j < 24; j++)
263 for (Int_t k = 0; k < 4; k++) vTRU[i][j][k]=0;
266 AliWarning("AliEMCALTriggerSTU::PatchGenerator(): STU region has been reset!");
268 // Fill the patch FOR at 'pos' w/ value 'val'
270 TIter NextPosition( lpos );
272 while ( TVector2 *pos = (TVector2*)NextPosition() )
275 for (Int_t i=(Int_t)(pos->X()*fSubRegionSize->X()); i<int((pos->X()+fPatchSize->X())*fSubRegionSize->X()); i++)
276 for (Int_t j=(Int_t)(pos->Y()*fSubRegionSize->Y()); j<int((pos->Y()+fPatchSize->Y())*fSubRegionSize->Y()); j++)
280 for (Int_t i=0; i<2; i++)
282 for (Int_t j=0; j<16; j++)
284 Int_t iTRU = ( !i ) ? 31 - j : 15 - j;
286 for (Int_t i1=24*i; i1<24*(i+1); i1++)
288 for (Int_t j1=4*j; j1<4*(j+1); j1++)
290 vTRU[iTRU][i1%24][j1%4] = fRegion[i1][j1];
296 gSystem->Exec(Form("mkdir -p GP"));
298 for (Int_t i=0; i<32; i++)
300 std::ofstream outfile(Form("GP/data_TRU%d.txt",i),std::ios::trunc);
302 for (Int_t j=0;j<96;j++)
317 outfile << vTRU[i][ietam][iphim] << endl;
325 Float_t AliEMCALTriggerSTU::ThresholdFromV0( L1TriggerType_t type )
327 // Compute threshold FIXME: need an access to the OCDB
328 // to get f(V0) parameters depending on trigger type
338 AliError("AliEMCALTriggerSTU::ThresholdFromV0(): Undefined trigger type, pls check!");
345 void AliEMCALTriggerSTU::Reset()