PlotElecCont.C

void PlotElecCont()
{
  float dummy;
  double fECont[2][9][6][12];

  Int_t fMarkerColor[6] = {2,95,209,226,4,221};
  Int_t fMarkerStyle[6][2] = {{24,20},{26,22},{25,21},{27,33},{28,34},{30,29}};
  Int_t fMarkerStyleb[2][2] = {{24,20},{27,33}};

  TCanvas c1("Comp","Comp",3200,1600);

  double z_range[12] = {.225,.275,.325,.375,.425,.475,.525,.575,.625,.675,.725,.8};

  c1.SetFillColor(0);

  c1.Divide(9,5,0,0);

  TGraph* H_acc[2][9][6];

  ifstream econt("electron.txt");
  for(int c=0; c<2; c++)
  {
    for(int i=0; i<9; i++)
    {
      for(int j=0; j<6; j++)
      {
        for(int k=0; k<12; k++)
        {
          econt >> fECont[c][i][j][k] >> dummy;
        }
      }
    }
  }
  econt.close();

  TLine l1(0.1,0.05,0.9,0.05);
  TLine l2(0.1,0.1,0.9,0.1);
  TLine l3(0.1,0.,0.9,0.);

  l1.SetLineStyle(3); l2.SetLineStyle(3);

  for(int c=0; c<2; c++)
  {
    for(int i=0; i<9; i++)
    {
      for(int j=0; j<5; j++)
      {
        std::vector<double> h_a;
        std::vector<double> h_err;
        std::vector<double> z_range_h;

        for(int k=0; k<12; k++)
        {
          fECont[c][i][j][k]=1-fECont[c][i][j][k];
          
          if((j==4 && k>4)
          || (i==0 && j==5 && k==3)
          || (i>4 && j==4 && k==5)
          || (i==0 && j==3 && k==9)
          || (j==3 && k==11)
          || (i==8 && j==4 && k==4)
          || (i==6 && j==3 && k==10)
          || (i==7 && j==3 && k>7)
          || (i==8 && j==3 && k>6)
          || (i>5 && j==5 && k==2)
          || (j==1 && k==3)
          || (i<7 && j==4 && k>5)
          || (i>6 && j==4)
          || (i>5 && j==3 && k>10)
          || (i>7 && j==3 && k>8)
          || (j==2 && k<2)
          || (j==1 && k<4)
          || (j==0 && k<7)
          || (i==8 && j==0)
          || (i==0 && j<3)
          || (i==1 && j==0))
          {
            fECont[c][i][j][k]=0;
          }
          h_a.push_back(fECont[c][i][j][k]);

          z_range_h.push_back(z_range[k]);
        }

        for(int k=12; k>0; k--)
        {
          if(!h_a[k-1]) {h_a.erase(h_a.begin()+k-1); z_range_h.erase(z_range_h.begin()+k-1);}
        }

        bool h_a_empty = 0;

        if(!(int(h_a.size()))) h_a_empty = 1;

        H_acc[c][i][j] = new TGraphErrors(int(h_a.size()),&(z_range_h[0]),&(h_a[0]));

        if(!c)
        {
          H_acc[c][i][j]->SetMarkerColor(fMarkerColor[4]);
        }
        else
        {
          H_acc[c][i][j]->SetMarkerColor(fMarkerColor[0]);
        }

        H_acc[c][i][j]->SetMarkerSize(1);

        H_acc[c][i][j]->SetMarkerStyle(fMarkerStyleb[0][c]);

        H_acc[c][i][j]->GetYaxis()->SetTitle("");

        H_acc[c][i][j]->GetXaxis()->SetTitle("");

        H_acc[c][i][j]->SetTitle("");

        if(!h_a_empty)
        {
          c1.cd(i+1+9*j);
          gPad->SetFillStyle(4000);
          if(H_acc[c][i][j])
          {
            if(!c)
            {
              H_acc[c][i][j]->Draw("SAMEPA");
              H_acc[c][i][j]->GetXaxis()->SetLimits(0.1,0.9);
              H_acc[c][i][j]->SetMinimum(-0.01);
              H_acc[c][i][j]->SetMaximum(0.09);
              H_acc[c][i][j]->GetXaxis()->SetLabelSize(0.06);
              H_acc[c][i][j]->GetYaxis()->SetLabelSize(0.06);
              H_acc[c][i][j]->SetTitle("");
              if(j==5) gPad->SetBottomMargin(.15);
              if(i==0) gPad->SetLeftMargin(.22);
              if(i==8 && j==5)
              {
                H_acc[c][i][j]->GetXaxis()->SetTitle("#font[ 12]{z}");
                H_acc[c][i][j]->GetXaxis()->SetTitleSize(0.08);
                H_acc[c][i][j]->GetXaxis()->SetTitleOffset(.8);
              }
              H_acc[c][i][j]->GetXaxis()->SetNdivisions(304,kTRUE);
              H_acc[c][i][j]->GetYaxis()->SetNdivisions(304,kTRUE);
              if(i==1 && j==0)
              {
                H_acc[c][i][j]->GetYaxis()->SetTitle("#font[12]{e}^{#font[ 12]{#pi}}");
                H_acc[c][i][j]->GetYaxis()->SetTitleSize(0.08);
              }
              H_acc[c][i][j]->Draw("SAMEP");
              l1.Draw("SAME"); l2.Draw("SAME"); l3.Draw("SAME");
              H_acc[c][i][j]->GetXaxis()->SetLimits(0.1,0.9);
              H_acc[c][i][j]->SetMinimum(-0.01);
              H_acc[c][i][j]->SetMaximum(0.09);
              c1.Range(0.1,0.,0.9,1.2);
            }
            else
            {
              H_acc[c][i][j]->Draw("SAMEP");
              l1.Draw("SAME"); l2.Draw("SAME"); l3.Draw("SAME");
              H_acc[c][i][j]->GetXaxis()->SetLimits(0.1,0.9);
              H_acc[c][i][j]->SetMinimum(-0.01);
              H_acc[c][i][j]->SetMaximum(0.09);
            }
          }
          c1.Update();
        }
      }
    }
  }

  c1.Print("Electron_contamination.pdf");
}