macros/scattered-electron_analysis/plot_efficiency_purity.C

0001 int plot_efficiency_purity()
0002 {
0003   vector< TString > v_cuts;
0004
0005   TString base_ep = "em_cluster_e/em_track_ptotal > ";
0006   TString base_prob = "em_cluster_prob >= ";
0007
0008   Double_t cut_ep = 0.90;
0009   Double_t cut_prob = 0.02;
0010   v_cuts.push_back(TString(base_ep) += cut_ep);
0011   v_cuts.push_back(TString(base_prob) += cut_prob);
0012
0013   for ( unsigned i = 0; i < v_cuts.size(); i++ )
0014     {
0015       cout << v_cuts.at(i) << endl;
0016     }
0017
0018   vector< float > v_momenta;
0019   v_momenta.push_back(1.0);
0020   v_momenta.push_back(2.0);
0021   v_momenta.push_back(5.0);
0022   v_momenta.push_back(10.0);
0023   v_momenta.push_back(20.0);
0024
0025   /* Loop over detectors */
0026   {
0027     //...
0028     /* loope over cuts */
0029     {
0030       /* loop over energies */
0031       {
0032     vector< float > v_rates_tp;
0033     vector< float > v_rates_fp;
0034
0035     for ( unsigned idx_p = 0; idx_p < v_momenta.size(); idx_p++ )
0036       {
0037         v_rates_tp.push_back( 1. - 0.1 * v_momenta.at(idx_p) );
0038         v_rates_fp.push_back( 0.1 * v_momenta.at(idx_p) );
0039       }
0040     TGraph* g_tp = new TGraph(v_momenta.size(), &(v_momenta[0]), &(v_rates_tp[0]));
0041     TGraph* g_fp = new TGraph(v_momenta.size(), &(v_momenta[0]), &(v_rates_fp[0]));
0042
0043     g_tp->SetMarkerColor(kGreen+2);
0044     g_fp->SetMarkerColor(kRed+2);
0045
0046     g_tp->SetLineColor(g_tp->GetMarkerColor());
0047     g_fp->SetLineColor(g_fp->GetMarkerColor());
0048
0049     g_tp->Draw("LPA");
0050     g_fp->Draw("LPsame");
0051       }
0052     }
0053   }
0054   return 0;
0055 }