[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / PiKaPr / TOF / pPb502 / macros / CommonDefs.C


#ifndef COMMONDEFS
#define COMMONDEFS

#define PROTONLEAD 1
#define LEADLEAD 0
#define PROTONPROTON 0

/* analysis dependent stuff */

/***** p-Pb *****/
#if PROTONLEAD

Double_t tofReso = 80.;
Double_t tofTail = 75.;

Double_t scaletexpreso[5] = {1., 1., 1., 1., 1.};
Double_t scaletimezerosigma = 1.;
Double_t forcetimezeroineff = 0.;
Double_t timezero_spread = 200.;

Int_t acceptEventType = 0;
Int_t centralityEstimator = 1;
const Int_t NcentralityBins = 7;
Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 20., 40., 60., 80., 100.};

Float_t rapidityShift = -0.465;
Float_t rapidityMinCut = 0.0;
Float_t rapidityMaxCut = 0.5;

/* TZERO corrections */
Double_t TZEROFILL_shift = 0.;
Double_t TZEROA_shift = 0.;
Double_t TZEROC_shift = 0.;
Double_t TZEROTOF_shift = 0.;

Double_t TZEROvertexCorr = 0.5;

Double_t TZEROA_sigma = 1000.;//247.187;
Double_t TZEROC_sigma = 1000.;//235.113;
Double_t TZEROTOF_sigma = 1000.;//234.423;

Double_t TZEROTOF_resoScaleFactor = 1.25;


Int_t multcentColor[7] = {
  kRed,
  kPink+1,
  kOrange+1,
  kYellow+1,
  kGreen+1,
  kAzure+1,
  kViolet+1,
};

/***** Pb-Pb *****/
#elif LEADLEAD

Float_t rapidityShift = 0.;
Float_t rapidityMinCut = -0.5;
Float_t rapidityMaxCut = 0.5;
Float_t rapidityCut = 0.5;

Int_t acceptEventType = 0;
Int_t centralityEstimator = AliAnalysisEvent::kCentEst_V0M;
const Int_t NcentralityBins = 10;
Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 20., 30., 40., 50., 60., 70., 80., 90.};

/* TZERO corrections */
Double_t TZEROFILL_shift = 0.;
Double_t TZEROA_shift = 0.;
Double_t TZEROC_shift = 0.;
Double_t TZEROTOF_shift = 0.;

Double_t TZEROvertexCorr = 0.5;

Double_t TZEROA_sigma = 1000.;//247.187;
Double_t TZEROC_sigma = 1000.;//235.113;
Double_t TZEROTOF_sigma = 1000.;//234.423;

Double_t TZEROTOF_resoScaleFactor = 1.25;


Int_t multcentColor[10] = {
  kRed,
  kOrange+1,
  kOrange,
  kYellow,
  kYellow+1,
  kGreen,
  kGreen+1,
  kCyan+1,
  kBlue,
  kMagenta,
  //  kMagenta+1  
};

/***** p-p *****/
#elif PROTONPROTON /* LHC10d, run 126088 */

Float_t rapidityShift = 0.;
Float_t rapidityMinCut = -0.5;
Float_t rapidityMaxCut = 0.5;
Float_t rapidityCut = 0.5;

Int_t acceptEventType = 1;
Int_t centralityEstimator = 999;
const Int_t NcentralityBins = 6;
Double_t centralityBin[NcentralityBins + 1] = {0., 5., 10., 15., 20., 25., 100};

/* TZERO corrections */
Double_t TZEROFILL_shift = 0.;//-30.6174;
Double_t TZEROA_shift = -6.5e6;// + -13.8783;
Double_t TZEROC_shift = -6.5e6;// + -21.7581;
Double_t TZEROTOF_shift = 0.;//-45.1779;

Double_t TZEROvertexCorr = 0.5;

Double_t TZEROA_sigma = 1000.;//247.187;
Double_t TZEROC_sigma = 1000.;//235.113;
Double_t TZEROTOF_sigma = 1000.;//234.423;

Double_t TZEROA_resolution[NcentralityBins] = {86.4697,
					       81.6549,
					       79.0132,
					       76.5572,
					       74.8259,
					       73.3834};

Double_t TZEROC_resolution[NcentralityBins] = {58.0423,
					       52.9241,
					       50.8822,
					       49.05,
					       48.6712,
					       46.5776};

Double_t TZEROAND_resolution[NcentralityBins] = {48.383,
						 51.9983,
						 51.2704,
						 50.1974,
						 51.0198,
						 48.6588};
						  
Double_t TOFTZEROADIFF[NcentralityBins] = {-20.2311,
					   -19.539,
					   -18.8685,
					   -19.5327,
					   -18.4323,
					   -18.0336};

Double_t TOFTZEROCDIFF[NcentralityBins] = {-16.0008,
					   -11.9249,
					   -9.29299,
					   -8.46063,
					   -8.1128,
					   -4.89226};

Double_t TOFTZEROTOFDIFF[NcentralityBins] = {-11.9418,
					     -12.4526,
					     -14.0921,
					     -18.1132,
					     -17.5909,
					     -18.9787};

Double_t TZEROTOF_resoScaleFactor = 1.25;

Double_t TZEROFILL_sigma = 1000.;

Int_t multcentColor[10] = {
  kRed,
  kOrange+1,
  kOrange,
  kYellow,
  kYellow+1,
  kGreen,
  kGreen+1,
  kCyan+1,
  kBlue,
  kMagenta,
  //  kMagenta+1  
};
#endif

/* commom binning and similar business */

const Char_t *t0FillOnlineFileName = "T0FillOnline.139465.extended.root";
Double_t t0Fill_offset = -1.26416e+04;

//const Char_t *enabledChannelsFileName = "enabledChannels.root";
const Char_t *enabledChannelsFileName = NULL;

const Int_t NptBins = 46;
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};

const Int_t NpBins = 46;
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};

const Int_t NmtBins = 46;
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};

Int_t particleColor[5] = {1, 1, 4, 8, 2};
Int_t chargeMarker[2] = {20, 25};

const Char_t *partChargeName[5][2] = {"e^{+}", "e^{-}", "#mu^{+}", "#mu^{-}", "#pi^{+}", "#pi^{-}", "K^{+}", "K^{-}", "p", "#bar{p}"};

const Double_t kEpsilon = 0.001;

enum ECharge_t {
  kPositive,
  kNegative,
  kNCharges
};
const Char_t *chargeName[kNCharges] = {
  "positive",
  "negative"
};

const Int_t NetaBins = 10;
Double_t etaMin = -1.;
Double_t etaMax = 1.;
Double_t etaStep = (etaMax - etaMin) / NetaBins;
Double_t etaBin[NetaBins + 1]; /* computed at run-time */

const Int_t NyBins = 20;
Double_t yMin = -1.;
Double_t yMax = 1.;
Double_t yStep = (yMax - yMin) / NyBins;
Double_t yBin[NyBins + 1]; /* computed at run-time */

const Int_t NphiBins = 10;
Double_t phiMin = 0.;
Double_t phiMax = 2. * TMath::Pi();
Double_t phiStep = (phiMax - phiMin) / NphiBins;
Double_t phiBin[NphiBins + 1]; /* computed at run-time */

const Int_t NptsubBins = 4;
Double_t ptsubBin[NptsubBins + 1] = {0.2, 0.5, 1.0, 1.5, 5.0};
Int_t ptsubBinMin[NptsubBins] = {0, 6, 16, 21};
Int_t ptsubBinMax[NptsubBins] = {5, 15, 20, 45};

const Int_t NdcaBins = 2000;
Double_t dcaBin[NdcaBins + 1];
Double_t dcaMin = -5., dcaMax = 5., dcaStep = (dcaMax - dcaMin) / NdcaBins;


#endif
Commit	Line	Data
59e49925	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