Makeing defualt plots - report from the raw data QA

(Marian)

diff --git a/TPC/CalibMacros/CalibQA.C b/TPC/CalibMacros/CalibQA.C
new file mode 100644 (file)
index 0000000..4178365
--- /dev/null
+++ b/TPC/CalibMacros/CalibQA.C
@@ -0,0 +1,111 @@
+/*
+  Make default plotsfrom AliTPCdataQA components:
+
+  aliroot -b -q  $ALICE_ROOT/TPC/CalibMacros/CalibQA.C\(121694\);
+  
+  .L $ALICE_ROOT/TPC/CalibMacros/CalibQA.C
+  Int_t run=121694;
+  CalibQA(run);
+*/
+
+TCut cutNoise="PadNoise.fElements<1.5&&abs(PadNoise.fElements/PadNoise_Median-1)<0.5";
+TCut cutTime="abs(TimePosition.fElements-TimePosition_Median)<100";
+TCut cutOccu="abs(NoThreshold.fElements/NoThreshold_Median-1)<0.9";
+TCut cutAmp="abs(MaxCharge.fElements/MaxCharge_Median-1)<0.99";
+TCut cutIROC="sector<36";
+TCut cutOROC="sector>=36";
+
+
+void CalibQA(Int_t run){
+  InitOCDB(run);
+  MakeTree();
+  TCanvas *canvas=0;
+  //
+  TPostScript *ps = new TPostScript("rawQA.ps", 112);  
+  ps->NewPage();
+  canvas=DrawOccupancy();
+  ps->NewPage();
+  canvas->Update();
+  ps->Close();
+  delete ps;
+
+}
+
+void InitOCDB(Int_t run){
+  gSystem->Load("libANALYSIS");
+  gSystem->Load("libTPCcalib");
+  gROOT->LoadMacro("$ALICE_ROOT/TPC/scripts/OCDBscan/ConfigOCDBLustre.C");
+  //gROOT->LoadMacro("$ALICE_ROOT/TPC/scripts/OCDBscan/ConfigOCDB.C");
+  gROOT->Macro("$ALICE_ROOT/TPC/scripts/OCDBscan/NimStyle.C");
+  ConfigOCDB(run);
+}
+
+void MakeTree(){
+  //
+  // make summary tree
+  //
+  AliTPCcalibDB::Instance()->UpdateNonRec();
+  AliTPCdataQA* dataQA =   AliTPCcalibDB::Instance()->GetDataQA();
+  AliTPCCalPad* gain   =   AliTPCcalibDB::Instance()->GetDedxGainFactor();
+  AliTPCCalPad * padNoise = AliTPCcalibDB::Instance()->GetPadNoise();
+  AliTPCPreprocessorOnline preprocesor;
+  gain->SetName("krGain");
+  preprocesor.AddComponent(gain);
+  preprocesor.AddComponent(dataQA->GetNPads());
+  preprocesor.AddComponent(dataQA->GetNTimeBins());
+  preprocesor.AddComponent(dataQA->GetMaxCharge());
+  preprocesor.AddComponent(dataQA->GetNoThreshold());
+  preprocesor.AddComponent(dataQA->GetNLocalMaxima());
+  preprocesor.AddComponent(dataQA->GetTimePosition());
+  preprocesor.AddComponent(padNoise);
+  preprocesor.DumpToFile("QA.root");
+}
+
+
+TCanvas * DrawOccupancy(){
+  TH1::AddDirectory(0);
+  gStyle->SetOptStat(0);
+  TFile f("QA.root");
+  TTree  * tree = (TTree*)f.Get("calPads");
+  TLegend *legend=0;
+  TProfile * phoccGainIROC=0;
+  TProfile * phoccGainOROC=0;
+  //
+  TCanvas * canvas = new TCanvas("occupancy","occupancy",700,700);
+  canvas->Divide(2,2);
+  canvas->cd(1);
+  tree->Draw("NoThreshold.fElements:gy.fElements:gx.fElements>>hisOccuA(250,-250,250,250,-250,250)",cutNoise+cutTime+cutOccu+"sector%36<18","profcolz");
+  canvas->cd(2);
+  tree->Draw("NoThreshold.fElements:gy.fElements:gx.fElements>>hisOccuC(250,-250,250,250,-250,250)",cutNoise+cutTime+cutOccu+"sector%36>=18","profcolz");
+
+  canvas->cd(3);
+  tree->Draw("NoThreshold.fElements:krGain.fElements>>hoccGainIROC(20,0.7,1.2)",cutNoise+cutTime+cutOccu+cutIROC,"prof");
+  legend = new TLegend(0.45,0.15,0.85,0.35, "Raw cluster occupancy");
+  phoccGainIROC = (TProfile*)(gROOT->FindObject("hoccGainIROC")->Clone());
+  phoccGainIROC->Draw();
+  phoccGainIROC->SetTitle("IROC");
+  phoccGainIROC->SetName("IROC");
+  phoccGainIROC->GetXaxis()->SetTitle("Krypton Amp (a.u.)");
+  legend->AddEntry(phoccGainIROC);
+  legend->Draw();
+
+  canvas->cd(4);
+  tree->Draw("NoThreshold.fElements:krGain.fElements>>hoccGainOROC(20,0.8,1.2)",cutNoise+cutTime+cutOccu+cutOROC+"lx.fElements<200","prof");
+  phoccGainOROC = (TProfile*)(gROOT->FindObject("hoccGainOROC")->Clone());
+  phoccGainOROC->Draw();
+  phoccGainOROC->SetTitle("OROC");
+  phoccGainOROC->SetName("OROC");
+  phoccGainOROC->GetXaxis()->SetTitle("Krypton Amp (a.u.)");
+  legend = new TLegend(0.45,0.15,0.85,0.35, "Raw cluster occupancy");
+  legend->AddEntry(phoccGainOROC);
+  legend->Draw();
+  return canvas;
+}
+
+TCanvas * DrawGain(){
+  //
+  // Compare the amplitude with krypton gain amplitude
+  // Similar filtering as in occupancy plot
+  //                                           
+  return 0;
+}