higgspolarization_incorrect_boost.C

Go to the documentation of this file.

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



void higgspolarization(char *filelistname)

{
  TString filename;
  TFile * f;


  ifstream filelist(filelistname);
  string line;


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

  double x[n],y[n],y1[n],ex[n],ey[n],ey1[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");
      thish->SetTitle("Cos(#theta) of muon1 in Z1 rest frame");
      thish->Rebin(30);
      TF1 *f2 = new TF1("f2","[0]*([1]*(3/4)*(1-x**2)+(3/8)*(1-[1])*(1+x**2))",-1,1);
      TF1 *f3 = new TF1("f3","[0]*( (1/(1+[1]))*(3/4)*(1-x**2) + (3/8)*([1]/(1+[1]))*(1+x**2) )",-1,1);
      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);
      
      
       f2->SetParNames("Normalization","alpha");
       f3->SetParNames("Normalization","R");
       f5->SetParNames("Normalization","Asymmetry");      
 
       f2->SetParameter(1,0.9);
       f3->SetParameter(1,0.1);
       f5->SetParameter(1,0.5);      

       f2->FixParameter(0,thish->Integral("bin width"));
       f3->FixParameter(0,thish->Integral("bin width"));
       f5->FixParameter(0,thish->Integral("bin width"));

       thish->Fit("f2","L+");
       thish->Fit("f3","L+");
       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("Asymmetry");
       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.876**4)/(x**4))*(8/((1-((2*(91.768**2))/(x**2)))**2))",150,500);
  
  //  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 versus Higgs mass");
  gr1->GetXaxis()->SetTitle("Higgs Mass (GeV)"); 
  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("Plot of Ratio of transverse to longitudinal polarization");
  gr2->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  gr2->GetYaxis()->SetTitle("Ratio");
  gr2->Draw("AP");
  f6->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","Tully Ratio",200,10,600,400);
  c4->cd();
  gr4->SetMarkerStyle(7);
  gr4->SetMarkerSize(3);
  gr4->SetTitle("Asymmetry");
  gr4->GetXaxis()->SetTitle("Higgs Mass (GeV)");
  gr4->Draw("AP");
  f7->Draw("CSAME");
}

---