{

gROOT->Reset();

TFile f("SDD_histos_test.root");

Int_t nbins = 18;
Float_t dmax = 36.;

a = local->ProjectionX();
t = local->ProjectionY();
local->Draw();

TF1 *faa = new TF1("faa","gaus",-100.,100.);
a->Fit("faa","R","Q",-100.,100.);
TF1 *fat = new TF1("fat","gaus",-50,50);
t->Fit("fat","R","Q",-50,50);

//c1->Size(21.,15.);
c1->SetFillColor(0);
c1->Clear();
c1->Divide(1,3);
c1->cd(1);
local->Draw();
c1->cd(2);
a->Draw();
c1->cd(3);
t->Draw();

c1->SaveAs("ITS_res_1.ps");

TH1F *anode_resolution = new TH1F("anoder","Anode resolution vs Drift Path",nbins,0.,dmax);
TH1D *anodes[nbins];
Float_t res_anodes[nbins];
Float_t errres_anodes[nbins];
Float_t dmin = 60.;

for(Int_t i=0;i<nbins;i++){
  if(i>6) dmin = 70.;
  if(i>11) dmin = 90.;
  TString *aa = new TString("aa_");
  Char_t ai[2];
  sprintf(ai,"%d",i+1);
  aa->Append(ai);
  anodes[i] = at->ProjectionY(aa->Data(),i,i+1);
  TF1 *fa = new TF1("fa","gaus",-1.*dmin,dmin);
  anodes[i]->Fit("fa","R","Q",-1.*dmin,dmin);
  res_anodes[i] = fa->GetParameter(2);
  Float_t RMS = anodes[i]->GetRMS();
  //if(res_anodes[i] > RMS) 
  res_anodes[i] = RMS;
  errres_anodes[i] = fa->GetParError(2);
  anode_resolution->Fill(i*dmax/nbins+dmax/(2*nbins),res_anodes[i]);
  anode_resolution->SetBinError(i+1,(Stat_t) errres_anodes[i]);
  anode_resolution->SetMarkerColor(2);
  anode_resolution->SetLineColor(2);
}
c1->Clear();
f->cd();

TH1F *time_resolution = new TH1F("timer","Time resolution vs Drift Path",nbins,0.,dmax);
TH1D *times[nbins];
Float_t res_times[nbins];
Float_t errres_times[nbins];
for(Int_t i=0;i<nbins;i++){
  TString *ta = new TString("tt_");
  Char_t ti[2];
  sprintf(ti,"%d",i+1);
  ta->Append(ti);
  times[i] = tt->ProjectionY(ta->Data(),i,i+1);
  TF1 *ft = new TF1("ft","gaus",-50,50);
  times[i]->Fit("ft","R","Q",-50,50);
  res_times[i] = ft->GetParameter(2);
  Float_t RMS = times[i]->GetRMS();
  //if(res_times[i] > RMS) 
  res_times[i] = RMS;
  errres_times[i] = ft->GetParError(2);
  time_resolution->Fill(i*dmax/nbins+dmax/(2*nbins),res_times[i]);
  time_resolution->SetBinError(i+1,(Stat_t) errres_times[i]);
  time_resolution->SetMarkerColor(6);
  time_resolution->SetLineColor(6);
}
c1->Clear();

Float_t x1 = 4000.;
Float_t x2 = x1;
Float_t y1 = 94;
Float_t y2 = 84;

TMarker *m1 = new TMarker(x1,y1,21);
TMarker *m2 = new TMarker(x2,y2,22);
Text_t *text1 = "Anode";
Text_t *text2 = "Time";
TText *t1 = new TText(x1 + 250,y1-2,text1);
TText *t2 = new TText(x2 + 250,y2-2,text2);
//text1->SetLineColor(2);
//text2->SetLineColor(4);

anode_resolution->SetMarkerStyle(21);
anode_resolution->SetMaximum(100.);
anode_resolution->SetMinimum(0.);
anode_resolution->Draw("Ep");
//m1->SetColor(2);
//m2->SetColor(4);
m1->Draw();
m2->Draw();
//t1->Draw();
//t2->Draw();
time_resolution->SetMarkerStyle(23);
time_resolution->Draw("Ep,same");
c1->SaveAs("ITS_res.ps");

c1->Clear();
/*
c1->Divide(1,2);
c1->cd(1);
dtrh->Draw();
dtrp->SetLineColor(2);
dtrp->Draw("SAME");

c1->cd(2);
dtrhall->Draw();
dtrpall->SetLineColor(1);
//dtrpall->SetLineStyle(2);
dtrpall->SetLineColor(2);
dtrpall->Draw("SAME");

c1->SaveAs("ITS_dtr_1.ps");
*/
//c1->Size(21.,14.);

c1->Clear();
at->Draw();
c1->SaveAs("ITS_at.ps");

c1->Clear();
tt->Draw();
c1->SaveAs("ITS_tt.ps");

/*
c1->cd(2);
dtrh->Draw();
dtrp->SetLineColor(1);
//dtrp->SetLineStyle(2);
dtrp->SetLineColor(2);
dtrp->Draw("SAME");

c1->SaveAs("ITS_dtr_2.ps");

effh->SetMinimum(0.9);
effh->Draw();
effp->SetLineColor(2);
effp->Draw("SAME");
c1->SaveAs("ITS_eff_1.ps");
*/

hit_vs_time->Draw();
rec_vs_time->SetLineColor(2);
rec_vs_time->Draw("SAME");
c1->SaveAs("ITS_eff_2.ps");

nanodes->ProfileX();
nanodes_pfx->Draw();
nsampls->ProfileX();
nsampls_pfx->Draw();
nanodes_pfx->SetLineColor(2);
nanodes_pfx->Draw("SAME");
c1->SaveAs("ITS_clsize.ps");

/*
occupancy3->Draw();
occupancy4->Draw("SAME");
c1->SaveAs("ITS_occ.ps");

xy3->ProjectionX();
xz3->ProjectionX();
xy3->ProjectionY();

xy3_py->Draw();
xy3_px->Draw("SAME");
xz3_px->Draw("SAME");

c1->SaveAs("ITS_xyz_3.ps");

xy4->ProjectionX();
xz4->ProjectionX();
xy4->ProjectionY();

xy4_py->Draw();
xy4_px->Draw("SAME");
xz4_px->Draw("SAME");

c1->SaveAs("ITS_xyz_4.ps");
*/
amplit->Draw();
c1->SaveAs("ITS_ampl_1.ps");

amplit->ProfileX();
amplit_pfx->Draw();
c1->SaveAs("ITS_ampl_2.ps");

TF1 *fhc = new TF1("fhc","gaus",-5,5);
hitpnt->Fit("fhc","R","Q",-5,5);
hitpnt->Draw();
c1->SaveAs("ITS_hitpnt.ps");

f.cd();
anoder->SetMinimum(0.);
anoder->Draw();
anoder->SetXTitle("Drift Path (mm)");
anoder->SetYTitle("Resolution (um)");

TMarker *mk2 = new TMarker(16.,62.,21);
mk2->SetMarkerColor(2);
mk2->Draw();
TMarker *mk3 = new TMarker(16.,54.,23);
mk3->SetMarkerColor(6);
mk3->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(2,90,text);
Text_t *text0 = "Spatial Resolution";
TText *t0 = new TText(2,82,text0);
Text_t *text2 = "Simulation (Anode)";
TText *t2 = new TText(18,60,text2);
Text_t *text4 = "Simulation (Time)";
TText *t4 = new TText(18,52,text4);

t4->Draw();
t3->Draw();
t0->Draw();
t2->Draw();

timer->Draw("SAME");
c1->SaveAs("ITS_res_check_294.ps");

f.cd();
nanodes_pfx->SetXTitle("Drift Time (ns)");
nanodes_pfx->SetYTitle("Anodes/Cluster");
nanodes_pfx->SetMaximum(3.);
nanodes_pfx->SetMarkerStyle(21);
nanodes_pfx->SetLineColor(2);
nanodes_pfx->SetMarkerColor(2);
nanodes_pfx->Draw();

TMarker *mk2 = new TMarker(4000.,0.3,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(200,2.7,text);
Text_t *text0 = "Average Number of Anodes/Cluster";
TText *t0 = new TText(200,2.4,text0);
Text_t *text2 = "Simulation";
TText *t2 = new TText(4200,0.2,text2);
t3->Draw();
t0->Draw();
t2->Draw();
c1->SaveAs("ITS_and_check_294.ps");

f.cd();
nsampls_pfx->SetXTitle("Drift Time (ns)");
nsampls_pfx->SetYTitle("Time bins/Anode/Cluster");
nsampls_pfx->SetMaximum(10.);
nsampls_pfx->SetMarkerStyle(21);
nsampls_pfx->SetLineColor(2);
nsampls_pfx->SetMarkerColor(2);
nsampls_pfx->Draw();
TMarker *mk2 = new TMarker(4000.,1.,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(200,9.,text);
Text_t *text0 = "Average Number of Time bins/Anode/Cluster";
TText *t0 = new TText(200,8.,text0);
Text_t *text2 = "Simulation";
TText *t2 = new TText(4200,0.8,text2);
t3->Draw();
t0->Draw();
t2->Draw();
c1->SaveAs("ITS_tim_check_294.ps");

f.cd();
amplit->ProfileX();
amplit_pfx->SetLineColor(2);
amplit_pfx->SetMaximum(500.);
amplit_pfx->SetMarkerStyle(21);
amplit_pfx->SetMarkerColor(2);
amplit_pfx->Draw();
TMarker *mk2 = new TMarker(500.,50,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(1000,450.,text);
Text_t *text0 = "Peak Amplitude";
TText *t0 = new TText(1000,400.,text0);
Text_t *text2 = "Simulation";
TText *t2 = new TText(700,40.,text2);
t3->Draw();
t0->Draw();
t2->Draw();
c1->SaveAs("ITS_amp_check_294.ps");

f.cd();
chp->ProfileX();
chp_pfx->SetLineColor(2);
chp_pfx->SetMaximum(2000.);
chp_pfx->SetMarkerStyle(21);
chp_pfx->SetMarkerColor(2);
chp_pfx->Draw();
TMarker *mk2 = new TMarker(500.,200.,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(300,1800.,text);
Text_t *text0 = "Total Charge";
TText *t0 = new TText(300,1600.,text0);
Text_t *text2 = "Simulation";
TText *t2 = new TText(700,150.,text2);
t3->Draw();
t0->Draw();
t2->Draw();
c1->SaveAs("ITS_cha_check_294.ps");

TH1F *eff = f.Get("rec_vs_time");
TH1F *efh = f.Get("hit_vs_time");
eff->Divide(efh);

TH1F *effn = new TH1F("effn","Efficiency vs. drift path (mm)",18,0.,36.);
effn->SetMinimum(0.9);

for(Int_t i=1;i<=36;i++) {
  Float_t cont = eff->GetBinContent(i);
  i++;
  cont += eff->GetBinContent(i);
  cont /= 2.;
  effn->SetBinContent(i/2,cont);
}

effn->SetXTitle("Drift Path (mm)");
effn->SetYTitle("Reconstruction Efficiency");
effn->SetMaximum(1.2);
effn->SetMinimum(0.6);
effn->SetMarkerStyle(21);
effn->SetLineColor(2);
effn->SetMarkerColor(2);
effn->Draw("p");
TMarker *mk2 = new TMarker(20.,0.7,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text = "294 um pitch detector";
TText *t3 = new TText(2.,1.13,text);
Text_t *text2 = "Simulation";
TText *t2 = new TText(22,0.685,text2);
t3->Draw();
t0->Draw();
t2->Draw();
c1->SaveAs("ITS_eff_check_294.ps");

f.cd();
TH1F *pa = (TH1F *) amplit_pfx;
TH1F *pc = (TH1F *) chp_pfx;

TH1F *sd = pa;
sd->Divide(pc);

sd->SetMaximum(0.4);
TH1F *nsd = new TH1F("nsd","nsd",28,0.,7000.);
for(Int_t i=1;i<=28;i++) {
  Float_t sc = sd->GetBinContent(i);
  nsd->SetBinContent(i,sc);
}
nsd->SetMarkerStyle(21);
nsd->SetMarkerColor(2);
nsd->SetLineColor(2);
nsd->SetMaximum(0.4);
nsd->Draw("p");
TMarker *mk2 = new TMarker(300.,0.05,21);
mk2->SetMarkerColor(2);
mk2->Draw();
Text_t *text0 = "294 um pitch detector";
TText *t0 = new TText(300.,0.35,text0);
Text_t *text = "Peak Amplitude / Total Charge";
TText *t3 = new TText(300.,0.3,text);
Text_t *text2 = "Simulation";
TText *t2 = new TText(500,0.04,text2);
t3->Draw();
t0->Draw();
t2->Draw();

c1->SaveAs("ITS_coa_check_294.ps");

}