PWGUD/dNdPt/macros/plots/readAliceYield.C

   1 Int_t readAliceYield()
   2 {
   3 // read ALICE INVARIANT YIELD data from file
   4 //
   5
   6 cout << endl;
   7 cout << "================================" << endl;
   8 cout << endl;
   9 cout << "read ALICE INVARIANT YIELD data from file" <<endl;
  10 cout << "Number of bins: " << binsYieldAlice <<endl;
  11 cout << "Filename:       " << filenameYieldAlice <<endl;
  12 cout << endl;
  13 cout << "================================" << endl;
  14 cout << endl;
  15
  16 ifstream fileYieldAlice;
  17 fileYieldAlice.open(filenameYieldAlice);
  18
  19 Int_t i = 0;
  20 while(!fileYieldAlice.eof()) {
  21     if(i == binsYieldAlice) break;
  22     fileYieldAlice >> centerPtYieldAlice[i] >>  yieldAlice[i] >> statYieldAlice[i] >> systYieldAlice[i];
  23     //systYieldAlice[i] = 0.15*yieldAlice[i]+1e-7; // sys error has to be provided, currently 10%
  24
  25     // the width of the pt bins, currently hardwired....
  26     widthPtYieldAlice[i] = 0.05;
  27     if (centerPtYieldAlice[i] > 1) widthPtYieldAlice[i] = 0.1;
  28     if (centerPtYieldAlice[i] > 2) widthPtYieldAlice[i] = 0.2;
  29     if (centerPtYieldAlice[i] > 4) widthPtYieldAlice[i] = 0.5;
  30     if (centerPtYieldAlice[i] > 7) widthPtYieldAlice[i] = 1.0;
  31
  32     ptYieldAlice[i]       = centerPtYieldAlice[i];
  33     errPtYieldAlice[i]    = widthPtYieldAlice[i] / 2.0;
  34     lowPtYieldAlice[i]    = centerPtYieldAlice[i] - errPtYieldAlice[i];
  35     highPtYieldAlice[i]   = centerPtYieldAlice[i] + errPtYieldAlice[i];
  36
  37     lowStatYieldAlice[i]  = yieldAlice[i] - statYieldAlice[i];
  38     highStatYieldAlice[i] = yieldAlice[i] + statYieldAlice[i];
  39     relStatYieldAlice[i]  = statYieldAlice[i] / yieldAlice[i];
  40     lowSystYieldAlice[i]  = yieldAlice[i] - systYieldAlice[i];
  41     highSystYieldAlice[i] = yieldAlice[i] + systYieldAlice[i];
  42     relSystYieldAlice[i]  = systYieldAlice[i] / yieldAlice[i];
  43     errYieldAlice[i]      = systYieldAlice[i] + statYieldAlice[i];
  44     lowErrYieldAlice[i]   = yieldAlice[i] - errYieldAlice[i];
  45     highErrYieldAlice[i]  = yieldAlice[i] + errYieldAlice[i];
  46     relErrYieldAlice[i]   = errYieldAlice[i] / yieldAlice[i];
  47     err2YieldAlice[i]     = sqrt(systYieldAlice[i]*systYieldAlice[i] + statYieldAlice[i]*statYieldAlice[i]);
  48     lowErr2YieldAlice[i]  = yieldAlice[i] - err2YieldAlice[i];
  49     highErr2YieldAlice[i] = yieldAlice[i] + errYieldAlice[i];
  50     relErr2YieldAlice[i]  = err2YieldAlice[i] / yieldAlice[i];
  51
  52
  53
  54
  55     ptYield2PiPtAlice[i]       = ptYieldAlice[i];
  56     centerPtYield2PiPtAlice[i] = centerPtYieldAlice[i];
  57     widthPtYield2PiPtAlice[i]  = widthPtYieldAlice[i];
  58     errPtYield2PiPtAlice[i]    = errPtYieldAlice[i];
  59     lowPtYield2PiPtAlice[i]    = lowPtYieldAlice[i];
  60     highPtYield2PiPtAlice[i]   = highPtYieldAlice[i];
  61
  62     yield2PiPtAlice[i]         = yieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  63     statYield2PiPtAlice[i]     = statYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  64     lowStatYield2PiPtAlice[i]  = lowStatYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  65     highStatYield2PiPtAlice[i] = highStatYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  66     relStatYield2PiPtAlice[i]  = relStatYieldAlice[i];
  67     systYield2PiPtAlice[i]     = systYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  68     lowSystYield2PiPtAlice[i]  = lowSystYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  69     highSystYield2PiPtAlice[i] = highSystYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  70     relSystYield2PiPtAlice[i]  = relSystYieldAlice[i];
  71     errYield2PiPtAlice[i]      = errYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  72     lowErrYield2PiPtAlice[i]   = lowErrYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  73     highErrYield2PiPtAlice[i]  = highErrYieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  74     relErrYield2PiPtAlice[i]   = relErrYieldAlice[i];
  75     err2Yield2PiPtAlice[i]     = err2YieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  76     lowErr2Yield2PiPtAlice[i]  = lowErr2YieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  77     highErr2Yield2PiPtAlice[i] = highErr2YieldAlice[i]*centerPtYieldAlice[i]*2*M_PI;
  78     relErr2Yield2PiPtAlice[i]  = relErr2YieldAlice[i];
  79
  80
  81
  82
  83
  84
  85     cout << "ptYieldAlice[" << i << "]          = " << ptYieldAlice[i] <<endl;
  86     cout << "   centerPtYieldAlice[" << i << "] = " << centerPtYieldAlice[i] <<endl;
  87     cout << "   widthPtYieldAlice[" << i << "]  = " << widthPtYieldAlice[i] <<endl;
  88     cout << "   errPtYieldAlice[" << i << "]    = " << errPtYieldAlice[i] <<endl;
  89     cout << "   lowPtYieldAlice[" << i << "]    = " << lowPtYieldAlice[i] <<endl;
  90     cout << "   highPtYieldAlice[" << i << "]   = " << highPtYieldAlice[i] <<endl;
  91     cout << "yieldAlice[" << i << "]            = " << yieldAlice[i] <<endl;
  92     cout << "   statYieldAlice[" << i << "]     = " << statYieldAlice[i] <<endl;
  93     cout << "   lowStatYieldAlice[" << i << "]  = " << lowStatYieldAlice[i] <<endl;
  94     cout << "   highStatYieldAlice[" << i << "] = " << highStatYieldAlice[i] <<endl;
  95     cout << "   relStatYieldAlice[" << i << "]  = " << relStatYieldAlice[i] <<endl;
  96     cout << "   systYieldAlice[" << i << "]     = " << systYieldAlice[i] <<endl;
  97     cout << "   lowSystYieldAlice[" << i << "]  = " << lowSystYieldAlice[i] <<endl;
  98     cout << "   highSystYieldAlice[" << i << "] = " << highSystYieldAlice[i] <<endl;
  99     cout << "   relSystYieldAlice[" << i << "]  = " << relSystYieldAlice[i] <<endl;
 100     cout << "errYieldAlice[" << i << "]         = " << errYieldAlice[i] <<endl;
 101     cout << "   lowErrYieldAlice[" << i << "]   = " << lowErrYieldAlice[i] <<endl;
 102     cout << "   highErrYieldAlice[" << i << "]  = " << highErrYieldAlice[i] <<endl;
 103     cout << "   relErrYieldAlice[" << i << "]   = " << relErrYieldAlice[i] <<endl;
 104     cout << "err2YieldAlice[" << i << "]        = " << err2YieldAlice[i] <<endl;
 105     cout << "   lowErr2YieldAlice[" << i << "]  = " << lowErr2YieldAlice[i] <<endl;
 106     cout << "   highErr2YieldAlice[" << i << "] = " << highErr2YieldAlice[i] <<endl;
 107     cout << "   relErr2YieldAlice[" << i << "]  = " << relErr2YieldAlice[i] <<endl;
 108     cout << endl;
 109     /*
 110     ptYieldAliceFit[i]      = ptYieldAlice[i];
 111     YieldYieldAliceFit[i]     = YieldYieldAlice[i]*ptYieldAlice[i];
 112     errYieldYieldAliceFit[i]  = errYieldYieldAlice[i];
 113     widthPtYieldAliceFit[i] = widthPtYieldAlice[i];
 114     errPtYieldAliceFit[i]  = 0; //errPtYieldAlice[i];
 115     lowPtYieldAliceFit[i]  = lowPtYieldAlice[i];
 116     highPtYieldAliceFit[i] = highPtYieldAlice[i];
 117     */
 118
 119    i++;
 120 } // while(!fileYieldAlice.eof())
 121 fileYieldAlice.close();
 122 //if (fileYieldAlice) { delete fileYieldAlice; }
 123 //fileYieldAlice=0;
 124
 125 cout << "================================" << endl;
 126 cout << endl;
 127 cout << "Finished reading ALICE INVARIANT YIELD data" <<endl;
 128 cout << "Number of bins read: " << i <<endl;
 129 cout << endl;
 130 cout << "================================" << endl;
 131 cout << endl;
 132
 133 return i;
 134 }