PWGUD/dNdPt/macros/plots/settings.C

   1 //------------------------------------------------------------------------------
   2 // settings.C
   3 //
   4 // setting for all data files and folder
   5 // plot options are not set exclusively here
   6 //------------------------------------------------------------------------------
   7
   8
   9 //
  10 // store output plots? (gif+eps) directory?
  11 //
  12 Bool_t SAVE_FIGURES = kTRUE;
  13 // Bool_t SAVE_FIGURES = kFALSE;
  14 TString outputDir = "/u/jacek/CMS_ALICE_UA1_Michael/figures/";
  15
  16 //
  17 // draw ALICE logo and/or PRELIMINARY tag?
  18 //
  19 Bool_t SHOW_LOGO = kTRUE;
  20 Bool_t SHOW_PRELIM = kTRUE;
  21 //Bool_t SHOW_LOGO = kFALSE;
  22 //Bool_t SHOW_PRELIM = kFALSE;
  23
  24 //
  25 // graphics and plot options
  26 //
  27 Int_t textFont = 43;
  28 Int_t titleFont = 43;
  29 Int_t labelFont = 43;
  30 Int_t font = 43;
  31
  32 Float_t titleFontSize = 24;
  33 Float_t markerSize = 1.25;
  34 Float_t legendTextSize = 20;
  35 Float_t labelSize = 24;
  36 Float_t titleSize = 24;
  37
  38 Float_t leftMargin = 0.16;
  39
  40 const Color_t colorAlice = kGray+2;
  41 const Color_t colorAliceFit = kBlack;
  42 const Color_t colorPhojet = kMagenta+2;
  43 const Color_t colorPythia109 = kCyan+3;
  44 const Color_t colorPythia306 = kRed+2;
  45 const Color_t colorPythia320 = kBlue+1;
  46 const Color_t colorAtlas = kBlue+1;
  47 const Color_t colorCms = kRed+2;
  48 const Color_t colorAliceErrors = kGray; //kSpring-8;
  49 const Color_t colorUa1 = kRed+2;
  50
  51
  52 //
  53 // various settings
  54 //
  55
  56 // limit for shifting pt bins
  57 const Double_t fit_limit = 1e-14;
  58 const Double_t fit_limit2 = 1e-14;
  59
  60
  61 // ua1 inelastic cross section (mb); used to calculate yield/
  62 const Double_t sigmaInelUa1 = 43.5;
  63 // yields for (h+ + h-)/2
  64 const Double_t avgToHadr = 2.0;
  65
  66 // eta range used for model comparison
  67 // -0.8 < eta < 0.8
  68 const Double_t etaRange = 1.6;
  69
  70 // pi
  71 const Double_t M_PI = 3.1415926535897932384626433832795028841971693993751;
  72
  73
  74 //
  75 // set filenames and number of data points
  76 //
  77
  78 // compatibility setting
  79 const Int_t binsAlice = 46;
  80
  81 // ALICE INEL data file (for comparison to models)
  82 Char_t* filenameInelAlice = "/u/jacek/alice/dNdPt/output/data_points/ALICE_Yield_INEL_900GeV.txt";
  83 const Int_t binsInelAlice = 46;                 // number of bins (= lines in data textfile)
  84
  85 // ALICE NSD data file (for comparison to ATLAS and CMS)
  86 Char_t* filenameNsdAlice = "/u/jacek/alice/dNdPt/output/data_points/ALICE_Yield_NSD_900GeV.txt";
  87 const Int_t binsNsdAlice = 46;                 // number of bins (= lines in data textfile)
  88
  89 // ALICE NSD Invariant Yield data file (for comparison to UA1)
  90 Char_t* filenameYieldAlice = "/u/jacek/alice/dNdPt/output/data_points/ALICE_InvYield_NSD_900GeV.txt";
  91 const Int_t binsYieldAlice = 46;                 // number of bins (= lines in data textfile)
  92
  93 // ATLAS NSD data file
  94 Char_t* filenameAtlas = "/u/mknichel/dNdPt/data/atlas_data.txt";
  95 const Int_t binsAtlas = 33;                 // number of bins (= lines in data textfile)
  96
  97 // CMS NSD data file
  98 Char_t* filenameCms = "/u/mknichel/dNdPt/data/cms_data.txt";
  99 const Int_t binsCms = 24;                   // number of bins (= lines in data textfile)
 100
 101 // UA1 invariant Cross section data file
 102 Char_t* filenameUa1 = "/u/mknichel/dNdPt/data/UA1data_0.9TeV_pt_CrossSec_errCrossSec.txt";
 103 const Int_t binsUa1 = 52;              // number of bins (= lines in data textfile)
 104
 105 // photjet root file
 106 Char_t* filenamePhojet = "/u/jacek/alice/dNdPt/output/fastsim/out_fastsim_900Phojet_EtaRange0.8.root";
 107 const Int_t binsPhojet = 54;   // number of bins
 108
 109 // Pythia D6T (109) root file
 110 Char_t* filenamePythia109 = "/u/jacek/alice/dNdPt/output/fastsim/out_fastsim_900Pythia109_EtaRange0.8.root";
 111 const Int_t binsPythia109 = 54;   // number of bins
 112
 113 // Pythia ATLAS-CSC (306) root file
 114 Char_t* filenamePythia306 = "/u/jacek/alice/dNdPt/output/fastsim/out_fastsim_900Pythia306_EtaRange0.8.root";
 115 const Int_t binsPythia306 = 54;   // number of bins
 116
 117 // Pythia Perugia0 (320) root file
 118 Char_t* filenamePythia320 = "/u/jacek/alice/dNdPt/output/fastsim/out_fastsim_900Pythia320_EtaRange0.8.root";
 119 const Int_t binsPythia320 = 54;   // number of bins
 120
 121 // ALICE Logo file
 122 Char_t* filenameAliceLogo = "/u/mknichel/dNdPt/LogoALICE-DEF-transp.png";
 123 TImage* logo = TImage::Open(filenameAliceLogo);
 124
 125 //
 126 // fit functions
 127 // these have to be set also in shiftPtxxx.C and makePlotsxxx.C
 128 //
 129 TF1* fitNsd = new TF1("fitNsd","[0]*(x*x/sqrt(0.14*0.14+x*x))*(1+x/[2])^((-1) * [1])",0.0,10.0);
 130 TF1* fitInel = new TF1("fitInel","[0]*(x*x/sqrt(0.14*0.14+x*x))*(1+x/[2])^((-1) * [1])",0.0,10.0);
 131 TF1* fitYield = new TF1("fitYield","[0]*(x*x/sqrt(0.14*0.14+x*x))*(1+x/[2])^((-1) * [1])",0.0,10.0);