|
| gStyle | SetOptStat (0000) |
| |
| gStyle | SetOptFit () |
| |
| gStyle | SetPalette (1) |
| |
| gROOT | SetStyle ("Plain") |
| |
| c1 | Divide (4, 3) |
| |
| | while (1) |
| |
| | fclose (fp) |
| |
| c1 | cd (1) |
| |
| h1 | Draw () |
| |
| h1 | GetXaxis () -> SetLabelSize(0.025) |
| |
| h1 | GetYaxis () -> SetLabelSize(0.025) |
| |
| h1 | SetLineColor (3) |
| |
| h1 | SetFillColor (3) |
| |
| h11 | SetLineColor (8) |
| |
| h11 | SetFillColor (8) |
| |
| h11 | Draw ("same") |
| |
| c1 | cd (5) |
| |
| h2 | SetLineColor (2) |
| |
| h2 | SetFillColor (2) |
| |
| h20 | SetLineColor (5) |
| |
| h20 | SetFillColor (5) |
| |
| c1 | cd (7) |
| |
| ntupleYX | Draw ("Y:X>>hist","vox","colz") |
| |
| gPad | SetLogz () |
| |
| hist | Draw ("colz") |
| |
| hist | GetZaxis () -> SetLabelSize(0.025) |
| |
| hist | SetTitle ("Cytoplasm intensity on transverse section") |
| |
| c1 | cd (3) |
| |
| hist | SetTitle ("Nucleus intensity on transverse section") |
| |
| | system ("rm -rf microbeam.root") |
| |
| | system ("hadd microbeam.root microbeam_*.root") |
| |
| TFile | f ("microbeam.root") |
| |
| c1 | cd (2) |
| |
| ntuple3 | Project ("h1","doseN") |
| |
| h1 | Scale (scale) |
| |
| c1 | cd (6) |
| |
| ntuple3 | Project ("h10","doseC") |
| |
| c1 | cd (9) |
| |
| ntuple0 | Project ("h2","sp") |
| |
| h2 | SetTitle ("dE/dx at cell entrance") |
| |
| h2 | SetFillColor (4) |
| |
| h2 | SetLineColor (4) |
| |
| h2 | Fit ("gaus") |
| |
| gaus | SetLineColor (6) |
| |
| ntuple2 | SetBranchAddress ("x",&x) |
| |
| ntuple2 | SetBranchAddress ("y",&y) |
| |
| ntuple2 | SetBranchAddress ("z",&z) |
| |
| | for (Int_t i=0;i< nentries;i++) |
| |
| c1 | cd (10) |
| |
| ntupleR | Draw ("X2:Z2","abs(X2)<50","surf3") |
| |
| c1 | cd (11) |
| |
| ntuple4 | Draw ("y*0.359060:x*0.359060>>hist","doseV","contz") |
| |
| hist | Draw ("contz") |
| |
| hist | SetTitle ("Mean energy deposit -transverse- (z axis in eV)") |
| |
| c1 | cd (12) |
| |
| ntuple4 | Draw ("x*0.359060:z*0.162810>>hist","doseV","contz") |
| |
| hist | SetTitle ("Mean energy deposit -longitudinal- (z axis in eV)") |
| |
| c1 | cd (4) |
| |
| ntuple1 | Project ("hx","x") |
| |
| h77 | SetTitle ("Beam X position on cell") |
| |
| c1 | cd (8) |
| |
| ntuple1 | Project ("hy","y") |
| |
| h88 | SetTitle ("Beam Y position on cell") |
| |
|
| Double_t | scale = 1/h1->Integral() |
| |
| | c1 = new TCanvas ("c1","",20,20,1200,900) |
| |
| jump | __pad0__ |
| |
| jump | r |
| |
| Float_t | xVox |
| |
| Float_t | yVox |
| |
| Float_t | zVox |
| |
| Float_t | tmp |
| |
| Float_t | den |
| |
| Float_t | dose |
| |
| Float_t | X |
| |
| Float_t | Y |
| |
| Float_t | Z |
| |
| Float_t | vox = 0 |
| |
| Float_t | mat = 0 |
| |
| Float_t | voxelSizeX |
| |
| Float_t | voxelSizeY |
| |
| Float_t | voxelSizeZ |
| |
| TH1F * | h1 = new TH1F("h1","Nucleus marker intensity",100,1,300) |
| |
| TH1F * | h11 = new TH1F("h11 ","",100,1,300) |
| |
| TH1F * | h2 = new TH1F("h2","Cytoplasm marker intensity",100,1,300) |
| |
| TH1F * | h20 = new TH1F("h20 ","",100,1,300) |
| |
| TNtuple * | ntupleYXN = new TNtuple("NUCLEUS","ntuple","Y:X:vox") |
| |
| TNtuple * | ntupleZX = new TNtuple("CYTOPLASM","ntuple","Z:X:vox") |
| |
| TNtuple * | ntupleYX = new TNtuple("CYTOPLASM","ntuple","Y:X:vox") |
| |
| | nlines =0 |
| |
| | ncols =0 |
| |
| TH2F * | hist = new TH2F("hist","hist",50,-20,20,50,-20,20) |
| |
| TNtuple * | ntuple0 = (TNtuple*)f->Get("ntuple0") |
| |
| TNtuple * | ntuple1 = (TNtuple*)f->Get("ntuple1") |
| |
| TNtuple * | ntuple2 = (TNtuple*)f->Get("ntuple2") |
| |
| TNtuple * | ntuple3 = (TNtuple*)f->Get("ntuple3") |
| |
| TNtuple * | ntuple4 = (TNtuple*)f->Get("ntuple4") |
| |
| TH1F * | h10 = new TH1F("h10","Dose distribution in Cytoplasm",100,0.001,.2) |
| |
| Float_t | d |
| |
| Double_t | Xc = -1295.59e3 - 955e3*sin(10*TMath::Pi()/180) |
| |
| Double_t | Zc = -1327e3 + 955e3*cos(10*TMath::Pi()/180) |
| |
| Double_t | X1 |
| |
| Double_t | Y1 |
| |
| Double_t | Z1 |
| |
| Double_t | X2 |
| |
| Double_t | Y2 |
| |
| Double_t | Z2 |
| |
| TNtuple * | ntupleR = new TNtuple("Rmax","ntuple","Z2:Y2:X2") |
| |
| Double_t | x |
| |
| Double_t | y |
| |
| Double_t | z |
| |
| Double_t | xx |
| |
| Double_t | zz |
| |
| Int_t | nentries = (Int_t)ntuple2->GetEntries() |
| |
| TH1F * | h77 = new TH1F("hx","h1",200,-10,10) |
| |
| TH1F * | h88 = new TH1F("hy","h1",200,-10,10) |
| |
