]>
Commit | Line | Data |
---|---|---|
8c823e5a | 1 | /************************************************************************** |
2 | * Copyright(c) 2007-2009, 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 | // ************************************************************* | |
17 | // Checks the quality assurance | |
18 | // by comparing with reference data | |
19 | // contained in a DB | |
20 | // ------------------------------------------------------------- | |
21 | // W. Ferrarese + P. Cerello Feb 2008 | |
22 | // INFN Torino | |
23 | ||
24 | // --- ROOT system --- | |
25 | #include <TH2D.h> | |
26 | #include <TTree.h> | |
27 | #include <TString.h> | |
28 | ||
29 | // --- Standard library --- | |
30 | ||
31 | // --- AliRoot header files --- | |
32 | #include "AliITSQASSDDataMakerRec.h" | |
33 | #include "AliLog.h" | |
34 | #include "AliQA.h" | |
35 | #include "AliQAChecker.h" | |
36 | #include "AliRawReader.h" | |
37 | #include "AliITSRawStreamSSD.h" | |
38 | #include "AliITSgeomTGeo.h" | |
39 | ||
40 | ||
41 | ClassImp(AliITSQASSDDataMakerRec) | |
42 | ||
43 | //____________________________________________________________________________ | |
44 | AliITSQASSDDataMakerRec::AliITSQASSDDataMakerRec(AliITSQADataMakerRec *aliITSQADataMakerRec, Bool_t kMode, Int_t ldc) : | |
45 | TObject(), | |
46 | fAliITSQADataMakerRec(aliITSQADataMakerRec), | |
47 | fkOnline(kMode), | |
48 | fLDC(ldc), | |
49 | fSSDhRaws(0), | |
50 | fSSDhRecs(0), | |
51 | fRawsOffset(0), | |
52 | fRecsOffset(0) | |
53 | { | |
54 | //ctor used to discriminate OnLine-Offline analysis | |
55 | //initilize the raw signal vs strip number histograms | |
56 | Int_t fLayer = 0,fLadder = 0, fModule = 0; | |
57 | TString fTitle; Int_t fHistCounter = 0; | |
58 | for(Int_t i = 500; i < fgkSSDMODULES + 500; i++) { | |
59 | AliITSgeomTGeo::GetModuleId(i,fLayer,fLadder,fModule); | |
60 | fTitle = "SSD_RawSignal_Layer"; fTitle += fLayer; | |
61 | fTitle += "_Ladder"; fTitle += fLadder; | |
62 | fTitle += "_Module"; fTitle += fModule; | |
63 | fHistSSDRawSignalModule[fHistCounter] = new TH1D(fTitle.Data(),fTitle.Data(),1540,0,1540); | |
64 | fHistCounter++; | |
65 | } | |
66 | } | |
67 | ||
68 | //____________________________________________________________________________ | |
69 | AliITSQASSDDataMakerRec::AliITSQASSDDataMakerRec(const AliITSQASSDDataMakerRec& qadm) : | |
70 | TObject(), | |
71 | fAliITSQADataMakerRec(qadm.fAliITSQADataMakerRec), | |
72 | fkOnline(qadm.fkOnline), | |
73 | fLDC(qadm.fLDC), | |
74 | fSSDhRaws(qadm.fSSDhRaws), | |
75 | fSSDhRecs(qadm.fSSDhRecs), | |
76 | fRawsOffset(qadm.fRawsOffset), | |
77 | fRecsOffset(qadm.fRecsOffset) | |
78 | { | |
79 | //copy ctor | |
80 | fAliITSQADataMakerRec->SetName((const char*)qadm.fAliITSQADataMakerRec->GetName()) ; | |
81 | fAliITSQADataMakerRec->SetTitle((const char*)qadm.fAliITSQADataMakerRec->GetTitle()); | |
82 | } | |
83 | ||
84 | //__________________________________________________________________ | |
85 | AliITSQASSDDataMakerRec& AliITSQASSDDataMakerRec::operator = (const AliITSQASSDDataMakerRec& qac ) | |
86 | { | |
87 | // Equal operator. | |
88 | this->~AliITSQASSDDataMakerRec(); | |
89 | new(this) AliITSQASSDDataMakerRec(qac); | |
90 | return *this; | |
91 | } | |
92 | ||
93 | //____________________________________________________________________________ | |
94 | void AliITSQASSDDataMakerRec::StartOfDetectorCycle() | |
95 | { | |
96 | //Detector specific actions at start of cycle | |
97 | AliDebug(1,"AliITSQADM::Start of SSD Cycle\n"); | |
98 | } | |
99 | ||
100 | //____________________________________________________________________________ | |
101 | void AliITSQASSDDataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX /*task*/, TObjArray* /*list*/) | |
102 | { | |
103 | // launch the QA checking | |
104 | AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n"); | |
105 | ||
106 | //SSD occupancy plots | |
107 | if(fkOnline) { | |
108 | Int_t fLayer = 0, fLadder = 0, fModule = 0; | |
109 | Int_t fOccPosition = 0; | |
110 | for(Int_t i = 0; i < fgkSSDMODULES; i++) { | |
111 | AliITSgeomTGeo::GetModuleId(i+500,fLayer,fLadder,fModule); | |
112 | fOccPosition = (fLayer == 5) ? fLadder : fLadder + 34; | |
113 | Double_t occupancy = GetSSDOccupancyRaws(fHistSSDRawSignalModule[i]); | |
114 | //if(occupancy > 0) { | |
115 | //AliInfo(Form("=========================================")); | |
116 | //AliInfo(Form("%s Position: %d - Initial histos: %d - Occupancy: %lf\n",fHistSSDRawSignalModule[i]->GetName(),fOccPosition,fInitialHistograms,occupancy)); | |
117 | //} | |
118 | fAliITSQADataMakerRec->GetRawsData(fRawsOffset + fOccPosition - 1)->Fill(fModule,occupancy); | |
119 | } | |
120 | }//online flag for SSD | |
121 | //AliQAChecker::Instance()->Run( AliQA::kITS , task, list); | |
122 | } | |
123 | ||
124 | //____________________________________________________________________________ | |
125 | void AliITSQASSDDataMakerRec::InitRaws() { | |
126 | ||
127 | // Initialization for RAW data - SSD - | |
128 | fRawsOffset = (fAliITSQADataMakerRec->fRawsQAList)->GetEntries(); | |
129 | ||
130 | const Int_t fSSDLADDERS = 72; | |
131 | if(fkOnline) { | |
132 | AliInfo("Book Online Histograms for SSD\n"); | |
133 | } | |
134 | else { | |
135 | AliInfo("Book Offline Histograms for SSD\n "); | |
136 | } | |
137 | AliInfo(Form("Number of histograms (SPD+SDD): %d\n",fRawsOffset)); | |
138 | TString fTitle = 0; | |
139 | //Int_t fLayer = 0,fLadder = 0, fModule = 0; | |
140 | //book online QA histos | |
141 | if(fkOnline) { | |
142 | TH1D *hOccupancyLadder[fSSDLADDERS]; | |
143 | for(Int_t iLayer = 5; iLayer < 7; iLayer++) { | |
144 | for(Int_t iLadder = 1; iLadder < AliITSgeomTGeo::GetNLadders(iLayer) + 1; iLadder++) { | |
145 | fTitle = "SSD_Occupancy_Layer"; fTitle += iLayer; | |
146 | fTitle += "_Ladder"; fTitle += iLadder; | |
147 | hOccupancyLadder[fSSDhRaws] = new TH1D(fTitle.Data(),fTitle.Data(),AliITSgeomTGeo::GetNDetectors(iLayer),0,AliITSgeomTGeo::GetNDetectors(iLayer)); | |
148 | hOccupancyLadder[fSSDhRaws]->GetXaxis()->SetTitleColor(1); | |
149 | hOccupancyLadder[fSSDhRaws]->GetXaxis()->SetTitle("Module number"); | |
150 | hOccupancyLadder[fSSDhRaws]->GetYaxis()->SetTitle("Occupancy [%]"); | |
151 | hOccupancyLadder[fSSDhRaws]->SetMarkerStyle(kFullCircle); | |
152 | fAliITSQADataMakerRec->Add2RawsList(hOccupancyLadder[fSSDhRaws], fRawsOffset + fSSDhRaws); | |
153 | fSSDhRaws++; | |
154 | }//ladder loop | |
155 | }//layer loop | |
156 | }//online flag | |
157 | AliDebug(1,Form("%d SSD Raws histograms booked\n",fSSDhRaws)); | |
158 | AliInfo(Form("Number of histograms (SPD+SDD+SSD): %d\n",fRawsOffset+fSSDhRaws)); | |
159 | ||
160 | ||
161 | } | |
162 | ||
163 | //____________________________________________________________________________ | |
164 | void AliITSQASSDDataMakerRec::MakeRaws(AliRawReader* rawReader) { | |
165 | // Fill QA for RAW - SSD - | |
166 | Int_t fStripNumber; | |
167 | Int_t fHistPosition; | |
168 | Int_t fLayer = 0,fLadder = 0, fModule = 0; | |
169 | ||
170 | rawReader->Select("ITSSSD",-1,-1); | |
171 | rawReader->Reset(); | |
172 | AliITSRawStreamSSD fSSDStream(rawReader); | |
173 | while (fSSDStream.Next()) { | |
174 | if(fSSDStream.GetModuleID() < 0) continue; | |
175 | AliITSgeomTGeo::GetModuleId(fSSDStream.GetModuleID(),fLayer,fLadder,fModule); | |
176 | fStripNumber = (fSSDStream.GetSideFlag() == 0) ? fSSDStream.GetStrip() : fSSDStream.GetStrip() + 768; | |
177 | fHistPosition = (fLayer == 5) ? ((fLadder - 1)*22 + fModule - 1) : ((fLadder - 1)*25 + fModule + 748 - 1); | |
178 | //AliInfo(Form("ModulePosition: %d - Layer: %d - Ladder: %d - Module: %d\n",fHistPosition,fLayer,fLadder,fModule)); | |
179 | fHistSSDRawSignalModule[fHistPosition]->Fill(fStripNumber,fSSDStream.GetSignal()); | |
180 | }//streamer loop | |
181 | } | |
182 | ||
183 | //____________________________________________________________________________ | |
184 | Double_t AliITSQASSDDataMakerRec::GetSSDOccupancyRaws(TH1 *lHisto) { | |
185 | // bo: TDC >0 or # of sigmas wrt noise ? | |
186 | Int_t lNumFiredBins = 0; | |
187 | for(Int_t iBin = 1; iBin < lHisto->GetNbinsX(); iBin++){ | |
188 | if (lHisto->GetBinContent(iBin)) | |
189 | lNumFiredBins++; | |
190 | } | |
191 | Double_t lOccupancy = (100.*lNumFiredBins)/lHisto->GetNbinsX(); // percentage | |
192 | //if(lOccupancy > 0) | |
193 | //AliInfo(Form("Fired strips: %d - Total strips: %d - Occupancy :%lf\n",lNumFiredBins,lHisto->GetNbinsX(),lOccupancy)); | |
194 | ||
195 | return lOccupancy; | |
196 | } | |
197 | ||
198 | ||
199 | //____________________________________________________________________________ | |
200 | void AliITSQASSDDataMakerRec::InitRecPoints() | |
201 | { | |
202 | // Initialization for RECPOINTS - SSD - | |
203 | fRecsOffset = (fAliITSQADataMakerRec->fRecPointsQAList)->GetEntries(); | |
204 | ||
205 | // custom code here | |
206 | //fSSDhRecs must be incremented by one unit every time a histogram is ADDED to the QA List | |
207 | ||
208 | AliDebug(1,Form("%d SSD Recs histograms booked\n",fSSDhRecs)); | |
209 | } | |
210 | ||
211 | //____________________________________________________________________________ | |
212 | void AliITSQASSDDataMakerRec::MakeRecPoints(TTree * /*clustersTree*/) | |
213 | { | |
214 | // Fill QA for recpoints - SSD - | |
215 | } | |
216 | ||
217 |