1 /**************************************************************************
2 * Copyright(c) 1998-2007, 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 **************************************************************************/
16 /* History of cvs commits:
19 * Revision 1.2 2007/10/16 14:36:39 pavlinov
20 * fixed code violation (almost)
22 * Revision 1.1 2007/09/12 11:19:24 pavlinov
23 * added pi0 calibration, linearity, shower profile
27 // Aug 1, 2007 - cells information in one place
28 #include "AliEMCALCellInfo.h"
29 #include "AliEMCALGeometry.h"
31 #include <TObjArray.h>
33 ClassImp(AliEMCALCellIndexes)
35 // ------------------------------------------------------------------------------
36 AliEMCALCellIndexes::AliEMCALCellIndexes() : fAbsId(-1), fNSupMod(-1), fNModule(-1), fNIphi(-1),
37 fNIeta(-1), fIPhi(-1), fIEta(-1), fIPhim(-1), fIEtam(-1)
42 ClassImp(AliEMCALCellInfo)
44 // ------------------------------------------------------------------------------
45 AliEMCALCellInfo::AliEMCALCellInfo() : TNamed("",""), fTable(0), fCurrentInd(0)
50 // ------------------------------------------------------------------------------
51 AliEMCALCellInfo::AliEMCALCellInfo(const AliEMCALCellInfo& info)
52 : TNamed(info.GetName(),info.GetTitle()),
53 fTable(info.fTable), fCurrentInd(info.fCurrentInd)
58 // ------------------------------------------------------------------------------
59 AliEMCALCellInfo::AliEMCALCellInfo(const char* name, const Int_t nrow) : TNamed(name,"table of cell information") , fTable(0), fCurrentInd(0)
61 fTable = new TObjArray(nrow);
64 // ------------------------------------------------------------------------------
65 void AliEMCALCellInfo::AddAt(AliEMCALCellIndexes* r)
67 (*fTable)[fCurrentInd] = new AliEMCALCellIndexes(*r);
71 // ------------------------------------------------------------------------------
72 AliEMCALCellInfo::~AliEMCALCellInfo()
80 // ------------------------------------------------------------------------------
81 AliEMCALCellIndexes* AliEMCALCellInfo::GetTable(Int_t i) const
84 return (AliEMCALCellIndexes*)fTable->At(i);
87 // ------------------------------------------------------------------------------
88 void AliEMCALCellInfo::PrintTable(int ind1, int ind2) const
91 printf(" %s : %s : #rows %i \n", GetName(), GetTitle(), fTable->GetSize());
92 if(ind1==-1 && ind2==-1) return;
93 printf(" fAbsId fNSupMod fNModule fNIphi fNIeta fIPhi fIEta fIPhim fIEtam\n");
94 if(ind1 < 0) ind1 = 0;
95 if(ind2 >= fTable->GetSize()) ind2 = fTable->GetSize();
96 for(int i=ind1; i<ind2; i++) {
97 AliEMCALCellIndexes* r = GetTable(i);
99 printf(" %5.5i %2.2i %3.3i %1.1i %1.1i %2.2i %2.2i %2.2i %2.2i\n",
100 r->fAbsId, r->fNSupMod, r->fNModule, r->fNIphi, r->fNIeta, r->fIPhi, r->fIEta, r->fIPhim,r->fIEtam);
104 // ------------------------------------------------------------------------------
105 AliEMCALCellInfo *AliEMCALCellInfo::GetTableForGeometry(const char* geoName)
108 if(geoName==0) return 0;
109 AliEMCALGeometry *g=AliEMCALGeometry::GetInstance(geoName);
110 return GetTableForGeometry(g);
113 // ------------------------------------------------------------------------------
114 AliEMCALCellInfo *AliEMCALCellInfo::GetTableForGeometry(AliEMCALGeometry *g)
118 AliEMCALCellInfo *t = new AliEMCALCellInfo("CellInfo",g->GetNCells());
120 for(Int_t absid=0; absid<g->GetNCells(); absid++){
121 AliEMCALCellIndexes r;
124 g->GetCellIndex(r.fAbsId, r.fNSupMod, r.fNModule, r.fNIphi, r.fNIeta);
125 g->GetCellPhiEtaIndexInSModule(r.fNSupMod,r.fNModule,r.fNIphi,r.fNIeta, r.fIPhi,r.fIEta);
126 g->GetModulePhiEtaIndexInSModule(r.fNSupMod,r.fNModule, r.fIPhim, r.fIEtam);