Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-05 08:13:11 

0001 #include <iostream>
0002 #include <TFile.h>
0003 #include <TH1D.h>
0004 #include "/sphenix/user/hanpuj/CaloDataAna24_skimmed/src/draw_template.C"
0005 #include "unfold_Def.h"
0006 
0007 void get_purityefficiency() {
0008   // Read Files
0009   TFile *f_in = new TFile("output_sim.root", "READ");
0010   if (!f_in) {
0011     std::cout << "Error: cannot open output_sim.root" << std::endl;
0012     return;
0013   }
0014 
0015   // Read histograms
0016   TH1D* h_truth_calib_dijet = (TH1D*)f_in->Get("h_truth_calib_dijet");
0017   TH1D* h_measure_calib_dijet = (TH1D*)f_in->Get("h_measure_calib_dijet");
0018   TH2D* h_respmatrix_calib_dijet = (TH2D*)f_in->Get("h_respmatrix_calib_dijet");
0019 
0020   TH1D* h_truth_calib_dijet_jesdown = (TH1D*)f_in->Get("h_truth_calib_dijet_jesdown");
0021   TH1D* h_measure_calib_dijet_jesdown = (TH1D*)f_in->Get("h_measure_calib_dijet_jesdown");
0022   TH2D* h_respmatrix_calib_dijet_jesdown = (TH2D*)f_in->Get("h_respmatrix_calib_dijet_jesdown");
0023 
0024   TH1D* h_truth_calib_dijet_jesup = (TH1D*)f_in->Get("h_truth_calib_dijet_jesup");
0025   TH1D* h_measure_calib_dijet_jesup = (TH1D*)f_in->Get("h_measure_calib_dijet_jesup");
0026   TH2D* h_respmatrix_calib_dijet_jesup = (TH2D*)f_in->Get("h_respmatrix_calib_dijet_jesup");
0027 
0028   TH1D* h_truth_calib_dijet_jerdown = (TH1D*)f_in->Get("h_truth_calib_dijet_jerdown");
0029   TH1D* h_measure_calib_dijet_jerdown = (TH1D*)f_in->Get("h_measure_calib_dijet_jerdown");
0030   TH2D* h_respmatrix_calib_dijet_jerdown = (TH2D*)f_in->Get("h_respmatrix_calib_dijet_jerdown");
0031 
0032   TH1D* h_truth_calib_dijet_jerup = (TH1D*)f_in->Get("h_truth_calib_dijet_jerup");
0033   TH1D* h_measure_calib_dijet_jerup = (TH1D*)f_in->Get("h_measure_calib_dijet_jerup");
0034   TH2D* h_respmatrix_calib_dijet_jerup = (TH2D*)f_in->Get("h_respmatrix_calib_dijet_jerup");
0035   
0036   TH1D* h_truth_calib_frac = (TH1D*)f_in->Get("h_truth_calib_frac");
0037   TH1D* h_measure_calib_frac = (TH1D*)f_in->Get("h_measure_calib_frac");
0038   TH2D* h_respmatrix_calib_frac = (TH2D*)f_in->Get("h_respmatrix_calib_frac");
0039 
0040   TH1D* h_truth_calib_frac_jesdown = (TH1D*)f_in->Get("h_truth_calib_frac_jesdown");
0041   TH1D* h_measure_calib_frac_jesdown = (TH1D*)f_in->Get("h_measure_calib_frac_jesdown");
0042   TH2D* h_respmatrix_calib_frac_jesdown = (TH2D*)f_in->Get("h_respmatrix_calib_frac_jesdown");
0043 
0044   TH1D* h_truth_calib_frac_jesup = (TH1D*)f_in->Get("h_truth_calib_frac_jesup");
0045   TH1D* h_measure_calib_frac_jesup = (TH1D*)f_in->Get("h_measure_calib_frac_jesup");
0046   TH2D* h_respmatrix_calib_frac_jesup = (TH2D*)f_in->Get("h_respmatrix_calib_frac_jesup");
0047 
0048   TH1D* h_truth_calib_frac_jerdown = (TH1D*)f_in->Get("h_truth_calib_frac_jerdown");
0049   TH1D* h_measure_calib_frac_jerdown = (TH1D*)f_in->Get("h_measure_calib_frac_jerdown");
0050   TH2D* h_respmatrix_calib_frac_jerdown = (TH2D*)f_in->Get("h_respmatrix_calib_frac_jerdown");
0051 
0052   TH1D* h_truth_calib_frac_jerup = (TH1D*)f_in->Get("h_truth_calib_frac_jerup");
0053   TH1D* h_measure_calib_frac_jerup = (TH1D*)f_in->Get("h_measure_calib_frac_jerup");
0054   TH2D* h_respmatrix_calib_frac_jerup = (TH2D*)f_in->Get("h_respmatrix_calib_frac_jerup");
0055 
0056   //TH1D* h_truth_calib_frac_half1 = (TH1D*)f_in->Get("h_truth_calib_frac_half1");
0057   //TH1D* h_measure_calib_frac_half1 = (TH1D*)f_in->Get("h_measure_calib_frac_half1");
0058   //TH2D* h_respmatrix_calib_frac_half1 = (TH2D*)f_in->Get("h_respmatrix_calib_frac_half1");
0059 
0060   // Form purity and efficiency histograms
0061   TH1D* h_purity_calib_dijet = (TH1D*)h_respmatrix_calib_dijet->ProjectionX("h_purity_calib_dijet");
0062   TH1D* h_efficiency_calib_dijet = (TH1D*)h_respmatrix_calib_dijet->ProjectionY("h_efficiency_calib_dijet");
0063 
0064   TH1D* h_purity_calib_dijet_jesdown = (TH1D*)h_respmatrix_calib_dijet_jesdown->ProjectionX("h_purity_calib_dijet_jesdown");
0065   TH1D* h_efficiency_calib_dijet_jesdown = (TH1D*)h_respmatrix_calib_dijet_jesdown->ProjectionY("h_efficiency_calib_dijet_jesdown");
0066 
0067   TH1D* h_purity_calib_dijet_jesup = (TH1D*)h_respmatrix_calib_dijet_jesup->ProjectionX("h_purity_calib_dijet_jesup");
0068   TH1D* h_efficiency_calib_dijet_jesup = (TH1D*)h_respmatrix_calib_dijet_jesup->ProjectionY("h_efficiency_calib_dijet_jesup");
0069 
0070   TH1D* h_purity_calib_dijet_jerdown = (TH1D*)h_respmatrix_calib_dijet_jerdown->ProjectionX("h_purity_calib_dijet_jerdown");
0071   TH1D* h_efficiency_calib_dijet_jerdown = (TH1D*)h_respmatrix_calib_dijet_jerdown->ProjectionY("h_efficiency_calib_dijet_jerdown");
0072 
0073   TH1D* h_purity_calib_dijet_jerup = (TH1D*)h_respmatrix_calib_dijet_jerup->ProjectionX("h_purity_calib_dijet_jerup");
0074   TH1D* h_efficiency_calib_dijet_jerup = (TH1D*)h_respmatrix_calib_dijet_jerup->ProjectionY("h_efficiency_calib_dijet_jerup");
0075 
0076   TH1D* h_purity_calib_frac = (TH1D*)h_respmatrix_calib_frac->ProjectionX("h_purity_calib_frac");
0077   TH1D* h_efficiency_calib_frac = (TH1D*)h_respmatrix_calib_frac->ProjectionY("h_efficiency_calib_frac");
0078 
0079   TH1D* h_purity_calib_frac_jesdown = (TH1D*)h_respmatrix_calib_frac_jesdown->ProjectionX("h_purity_calib_frac_jesdown");
0080   TH1D* h_efficiency_calib_frac_jesdown = (TH1D*)h_respmatrix_calib_frac_jesdown->ProjectionY("h_efficiency_calib_frac_jesdown");
0081 
0082   TH1D* h_purity_calib_frac_jesup = (TH1D*)h_respmatrix_calib_frac_jesup->ProjectionX("h_purity_calib_frac_jesup");
0083   TH1D* h_efficiency_calib_frac_jesup = (TH1D*)h_respmatrix_calib_frac_jesup->ProjectionY("h_efficiency_calib_frac_jesup");
0084 
0085   TH1D* h_purity_calib_frac_jerdown = (TH1D*)h_respmatrix_calib_frac_jerdown->ProjectionX("h_purity_calib_frac_jerdown");
0086   TH1D* h_efficiency_calib_frac_jerdown = (TH1D*)h_respmatrix_calib_frac_jerdown->ProjectionY("h_efficiency_calib_frac_jerdown");
0087 
0088   TH1D* h_purity_calib_frac_jerup = (TH1D*)h_respmatrix_calib_frac_jerup->ProjectionX("h_purity_calib_frac_jerup");
0089   TH1D* h_efficiency_calib_frac_jerup = (TH1D*)h_respmatrix_calib_frac_jerup->ProjectionY("h_efficiency_calib_frac_jerup");
0090 
0091   //TH1D* h_purity_calib_frac_half1 = (TH1D*)h_respmatrix_calib_frac_half1->ProjectionX("h_purity_calib_frac_half1");
0092   //TH1D* h_efficiency_calib_frac_half1 = (TH1D*)h_respmatrix_calib_frac_half1->ProjectionY("h_efficiency_calib_frac_half1");
0093 
0094   // Calculate purity and efficiency
0095   h_purity_calib_dijet->Divide(h_purity_calib_dijet, h_measure_calib_dijet, 1, 1, "B");
0096   h_efficiency_calib_dijet->Divide(h_efficiency_calib_dijet, h_truth_calib_dijet, 1, 1, "B");
0097 
0098   h_purity_calib_dijet_jesdown->Divide(h_purity_calib_dijet_jesdown, h_measure_calib_dijet_jesdown, 1, 1, "B");
0099   h_efficiency_calib_dijet_jesdown->Divide(h_efficiency_calib_dijet_jesdown, h_truth_calib_dijet_jesdown, 1, 1, "B");
0100 
0101   h_purity_calib_dijet_jesup->Divide(h_purity_calib_dijet_jesup, h_measure_calib_dijet_jesup, 1, 1, "B");
0102   h_efficiency_calib_dijet_jesup->Divide(h_efficiency_calib_dijet_jesup, h_truth_calib_dijet_jesup, 1, 1, "B");
0103 
0104   h_purity_calib_dijet_jerdown->Divide(h_purity_calib_dijet_jerdown, h_measure_calib_dijet_jerdown, 1, 1, "B");
0105   h_efficiency_calib_dijet_jerdown->Divide(h_efficiency_calib_dijet_jerdown, h_truth_calib_dijet_jerdown, 1, 1, "B");
0106 
0107   h_purity_calib_dijet_jerup->Divide(h_purity_calib_dijet_jerup, h_measure_calib_dijet_jerup, 1, 1, "B");
0108   h_efficiency_calib_dijet_jerup->Divide(h_efficiency_calib_dijet_jerup, h_truth_calib_dijet_jerup, 1, 1, "B");
0109 
0110   h_purity_calib_frac->Divide(h_purity_calib_frac, h_measure_calib_frac, 1, 1, "B");
0111   h_efficiency_calib_frac->Divide(h_efficiency_calib_frac, h_truth_calib_frac, 1, 1, "B");
0112 
0113   h_purity_calib_frac_jesdown->Divide(h_purity_calib_frac_jesdown, h_measure_calib_frac_jesdown, 1, 1, "B");
0114   h_efficiency_calib_frac_jesdown->Divide(h_efficiency_calib_frac_jesdown, h_truth_calib_frac_jesdown, 1, 1, "B");
0115 
0116   h_purity_calib_frac_jesup->Divide(h_purity_calib_frac_jesup, h_measure_calib_frac_jesup, 1, 1, "B");
0117   h_efficiency_calib_frac_jesup->Divide(h_efficiency_calib_frac_jesup, h_truth_calib_frac_jesup, 1, 1, "B");
0118 
0119   h_purity_calib_frac_jerdown->Divide(h_purity_calib_frac_jerdown, h_measure_calib_frac_jerdown, 1, 1, "B");
0120   h_efficiency_calib_frac_jerdown->Divide(h_efficiency_calib_frac_jerdown, h_truth_calib_frac_jerdown, 1, 1, "B");
0121 
0122   h_purity_calib_frac_jerup->Divide(h_purity_calib_frac_jerup, h_measure_calib_frac_jerup, 1, 1, "B");
0123   h_efficiency_calib_frac_jerup->Divide(h_efficiency_calib_frac_jerup, h_truth_calib_frac_jerup, 1, 1, "B");
0124 
0125   //h_purity_calib_frac_half1->Divide(h_purity_calib_frac_half1, h_measure_calib_frac_half1, 1, 1, "B");
0126   //h_efficiency_calib_frac_half1->Divide(h_efficiency_calib_frac_half1, h_truth_calib_frac_half1, 1, 1, "B");
0127 
0128   // Write histograms
0129   TFile *f_out = new TFile("output_purityefficiency.root", "RECREATE");
0130   f_out->cd();
0131   h_purity_calib_dijet->Write();
0132   h_efficiency_calib_dijet->Write();
0133   
0134   h_purity_calib_dijet_jesdown->Write();
0135   h_efficiency_calib_dijet_jesdown->Write();
0136 
0137   h_purity_calib_dijet_jesup->Write();
0138   h_efficiency_calib_dijet_jesup->Write();
0139 
0140   h_purity_calib_dijet_jerdown->Write();
0141   h_efficiency_calib_dijet_jerdown->Write();
0142 
0143   h_purity_calib_dijet_jerup->Write();
0144   h_efficiency_calib_dijet_jerup->Write();
0145 
0146   h_purity_calib_frac->Write();
0147   h_efficiency_calib_frac->Write();
0148 
0149   h_purity_calib_frac_jesdown->Write();
0150   h_efficiency_calib_frac_jesdown->Write();
0151 
0152   h_purity_calib_frac_jesup->Write();
0153   h_efficiency_calib_frac_jesup->Write();
0154 
0155   h_purity_calib_frac_jerdown->Write();
0156   h_efficiency_calib_frac_jerdown->Write();
0157 
0158   h_purity_calib_frac_jerup->Write();
0159   h_efficiency_calib_frac_jerup->Write();
0160 
0161   //h_purity_calib_frac_half1->Write();
0162   //h_efficiency_calib_frac_half1->Write();
0163   f_out->Close();
0164 
0165   // Apply purity correction.
0166   TFile* f_data = new TFile("output_data.root", "READ");
0167   TFile* f_dataout = new TFile("output_data_puritycorr.root", "RECREATE");
0168 
0169   // Read histograms
0170   TH1D* h_calibjet_pt_dijet_eff = (TH1D*)f_data->Get("h_calibjet_pt_dijet_eff");
0171   TH1D* h_calibjet_pt_dijet_effdown = (TH1D*)f_data->Get("h_calibjet_pt_dijet_effdown");
0172   TH1D* h_calibjet_pt_dijet_effup = (TH1D*)f_data->Get("h_calibjet_pt_dijet_effup");
0173   TH1D* h_calibjet_pt_frac_eff = (TH1D*)f_data->Get("h_calibjet_pt_frac_eff");
0174   TH1D* h_calibjet_pt_frac_effdown = (TH1D*)f_data->Get("h_calibjet_pt_frac_effdown");
0175   TH1D* h_calibjet_pt_frac_effup = (TH1D*)f_data->Get("h_calibjet_pt_frac_effup");
0176 
0177   // Histograms for purity-corrected jet pT.
0178   TH1D* h_calibjet_pt_puritycorr_dijet = (TH1D*)h_calibjet_pt_dijet_eff->Clone("h_calibjet_pt_puritycorr_dijet");
0179   TH1D* h_calibjet_pt_puritycorr_dijet_effdown = (TH1D*)h_calibjet_pt_dijet_effdown->Clone("h_calibjet_pt_puritycorr_dijet_effdown");
0180   TH1D* h_calibjet_pt_puritycorr_dijet_effup = (TH1D*)h_calibjet_pt_dijet_effup->Clone("h_calibjet_pt_puritycorr_dijet_effup");
0181   TH1D* h_calibjet_pt_puritycorr_dijet_jesdown = (TH1D*)h_calibjet_pt_dijet_eff->Clone("h_calibjet_pt_puritycorr_dijet_jesdown");
0182   TH1D* h_calibjet_pt_puritycorr_dijet_jesup = (TH1D*)h_calibjet_pt_dijet_eff->Clone("h_calibjet_pt_puritycorr_dijet_jesup");
0183   TH1D* h_calibjet_pt_puritycorr_dijet_jerdown = (TH1D*)h_calibjet_pt_dijet_eff->Clone("h_calibjet_pt_puritycorr_dijet_jerdown");
0184   TH1D* h_calibjet_pt_puritycorr_dijet_jerup = (TH1D*)h_calibjet_pt_dijet_eff->Clone("h_calibjet_pt_puritycorr_dijet_jerup");
0185 
0186   TH1D* h_calibjet_pt_puritycorr_frac = (TH1D*)h_calibjet_pt_frac_eff->Clone("h_calibjet_pt_puritycorr_frac");
0187   TH1D* h_calibjet_pt_puritycorr_frac_effdown = (TH1D*)h_calibjet_pt_frac_effdown->Clone("h_calibjet_pt_puritycorr_frac_effdown");
0188   TH1D* h_calibjet_pt_puritycorr_frac_effup = (TH1D*)h_calibjet_pt_frac_effup->Clone("h_calibjet_pt_puritycorr_frac_effup");
0189   TH1D* h_calibjet_pt_puritycorr_frac_jesdown = (TH1D*)h_calibjet_pt_frac_eff->Clone("h_calibjet_pt_puritycorr_frac_jesdown");
0190   TH1D* h_calibjet_pt_puritycorr_frac_jesup = (TH1D*)h_calibjet_pt_frac_eff->Clone("h_calibjet_pt_puritycorr_frac_jesup");
0191   TH1D* h_calibjet_pt_puritycorr_frac_jerdown = (TH1D*)h_calibjet_pt_frac_eff->Clone("h_calibjet_pt_puritycorr_frac_jerdown");
0192   TH1D* h_calibjet_pt_puritycorr_frac_jerup = (TH1D*)h_calibjet_pt_frac_eff->Clone("h_calibjet_pt_puritycorr_frac_jerup");
0193 
0194   h_calibjet_pt_puritycorr_dijet->Divide(h_purity_calib_dijet);
0195   h_calibjet_pt_puritycorr_dijet_effdown->Divide(h_purity_calib_dijet);
0196   h_calibjet_pt_puritycorr_dijet_effup->Divide(h_purity_calib_dijet);
0197   h_calibjet_pt_puritycorr_dijet_jesdown->Divide(h_purity_calib_dijet_jesdown);
0198   h_calibjet_pt_puritycorr_dijet_jesup->Divide(h_purity_calib_dijet_jesup);
0199   h_calibjet_pt_puritycorr_dijet_jerdown->Divide(h_purity_calib_dijet_jerdown);
0200   h_calibjet_pt_puritycorr_dijet_jerup->Divide(h_purity_calib_dijet_jerup);
0201 
0202   h_calibjet_pt_puritycorr_frac->Divide(h_purity_calib_frac);
0203   h_calibjet_pt_puritycorr_frac_effdown->Divide(h_purity_calib_frac);
0204   h_calibjet_pt_puritycorr_frac_effup->Divide(h_purity_calib_frac);
0205   h_calibjet_pt_puritycorr_frac_jesdown->Divide(h_purity_calib_frac_jesdown);
0206   h_calibjet_pt_puritycorr_frac_jesup->Divide(h_purity_calib_frac_jesup);
0207   h_calibjet_pt_puritycorr_frac_jerdown->Divide(h_purity_calib_frac_jerdown);
0208   h_calibjet_pt_puritycorr_frac_jerup->Divide(h_purity_calib_frac_jerup);
0209 
0210   // Write histograms
0211   f_dataout->cd();
0212   h_calibjet_pt_puritycorr_dijet->Write();
0213   h_calibjet_pt_puritycorr_dijet_effdown->Write();
0214   h_calibjet_pt_puritycorr_dijet_effup->Write();
0215   h_calibjet_pt_puritycorr_dijet_jesdown->Write();
0216   h_calibjet_pt_puritycorr_dijet_jesup->Write();
0217   h_calibjet_pt_puritycorr_dijet_jerdown->Write();
0218   h_calibjet_pt_puritycorr_dijet_jerup->Write();
0219 
0220   h_calibjet_pt_puritycorr_frac->Write();
0221   h_calibjet_pt_puritycorr_frac_effdown->Write();
0222   h_calibjet_pt_puritycorr_frac_effup->Write();
0223   h_calibjet_pt_puritycorr_frac_jesdown->Write();
0224   h_calibjet_pt_puritycorr_frac_jesup->Write();
0225   h_calibjet_pt_puritycorr_frac_jerdown->Write();
0226   h_calibjet_pt_puritycorr_frac_jerup->Write();
0227   f_dataout->Close();
0228 
0229   // Draw purity and efficiency.
0230   TH1D* h_purity_calib_dijet_proj = (TH1D*)h_respmatrix_calib_dijet->ProjectionX("h_purity_calib_dijet_proj");
0231   TH1D* h_efficiency_calib_dijet_proj = (TH1D*)h_respmatrix_calib_dijet->ProjectionY("h_efficiency_calib_dijet_proj");
0232 
0233   TH1D* h_purity_calib_dijet_jesdown_proj = (TH1D*)h_respmatrix_calib_dijet_jesdown->ProjectionX("h_purity_calib_dijet_jesdown_proj");
0234   TH1D* h_efficiency_calib_dijet_jesdown_proj = (TH1D*)h_respmatrix_calib_dijet_jesdown->ProjectionY("h_efficiency_calib_dijet_jesdown_proj");
0235 
0236   TH1D* h_purity_calib_dijet_jesup_proj = (TH1D*)h_respmatrix_calib_dijet_jesup->ProjectionX("h_purity_calib_dijet_jesup_proj");
0237   TH1D* h_efficiency_calib_dijet_jesup_proj = (TH1D*)h_respmatrix_calib_dijet_jesup->ProjectionY("h_efficiency_calib_dijet_jesup_proj");
0238 
0239   TH1D* h_purity_calib_dijet_jerdown_proj = (TH1D*)h_respmatrix_calib_dijet_jerdown->ProjectionX("h_purity_calib_dijet_jerdown_proj");
0240   TH1D* h_efficiency_calib_dijet_jerdown_proj = (TH1D*)h_respmatrix_calib_dijet_jerdown->ProjectionY("h_efficiency_calib_dijet_jerdown_proj");
0241 
0242   TH1D* h_purity_calib_dijet_jerup_proj = (TH1D*)h_respmatrix_calib_dijet_jerup->ProjectionX("h_purity_calib_dijet_jerup_proj");
0243   TH1D* h_efficiency_calib_dijet_jerup_proj = (TH1D*)h_respmatrix_calib_dijet_jerup->ProjectionY("h_efficiency_calib_dijet_jerup_proj");
0244 
0245   TH1D* h_purity_calib_frac_proj = (TH1D*)h_respmatrix_calib_frac->ProjectionX("h_purity_calib_frac_proj");
0246   TH1D* h_efficiency_calib_frac_proj = (TH1D*)h_respmatrix_calib_frac->ProjectionY("h_efficiency_calib_frac_proj");
0247 
0248   TH1D* h_purity_calib_frac_jesdown_proj = (TH1D*)h_respmatrix_calib_frac_jesdown->ProjectionX("h_purity_calib_frac_jesdown_proj");
0249   TH1D* h_efficiency_calib_frac_jesdown_proj = (TH1D*)h_respmatrix_calib_frac_jesdown->ProjectionY("h_efficiency_calib_frac_jesdown_proj");
0250 
0251   TH1D* h_purity_calib_frac_jesup_proj = (TH1D*)h_respmatrix_calib_frac_jesup->ProjectionX("h_purity_calib_frac_jesup_proj");
0252   TH1D* h_efficiency_calib_frac_jesup_proj = (TH1D*)h_respmatrix_calib_frac_jesup->ProjectionY("h_efficiency_calib_frac_jesup_proj");
0253 
0254   TH1D* h_purity_calib_frac_jerdown_proj = (TH1D*)h_respmatrix_calib_frac_jerdown->ProjectionX("h_purity_calib_frac_jerdown_proj");
0255   TH1D* h_efficiency_calib_frac_jerdown_proj = (TH1D*)h_respmatrix_calib_frac_jerdown->ProjectionY("h_efficiency_calib_frac_jerdown_proj");
0256 
0257   TH1D* h_purity_calib_frac_jerup_proj = (TH1D*)h_respmatrix_calib_frac_jerup->ProjectionX("h_purity_calib_frac_jerup_proj");
0258   TH1D* h_efficiency_calib_frac_jerup_proj = (TH1D*)h_respmatrix_calib_frac_jerup->ProjectionY("h_efficiency_calib_frac_jerup_proj");
0259 
0260   std::vector<TH1F*> h_input;
0261   std::vector<int> color;
0262   std::vector<int> markerstyle;
0263   std::vector<std::string> text;
0264   std::vector<std::string> legend;
0265 
0266   h_input.push_back((TH1F*)h_measure_calib_dijet);
0267   h_input.push_back((TH1F*)h_purity_calib_dijet_proj);
0268   color.push_back(kBlack);
0269   color.push_back(kRed);
0270   markerstyle.push_back(20);
0271   markerstyle.push_back(20);
0272   text.push_back("#bf{#it{sPHENIX}} Simulation");
0273   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0274   text.push_back("anti-k_{T}#it{R} = 0.4");
0275   text.push_back("|#eta^{jet}| < 0.7");
0276   text.push_back("#DeltaR_{matching} < 0.3");
0277   legend.push_back("Reco jet spectrum");
0278   legend.push_back("Matched reco jet spectrum");
0279   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0280                               false, 10, true,
0281                               true, 15, 72, false,
0282                               false, 0, 0.5, true,
0283                               true, 0.75, 1.1,
0284                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0285                               true, text, 0.55, 0.9, 0.04,
0286                               true, legend, 0.58, 0.65, 0.04,
0287                               "figure/purity_calib_dijet.pdf");
0288   h_input.clear();
0289   color.clear();
0290   markerstyle.clear();
0291   text.clear();
0292   legend.clear();
0293 
0294   h_input.push_back((TH1F*)h_measure_calib_dijet_jesdown);
0295   h_input.push_back((TH1F*)h_purity_calib_dijet_jesdown_proj);
0296   color.push_back(kBlack);
0297   color.push_back(kRed);
0298   markerstyle.push_back(20);
0299   markerstyle.push_back(20);
0300   text.push_back("#bf{#it{sPHENIX}} Simulation");
0301   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0302   text.push_back("anti-k_{T}#it{R} = 0.4");
0303   text.push_back("|#eta^{jet}| < 0.7");
0304   text.push_back("#DeltaR_{matching} < 0.3");
0305   text.push_back("With JES calibration");
0306   legend.push_back("Reco jet spectrum");
0307   legend.push_back("Matched reco jet spectrum");
0308   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0309                               false, 10, true,
0310                               true, 15, 72, false,
0311                               false, 0, 0.5, true,
0312                               true, 0.75, 1.1,
0313                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0314                               true, text, 0.55, 0.9, 0.04,
0315                               true, legend, 0.58, 0.6, 0.04,
0316                               "figure/purity_calib_dijet_jesdown.pdf");
0317   h_input.clear();
0318   color.clear();
0319   markerstyle.clear();
0320   text.clear();
0321   legend.clear();
0322 
0323   h_input.push_back((TH1F*)h_measure_calib_dijet_jesup);
0324   h_input.push_back((TH1F*)h_purity_calib_dijet_jesup_proj);
0325   color.push_back(kBlack);
0326   color.push_back(kRed);
0327   markerstyle.push_back(20);
0328   markerstyle.push_back(20);
0329   text.push_back("#bf{#it{sPHENIX}} Simulation");
0330   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0331   text.push_back("anti-k_{T}#it{R} = 0.4");
0332   text.push_back("|#eta^{jet}| < 0.7");
0333   text.push_back("#DeltaR_{matching} < 0.3");
0334   legend.push_back("Truth jet spectrum");
0335   legend.push_back("Matched truth jet spectrum");
0336   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0337                               false, 10, true,
0338                               true, 15, 72, false,
0339                               false, 0, 0.5, true,
0340                               true, 0.75, 1.1,
0341                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0342                               true, text, 0.55, 0.9, 0.04,
0343                               true, legend, 0.58, 0.65, 0.04,
0344                               "figure/purity_calib_dijet_jesup.pdf");
0345   h_input.clear();
0346   color.clear();
0347   markerstyle.clear();
0348   text.clear();
0349   legend.clear();
0350 
0351   h_input.push_back((TH1F*)h_measure_calib_dijet_jerdown);
0352   h_input.push_back((TH1F*)h_purity_calib_dijet_jerdown_proj);
0353   color.push_back(kBlack);
0354   color.push_back(kRed);
0355   markerstyle.push_back(20);
0356   markerstyle.push_back(20);
0357   text.push_back("#bf{#it{sPHENIX}} Simulation");
0358   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0359   text.push_back("anti-k_{T}#it{R} = 0.4");
0360   text.push_back("|#eta^{jet}| < 0.7");
0361   text.push_back("#DeltaR_{matching} < 0.3");
0362   text.push_back("With JES calibration");
0363   legend.push_back("Truth jet spectrum");
0364   legend.push_back("Matched truth jet spectrum");
0365   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0366                               false, 10, true,
0367                               true, 15, 72, false,
0368                               false, 0, 0.5, true,
0369                               true, 0.75, 1.1,
0370                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0371                               true, text, 0.55, 0.9, 0.04,
0372                               true, legend, 0.58, 0.6, 0.04,
0373                               "figure/purity_calib_dijet_jerdown.pdf");
0374   h_input.clear();
0375   color.clear();
0376   markerstyle.clear();
0377   text.clear();
0378   legend.clear();
0379 
0380   h_input.push_back((TH1F*)h_measure_calib_dijet_jerup);
0381   h_input.push_back((TH1F*)h_purity_calib_dijet_jerup_proj);
0382   color.push_back(kBlack);
0383   color.push_back(kRed);
0384   markerstyle.push_back(20);
0385   markerstyle.push_back(20);
0386   text.push_back("#bf{#it{sPHENIX}} Simulation");
0387   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0388   text.push_back("anti-k_{T}#it{R} = 0.4");
0389   text.push_back("|#eta^{jet}| < 0.7");
0390   text.push_back("#DeltaR_{matching} < 0.3");
0391   text.push_back("With JES calibration");
0392   legend.push_back("Truth jet spectrum");
0393   legend.push_back("Matched truth jet spectrum");
0394   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0395                               false, 10, true,
0396                               true, 15, 72, false,
0397                               false, 0, 0.5, true,
0398                               true, 0.75, 1.1,
0399                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0400                               true, text, 0.55, 0.9, 0.04,
0401                               true, legend, 0.58, 0.6, 0.04,
0402                               "figure/purity_calib_dijet_jerup.pdf");
0403   h_input.clear();
0404   color.clear();
0405   markerstyle.clear();
0406   text.clear();
0407   legend.clear();
0408 
0409   h_input.push_back((TH1F*)h_measure_calib_frac);
0410   h_input.push_back((TH1F*)h_purity_calib_frac_proj);
0411   color.push_back(kBlack);
0412   color.push_back(kRed);
0413   markerstyle.push_back(20);
0414   markerstyle.push_back(20);
0415   text.push_back("#bf{#it{sPHENIX}} Simulation");
0416   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0417   text.push_back("anti-k_{T}#it{R} = 0.4");
0418   text.push_back("|#eta^{jet}| < 0.7");
0419   text.push_back("#DeltaR_{matching} < 0.3");
0420   legend.push_back("Reco jet spectrum");
0421   legend.push_back("Matched reco jet spectrum");
0422   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0423                               false, 10, true,
0424                               true, 15, 72, false,
0425                               false, 0, 0.5, true,
0426                               true, 0.75, 1.1,
0427                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0428                               true, text, 0.55, 0.9, 0.04,
0429                               true, legend, 0.58, 0.65, 0.04,
0430                               "figure/purity_calib_frac.pdf");
0431   h_input.clear();
0432   color.clear();
0433   markerstyle.clear();
0434   text.clear();
0435   legend.clear();
0436 
0437   h_input.push_back((TH1F*)h_measure_calib_frac_jesdown);
0438   h_input.push_back((TH1F*)h_purity_calib_frac_jesdown_proj);
0439   color.push_back(kBlack);
0440   color.push_back(kRed);
0441   markerstyle.push_back(20);
0442   markerstyle.push_back(20);
0443   text.push_back("#bf{#it{sPHENIX}} Simulation");
0444   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0445   text.push_back("anti-k_{T}#it{R} = 0.4");
0446   text.push_back("|#eta^{jet}| < 0.7");
0447   text.push_back("#DeltaR_{matching} < 0.3");
0448   text.push_back("With JES calibration");
0449   legend.push_back("Reco jet spectrum");
0450   legend.push_back("Matched reco jet spectrum");
0451   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0452                               false, 10, true,
0453                               true, 15, 72, false,
0454                               false, 0, 0.5, true,
0455                               true, 0.75, 1.1,
0456                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0457                               true, text, 0.55, 0.9, 0.04,
0458                               true, legend, 0.58, 0.6, 0.04,
0459                               "figure/purity_calib_frac_jesdown.pdf");
0460   h_input.clear();
0461   color.clear();
0462   markerstyle.clear();
0463   text.clear();
0464   legend.clear();
0465 
0466   h_input.push_back((TH1F*)h_measure_calib_frac_jesup);
0467   h_input.push_back((TH1F*)h_purity_calib_frac_jesup_proj);
0468   color.push_back(kBlack);
0469   color.push_back(kRed);
0470   markerstyle.push_back(20);
0471   markerstyle.push_back(20);
0472   text.push_back("#bf{#it{sPHENIX}} Simulation");
0473   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0474   text.push_back("anti-k_{T}#it{R} = 0.4");
0475   text.push_back("|#eta^{jet}| < 0.7");
0476   text.push_back("#DeltaR_{matching} < 0.3");
0477   legend.push_back("Truth jet spectrum");
0478   legend.push_back("Matched truth jet spectrum");
0479   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0480                               false, 10, true,
0481                               true, 15, 72, false,
0482                               false, 0, 0.5, true,
0483                               true, 0.75, 1.1,
0484                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0485                               true, text, 0.55, 0.9, 0.04,
0486                               true, legend, 0.58, 0.65, 0.04,
0487                               "figure/purity_calib_frac_jesup.pdf");
0488   h_input.clear();
0489   color.clear();
0490   markerstyle.clear();
0491   text.clear();
0492   legend.clear();
0493 
0494   h_input.push_back((TH1F*)h_measure_calib_frac_jerdown);
0495   h_input.push_back((TH1F*)h_purity_calib_frac_jerdown_proj);
0496   color.push_back(kBlack);
0497   color.push_back(kRed);
0498   markerstyle.push_back(20);
0499   markerstyle.push_back(20);
0500   text.push_back("#bf{#it{sPHENIX}} Simulation");
0501   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0502   text.push_back("anti-k_{T}#it{R} = 0.4");
0503   text.push_back("|#eta^{jet}| < 0.7");
0504   text.push_back("#DeltaR_{matching} < 0.3");
0505   text.push_back("With JES calibration");
0506   legend.push_back("Truth jet spectrum");
0507   legend.push_back("Matched truth jet spectrum");
0508   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0509                               false, 10, true,
0510                               true, 15, 72, false,
0511                               false, 0, 0.5, true,
0512                               true, 0.75, 1.1,
0513                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0514                               true, text, 0.55, 0.9, 0.04,
0515                               true, legend, 0.58, 0.6, 0.04,
0516                               "figure/purity_calib_frac_jerdown.pdf");
0517   h_input.clear();
0518   color.clear();
0519   markerstyle.clear();
0520   text.clear();
0521   legend.clear();
0522 
0523   h_input.push_back((TH1F*)h_measure_calib_frac_jerup);
0524   h_input.push_back((TH1F*)h_purity_calib_frac_jerup_proj);
0525   color.push_back(kBlack);
0526   color.push_back(kRed);
0527   markerstyle.push_back(20);
0528   markerstyle.push_back(20);
0529   text.push_back("#bf{#it{sPHENIX}} Simulation");
0530   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0531   text.push_back("anti-k_{T}#it{R} = 0.4");
0532   text.push_back("|#eta^{jet}| < 0.7");
0533   text.push_back("#DeltaR_{matching} < 0.3");
0534   text.push_back("With JES calibration");
0535   legend.push_back("Truth jet spectrum");
0536   legend.push_back("Matched truth jet spectrum");
0537   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0538                               false, 10, true,
0539                               true, 15, 72, false,
0540                               false, 0, 0.5, true,
0541                               true, 0.75, 1.1,
0542                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Purity", 1,
0543                               true, text, 0.55, 0.9, 0.04,
0544                               true, legend, 0.58, 0.6, 0.04,
0545                               "figure/purity_calib_frac_jerup.pdf");
0546   h_input.clear();
0547   color.clear();
0548   markerstyle.clear();
0549   text.clear();
0550   legend.clear();
0551 
0552   h_input.push_back((TH1F*)h_truth_calib_dijet);
0553   h_input.push_back((TH1F*)h_efficiency_calib_dijet_proj);
0554   color.push_back(kBlack);
0555   color.push_back(kRed);
0556   markerstyle.push_back(20);
0557   markerstyle.push_back(20);
0558   text.push_back("#bf{#it{sPHENIX}} Simulation");
0559   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0560   text.push_back("anti-k_{T}#it{R} = 0.4");
0561   text.push_back("|#eta^{jet}| < 0.7");
0562   text.push_back("#DeltaR_{matching} < 0.3");
0563   legend.push_back("Reco jet spectrum");
0564   legend.push_back("Matched reco jet spectrum");
0565   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0566                               false, 10, true,
0567                               true, 15, 72, false,
0568                               true, 1e-9, 1e-1, true,
0569                               true, 0., 1.2,
0570                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0571                               true, text, 0.55, 0.9, 0.04,
0572                               true, legend, 0.58, 0.65, 0.04,
0573                               "figure/efficiency_calib_dijet.pdf");
0574   h_input.clear();
0575   color.clear();
0576   markerstyle.clear();
0577   text.clear();
0578   legend.clear();
0579 
0580   h_input.push_back((TH1F*)h_truth_calib_dijet_jesdown);
0581   h_input.push_back((TH1F*)h_efficiency_calib_dijet_jesdown_proj);
0582   color.push_back(kBlack);
0583   color.push_back(kRed);
0584   markerstyle.push_back(20);
0585   markerstyle.push_back(20);
0586   text.push_back("#bf{#it{sPHENIX}} Simulation");
0587   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0588   text.push_back("anti-k_{T}#it{R} = 0.4");
0589   text.push_back("|#eta^{jet}| < 0.7");
0590   text.push_back("#DeltaR_{matching} < 0.3");
0591   text.push_back("With JES calibration");
0592   legend.push_back("Reco jet spectrum");
0593   legend.push_back("Matched reco jet spectrum");
0594   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0595                               false, 10, true,
0596                               true, 15, 72, false,
0597                               true, 1e-9, 1e-1, true,
0598                               true, 0., 1.2,
0599                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0600                               true, text, 0.55, 0.9, 0.04,
0601                               true, legend, 0.58, 0.6, 0.04,
0602                               "figure/efficiency_calib_dijet_jesdown.pdf");
0603   h_input.clear();
0604   color.clear();
0605   markerstyle.clear();
0606   text.clear();
0607   legend.clear();
0608 
0609   h_input.push_back((TH1F*)h_truth_calib_dijet_jesup);
0610   h_input.push_back((TH1F*)h_efficiency_calib_dijet_jesup_proj);
0611   color.push_back(kBlack);
0612   color.push_back(kRed);
0613   markerstyle.push_back(20);
0614   markerstyle.push_back(20);
0615   text.push_back("#bf{#it{sPHENIX}} Simulation");
0616   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0617   text.push_back("anti-k_{T}#it{R} = 0.4");
0618   text.push_back("|#eta^{jet}| < 0.7");
0619   text.push_back("#DeltaR_{matching} < 0.3");
0620   legend.push_back("Truth jet spectrum");
0621   legend.push_back("Matched truth jet spectrum");
0622   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0623                               false, 10, true,
0624                               true, 15, 72, false,
0625                               true, 1e-9, 1e-1, true,
0626                               true, 0, 1.2,
0627                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0628                               true, text, 0.55, 0.9, 0.04,
0629                               true, legend, 0.58, 0.65, 0.04,
0630                               "figure/efficiency_calib_dijet_jesup.pdf");
0631   h_input.clear();
0632   color.clear();
0633   markerstyle.clear();
0634   text.clear();
0635   legend.clear();
0636 
0637   h_input.push_back((TH1F*)h_truth_calib_dijet_jerdown);
0638   h_input.push_back((TH1F*)h_efficiency_calib_dijet_jerdown_proj);
0639   color.push_back(kBlack);
0640   color.push_back(kRed);
0641   markerstyle.push_back(20);
0642   markerstyle.push_back(20);
0643   text.push_back("#bf{#it{sPHENIX}} Simulation");
0644   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0645   text.push_back("anti-k_{T}#it{R} = 0.4");
0646   text.push_back("|#eta^{jet}| < 0.7");
0647   text.push_back("#DeltaR_{matching} < 0.3");
0648   text.push_back("With JES calibration");
0649   legend.push_back("Truth jet spectrum");
0650   legend.push_back("Matched truth jet spectrum");
0651   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0652                               false, 10, true,
0653                               true, 15, 72, false,
0654                               true, 1e-9, 1e-1, true,
0655                               true, 0, 1.2,
0656                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0657                               true, text, 0.55, 0.9, 0.04,
0658                               true, legend, 0.58, 0.6, 0.04,
0659                               "figure/efficiency_calib_dijet_jerdown.pdf");
0660   h_input.clear();
0661   color.clear();
0662   markerstyle.clear();
0663   text.clear();
0664   legend.clear();
0665 
0666   h_input.push_back((TH1F*)h_truth_calib_dijet_jerup);
0667   h_input.push_back((TH1F*)h_efficiency_calib_dijet_jerup_proj);
0668   color.push_back(kBlack);
0669   color.push_back(kRed);
0670   markerstyle.push_back(20);
0671   markerstyle.push_back(20);
0672   text.push_back("#bf{#it{sPHENIX}} Simulation");
0673   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0674   text.push_back("anti-k_{T}#it{R} = 0.4");
0675   text.push_back("|#eta^{jet}| < 0.7");
0676   text.push_back("#DeltaR_{matching} < 0.3");
0677   text.push_back("With JES calibration");
0678   legend.push_back("Truth jet spectrum");
0679   legend.push_back("Matched truth jet spectrum");
0680   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0681                               false, 10, true,
0682                               true, 15, 72, false,
0683                               true, 1e-9, 1e-1, true,
0684                               true, 0, 1.2,
0685                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0686                               true, text, 0.55, 0.9, 0.04,
0687                               true, legend, 0.58, 0.6, 0.04,
0688                               "figure/efficiency_calib_dijet_jerup.pdf");
0689   h_input.clear();
0690   color.clear();
0691   markerstyle.clear();
0692   text.clear();
0693   legend.clear();
0694 
0695   h_input.push_back((TH1F*)h_truth_calib_frac);
0696   h_input.push_back((TH1F*)h_efficiency_calib_frac_proj);
0697   color.push_back(kBlack);
0698   color.push_back(kRed);
0699   markerstyle.push_back(20);
0700   markerstyle.push_back(20);
0701   text.push_back("#bf{#it{sPHENIX}} Simulation");
0702   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0703   text.push_back("anti-k_{T}#it{R} = 0.4");
0704   text.push_back("|#eta^{jet}| < 0.7");
0705   text.push_back("#DeltaR_{matching} < 0.3");
0706   legend.push_back("Reco jet spectrum");
0707   legend.push_back("Matched reco jet spectrum");
0708   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0709                               false, 10, true,
0710                               true, 15, 72, false,
0711                               true, 1e-9, 1e-1, true,
0712                               true, 0., 1.2,
0713                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0714                               true, text, 0.55, 0.9, 0.04,
0715                               true, legend, 0.58, 0.65, 0.04,
0716                               "figure/efficiency_calib_frac.pdf");
0717   h_input.clear();
0718   color.clear();
0719   markerstyle.clear();
0720   text.clear();
0721   legend.clear();
0722 
0723   h_input.push_back((TH1F*)h_truth_calib_frac_jesdown);
0724   h_input.push_back((TH1F*)h_efficiency_calib_frac_jesdown_proj);
0725   color.push_back(kBlack);
0726   color.push_back(kRed);
0727   markerstyle.push_back(20);
0728   markerstyle.push_back(20);
0729   text.push_back("#bf{#it{sPHENIX}} Simulation");
0730   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0731   text.push_back("anti-k_{T}#it{R} = 0.4");
0732   text.push_back("|#eta^{jet}| < 0.7");
0733   text.push_back("#DeltaR_{matching} < 0.3");
0734   text.push_back("With JES calibration");
0735   legend.push_back("Reco jet spectrum");
0736   legend.push_back("Matched reco jet spectrum");
0737   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0738                               false, 10, true,
0739                               true, 15, 72, false,
0740                               true, 1e-9, 1e-1, true,
0741                               true, 0., 1.2,
0742                               true, "p_{T}^{reco jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0743                               true, text, 0.55, 0.9, 0.04,
0744                               true, legend, 0.58, 0.6, 0.04,
0745                               "figure/efficiency_calib_frac_jesdown.pdf");
0746   h_input.clear();
0747   color.clear();
0748   markerstyle.clear();
0749   text.clear();
0750   legend.clear();
0751 
0752   h_input.push_back((TH1F*)h_truth_calib_frac_jesup);
0753   h_input.push_back((TH1F*)h_efficiency_calib_frac_jesup_proj);
0754   color.push_back(kBlack);
0755   color.push_back(kRed);
0756   markerstyle.push_back(20);
0757   markerstyle.push_back(20);
0758   text.push_back("#bf{#it{sPHENIX}} Simulation");
0759   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0760   text.push_back("anti-k_{T}#it{R} = 0.4");
0761   text.push_back("|#eta^{jet}| < 0.7");
0762   text.push_back("#DeltaR_{matching} < 0.3");
0763   legend.push_back("Truth jet spectrum");
0764   legend.push_back("Matched truth jet spectrum");
0765   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0766                               false, 10, true,
0767                               true, 15, 72, false,
0768                               true, 1e-9, 1e-1, true,
0769                               true, 0, 1.2,
0770                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0771                               true, text, 0.55, 0.9, 0.04,
0772                               true, legend, 0.58, 0.65, 0.04,
0773                               "figure/efficiency_calib_frac_jesup.pdf");
0774   h_input.clear();
0775   color.clear();
0776   markerstyle.clear();
0777   text.clear();
0778   legend.clear();
0779 
0780   h_input.push_back((TH1F*)h_truth_calib_frac_jerdown);
0781   h_input.push_back((TH1F*)h_efficiency_calib_frac_jerdown_proj);
0782   color.push_back(kBlack);
0783   color.push_back(kRed);
0784   markerstyle.push_back(20);
0785   markerstyle.push_back(20);
0786   text.push_back("#bf{#it{sPHENIX}} Simulation");
0787   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0788   text.push_back("anti-k_{T}#it{R} = 0.4");
0789   text.push_back("|#eta^{jet}| < 0.7");
0790   text.push_back("#DeltaR_{matching} < 0.3");
0791   text.push_back("With JES calibration");
0792   legend.push_back("Truth jet spectrum");
0793   legend.push_back("Matched truth jet spectrum");
0794   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0795                               false, 10, true,
0796                               true, 15, 72, false,
0797                               true, 1e-9, 1e-1, true,
0798                               true, 0, 1.2,
0799                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0800                               true, text, 0.55, 0.9, 0.04,
0801                               true, legend, 0.58, 0.6, 0.04,
0802                               "figure/efficiency_calib_frac_jerdown.pdf");
0803   h_input.clear();
0804   color.clear();
0805   markerstyle.clear();
0806   text.clear();
0807   legend.clear();
0808 
0809   h_input.push_back((TH1F*)h_truth_calib_frac_jerup);
0810   h_input.push_back((TH1F*)h_efficiency_calib_frac_jerup_proj);
0811   color.push_back(kBlack);
0812   color.push_back(kRed);
0813   markerstyle.push_back(20);
0814   markerstyle.push_back(20);
0815   text.push_back("#bf{#it{sPHENIX}} Simulation");
0816   text.push_back("PYTHIA8 p+p#sqrt{s} = 200 GeV");
0817   text.push_back("anti-k_{T}#it{R} = 0.4");
0818   text.push_back("|#eta^{jet}| < 0.7");
0819   text.push_back("#DeltaR_{matching} < 0.3");
0820   text.push_back("With JES calibration");
0821   legend.push_back("Truth jet spectrum");
0822   legend.push_back("Matched truth jet spectrum");
0823   draw_1D_multiple_plot_ratio(h_input, color, markerstyle,
0824                               false, 10, true,
0825                               true, 15, 72, false,
0826                               true, 1e-9, 1e-1, true,
0827                               true, 0, 1.2,
0828                               true, "p_{T}^{truth jet} [GeV]", "Arbitrary Unit", "Efficiency", 1,
0829                               true, text, 0.55, 0.9, 0.04,
0830                               true, legend, 0.58, 0.6, 0.04,
0831                               "figure/efficiency_calib_frac_jerup.pdf");
0832   h_input.clear();
0833   color.clear();
0834   markerstyle.clear();
0835   text.clear();
0836   legend.clear();
0837 }
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!