+
+
+
+ if(fSubDetector == 0){
+
+ // Check id histograms already created for this Event Specie
+ AliITSRecPointContainer* rpcont=AliITSRecPointContainer::Instance();
+ TClonesArray *recpoints =NULL;
+ if(fkOnline){
+ recpoints= rpcont->FetchClusters(0,clustersTree,GetEventNumber());
+ }
+ else{
+ recpoints= rpcont->FetchClusters(0,clustersTree);
+ }
+ if(!rpcont->GetStatusOK()){
+ AliError("cannot access to ITS recpoints");
+ return;
+ }
+
+ Int_t offset = fRecPointsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries();
+ Float_t cluGlo[3] = {0.,0.,0.};
+ Int_t lay, lad, det;
+ // Fill QA for recpoints
+ for(Int_t module=0; module<rpcont->GetNumberOfModules();module++){
+ // AliInfo(Form("Module %d\n",module));
+ recpoints = rpcont->UncheckedGetClusters(module);
+ AliITSgeomTGeo::GetModuleId(module, lay, lad, det);
+ for(Int_t j=0;j<recpoints->GetEntries();j++){
+ AliITSRecPoint *rcp = (AliITSRecPoint*)recpoints->At(j);
+ //Check id histograms already created for this Event Specie
+ rcp->GetGlobalXYZ(cluGlo);
+ Double_t rad=TMath::Sqrt(cluGlo[0]*cluGlo[0]+cluGlo[1]*cluGlo[1]+cluGlo[2]*cluGlo[2]);
+ Double_t phi= TMath::Pi() + TMath::ATan2(-cluGlo[1],-cluGlo[0]);
+ Double_t theta = TMath::ACos(cluGlo[2]/rad);
+ Double_t eta = 100.;
+ if(AreEqual(rad,0.) == kFALSE) {
+ if(theta<=1.e-14){ eta=30.; }
+ else { eta = -TMath::Log(TMath::Tan(theta/2.));}
+ }
+ // printf("=========================>hlt rcp->GetLayer() = %d \n",rcp->GetLayer());
+ (GetRecPointsData( rcp->GetLayer() + offset - 6))->Fill(eta,phi);
+ }
+ }
+ }
+
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerRec::FillRecPoint(AliITSRecPoint rcp)
+{
+
+ // Fill QA for recpoints
+ Float_t cluGlo[3] = {0.,0.,0.};
+ Int_t offset = fRecPointsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries();
+ // Check id histograms already created for this Event Specie
+ rcp.GetGlobalXYZ(cluGlo);
+ Double_t rad=TMath::Sqrt(cluGlo[0]*cluGlo[0]+cluGlo[1]*cluGlo[1]+cluGlo[2]*cluGlo[2]);
+ Double_t phi= TMath::Pi() + TMath::ATan2(-cluGlo[1],-cluGlo[0]);
+ Double_t theta = TMath::ACos(cluGlo[2]/rad);
+ Double_t eta = 100.;
+ if(AreEqual(rad,0.)==kFALSE) {
+ if(theta<=1.e-14){eta=30.;}
+ else {eta = -TMath::Log(TMath::Tan(theta/2.));}
+ }
+ (GetRecPointsData( rcp.GetLayer() + offset - 6))->Fill(eta,phi);
+