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 **************************************************************************/
18 Revision 1.3 2007/09/12 17:44:22 pavlinov
19 fixed compilation problem under SuSe Linux
21 Revision 1.2 2007/09/11 19:38:15 pavlinov
22 added pi0 calibration, linearity, shower profile
26 //_________________________________________________________________________
27 // Cell folder which will keep all information
28 // about cell(tower) itself
29 // Initial version was created with TDataSet staf
30 // TObjectSet -> TFolder; Sep 6, 2007
32 //*-- Author: Aleksei Pavlinov (WSU, Detroit, USA)
34 #include "AliEMCALCell.h"
35 #include "AliEMCALHistoUtilities.h"
36 #include "AliEMCALGeometry.h"
37 #include "AliEMCALFolder.h"
38 #include "AliEMCALSuperModule.h"
39 #include "AliEMCALCalibData.h"
40 #include "AliEMCALRecPointsQaESDSelector.h"
41 #include "AliEMCALCalibCoefs.h"
52 typedef AliEMCALHistoUtilities u;
54 ClassImp(AliEMCALCell)
56 //______________________________________________________________
57 AliEMCALCell::AliEMCALCell() :
60 fAbsId(0),fSupMod(0),fModule(0),fPhi(0),fEta(0),fPhiCell(0),fEtaCell(0),fCcIn(0),fCcOut(0),
66 //______________________________________________________________
67 AliEMCALCell::AliEMCALCell(const AliEMCALCell& cell) :
68 TFolder(cell.GetName(),cell.GetTitle()),
69 fParent(cell.fParent),fLh(cell.fLh),
70 fAbsId(cell.fAbsId),fSupMod(cell.fSupMod),
71 fModule(cell.fModule),fPhi(cell.fPhi),
72 fEta(cell.fEta),fPhiCell(cell.fPhiCell),
73 fEtaCell(cell.fEtaCell),fCcIn(cell.fCcIn),
74 fCcOut(cell.fCcOut),fFun(cell.fFun)
79 //______________________________________________________________
80 AliEMCALCell::AliEMCALCell(const Int_t absId, const char* title) :
81 TFolder(Form("Cell%4.4i",absId),title),
83 fAbsId(absId),fSupMod(0),fModule(0),fPhi(0),fEta(0),fPhiCell(0),fEtaCell(0),fCcIn(0),fCcOut(0),
87 AliEMCALGeometry *g = AliEMCALGeometry::GetInstance();
88 g->GetCellIndex(fAbsId, fSupMod, fModule, fPhi, fEta);
89 g->GetCellPhiEtaIndexInSModule(fSupMod, fModule, fPhi, fEta, fPhiCell, fEtaCell);
93 //______________________________________________________________
94 AliEMCALCell::~AliEMCALCell()
99 //-------------------------------------------------------------------------------------
100 void AliEMCALCell::SetCCfromDB(AliEMCALCalibData *ccDb)
103 if(ccDb == 0) return;
104 // fADCchannelEC = fCalibData->GetADCchannel(iSupMod,ieta,iphi);
105 // fetaCel-column; fPhiCell- row
106 fCcIn = ccDb->GetADCchannel(fSupMod, fEtaCell, fPhiCell); // in GeV
108 TH1* h = (TH1*)GetHists()->At(0);
109 u::AddToNameAndTitle(h, 0, Form(", cc %5.2f MeV", fCcIn*1.e+3));
112 //______________________________________________________________
113 void AliEMCALCell::SetCCfromCCTable(AliEMCALCalibCoefs *t)
123 AliEMCALCalibCoef *r = t->GetTable(fAbsId);
124 if(r && r->fAbsId == fAbsId) {
126 } else { // something wrong
127 if(r) printf(" fAbsId %i : r->absId %i \n", fAbsId, r->fAbsId);
131 TH1* h = (TH1*)GetHists()->At(0);
132 u::AddToNameAndTitle(h, 0, Form(", cc %5.2f MeV", fCcIn*1.e+3));
135 //______________________________________________________________
136 void AliEMCALCell::FillEffMass(const Double_t mgg)
138 u::FillH1(GetHists(), 0, mgg);
141 //______________________________________________________________
142 void AliEMCALCell::FillCellNtuple(TNtuple *nt)
145 nt->Fill(fAbsId,fSupMod,fModule,fPhi,fEta,fPhiCell,fEtaCell,fCcIn,fCcOut);
148 //______________________________________________________________
149 void AliEMCALCell::FitHist(TH1* h, const char* name, const char* opt)
152 TString optFit(""), sopt(opt);
155 printf("<I> AliEMCALCell::FitHist : |%s| is started : opt %s\n", h->GetName(), opt);
156 TString tit(h->GetTitle());
158 TF1 *gausPol2 = 0, *g=0, *bg=0;
159 if(h->GetListOfFunctions()->GetSize() == 0 || 1) {
160 g = u::Gausi(name, 0.0, 0.4, h); // gaus estimation
162 g->SetParLimits(0, h->Integral()/20., h->Integral());
163 g->SetParLimits(1, 0.08, 0.20);
164 g->SetParLimits(2, 0.001, 0.02);
166 g->SetParameter(0, 1200.);
167 g->SetParameter(1, h->GetBinLowEdge(h->GetMaximumBin()));
168 g->SetParameter(2, 0.010);
170 g->SetLineColor(kRed);
174 h->Fit(g, optFit.Data(),"", 0.001, 0.3);
177 h->Fit(g, optFit.Data(),"", 0.001, 0.3);
179 bg = new TF1(Form("bg%s",name), "pol2", 0.0, 0.3);
181 h->Fit(bg, optFit.Data(),"", 0.0, 0.3);
183 gausPol2 = u::GausiPol2(name, 0.00, 0.3, g, bg);
185 h->Fit(gausPol2, optFit.Data(),"", 0.03, 0.28);
189 optFit = "0IME+"; // no drwaing at all
190 if(tit.Contains("SM") || sopt.Contains("DRAW")) optFit = "IME+";
192 gausPol2 = (TF1*)h->GetListOfFunctions()->At(0);
194 printf("<I> Function is defined alredy : %s optFit %s \n", gausPol2->GetTitle(), optFit.Data());
197 h->Fit(gausPol2, optFit.Data(),"", 0.01, 0.28);
199 if(optFit.Contains("0") == 0) {
200 gStyle->SetOptFit(111);
203 printf("<I> AliEMCALCell::FitHist : |%s| is ended \n\n", h->GetName());
206 //______________________________________________________________
207 void AliEMCALCell::FitEffMassHist(const char* opt)
210 static Double_t mPI0 = 0.13498; // mass of pi0
211 static Double_t mPI02 = mPI0*mPI0; // mass**2
213 AliEMCALFolder* emcal = AliEMCALRecPointsQaESDSelector::GetEmcalFolder();
214 Int_t it = emcal->GetIterationNumber();
216 TH1* h = (TH1*)GetHists()->At(0);
218 FitHist(h, GetName(), opt);
220 fFun = (TF1*)h->GetListOfFunctions()->At(0);
222 Double_t mpi = fFun->GetParameter(1), mpi2 = mpi*mpi;
223 Double_t ccTmp = fCcIn * mPI02 / mpi2;
224 if(it<=1) { // Jul 16, 2007
227 fCcOut = (ccTmp + fCcIn)/2.;
229 printf(" fFun %s | %s : iet %i\n", fFun->GetName(), fFun->GetTitle(), it);
231 printf(" %s | fCcIn %6.5f -> % 6.5f <- fCcOut \n", GetTitle(), fCcIn , fCcOut);
234 //______________________________________________________________
235 void AliEMCALCell::PrintInfo()
238 printf(" %s %s \n", GetName(), GetTitle());
239 if(fLh == 0 ) return;
240 TH1* h = (TH1*)GetHists()->At(0);
241 TF1 *f = (TF1*)h->GetListOfFunctions()->At(0);
242 if(fFun) printf(" fFun : %s | %s \n", fFun->GetName(), fFun->GetTitle());
247 //______________________________________________________________
248 TList* AliEMCALCell::BookHists()
252 TH1::AddDirectory(1);
254 AliEMCALFolder* emcal = AliEMCALRecPointsQaESDSelector::GetEmcalFolder();
255 Int_t it = emcal->GetIterationNumber();
257 new TH1F("01_EffMass", "effective mass of #gamma,#gamma(m_{#pi^{0}}=134.98 MeV) ", 60,0.0,0.3);
259 TList *l = AliEMCALHistoUtilities::MoveHistsToList(Form("HistsOfCell%4.4i",fAbsId), kFALSE);
260 AliEMCALHistoUtilities::AddToNameAndTitleToList(l, Form("4.4i_It%i",fAbsId,it),
261 Form(" Cell %4.4i, iter. %i",fAbsId, it));
263 TH1::AddDirectory(0);