LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Continued to adapt to root6

Browse Source
This commit is contained in:
prokscha
2017-03-11 18:45:59 +01:00
parent 66af7291e5
commit b613c0dde1
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
run/macros/test1_tr.C
View File

@ -14,7 +14,7 @@ void sr_transm(const char* fname, int saveID)
TFile* f1 = new TFile(fname);
t1 = (TTree*) gDirectory->FindObjectAny("t1");
sprintf(cutStr,"save_detID== %3i &&save_particleID==-13",saveID);
sprintf(cutStr,"save_detID== %3i &&save_particleID==-13 && save_pz>0",saveID);
gStyle->SetPalette(1,0);
TCanvas* c1 = new TCanvas("c1","Polarization and Transmission",60,40,600,800);
c1->Divide(2,3);
@ -318,9 +318,9 @@ void Beam_Env(const char* fname, Int_t iSlice, Int_t fSlice)
// Resulting beam spot
Cond.Form("save_detID==%d&&save_particleID==-13",i);
// Cond.Form("save_detID==%d&& (save_particleID!=22) ",i);
t1->Draw("save_z:save_y:save_x>>bspot_hist(64,-40,40,64,-40,40,100,-2600,300)",Cond);
t1->Draw("save_z:save_y:save_x>>bspot_hist(64,-80,80,64,-80,80,100,-2600,300)",Cond);
// t1->Draw("save_z:save_y:save_x>>bspot_hist(64,-40,40,64,-40,40,100,-290,290)",Cond);
bspot_hist = (TH2D*) gDirectory->Get("bspot_hist");
// Get properties of slice
zSlice=bspot_hist->GetMean(3);