TPC/CalibMacros/distortionMapDraw.C

   1 //
   2 // This is an example macro for drawing of the distortion maps which were created during the CPass0
   3 // respectivally,  CPass1
   4 //
   5 // Macro to be extended - resulting QA plots to be published on the web page
   6 //
   7 void distortionMapDraw(){
   8   //
   9   // 1.) connect the grid
  10   //
  11   TGrid::Connect("alien");
  12   //
  13   // 2.) Open calibration of interest. e.g
  14   //
  15   TFile *f  = TFile::Open("alien:///alice/data/2011/LHC11h/000170593/cpass0_HLT/OCDB/meanITSVertex.root");
  16   //
  17   // 3.) Get the content - distortion map
  18   //   TPC=Vertex, TPC-ITS, TPC-TRD, TPC-TOF
  19   //   f.ls()
  20   //   TAlienFile**            alien:///alice/data/2011/LHC11h/000170593/cpass0_HLT/OCDB/meanITSVertex.root             //   TAlienFile*            alien:///alice/data/2011/LHC11h/000170593/cpass0_HLT/OCDB/meanITSVertex.root
  21   //   KEY: TTree    ITSdy;1 ITSdy
  22   //   KEY: TTree    ITSdz;1 ITSdz
  23   //   KEY: TTree    ITSdsnp;1       ITSdsnp
  24   //   KEY: TTree    Vertexdy;1      Vertexdy
  25   //   KEY: TTree    Vertexdsnp;1    Vertexdsnp
  26   //   KEY: TTree    Vertexdz;1      Vertexdz
  27   //   KEY: TTree    TOFdy;1 TOFdy
  28   //   KEY: TTree    TRDdy;1 TRDdy
  29   //
  30   // 4.) To see the variables of the distrotion tree
  31   //     ITSdy.Print()
  32   //
  33   // 5. Example draw - Note that  for the example run we show residuals obtained in  CPass0
  34   //    before run dependent alignment and ExB twist
  35   //
  36   ITSdy->SetMarkerStyle(25);
  37   Vertexdy->SetMarkerStyle(25);
  38   TOFdy->SetMarkerStyle(25);
  39   // Draw mean residual
  40   ITSdy->Draw("mean","entries>100&&theta<0&&abs(snp)<0.1","");
  41   //
  42   // Draw mean residual/distortion as function of sector position - theta in color code
  43   //
  44   ITSdy->Draw("mean:sector:abs(theta)","entries>100&&theta<0","colz");
  45   Vertexdy->Draw("mean:sector:abs(theta)","entries>100&&theta<0","colz")
  46   TOFdy->Draw("mean:sector:abs(theta)","entries>100&&theta<0","colz");
  47 }