ITS/UPGRADE/testDetector.C

   1
   2 // You have to load the class before ... ;-)
   3 // .L DetectorK.cxx++
   4
   5 void standardPlots() {
   6
   7   DetectorK its("ALICE","ITS");
   8
   9   its.MakeAliceCurrent(0,0);
  10
  11   // its.SetRadius("bpipe",2.1);
  12   // its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
  13
  14   its.PrintLayout();
  15   its.SolveViaBilloir();
  16
  17   its.MakeStandardPlots(0,2);
  18
  19   its.AddTPC(0.1,0.1);
  20   its.SolveViaBilloir();
  21
  22   its.MakeStandardPlots(1,1);
  23
  24 }
  25
  26
  27 void particleDependendResolution() {
  28 // particle dependency on resolution
  29
  30   // .L Detector.cxx++
  31
  32   Detector its("ALICE","ITS");
  33
  34   its.MakeAliceCurrent();
  35   its.PrintLayout();
  36   its.PlotLayout();
  37
  38   its.SolveViaBilloir(0);
  39
  40   its.SetRadius("bpipe",2.1);
  41   its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
  42
  43   TCanvas *c1 = new TCanvas("c1","c1");
  44
  45   c1->Divide(2,1);
  46   c1->cd(1); gPad->SetGridx();   gPad->SetGridy();
  47   gPad->SetLogx(); //gPad->SetLogy();
  48   c1->cd(2); gPad->SetGridx();   gPad->SetGridy();
  49   gPad->SetLogx(); //gPad->SetLogy();
  50
  51
  52   // compare to telescope equation ?
  53   //  c1->cd(1); its.GetGraphPointingResolutionTeleEqu(0,1)->Draw("AC");
  54   //  c1->cd(2); its.GetGraphPointingResolutionTeleEqu(1,1)->Draw("AC");
  55
  56   its.SetParticleMass(0.140); // pion
  57   its.SolveViaBilloir(0,0);
  58   c1->cd(1); its.GetGraphPointingResolution(0,1)->Draw("AC");
  59   c1->cd(2); its.GetGraphPointingResolution(1,1)->Draw("AC");
  60
  61   its.SetParticleMass(0.498); // kaon
  62   its.SolveViaBilloir(0,0);
  63   c1->cd(1); its.GetGraphPointingResolution(0,2)->Draw("C");
  64   c1->cd(2); its.GetGraphPointingResolution(1,2)->Draw("C");
  65
  66   its.SetParticleMass(0.00051); // electron
  67   its.SolveViaBilloir(0,0);
  68   c1->cd(1); its.GetGraphPointingResolution(0,3)->Draw("C");
  69   c1->cd(2); its.GetGraphPointingResolution(1,3)->Draw("C");
  70
  71   its.SetParticleMass(0.938); // proton
  72   its.SolveViaBilloir(0,0);
  73   c1->cd(1); its.GetGraphPointingResolution(0,4)->Draw("C");
  74   c1->cd(2); its.GetGraphPointingResolution(1,4)->Draw("C");
  75
  76
  77
  78 }