dna/wholeNuclearDNA/plot.C

0001 // *********************************************************************
0002 // To execute this macro under ROOT,
0003 //   1 - launch ROOT (usually type 'root' at your machine's prompt)
0004 //   2 - type '.X plot.C' at the ROOT session prompt
0005 // This macro needs the output ROOT file
0006 // *********************************************************************
0007 {
0008 gROOT->Reset();
0009 gStyle->SetPalette(1);
0010 gROOT->SetStyle("Plain");
0011 Double_t scale;
0012
0013 c1 = new TCanvas ("c1","",20,20,1000,500);
0014
0015 // The command below stands for:
0016 // IF 1 THEN 2
0017 system ("ls wholeNuclearDNA_* 2> /dev/null && hadd -f wholeNuclearDNA.root wholeNuclearDNA_*.root");
0018
0019 TFile f("wholeNuclearDNA.root");
0020
0021 TNtuple* ntuple;
0022 ntuple = (TNtuple*)f.Get("ntuple");
0023
0024 c1->cd(1);
0025 gStyle->SetOptStat(000000);
0026
0027 ntuple->SetMarkerColor(2);
0028 ntuple->SetMarkerStyle(20);
0029 ntuple->SetMarkerSize(1);
0030 ntuple->Draw("x/1000:y/1000:z/1000","flagVolume==1","");
0031
0032 ntuple->SetMarkerColor(4);
0033 ntuple->SetMarkerStyle(24);
0034 ntuple->Draw("x/1000:y/1000:z/1000","flagVolume==2","same");
0035
0036 }