#include "setstyle.C"

void vdc_draw(int runnumber,int choice){
  setstyle();

  if(runnumber<0){
    TChain *T = new TChain("T");
    if(choice==0){
      T->Add("/taidisk/L_det/e07006_L_det_1231.root");
      T->Add("/taidisk/L_det/e07006_L_det_1235.root");
      T->Add("/taidisk/L_det/e07006_L_det_1236.root");
      T->Add("/taidisk/L_det/e07006_L_det_1237.root");
      T->Add("/taidisk/L_det/e07006_L_det_1241.root");
      T->Add("/taidisk/L_det/e07006_L_det_1241_1.root");
      T->Add("/taidisk/L_det/e07006_L_det_1241_2.root");
      T->Add("/taidisk/L_det/e07006_L_det_1242.root");
      T->Add("/taidisk/L_det/e07006_L_det_1243.root");
    }
    if(choice==1){
      T->Add("/taidisk/L_detCheck/e07006_L_det_1231.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1235.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1236.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1237.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1241.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1242.root");
      T->Add("/taidisk/L_detCheck/e07006_L_det_1243.root");
    }
  }
  else{

    if(choice==0) TFile *file = new TFile(Form("/taidisk/L_det/e07006_L_det_%d.root",runnumber));
    if(choice==1) TFile *file = new TFile(Form("/taidisk/L_detCheck/e07006_L_det_%d.root",runnumber));
    TTree *T = (TTree*)file->Get("T");
  }


  if(choice==0) TCanvas *c = new TCanvas("vdc_left_bf",Form("vdc left before %d",runnumber),1200,1000);
  if(choice==1) TCanvas *c = new TCanvas("vdc_left_af",Form("vdc left after %d",runnumber),1200,1000);
  //c->Divide(2,3);

  //c->cd(1);
   T->Draw("L.vdc.u1.time*10^9>>ax(500,-100,400)");
   int ymax = ax->GetMaximum();
   printf("%d\n",ymax);
   TLine *lzero = new TLine(0,0,0,ymax);
   lzero->SetLineColor(3);
   lzero->Draw("same");

   c->Clear();
   c->Divide(2,2);

   c->cd(1); T->Draw("L.vdc.u1.time*10^9>>u1g(500,-100,400)"); lzero->Draw("same"); c->Update();
   printf("lost:total = %2.3e\n",u1g->Integral(0,100)/u1g->Integral(0,500));
   c->cd(2); T->Draw("L.vdc.v1.time*10^9>>v1g(500,-100,400)"); lzero->Draw("same"); c->Update();
   printf("lost:total = %2.3e\n",v1g->Integral(0,100)/v1g->Integral(0,500));
   c->cd(3); T->Draw("L.vdc.u2.time*10^9>>u2g(500,-100,400)"); lzero->Draw("same"); c->Update();
   printf("lost:total = %2.3e\n",u2g->Integral(0,100)/u2g->Integral(0,500));
   c->cd(4); T->Draw("L.vdc.v2.time*10^9>>v2g(500,-100,400)"); lzero->Draw("same"); c->Update();
   printf("lost:total = %2.3e\n",v2g->Integral(0,100)/v2g->Integral(0,500));

}

void vdc_draw(int runnumber){
  vdc_draw(runnumber,0);
  vdc_draw(runnumber,1);
}