extract_zdc_histo_from_dst_old_release.C 1.8 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859
  1. /*************************************************************************************
  2. * Example by: Elena Litvinenko
  3. * e-mail: litvin@nf.jinr.ru
  4. * Version: 11-Dec-2012
  5. ************************************************************************************/
  6. void extract_zdc_histo_from_dst (const char *filename="dst_auau.root")
  7. {
  8. Double_t x,y,ee;
  9. TFile *f=new TFile (filename);
  10. TH2F* h_zdc_1_orig=(TH2F*)f->Get("HistZdc1En");
  11. TH2F* h_zdc_2_orig=(TH2F*)f->Get("HistZdc2En");
  12. TChain chain ("cbmsim");
  13. chain.Add(filename);
  14. TChain *t=&chain;
  15. Int_t n_events = t->GetEntries();
  16. TClonesArray *ELossZdc1Histo_arr=0;
  17. t->SetBranchAddress("ELossZdc1Histo",&ELossZdc1Histo_arr);
  18. TClonesArray *ELossZdc2Histo_arr=0;
  19. t->SetBranchAddress("ELossZdc2Histo",&ELossZdc2Histo_arr);
  20. TString main_title_zdc1=h_zdc_1_orig->GetTitle();
  21. TString main_title_zdc2=h_zdc_2_orig->GetTitle();
  22. for (int i=1;i<4;i++) {
  23. t->GetEntry(i);
  24. TH2F *h_zdc_1 = h_zdc_1_orig->Clone("h_zdc_1");
  25. TH2F *h_zdc_2 = h_zdc_2_orig->Clone("h_zdc_2");
  26. TVectorT<float> &v1=*((TVectorT<float>*)ELossZdc1Histo_arr->At(0));
  27. h_zdc_1->Set(v1.GetNoElements(),v1.GetMatrixArray());
  28. h_zdc_1->SetEntries(v1.GetNoElements());
  29. h_zdc_1->SetTitle(Form("%s (event %d)",main_title_zdc1.Data(),i));
  30. cout << h_zdc_1->GetTitle() << " " << h_zdc_1->GetSum() << endl;
  31. if (i<8)
  32. new TCanvas();
  33. h_zdc_1->Draw("colz");
  34. TVectorT<float> &v2=*((TVectorT<float>*)ELossZdc2Histo_arr->First());
  35. h_zdc_2->Set(v2.GetNoElements(),v2.GetMatrixArray());
  36. h_zdc_2->SetEntries(v2.GetNoElements());
  37. h_zdc_2->SetTitle(Form("%s (event %d)",main_title_zdc2.Data(),i));
  38. cout << h_zdc_2->GetTitle() << " " << h_zdc_2->GetSum() << endl;
  39. if (i<8)
  40. new TCanvas();
  41. h_zdc_2->Draw("colz");
  42. }
  43. }