if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeRecPoints(clustersTree);
+
+
+
// Check id histograms already created for this Event Specie
TBranch *branchRecP = clustersTree->GetBranch("ITSRecPoints");
if (!branchRecP) {
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
+ //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]);
Double_t phi= TMath::Pi() + TMath::ATan2(-cluGlo[1],-cluGlo[0]);
Double_t theta = TMath::ACos(cluGlo[2]/rad);
Double_t eta = 100.;
- if(rad != 0) eta = -TMath::Log(TMath::Tan(theta/2.));
+ if(IsEqual(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);
}
}
Double_t phi= TMath::Pi() + TMath::ATan2(-cluGlo[1],-cluGlo[0]);
Double_t theta = TMath::ACos(cluGlo[2]/rad);
Double_t eta = 100.;
+ if(IsEqual(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);
+ /*
if(rad != 0) eta = -TMath::Log(TMath::Tan(theta/2.));
(GetRecPointsData( rcp.GetLayer() + offset - 6))->Fill(eta,phi);
+ */
}
//____________________________________________________________________________
}
//return offset;
}
+
+//____________________________________________________________________
+
+Bool_t AliITSQADataMakerRec::IsEqual(Double_t a1,Double_t a2)
+{
+ const Double_t kEpsilon= 0.000000001;
+ return TMath::Abs(a1-a2)<=kEpsilon*TMath::Abs(a1);
+}
+
virtual ~AliITSQADataMakerRec(); // dtor
Short_t GetSubDet(){return fSubDetector;};
Int_t GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task);
- TH2F *GetITSGlobalHisto(Int_t layer);
+ TH2F *GetITSGlobalHisto(Int_t layer);
+ static Bool_t IsEqual(Double_t a1, Double_t a2);
private:
//____________________________________________________________________________
void AliITSQASDDDataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t /*task*/, TObjArray* /*list*/)
{
+
// launch the QA checking
if(fkOnline) {
AnalyseBNG(); // Analyse Baseline, Noise, Gain
AnalyseINJ(); // Analyse Injectors
}
+
AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
}
TH2D *hphil3 = new TH2D("SDDphizL3","SDD #varphiz Layer3 ",6,0.5,6.5,14,0.5,14.5);
hphil3->GetXaxis()->SetTitle("z[#Module L3 ]");
hphil3->GetYaxis()->SetTitle("#varphi[#Ladder L3]");
- rv = fAliITSQADataMakerRec->Add2RawsList((new TH2D(*hphil3)),1+fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()], expert, image, saveCorr);
+ rv = fAliITSQADataMakerRec->Add2RawsList((new TH2D(*hphil3)),1+fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()], !expert, image, saveCorr);
delete hphil3;
fSDDhRawsTask++;
TH2D *hphil4 = new TH2D("SDDphizL4","SDD #varphiz Layer4 ",8,0.5,8.5,22,0.5,22.5);
hphil4->GetXaxis()->SetTitle("z[#Module L4]");
hphil4->GetYaxis()->SetTitle("#varphi[#Ladder L4]");
- rv = fAliITSQADataMakerRec->Add2RawsList((new TH2D(*hphil4)),2+fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()], expert, image, saveCorr);
+ rv = fAliITSQADataMakerRec->Add2RawsList((new TH2D(*hphil4)),2+fGenRawsOffset[fAliITSQADataMakerRec->GetEventSpecie()], !expert, image, saveCorr);
delete hphil4;
fSDDhRawsTask++;