16 gSystem->Exec( Form(
"mkdir %s", (name+
"/"+surface_info.
name).c_str()) );
19 if(surface_info.
type == 2){
21 TText *vol =
new TText(.1,.95,surface_info.
name.c_str());
23 TText *
surface =
new TText(.1,.9,surface_info.
id.c_str());
25 TText *surface_z =
new TText(.1,.85,(
"Z = " +
to_string(surface_info.
pos)).c_str() );
28 TCanvas *c1 =
new TCanvas(
"c1",
"mat_X0",1200,1200);
29 c1->SetRightMargin(0.14);
30 c1->SetTopMargin(0.14);
31 c1->SetLeftMargin(0.14);
32 c1->SetBottomMargin(0.14);
33 Map_prop[0]->Divide(Map_geant[0]);
34 Map_prop[0]->GetZaxis()->SetTitle(
"X0 ratio");
35 Map_prop[0]->SetMaximum(2.);
36 Map_prop[0]->Draw(
"COLZ");
40 c1->Print( (out_name+
"_X0.pdf").c_str());
51 if(surface_info.
type == 1){
53 TText *vol =
new TText(.1,.95,surface_info.
name.c_str());
55 TText *
surface =
new TText(.1,.9,surface_info.
id.c_str());
57 TText *surface_r =
new TText(.1,.85,(
"R = " +
to_string(surface_info.
pos)).c_str() );
60 TCanvas *c1 =
new TCanvas(
"c1",
"mat_X0",1200,1200);
61 c1->SetRightMargin(0.14);
62 c1->SetTopMargin(0.14);
63 c1->SetLeftMargin(0.14);
64 c1->SetBottomMargin(0.14);
65 Map_prop[0]->Divide(Map_geant[0]);
66 Map_prop[0]->GetZaxis()->SetTitle(
"X0 ratio");
67 Map_prop[0]->SetMaximum(2.);
68 Map_prop[0]->Draw(
"COLZ");
72 c1->Print( (out_name+
"_X0.pdf").c_str());
95 gStyle->SetOptStat(0);
96 gStyle->SetOptTitle(0);
98 std::map<uint64_t,std::vector<TH2F*>> surface_hist_prop;
99 std::map<uint64_t,sinfo> surface_info_prop;
101 std::map<uint64_t,std::vector<TH2F*>> surface_hist_geant;
102 std::map<uint64_t,sinfo> surface_info_geant;
104 Fill(surface_hist_prop, surface_info_prop, input_file_prop, nbprocess);
105 Fill(surface_hist_geant, surface_info_geant, input_file_geant, nbprocess);
107 for (
auto hist_it = surface_hist_prop.begin(); hist_it != surface_hist_prop.end(); hist_it++){
108 if(name_prop !=
"")
plot(hist_it->second, surface_info_prop[hist_it->first], name_prop);
109 if(name_geant !=
"")
plot(surface_hist_geant[hist_it->first], surface_info_geant[hist_it->first], name_geant);
110 plot_ratio(hist_it->second,surface_hist_geant[hist_it->first], surface_info_prop[hist_it->first],
name);
112 for (
auto hist : hist_it->second){
115 hist_it->second.clear();
116 for (
auto hist : surface_hist_geant[hist_it->first]){
119 surface_hist_geant[hist_it->first].clear();