ITS/AnalyzeSDDInjectorsAllMod.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include <TH2F.h>
   3 #include <TCanvas.h>
   4 #include <TGraphErrors.h>
   5 #include <TStopwatch.h>
   6 #include <TStyle.h>
   7 #include <TLatex.h>
   8 #include <TFile.h>
   9 #include <TMath.h>
  10 #include <TGrid.h>
  11 #include <TF1.h>
  12 #include <TLine.h>
  13 #include "AliRawReader.h"
  14 #include "AliRawReaderDate.h"
  15 #include "AliRawReaderRoot.h"
  16 #include "AliITSOnlineSDDInjectors.h"
  17 #include "AliITSRawStreamSDD.h"
  18 #include "AliITSDDLModuleMapSDD.h"
  19 #endif
  20
  21 // Macro for the analysis of PULSER runs (equivalent to ITSSDDINJda.cxx)
  22 // Two functions named AnalyzeSDDInjectorsAllModules:
  23 // The first is for analyzing a local raw data file and takes as agrument the file name.
  24 // The second is for running on ALIEN
  25 // All DDLs are analyzed, the argument nDDL selects the DDL to be plotted
  26 // Origin: F. Prino (prino@to.infn.it)
  27
  28
  29 void AnalyzeSDDInjectorsAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=15,Int_t jpad=16, Int_t statuscut=7){
  30
  31   const Int_t kTotDDL=24;
  32   const Int_t kModPerDDL=12;
  33   const Int_t kSides=2;
  34
  35   AliITSDDLModuleMapSDD* dmap=new AliITSDDLModuleMapSDD();
  36   dmap->SetJun08Map();
  37
  38   TH2F** histo = new TH2F*[kTotDDL*kModPerDDL*kSides];
  39   Int_t nWrittenEv[kTotDDL*kModPerDDL*kSides];
  40   TGraphErrors** gvel = new TGraphErrors*[kTotDDL*kModPerDDL*kSides];
  41   AliITSOnlineSDDInjectors **anal=new AliITSOnlineSDDInjectors*[kTotDDL*kModPerDDL*kSides];
  42   TH1F** hvdriftl=new TH1F*[260];
  43   TH1F** hvdriftr=new TH1F*[260];
  44   Char_t hisnam[20];
  45   for(Int_t idet=0; idet<260;idet++){
  46     sprintf(hisnam,"vdriftl%03d",idet);
  47     hvdriftl[idet]=new TH1F(hisnam,"",500,5.5,8.0);
  48     sprintf(hisnam,"vdriftr%03d",idet);
  49     hvdriftr[idet]=new TH1F(hisnam,"",500,5.5,8.0);
  50   }
  51   for(Int_t iddl=0; iddl<kTotDDL;iddl++){
  52     for(Int_t imod=0; imod<kModPerDDL;imod++){
  53       for(Int_t isid=0;isid<kSides;isid++){
  54         Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
  55         sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
  56         histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
  57         anal[index]=new AliITSOnlineSDDInjectors(iddl,imod,isid);
  58 /* Uncomment these lines for analysis of runs with 40 MHz sapling */
  59 //      anal[index]->SetInjLineRange(0,20,50);
  60 //      anal[index]->SetInjLineRange(1,90,160);
  61 //      anal[index]->SetInjLineRange(2,170,240);
  62 //      anal[index]->SetTimeStep(25.);
  63 /* END of lines to be uncommented */
  64         nWrittenEv[index]=0;
  65       }
  66     }
  67   }
  68   TGraphErrors *gvvsmod0=new TGraphErrors(0);
  69   TGraphErrors *gvvsmod1=new TGraphErrors(0);
  70   TGraphErrors *gtvsmod0=new TGraphErrors(0);
  71   TGraphErrors *gtvsmod1=new TGraphErrors(0);
  72   Float_t gvmin=6.0, gvmax=7.5;
  73   Float_t gtmin=288., gtmax=308.;
  74   TH1F* hanst=new TH1F("hanst","",8,-0.5,7.5);
  75   TH1F* hpad7l=new TH1F("hpad7l","",33,-0.5,32.5);
  76   TH1F* hpad7r=new TH1F("hpad7r","",33,-0.5,32.5);
  77
  78   TCanvas* c0 = new TCanvas("c0","Event display",900,900);
  79   gStyle->SetPalette(1);
  80   TCanvas* c1 = new TCanvas("c1","Drift Speed vs. anode",900,900);
  81   Char_t text[50];
  82
  83   Int_t iev=firstEv;
  84   AliRawReader *rd;
  85   if(strstr(datafil,".root")!=0){
  86     rd=new AliRawReaderRoot(datafil,iev);
  87   }else{
  88     rd=new AliRawReaderDate(datafil,iev);
  89   }
  90
  91   Char_t gname[15];
  92   TF1 *funz=new TF1("funz","[0]+[1]*x+[2]*x*x+[3]*x*x*x",0.,255.);
  93   TLatex *t0=new TLatex();
  94   t0->SetNDC();
  95   t0->SetTextSize(0.06);
  96   t0->SetTextColor(4);
  97   Int_t readEv=0;
  98   do{
  99     c0->Clear();
 100     c0->Divide(4,6,0.001,0.001);
 101     c1->Clear();
 102     c1->Divide(4,6,0.001,0.001);
 103     printf("Event # %d\n",iev);
 104     UInt_t timeSt=rd->GetTimestamp();
 105     rd->Reset();
 106     for(Int_t iddl=0; iddl<kTotDDL;iddl++){
 107       for(Int_t imod=0; imod<kModPerDDL;imod++){
 108         for(Int_t isid=0;isid<kSides;isid++){
 109           Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
 110           histo[index]->Reset();
 111         }
 112       }
 113     }
 114
 115     AliITSRawStreamSDD s(rd);
 116     while(s.Next()){
 117       Int_t iDDL=rd->GetDDLID();
 118       Int_t iCarlos=s.GetCarlosId();
 119       if(iDDL>=0 && iDDL<kTotDDL && s.IsCompletedModule()==kFALSE){
 120         Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s.GetChannel();
 121         histo[index]->Fill(s.GetCoord2(),s.GetCoord1(),s.GetSignal());
 122       }
 123     }
 124
 125     for(Int_t iddl=0; iddl<kTotDDL;iddl++){
 126       for(Int_t imod=0; imod<kModPerDDL;imod++){
 127         for(Int_t isid=0;isid<kSides;isid++){
 128           Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
 129           anal[index]->AnalyzeEvent(histo[index]);
 130           anal[index]->WriteToASCII(iev,timeSt,nWrittenEv[index]);
 131           nWrittenEv[index]++;
 132           Int_t iMod=dmap->GetModuleNumber(iddl,imod);
 133           if(iMod!=-1){
 134             for(Int_t ipad=0;ipad<33;ipad++){
 135               Int_t st=anal[index]->GetInjPadStatus(ipad);
 136               hanst->Fill(st);
 137               if(anal[index]->GetInjPadStatus(ipad)>=statuscut){
 138                 if(isid==0) hpad7l->Fill(ipad);
 139                 if(isid==1) hpad7r->Fill(ipad);
 140               }
 141             }
 142             if(anal[index]->GetInjPadStatus(jpad)>=statuscut){
 143               Float_t vel=anal[index]->GetDriftSpeed(jpad);
 144               if(isid==0) hvdriftl[iMod-240]->Fill(vel);
 145               if(isid==1) hvdriftr[iMod-240]->Fill(vel);
 146             }
 147           }
 148           if(iddl==nDDL){
 149             Int_t index2=kSides*imod+isid;
 150             c0->cd(index2+1);
 151             histo[index]->SetMaximum(100.);
 152             histo[index]->DrawCopy("colz");
 153             sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid);
 154             t0->DrawLatex(0.15,0.92,text);
 155             c0->Update();
 156             c1->cd(index2+1);
 157             gvel[index]=anal[index]->GetDriftSpeedGraph();
 158             gvel[index]->SetMarkerStyle(20);
 159             gvel[index]->SetTitle("");
 160             sprintf(gname,"gvel%dev%d",index,iev);
 161             gvel[index]->SetName(gname);
 162             //      gvel[index]->SetMinimum(0);
 163             //gvel[index]->SetMaximum(200);
 164
 165             gvel[index]->GetXaxis()->SetLimits(0,256);
 166             gvel[index]->GetXaxis()->SetTitle("Anode");
 167             gvel[index]->GetYaxis()->SetTitle("Drift vel.");
 168             gvel[index]->GetXaxis()->SetTitleSize(0.07);
 169             gvel[index]->GetYaxis()->SetTitleSize(0.07);
 170             gvel[index]->GetXaxis()->SetTitleOffset(0.6);
 171             gvel[index]->GetYaxis()->SetTitleOffset(0.6);
 172             if(gvel[index]->GetN()>0) gvel[index]->Draw("AP");
 173             Float_t *param=anal[index]->GetDriftSpeedFitParam();
 174             funz->SetParameters(param[0],param[1],param[2],param[3]);
 175             funz->SetLineColor(2);
 176             funz->DrawCopy("LSAME");
 177             t0->DrawLatex(0.15,0.92,text);
 178             c1->Update();
 179           }
 180         }
 181       }
 182     }
 183     iev++;
 184     readEv++;
 185     printf(" --- OK\n");
 186   }while(rd->NextEvent()&&iev<=lastEv);
 187   printf("Total number of events = %d\n",readEv);
 188   Float_t nfac=1./(Float_t)readEv/33./520.;
 189   hanst->Scale(nfac);
 190   nfac=1./(Float_t)readEv;
 191   hpad7l->Scale(nfac);
 192   hpad7r->Scale(nfac);
 193
 194   Int_t ipt0=0, ipt1=0;
 195   Float_t Edrift=(1800-45)/291/0.012;
 196   TFile *outfil=new TFile("DriftSpeedHistos.root","recreate");
 197   for(Int_t iMod=0; iMod<260; iMod++){
 198     outfil->cd();
 199     hvdriftl[iMod]->Write();
 200     hvdriftr[iMod]->Write();
 201     Float_t modid=iMod+240;
 202     if(hvdriftl[iMod]->GetEntries()>0){
 203       Float_t avevell=hvdriftl[iMod]->GetMean();
 204       Float_t rmsvell=hvdriftl[iMod]->GetRMS();
 205       if(avevell > 5.5 && avevell < 8.5){
 206         gvvsmod0->SetPoint(ipt0,modid,avevell);
 207         gvvsmod0->SetPointError(ipt0,0,rmsvell);
 208         Float_t mob=avevell*1.E5/Edrift;
 209         Float_t temper=293.15*TMath::Power((mob/1350.),-1/2.4);
 210         gtvsmod0->SetPoint(ipt0,modid,temper);
 211         ++ipt0;
 212       }
 213     }
 214     if(hvdriftr[iMod]->GetEntries()>0){
 215       Float_t avevelr=hvdriftr[iMod]->GetMean();
 216       Float_t rmsvelr=hvdriftr[iMod]->GetRMS();
 217       if(avevelr > 5.5 && avevelr < 8.5){
 218         gvvsmod1->SetPoint(ipt1,modid,avevelr);
 219         gvvsmod1->SetPointError(ipt1,0,rmsvelr);
 220         Float_t mob=avevelr*1.E5/Edrift;
 221         Float_t temper=293.15*TMath::Power((mob/1350.),-1./2.4);
 222         gtvsmod1->SetPoint(ipt1,modid,temper);
 223         ++ipt1;
 224       }
 225     }
 226   }
 227   gvvsmod0->SetName("gvvsmod0");
 228   gvvsmod1->SetName("gvvsmod1");
 229   gtvsmod0->SetName("gtvsmod0");
 230   gtvsmod1->SetName("gtvsmod1");
 231   outfil->cd();
 232   gvvsmod0->Write();
 233   gvvsmod1->Write();
 234   gtvsmod0->Write();
 235   gtvsmod1->Write();
 236   outfil->Close();
 237
 238   TCanvas* c8=new TCanvas("c8","Drift Speed vs. mod");
 239   gvvsmod0->SetTitle("");
 240   gvvsmod1->SetTitle("");
 241   gvvsmod0->SetMarkerStyle(20);
 242   gvvsmod1->SetMarkerStyle(21);
 243   gvvsmod1->SetMarkerColor(2);
 244   gvvsmod0->Draw("AP");
 245   gvvsmod0->SetMinimum(gvmin);
 246   gvvsmod0->SetMaximum(gvmax);
 247   gvvsmod0->GetXaxis()->SetTitle("Module Number");
 248   Char_t title[25];
 249   sprintf(title,"Vdrift at injector pad %d",jpad);
 250   gvvsmod0->GetYaxis()->SetTitle(title);
 251   gvvsmod1->Draw("PSAME");
 252   TLatex* tleft=new TLatex(0.7,0.82,"Side 0");
 253   tleft->SetNDC();
 254   tleft->SetTextColor(1);
 255   tleft->Draw();
 256   TLatex* tright=new TLatex(0.7,0.75,"Side 1");
 257   tright->SetNDC();
 258   tright->SetTextColor(2);
 259   tright->Draw();
 260
 261   TLine *lin=new TLine(323,gvmin,323,gvmax);
 262   lin->SetLineColor(4);
 263   lin->Draw();
 264   c8->Update();
 265   c8->SaveAs("VdriftVsMod.gif");
 266
 267   TCanvas* c8t=new TCanvas("c8t","Temeprature vs. mod");
 268   gtvsmod0->SetTitle("");
 269   gtvsmod1->SetTitle("");
 270   gtvsmod0->SetMarkerStyle(20);
 271   gtvsmod1->SetMarkerStyle(21);
 272   gtvsmod1->SetMarkerColor(2);
 273   gtvsmod0->Draw("AP");
 274   gtvsmod0->SetMinimum(gtmin);
 275   gtvsmod0->SetMaximum(gtmax);
 276   gtvsmod0->GetXaxis()->SetTitle("Module Number");
 277   sprintf(title,"Estimated Temperature (K)");
 278   gtvsmod0->GetYaxis()->SetTitle(title);
 279   gtvsmod1->Draw("PSAME");
 280   tleft->Draw();
 281   tright->Draw();
 282   TLine *lint=new TLine(323,gtmin,323,gtmax);
 283   lint->SetLineColor(4);
 284   lint->Draw();
 285   c8t->Update();
 286   c8t->SaveAs("TempVsMod.gif");
 287
 288   TCanvas* c9=new TCanvas("c9","Injector status");
 289   hanst->SetStats(0);
 290   hanst->Draw();
 291   hanst->GetXaxis()->SetTitle("Injector pad status");
 292   hanst->GetXaxis()->CenterTitle();
 293   c9->SaveAs("InjStatus.gif");
 294
 295 //   TCanvas* c10=new TCanvas("c10","Pad status 7",1200,600);
 296 //   hpad7l->SetStats(0);
 297 //   hpad7r->SetStats(0);
 298 //   c10->Divide(2,1);
 299 //   c10->cd(1);
 300 //   hpad7l->Draw();
 301 //   hpad7l->GetXaxis()->SetTitle("Side Left -- Pad number");
 302 //   hpad7l->GetXaxis()->CenterTitle();
 303 //   hpad7l->GetYaxis()->SetTitle("Number of status 7");
 304 //   c10->cd(2);
 305 //   hpad7r->Draw();
 306 //   hpad7r->GetXaxis()->SetTitle("Side Right -- Pad number");
 307 //   hpad7r->GetXaxis()->CenterTitle();
 308 //   hpad7r->GetYaxis()->SetTitle("Number of status 7");
 309 //   printf("Side 0, maximum pad=%d\n",hpad7l->GetMaximumBin());
 310 //   printf("Side 1, maximum pad=%d\n",hpad7r->GetMaximumBin());
 311
 312
 313 }
 314
 315 void AnalyzeSDDInjectorsAllMod(Int_t nrun, Int_t n2, Char_t* dir="LHC08d_SDD", Int_t nDDL=0, Int_t firstEv=15, Int_t lastEv=15){
 316   TGrid::Connect("alien:",0,0,"t");
 317   Char_t filnam[200];
 318   sprintf(filnam,"alien:///alice/data/2008/%s/%09d/raw/08%09d%03d.10.root",dir,nrun,nrun,n2);
 319   printf("Open file %s\n",filnam);
 320   AnalyzeSDDInjectorsAllMod(filnam,nDDL,firstEv,lastEv);
 321 }
 322
 323
 324