[u/mrichter/AliRoot.git] / PWG2 / EVCHAR / combscalingfactor.C

#include "Riostream.h"
#include "TROOT.h"
#include "TStyle.h"
#include "TMath.h"
#include "TH1.h"
#include "TH2.h"
#include "TF1.h" 
#include "TFile.h"
#include "TList.h"
#include "TTree.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TCanvas.h"
void combscalingfactor(Char_t *fileall, Char_t *filecomb, Bool_t scalebkg = kFALSE, Float_t scfac=1.) {

 //  gStyle->SetOptStat(kFALSE);
  gStyle->SetOptTitle(kFALSE);
//  gStyle->SetOptLogy(kFALSE);
  gStyle->SetFrameLineWidth(2.5);
  gStyle->SetCanvasColor(0);
  gStyle->SetPadColor(0);
  gStyle->SetHistLineWidth(2.5);
  gStyle->SetLabelSize(0.05, "x");
  gStyle->SetLabelSize(0.05, "y");
  gStyle->SetTitleSize(0.05, "x");
  gStyle->SetTitleSize(0.05, "y");
  gStyle->SetTitleOffset(1.1, "x");
  gStyle->SetPadBottomMargin(0.14);
  gStyle->SetFrameFillColor(0);
  gStyle->SetFillColor(0);
//  gROOT->ForceStyle();


 TFile *fall = new TFile(fileall);
 TFile *fcomb = new TFile(filecomb);

 TList *lall = (TList*) fall->Get("cOutput");
 TList *lcomb = (TList*) fcomb->Get("cOutput");

 TH2F *hall = (TH2F*) lall->FindObject("fHistSPDdePhideTheta");
 TH2F *hcomb = (TH2F*) lcomb->FindObject("fHistSPDdePhideTheta");

 TH1D* hproall = new TH1D("_xall","",2000,-1.,1.);
 hproall->Sumw2(); 
 hproall = hall->ProjectionX("_xall",0,-1,"e");

 TH1D* hprocomb = new TH1D("_xcomb","",2000,-1.,1.);
 hprocomb->Sumw2(); 
 hprocomb = hcomb->ProjectionX("_xcomb",0,-1,"e");


 TH1D* hratio = new TH1D("ratio","",2000,-1.,1.); 
 hratio->Divide(hproall,hprocomb,1.,1.); 

 TF1 *fleft = new TF1("fleft","pol0",-.8,-0.6); //check values
 TF1 *fright = new TF1("fright","pol0",0.6,.8);

 hratio->Fit("fleft","R");
 hratio->Fit("fright","R");

/*
 // Projection on DeltaTheta 
 TH1D* hproall = new TH1D("_xall","",1000,-.25,.25);
 hproall->Sumw2(); 
 hproall = hall->ProjectionY("_xall",0,-1,"e");

 TH1D* hprocomb = new TH1D("_xcomb","",1000,-.25,.25);
 hprocomb->Sumw2(); 
 hprocomb = hcomb->ProjectionY("_xcomb",0,-1,"e");

 TH1D* hratio = new TH1D("ratio","",1000,-.25,.25); 
 hratio->Divide(hproall,hprocomb,1.,1.);
 
 TF1 *fleft = new TF1("fleft","pol0",-.8,-0.6); //check values
 TF1 *fright = new TF1("fright","pol0",0.6,.8);
 hratio->Fit("fleft","R");
 hratio->Fit("fright","R");
*/

  TFile *fout= new TFile("scaling.root","RECREATE");

  new TCanvas();
  hproall->Draw("histo");
  
  if (scalebkg) hprocomb->Scale(scfac);
  hprocomb->Draw("histo,same");

  cout<<"Percentage of background"<< hprocomb->Integral(1,1000)/hproall->Integral(1,1000)<<endl;
  cout<<"Percentage of background in |#D#phi|<0.08"<< hprocomb->Integral(921,1080)/hproall->Integral(921,1080)<<endl; 
  hproall->Write();
  hprocomb->SetLineColor(kBlue);
  hprocomb->Write(); 
  fout->Write();
  fout->Close();

}
Commit	Line	Data
c6788f58	1	#include "Riostream.h"
	2	#include "TROOT.h"
	3	#include "TStyle.h"
	4	#include "TMath.h"
	5	#include "TH1.h"
	6	#include "TH2.h"
	7	#include "TF1.h"
	8	#include "TFile.h"
	9	#include "TList.h"
	10	#include "TTree.h"
	11	#include "TLegend.h"
	12	#include "TLegendEntry.h"
	13	#include "TCanvas.h"
	14	void combscalingfactor(Char_t fileall, Char_t filecomb, Bool_t scalebkg = kFALSE, Float_t scfac=1.) {
	15
	16	// gStyle->SetOptStat(kFALSE);
	17	gStyle->SetOptTitle(kFALSE);
	18	// gStyle->SetOptLogy(kFALSE);
	19	gStyle->SetFrameLineWidth(2.5);
	20	gStyle->SetCanvasColor(0);
	21	gStyle->SetPadColor(0);
	22	gStyle->SetHistLineWidth(2.5);
	23	gStyle->SetLabelSize(0.05, "x");
	24	gStyle->SetLabelSize(0.05, "y");
	25	gStyle->SetTitleSize(0.05, "x");
	26	gStyle->SetTitleSize(0.05, "y");
	27	gStyle->SetTitleOffset(1.1, "x");
	28	gStyle->SetPadBottomMargin(0.14);
	29	gStyle->SetFrameFillColor(0);
	30	gStyle->SetFillColor(0);
	31	// gROOT->ForceStyle();
	32
	33
	34
	35	TFile *fall = new TFile(fileall);
	36	TFile *fcomb = new TFile(filecomb);
	37
	38	TList lall = (TList) fall->Get("cOutput");
	39	TList lcomb = (TList) fcomb->Get("cOutput");
	40
	41	TH2F hall = (TH2F) lall->FindObject("fHistSPDdePhideTheta");
	42	TH2F hcomb = (TH2F) lcomb->FindObject("fHistSPDdePhideTheta");
	43
	44	TH1D* hproall = new TH1D("_xall","",2000,-1.,1.);
	45	hproall->Sumw2();
	46	hproall = hall->ProjectionX("_xall",0,-1,"e");
	47
	48	TH1D* hprocomb = new TH1D("_xcomb","",2000,-1.,1.);
	49	hprocomb->Sumw2();
	50	hprocomb = hcomb->ProjectionX("_xcomb",0,-1,"e");
	51
	52
	53	TH1D* hratio = new TH1D("ratio","",2000,-1.,1.);
	54	hratio->Divide(hproall,hprocomb,1.,1.);
	55
	56	TF1 *fleft = new TF1("fleft","pol0",-.8,-0.6); //check values
	57	TF1 *fright = new TF1("fright","pol0",0.6,.8);
	58
	59	hratio->Fit("fleft","R");
	60	hratio->Fit("fright","R");
	61
	62	/*
	63	// Projection on DeltaTheta
	64	TH1D* hproall = new TH1D("_xall","",1000,-.25,.25);
65	hproall->Sumw2();
66	hproall = hall->ProjectionY("_xall",0,-1,"e");
67
68	TH1D* hprocomb = new TH1D("_xcomb","",1000,-.25,.25);
69	hprocomb->Sumw2();
70	hprocomb = hcomb->ProjectionY("_xcomb",0,-1,"e");
71
72	TH1D* hratio = new TH1D("ratio","",1000,-.25,.25);
73	hratio->Divide(hproall,hprocomb,1.,1.);
74
75	TF1 *fleft = new TF1("fleft","pol0",-.8,-0.6); //check values
76	TF1 *fright = new TF1("fright","pol0",0.6,.8);
77	hratio->Fit("fleft","R");
78	hratio->Fit("fright","R");
79	*/
80
81	TFile *fout= new TFile("scaling.root","RECREATE");
82
83	new TCanvas();
84	hproall->Draw("histo");
85
86	if (scalebkg) hprocomb->Scale(scfac);
87	hprocomb->Draw("histo,same");
88
89	cout<<"Percentage of background"<< hprocomb->Integral(1,1000)/hproall->Integral(1,1000)<<endl;
90	cout<<"Percentage of background in \|#D#phi\|<0.08"<< hprocomb->Integral(921,1080)/hproall->Integral(921,1080)<<endl;
91	hproall->Write();
92	hprocomb->SetLineColor(kBlue);
93	hprocomb->Write();
94	fout->Write();
95	fout->Close();
96
97	}