Esempio n. 1
0
float error(TFile * theFile){
  if(!theFile) return 0;
  TH1F * denominator = (TH1F*) (theFile->Get("higgsEta_7"));
  if(!denominator)return 0;
  float eff = efficiency(theFile);
  return sqrt(eff*(1-eff)/denominator->Integral()); 
}
Esempio n. 2
0
Double_t eff(Double_t *x, Double_t *par) {
  Double_t xx = x[0];
  Double_t m0 = par[0];
  Double_t sigma = par[1];
  Double_t alpha = par[2];
  Double_t n = par[3];
  Double_t norm = par[4];
 
  return efficiency(xx,m0,sigma,alpha,n,norm);
}
Esempio n. 3
0
double ratioEfficiencyElec::mcEfficiency(double pt, bool EB) const
{
  double ratio_ = 0 ;
  if (pt<15) return ratio_ ;

  double meanMC,sigmaMC,alphaMC,nMC,normMC ;

  if (EB) {
    meanMC=20.58604584 ; sigmaMC=-1.89456806 ; alphaMC=3.69311772 ; nMC=1.05480046 ; normMC=1.28655181 ;
  } else {
    meanMC=20.15425918 ; sigmaMC=0.75449122 ; alphaMC=1.06027513   ; nMC=1.01106686  ; normMC=7.01956561 ;
  }
  double effMC = efficiency(pt,meanMC,sigmaMC,alphaMC,nMC,normMC) ;

  ratio_ = effMC ;
  return ratio_ ;
}
Esempio n. 4
0
double ratioEfficiencyMu::mcEfficiency(double pt, bool EB) const
{
  double ratio_ = 0 ;
  if (pt<15) return ratio_ ;

  double meanMC,sigmaMC,alphaMC,nMC,normMC ;

  if (EB) {
    meanMC=16.99389526 ; sigmaMC=-0.04080190 ; alphaMC=0.00794730 ; nMC=1.60377906 ; normMC=0.99626161;
  } else {
    meanMC=16.99065795 ; sigmaMC=-0.11993730 ; alphaMC=0.01384991 ; nMC=2.38867304 ; normMC=0.86552275;
  }
  double effMC = efficiency(pt,meanMC,sigmaMC,alphaMC,nMC,normMC) ;

  ratio_ = effMC ;
  return ratio_ ;
}
Esempio n. 5
0
double Recipe::originalGravity() const
{
    double originalGravity = 0;

    foreach(RecipeIngredient* ingredient, ingredients()) {
        GrainIngredient *grainIngredient =
                qobject_cast<GrainIngredient *>(ingredient->ingredient());
        if(grainIngredient) {
            double specificGravity = (grainIngredient->specificGravity() - 1);
            specificGravity *= ingredient->quantity().valueToPound();

            if(grainIngredient->extract()) {
                originalGravity += specificGravity;
            } else {
                originalGravity += specificGravity * efficiency();
            }

            continue;
        }
    }
Esempio n. 6
0
/*	================== getOption =================
 This function reads in the user's desired chose
 of operation. Calls upon other functions to
 perform the procedure.
 Pre		pHeader - pointer to HEAD structure
 Post
 Return
 */
void getOption (HEAD* pHeader)
{
	//	Local Declarations
    char command;
    DATA target;
    DATA* airport = NULL;
	int i;
    
	//	Statements
    while ((command = menu()) != 'Q') {
        switch (command)
        {
            case 'A':
                if (addAirport(pHeader))
                {
                    while (checkHash(pHeader->pHash) == 1) {
                        pHeader->pHash = upsizeHash(pHeader->pHash);
                    }
					printf ("\n Succesfully added data.\n\n");
                }
                break;
            case 'D':
                printf("Enter the airport code: ");
                scanf(" %s", target.arpCode);
                
				if (deleteHash (pHeader, target))
                {
                    while (checkHash(pHeader->pHash) == -1) {
                        pHeader->pHash = downsizeHash(pHeader->pHash);
                    }
					printf ("\n Succesfully deleted data.\n\n");
                }
                
                break;
            case 'F':
                printf("Enter the airport code: ");
                scanf(" %s", target.arpCode);
				// fix sensitive input cases
				for (i = 0; i < strlen(target.arpCode); i++) {
					target.arpCode[i] = toupper(target.arpCode[i]);
				}
                airport = findHash(pHeader->pHash, &target);
                if (airport != NULL) {
                    processScreen(airport);
                }
                else printf("No airport exists\n");
                
                break;
            case 'L':
                printHash(pHeader->pHash);
                break;
            case 'K':
                BST_Traverse(pHeader->pTree, processScreen);
                break;
            case 'P':
                printTree(pHeader->pTree->root, 0);
				printf("\n");
                break;
            case 'W':
				outputFile (pHeader->pHash);
                break;
            case 'E':
				efficiency(pHeader->pHash);
                break;
			case 'H':
				pHeader->pHash = hashDemo(pHeader->pHash);
				break;
            default:
                printf("Invalid choice. Choose again\n");
                break;
        }
    }
    return;
}	// getOption
 Quantity ZoneHVACBaseboardConvectiveElectric_Impl::getEfficiency(bool returnIP) const {
   OptionalDouble value = efficiency();
   OSOptionalQuantity result = getQuantityFromDouble(OS_ZoneHVAC_Baseboard_Convective_ElectricFields::Efficiency, value, returnIP);
   OS_ASSERT(result.isSet());
   return result.get();
 }
Esempio n. 8
0
void higgs_hlt87(){

  float x[16], ex[16];
  float x2[18], ex2[18];
  float x3[19], ex3[19];

  TFile *_fileggh[19];
  _fileggh[0] = TFile::Open("ggh_m60_8TeV/output.root"); 
  _fileggh[1] = TFile::Open("ggh_m65_8TeV/output.root"); 
  _fileggh[2] = TFile::Open("ggh_m70_8TeV/output.root"); 
  _fileggh[3] = TFile::Open("ggh_m75_8TeV/output.root"); 
  _fileggh[4] = TFile::Open("ggh_m80_8TeV/output.root"); 
  _fileggh[5] = TFile::Open("ggh_m85_8TeV/output.root"); 
  _fileggh[6] = TFile::Open("ggh_m90_8TeV/output.root");
  _fileggh[7] = TFile::Open("ggh_m95_8TeV/output.root");
  _fileggh[8] = TFile::Open("ggh_m100_8TeV/output.root");
  _fileggh[9] = TFile::Open("ggh_m105_8TeV/output.root");
  _fileggh[10] = TFile::Open("ggh_m110_8TeV/output.root");
  _fileggh[11] = TFile::Open("ggh_m115_8TeV/output.root");
  _fileggh[12] = TFile::Open("ggh_m120_8TeV/output.root");
  _fileggh[13] = TFile::Open("ggh_m125_8TeV/output.root");
  _fileggh[14] = TFile::Open("ggh_m130_8TeV/output.root");
  _fileggh[15] = TFile::Open("ggh_m135_8TeV/output.root");
  _fileggh[16] = TFile::Open("ggh_m140_8TeV/output.root");
  _fileggh[17] = TFile::Open("ggh_m145_8TeV/output.root");
  _fileggh[18] = TFile::Open("ggh_m150_8TeV/output.root");
  
  TFile *_filetth[19];
  _filetth[0] = TFile::Open("tth_m60_8TeV/output.root"); 
  _filetth[1] = TFile::Open("tth_m65_8TeV/output.root"); 
  _filetth[2] = TFile::Open("tth_m70_8TeV/output.root"); 
  _filetth[3] = TFile::Open("tth_m75_8TeV/output.root"); 
  _filetth[4] = TFile::Open("tth_m80_8TeV/output.root"); 
  _filetth[5] = TFile::Open("tth_m85_8TeV/output.root"); 
  _filetth[6] = TFile::Open("tth_m90_8TeV/output.root");
  _filetth[7] = TFile::Open("tth_m95_8TeV/output.root");
  _filetth[8] = TFile::Open("tth_m100_8TeV/output.root");
  _filetth[9] = TFile::Open("tth_m105_8TeV/output.root");
  _filetth[10] = TFile::Open("tth_m110_8TeV/output.root");
  _filetth[11] = TFile::Open("tth_m115_8TeV/output.root");
  _filetth[12] = TFile::Open("tth_m120_8TeV/output.root");
  _filetth[13] = TFile::Open("tth_m125_8TeV/output.root");
  _filetth[14] = TFile::Open("tth_m130_8TeV/output.root");
  _filetth[15] = TFile::Open("tth_m135_8TeV/output.root");
  _filetth[16] = TFile::Open("tth_m140_8TeV/output.root");
  _filetth[17] = TFile::Open("tth_m145_8TeV/output.root");
  _filetth[18] = TFile::Open("tth_m150_8TeV/output.root");

  TFile *_filevbf[19];
  _filevbf[0] = TFile::Open("vbf_m60_8TeV/output.root"); 
  _filevbf[1] = TFile::Open("vbf_m65_8TeV/output.root"); 
  _filevbf[2] = TFile::Open("vbf_m70_8TeV/output.root"); 
  _filevbf[3] = TFile::Open("vbf_m75_8TeV/output.root"); 
  _filevbf[4] = TFile::Open("vbf_m80_8TeV/output.root"); 
  _filevbf[5] = TFile::Open("vbf_m85_8TeV/output.root"); 
  _filevbf[6] = TFile::Open("vbf_m90_8TeV/output.root");
  _filevbf[7] = TFile::Open("vbf_m95_8TeV/output.root");
  _filevbf[8] = TFile::Open("vbf_m100_8TeV/output.root");
  _filevbf[9] = TFile::Open("vbf_m105_8TeV/output.root");
  _filevbf[10] = TFile::Open("vbf_m110_8TeV/output.root");
  _filevbf[11] = TFile::Open("vbf_m115_8TeV/output.root");
  _filevbf[12] = TFile::Open("vbf_m120_8TeV/output.root");
  _filevbf[13] = TFile::Open("vbf_m125_8TeV/output.root");
  _filevbf[14] = TFile::Open("vbf_m130_8TeV/output.root");
  _filevbf[15] = TFile::Open("vbf_m135_8TeV/output.root");
  _filevbf[16] = TFile::Open("vbf_m140_8TeV/output.root");
  _filevbf[17] = TFile::Open("vbf_m145_8TeV/output.root");
  _filevbf[18] = TFile::Open("vbf_m150_8TeV/output.root");

  TFile *_filewzh[19];
  _filewzh[0] = TFile::Open("wzh_m60_8TeV/output.root"); 
  _filewzh[1] = TFile::Open("wzh_m65_8TeV/output.root"); 
  _filewzh[2] = TFile::Open("wzh_m70_8TeV/output.root"); 
  _filewzh[3] = TFile::Open("wzh_m75_8TeV/output.root"); 
  _filewzh[4] = TFile::Open("wzh_m80_8TeV/output.root"); 
  _filewzh[5] = TFile::Open("wzh_m85_8TeV/output.root"); 
  _filewzh[6] = TFile::Open("wzh_m90_8TeV/output.root");
  _filewzh[7] = TFile::Open("wzh_m95_8TeV/output.root");
  _filewzh[8] = TFile::Open("wzh_m100_8TeV/output.root");
  _filewzh[9] = TFile::Open("wzh_m105_8TeV/output.root");
  _filewzh[10] = TFile::Open("wzh_m110_8TeV/output.root");
  _filewzh[11] = TFile::Open("wzh_m115_8TeV/output.root");
  _filewzh[12] = TFile::Open("wzh_m120_8TeV/output.root");
  _filewzh[13] = TFile::Open("wzh_m125_8TeV/output.root");
  _filewzh[14] = TFile::Open("wzh_m130_8TeV/output.root");
  _filewzh[15] = TFile::Open("wzh_m135_8TeV/output.root");
  _filewzh[16] = TFile::Open("wzh_m140_8TeV/output.root");
  _filewzh[17] = TFile::Open("wzh_m140_8TeV/output.root");
  _filewzh[18] = TFile::Open("wzh_m145_8TeV/output.root");

  float avEFFggh[19];
  float avEFFtth[19];
  float avEFFwzh[19];
  float avEFFvbf[19];
  float avERRggh[19];
  float avERRtth[19];
  float avERRwzh[19];
  float avERRvbf[19];

  float mass_value=55.;
  
  for(int k=0;k<19;k++){
    mass_value=mass_value+5.;
    x3[k]  = mass_value;
    ex3[k] = 0;
    avEFFggh[k] = efficiency(_fileggh[k]);
    avEFFtth[k] = efficiency(_filetth[k]);
    avEFFvbf[k] = efficiency(_filevbf[k]);
    avEFFwzh[k] = efficiency(_filewzh[k]);

    avERRggh[k] = error(_fileggh[k]);
    avERRwzh[k] = error(_filewzh[k]);
    avERRtth[k] = error(_filetth[k]);
    avERRvbf[k] = error(_filevbf[k]);
  }
 
  cout << x3[0] <<" : " << avEFFggh[0] << "------" << avERRggh[0] << endl;
  cout << x3[1] <<" : " << avEFFggh[1] << "------" << avERRggh[1] << endl;
  cout << x3[2] <<" : " << avEFFggh[2] << "------" << avERRggh[2] << endl;
  cout << x3[3] <<" : " << avEFFggh[3] << "------" << avERRggh[3] << endl;
    
  TGraphErrors * ggh = new TGraphErrors(19,x3,avEFFggh,ex3,avERRggh);
  TGraphErrors * vbf = new TGraphErrors(19,x3,avEFFvbf,ex3,avERRvbf);
  TGraphErrors * wzh = new TGraphErrors(19,x3,avEFFwzh,ex3,avERRwzh);
  TGraphErrors * tth = new TGraphErrors(19,x3,avEFFtth,ex3,avERRtth);

  ggh->SetMarkerStyle(22);ggh->SetMarkerSize(2.0); 
  vbf->SetMarkerStyle(23);vbf->SetMarkerSize(2.0); 
  wzh->SetMarkerStyle(20);wzh->SetMarkerSize(2.0); 
  tth->SetMarkerStyle(21);tth->SetMarkerSize(2.0); 

  ggh->SetMarkerColor(kRed);
  vbf->SetMarkerColor(kBlue);
  wzh->SetMarkerColor(kGreen);
  tth->SetMarkerColor(kBlack);

  tth->GetXaxis()->SetTitle("M_{H} (GeV)");  tth->GetYaxis()->SetTitle("#epsilon");
  ggh->GetXaxis()->SetTitle("M_{H} (GeV)");  ggh->GetYaxis()->SetTitle("#epsilon");
  vbf->GetXaxis()->SetTitle("M_{H} (GeV)");  vbf->GetYaxis()->SetTitle("#epsilon");
  wzh->GetXaxis()->SetTitle("M_{H} (GeV)");  wzh->GetYaxis()->SetTitle("#epsilon");
  
  TCanvas *C1 = new TCanvas("C1","Global Efficiency",200,10,1200,800);
  C1->SetFillColor(10);
  C1->Draw();
  C1->cd();

  float min=0.0;
  float max=1.0;
  ggh->SetMinimum(min);ggh->SetMaximum(max);
  vbf->SetMinimum(min);vbf->SetMaximum(max);
  wzh->SetMinimum(min);wzh->SetMaximum(max);
  tth->SetMinimum(min);tth->SetMaximum(max);
  
  ggh->Draw("AP"); 
  vbf->Draw("Psame"); 
  wzh->Draw("Psame"); 
  tth->Draw("Psame"); 

  C1->Update();

  TLegend *leg = new TLegend(0.6,0.35,0.7,0.2);
  //leg->AddEntry(glui,"gluino split-SUSY","lp");                                                                                                               
  //leg->AddEntry(stop,"stop MSSM","lp");                                                                                                                       
  leg->AddEntry(ggh,"ggh","p");
  leg->AddEntry(vbf,"vbf","p");
  leg->AddEntry(tth,"tth","p");
  leg->AddEntry(wzh,"wzh","p");
  leg->SetFillColor(0);
  leg->Draw("same");
  
  C1->SaveAs("summary87.png");

  min=0.9;
  max=1.0;
  ggh->SetMinimum(min);ggh->SetMaximum(max);
  vbf->SetMinimum(min);vbf->SetMaximum(max);
  wzh->SetMinimum(min);wzh->SetMaximum(max);
  tth->SetMinimum(min);tth->SetMaximum(max);
  
  ggh->Draw("AP"); 
  vbf->Draw("Psame"); 
  wzh->Draw("Psame"); 
  tth->Draw("Psame"); 

  C1->Update();

  TLegend *leg = new TLegend(0.6,0.35,0.7,0.2);
  //leg->AddEntry(glui,"gluino split-SUSY","lp");                                                                                                               
  //leg->AddEntry(stop,"stop MSSM","lp");                                                                                                                       
  leg->AddEntry(ggh,"ggh","p");
  leg->AddEntry(vbf,"vbf","p");
  leg->AddEntry(tth,"tth","p");
  leg->AddEntry(wzh,"wzh","p");
  leg->SetFillColor(0);
  leg->Draw("same");
  
  C1->SaveAs("summary87_0p9.png");

  min=0.7;
  max=1.0;
  
  ggh->SetMinimum(min);ggh->SetMaximum(max);
  vbf->SetMinimum(min);vbf->SetMaximum(max);
  wzh->SetMinimum(min);wzh->SetMaximum(max);
  tth->SetMinimum(min);tth->SetMaximum(max);
  
  ggh->Draw("AP"); 
  vbf->Draw("Psame"); 
  wzh->Draw("Psame"); 
  tth->Draw("Psame"); 

  C1->Update();

  TLegend *leg = new TLegend(0.6,0.35,0.7,0.2);
  leg->AddEntry(ggh,"ggh","p");
  leg->AddEntry(vbf,"vbf","p");
  leg->AddEntry(tth,"tth","p");
  leg->AddEntry(wzh,"wzh","p");
  leg->SetFillColor(0);
  leg->Draw("same");
  
  C1->SaveAs("summary87_0p7.png");

  exit(0);
}