ITS/UPGRADE/testITSU/testDetectorUp.C

   1
   2 // You have to load the class before ... ;-)
   3 // .L DetectorK.cxx++
   4
   5 //void standardPlots() {
   6 void testDetectorUp() {
   7
   8
   9   DetectorK its("ALICE","ITS");
  10
  11   its.AddLayer((char*)"bpipe",2.0.0,0.0022);
  12   its.AddLayer((char*)"vertex",     0,     0); // dummy vertex for matrix calculation
  13   // new ideal Pixel properties?
  14   Double_t x0IB     = 0.003;
  15   Double_t x0OB     = 0.008;
  16   Double_t resRPhiIB     = 0.0004;
  17   Double_t resZIB        = 0.0004;
  18   Double_t resRPhiOB     = 0.0004;
  19   Double_t resZOB        = 0.0004;
  20   Double_t eff           = 0.95;
  21   //
  22   /*
  23   its.AddLayer((char*)"ddd1",  2.32 ,  x0IB, resRPhiIB, resZIB,eff);
  24   its.AddLayer((char*)"ddd2",  3.13 ,  x0IB, resRPhiIB, resZIB,eff);
  25   its.AddLayer((char*)"ddd3",  3.91 ,  x0IB, resRPhiIB, resZIB,eff);
  26   its.AddLayer((char*)"ddd4",  19.41,  x0OB, resRPhiOB, resZOB,eff);
  27   its.AddLayer((char*)"ddd5",  24.71 ,  x0OB, resRPhiOB, resZOB,eff);
  28   its.AddLayer((char*)"ddd6",  35.33 ,  x0OB, resRPhiOB, resZOB,eff);
  29   its.AddLayer((char*)"ddd7",  40.53 ,  x0OB, resRPhiOB, resZOB,eff);
  30   */
  31   //
  32   its.AddLayer((char*)"ddd1",  2.32 ,  x0IB, resRPhiIB, resZIB,eff);
  33   its.AddLayer((char*)"ddd2",  3.13 ,  x0IB, resRPhiIB, resZIB,eff);
  34   its.AddLayer((char*)"ddd3",  3.91 ,  x0IB, resRPhiIB, resZIB,eff);
  35   //  its.AddLayer((char*)"ddd4",  19.41,  x0OB, resRPhiOB, resZOB,eff);
  36   //  its.AddLayer((char*)"ddd5",  24.71 ,  x0OB, resRPhiOB, resZOB,eff);
  37   its.AddLayer((char*)"ddd4",  5.,  x0OB, resRPhiOB, resZOB,eff);
  38   its.AddLayer((char*)"ddd5",  33. ,  x0OB, resRPhiOB, resZOB,eff);
  39   its.AddLayer((char*)"ddd6",  35.33 ,  x0OB, resRPhiOB, resZOB,eff);
  40   its.AddLayer((char*)"ddd7",  40.53 ,  x0OB, resRPhiOB, resZOB,eff);
  41
  42
  43
  44   its.SetAtLeastCorr(5);
  45   its.SetAtLeastFake(1);
  46   //
  47   its.PrintLayout();
  48   its.SolveViaBilloir(0);
  49
  50   its.MakeStandardPlots(0,2,1,kTRUE);
  51   //  return;
  52   its.AddTPC(0.1,0.1);
  53   its.SolveViaBilloir(0);
  54
  55   its.MakeStandardPlots(1,1,1,kTRUE);
  56
  57 }
  58
  59
  60 void particleDependendResolution() {
  61 // particle dependency on resolution
  62
  63   // .L Detector.cxx++
  64
  65   Detector its("ALICE","ITS");
  66
  67   its.MakeAliceCurrent();
  68   its.PrintLayout();
  69   its.PlotLayout();
  70
  71   its.SolveViaBilloir(0);
  72
  73   its.SetRadius("bpipe",2.1);
  74   its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
  75
  76   TCanvas *c1 = new TCanvas("c1","c1");
  77
  78   c1->Divide(2,1);
  79   c1->cd(1); gPad->SetGridx();   gPad->SetGridy();
  80   gPad->SetLogx(); //gPad->SetLogy();
  81   c1->cd(2); gPad->SetGridx();   gPad->SetGridy();
  82   gPad->SetLogx(); //gPad->SetLogy();
  83
  84
  85   // compare to telescope equation ?
  86   //  c1->cd(1); its.GetGraphPointingResolutionTeleEqu(0,1)->Draw("AC");
  87   //  c1->cd(2); its.GetGraphPointingResolutionTeleEqu(1,1)->Draw("AC");
  88
  89   its.SetParticleMass(0.140); // pion
  90   its.SolveViaBilloir(0,0);
  91   c1->cd(1); its.GetGraphPointingResolution(0,1)->Draw("AC");
  92   c1->cd(2); its.GetGraphPointingResolution(1,1)->Draw("AC");
  93
  94   its.SetParticleMass(0.498); // kaon
  95   its.SolveViaBilloir(0,0);
  96   c1->cd(1); its.GetGraphPointingResolution(0,2)->Draw("C");
  97   c1->cd(2); its.GetGraphPointingResolution(1,2)->Draw("C");
  98
  99   its.SetParticleMass(0.00051); // electron
 100   its.SolveViaBilloir(0,0);
 101   c1->cd(1); its.GetGraphPointingResolution(0,3)->Draw("C");
 102   c1->cd(2); its.GetGraphPointingResolution(1,3)->Draw("C");
 103
 104   its.SetParticleMass(0.938); // proton
 105   its.SolveViaBilloir(0,0);
 106   c1->cd(1); its.GetGraphPointingResolution(0,4)->Draw("C");
 107   c1->cd(2); its.GetGraphPointingResolution(1,4)->Draw("C");
 108
 109
 110
 111 }
 112