* provided "as is" without express or implied warranty. *\r
**************************************************************************/\r
\r
-/* $Id: AliITSQADataMakerRec.cxx 23528 2008-01-24 23:46:36Z masera $ */\r
+/* $Id$ */\r
\r
// *************************************************************\r
// Checks the quality assurance \r
// INFN Torino\r
\r
// --- ROOT system ---\r
+#include <TProfile2D.h>\r
#include <TH2D.h>\r
#include <TBranch.h>\r
#include <TTree.h>\r
fAliITSQADataMakerRec(aliITSQADataMakerRec),\r
fkOnline(kMode),\r
fLDC(ldc),\r
-fSDDhRaws(0),\r
-fSDDhRecs(0),\r
-fRawsOffset(0),\r
-fRecsOffset(0),\r
-fSDDDDLModuleMap(0)\r
+fSDDhTask(0),\r
+fGenOffset(0),\r
+fTimeBinSize(1),\r
+fDDLModuleMap(0)\r
{\r
//ctor used to discriminate OnLine-Offline analysis\r
if(fLDC < 0 || fLDC > 4) {\r
AliError("Error: LDC number out of range; return\n");\r
}\r
- for(Int_t i=0;i<2*fgknSDDmodules;i++){\r
- fModuleChargeMap[i] = NULL;\r
- fmonoD[i] = NULL;\r
- }\r
+ //fDDLModuleMap=NULL;\r
}\r
\r
//____________________________________________________________________________ \r
fAliITSQADataMakerRec(qadm.fAliITSQADataMakerRec),\r
fkOnline(qadm.fkOnline),\r
fLDC(qadm.fLDC),\r
-fSDDhRaws(qadm.fSDDhRaws),\r
-fSDDhRecs(qadm.fSDDhRecs),\r
-fRawsOffset(qadm.fRawsOffset),\r
-fRecsOffset(qadm.fRecsOffset),\r
-fSDDDDLModuleMap(0)\r
+fSDDhTask(qadm.fSDDhTask),\r
+fGenOffset(qadm.fGenOffset),\r
+fTimeBinSize(1),\r
+fDDLModuleMap(0)\r
{\r
//copy ctor \r
fAliITSQADataMakerRec->SetName((const char*)qadm.fAliITSQADataMakerRec->GetName()) ; \r
fAliITSQADataMakerRec->SetTitle((const char*)qadm.fAliITSQADataMakerRec->GetTitle());\r
+ fDDLModuleMap=NULL;\r
}\r
\r
//____________________________________________________________________________ \r
AliITSQASDDDataMakerRec::~AliITSQASDDDataMakerRec(){\r
// destructor\r
- for(Int_t i=0;i<2*fgknSDDmodules;i++){\r
- if(fModuleChargeMap[i]) delete fModuleChargeMap[i];\r
- if(fmonoD[i]) delete fmonoD[i];\r
- }\r
- if(fSDDDDLModuleMap) delete fSDDDDLModuleMap; \r
+ // if(fDDLModuleMap) delete fDDLModuleMap;\r
}\r
//__________________________________________________________________\r
AliITSQASDDDataMakerRec& AliITSQASDDDataMakerRec::operator = (const AliITSQASDDDataMakerRec& qac )\r
}\r
\r
//____________________________________________________________________________ \r
-void AliITSQASDDDataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX /*task*/, TObjArray* /*list*/)\r
+void AliITSQASDDDataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t /*task*/, TObjArray* /*list*/)\r
{\r
// launch the QA checking\r
AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n"); \r
- \r
- //AliQAChecker::Instance()->Run( AliQA::kITS , task, list);\r
}\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerRec::InitRaws()\r
{ \r
// Initialization for RAW data - SDD -\r
- fRawsOffset = (fAliITSQADataMakerRec->fRawsQAList)->GetEntries();\r
+ fGenOffset = (fAliITSQADataMakerRec->fRawsQAList)->GetEntries();\r
\r
AliCDBEntry *ddlMapSDD = AliCDBManager::Instance()->Get("ITS/Calib/DDLMapSDD");\r
+ Bool_t cacheStatus = AliCDBManager::Instance()->GetCacheFlag();\r
+\r
if( !ddlMapSDD){\r
- AliFatal("Calibration object retrieval failed! ");\r
+ AliError("Calibration object retrieval failed! SDD will not be processed");\r
+ fDDLModuleMap = NULL;\r
+ return;\r
} \r
- AliITSDDLModuleMapSDD *ddlsdd=(AliITSDDLModuleMapSDD*)ddlMapSDD->GetObject();\r
+ fDDLModuleMap = (AliITSDDLModuleMapSDD*)ddlMapSDD->GetObject();\r
+ if(!cacheStatus)ddlMapSDD->SetObject(NULL);\r
ddlMapSDD->SetOwner(kTRUE);\r
- fSDDDDLModuleMap = ddlsdd;\r
+ if(!cacheStatus)delete ddlMapSDD;\r
\r
Int_t lay, lad, det;\r
Int_t LAY = -1; //, LAD = -1;\r
else {\r
AliInfo("Book Offline Histograms for SDD\n ");\r
}\r
- TH1F *h0 = new TH1F("ModPattern","Modules pattern",fgknSDDmodules,-0.5,259.5);\r
+ TH1D *h0 = new TH1D("ModPattern","HW Modules pattern",fgknSDDmodules,-0.5,259.5);\r
h0->GetXaxis()->SetTitle("Module Number");\r
h0->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h0,0+fRawsOffset);\r
- fSDDhRaws++;\r
+ fAliITSQADataMakerRec->Add2RawsList((new TH1D(*h0)),0+fGenOffset);\r
+ delete h0;\r
+ fSDDhTask++;\r
if(fLDC==0 || fLDC==1 || fLDC==2){\r
- TH1F *h1 = new TH1F("LadPatternL3","Ladder pattern L3",14,0.5,14.5); \r
+ TH1D *h1 = new TH1D("LadPatternL3","Ladder pattern L3",14,0.5,14.5); \r
h1->GetXaxis()->SetTitle("Ladder Number on Lay3");\r
h1->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h1,1+fRawsOffset);\r
- fSDDhRaws++;\r
+ fAliITSQADataMakerRec->Add2RawsList((new TH1D(*h1)),1+fGenOffset);\r
+ delete h1;\r
+ fSDDhTask++;\r
} \r
if(fLDC==0 || fLDC==3 || fLDC==4){\r
- TH1F *h2 = new TH1F("LadPatternL4","Ladder pattern L4",22,0.5,22.5); \r
+ TH1D *h2 = new TH1D("LadPatternL4","Ladder pattern L4",22,0.5,22.5); \r
h2->GetXaxis()->SetTitle("Ladder Number on Lay4");\r
h2->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h2,2+fRawsOffset);\r
- fSDDhRaws++;\r
+ fAliITSQADataMakerRec->Add2RawsList((new TH1D(*h2)),2+fGenOffset);\r
+ delete h2;\r
+ fSDDhTask++;\r
}\r
if(fLDC==0 || fLDC==1 || fLDC==2){\r
- TH1D *h3[fgkLADDonLAY3] ; \r
- for(Int_t i=1; i<=fgkLADDonLAY3; i++) {\r
+ for(Int_t i=1; i<=fgkLADDonLAY3; i++) {\r
sprintf(hname0,"ModPattern_L3_%d",i);\r
- h3[i-1] = new TH1D(hname0,hname0,6,0.5,6.5);\r
- h3[i-1]->GetXaxis()->SetTitle("Channel Number");\r
- h3[i-1]->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h3[i-1],i-1+3+fRawsOffset);\r
- fSDDhRaws++;\r
+ TH1D *h3 = new TH1D(hname0,hname0,6,0.5,6.5);\r
+ h3->GetXaxis()->SetTitle("Module Number");\r
+ h3->GetYaxis()->SetTitle("Counts");\r
+ fAliITSQADataMakerRec->Add2RawsList((new TH1D(*h3)),i-1+3+fGenOffset);\r
+ delete h3;\r
+ fSDDhTask++;\r
}\r
}\r
if(fLDC==0 || fLDC==3 || fLDC==4){\r
- TH1D *h4[fgkLADDonLAY4] ; \r
for(Int_t i=1; i<=fgkLADDonLAY4; i++) {\r
sprintf(hname0,"ModPattern_L4_%d",i);\r
- h4[i-1] = new TH1D(hname0,hname0,8,0.5,8.5);\r
- h4[i-1]->GetXaxis()->SetTitle("Channel Number");\r
- h4[i-1]->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h4[i-1],i-1+17+fRawsOffset);\r
- fSDDhRaws++;\r
+ TH1D *h4 = new TH1D(hname0,hname0,8,0.5,8.5);\r
+ h4->GetXaxis()->SetTitle("Module Number");\r
+ h4->GetYaxis()->SetTitle("Counts");\r
+ fAliITSQADataMakerRec->Add2RawsList((new TH1D(*h4)),i-1+17+fGenOffset);\r
+ delete h4;\r
+ fSDDhTask++;\r
}\r
}\r
\r
Int_t indexlast2 = 0;\r
\r
if(fkOnline) {\r
+ fTimeBinSize = 4;\r
indexlast = 0;\r
- TH1D *h5[fgknSDDmodules*2] ; \r
index1 = 0;\r
-\r
-/*\r
- if(fLDC == 0 || fLDC == 1 || fLDC == 4){\r
- \r
- for(Int_t lad=0; lad<fgkLADDonLAY3; lad++){\r
- for(Int_t ichan=0; det<3; det++){\r
- for(Int_t iside=0; iside<fgknSide; iside++){\r
- sprintf(hname0,"ProjYMap_L%d_%d_%d_%d",lay,lad,det+sideflag*3,iside); //need to define "lay"\r
- h5[index1] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
- h5[index1]->GetXaxis()->SetTitle("Anode Number");\r
- h5[index1]->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h5[index1],index1+39+fRawsOffset);\r
- fSDDhRaws++;\r
- index1++;\r
- indexlast = index1+39;\r
- }\r
- } \r
- }\r
- }\r
-*/ \r
- for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
- for(Int_t iside=0; iside<fgknSide; iside++){\r
- AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
- //if(fLDC == 0 || ((fLDC == 1 || fLDC == 2) && lay == 2) || ((fLDC == 3 || fLDC == 4) && lay == 3)\r
- sprintf(hname0,"ProjYMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
- h5[index1] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
- h5[index1]->GetXaxis()->SetTitle("Anode Number");\r
- h5[index1]->GetYaxis()->SetTitle("Counts");\r
- fAliITSQADataMakerRec->Add2RawsList(h5[index1],index1+39+fRawsOffset);\r
- fSDDhRaws++;\r
- index1++;\r
- indexlast = index1+39;\r
- }\r
- }\r
- }\r
- \r
- TH1D *h6[fgknSDDmodules*8] ; \r
- Int_t indextot = 0;\r
- indexlast1 = 0;\r
- for(Int_t htype=0; htype<4; htype++){\r
- for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) { \r
- for(Int_t iside=0; iside<fgknSide; iside++){\r
- AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
- if(htype == 0) sprintf(hname0,"CountsVSAnode_L%d_%d_%d_%d",lay,lad,det,iside);\r
- if(htype == 1) sprintf(hname0,"ChargeVSAnode_L%d_%d_%d_%d",lay,lad,det,iside);\r
- if(htype == 2) sprintf(hname0,"CountsVSTbin_L%d_%d_%d_%d",lay,lad,det,iside);\r
- if(htype == 3) sprintf(hname0,"ChargeVSTbin_L%d_%d_%d_%d",lay,lad,det,iside);\r
- h6[indextot] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
- if(htype == 0){\r
- h6[indextot]->GetXaxis()->SetTitle("Anode Number");\r
- h6[indextot]->GetYaxis()->SetTitle("Counts");\r
- }\r
- if(htype == 1){\r
- h6[indextot]->GetXaxis()->SetTitle("Anode Number");\r
- h6[indextot]->GetYaxis()->SetTitle("Charge (ADC)");\r
- }\r
- if(htype == 2){\r
- h6[indextot]->GetXaxis()->SetTitle("Time bin");\r
- h6[indextot]->GetYaxis()->SetTitle("Counts");\r
- }\r
- if(htype == 3){\r
- h6[indextot]->GetXaxis()->SetTitle("Time bin");\r
- h6[indextot]->GetYaxis()->SetTitle("Charge (ADC)");\r
- }\r
- fAliITSQADataMakerRec->Add2RawsList(h6[indextot],indexlast + indextot + fRawsOffset);\r
- fSDDhRaws++;\r
- indextot++;\r
- indexlast1 = indexlast + indextot;\r
- }\r
- }\r
- }\r
- }\r
-\r
-\r
+ indexlast1 = fSDDhTask;\r
indexlast2 = 0;\r
- index1 = 0;\r
char *hname[3];\r
for(Int_t i=0; i<3; i++) hname[i]= new char[50];\r
for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
- for(Int_t iside=0;iside<fgknSide;iside++){\r
- AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
- sprintf(hname[0],"chargeMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
- sprintf(hname[1],"TotalCharge_L%d_%d_%d_%d",lay,lad,det,iside);\r
- sprintf(hname[2],"hmonoDMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
- fModuleChargeMap[index1] = new TH2D(hname[0],hname[1],256,-0.5,255.5,256,-0.5,255.5);\r
- fModuleChargeMap[index1]->GetXaxis()->SetTitle("Time Bin");\r
- fModuleChargeMap[index1]->GetYaxis()->SetTitle("Anode");\r
- fAliITSQADataMakerRec->Add2RawsList(fModuleChargeMap[index1],indexlast1 + index1 + fRawsOffset);\r
- fSDDhRaws++;\r
- fmonoD[index1] = new TH1D(hname[2],hname[2],256,-0.5,255.5);\r
- index1++; \r
- indexlast2 = indexlast1 + index1;\r
- }\r
- }\r
- }\r
- \r
- \r
- TH2D *h7[fgknSDDmodules*2] ;\r
- index1 = 0;\r
- for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
- for(Int_t iside=0; iside<fgknSide; iside++){\r
- AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
- sprintf(hname0,"Anode_vs_Charge_L%d_%d_%d_%d",lay,lad,det,iside);\r
- h7[index1] = new TH2D(hname0,hname0,fgknSDDmodules*2,-0.5,-0.5+fgknSDDmodules*2,256,0.5,256.5);\r
- h7[index1]->GetXaxis()->SetTitle("Charge Value (ADC)");\r
- h7[index1]->GetYaxis()->SetTitle("Anode");\r
- fAliITSQADataMakerRec->Add2RawsList(h7[index1],indexlast2 + index1 + fRawsOffset);\r
- fSDDhRaws++;\r
- index1++;\r
- }\r
- }\r
- }\r
- /*\r
- TH2F *h8[3];\r
- for(Int_t htype =0; htype<3; htype++){\r
- if(htype == 0) sprintf(hname0,"BLcompare");\r
- if(htype == 1) sprintf(hname0,"NoiseCompare");\r
- if(htype == 2) sprintf(hname0,"VdriftCompare");\r
- h8[htype]= new TH2F(hname0,hname0,520,0,260,100,-0.5,99.5);\r
- if(htype == 0){\r
- h8[htype]->GetXaxis()->SetTitle("Module Number");\r
- h8[htype]->GetYaxis()->SetTitle("Time #");\r
- }\r
- if(htype == 1){\r
- h8[htype]->GetXaxis()->SetTitle("Module Number");\r
- h8[htype]->GetYaxis()->SetTitle("Time #)");\r
- }\r
- if(htype == 2){\r
- h8[htype]->GetXaxis()->SetTitle("Module Number");\r
- h8[htype]->GetYaxis()->SetTitle("Time #");\r
- }\r
+ for(Int_t iside=0;iside<fgknSide;iside++){\r
+ AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
+ sprintf(hname[0],"chargeMapFSE_L%d_%d_%d_%d",lay,lad,det,iside);\r
+ sprintf(hname[1],"ChargeMapForSingleEvent_L%d_%d_%d_%d",lay,lad,det,iside);\r
+ sprintf(hname[2],"hmonoDMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
+ TProfile2D *fModuleChargeMapFSE = new TProfile2D(hname[0],hname[1],256/fTimeBinSize,-0.5,255.5,256,-0.5,255.5);\r
+ fModuleChargeMapFSE->GetXaxis()->SetTitle("Time Bin");\r
+ fModuleChargeMapFSE->GetYaxis()->SetTitle("Anode");\r
+ fAliITSQADataMakerRec->Add2RawsList((new TProfile2D(*fModuleChargeMapFSE)),indexlast1 + index1 + fGenOffset);\r
+ delete fModuleChargeMapFSE;\r
+ \r
+ fSDDhTask++;\r
+ index1++; \r
+ indexlast2 = indexlast1 + index1;\r
+ }\r
+ }\r
\r
- fAliITSQADataMakerRec->Add2RawsList(h8[htype],indexlast2 + index1 + htype+1 + fRawsOffset);\r
- fSDDhRaws++;\r
- printf("added histos %c, at position %d \n", hname0, indexlast2 + index1 + htype+1 + fRawsOffset);\r
- }\r
- */\r
+ for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
+ for(Int_t iside=0;iside<fgknSide;iside++){\r
+ AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
+ sprintf(hname[0],"chargeMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
+ sprintf(hname[1],"ChargeMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
+ TProfile2D *fModuleChargeMap = new TProfile2D(hname[0],hname[1],256/fTimeBinSize,-0.5,255.5,256,-0.5,255.5);\r
+ fModuleChargeMap->GetXaxis()->SetTitle("Time Bin");\r
+ fModuleChargeMap->GetYaxis()->SetTitle("Anode");\r
+ fAliITSQADataMakerRec->Add2RawsList((new TProfile2D(*fModuleChargeMap)),indexlast1 + index1 + fGenOffset);\r
+ delete fModuleChargeMap;\r
+\r
+ fSDDhTask++;\r
+ index1++; \r
+ indexlast2 = indexlast1 + index1;\r
+ }\r
+ }\r
+ \r
} // kONLINE\r
\r
-\r
- AliDebug(1,Form("%d SDD Raws histograms booked\n",fSDDhRaws));\r
+ AliDebug(1,Form("%d SDD Raws histograms booked\n",fSDDhTask));\r
}\r
\r
\r
void AliITSQASDDDataMakerRec::MakeRaws(AliRawReader* rawReader)\r
{ \r
// Fill QA for RAW - SDD -\r
+ if(!fDDLModuleMap){\r
+ AliError("SDD DDL module map not available - skipping SDD QA");\r
+ return;\r
+ }\r
if(rawReader->GetType() != 7) return; // skips non physical triggers\r
AliDebug(1,"entering MakeRaws\n"); \r
rawReader->SelectEquipment(17,fgkeqOffset,fgkeqOffset + fgkDDLidRange); \r
\r
- /*\r
- if(rawReader->GetEventId()!=fEvtId){\r
- TFile *DAoutput = new TFile::Open(filename);\r
- TH1F *BLhisto;\r
- for(Int_t imod=0; imod<nSDDmodules; imod++){ \r
- BLhisto = (TH1F*)DAoutput->Get("BLhistoname[imod]");\r
- mean = BLhisto->take mean;\r
- fAliITSQADataMakerRec->GetRawsData(i+1887)->Fill(mean, imod);\r
- Noisehisto....;\r
- Vdrifthisto...;\r
- <Q>histo....;\r
- }\r
- }\r
- fEvtId==rawReader->GetEventId();\r
- */\r
-\r
rawReader->Reset(); \r
AliITSRawStreamSDD s(rawReader); \r
- s.SetDDLModuleMap(fSDDDDLModuleMap);\r
+ s.SetDDLModuleMap(fDDLModuleMap);\r
Int_t lay, lad, det; \r
\r
Int_t index=0;\r
if(fkOnline) {\r
for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
for(Int_t iside=0;iside<fgknSide;iside++) {\r
- if(fSDDhRaws > 39+12 * 132 + index) fAliITSQADataMakerRec->GetRawsData(39+12 * 132 + index +fRawsOffset)->Reset();\r
+ if(fSDDhTask > 39 + index) fAliITSQADataMakerRec->GetRawsData(39 + index +fGenOffset)->Reset();\r
index++;\r
}\r
- }\r
}\r
}\r
Int_t cnt = 0;\r
Int_t ildcID = -1;\r
Int_t iddl = -1;\r
Int_t isddmod = -1;\r
- Int_t coord1, coord2, signal, moduleSDD, ioffset, iorder, activeModule, index1 ,nBins , iside, bin ;\r
+ Int_t coord1, coord2, signal, moduleSDD, ioffset, iorder, activeModule, index1;\r
while(s.Next()) {\r
ildcID = rawReader->GetLDCId();\r
iddl = rawReader->GetDDLID() - fgkDDLIDshift;\r
AliDebug(1,Form( "Module SDD = %d, resetting it to 1 \n",moduleSDD));\r
moduleSDD = 1;\r
}\r
- fAliITSQADataMakerRec->GetRawsData(0 +fRawsOffset)->Fill(moduleSDD); \r
+ fAliITSQADataMakerRec->GetRawsData(0 +fGenOffset)->Fill(moduleSDD); \r
\r
AliITSgeomTGeo::GetModuleId(isddmod, lay, lad, det);\r
ioffset = 3;\r
ioffset += 14;\r
iorder = 2; \r
} \r
- fAliITSQADataMakerRec->GetRawsData(iorder +fRawsOffset)->Fill(lad);\r
- fAliITSQADataMakerRec->GetRawsData(ioffset+lad-1 +fRawsOffset)->Fill(det); //-1 because ladder# starts from 1 \r
+ fAliITSQADataMakerRec->GetRawsData(iorder +fGenOffset)->Fill(lad);\r
+ fAliITSQADataMakerRec->GetRawsData(ioffset+lad-1 +fGenOffset)->Fill(det); //-1 because ladder# starts from 1 \r
\r
Short_t iside = s.GetChannel();\r
activeModule = moduleSDD;\r
- if(moduleSDD > 35) activeModule -= 48;\r
index1 = activeModule * 2 + iside;\r
\r
if(index1<0){\r
}\r
\r
if(fkOnline) {\r
- if(fSDDhRaws > 39+12 * 132 + index1) {\r
- \r
- fAliITSQADataMakerRec->GetRawsData(39+ 2 * 132 + index1 +fRawsOffset)->Fill(coord1); \r
- fAliITSQADataMakerRec->GetRawsData(39+ 4 * 132 + index1 +fRawsOffset)->Fill(coord1,signal); \r
- fAliITSQADataMakerRec->GetRawsData(39+ 6 * 132 + index1 +fRawsOffset)->Fill(coord2); \r
- fAliITSQADataMakerRec->GetRawsData(39+ 8 * 132 + index1 +fRawsOffset)->Fill(coord2,signal); \r
- fAliITSQADataMakerRec->GetRawsData(39+12 * 132 + index1 +fRawsOffset)->Fill(signal,coord1);\r
- ((TH2D *)(fAliITSQADataMakerRec->GetRawsData(39+10 * 132 + index1 +fRawsOffset)))->Fill(coord2, coord1, signal);\r
-\r
+ if(fSDDhTask > 39 + index1) {\r
+ ((TProfile2D *)(fAliITSQADataMakerRec->GetRawsData(39 + index1 +fGenOffset)))->Fill(coord2, coord1, signal);\r
+ ((TProfile2D *)(fAliITSQADataMakerRec->GetRawsData(39 + index1 + 260*2 +fGenOffset)))->Fill(coord2, coord1, signal);\r
}\r
}\r
cnt++;\r
if(!(cnt%10000)) AliDebug(1,Form(" %d raw digits read",cnt));\r
}\r
- AliDebug(1,Form("Event completed, %d raw digits read",cnt)); \r
-\r
- if(fkOnline) {\r
- nBins = 256;\r
- for(moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
- if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
- for(iside=0; iside<fgknSide; iside++){\r
- activeModule = moduleSDD;\r
- if(moduleSDD > 35) activeModule -= 48;\r
- index1 = activeModule * 2 + iside;\r
- if(fSDDhRaws > 39 + 2 * 132 + index1) {\r
- fmonoD[index1] = ((TH2D *) (fAliITSQADataMakerRec->GetRawsData(39+10 * 132 + index1 +fRawsOffset)))->ProjectionY();\r
- for(bin=0; bin<nBins; bin++) fAliITSQADataMakerRec->GetRawsData(index1+39 +fRawsOffset)->Fill(bin,fmonoD[index1]->GetBinContent(bin+1) );\r
- } \r
- }\r
- }\r
- }\r
- for(Int_t i=0; i<fSDDhRaws; i++){\r
- Int_t entries = static_cast<Int_t>(fAliITSQADataMakerRec->GetRawsData(i +fRawsOffset)->GetEntries());\r
- if(entries != 0)\r
- AliDebug(1,Form("histo %d, name %s , entries %d ",i,fAliITSQADataMakerRec->GetRawsData(i +fRawsOffset)->GetName(),entries));\r
- }\r
- }\r
-}\r
+ AliDebug(1,Form("Event completed, %d raw digits read",cnt)); \r
+ }\r
\r
//____________________________________________________________________________ \r
void AliITSQASDDDataMakerRec::InitRecPoints()\r
{\r
// Initialization for RECPOINTS - SDD -\r
- fRecsOffset = (fAliITSQADataMakerRec->fRecPointsQAList)->GetEntries();\r
+ fGenOffset = (fAliITSQADataMakerRec->fRecPointsQAList)->GetEntries();\r
\r
- TH1F * h0 = new TH1F("Lay3TotCh","Layer 3 total charge",1000,-0.5, 499.5);\r
+ TH1F *h0 = new TH1F("Lay3TotCh","Layer 3 total charge",1000,-0.5, 499.5);\r
h0->GetXaxis()->SetTitle("ADC value");\r
h0->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h0, 0 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h0)), 0 +fGenOffset);\r
+ delete h0;\r
+ fSDDhTask++;\r
\r
- TH1F * h1 = new TH1F("Lay4TotCh","Layer 4 total charge",1000,-0.5, 499.5);\r
+ TH1F *h1 = new TH1F("Lay4TotCh","Layer 4 total charge",1000,-0.5, 499.5);\r
h1->GetXaxis()->SetTitle("ADC value");\r
h1->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h1, 1 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h1)), 1 +fGenOffset);\r
+ delete h1;\r
+ fSDDhTask++;\r
\r
\r
char hisnam[50];\r
- TH1F *h2[3]; \r
for(Int_t i=1; i<=3; i++){\r
sprintf(hisnam,"Charge_L3_Strip%d",i);\r
- h2[i] = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
- fAliITSQADataMakerRec->Add2RecPointsList(h2[i],i+1 +fRecsOffset);\r
- fSDDhRecs++;\r
+ TH1F *h2 = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h2)),i+1 +fGenOffset);\r
+ delete h2;\r
+ fSDDhTask++;\r
}\r
\r
- TH1F *h3[4]; \r
for(Int_t i=1; i<=4; i++){\r
sprintf(hisnam,"Charge_L4_Strip%d",i);\r
- h3[i] = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
- fAliITSQADataMakerRec->Add2RecPointsList(h3[i],i+4 +fRecsOffset);\r
- fSDDhRecs++;\r
+ TH1F *h3 = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h3)),i+4 +fGenOffset);\r
+ delete h3;\r
+ fSDDhTask++;\r
}\r
\r
TH1F *h4 = new TH1F("ModPatternRP","Modules pattern RP",fgknSDDmodules,239.5,499.5); \r
h4->GetXaxis()->SetTitle("Module number");\r
h4->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h4,9 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h4)),9 +fGenOffset);\r
+ delete h4;\r
+ fSDDhTask++;\r
TH1F *h5 = new TH1F("ModPatternL3 RP","Ladder pattern L3 RP",14,0.5,14.5); \r
h5->GetXaxis()->SetTitle("Ladder #, Layer 3");\r
h5->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h5,10 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h5)),10 +fGenOffset);\r
+ delete h5;\r
+ fSDDhTask++;\r
TH1F *h6 = new TH1F("ModPatternL4 RP","Ladder pattern L4 RP",22,0.5,22.5); \r
h6->GetXaxis()->SetTitle("Ladder #, Layer 4");\r
h6->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h6,11 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h6)),11 +fGenOffset);\r
+ delete h6;\r
+ fSDDhTask++;\r
TH2F *h7 = new TH2F("Local Coord Distrib","Local Coord Distrib",1000,-4,4,1000,-4,4);\r
h7->GetXaxis()->SetTitle("X local coord, drift, cm");\r
h7->GetYaxis()->SetTitle("Z local coord, anode, cm");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h7,12 +fRecsOffset);\r
- fSDDhRecs++;\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH2F(*h7)),12 +fGenOffset);\r
+ delete h7;\r
+ fSDDhTask++;\r
TH2F *h8 = new TH2F("Global Coord Distrib","Global Coord Distrib",6000,-30,30,6000,-30,30);\r
h8->GetYaxis()->SetTitle("Y glob coord, cm");\r
h8->GetXaxis()->SetTitle("X glob coord, cm");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h8,13 +fRecsOffset);\r
- fSDDhRecs++;\r
- \r
- TH1F *h9[2]; \r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH2F(*h8)),13 +fGenOffset);\r
+ delete h8;\r
+ fSDDhTask++;\r
+ \r
for(Int_t iLay=0; iLay<=1; iLay++){\r
sprintf(hisnam,"hr_Layer%d",iLay+3);\r
- h9[iLay] = new TH1F(hisnam,hisnam,100,10.,30.);\r
- h9[iLay]->GetXaxis()->SetTitle("r (cm)");\r
- h9[iLay]->GetXaxis()->CenterTitle();\r
- h9[iLay]->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h9[iLay],iLay+14 +fRecsOffset);\r
- fSDDhRecs++;\r
+ TH1F *h9 = new TH1F(hisnam,hisnam,100,10.,30.);\r
+ h9->GetXaxis()->SetTitle("r (cm)");\r
+ h9->GetXaxis()->CenterTitle();\r
+ h9->GetYaxis()->SetTitle("Entries");\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h9)),iLay+14 +fGenOffset);\r
+ delete h9;\r
+ fSDDhTask++;\r
}\r
- TH1F *h10[2]; \r
+\r
for(Int_t iLay=0; iLay<=1; iLay++){\r
sprintf(hisnam,"hphi_Layer%d",iLay+3);\r
- h10[iLay] = new TH1F(hisnam,hisnam,100,-TMath::Pi(),TMath::Pi());\r
- h10[iLay]->GetXaxis()->SetTitle("#varphi (rad)");\r
- h10[iLay]->GetXaxis()->CenterTitle();\r
- h10[iLay]->GetYaxis()->SetTitle("Entries");\r
- fAliITSQADataMakerRec->Add2RecPointsList(h10[iLay],iLay+16 +fRecsOffset);\r
- fSDDhRecs++;\r
+ TH1F *h10 = new TH1F(hisnam,hisnam,100,-TMath::Pi(),TMath::Pi());\r
+ h10->GetXaxis()->SetTitle("#varphi (rad)");\r
+ h10->GetXaxis()->CenterTitle();\r
+ h10->GetYaxis()->SetTitle("Entries");\r
+ fAliITSQADataMakerRec->Add2RecPointsList((new TH1F(*h10)),iLay+16 +fGenOffset);\r
+ delete h10;\r
+ fSDDhTask++;\r
}\r
\r
- AliDebug(1,Form("%d SDD Recs histograms booked\n",fSDDhRecs));\r
+ AliDebug(1,Form("%d SDD Recs histograms booked\n",fSDDhTask));\r
}\r
\r
//____________________________________________________________________________ \r
AliError("can't get the branch with the ITS clusters !");\r
return;\r
}\r
- TClonesArray * recpoints = new TClonesArray("AliITSRecPoint") ;\r
+ static TClonesArray statRecpoints("AliITSRecPoint") ;\r
+ TClonesArray *recpoints = &statRecpoints;\r
branchRecP->SetAddress(&recpoints);\r
Int_t npoints = 0; \r
Float_t cluglo[3]={0.,0.,0.}; \r
\r
for(Int_t j=0;j<recpoints->GetEntries();j++){\r
AliITSRecPoint *recp = (AliITSRecPoint*)recpoints->At(j); \r
- fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
+ fAliITSQADataMakerRec->GetRecPointsData(9 +fGenOffset)->Fill(module);\r
recp->GetGlobalXYZ(cluglo);\r
Float_t rad=TMath::Sqrt(cluglo[0]*cluglo[0]+cluglo[1]*cluglo[1]); \r
Float_t phi=TMath::ATan2(cluglo[1],cluglo[0]);\r
if(recp->GetLayer() ==2) {\r
- fAliITSQADataMakerRec->GetRecPointsData(0 +fRecsOffset)->Fill(recp->GetQ()) ;\r
- fAliITSQADataMakerRec->GetRecPointsData(10 +fRecsOffset)->Fill(lad);\r
- fAliITSQADataMakerRec->GetRecPointsData(14 +fRecsOffset)->Fill(rad);\r
- fAliITSQADataMakerRec->GetRecPointsData(16 +fRecsOffset)->Fill(phi);\r
- fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
- fAliITSQADataMakerRec->GetRecPointsData(12 +fRecsOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
- fAliITSQADataMakerRec->GetRecPointsData(13 +fRecsOffset)->Fill(cluglo[0],cluglo[1]);\r
+ fAliITSQADataMakerRec->GetRecPointsData(0 +fGenOffset)->Fill(recp->GetQ()) ;\r
+ fAliITSQADataMakerRec->GetRecPointsData(10 +fGenOffset)->Fill(lad);\r
+ fAliITSQADataMakerRec->GetRecPointsData(14 +fGenOffset)->Fill(rad);\r
+ fAliITSQADataMakerRec->GetRecPointsData(16 +fGenOffset)->Fill(phi);\r
+ fAliITSQADataMakerRec->GetRecPointsData(9 +fGenOffset)->Fill(module);\r
+ fAliITSQADataMakerRec->GetRecPointsData(12 +fGenOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
+ fAliITSQADataMakerRec->GetRecPointsData(13 +fGenOffset)->Fill(cluglo[0],cluglo[1]);\r
}\r
else if(recp->GetLayer() ==3) {\r
- fAliITSQADataMakerRec->GetRecPointsData(1 +fRecsOffset)->Fill(recp->GetQ()) ;\r
- fAliITSQADataMakerRec->GetRecPointsData(11 +fRecsOffset)->Fill(lad);\r
- fAliITSQADataMakerRec->GetRecPointsData(15 +fRecsOffset)->Fill(rad);\r
- fAliITSQADataMakerRec->GetRecPointsData(17 +fRecsOffset)->Fill(phi);\r
- fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
- fAliITSQADataMakerRec->GetRecPointsData(12 +fRecsOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
- fAliITSQADataMakerRec->GetRecPointsData(13 +fRecsOffset)->Fill(cluglo[0],cluglo[1]);\r
+ fAliITSQADataMakerRec->GetRecPointsData(1 +fGenOffset)->Fill(recp->GetQ()) ;\r
+ fAliITSQADataMakerRec->GetRecPointsData(11 +fGenOffset)->Fill(lad);\r
+ fAliITSQADataMakerRec->GetRecPointsData(15 +fGenOffset)->Fill(rad);\r
+ fAliITSQADataMakerRec->GetRecPointsData(17 +fGenOffset)->Fill(phi);\r
+ fAliITSQADataMakerRec->GetRecPointsData(9 +fGenOffset)->Fill(module);\r
+ fAliITSQADataMakerRec->GetRecPointsData(12 +fGenOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
+ fAliITSQADataMakerRec->GetRecPointsData(13 +fGenOffset)->Fill(cluglo[0],cluglo[1]);\r
}\r
}\r
}\r
- recpoints->Delete();\r
- delete recpoints;\r
-\r
+ statRecpoints.Clear();\r
}\r
\r