1 #include "AliPHOSDA2.h"
6 //----------------------------------------------------------------
7 AliPHOSDA2::AliPHOSDA2(int module) : TNamed(),
8 fHistoFile(0),fFiredCells(0),fMod(module)
11 // Create AliPHOSDA2 ("Bad channels finder") object.
12 // module is the PHOS module number (0..4).
13 // Quality histogram names: module_iX_iZ_gain.
14 // Root file name: PHOS_ModuleX_BCM.root, where X - module number.
17 sprintf(name,"PHOS_Module%d_BCM",fMod);
20 SetTitle("Detector Algorithm to check for PHOS channels quality");
23 sprintf(rootname,"%s.root",GetName());
25 fHistoFile = new TFile(rootname,"recreate"); // new file!
27 for(Int_t iX=0; iX<64; iX++) {
28 for(Int_t iZ=0; iZ<56; iZ++) {
29 for(Int_t iGain=0; iGain<2; iGain++) {
30 fHQuality[iX][iZ][iGain] = 0;
38 fFiredCells = new TH1I("fFiredCells","Number of fired cells per event",100,0,1000);
42 //----------------------------------------------------------------
43 AliPHOSDA2::AliPHOSDA2(Int_t module, TObjArray* oldHistos) : TNamed(),
44 fHistoFile(0),fFiredCells(0),fMod(module)
47 // Create AliPHOSDA2 ("Bad channels finder") object.
48 // module is the PHOS module number (0..4).
49 // Quality histogram names: module_iX_iZ_gain.
50 // Read histograms from array oldHistos (if any).
51 // Do not produce an output file!
54 sprintf(name,"PHOS_Module%d_BCM",fMod);
57 SetTitle("Detector Algorithm to check for PHOS channels quality");
62 for(Int_t iX=0; iX<64; iX++) {
63 for(Int_t iZ=0; iZ<56; iZ++) {
64 for(Int_t iGain=0; iGain<2; iGain++) {
65 sprintf(hname,"%d_%d_%d_%d",fMod,iX,iZ,iGain);
67 hist1 = (TH1F*)oldHistos->FindObject(hname);
68 if(hist1) fHQuality[iX][iZ][iGain] = hist1;
70 fHQuality[iX][iZ][iGain] = 0;
78 fFiredCells = new TH1I("fFiredCells","Number of fired cells per event",100,0,1000);
82 //-------------------------------------------------------------------
83 AliPHOSDA2::AliPHOSDA2(const AliPHOSDA2& da) : TNamed(da),
84 fHistoFile(0),fFiredCells(0),fMod(da.fMod)
91 for(Int_t iX=0; iX<64; iX++) {
92 for(Int_t iZ=0; iZ<56; iZ++) {
93 for(Int_t iGain=0; iGain<2; iGain++) {
95 sprintf(hname,"%d_%d_%d_%d",fMod,iX,iZ,iGain);
96 hist1 = (TH1F*)da.fHistoFile->Get(hname);
97 if(hist1) fHQuality[iX][iZ][iGain] = new TH1F(*hist1);
99 fHQuality[iX][iZ][iGain] = 0;
105 fMaps[0] = new TH2F(*da.fMaps[0]);
110 fMaps[1] = new TH2F(*da.fMaps[1]);
114 fHistoFile = new TFile(da.GetName(),"recreate");
115 fFiredCells = new TH1I(*da.fFiredCells);
119 //-------------------------------------------------------------------
120 AliPHOSDA2& AliPHOSDA2::operator= (const AliPHOSDA2& da)
122 //Assignment operator.
126 TString oldname(fHistoFile->GetName());
127 TString newname(da.fHistoFile->GetName());
129 if(oldname != newname) {
131 fHistoFile = new TFile(da.fHistoFile->GetName(),"update");
136 SetName(da.GetName());
137 SetTitle(da.GetTitle());
139 for(Int_t iX=0; iX<64; iX++) {
140 for(Int_t iZ=0; iZ<56; iZ++) {
141 for(Int_t iGain=0; iGain<2; iGain++) {
142 if (fHQuality[iX][iZ][iGain]) delete fHQuality[iX][iZ][iGain];
143 fHQuality[iX][iZ][iGain] = da.fHQuality[iX][iZ][iGain];
150 fMaps[0] = da.fMaps[0];
155 fMaps[1] = da.fMaps[1];
160 fFiredCells = da.fFiredCells;
169 //-------------------------------------------------------------------
170 AliPHOSDA2::~AliPHOSDA2()
175 if(fHistoFile) delete fHistoFile;
179 //-------------------------------------------------------------------
180 void AliPHOSDA2::FillQualityHistograms(Float_t quality[64][56][2])
182 // Fills quality histograms.
183 // _By definition_, qood quality is 0<quality<1,
184 // all outside that is a bad quality.
185 // If no quality value read for particular channel,
186 // the correspondent array entry should be filled by zero.
187 // WARNING: this function should be called once per event!
192 for(Int_t iX=0; iX<64; iX++) {
193 for (Int_t iZ=0; iZ<56; iZ++) {
195 for(Int_t iGain=0; iGain<2; iGain++) {
196 if(!quality[iX][iZ][iGain]) continue;
198 if(fHQuality[iX][iZ][iGain])
199 fHQuality[iX][iZ][iGain]->Fill(quality[iX][iZ][iGain]);
201 sprintf(hname,"%d_%d_%d_%d",fMod,iX,iZ,iGain);
202 sprintf(htitl,"Quality for crystal %d_%d_%d and gain %d",fMod,iX,iZ,iGain);
203 fHQuality[iX][iZ][iGain] = new TH1F(hname,htitl,100,1.e-6,10.);
204 fHQuality[iX][iZ][iGain]->Fill(quality[iX][iZ][iGain]);
213 //-------------------------------------------------------------------
214 void AliPHOSDA2::FillFiredCellsHistogram(Int_t nCells)
216 fFiredCells->Fill(nCells);
219 //-------------------------------------------------------------------
220 void AliPHOSDA2::UpdateHistoFile()
222 // Write histograms to file
224 if(!fHistoFile) return;
225 if(!fHistoFile->IsOpen()) return;
232 sprintf(titl,"Quality map for Low gain");
233 fMaps[0] = new TH2F("gmaplow", titl, 64,0.,64.,56,0.,56.);
239 sprintf(titl,"Quality map for High gain");
240 fMaps[1] = new TH2F("gmaphigh", titl, 64,0.,64.,56,0.,56.);
245 for(Int_t iX=0; iX<64; iX++) {
246 for(Int_t iZ=0; iZ<56; iZ++) {
248 for(Int_t iGain=0; iGain<2; iGain++) {
249 hist1 = fHQuality[iX][iZ][iGain];
251 hist1->Write(hist1->GetName(),TObject::kWriteDelete);
252 Double_t mean = hist1->GetMean();
253 fMaps[iGain]->SetBinContent(iX+1,iZ+1,mean);
260 fMaps[0]->Write(fMaps[0]->GetName(),TObject::kWriteDelete);
261 fMaps[1]->Write(fMaps[1]->GetName(),TObject::kWriteDelete);
263 fFiredCells->Write();