higgspolarization.C

Go to the documentation of this file.

double beta(double x)
{
  double gamma = ((91.1876**4)/(x**4))*(8/((1-((2*(91.1876**2))/(x**2)))**2));
  return (-gamma+2)/(2+gamma);
 
}

double R_L(double x)
{
  double z = (91.1876*91.1876)/(x*x);
  
  double a = ((3-16*z+20*z*z)/(sqrt(4*z-1)))*acos((3*z-1)/(2*z*sqrt(z)));

  double b = ((1-z)/(2*z))*(2-13*z+15*z*z);

  double c = 0.5*(3-10*z+4*z*z)*log(z);

  return (a-b-c);

}



double R_T(double x)
{

  double y = (91.1876*91.1876)/(x*x);

  double s = ((3*(1-8*y+20*y*y))/(sqrt(4*y-1)))*acos((3*y-1)/(2*y*sqrt(y)));

  double t = ((1-y)/(2*y))*(2-13*y+47*y*y);

  double u = 1.5*(1-6*y+4*y*y)*log(y);

  return (s-t-u);
}





void higgspolarization(char *filelistname)

{
  TString filename;
  TFile * f;


  ifstream filelist(filelistname);
  string line;


  const int n = 37; //how large the arrray is that it creates

  double x[n],y[n],y1[n],y2[n],ex[n],ey[n],ey1[n],ey2[n],y3[n],ey3[n];


  int ipoint = 0;

  filelist >> line;
  double higgsmass = 140;
  
  //loop over the files from file list
 
  while (!filelist.eof())
    {
      filename = line;
      //open the file, now pointed by f


      f = TFile::Open(filename);
      
      TH1F * thish = (TH1F*) f->Get("Cosine1graph");

      TF1 * f2 = (TF1*)thish->GetListOfFunctions()->At(0);
      TF1 * f3 = (TF1*)thish->GetListOfFunctions()->At(1);
      TF1 * f4 = (TF1*)thish->GetListOfFunctions()->At(2);
      TF1 *f5 = new TF1("f5","[0]*( ((-[1]-1)/([1]-3)) * (3/4)*(1-x**2) + (3/8)*((2*[1]-2)/([1]-3))*(1+x**2) )",-1,1);
      
      f5->SetParNames("Normalization","Correct Tully Ratio");
      f5->SetParameter(1,0.5);
      f5->FixParameter(0,thish->Integral("bin width"));
      thish->Fit("f5","+L");
      
      float alpha =  f2->GetParameter("alpha");
      float alpha_error = f2->GetParError(1);
      float ratio = f3->GetParameter("R");
      float ratio_error = f3->GetParError(1);
      float tully_ratio = f4->GetParameter("tullyR");
      float tully_ratio_error = f4->GetParError(1);
      float correct_tratio = f5->GetParameter("Correct Tully Ratio");
      float correct_tratio_error = f5->GetParError(1);
      
      f->Close();



     
      x[ipoint] = higgsmass;
      y[ipoint] = alpha;
      ey[ipoint] = alpha_error;
      y1[ipoint] = ratio;
      ey1[ipoint] = ratio_error;
      y2[ipoint] = tully_ratio;
      ey2[ipoint] = tully_ratio_error;
      y3[ipoint] = correct_tratio;
      ey3[ipoint] = correct_tratio_error;

      ipoint++;
      higgsmass+=10;

      filelist >> line; 

    }
 
  TGraphErrors *gr1 = new TGraphErrors(n,x,y,ex,ey);
  
  TGraphErrors *gr2 = new TGraphErrors(n,x,y1,ex,ey1);
  TF1 *f6 = new TF1("f6","((91.1876**4)/(x**4))*(8/((1-((2*(91.1876**2))/(x**2)))**2))",182.3752,500);
  TF1 *f8 = new TF1("f6","1/(-(R_L(x)/R_T(x))/(R_L(x)/R_T(x)-1))",100,182.3752);

  TGraphErrors *gr3 = new TGraphErrors(n,x,y2,ex,ey2);

  TGraphErrors *gr4 = new TGraphErrors(n,x,y3,ex,ey3);
  TF1 *f7 = new TF1("f7","beta(x)",150,500);

  TCanvas *c1 = new TCanvas("c1","alpha",200,10,600,400);

  gr1->SetLineColor(4);
  gr1->SetMarkerStyle(7);
  gr1->SetMarkerSize(3);
  gr1->SetMarkerColor(4);
  gr1->SetTitle("#alpha vs. Higgs mass for incorrect boost 1");
  gr1->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  gr1->GetXaxis()->CenterTitle();
  gr1->GetYaxis()->SetTitle("#alpha"); 
  gr1->GetYaxis()->CenterTitle();
  gr1->Draw("AP");
  
  TCanvas *c2 = new TCanvas("c2","Ratio",200,10,600,400);
  
  c2->cd();
  gr2->SetLineColor(2);
  gr2->SetMarkerStyle(7);
  gr2->SetMarkerSize(3);
  gr2->SetMarkerColor(2);
  gr2->SetTitle("#Gamma_{T}/#Gamma_{L} vs. Higgs Mass (GeV) for incorrect boost 1");
  gr2->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  gr2->GetXaxis()->CenterTitle();
  gr2->GetYaxis()->SetTitle("#Gamma_{T}/#Gamma_{L}");
  gr2->GetYaxis()->CenterTitle();
  gr2->Draw("AP");
  f6->Draw("CSAME");
  f8->Draw("CSAME");

  // TCanvas *c3 = new TCanvas("c3","Tully Ratio",200,10,600,400);

  // c3->cd();
  //gr3->SetMarkerStyle(7);
  //gr3->SetMarkerSize(3);
  //gr3->SetTitle("Tully Ratio");
  //gr3->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  //gr3->Draw("AP");

  TCanvas *c4 = new TCanvas("c4","Correct Tully Ratio",200,10,600,400);
  c4->cd();
  gr4->SetMarkerStyle(7);
  gr4->SetMarkerSize(3);
  gr4->SetTitle("Tully Ratio for incorrect boost 1");
  gr4->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  gr4->GetXaxis()->CenterTitle();
  gr4->GetYaxis()->SetTitle("R = #frac{L-T}{L+T}");
  gr4->GetYaxis()->CenterTitle();
  gr4->Draw("AP");
  f7->Draw("CSAME");
}

---