Load pythia libraries.
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALCellInfo.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-2007, 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 /* History of cvs commits:
17  *
18 * $Log$
19 * Revision 1.2  2007/10/16 14:36:39  pavlinov
20 * fixed code violation (almost)
21 *
22 * Revision 1.1  2007/09/12 11:19:24  pavlinov
23 * added pi0 calibration, linearity, shower profile
24 *
25 */
26
27 // Aug 1, 2007 - cells information in one place
28 #include "AliEMCALCellInfo.h"
29 #include "AliEMCALGeometry.h"
30
31 #include <TObjArray.h>
32  
33 ClassImp(AliEMCALCellIndexes) 
34
35 // ------------------------------------------------------------------------------
36 AliEMCALCellIndexes::AliEMCALCellIndexes() : fAbsId(-1), fNSupMod(-1), fNModule(-1), fNIphi(-1), 
37 fNIeta(-1), fIPhi(-1), fIEta(-1), fIPhim(-1), fIEtam(-1)   
38 {
39 }
40
41
42 ClassImp(AliEMCALCellInfo)
43
44 // ------------------------------------------------------------------------------
45 AliEMCALCellInfo::AliEMCALCellInfo() : TNamed("",""), fTable(0), fCurrentInd(0)
46 {
47   //default constructor
48 }
49
50 // ------------------------------------------------------------------------------
51 AliEMCALCellInfo::AliEMCALCellInfo(const AliEMCALCellInfo& info) 
52   : TNamed(info.GetName(),info.GetTitle()), 
53     fTable(info.fTable), fCurrentInd(info.fCurrentInd)
54 {
55   //copy constructor
56 }
57
58 // ------------------------------------------------------------------------------
59 AliEMCALCellInfo::AliEMCALCellInfo(const char* name, const Int_t nrow) : TNamed(name,"table of cell information") , fTable(0), fCurrentInd(0)
60 {
61   fTable = new TObjArray(nrow);
62 }
63
64 // ------------------------------------------------------------------------------
65 void AliEMCALCellInfo::AddAt(AliEMCALCellIndexes* r)
66 {
67   (*fTable)[fCurrentInd] = new AliEMCALCellIndexes(*r);
68   fCurrentInd++;
69 }
70
71 // ------------------------------------------------------------------------------
72 AliEMCALCellInfo::~AliEMCALCellInfo()
73 {
74   if(fTable) {
75     fTable->Delete();
76     delete fTable;
77   }
78 }
79
80 // ------------------------------------------------------------------------------
81 AliEMCALCellIndexes* AliEMCALCellInfo::GetTable(Int_t i) const
82 {
83   // Oct 16, 2007
84   return (AliEMCALCellIndexes*)fTable->At(i);
85 }
86
87 // ------------------------------------------------------------------------------
88 void AliEMCALCellInfo::PrintTable(int ind1, int ind2) const
89 {
90   // Oct 16, 2007
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);
98     if(r==0) break;
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);
101   }
102 }
103
104 // ------------------------------------------------------------------------------
105 AliEMCALCellInfo *AliEMCALCellInfo::GetTableForGeometry(const char* geoName)
106 {
107   // Oct 16, 2007
108   if(geoName==0) return 0;
109   AliEMCALGeometry *g=AliEMCALGeometry::GetInstance(geoName);  
110   return GetTableForGeometry(g);
111 }
112
113 // ------------------------------------------------------------------------------
114 AliEMCALCellInfo *AliEMCALCellInfo::GetTableForGeometry(AliEMCALGeometry *g)
115 {
116   // Oct 16, 2007
117   if(g==0) return 0;
118   AliEMCALCellInfo *t = new AliEMCALCellInfo("CellInfo",g->GetNCells());
119
120   for(Int_t absid=0; absid<g->GetNCells(); absid++){
121     AliEMCALCellIndexes r;
122     r.fAbsId   = absid;
123
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); 
127
128     t->AddAt(&r);
129   }
130   return t;
131 }