Minor updates in macros used to plot SDD calib. quantities vs. time

[u/mrichter/AliRoot.git] / ITS / PlotDriftSpeedSDDVsTime.C

diff --git a/ITS/PlotDriftSpeedSDDVsTime.C b/ITS/PlotDriftSpeedSDDVsTime.C
index 58ee8d4e80939503d218d5d50e14a7ec355850d2..259050e6b9f093255fcabf2c56e13a05eee5fd44 100644 (file)
--- a/ITS/PlotDriftSpeedSDDVsTime.C
+++ b/ITS/PlotDriftSpeedSDDVsTime.C
@@ -2,10 +2,12 @@
 #include <Riostream.h>
 #include <TFile.h>
 #include <TH1F.h>
+#include <TH2F.h>
 #include <TGraph.h>
 #include <TStyle.h>
 #include <TLegend.h>
 #include <TLegendEntry.h>
+#include <TLine.h>
 #include <TSystem.h>
 #include <TGrid.h>
 #include <TMath.h>
@@ -39,6 +41,15 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
   Char_t filnam[200],filnamalien[200];
   TGraphErrors** gvdrvstime=new TGraphErrors*[520];
   TGraphErrors** gvdrvsrun=new TGraphErrors*[520];
+  TGraph* gGoodInjVsRun=new TGraph(0);
+  gGoodInjVsRun->SetName("gGoodInjVsRun");
+  TGraph* gAverSpeedVsRun=new TGraph(0);
+  gAverSpeedVsRun->SetName("gAverSpeedVsRun");
+  TGraph* gGoodInjVsTime=new TGraph(0);
+  gGoodInjVsTime->SetName("gGoodInjVsTime");
+  TGraph* gAverSpeedVsTime=new TGraph(0);
+  gAverSpeedVsTime->SetName("gAverSpeedVsIime");
+  
   for(Int_t iMod=0; iMod<260;iMod++){
     gvdrvstime[iMod]=new TGraphErrors(0);    
     gvdrvstime[iMod]->SetTitle(Form("Module %d",iMod+240));
@@ -67,27 +78,48 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
     AliCDBEntry *ent=(AliCDBEntry*)f->Get("AliCDBEntry");
     TObjArray *drspSDD = (TObjArray *)ent->GetObject();
     
+    Int_t iGoodInj=0;
+    Int_t iAverSpeed=0;
+    Float_t timeday=0;
+    
     AliITSDriftSpeedArraySDD *vdriftarr;
+    AliITSDriftSpeedArraySDD *vdriftarr0;
+    AliITSDriftSpeedArraySDD *vdriftarr1;
+
     for(Int_t iMod=0; iMod<260;iMod++){
       Int_t index=-1;
       if(anode<256) index=2*iMod;
       else index=2*iMod+1;
       vdriftarr=(AliITSDriftSpeedArraySDD*)drspSDD->At(index);
+         
+      Int_t i0=2*iMod;
+      Int_t i1=1+2*iMod;
+      vdriftarr0=(AliITSDriftSpeedArraySDD*)drspSDD->At(i0);
+      vdriftarr1=(AliITSDriftSpeedArraySDD*)drspSDD->At(i1);
+         
+      Int_t statusInj0=vdriftarr0->GetInjectorStatus();
+      Int_t statusInj1=vdriftarr1->GetInjectorStatus();
+      if(statusInj0>0) iGoodInj++;
+      else iAverSpeed++;
+      if(statusInj1>0) iGoodInj++;
+      else iAverSpeed++;
+         
       Int_t iAn=anode;
       if(anode>256) iAn=anode-256;
       Float_t vdrift=vdriftarr->GetDriftSpeed(0,iAn);
       if(vdrift<4. || vdrift > 8.) continue;
-      UInt_t timest=vdriftarr->GetTimestamp(0);
-      if(timest==0) continue;
+      if(statusInj0==0) continue;
       Int_t npt=gvdrvsrun[iMod]->GetN();
       gvdrvsrun[iMod]->SetPoint(npt,(Float_t)nrun,vdrift);
       gvdrvsrun[iMod]->SetPointError(npt,0,errSpeed[iMod]);
       
+      UInt_t timest=vdriftarr->GetTimestamp(0);
+      if(timest==0) continue;
       Float_t timeZero;
       if(year==2009) timeZero=1247762992;
       else timeZero=1262300400;
       if(timest<timeZero) continue;
-      Float_t timeday=float(timest-timeZero)/60./60./24.;
+      timeday=float(timest-timeZero)/60./60./24.;
       Float_t mob=vdrift*1.E5/Edrift;  
       Float_t temper=293.15*TMath::Power((mob/1350.),-1/2.4); 
       if(iMod==497-240) printf("Run %s   Time %d Day %f Speed=%f Temp=%f\n",filnam,timest,timeday,vdrift,temper);
@@ -95,12 +127,22 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
       gvdrvstime[iMod]->SetPoint(npt,timeday,vdrift);
       gvdrvstime[iMod]->SetPointError(npt,0,errSpeed[iMod]);
     }
+    Int_t npt=gGoodInjVsRun->GetN();
+    gGoodInjVsRun->SetPoint(npt,(Float_t)nrun,iGoodInj);
+    gAverSpeedVsRun->SetPoint(npt,(Float_t)nrun,iAverSpeed);
+       
+    npt=gGoodInjVsTime->GetN();
+    gGoodInjVsTime->SetPoint(npt,timeday,iGoodInj);
+    gAverSpeedVsTime->SetPoint(npt,timeday,iAverSpeed);
+       
+    printf("Number of half-modules with drift speed from injectors = %d\n",iGoodInj);
+    printf("Number of half-modules with average drift speed        = %d\n",iAverSpeed);
     f->Close();
   }
 
   Int_t mod1=244-240;
   Int_t mod2=277-240;
-//   Int_t mod1=268-240;
+  //  Int_t mod2=259-240;
 //   Int_t mod2=274-240;
   Int_t mod3=327-240;
    Int_t mod4=453-240;
@@ -120,10 +162,14 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
     gvdrvstime[iMod]->Write();
     gvdrvsrun[iMod]->Write();
   }
+  gGoodInjVsRun->Write();
+  gGoodInjVsTime->Write();
+  gAverSpeedVsRun->Write();
+  gAverSpeedVsTime->Write();
   ofil->Close();
 
   gStyle->SetOptTitle(0);
-  TCanvas* c0=new TCanvas();
+  TCanvas* c0=new TCanvas("c0","Vdrift vs. time");
   c0->SetGridx();
   c0->SetGridy();
   gvdrvstime[mod1]->SetMarkerStyle(20);
@@ -163,7 +209,7 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
   lent->SetTextColor(4);
   leg->Draw();
 
-  TCanvas* c1=new TCanvas();
+  TCanvas* c1=new TCanvas("c1","Vdrift vs. run");
   c1->SetGridx();
   c1->SetGridy();
   gvdrvsrun[mod1]->SetMarkerStyle(20);
@@ -252,7 +298,24 @@ void PlotDriftSpeedSDDVsTime(Int_t year=2010, Int_t firstRun=62840,
     lin->SetY2(i+0.5);
     lin->DrawClone();
   }
+
+  TCanvas* c4=new TCanvas("c4","GoodMod vs. run");
+  c4->SetGridx();
+  c4->SetGridy();
+  gGoodInjVsRun->SetMarkerStyle(20);
+  gGoodInjVsRun->SetMarkerColor(2);
+  gGoodInjVsRun->GetXaxis()->SetTitle("Run number");
+  gGoodInjVsRun->GetYaxis()->SetTitle("Half-modules with drift speed from injectors");
+  gGoodInjVsRun->Draw("AP");
   
+  TCanvas* c5=new TCanvas("c5","GoodMod vs. time");
+  c5->SetGridx();
+  c5->SetGridy();
+  gGoodInjVsTime->SetMarkerStyle(20);
+  gGoodInjVsTime->SetMarkerColor(2);
+  gGoodInjVsTime->GetXaxis()->SetTitle(title);
+  gGoodInjVsTime->GetYaxis()->SetTitle("Half-modules with drift speed from injectors");
+  gGoodInjVsTime->Draw("AP");
 }
 
 void FillErrors(Float_t errSpeed[260]){