[u/mrichter/AliRoot.git] / ITS / UPGRADE / testDetector.C


// You have to load the class before ... ;-)
// .L DetectorK.cxx++

void standardPlots() {

  DetectorK its("ALICE","ITS");

  its.MakeAliceCurrent(0,0); 

  // its.SetRadius("bpipe",2.1);
  // its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);

  its.PrintLayout();
  its.SolveViaBilloir();
 
  its.MakeStandardPlots(0,2);

  its.AddTPC(0.1,0.1);
  its.SolveViaBilloir();
 
  its.MakeStandardPlots(1,1);

}


void particleDependendResolution() { 
// particle dependency on resolution

  // .L Detector.cxx++

  Detector its("ALICE","ITS");

  its.MakeAliceCurrent(); 
  its.PrintLayout();
  its.PlotLayout();

  its.SolveViaBilloir(0);
 
  its.SetRadius("bpipe",2.1);
  its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);

  TCanvas *c1 = new TCanvas("c1","c1");
 
  c1->Divide(2,1);
  c1->cd(1); gPad->SetGridx();   gPad->SetGridy(); 
  gPad->SetLogx(); //gPad->SetLogy();
  c1->cd(2); gPad->SetGridx();   gPad->SetGridy(); 
  gPad->SetLogx(); //gPad->SetLogy();

  
  // compare to telescope equation ?
  //  c1->cd(1); its.GetGraphPointingResolutionTeleEqu(0,1)->Draw("AC");
  //  c1->cd(2); its.GetGraphPointingResolutionTeleEqu(1,1)->Draw("AC");

  its.SetParticleMass(0.140); // pion  
  its.SolveViaBilloir(0,0);
  c1->cd(1); its.GetGraphPointingResolution(0,1)->Draw("AC");
  c1->cd(2); its.GetGraphPointingResolution(1,1)->Draw("AC");
 
  its.SetParticleMass(0.498); // kaon  
  its.SolveViaBilloir(0,0);
  c1->cd(1); its.GetGraphPointingResolution(0,2)->Draw("C");
  c1->cd(2); its.GetGraphPointingResolution(1,2)->Draw("C");

  its.SetParticleMass(0.00051); // electron  
  its.SolveViaBilloir(0,0);
  c1->cd(1); its.GetGraphPointingResolution(0,3)->Draw("C");
  c1->cd(2); its.GetGraphPointingResolution(1,3)->Draw("C");

  its.SetParticleMass(0.938); // proton  
  its.SolveViaBilloir(0,0);
  c1->cd(1); its.GetGraphPointingResolution(0,4)->Draw("C");
  c1->cd(2); its.GetGraphPointingResolution(1,4)->Draw("C");


}
Commit	Line	Data
54fc064a	1
54fc064a	2	// You have to load the class before ... ;-)
caadc8fb	3	// .L DetectorK.cxx++
54fc064a	4
	5	void standardPlots() {
	6
caadc8fb	7	DetectorK its("ALICE","ITS");
54fc064a	8
c40d3fcb	9	its.MakeAliceCurrent(0,0);
54fc064a	10
c40d3fcb	11	// its.SetRadius("bpipe",2.1);
c40d3fcb	12	// its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
54fc064a	13
	14	its.PrintLayout();
	15	its.SolveViaBilloir();
	16
c40d3fcb	17	its.MakeStandardPlots(0,2);
54fc064a	18
c40d3fcb	19	its.AddTPC(0.1,0.1);
c40d3fcb	20	its.SolveViaBilloir();
54fc064a	21
c40d3fcb	22	its.MakeStandardPlots(1,1);
54fc064a	23
	24	}
	25
	26
	27	void particleDependendResolution() {
	28	// particle dependency on resolution
	29
c40d3fcb	30	// .L Detector.cxx++
c40d3fcb	31
54fc064a	32	Detector its("ALICE","ITS");
	33
	34	its.MakeAliceCurrent();
54fc064a	35	its.PrintLayout();
c40d3fcb	36	its.PlotLayout();
c40d3fcb	37
54fc064a	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	}