PWG2/EVCHAR/combscalingfactor.C

   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 }