Flexible pt range for the efficiency histogramming

[u/mrichter/AliRoot.git] / TOF / CheckEnabledChannels.C

diff --git a/TOF/CheckEnabledChannels.C b/TOF/CheckEnabledChannels.C
index 63e0f8365e2b2032226ce0d8253d3909543e30d1..09ebe23c812aaa0073d943428aa9f96d26cb47bf 100644 (file)
--- a/TOF/CheckEnabledChannels.C
+++ b/TOF/CheckEnabledChannels.C
@@ -1,4 +1,50 @@
-CheckEnabledChannels(Int_t run, const Char_t *dbString)
+CheckEnabledChannels(const Char_t *runlist)
+{
+
+  ifstream is(runlist);
+  Char_t buf[4096];
+  Int_t run[1024];
+  Int_t nrun = 0;
+  while(!is.eof()) {
+    is.getline(buf, 4096);
+    if (is.eof()) break;
+    run[nrun] = atoi(buf);
+    printf("added run number %d\n", run[nrun]);
+    nrun++;
+  }
+  printf("%d runs added\n", nrun);
+  is.close();
+
+  TH1F *hActive = new TH1F("hActive", "active channels;run;fraction", nrun, 0, nrun);
+  TH1F *hReadout = new TH1F("hReadout", "good readout;run;fraction", nrun, 0, nrun);
+  for (Int_t irun = 0; irun < nrun; irun++) {
+    hr = CheckEnabledChannels(run[irun], kTRUE);
+    ha = CheckEnabledChannels(run[irun], kFALSE);
+    hReadout->SetBinContent(irun + 1, hr->Integral());
+    hActive->SetBinContent(irun + 1, ha->Integral());
+    hReadout->GetXaxis()->SetBinLabel(irun + 1, Form("%d", run[irun]));
+    delete hr; delete ha;
+  }
+  
+  hReadout->SetMarkerStyle(20);
+  hReadout->SetMarkerColor(4);
+  hActive->SetMarkerStyle(25);
+  hActive->SetMarkerColor(2);
+  hReadout->Sumw2();
+  hActive->Sumw2();
+  hReadout->Divide(hReadout, hActive, 1., 1., "B");
+  hActive->Scale(1. / 152928.);
+  hReadout->SetMinimum(0.);
+  hReadout->SetMaximum(1.);
+  hReadout->Draw("E");
+  hActive->Draw("E, same");
+  TLegend *l = gPad->BuildLegend();
+  l->SetFillStyle(0);
+
+}
+
+TH1F *
+CheckEnabledChannels(Int_t run, Bool_t checkROEff = kTRUE, const Char_t *dbString = "raw://")
 {
 
   /* init */
@@ -9,6 +55,7 @@ CheckEnabledChannels(Int_t run, const Char_t *dbString)
   calib.Init();
 
   TH2F *hEnabledMap = new TH2F("hEnabledMap", "Enabled channel map;sector;strip", 72, 0., 18., 91, 0., 91.);
+  TH1F *hEnabledFlag = new TH1F("hEnabledFlag", "Enabled channel flag;index;flag", 157248, 0., 157248.);
 
   AliTOFcalibHisto calibhisto;
   calibhisto.LoadCalibHisto();
@@ -18,7 +65,7 @@ CheckEnabledChannels(Int_t run, const Char_t *dbString)
   Float_t hitmapx, hitmapy;
   /* loop over channels */
   for (Int_t ich = 0; ich < 157248; ich++) {
-    if (!calib.IsChannelEnabled(ich, kTRUE, kTRUE)) continue;
+    if (!calib.IsChannelEnabled(ich, checkROEff)) continue;
     sector = calibhisto.GetCalibMap(AliTOFcalibHisto::kSector, ich);
     sectorStrip = calibhisto.GetCalibMap(AliTOFcalibHisto::kSectorStrip, ich);
     padx = calibhisto.GetCalibMap(AliTOFcalibHisto::kPadX, ich);
@@ -26,11 +73,10 @@ CheckEnabledChannels(Int_t run, const Char_t *dbString)
     hitmapx = sector + ((Double_t)(3 - fea) + 0.5) / 4.;
     hitmapy = sectorStrip;
     hEnabledMap->Fill(hitmapx, hitmapy);
+    hEnabledFlag->SetBinContent(ich + 1, 1);
   }
   
   hEnabledMap->DrawCopy("colz");
-  TFile *fileout = TFile::Open("CheckEnabledChannels.root", "RECREATE");
-  hEnabledMap->Write();
-  fileout->Close();
+  return hEnabledFlag;
 
 }