- switch(diaglevel)
- {
- case 1:
-
- TCanvas *c1 = new TCanvas("c1","Alice RICH pad hits",50,10,1200,700);
- c1->Divide(4,2);
- c1->cd(1);
- hc1->SetXTitle("ix (npads)");
- hc1->Draw();
- c1->cd(2);
- hc2->SetXTitle("ix (npads)");
- hc2->Draw();
- c1->cd(3);
- hc3->SetXTitle("ix (npads)");
- hc3->Draw();
- c1->cd(4);
- hc4->SetXTitle("ix (npads)");
- hc4->Draw();
- c1->cd(5);
- hc5->SetXTitle("ix (npads)");
- hc5->Draw();
- c1->cd(6);
- hc6->SetXTitle("ix (npads)");
- hc6->Draw();
- c1->cd(7);
- hc7->SetXTitle("ix (npads)");
- hc7->Draw();
-//
- TCanvas *c4 = new TCanvas("c4","Hits per type",400,10,600,700);
- c4->Divide(2,2);
-
- c4->cd(1);
- feedback->SetFillColor(42);
- feedback->SetXTitle("x (pads)");
- feedback->SetYTitle("y (pads)");
- feedback->Draw();
-
- c4->cd(2);
- mip->SetFillColor(42);
- mip->SetXTitle("x (pads)");
- mip->SetYTitle("y (pads)");
- mip->Draw();
-
- c4->cd(3);
- cerenkov->SetFillColor(42);
- cerenkov->SetXTitle("x (pads)");
- cerenkov->SetYTitle("y (pads)");
- cerenkov->Draw();
-
- c4->cd(4);
- h->SetFillColor(42);
- h->SetXTitle("x (cm)");
- h->SetYTitle("y (cm)");
- h->Draw();
-
- TCanvas *c10 = new TCanvas("c10","Hits distribution",400,10,600,350);
- c10->Divide(2,1);
-
- c10->cd(1);
- hitsX->SetFillColor(42);
- hitsX->SetXTitle("(GeV)");
- hitsX->Draw();
-
- c10->cd(2);
- hitsY->SetFillColor(42);
- hitsY->SetXTitle("(GeV)");
- hitsY->Draw();
-
-
- break;
-//
- case 2:
-
- /*TCanvas *c6 = new TCanvas("c6","Photon Spectra",50,10,600,350);
- c6->Divide(2,1);
-
- c6->cd(1);
- phspectra2->SetFillColor(42);
- phspectra2->SetXTitle("energy (eV)");
- phspectra2->Draw();
- c6->cd(2);
- phspectra1->SetFillColor(42);
- phspectra1->SetXTitle("energy (eV)");
- phspectra1->Draw();*/
-
- //TCanvas *c9 = new TCanvas("c9","Particles Spectra",400,10,600,700);
- TCanvas *c9 = new TCanvas("c9","Pion Spectra",400,10,600,700);
- //c9->Divide(2,2);
-
- //c9->cd(1);
- pionspectra1->SetFillColor(42);
- pionspectra1->SetXTitle("log(GeV)");
- pionspectra2->SetFillColor(46);
- pionspectra2->SetXTitle("log(GeV)");
- pionspectra3->SetFillColor(10);
- pionspectra3->SetXTitle("log(GeV)");
- //c9->SetLogx();
- pionspectra1->Draw();
- pionspectra2->Draw("same");
- pionspectra3->Draw("same");
-
- //c9->cd(2);
-
- TCanvas *c10 = new TCanvas("c10","Proton Spectra",400,10,600,700);
- protonspectra1->SetFillColor(42);
- protonspectra1->SetXTitle("log(GeV)");
- protonspectra2->SetFillColor(46);
- protonspectra2->SetXTitle("log(GeV)");
- protonspectra3->SetFillColor(10);
- protonspectra3->SetXTitle("log(GeV)");
- //c10->SetLogx();
- protonspectra1->Draw();
- protonspectra2->Draw("same");
- protonspectra3->Draw("same");
-
- //c9->cd(3);
- TCanvas *c11 = new TCanvas("c11","Kaon Spectra",400,10,600,700);
- kaonspectra1->SetFillColor(42);
- kaonspectra1->SetXTitle("log(GeV)");
- kaonspectra2->SetFillColor(46);
- kaonspectra2->SetXTitle("log(GeV)");
- kaonspectra3->SetFillColor(10);
- kaonspectra3->SetXTitle("log(GeV)");
- //c11->SetLogx();
- kaonspectra1->Draw();
- kaonspectra2->Draw("same");
- kaonspectra3->Draw("same");
-
- //c9->cd(4);
- TCanvas *c12 = new TCanvas("c12","Charged Particles Spectra",400,10,600,700);
- chargedspectra1->SetFillColor(42);
- chargedspectra1->SetXTitle("log(GeV)");
- chargedspectra2->SetFillColor(46);
- chargedspectra2->SetXTitle("log(GeV)");
- chargedspectra3->SetFillColor(10);
- chargedspectra3->SetXTitle("log(GeV)");
- //c12->SetLogx();
- chargedspectra1->Draw();
- chargedspectra2->Draw("same");
- chargedspectra3->Draw("same");
-
- //TCanvas *c16 = new TCanvas("c16","Particles Spectra II",400,10,600,700);
- //c16->Divide(2,2);
-
- //c16->cd(1);
- TCanvas *c13 = new TCanvas("c13","Electron Spectra",400,10,600,700);
- electronspectra1->SetFillColor(42);
- electronspectra1->SetXTitle("log(GeV)");
- electronspectra2->SetFillColor(46);
- electronspectra2->SetXTitle("log(GeV)");
- electronspectra3->SetFillColor(10);
- electronspectra3->SetXTitle("log(GeV)");
- //c13->SetLogx();
- electronspectra1->Draw();
- electronspectra2->Draw("same");
- electronspectra3->Draw("same");
-
- //c16->cd(2);
- TCanvas *c14 = new TCanvas("c14","Muon Spectra",400,10,600,700);
- muonspectra1->SetFillColor(42);
- muonspectra1->SetXTitle("log(GeV)");
- muonspectra2->SetFillColor(46);
- muonspectra2->SetXTitle("log(GeV)");
- muonspectra3->SetFillColor(10);
- muonspectra3->SetXTitle("log(GeV)");
- //c14->SetLogx();
- muonspectra1->Draw();
- muonspectra2->Draw("same");
- muonspectra3->Draw("same");
-
- //c16->cd(4);
- TCanvas *c16 = new TCanvas("c16","Neutron Spectra",400,10,600,700);
- neutronspectra1->SetFillColor(42);
- neutronspectra1->SetXTitle("log(GeV)");
- neutronspectra2->SetFillColor(46);
- neutronspectra2->SetXTitle("log(GeV)");
- neutronspectra3->SetFillColor(10);
- neutronspectra3->SetXTitle("log(GeV)");
- //c16->SetLogx();
- neutronspectra1->Draw();
- neutronspectra2->Draw("same");
- neutronspectra3->Draw("same");
-
- TCanvas *c15 = new TCanvas("c15","Mothers Production Vertices",500,100,800,800);
- production->SetFillColor(42);
- production->SetXTitle("z (m)");
- production->SetYTitle("R (m)");
- production->Draw();
-
- break;
-
- case 3:
-
- if (nrawclusters) {
- TCanvas *c3=new TCanvas("c3","Clusters Statistics",400,10,600,700);
- c3->Divide(2,2);
-
- c3->cd(1);
- Clcharge->SetFillColor(42);
- Clcharge->SetXTitle("ADC units");
- Clcharge->Draw();
-
- c3->cd(2);
- padnumber->SetFillColor(42);
- padnumber->SetXTitle("(counts)");
- padnumber->Draw();
-
- c3->cd(3);
- clusev->SetFillColor(42);
- clusev->SetXTitle("(counts)");
- clusev->Draw();
- }
-
- if (nev<1)
- {
- TCanvas *c11 = new TCanvas("c11","Cherenkov per Mip",400,10,600,700);
- mother->SetFillColor(42);
- mother->SetXTitle("counts");
- mother->Draw();
- }
-
-
- TCanvas *c5 = new TCanvas("c5","Ring Statistics",50,10,600,350);
- c5->Divide(2,1);
-
- c5->cd(1);
- ckovangle->SetFillColor(42);
- ckovangle->SetXTitle("angle (radians)");
- ckovangle->Draw();
-
- c5->cd(2);
- radius->SetFillColor(42);
- radius->SetXTitle("radius (cm)");
- radius->Draw();
-
- TCanvas *c7 = new TCanvas("c7","Production Statistics",400,10,600,700);
- c7->Divide(2,2);
-
- c7->cd(1);
- totalphotons->SetFillColor(42);
- totalphotons->SetXTitle("Photons (counts)");
- totalphotons->Draw();
-
- c7->cd(2);
- photev->SetFillColor(42);
- photev->SetXTitle("(counts)");
- photev->Draw();
-
- c7->cd(3);
- feedev->SetFillColor(42);
- feedev->SetXTitle("(counts)");
- feedev->Draw();
-
- c7->cd(4);
- padsev->SetFillColor(42);
- padsev->SetXTitle("(counts)");
- padsev->Draw();
-
- break;
-
- }
-
- /*
- TCanvas *c8 = new TCanvas("c25","Number of pads per event inside MIP region",400,10,600,700);
- padsmip->SetFillColor(42);
- padsmip->SetXTitle("(counts)");
- padsmip->Draw();
- */
-
-
- /*TCanvas *c8 = new TCanvas("c8","Number of pads per event inside MIP region",400,10,600,700);
- hckphi->SetFillColor(42);
- hckphi->SetXTitle("phi");
- hckphi->Draw();*/
-