Просмотр исходного кода

Small change in Draw macro: now all TGraph/TGraphErrors are saved in the graphs_e895.root file

PeterParfenov 1 год назад
Родитель
Сommit
aa47566dd2
1 измененных файлов с 135 добавлено и 0 удалено
  1. 135 0
      Draw_e895.C

+ 135 - 0
Draw_e895.C

@@ -84,6 +84,8 @@ void Draw_e895()
   const std::vector<double> v4_pt_bins= {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0};
 
   /////////////////////////////////////////////////////////////////////////////////////////////////
+  TFile *fo = new TFile("./graphs_e895.root", "recreate");
+
   TGraph *gr_pub_v1_y[nene][npid];
   TGraph *gr_pub_v2_pt[nene][npid];
   TGraph *gr_pub_v2_b[nene][npid];
@@ -342,6 +344,24 @@ void Draw_e895()
   {
     for (int ipid=0; ipid<npid; ipid++)
     {
+      gr_pub_v1_y[iene][ipid]->SetName(Form("published_v1_y_ene%i_ipid%i", iene, ipid));
+      gr_pub_v2_pt[iene][ipid]->SetName(Form("published_v2_pt_ene%i_ipid%i", iene, ipid));
+      gr_pub_v2_b[iene][ipid]->SetName(Form("published_v2_b_ene%i_ipid%i", iene, ipid));
+      gr_pub_ptcut_v2_b[iene][ipid]->SetName(Form("published_v2_b_withptcut_ene%i_ipid%i", iene, ipid));
+      gr_pub_ptcut_v2_Mult[iene][ipid]->SetName(Form("published_v1_mult_withptcut_ene%i_ipid%i", iene, ipid));
+
+      gr_pub_v1_y[iene][ipid]->SetTitle(Form("Published v_{1}(y) for ene %i and ipid %i;y;v_{1}", iene, ipid));
+      gr_pub_v2_pt[iene][ipid]->SetTitle(Form("Published v_{2}(p_{T}) for ene %i and ipid %i;p_{T}, GeV/c;v_{2}", iene, ipid));
+      gr_pub_v2_b[iene][ipid]->SetTitle(Form("Published v_{2}(b) for ene %i and ipid %i;b, fm;v_{2}", iene, ipid));
+      gr_pub_ptcut_v2_b[iene][ipid]->SetTitle(Form("Published v_{2}(b) with p_{T}-cut for ene %i and ipid %i;b, fm;v_{2}", iene, ipid));
+      gr_pub_ptcut_v2_Mult[iene][ipid]->SetTitle(Form("Published v_{2}(N_{ch}) with p_{T}-cut for ene %i and ipid %i;N_{ch};v_{2}", iene, ipid));
+    }
+  }
+
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
       gr_pub_v1_y[iene][ipid]->SetMarkerStyle(vMarkerStyle.at(0));
       gr_pub_v1_y[iene][ipid]->SetMarkerColor(vMarkerColor.at(0));
       gr_pub_v1_y[iene][ipid]->SetLineColor(vMarkerColor.at(0));
@@ -1021,6 +1041,121 @@ void Draw_e895()
       canv_v4_Mfilt[iene][ipid]->SaveAs(Form("./pics/c_v4_Mfilt_ien%i_ipid%i.png", iene, ipid));
     }
   }
+
+  fo->cd();
+
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      if (gr_pub_v1_y[iene][ipid]->GetN()>0) gr_pub_v1_y[iene][ipid]->Write();
+    }
+  }
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      if (gr_pub_v2_pt[iene][ipid]->GetN()>0) gr_pub_v2_pt[iene][ipid]->Write();
+      if (gr_pub_v2_b[iene][ipid]->GetN()>0) gr_pub_v2_b[iene][ipid]->Write();
+      if (gr_pub_ptcut_v2_b[iene][ipid]->GetN()>0) gr_pub_ptcut_v2_b[iene][ipid]->Write();
+      if (gr_pub_ptcut_v2_Mult[iene][ipid]->GetN()>0) gr_pub_ptcut_v2_Mult[iene][ipid]->Write();
+    }
+  }
+
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v1_pt[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v1_ycm[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v1_b[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v1_Mcent[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v1_Mfilt[iene][ipid]->Write();
+    }
+  }
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v2_pt[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v2_ycm[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v2_b[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v2_Mcent[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v2_Mfilt[iene][ipid]->Write();
+    }
+  }
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v3_pt[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v3_ycm[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v3_b[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v3_Mcent[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v3_Mfilt[iene][ipid]->Write();
+    }
+  }
+  for (int iene=0; iene<nene; iene++)
+  {
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v4_pt[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v4_ycm[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v4_b[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v4_Mcent[iene][ipid]->Write();
+    }
+    for (int ipid=0; ipid<npid; ipid++)
+    {
+      gr_v4_Mfilt[iene][ipid]->Write();
+    }
+  }
+
+  fo->Close();
 }