[u/mrichter/AliRoot.git] / ITS / ITSPlotTracks.C

#include <iostream.h>
#include <fstream.h>

void ITSPlotTracks(){  
  
  ifstream in ("AliITSTra.out");
  
  TVector DataOut(10);
  
  ///////////////////////////////// Histograms definition ///////////////////////////////////////////

  
  TH1F *hp=new TH1F("hp","PHI resolution",50,-15.,15.); hp->SetFillColor(4);     
  TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-15.,15.); hl->SetFillColor(4);  
  
 //TH1F *hp=new TH1F("hp","PHI resolution",100,-50.,50.); hp->SetFillColor(4);     
 // TH1F *hl=new TH1F("hl","LAMBDA resolution",100,-50.,50.); hl->SetFillColor(4); 
      
 // TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-20.,20.);
 TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-10.,10.); 
  hpt->SetFillColor(2); 
  TH1F *hd=new TH1F("hd","Impact parameter distribution ",100,0,2000); 
  hd->SetFillColor(6);
  
  TH1F *hdr=new TH1F("hdr","Dr ",50,-1000,1000);  
  hdr->SetFillColor(kGreen);
  TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);  
  hdz->SetFillColor(kBlue);
  
  /*
  TH1F *hdr=new TH1F("hdr","Dr ",50,-100,100);  
  hdr->SetFillColor(kGreen);
  TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);  
  hdz->SetFillColor(kBlue);  
  */ 

  TH1F *hgood=new TH1F("hgood","Good tracks",10,0,2);
  TH1F *hfound=new TH1F("hfound","Found tracks",10,0,2); 
  TH1F *hfake=new TH1F("hfake","Fake tracks",10,0,2);    
  TH1F *hg=new TH1F("hg","",10,0,2); //efficiency for good tracks  
  hg->SetLineColor(4); hg->SetLineWidth(2);
  TH1F *hf=new TH1F("hf","Efficiency for fake tracks",10,0,2);
  hf->SetFillColor(1); hf->SetFillStyle(3013); hf->SetLineWidth(2);

  /////////////////////////////////////////////////////////////////////////////////////////////////// 

  ifstream in1 ("AliITSTrag.out");
  Double_t ptg; 
  for(;;) {
    in1 >> ptg; 
    if( in1.eof() ) break;
    hgood->Fill(ptg);
  }
  in1.close();
    
  for (;;){    
    for (int r=0; r<9; r++) in>>DataOut(r);
    if( in.eof() ) break;

    Double_t ptg=DataOut(0); Double_t labITS=DataOut(1); Double_t labTPC=DataOut(2); Double_t ptperc=DataOut(3);	  
    Double_t deltalam=DataOut(4); Double_t deltaphi=DataOut(5);
    Double_t Dtot=DataOut(6); Double_t Dr=DataOut(7); Double_t Dz=DataOut(8);	

	  if(labITS>=0) hfound->Fill(ptg); else 	{   hfake->Fill(ptg);}

	  if(labITS>=0 ) {	  
      hpt->Fill(ptperc);
      hl->Fill(deltalam);
      hp->Fill(deltaphi);
      hd->Fill(Dtot);
      hdr->Fill(Dr);
      hdz->Fill(Dz);
   }    
  }
 

  in.close();
  Stat_t ngood=hgood->GetEntries(); cerr<<"Good tracks "<<ngood<<endl;  
  Stat_t nfound=hfound->GetEntries(); cerr<<"Found tracks "<<nfound<<endl;
  Stat_t nfake=hfake->GetEntries(); cerr<<"Fake tracks "<<nfake<<endl;       
  gStyle->SetOptStat(111110);
  gStyle->SetOptFit(1);	   
  TCanvas *c1=new TCanvas("c1","",0,0,700,700);
  TPad *p1=new TPad("p1","",0,0.5,0.5,1); p1->Draw(); hp->SetXTitle("(mrad)");
  p1->cd(); hp->Draw();  hp->Fit("gaus"); c1->cd();
  TPad *p2=new TPad("p2","",0.5,0.5,1,1); p2->Draw(); hl->SetXTitle("(mrad)");
  p2->cd(); hl->Draw(); hl->Fit("gaus"); c1->cd();
  TPad *p3=new TPad("p3","",0,0,0.5,0.5); p3->Draw(); hpt->SetXTitle("(%)");
  p3->cd(); hpt->Draw(); hpt->Fit("gaus"); c1->cd();
  TPad *p4=new TPad("p4","",0.5,0,1,0.5); p4->Draw(); hd->SetXTitle("(micron)");
  p4->cd(); hd->Draw(); c1->cd();
   
  TCanvas *c3=new TCanvas("c3","",200,200,800,500);
  TPad *p7=new TPad("p7","",0,0,0.333,1); p7->Draw(); p7->cd(); hfound->Draw(); c3->cd(); 
  TPad *p8=new TPad("p8","",0.333,0,0.666,1); p8->Draw(); p8->cd(); hfake->Draw(); c3->cd();
  TPad *p9=new TPad("p9","",0.666,0,1,1); p9->Draw(); p9->cd(); hgood->Draw(); c3->cd();
     
  TCanvas *c4=new TCanvas("c4","",300,300,800,500);
  hg->Divide(hfound,hgood,1,1.,"b");
  hf->Divide(hfake,hgood,1,1.,"b");
  hg->SetMaximum(1.4);
  hg->SetYTitle("Tracking efficiency");
  hg->SetXTitle("Pt (GeV/c)");
  hg->Draw();
  TLine *line1 = new TLine(0,1.0,2,1.0); line1->SetLineStyle(4);
  line1->Draw("same");
  TLine *line2 = new TLine(0,0.9,2,0.9); line2->SetLineStyle(4);
  line2->Draw("same");
   
  hf->SetFillColor(1);
  hf->SetFillStyle(3013);
  hf->SetLineColor(2);
  hf->SetLineWidth(2);
  hf->Draw("histsame");
  TText *text = new TText(0.461176,0.248448,"Fake tracks");
  text->SetTextSize(0.05);
  text->Draw();
  text = new TText(0.453919,1.11408,"Good tracks");
  text->SetTextSize(0.05);
  text->Draw();
         
  TCanvas *c2=new TCanvas("c2","",100,100,700,400);
  TPad *p5=new TPad("p5","",0,0,0.5,1); p5->Draw(); hdr->SetXTitle("(micron)");
  p5->cd(); hdr->Draw();  hdr->Fit("gaus"); c2->cd();
  TPad *p6=new TPad("p6","",0.5,0,1,1); p6->Draw(); hdz->SetXTitle("(micron)");
  p6->cd(); hdz->Draw(); hdz->Fit("gaus"); c2->cd();
   
}
Commit	Line	Data
d1c69928	1	#include <iostream.h>
	2	#include <fstream.h>
	3
	4	void ITSPlotTracks(){
	5
	6	ifstream in ("AliITSTra.out");
	7
	8	TVector DataOut(10);
	9
	10	///////////////////////////////// Histograms definition ///////////////////////////////////////////
	11
	12
	13	TH1F *hp=new TH1F("hp","PHI resolution",50,-15.,15.); hp->SetFillColor(4);
	14	TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-15.,15.); hl->SetFillColor(4);
	15
	16	//TH1F *hp=new TH1F("hp","PHI resolution",100,-50.,50.); hp->SetFillColor(4);
	17	// TH1F *hl=new TH1F("hl","LAMBDA resolution",100,-50.,50.); hl->SetFillColor(4);
	18
	19	// TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-20.,20.);
	20	TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-10.,10.);
	21	hpt->SetFillColor(2);
	22	TH1F *hd=new TH1F("hd","Impact parameter distribution ",100,0,2000);
	23	hd->SetFillColor(6);
	24
	25	TH1F *hdr=new TH1F("hdr","Dr ",50,-1000,1000);
	26	hdr->SetFillColor(kGreen);
	27	TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);
	28	hdz->SetFillColor(kBlue);
	29
	30	/*
	31	TH1F *hdr=new TH1F("hdr","Dr ",50,-100,100);
	32	hdr->SetFillColor(kGreen);
	33	TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);
	34	hdz->SetFillColor(kBlue);
	35	*/
	36
	37	TH1F *hgood=new TH1F("hgood","Good tracks",10,0,2);
	38	TH1F *hfound=new TH1F("hfound","Found tracks",10,0,2);
	39	TH1F *hfake=new TH1F("hfake","Fake tracks",10,0,2);
	40	TH1F *hg=new TH1F("hg","",10,0,2); //efficiency for good tracks
	41	hg->SetLineColor(4); hg->SetLineWidth(2);
	42	TH1F *hf=new TH1F("hf","Efficiency for fake tracks",10,0,2);
	43	hf->SetFillColor(1); hf->SetFillStyle(3013); hf->SetLineWidth(2);
	44
	45	///////////////////////////////////////////////////////////////////////////////////////////////////
	46
	47	ifstream in1 ("AliITSTrag.out");
	48	Double_t ptg;
	49	for(;;) {
	50	in1 >> ptg;
	51	if( in1.eof() ) break;
	52	hgood->Fill(ptg);
	53	}
	54	in1.close();
	55
	56	for (;;){
	57	for (int r=0; r<9; r++) in>>DataOut(r);
	58	if( in.eof() ) break;
	59
	60	Double_t ptg=DataOut(0); Double_t labITS=DataOut(1); Double_t labTPC=DataOut(2); Double_t ptperc=DataOut(3);
	61	Double_t deltalam=DataOut(4); Double_t deltaphi=DataOut(5);
	62	Double_t Dtot=DataOut(6); Double_t Dr=DataOut(7); Double_t Dz=DataOut(8);
	63
	64	if(labITS>=0) hfound->Fill(ptg); else { hfake->Fill(ptg);}
65
66	if(labITS>=0 ) {
67	hpt->Fill(ptperc);
68	hl->Fill(deltalam);
69	hp->Fill(deltaphi);
70	hd->Fill(Dtot);
71	hdr->Fill(Dr);
72	hdz->Fill(Dz);
73	}
74	}
75
76
77	in.close();
78	Stat_t ngood=hgood->GetEntries(); cerr<<"Good tracks "<<ngood<<endl;
79	Stat_t nfound=hfound->GetEntries(); cerr<<"Found tracks "<<nfound<<endl;
80	Stat_t nfake=hfake->GetEntries(); cerr<<"Fake tracks "<<nfake<<endl;
81	gStyle->SetOptStat(111110);
82	gStyle->SetOptFit(1);
83	TCanvas *c1=new TCanvas("c1","",0,0,700,700);
84	TPad *p1=new TPad("p1","",0,0.5,0.5,1); p1->Draw(); hp->SetXTitle("(mrad)");
85	p1->cd(); hp->Draw(); hp->Fit("gaus"); c1->cd();
86	TPad *p2=new TPad("p2","",0.5,0.5,1,1); p2->Draw(); hl->SetXTitle("(mrad)");
87	p2->cd(); hl->Draw(); hl->Fit("gaus"); c1->cd();
88	TPad *p3=new TPad("p3","",0,0,0.5,0.5); p3->Draw(); hpt->SetXTitle("(%)");
89	p3->cd(); hpt->Draw(); hpt->Fit("gaus"); c1->cd();
90	TPad *p4=new TPad("p4","",0.5,0,1,0.5); p4->Draw(); hd->SetXTitle("(micron)");
91	p4->cd(); hd->Draw(); c1->cd();
92
93	TCanvas *c3=new TCanvas("c3","",200,200,800,500);
94	TPad *p7=new TPad("p7","",0,0,0.333,1); p7->Draw(); p7->cd(); hfound->Draw(); c3->cd();
95	TPad *p8=new TPad("p8","",0.333,0,0.666,1); p8->Draw(); p8->cd(); hfake->Draw(); c3->cd();
96	TPad *p9=new TPad("p9","",0.666,0,1,1); p9->Draw(); p9->cd(); hgood->Draw(); c3->cd();
97
98	TCanvas *c4=new TCanvas("c4","",300,300,800,500);
99	hg->Divide(hfound,hgood,1,1.,"b");
100	hf->Divide(hfake,hgood,1,1.,"b");
101	hg->SetMaximum(1.4);
102	hg->SetYTitle("Tracking efficiency");
103	hg->SetXTitle("Pt (GeV/c)");
104	hg->Draw();
105	TLine *line1 = new TLine(0,1.0,2,1.0); line1->SetLineStyle(4);
106	line1->Draw("same");
107	TLine *line2 = new TLine(0,0.9,2,0.9); line2->SetLineStyle(4);
108	line2->Draw("same");
109
110	hf->SetFillColor(1);
111	hf->SetFillStyle(3013);
112	hf->SetLineColor(2);
113	hf->SetLineWidth(2);
114	hf->Draw("histsame");
115	TText *text = new TText(0.461176,0.248448,"Fake tracks");
116	text->SetTextSize(0.05);
117	text->Draw();
118	text = new TText(0.453919,1.11408,"Good tracks");
119	text->SetTextSize(0.05);
120	text->Draw();
121
122	TCanvas *c2=new TCanvas("c2","",100,100,700,400);
123	TPad *p5=new TPad("p5","",0,0,0.5,1); p5->Draw(); hdr->SetXTitle("(micron)");
124	p5->cd(); hdr->Draw(); hdr->Fit("gaus"); c2->cd();
125	TPad *p6=new TPad("p6","",0.5,0,1,1); p6->Draw(); hdz->SetXTitle("(micron)");
126	p6->cd(); hdz->Draw(); hdz->Fit("gaus"); c2->cd();
127
128	}