PWG0/dNdPt/macros/plots/readPhojet.C

   1 Int_t readPhojet()
   2 {
   3 // read PHOJET INEL data from file
   4
   5 cout << endl;
   6 cout << "================================" << endl;
   7 cout << endl;
   8 cout << "read PHOJET INEL data from file" <<endl;
   9 cout << "Number of bins: " << binsPhojet <<endl;
  10 cout << "Filename:       " << filenamePhojet <<endl;
  11 cout << endl;
  12 cout << "================================" << endl;
  13 cout << endl;
  14
  15 TFile* filePhojet = new TFile(filenamePhojet,"OPEN");
  16 Double_t nEventsPhojet = ((TH1F*)filePhojet->Get("eventsINEL"))->GetBinContent(1);
  17
  18 TH1F* histPhojet = (TH1F*)filePhojet->Get("ptINEL");
  19
  20 int bins = 0;
  21 for (int i=0; i < histPhojet->GetNbinsX(); i++) {
  22     if (histPhojet->GetBinContent(i) > 0 ) {
  23         centerPtPhojet[bins]    = histPhojet->GetBinCenter(i);
  24         ptPhojet[bins]          = centerPtPhojet[bins];
  25         widthPtPhojet[bins]     = histPhojet->GetBinWidth(i);
  26         errPtPhojet[bins]       = widthPtPhojet[bins] / 2.0;
  27         lowPtPhojet[bins]       = centerPtPhojet[bins] - errPtPhojet[bins];
  28         highPtPhojet[bins]      = centerPtPhojet[bins] + errPtPhojet[bins];
  29
  30         inelPhojet[bins]        = histPhojet->GetBinContent(i) / (nEventsPhojet * etaRange * 2 * M_PI * ptPhojet[bins]);
  31         errInelPhojet[bins]     = histPhojet->GetBinError(i) / (nEventsPhojet * etaRange * 2 * M_PI * ptPhojet[bins]);
  32         lowErrInelPhojet[bins]  = inelPhojet[bins] - errInelPhojet[bins];
  33         highErrInelPhojet[bins] = inelPhojet[bins] + errInelPhojet[bins];
  34         relErrInelPhojet[bins]  = errInelPhojet[bins] / inelPhojet[bins];
  35
  36         centerPt2PiPtPhojet[bins]    = histPhojet->GetBinCenter(i);
  37         pt2PiPtPhojet[bins]          = centerPt2PiPtPhojet[bins];
  38         widthPt2PiPtPhojet[bins]     = histPhojet->GetBinWidth(i);
  39         errPt2PiPtPhojet[bins]       = widthPt2PiPtPhojet[bins] / 2.0;
  40         lowPt2PiPtPhojet[bins]       = centerPt2PiPtPhojet[bins] - errPt2PiPtPhojet[bins];
  41         highPt2PiPtPhojet[bins]      = centerPt2PiPtPhojet[bins] + errPt2PiPtPhojet[bins];
  42         inel2PiPtPhojet[bins]        = histPhojet->GetBinContent(i) / (nEventsPhojet * etaRange);
  43         errInel2PiPtPhojet[bins]     = histPhojet->GetBinError(i) / (nEventsPhojet * etaRange);
  44         lowErrInel2PiPtPhojet[bins]  = inel2PiPtPhojet[bins] - errInel2PiPtPhojet[bins];
  45         highErrInel2PiPtPhojet[bins] = inel2PiPtPhojet[bins] + errInel2PiPtPhojet[bins];
  46         relErrInel2PiPtPhojet[bins]  = errInel2PiPtPhojet[bins] / inel2PiPtPhojet[bins];
  47
  48         cout << "ptPhojet[" << bins << "]              = " << ptPhojet[bins] <<endl;
  49         cout << "   centerPtPhojet[" << bins << "]     = " << centerPtPhojet[bins] <<endl;
  50         cout << "   widthPtPhojet[" << bins << "]      = " << widthPtPhojet[bins] <<endl;
  51         cout << "   errPtPhojet[" << bins << "]        = " << errPtPhojet[bins] <<endl;
  52         cout << "   lowPtPhojet[" << bins << "]        = " << lowPtPhojet[bins] <<endl;
  53         cout << "   highPtPhojet[" << bins << "]       = " << highPtPhojet[bins] <<endl;
  54         cout << "inelPhojet[" << bins << "]            = " << inelPhojet[bins] <<endl;
  55         cout << "errInelPhojet[" << bins << "]         = " << errInelPhojet[bins] <<endl;
  56         cout << "   lowErrInelPhojet[" << bins << "]   = " << lowErrInelPhojet[bins] <<endl;
  57         cout << "   highErrInelPhojet[" << bins << "]  = " << highErrInelPhojet[bins] <<endl;
  58         cout << "   relErrInelPhojet[" << bins << "]   = " << relErrInelPhojet[bins] <<endl;
  59         cout << endl;
  60         bins++;
  61     }
  62 }
  63
  64 filePhojet->Close();
  65
  66 cout << "================================" << endl;
  67 cout << endl;
  68 cout << "Finished reading PHOJET INEL data" <<endl;
  69 cout << "Number of bins read: " << bins <<endl;
  70 cout << endl;
  71 cout << "================================" << endl;
  72 cout << endl;
  73
  74 return bins;
  75 }