PWGLF/SPECTRA/PiKaPr/TOF/pPb502/macros/CommonDefs.C

   1
   2 #ifndef COMMONDEFS
   3 #define COMMONDEFS
   4
   5 #define PROTONLEAD 1
   6 #define LEADLEAD 0
   7 #define PROTONPROTON 0
   8
   9 /* analysis dependent stuff */
  10
  11 /***** p-Pb *****/
  12 #if PROTONLEAD
  13
  14 Double_t tofReso = 80.;
  15 Double_t tofTail = 75.;
  16
  17 Double_t scaletexpreso[5] = {1., 1., 1., 1., 1.};
  18 Double_t scaletimezerosigma = 1.;
  19 Double_t forcetimezeroineff = 0.;
  20 Double_t timezero_spread = 200.;
  21
  22 Int_t acceptEventType = 0;
  23 Int_t centralityEstimator = 1;
  24 const Int_t NcentralityBins = 7;
  25 Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 20., 40., 60., 80., 100.};
  26
  27 Float_t rapidityShift = -0.465;
  28 Float_t rapidityMinCut = 0.0;
  29 Float_t rapidityMaxCut = 0.5;
  30
  31 /* TZERO corrections */
  32 Double_t TZEROFILL_shift = 0.;
  33 Double_t TZEROA_shift = 0.;
  34 Double_t TZEROC_shift = 0.;
  35 Double_t TZEROTOF_shift = 0.;
  36
  37 Double_t TZEROvertexCorr = 0.5;
  38
  39 Double_t TZEROA_sigma = 1000.;//247.187;
  40 Double_t TZEROC_sigma = 1000.;//235.113;
  41 Double_t TZEROTOF_sigma = 1000.;//234.423;
  42
  43 Double_t TZEROTOF_resoScaleFactor = 1.25;
  44
  45
  46 Int_t multcentColor[7] = {
  47   kRed,
  48   kPink+1,
  49   kOrange+1,
  50   kYellow+1,
  51   kGreen+1,
  52   kAzure+1,
  53   kViolet+1,
  54 };
  55
  56 /***** Pb-Pb *****/
  57 #elif LEADLEAD
  58
  59 Float_t rapidityShift = 0.;
  60 Float_t rapidityMinCut = -0.5;
  61 Float_t rapidityMaxCut = 0.5;
  62 Float_t rapidityCut = 0.5;
  63
  64 Int_t acceptEventType = 0;
  65 Int_t centralityEstimator = AliAnalysisEvent::kCentEst_V0M;
  66 const Int_t NcentralityBins = 10;
  67 Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 20., 30., 40., 50., 60., 70., 80., 90.};
  68
  69 /* TZERO corrections */
  70 Double_t TZEROFILL_shift = 0.;
  71 Double_t TZEROA_shift = 0.;
  72 Double_t TZEROC_shift = 0.;
  73 Double_t TZEROTOF_shift = 0.;
  74
  75 Double_t TZEROvertexCorr = 0.5;
  76
  77 Double_t TZEROA_sigma = 1000.;//247.187;
  78 Double_t TZEROC_sigma = 1000.;//235.113;
  79 Double_t TZEROTOF_sigma = 1000.;//234.423;
  80
  81 Double_t TZEROTOF_resoScaleFactor = 1.25;
  82
  83
  84 Int_t multcentColor[10] = {
  85   kRed,
  86   kOrange+1,
  87   kOrange,
  88   kYellow,
  89   kYellow+1,
  90   kGreen,
  91   kGreen+1,
  92   kCyan+1,
  93   kBlue,
  94   kMagenta,
  95   //  kMagenta+1
  96 };
  97
  98 /***** p-p *****/
  99 #elif PROTONPROTON /* LHC10d, run 126088 */
 100
 101 Float_t rapidityShift = 0.;
 102 Float_t rapidityMinCut = -0.5;
 103 Float_t rapidityMaxCut = 0.5;
 104 Float_t rapidityCut = 0.5;
 105
 106 Int_t acceptEventType = 1;
 107 Int_t centralityEstimator = 999;
 108 const Int_t NcentralityBins = 6;
 109 Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 15., 20., 25., 100};
 110
 111 /* TZERO corrections */
 112 Double_t TZEROFILL_shift = 0.;//-30.6174;
 113 Double_t TZEROA_shift = -6.5e6;// + -13.8783;
 114 Double_t TZEROC_shift = -6.5e6;// + -21.7581;
 115 Double_t TZEROTOF_shift = 0.;//-45.1779;
 116
 117 Double_t TZEROvertexCorr = 0.5;
 118
 119 Double_t TZEROA_sigma = 1000.;//247.187;
 120 Double_t TZEROC_sigma = 1000.;//235.113;
 121 Double_t TZEROTOF_sigma = 1000.;//234.423;
 122
 123 Double_t TZEROA_resolution[NcentralityBins] = {86.4697,
 124                                                81.6549,
 125                                                79.0132,
 126                                                76.5572,
 127                                                74.8259,
 128                                                73.3834};
 129
 130 Double_t TZEROC_resolution[NcentralityBins] = {58.0423,
 131                                                52.9241,
 132                                                50.8822,
 133                                                49.05,
 134                                                48.6712,
 135                                                46.5776};
 136
 137 Double_t TZEROAND_resolution[NcentralityBins] = {48.383,
 138                                                  51.9983,
 139                                                  51.2704,
 140                                                  50.1974,
 141                                                  51.0198,
 142                                                  48.6588};
 143
 144 Double_t TOFTZEROADIFF[NcentralityBins] = {-20.2311,
 145                                            -19.539,
 146                                            -18.8685,
 147                                            -19.5327,
 148                                            -18.4323,
 149                                            -18.0336};
 150
 151 Double_t TOFTZEROCDIFF[NcentralityBins] = {-16.0008,
 152                                            -11.9249,
 153                                            -9.29299,
 154                                            -8.46063,
 155                                            -8.1128,
 156                                            -4.89226};
 157
 158 Double_t TOFTZEROTOFDIFF[NcentralityBins] = {-11.9418,
 159                                              -12.4526,
 160                                              -14.0921,
 161                                              -18.1132,
 162                                              -17.5909,
 163                                              -18.9787};
 164
 165 Double_t TZEROTOF_resoScaleFactor = 1.25;
 166
 167 Double_t TZEROFILL_sigma = 1000.;
 168
 169 Int_t multcentColor[10] = {
 170   kRed,
 171   kOrange+1,
 172   kOrange,
 173   kYellow,
 174   kYellow+1,
 175   kGreen,
 176   kGreen+1,
 177   kCyan+1,
 178   kBlue,
 179   kMagenta,
 180   //  kMagenta+1
 181 };
 182 #endif
 183
 184 /* commom binning and similar business */
 185
 186 const Char_t *t0FillOnlineFileName = "T0FillOnline.139465.extended.root";
 187 Double_t t0Fill_offset = -1.26416e+04;
 188
 189 //const Char_t *enabledChannelsFileName = "enabledChannels.root";
 190 const Char_t *enabledChannelsFileName = NULL;
 191
 192 const Int_t NptBins = 46;
 193 Double_t ptBin[NptBins + 1] = {0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0};
 194
 195 const Int_t NpBins = 46;
 196 Double_t pBin[NpBins + 1] = {0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0};
 197
 198 const Int_t NmtBins = 46;
 199 Double_t mtBin[NmtBins + 1] = {0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0};
 200
 201 Int_t particleColor[5] = {1, 1, 4, 8, 2};
 202 Int_t chargeMarker[2] = {20, 25};
 203
 204 const Char_t *partChargeName[5][2] = {"e^{+}", "e^{-}", "#mu^{+}", "#mu^{-}", "#pi^{+}", "#pi^{-}", "K^{+}", "K^{-}", "p", "#bar{p}"};
 205
 206 const Double_t kEpsilon = 0.001;
 207
 208 enum ECharge_t {
 209   kPositive,
 210   kNegative,
 211   kNCharges
 212 };
 213 const Char_t *chargeName[kNCharges] = {
 214   "positive",
 215   "negative"
 216 };
 217
 218 const Int_t NetaBins = 10;
 219 Double_t etaMin = -1.;
 220 Double_t etaMax = 1.;
 221 Double_t etaStep = (etaMax - etaMin) / NetaBins;
 222 Double_t etaBin[NetaBins + 1]; /* computed at run-time */
 223
 224 const Int_t NyBins = 20;
 225 Double_t yMin = -1.;
 226 Double_t yMax = 1.;
 227 Double_t yStep = (yMax - yMin) / NyBins;
 228 Double_t yBin[NyBins + 1]; /* computed at run-time */
 229
 230 const Int_t NphiBins = 10;
 231 Double_t phiMin = 0.;
 232 Double_t phiMax = 2. * TMath::Pi();
 233 Double_t phiStep = (phiMax - phiMin) / NphiBins;
 234 Double_t phiBin[NphiBins + 1]; /* computed at run-time */
 235
 236 const Int_t NptsubBins = 4;
 237 Double_t ptsubBin[NptsubBins + 1] = {0.2, 0.5, 1.0, 1.5, 5.0};
 238 Int_t ptsubBinMin[NptsubBins] = {0, 6, 16, 21};
 239 Int_t ptsubBinMax[NptsubBins] = {5, 15, 20, 45};
 240
 241 const Int_t NdcaBins = 2000;
 242 Double_t dcaBin[NdcaBins + 1];
 243 Double_t dcaMin = -5., dcaMax = 5., dcaStep = (dcaMax - dcaMin) / NdcaBins;
 244
 245
 246
 247 #endif