{	
		TString inFile="out_test_2pi_tpc_a.root";
	

	TFile *f = TFile::Open(inFile);
	
	
	gStyle->SetOptStat(kFALSE);
	gStyle->SetPalette(1);

	TCanvas *c4=new TCanvas("c4","",800,800);
	c4->Divide(1,3);
	c4->cd(1);
	nt->Draw("mc_p>>hMcP(100,0,5)");
	nt->Draw("mc_p>>hMcAccTpcP(100,0,5)","mc_tpc>0","E1");
	hMcAccTpcP->Divide(hMcP);
	hMcAccTpcP->SetTitle("Geometrical Acceptance");
	hMcAccTpcP->GetXaxis()->SetTitle("Momentum (GeV)");
	hMcAccTpcP->GetYaxis()->SetTitle("Acceptance");
	c4->cd(2);
	nt->Draw("mc_theta>>hMcTheta(50,0,180)");
	nt->Draw("mc_theta>>hMcAccTpcTheta(50,0,180)","mc_tpc>0","E1");
	hMcAccTpcTheta->Divide(hMcTheta);
	hMcAccTpcTheta->SetTitle("Geometrical Acceptance");
	hMcAccTpcTheta->GetXaxis()->SetTitle("#theta (degree)");
	hMcAccTpcTheta->GetYaxis()->SetTitle("Acceptance");
	c4->cd(3);
	nt->Draw("mc_phi>>hMcPhi(50,-200,200)");
	nt->Draw("mc_phi>>hMcAccTpcPhi(50,-200,200)","mc_tpc>0","E1");
	hMcAccTpcPhi->Divide(hMcPhi);
	hMcAccTpcPhi->SetTitle("Geometrical Acceptance");
	hMcAccTpcPhi->GetXaxis()->SetTitle("#phi (degree)");
	hMcAccTpcPhi->GetYaxis()->SetTitle("Acceptance");
	
	
	TCanvas *c5=new TCanvas("c5","",800,800);
	float numero=ntEvt->GetEntries(); 
	c5->cd();
	TH1F *f1=new TH1F("f1","",100,0,5);
	ntEvt.Draw("acc_tpc>>f1");
	f1->SetTitle("Geometrical Global Acceptance");
	f1->GetXaxis()->SetTitle("Primary Tracks");
	f1->GetYaxis()->SetTitle("Tpc Acceptance");
	f1->GetYaxis()->SetTitleOffset(1.5);
	f1->Draw();
	f1->Scale(1/numero);


	TCanvas *c6=new TCanvas("c6","",800,800);
	c6->Divide(1,3);
	c6->cd(1);
	nt->Draw("mc_p>>hMcP(100,0,5)");
	nt->Draw("mc_p>>hRecoEffP(100,0,5)","mult>0","E1");
	hRecoEffP->Divide(hMcP);
	hRecoEffP->SetTitle("Pattern Recognition Efficiency");
	hRecoEffP->GetXaxis()->SetTitle("Momentum (GeV)");
	hRecoEffP->GetYaxis()->SetTitle("Efficiency");
	hRecoEffP->Draw();
	c6->cd(2);
	nt->Draw("mc_theta>>hMcTheta(50,0,180)");
	nt->Draw("mc_theta>>hRecoEffTheta(50,0,180)","mult>0","E1");
	hRecoEffTheta->Divide(hMcTheta);
	hRecoEffTheta->SetTitle("Pattern Recognition Efficiency");
	hRecoEffTheta->GetXaxis()->SetTitle("#theta (degree)");
	hRecoEffTheta->GetYaxis()->SetTitle("Efficiency");
	hRecoEffTheta->Draw();
	c6->cd(3);
	nt->Draw("mc_phi>>hMcPhi(50,-200,200)");
	nt->Draw("mc_phi>>hRecoEffPhi(50,-200,200)","mult>0","E1");
	hRecoEffPhi->Divide(hMcPhi);
	hRecoEffPhi->SetTitle("Pattern Recognition Efficiency");
	hRecoEffPhi->GetXaxis()->SetTitle("#phi (degree)");
	hRecoEffPhi->GetYaxis()->SetTitle("Efficiency");
	hRecoEffPhi->Draw();
	
	TCanvas *c7=new TCanvas("c7","",800,800);

	c7->cd();
	TH1F *f1=new TH1F("f1","",100,0,5);
	ntEvt.Draw("cand>>f1");
	f1->SetTitle("Number of candidates");
	f1->GetXaxis()->SetTitle("Candidates");
	f1->GetYaxis()->SetTitle("1/evt");
	f1->GetYaxis()->SetTitleOffset(1.5);
	f1->Draw();
	f1->Scale(1/numero);

	TCanvas *c8=new TCanvas("c8","",800,800);

	c8->cd();
	TH1F *f1=new TH1F("f1","",100,0,5);
	ntEvt.Draw("eff>>f1");
	f1->SetTitle("Number of reconstructed primary tracks");
	f1->GetXaxis()->SetTitle("Primary Tracks");
	f1->GetYaxis()->SetTitle("Efficiency");
	f1->GetYaxis()->SetTitleOffset(1.5);
	f1->Draw();
	f1->Scale(1/numero);

	TCanvas *c9=new TCanvas("c9","",800,800);
        c9->Divide(3,2);
        c9->cd(1);
        nt->Draw("(mc_p-p):mc_p>>hresp_p(100,0,5,100,-0.5,0.5)","mult>0 && (stt>0||tpc>0)","colz");
        hresp_p->SetTitle("Resolution");
        hresp_p->GetXaxis()->SetTitle("MC Momentum (GeV)");
        hresp_p->GetYaxis()->SetTitle("(MC Mom - Reco Mom )");
        hresp_p->GetYaxis()->SetTitleOffset(1.5);

        c9->cd(2);
        nt->Draw("(mc_theta-theta):mc_theta>>hrestheta_theta(100,0,100,100,-2,2)","mult>0 && (stt>0||tpc>0)","colz");
        hrestheta_theta->SetTitle("Resolution");
        hrestheta_theta->GetXaxis()->SetTitle("MC #theta (degree)");
        hrestheta_theta->GetYaxis()->SetTitle("Mc #theta - Reco #theta ");
        hrestheta_theta->GetYaxis()->SetTitleOffset(1.5);

        c9->cd(3);
        nt->Draw("(mc_p-p):mc_theta>>hresp_theta(100,0,100,100,-0.5,0.5)","mult>0 &&  (stt>0||tpc>0)","colz");
        hresp_theta->SetTitle("Resolution");
        hresp_theta->GetXaxis()->SetTitle("MC #theta (degree)");
        hresp_theta->GetYaxis()->SetTitle("MC Mom - Reco Mom");
        hresp_theta->GetYaxis()->SetTitleOffset(1.5);

        c9->cd(4);
        nt->Draw("(mc_theta-theta):mc_p>>hrestheta_p(100,0,5,100,-2,2)","mult>0 &&  (stt>0||tpc>0)","colz");
        hrestheta_p->SetTitle("Resolution");
        hrestheta_p->GetXaxis()->SetTitle("MC Momentum (GeV)");
        hrestheta_p->GetYaxis()->SetTitle("Mc #theta - Reco #theta ");
        hrestheta_p->GetYaxis()->SetTitleOffset(1.5);

        c9->cd(5);
        nt->Draw("(mc_phi-phi):mc_phi>>hresphi_phi(100,-200,200,100,-2,2)","mult>0 &&  (stt>0||tpc>0)","colz");
        hresphi_phi->SetTitle("Resolution");
        hresphi_phi->GetXaxis()->SetTitle("MC Phi (degree)");
        hresphi_phi->GetYaxis()->SetTitle("Mc #phi - Reco #phi ");
        hresphi_phi->GetYaxis()->SetTitleOffset(1.5);

        c9->Update();

	TCanvas *c20=new TCanvas("c20","",800,800);
        c20->cd();
        TH1F *hist=new TH1F("hist","Total Momentum distribution",100,0,10);
        hist->GetXaxis()->SetTitle("MC Momentum (GeV)");
        nt->Draw("p>>hist","mult>0 && (stt>0 || tpc>0)");

        TCanvas *c21=new TCanvas("c21","",800,800);
        c21->cd();
        TH1F *hist2=new TH1F("hist2","#theta Angle Distribution",100,0,180);
        hist2->GetXaxis()->SetTitle("Theta (degree)");
        nt->Draw("theta>>hist2","mult>0 && (stt>0 || tpc>0)");

        TCanvas *c22=new TCanvas("c22","",800,800);
        c22->cd();
        TH1F *hist3=new TH1F("hist3","(Reconstructed Momentum - MC Momentum )/MC Momentum",100,-1,1);
        hist3->SetStats(kTRUE);
        hist3->GetXaxis()->SetTitle("(Reco-MC)/MC");
        nt->Draw("(mc_p-p)/mc_p>>hist3","mult>0 && (stt>0 || tpc>0)");
        TF1 *f10 = new TF1("f10","gaus",-0.09,0.09);
        hist3->Fit(f10,"R");

		
}