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

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TF1.h>
#include <TFile.h>
#include <TLine.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TGraph.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TString.h>
#include <TGrid.h>
#include <TCanvas.h>
#include <TLatex.h>
#include <TObjArray.h>
#include "AliCDBEntry.h"
#include "AliITSDriftSpeedArraySDD.h"
#include "AliITSDriftSpeedSDD.h"
#include "AliITSgeomTGeo.h"
#endif

// Macro to plot the calibration parameters from the OCDB file 
// created from an INJECTOR run (OCDB/ITS/Calib/DriftSpeedSDD)
// Two methods ShowDriftSpeedSDD:
//  - the first takes the name of the file to be displayed
//  - the second builds the alien path+name from run number and file version
//
// Origin: F. Prino (prino@to.infn.it)

Bool_t kNoDraw = kFALSE; // set to kTRUE to eliminate module dependent plots


void ShowDriftSpeedSDD(Char_t filnam[150]="$ALICE_ROOT/ITS/Calib/DriftSpeedSDD/Run0_9999999_v0_s0.root", Int_t firstmod=0, Int_t lastmod=260,Int_t nrun=0){
  TFile *f= TFile::Open(filnam);
  AliCDBEntry *ent=(AliCDBEntry*)f->Get("AliCDBEntry");
  TObjArray *drspSDD = (TObjArray *)ent->GetObject();
  AliITSDriftSpeedArraySDD *vdriftarr0;
  AliITSDriftSpeedArraySDD *vdriftarr1;

  TH2F* hlay3=new TH2F("hlay3","Injector Status Layer 3",12,-0.5,5.5,14,-0.5,13.5);
  hlay3->GetXaxis()->SetTitle("Detector");
  hlay3->GetYaxis()->SetTitle("Ladder");
  hlay3->GetXaxis()->SetTickLength(0);
  hlay3->GetYaxis()->SetTickLength(0);
  hlay3->SetStats(0);
  hlay3->SetMinimum(-0.01);
  hlay3->SetMaximum(7.);
  TH2F* hlay4=new TH2F("hlay4","Injector Status Layer 4",16,-0.5,7.5,22,-0.5,21.5);
  hlay4->GetXaxis()->SetTitle("Detector");
  hlay4->GetYaxis()->SetTitle("Ladder");
  hlay4->GetXaxis()->SetTickLength(0);
  hlay4->GetYaxis()->SetTickLength(0);
  hlay4->GetYaxis()->SetTitle("Ladder");
  hlay4->SetStats(0);
  hlay4->SetMinimum(-0.01);
  hlay4->SetMaximum(7.);

  TGraph **gvdr0=new TGraph*[260];
  TGraph **gvdr1=new TGraph*[260];
  TCanvas *c0=0x0;
  if(!kNoDraw)c0=new TCanvas("c0","Module Drift Speed",700,1000);

  TGraph *tempvsmod0=new TGraph(0);
  TGraph *tempvsmod1=new TGraph(0);
  TGraph *vvsmod0=new TGraph(0);
  TGraph *vvsmod1=new TGraph(0);
  TGraph *poldegvsmod0=new TGraph(0); 
  TGraph *poldegvsmod1=new TGraph(0); 
  TGraph *anmaxvsmod0=new TGraph(0); 
  TGraph *anmaxvsmod1=new TGraph(0); 
  TGraph *dvcevsmod0=new TGraph(0);
  TGraph *dvcevsmod1=new TGraph(0);
  TGraph *dveevsmod0=new TGraph(0);
  TGraph *dveevsmod1=new TGraph(0);

  char tit0[100];
  sprintf(tit0,"Temperature vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Temperature vs. mod. number - Run %d",nrun);
  tempvsmod0->SetTitle(tit0);
  tempvsmod1->SetTitle(tit0);

  sprintf(tit0,"Drift Speed vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Drift Speed vs. mod. number - Run %d",nrun);
  vvsmod0->SetTitle(tit0);
  vvsmod1->SetTitle(tit0);

  sprintf(tit0,"Degree of poly fit vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Degree of poly fit vs. mod. number - Run %d",nrun);
  poldegvsmod0->SetTitle(tit0);
  poldegvsmod1->SetTitle(tit0);

  sprintf(tit0,"Anode with max. vdrift vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Anode with max. vdrift vs. mod. number - Run %d",nrun);
  anmaxvsmod0->SetTitle(tit0);
  anmaxvsmod1->SetTitle(tit0);

  sprintf(tit0,"Delta Vdrift 128-0 vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Delta Vdrift 128-0 vs. mod. number - Run %d",nrun);
  dvcevsmod0->SetTitle(tit0);
  dvcevsmod1->SetTitle(tit0);

  sprintf(tit0,"Delta Vdrift 256-0 vs. mod. number");
  if(nrun!=0)sprintf(tit0,"Delta Vdrift 256-0 vs. mod. number - Run %d",nrun);
  dveevsmod0->SetTitle(tit0);
  dveevsmod1->SetTitle(tit0);

  TF1* fPoly=new TF1("fPoly","pol3",0.,256.);
  Char_t tit[100];
  TString psnm0 = "vdriftSDD.ps[";
  TString psnm1 = "vdriftSDD.ps";
  TString psnm2 = "vdriftSDD.ps]";
  if(!kNoDraw) c0->Print(psnm0.Data());
  Int_t cntpad = 0;
  Int_t iGoodInj=0;
  Int_t iAverSpeed=0;
  TLatex* tleft=new TLatex(0.2,0.82,"Side 0");
  tleft->SetNDC();
  tleft->SetTextColor(1);
  TLatex* tright=new TLatex(0.2,0.75,"Side 1");
  tright->SetNDC();
  tright->SetTextColor(2);

  for(Int_t i=firstmod; i<lastmod; i++){
    Int_t iMod=i+240;
    Int_t lay,lad,det;
    AliITSgeomTGeo::GetModuleId(iMod,lay,lad,det);
    Int_t i0=2*i;
    Int_t i1=1+2*i;
    vdriftarr0=(AliITSDriftSpeedArraySDD*)drspSDD->At(i0);
    vdriftarr1=(AliITSDriftSpeedArraySDD*)drspSDD->At(i1);
    AliITSDriftSpeedSDD* vdrift0=0x0;
    if(vdriftarr0) vdrift0=vdriftarr0->GetDriftSpeedObject(0);
    AliITSDriftSpeedSDD* vdrift1=0x0;
    if(vdriftarr1) vdrift1=vdriftarr1->GetDriftSpeedObject(0);

    gvdr0[i]=new TGraph(0);
    gvdr1[i]=new TGraph(0);
    gvdr0[i]->SetMarkerStyle(7);
    gvdr1[i]->SetMarkerStyle(7);
    gvdr0[i]->SetMarkerColor(kBlack);
    gvdr0[i]->SetLineColor(kBlack);
    gvdr0[i]->SetMinimum(5.);
    gvdr1[i]->SetMinimum(5.);
    gvdr0[i]->SetMaximum(7.5);
    gvdr1[i]->SetMaximum(7.5);
    sprintf(tit,"Mod %d\n",iMod);
    gvdr0[i]->SetTitle(tit);
    gvdr1[i]->SetTitle(tit);
    gvdr1[i]->SetMarkerColor(kRed);
    gvdr1[i]->SetLineColor(kRed);

    for(Int_t iAn=0; iAn<256; iAn++){
      Float_t vel0=0;
      if(vdrift0) vel0=vdrift0->GetDriftSpeedAtAnode(iAn);
      Float_t vel1=0;
      if(vdrift1) vel1=vdrift1->GetDriftSpeedAtAnode(iAn);
      gvdr0[i]->SetPoint(iAn,(Float_t)iAn,vel0);
      gvdr1[i]->SetPoint(iAn,(Float_t)iAn,vel1);
    }
    Int_t statusInj0=vdriftarr0->GetInjectorStatus();
    Int_t statusInj1=vdriftarr1->GetInjectorStatus();
    if(statusInj0>0) iGoodInj++;
    else iAverSpeed++;
    if(statusInj1>0) iGoodInj++;
    else iAverSpeed++;

    printf(" Mod. %d \tStatusLR=%X %X \t v(an 128l)= %f",iMod,statusInj0,statusInj1,vdriftarr0->GetDriftSpeed(0,128));
    printf("        \t v(an 128r)= %f  Degree=%d %d\n",vdriftarr1->GetDriftSpeed(0,128),vdrift0->GetDegreeofPoly(),vdrift1->GetDegreeofPoly());

    Int_t n7=(statusInj0&(0x1F<<25))>>25;
    Int_t n6=(statusInj0&(0x1F<<20))>>20;
    Int_t n5=(statusInj0&(0x1F<<15))>>15;
    Int_t n4=(statusInj0&(0x1F<<5))>>10;
    Int_t n3=(statusInj0&(0x1F<<5))>>5;
    Int_t n2=statusInj0&0x1F;
    Float_t aveStatus0=(7.*n7+6.*n6+5.*n5+4.*n4+3.*n3+2.*n2)/32.;
    n7=(statusInj1&(0x1F<<25))>>25;
    n6=(statusInj1&(0x1F<<20))>>20;
    n5=(statusInj1&(0x1F<<15))>>15;
    n4=(statusInj1&(0x1F<<5))>>10;
    n3=(statusInj1&(0x1F<<5))>>5;
    n2=statusInj1&0x1F;
    Float_t aveStatus1=(7.*n7+6.*n6+5.*n5+4.*n4+3.*n3+2.*n2)/32.;

    Int_t index=1+(det-1)*2;
    if(lay==3){ 
      hlay3->SetBinContent(index,lad,aveStatus0);
      hlay3->SetBinContent(index+1,lad,aveStatus1);
    }
    if(lay==4){ 
      hlay4->SetBinContent(index,lad,aveStatus0);
      hlay4->SetBinContent(index+1,lad,aveStatus1);
    }


    if(!kNoDraw){
      if (i%12==0 ) {
	c0->cd();
	c0->Modified();
	c0->Update();
	if (i) c0->Print(psnm1.Data());
	c0->Clear();
	c0->Divide(3,4);
	cntpad = 0;
      }
      c0->cd(++cntpad);
      gvdr0[i]->Draw("AP");
      gvdr0[i]->GetXaxis()->SetTitle("Anode");
      gvdr0[i]->GetYaxis()->SetTitle("Vdrift (#mum/ns)");      
      gvdr1[i]->Draw("P same");
      gvdr1[i]->GetXaxis()->SetTitle("Anode");
      gvdr1[i]->GetYaxis()->SetTitle("Vdrift (#mum/ns)");
      tleft->Draw();
      tright->Draw();
    }

    Float_t vel0=0;
    Float_t pd0=0;
    if(vdrift0){ 
      vel0=vdrift0->GetDriftSpeedAtAnode(128);
     pd0=vdrift0->GetDegreeofPoly();
    }
    Float_t vel1=0;
    Float_t pd1=0;
    if(vdrift1){ 
      vel1=vdrift1->GetDriftSpeedAtAnode(128);
      pd1=vdrift1->GetDegreeofPoly();
    }
    Float_t Edrift=(1800-45)/291/0.012;  
    Float_t mob0=vel0*1.E5/Edrift;  
    Float_t temper0=293.15*TMath::Power((mob0/1350.),-1/2.4); 
    Float_t mob1=vel1*1.E5/Edrift;  
    Float_t temper1=293.15*TMath::Power((mob1/1350.),-1/2.4); 
    tempvsmod0->SetPoint(tempvsmod0->GetN(),(Float_t)iMod,temper0);
    tempvsmod1->SetPoint(tempvsmod1->GetN(),(Float_t)iMod,temper1);
    vvsmod0->SetPoint(vvsmod0->GetN(),(Float_t)iMod,vel0);
    vvsmod1->SetPoint(vvsmod1->GetN(),(Float_t)iMod,vel1);
    poldegvsmod0->SetPoint(poldegvsmod0->GetN(),(Float_t)iMod,pd0);
    poldegvsmod1->SetPoint(poldegvsmod1->GetN(),(Float_t)iMod,pd1);

    for(Int_t ipar=0; ipar<=vdrift0->GetDegreeofPoly(); ipar++){
      fPoly->SetParameter(ipar,vdrift0->GetDriftSpeedParameter(ipar));
    }
    if(vdrift0->GetDegreeofPoly()<3){
      for(Int_t ipar=vdrift0->GetDegreeofPoly()+1; ipar<=3; ipar++) fPoly->SetParameter(ipar,0.);
    }

    anmaxvsmod0->SetPoint(anmaxvsmod0->GetN(),(Float_t)iMod,fPoly->GetMaximumX(0.,256.));
    dvcevsmod0->SetPoint(dvcevsmod0->GetN(),(Float_t)iMod,fPoly->Eval(128)-fPoly->Eval(0));
    dveevsmod0->SetPoint(dveevsmod0->GetN(),(Float_t)iMod,fPoly->Eval(256)-fPoly->Eval(0));
    
    for(Int_t ipar=0; ipar<=vdrift1->GetDegreeofPoly(); ipar++){
      fPoly->SetParameter(ipar,vdrift1->GetDriftSpeedParameter(ipar));
    }
    if(vdrift1->GetDegreeofPoly()<3){
      for(Int_t ipar=vdrift1->GetDegreeofPoly()+1; ipar<=3; ipar++) fPoly->SetParameter(ipar,0.);
    }
    anmaxvsmod1->SetPoint(anmaxvsmod1->GetN(),(Float_t)iMod,fPoly->GetMaximumX(0.,256.));
    dvcevsmod1->SetPoint(dvcevsmod1->GetN(),(Float_t)iMod,fPoly->Eval(128)-fPoly->Eval(0));
    dveevsmod1->SetPoint(dveevsmod1->GetN(),(Float_t)iMod,fPoly->Eval(256)-fPoly->Eval(0));
    //    getchar();
  }

  if(!kNoDraw){
    c0->cd();
    c0->Modified();
    c0->Update();
    c0->Print(psnm2.Data());
  }
  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);

  gStyle->SetPalette(59);

  TCanvas* clay=new TCanvas("clay","Injector Status",900,600);
  clay->Divide(2,1);
  clay->cd(1);
  hlay3->Draw("colz");
  TLine** linv3=new TLine*[5];
  for(Int_t i=0;i<5;i++){
    linv3[i]=new TLine(i+0.5,-0.5,i+0.5,13.5);
    linv3[i]->SetLineColor(kGray+1);
    linv3[i]->Draw();
  }
  TLine** linh3=new TLine*[13];
  for(Int_t i=0;i<13;i++){
    linh3[i]=new TLine(-0.5,i+0.5,5.5,i+0.5);
    linh3[i]->SetLineColor(kGray+1);
    linh3[i]->Draw();
  }
  clay->cd(2);
  hlay4->Draw("colz");
  TLine** linv4=new TLine*[7];
  for(Int_t i=0;i<7;i++){
    linv4[i]=new TLine(i+0.5,-0.5,i+0.5,21.5);
    linv4[i]->SetLineColor(kGray+1);
    linv4[i]->Draw();
  }
  TLine** linh4=new TLine*[21];
  for(Int_t i=0;i<21;i++){
    linh4[i]=new TLine(-0.5,i+0.5,7.5,i+0.5);
    linh4[i]->SetLineColor(kGray+1);
    linh4[i]->Draw();
  }
  clay->Modified();


  TCanvas* c2;
  c2=new TCanvas("c2","Vdrift vs. mod",1000,700);
  vvsmod0->SetMarkerStyle(20);
  vvsmod0->Draw("AP");
  vvsmod0->GetXaxis()->SetTitle("Module Number");
  vvsmod0->GetYaxis()->SetTitle("Vdrift (#mum/ns)");
  vvsmod1->SetMarkerStyle(21);
  vvsmod1->SetMarkerColor(2);
  vvsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();

  TCanvas* c2t;
  c2t=new TCanvas("c2t","Temper vs. mod",1000,700);
  tempvsmod0->SetMarkerStyle(20);
  tempvsmod0->Draw("AP");
  tempvsmod0->GetXaxis()->SetTitle("Module Number");
  tempvsmod0->GetYaxis()->SetTitle("Temperature (K)");
  tempvsmod1->SetMarkerStyle(21);
  tempvsmod1->SetMarkerColor(2);
  tempvsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();

  TCanvas* c3;
  c3=new TCanvas("c3","Params vs. mod",900,900);
  c3->Divide(2,2);
  
  c3->cd(1);
  gPad->SetLeftMargin(0.14);
  poldegvsmod0->SetMarkerStyle(20);
  poldegvsmod0->Draw("AP");
  poldegvsmod0->GetXaxis()->SetTitle("Module Number");
  poldegvsmod0->GetYaxis()->SetTitle("Degree of Polynomial fit");
  poldegvsmod0->GetYaxis()->SetTitleOffset(1.4);
  poldegvsmod1->SetMarkerStyle(21);
  poldegvsmod1->SetMarkerColor(2);
  poldegvsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();
  c3->cd(2);
  gPad->SetLeftMargin(0.14);
  anmaxvsmod0->SetMarkerStyle(20);
  anmaxvsmod0->Draw("AP");
  anmaxvsmod0->GetXaxis()->SetTitle("Module Number");
  anmaxvsmod0->GetYaxis()->SetTitle("Anode with max. drift speed");
  anmaxvsmod0->GetYaxis()->SetTitleOffset(1.4);
  anmaxvsmod1->SetMarkerStyle(21);
  anmaxvsmod1->SetMarkerColor(2);
  anmaxvsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();
  c3->cd(3);
  gPad->SetLeftMargin(0.14);
  dvcevsmod0->SetMarkerStyle(20);
  dvcevsmod0->Draw("AP");
  dvcevsmod0->GetXaxis()->SetTitle("Module Number");
  dvcevsmod0->GetYaxis()->SetTitle("vdrift(anode128)-vdrift(anode0)");
  dvcevsmod0->GetYaxis()->SetTitleOffset(1.4);
  dvcevsmod1->SetMarkerStyle(21);
  dvcevsmod1->SetMarkerColor(2);
  dvcevsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();
  c3->cd(4);
  gPad->SetLeftMargin(0.14);
  dveevsmod0->SetMarkerStyle(20);
  dveevsmod0->Draw("AP");
  dveevsmod0->GetYaxis()->SetTitleOffset(1.4);
  dveevsmod0->GetXaxis()->SetTitle("Module Number");
  dveevsmod0->GetYaxis()->SetTitle("vdrift(anode256)-vdrift(anode0)");
  dveevsmod1->SetMarkerStyle(21);
  dveevsmod1->SetMarkerColor(2);
  dveevsmod1->Draw("SAMEP");
  tleft->Draw();
  tright->Draw();

  
}


void ShowDriftSpeedSDD(Int_t nrun, Int_t year=2010, Int_t nv=-1){
  TGrid::Connect("alien:",0,0,"t");
  TString cmd=Form("gbbox find \"/alice/data/%d/OCDB/ITS/Calib/DriftSpeedSDD\" \"Run%d*.root\" > run.txt",year,nrun);
  if(nv>0){
    cmd.Form("gbbox find \"/alice/data/%d/OCDB/ITS/Calib/DriftSpeedSDD\" \"Run%d_999999999_v%d_s0.root\" > run.txt",year,nrun,nv);  
  }
  gSystem->Exec(cmd.Data());
  
  Char_t filnam[200],filnamalien[200];
  FILE* runtxt=fopen("run.txt","r");
  fscanf(runtxt,"%s\n",filnam);    
  if(!strstr(filnam,"/alice/data/")){
    printf("Bad run number\n");
    gSystem->Exec("rm run.txt");
    return;
  }  
  sprintf(filnamalien,"alien://%s",filnam);
  
  printf("Open file: %s\n",filnamalien);
  ShowDriftSpeedSDD(filnamalien,0,260,nrun);
  fclose(runtxt);
  gSystem->Exec("rm run.txt");
  
}
Commit	Line	Data
40b3a8c7	1	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
2420973e	2	#include <TF1.h>
40b3a8c7	3	#include <TFile.h>
aa5cbbc5	4	#include <TLine.h>
40b3a8c7	5	#include <TH1F.h>
aa5cbbc5	6	#include <TH2F.h>
40b3a8c7	7	#include <TGraph.h>
40b3a8c7	8	#include <TStyle.h>
551eb942	9	#include <TSystem.h>
551eb942	10	#include <TString.h>
40b3a8c7	11	#include <TGrid.h>
	12	#include <TCanvas.h>
	13	#include <TLatex.h>
	14	#include <TObjArray.h>
	15	#include "AliCDBEntry.h"
	16	#include "AliITSDriftSpeedArraySDD.h"
2420973e	17	#include "AliITSDriftSpeedSDD.h"
aa5cbbc5	18	#include "AliITSgeomTGeo.h"
40b3a8c7	19	#endif
	20
	21	// Macro to plot the calibration parameters from the OCDB file
	22	// created from an INJECTOR run (OCDB/ITS/Calib/DriftSpeedSDD)
	23	// Two methods ShowDriftSpeedSDD:
	24	// - the first takes the name of the file to be displayed
	25	// - the second builds the alien path+name from run number and file version
	26	//
	27	// Origin: F. Prino (prino@to.infn.it)
	28
b70dd7af	29	Bool_t kNoDraw = kFALSE; // set to kTRUE to eliminate module dependent plots
b70dd7af	30
e9b71015	31
b70dd7af	32	void ShowDriftSpeedSDD(Char_t filnam[150]="$ALICE_ROOT/ITS/Calib/DriftSpeedSDD/Run0_9999999_v0_s0.root", Int_t firstmod=0, Int_t lastmod=260,Int_t nrun=0){
40b3a8c7	33	TFile *f= TFile::Open(filnam);
	34	AliCDBEntry ent=(AliCDBEntry)f->Get("AliCDBEntry");
	35	TObjArray drspSDD = (TObjArray )ent->GetObject();
	36	AliITSDriftSpeedArraySDD *vdriftarr0;
	37	AliITSDriftSpeedArraySDD *vdriftarr1;
aa5cbbc5	38
	39	TH2F* hlay3=new TH2F("hlay3","Injector Status Layer 3",12,-0.5,5.5,14,-0.5,13.5);
	40	hlay3->GetXaxis()->SetTitle("Detector");
	41	hlay3->GetYaxis()->SetTitle("Ladder");
	42	hlay3->GetXaxis()->SetTickLength(0);
	43	hlay3->GetYaxis()->SetTickLength(0);
	44	hlay3->SetStats(0);
	45	hlay3->SetMinimum(-0.01);
	46	hlay3->SetMaximum(7.);
	47	TH2F* hlay4=new TH2F("hlay4","Injector Status Layer 4",16,-0.5,7.5,22,-0.5,21.5);
	48	hlay4->GetXaxis()->SetTitle("Detector");
	49	hlay4->GetYaxis()->SetTitle("Ladder");
	50	hlay4->GetXaxis()->SetTickLength(0);
	51	hlay4->GetYaxis()->SetTickLength(0);
	52	hlay4->GetYaxis()->SetTitle("Ladder");
	53	hlay4->SetStats(0);
	54	hlay4->SetMinimum(-0.01);
	55	hlay4->SetMaximum(7.);
	56
40b3a8c7	57	TGraph *gvdr0=new TGraph[260];
40b3a8c7	58	TGraph *gvdr1=new TGraph[260];
e9b71015	59	TCanvas *c0=0x0;
e9b71015	60	if(!kNoDraw)c0=new TCanvas("c0","Module Drift Speed",700,1000);
40b3a8c7	61
813d0163	62	TGraph *tempvsmod0=new TGraph(0);
813d0163	63	TGraph *tempvsmod1=new TGraph(0);
40b3a8c7	64	TGraph *vvsmod0=new TGraph(0);
40b3a8c7	65	TGraph *vvsmod1=new TGraph(0);
915c0469	66	TGraph *poldegvsmod0=new TGraph(0);
	67	TGraph *poldegvsmod1=new TGraph(0);
	68	TGraph *anmaxvsmod0=new TGraph(0);
	69	TGraph *anmaxvsmod1=new TGraph(0);
	70	TGraph *dvcevsmod0=new TGraph(0);
	71	TGraph *dvcevsmod1=new TGraph(0);
	72	TGraph *dveevsmod0=new TGraph(0);
	73	TGraph *dveevsmod1=new TGraph(0);
	74
b70dd7af	75	char tit0[100];
813d0163	76	sprintf(tit0,"Temperature vs. mod. number");
	77	if(nrun!=0)sprintf(tit0,"Temperature vs. mod. number - Run %d",nrun);
	78	tempvsmod0->SetTitle(tit0);
	79	tempvsmod1->SetTitle(tit0);
	80
b70dd7af	81	sprintf(tit0,"Drift Speed vs. mod. number");
	82	if(nrun!=0)sprintf(tit0,"Drift Speed vs. mod. number - Run %d",nrun);
	83	vvsmod0->SetTitle(tit0);
	84	vvsmod1->SetTitle(tit0);
915c0469	85
	86	sprintf(tit0,"Degree of poly fit vs. mod. number");
	87	if(nrun!=0)sprintf(tit0,"Degree of poly fit vs. mod. number - Run %d",nrun);
	88	poldegvsmod0->SetTitle(tit0);
	89	poldegvsmod1->SetTitle(tit0);
	90
	91	sprintf(tit0,"Anode with max. vdrift vs. mod. number");
	92	if(nrun!=0)sprintf(tit0,"Anode with max. vdrift vs. mod. number - Run %d",nrun);
	93	anmaxvsmod0->SetTitle(tit0);
	94	anmaxvsmod1->SetTitle(tit0);
	95
	96	sprintf(tit0,"Delta Vdrift 128-0 vs. mod. number");
	97	if(nrun!=0)sprintf(tit0,"Delta Vdrift 128-0 vs. mod. number - Run %d",nrun);
	98	dvcevsmod0->SetTitle(tit0);
	99	dvcevsmod1->SetTitle(tit0);
	100
	101	sprintf(tit0,"Delta Vdrift 256-0 vs. mod. number");
	102	if(nrun!=0)sprintf(tit0,"Delta Vdrift 256-0 vs. mod. number - Run %d",nrun);
	103	dveevsmod0->SetTitle(tit0);
	104	dveevsmod1->SetTitle(tit0);
	105
49c59b1d	106	TF1* fPoly=new TF1("fPoly","pol3",0.,256.);
40b3a8c7	107	Char_t tit[100];
e9b71015	108	TString psnm0 = "vdriftSDD.ps[";
	109	TString psnm1 = "vdriftSDD.ps";
	110	TString psnm2 = "vdriftSDD.ps]";
	111	if(!kNoDraw) c0->Print(psnm0.Data());
	112	Int_t cntpad = 0;
915c0469	113	Int_t iGoodInj=0;
915c0469	114	Int_t iAverSpeed=0;
e9b71015	115	TLatex* tleft=new TLatex(0.2,0.82,"Side 0");
	116	tleft->SetNDC();
	117	tleft->SetTextColor(1);
	118	TLatex* tright=new TLatex(0.2,0.75,"Side 1");
	119	tright->SetNDC();
	120	tright->SetTextColor(2);
	121
40b3a8c7	122	for(Int_t i=firstmod; i<lastmod; i++){
266a4b47	123	Int_t iMod=i+240;
aa5cbbc5	124	Int_t lay,lad,det;
aa5cbbc5	125	AliITSgeomTGeo::GetModuleId(iMod,lay,lad,det);
40b3a8c7	126	Int_t i0=2*i;
	127	Int_t i1=1+2*i;
	128	vdriftarr0=(AliITSDriftSpeedArraySDD*)drspSDD->At(i0);
	129	vdriftarr1=(AliITSDriftSpeedArraySDD*)drspSDD->At(i1);
915c0469	130	AliITSDriftSpeedSDD* vdrift0=0x0;
	131	if(vdriftarr0) vdrift0=vdriftarr0->GetDriftSpeedObject(0);
	132	AliITSDriftSpeedSDD* vdrift1=0x0;
	133	if(vdriftarr1) vdrift1=vdriftarr1->GetDriftSpeedObject(0);
	134
40b3a8c7	135	gvdr0[i]=new TGraph(0);
	136	gvdr1[i]=new TGraph(0);
	137	gvdr0[i]->SetMarkerStyle(7);
	138	gvdr1[i]->SetMarkerStyle(7);
e9b71015	139	gvdr0[i]->SetMarkerColor(kBlack);
e9b71015	140	gvdr0[i]->SetLineColor(kBlack);
40b3a8c7	141	gvdr0[i]->SetMinimum(5.);
40b3a8c7	142	gvdr1[i]->SetMinimum(5.);
e9b71015	143	gvdr0[i]->SetMaximum(7.5);
e9b71015	144	gvdr1[i]->SetMaximum(7.5);
266a4b47	145	sprintf(tit,"Mod %d\n",iMod);
40b3a8c7	146	gvdr0[i]->SetTitle(tit);
40b3a8c7	147	gvdr1[i]->SetTitle(tit);
e9b71015	148	gvdr1[i]->SetMarkerColor(kRed);
e9b71015	149	gvdr1[i]->SetLineColor(kRed);
40b3a8c7	150
	151	for(Int_t iAn=0; iAn<256; iAn++){
	152	Float_t vel0=0;
915c0469	153	if(vdrift0) vel0=vdrift0->GetDriftSpeedAtAnode(iAn);
40b3a8c7	154	Float_t vel1=0;
915c0469	155	if(vdrift1) vel1=vdrift1->GetDriftSpeedAtAnode(iAn);
40b3a8c7	156	gvdr0[i]->SetPoint(iAn,(Float_t)iAn,vel0);
	157	gvdr1[i]->SetPoint(iAn,(Float_t)iAn,vel1);
	158	}
aa5cbbc5	159	Int_t statusInj0=vdriftarr0->GetInjectorStatus();
	160	Int_t statusInj1=vdriftarr1->GetInjectorStatus();
	161	if(statusInj0>0) iGoodInj++;
915c0469	162	else iAverSpeed++;
aa5cbbc5	163	if(statusInj1>0) iGoodInj++;
915c0469	164	else iAverSpeed++;
915c0469	165
aa5cbbc5	166	printf(" Mod. %d \tStatusLR=%X %X \t v(an 128l)= %f",iMod,statusInj0,statusInj1,vdriftarr0->GetDriftSpeed(0,128));
915c0469	167	printf(" \t v(an 128r)= %f Degree=%d %d\n",vdriftarr1->GetDriftSpeed(0,128),vdrift0->GetDegreeofPoly(),vdrift1->GetDegreeofPoly());
2420973e	168
aa5cbbc5	169	Int_t n7=(statusInj0&(0x1F<<25))>>25;
	170	Int_t n6=(statusInj0&(0x1F<<20))>>20;
	171	Int_t n5=(statusInj0&(0x1F<<15))>>15;
	172	Int_t n4=(statusInj0&(0x1F<<5))>>10;
	173	Int_t n3=(statusInj0&(0x1F<<5))>>5;
	174	Int_t n2=statusInj0&0x1F;
	175	Float_t aveStatus0=(7.n7+6.n6+5.n5+4.n4+3.n3+2.n2)/32.;
	176	n7=(statusInj1&(0x1F<<25))>>25;
	177	n6=(statusInj1&(0x1F<<20))>>20;
	178	n5=(statusInj1&(0x1F<<15))>>15;
	179	n4=(statusInj1&(0x1F<<5))>>10;
	180	n3=(statusInj1&(0x1F<<5))>>5;
	181	n2=statusInj1&0x1F;
	182	Float_t aveStatus1=(7.n7+6.n6+5.n5+4.n4+3.n3+2.n2)/32.;
	183
	184	Int_t index=1+(det-1)*2;
	185	if(lay==3){
	186	hlay3->SetBinContent(index,lad,aveStatus0);
	187	hlay3->SetBinContent(index+1,lad,aveStatus1);
	188	}
	189	if(lay==4){
	190	hlay4->SetBinContent(index,lad,aveStatus0);
	191	hlay4->SetBinContent(index+1,lad,aveStatus1);
	192	}
	193
	194
2420973e	195	if(!kNoDraw){
e9b71015	196	if (i%12==0 ) {
	197	c0->cd();
	198	c0->Modified();
	199	c0->Update();
	200	if (i) c0->Print(psnm1.Data());
	201	c0->Clear();
	202	c0->Divide(3,4);
	203	cntpad = 0;
	204	}
	205	c0->cd(++cntpad);
2420973e	206	gvdr0[i]->Draw("AP");
2420973e	207	gvdr0[i]->GetXaxis()->SetTitle("Anode");
e9b71015	208	gvdr0[i]->GetYaxis()->SetTitle("Vdrift (#mum/ns)");
e9b71015	209	gvdr1[i]->Draw("P same");
2420973e	210	gvdr1[i]->GetXaxis()->SetTitle("Anode");
2420973e	211	gvdr1[i]->GetYaxis()->SetTitle("Vdrift (#mum/ns)");
e9b71015	212	tleft->Draw();
e9b71015	213	tright->Draw();
2420973e	214	}
2420973e	215
40b3a8c7	216	Float_t vel0=0;
915c0469	217	Float_t pd0=0;
	218	if(vdrift0){
	219	vel0=vdrift0->GetDriftSpeedAtAnode(128);
e9b71015	220	pd0=vdrift0->GetDegreeofPoly();
915c0469	221	}
40b3a8c7	222	Float_t vel1=0;
915c0469	223	Float_t pd1=0;
	224	if(vdrift1){
	225	vel1=vdrift1->GetDriftSpeedAtAnode(128);
	226	pd1=vdrift1->GetDegreeofPoly();
	227	}
813d0163	228	Float_t Edrift=(1800-45)/291/0.012;
	229	Float_t mob0=vel0*1.E5/Edrift;
	230	Float_t temper0=293.15*TMath::Power((mob0/1350.),-1/2.4);
	231	Float_t mob1=vel1*1.E5/Edrift;
	232	Float_t temper1=293.15*TMath::Power((mob1/1350.),-1/2.4);
	233	tempvsmod0->SetPoint(tempvsmod0->GetN(),(Float_t)iMod,temper0);
	234	tempvsmod1->SetPoint(tempvsmod1->GetN(),(Float_t)iMod,temper1);
266a4b47	235	vvsmod0->SetPoint(vvsmod0->GetN(),(Float_t)iMod,vel0);
266a4b47	236	vvsmod1->SetPoint(vvsmod1->GetN(),(Float_t)iMod,vel1);
915c0469	237	poldegvsmod0->SetPoint(poldegvsmod0->GetN(),(Float_t)iMod,pd0);
	238	poldegvsmod1->SetPoint(poldegvsmod1->GetN(),(Float_t)iMod,pd1);
	239
	240	for(Int_t ipar=0; ipar<=vdrift0->GetDegreeofPoly(); ipar++){
	241	fPoly->SetParameter(ipar,vdrift0->GetDriftSpeedParameter(ipar));
	242	}
	243	if(vdrift0->GetDegreeofPoly()<3){
	244	for(Int_t ipar=vdrift0->GetDegreeofPoly()+1; ipar<=3; ipar++) fPoly->SetParameter(ipar,0.);
	245	}
	246
	247	anmaxvsmod0->SetPoint(anmaxvsmod0->GetN(),(Float_t)iMod,fPoly->GetMaximumX(0.,256.));
	248	dvcevsmod0->SetPoint(dvcevsmod0->GetN(),(Float_t)iMod,fPoly->Eval(128)-fPoly->Eval(0));
	249	dveevsmod0->SetPoint(dveevsmod0->GetN(),(Float_t)iMod,fPoly->Eval(256)-fPoly->Eval(0));
40b3a8c7	250
915c0469	251	for(Int_t ipar=0; ipar<=vdrift1->GetDegreeofPoly(); ipar++){
	252	fPoly->SetParameter(ipar,vdrift1->GetDriftSpeedParameter(ipar));
	253	}
	254	if(vdrift1->GetDegreeofPoly()<3){
	255	for(Int_t ipar=vdrift1->GetDegreeofPoly()+1; ipar<=3; ipar++) fPoly->SetParameter(ipar,0.);
	256	}
	257	anmaxvsmod1->SetPoint(anmaxvsmod1->GetN(),(Float_t)iMod,fPoly->GetMaximumX(0.,256.));
	258	dvcevsmod1->SetPoint(dvcevsmod1->GetN(),(Float_t)iMod,fPoly->Eval(128)-fPoly->Eval(0));
	259	dveevsmod1->SetPoint(dveevsmod1->GetN(),(Float_t)iMod,fPoly->Eval(256)-fPoly->Eval(0));
40b3a8c7	260	// getchar();
40b3a8c7	261	}
915c0469	262
e9b71015	263	if(!kNoDraw){
	264	c0->cd();
	265	c0->Modified();
	266	c0->Update();
	267	c0->Print(psnm2.Data());
	268	}
915c0469	269	printf("Number of half-modules with drift speed from injectors = %d\n",iGoodInj);
	270	printf("Number of half-modules with average drift speed = %d\n",iAverSpeed);
	271
aa5cbbc5	272	gStyle->SetPalette(59);
	273
	274	TCanvas* clay=new TCanvas("clay","Injector Status",900,600);
	275	clay->Divide(2,1);
	276	clay->cd(1);
	277	hlay3->Draw("colz");
	278	TLine** linv3=new TLine*[5];
	279	for(Int_t i=0;i<5;i++){
	280	linv3[i]=new TLine(i+0.5,-0.5,i+0.5,13.5);
	281	linv3[i]->SetLineColor(kGray+1);
	282	linv3[i]->Draw();
	283	}
	284	TLine** linh3=new TLine*[13];
	285	for(Int_t i=0;i<13;i++){
	286	linh3[i]=new TLine(-0.5,i+0.5,5.5,i+0.5);
	287	linh3[i]->SetLineColor(kGray+1);
	288	linh3[i]->Draw();
	289	}
	290	clay->cd(2);
	291	hlay4->Draw("colz");
	292	TLine** linv4=new TLine*[7];
	293	for(Int_t i=0;i<7;i++){
	294	linv4[i]=new TLine(i+0.5,-0.5,i+0.5,21.5);
	295	linv4[i]->SetLineColor(kGray+1);
	296	linv4[i]->Draw();
	297	}
	298	TLine** linh4=new TLine*[21];
	299	for(Int_t i=0;i<21;i++){
	300	linh4[i]=new TLine(-0.5,i+0.5,7.5,i+0.5);
	301	linh4[i]->SetLineColor(kGray+1);
	302	linh4[i]->Draw();
	303	}
	304	clay->Modified();
	305
	306
40b3a8c7	307	TCanvas* c2;
aa5cbbc5	308	c2=new TCanvas("c2","Vdrift vs. mod",1000,700);
40b3a8c7	309	vvsmod0->SetMarkerStyle(20);
	310	vvsmod0->Draw("AP");
	311	vvsmod0->GetXaxis()->SetTitle("Module Number");
	312	vvsmod0->GetYaxis()->SetTitle("Vdrift (#mum/ns)");
	313	vvsmod1->SetMarkerStyle(21);
	314	vvsmod1->SetMarkerColor(2);
	315	vvsmod1->Draw("SAMEP");
40b3a8c7	316	tleft->Draw();
40b3a8c7	317	tright->Draw();
40b3a8c7	318
813d0163	319	TCanvas* c2t;
	320	c2t=new TCanvas("c2t","Temper vs. mod",1000,700);
	321	tempvsmod0->SetMarkerStyle(20);
	322	tempvsmod0->Draw("AP");
	323	tempvsmod0->GetXaxis()->SetTitle("Module Number");
	324	tempvsmod0->GetYaxis()->SetTitle("Temperature (K)");
	325	tempvsmod1->SetMarkerStyle(21);
	326	tempvsmod1->SetMarkerColor(2);
	327	tempvsmod1->Draw("SAMEP");
	328	tleft->Draw();
	329	tright->Draw();
	330
915c0469	331	TCanvas* c3;
aa5cbbc5	332	c3=new TCanvas("c3","Params vs. mod",900,900);
915c0469	333	c3->Divide(2,2);
	334
	335	c3->cd(1);
	336	gPad->SetLeftMargin(0.14);
	337	poldegvsmod0->SetMarkerStyle(20);
	338	poldegvsmod0->Draw("AP");
	339	poldegvsmod0->GetXaxis()->SetTitle("Module Number");
	340	poldegvsmod0->GetYaxis()->SetTitle("Degree of Polynomial fit");
	341	poldegvsmod0->GetYaxis()->SetTitleOffset(1.4);
	342	poldegvsmod1->SetMarkerStyle(21);
	343	poldegvsmod1->SetMarkerColor(2);
	344	poldegvsmod1->Draw("SAMEP");
	345	tleft->Draw();
	346	tright->Draw();
	347	c3->cd(2);
	348	gPad->SetLeftMargin(0.14);
	349	anmaxvsmod0->SetMarkerStyle(20);
	350	anmaxvsmod0->Draw("AP");
	351	anmaxvsmod0->GetXaxis()->SetTitle("Module Number");
	352	anmaxvsmod0->GetYaxis()->SetTitle("Anode with max. drift speed");
	353	anmaxvsmod0->GetYaxis()->SetTitleOffset(1.4);
	354	anmaxvsmod1->SetMarkerStyle(21);
	355	anmaxvsmod1->SetMarkerColor(2);
	356	anmaxvsmod1->Draw("SAMEP");
	357	tleft->Draw();
	358	tright->Draw();
	359	c3->cd(3);
	360	gPad->SetLeftMargin(0.14);
	361	dvcevsmod0->SetMarkerStyle(20);
	362	dvcevsmod0->Draw("AP");
	363	dvcevsmod0->GetXaxis()->SetTitle("Module Number");
	364	dvcevsmod0->GetYaxis()->SetTitle("vdrift(anode128)-vdrift(anode0)");
	365	dvcevsmod0->GetYaxis()->SetTitleOffset(1.4);
	366	dvcevsmod1->SetMarkerStyle(21);
	367	dvcevsmod1->SetMarkerColor(2);
	368	dvcevsmod1->Draw("SAMEP");
	369	tleft->Draw();
	370	tright->Draw();
	371	c3->cd(4);
	372	gPad->SetLeftMargin(0.14);
	373	dveevsmod0->SetMarkerStyle(20);
	374	dveevsmod0->Draw("AP");
	375	dveevsmod0->GetYaxis()->SetTitleOffset(1.4);
	376	dveevsmod0->GetXaxis()->SetTitle("Module Number");
	377	dveevsmod0->GetYaxis()->SetTitle("vdrift(anode256)-vdrift(anode0)");
	378	dveevsmod1->SetMarkerStyle(21);
	379	dveevsmod1->SetMarkerColor(2);
	380	dveevsmod1->Draw("SAMEP");
	381	tleft->Draw();
	382	tright->Draw();
	383
	384
40b3a8c7	385	}
	386
	387
	388
e9b71015	389	void ShowDriftSpeedSDD(Int_t nrun, Int_t year=2010, Int_t nv=-1){
40b3a8c7	390	TGrid::Connect("alien:",0,0,"t");
551eb942	391	TString cmd=Form("gbbox find \"/alice/data/%d/OCDB/ITS/Calib/DriftSpeedSDD\" \"Run%d*.root\" > run.txt",year,nrun);
e9b71015	392	if(nv>0){
	393	cmd.Form("gbbox find \"/alice/data/%d/OCDB/ITS/Calib/DriftSpeedSDD\" \"Run%d_999999999_v%d_s0.root\" > run.txt",year,nrun,nv);
	394	}
551eb942	395	gSystem->Exec(cmd.Data());
	396
	397	Char_t filnam[200],filnamalien[200];
	398	FILE* runtxt=fopen("run.txt","r");
	399	fscanf(runtxt,"%s\n",filnam);
	400	if(!strstr(filnam,"/alice/data/")){
	401	printf("Bad run number\n");
	402	gSystem->Exec("rm run.txt");
	403	return;
	404	}
	405	sprintf(filnamalien,"alien://%s",filnam);
	406
	407	printf("Open file: %s\n",filnamalien);
b70dd7af	408	ShowDriftSpeedSDD(filnamalien,0,260,nrun);
551eb942	409	fclose(runtxt);
	410	gSystem->Exec("rm run.txt");
	411
40b3a8c7	412	}