]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/GammaConv/AliAnalysisTaskPi0v2.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
including Eta and Pt cuts for AOD analysis (for ESDs this is done in ESDTrackCuts)
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / AliAnalysisTaskPi0v2.cxx
CommitLineData
f7cf81bd 1#include <exception>
2#include "TRandom3.h"
3#include "TChain.h"
4#include "AliLog.h"
5#include "AliAnalysisTask.h"
6#include "AliAnalysisManager.h"
7#include "AliAnalysisTaskPi0v2.h"
8
9#include "AliESDEvent.h"
10#include "AliAODEvent.h"
11#include "AliCentrality.h"
12#include <iostream>
13
14#include "TFile.h"
15#include "AliOADBContainer.h"
16
17// Author Daniel Lohner (Daniel.Lohner@cern.ch)
18
19using namespace std;
20
21ClassImp(AliAnalysisTaskPi0v2)
22
23
24//________________________________________________________________________
25 AliAnalysisTaskPi0v2::AliAnalysisTaskPi0v2(const char *name) : AliAnalysisTaskSE(name),
26 fV0Reader(NULL),
27 fConversionSelection(NULL),
28 fConversionGammas(NULL),
29 fNCentralityBins(1),
30 fCentralityBins(NULL),
31 fCentrality(-1),
32 fCentralityBin(0),
33 fNRadialBins(1),
34 fRadialBins(NULL),
35 fNBinsPhi(6),
36 fEP(NULL),
37 fWeightMultiplicity(kTRUE),
38 fEtaMax(0.75),
39 fEtaGap(1),
40 fRPTPCEtaA(0),
41 fRPTPCEtaC(0),
42 fRPV0A(0),
43 fRPV0C(0),
44 fNCuts(0),
45 fCutList(NULL),
46 fConversionCuts(NULL),
47 fRandomizer(NULL),
48 fOutputList(NULL),
49 fMesonPDGCode(kPi0),
50 fInvMassRange(NULL),
51 fDeltaPsiRP(0),
52 fRunNumber(0),
53 fRunIndex(0),
54 fNEPMethods(knEPMethod),
55 fFillQA(kTRUE),
2bdd97ae 56
57 hNEvents(NULL),
58 hRPTPC(NULL),
59 hRPV0A(NULL),
60 hRPV0C(NULL),
61 hRPTPCAC(NULL),
62 hRPV0ATPC(NULL),
63 hRPV0CTPC(NULL),
64 hRPV0AC(NULL),
65 hCos2TPC(NULL),
66 hCos2V0ATPC(NULL),
67 hCos2V0CTPC(NULL),
68 hCos2V0AC(NULL),
69 hRPTPCEtaA(NULL),
70 hRPTPCEtaC(NULL),
71 hRPTPCEtaAC(NULL),
72 hCos2TPCEta(NULL),
73 hGammaMultCent(NULL),
74 hGammaPhi(NULL),
75 hMultChargedvsNGamma(NULL),
76 hMultChargedvsVZERO(NULL),
77 hMultChargedvsSPD(NULL),
78 hGammadNdPhi(NULL),
79 hGammaMultdPhiTRUE(NULL),
80 hGammaMultdPhiRECOTRUE(NULL),
81 hGammaMultTRUE(NULL),
82 hGammaMultRECOTRUE(NULL),
83 hGammaMultdPhi(NULL),
84 hGammaMult(NULL),
85 hGamma(NULL),
86 hCharged(NULL),
87 hPi0(NULL),
88 hPi0BG(NULL),
89
f7cf81bd 90 fMultV0(0x0),
91 fV0Cpol(0.),
92 fV0Apol(0.)
2bdd97ae 93
f7cf81bd 94{
95 fInvMassRange=new Double_t[2];
96 fInvMassRange[0]=0.05;
97 fInvMassRange[1]=0.3;
98
99 fRandomizer=new TRandom3();
100 fRandomizer->SetSeed(0);
101
102 // Define input and output slots here
103 DefineInput(0, TChain::Class());
104 DefineOutput(1, TList::Class());
105}
106
107//________________________________________________________________________
108AliAnalysisTaskPi0v2::~AliAnalysisTaskPi0v2(){
109
110 if(fRandomizer){
111 delete fRandomizer;
112 fRandomizer=NULL;
113 }
114 if(fCentralityBins){
115 delete fCentralityBins;
116 fCentralityBins=NULL;
117 }
118 if(fInvMassRange){
119 delete fInvMassRange;
120 fInvMassRange=NULL;
121 }
122 if(fRadialBins){
123 delete fRadialBins;
124 fRadialBins=NULL;
125 }
126 if(fCutList){
127 delete fCutList;
128 fCutList=NULL;
129 }
130 if(fConversionSelection){
131 delete fConversionSelection;
132 fConversionSelection=NULL;
133 }
134}
135
136//________________________________________________________________________
137void AliAnalysisTaskPi0v2::UserCreateOutputObjects()
138{
139 OpenFile(1);
140
141 // GetConversionCuts
142 fConversionCuts=fV0Reader->GetConversionCuts();
143
144 // Flags
145
146 Bool_t IsMC=AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler();
147 Bool_t IsHeavyIon=fConversionCuts->IsHeavyIon();
148
149 if(!IsHeavyIon||IsMC)fNEPMethods=1;
150
151 if(!fCentralityBins){
152 fCentralityBins=new Double_t[fNCentralityBins+1];
153 fCentralityBins[0]=-0.5;
154 fCentralityBins[1]=0.5;
155 }
156
157 // Create histograms
158
159 if(fOutputList != NULL){
160 delete fOutputList;
161 fOutputList = NULL;
162 }
163 if(fOutputList == NULL){
164 fOutputList = new TList();
165 fOutputList->SetOwner(kTRUE);
166 }
167
168 Int_t kGCnXBinsSpectra = Int_t((fInvMassRange[1]-fInvMassRange[0])*500); //500 for range 0 - 1
169 Double_t kGCfirstXBinSpectra = fInvMassRange[0];
170 Double_t kGClastXBinSpectra = fInvMassRange[1];
171
172 Int_t nbinspi0[knbinsPi0]={kGCnYBinsSpectra,kGCnXBinsSpectra,6,fNCentralityBins,fNEPMethods};
173 Double_t minpi0[knbinsPi0]={kGCfirstYBinSpectra,kGCfirstXBinSpectra,0,-0.5,-0.5};
174 Double_t maxpi0[knbinsPi0]={kGClastYBinSpectra,kGClastXBinSpectra,TMath::Pi()/2,fNCentralityBins-0.5,fNEPMethods-0.5};
175 const char *binpi0[knbinsPi0]={"pt","mass","dPhi","centr","EPm"};
176
177 Int_t nbinsg[knbinsGamma]={kGCnYBinsSpectra,fNRadialBins,6,fNCentralityBins,fNEPMethods};
178 Double_t ming[knbinsGamma]={kGCfirstYBinSpectra,-0.5,0,-0.5,-0.5};
179 Double_t maxg[knbinsGamma]={kGClastYBinSpectra,fNRadialBins-0.5,TMath::Pi()/2,fNCentralityBins-0.5,fNEPMethods-0.5};
180 const char *bingamma[knbinsGamma]={"pt","R","dPhi","centr","EPm"};
181
182 // Define Binning
183
184 if(!IsMC){
185
186 hPi0=new THnSparseF*[fNCuts];
187 hPi0BG=new THnSparseF*[fNCuts];
188 hGamma=new THnSparseF*[fNCuts];
189
190 // Photon Cuts
191 for(Int_t iCut=0;iCut<fNCuts;iCut++){
192 AliConversionCuts *fPhotonv2Cuts=(AliConversionCuts*)fCutList->At(iCut);
193
194 TList *fCutOutputList=new TList();
195 fCutOutputList->SetName(fPhotonv2Cuts->GetCutNumber().Data());
196 fCutOutputList->SetOwner(kTRUE);
197 fOutputList->Add(fCutOutputList);
198
199 /*fPhotonv2Cuts->SetFillCutHistograms("",kFALSE);
200 TList *fCutQAList=fPhotonv2Cuts->GetCutHistograms();
201 if(fCutQAList){
202 fCutQAList->SetOwner(kTRUE);
203 fCutOutputList->Add(fCutQAList);
204 } */
205
206 hPi0[iCut]=new THnSparseF("Pi0_Sparse","Pi0_Sparse",knbinsPi0,nbinspi0,minpi0,maxpi0);
207 for(Int_t i=0;i<knbinsPi0;i++) hPi0[iCut]->GetAxis(i)->SetName(binpi0[i]);
208 hPi0[iCut]->Sumw2();
209 fCutOutputList->Add(hPi0[iCut]);
210
211 hPi0BG[iCut]=new THnSparseF("Pi0BG_Sparse","Pi0BG_Sparse",knbinsPi0,nbinspi0,minpi0,maxpi0);
212 for(Int_t i=0;i<knbinsPi0;i++) hPi0BG[iCut]->GetAxis(i)->SetName(binpi0[i]);
213 hPi0BG[iCut]->Sumw2();
214 fCutOutputList->Add(hPi0BG[iCut]);
215
216 // Gamma
217
218 hGamma[iCut]=new THnSparseF("Gamma_Sparse","Gamma_Sparse",knbinsGamma,nbinsg,ming,maxg);
219 for(Int_t i=0;i<knbinsGamma;i++) hGamma[iCut]->GetAxis(i)->SetName(bingamma[i]);
220 hGamma[iCut]->Sumw2();
221 fCutOutputList->Add(hGamma[iCut]);
222
223 }
224 }
225
226 if(IsHeavyIon&&!IsMC){
227
228 // RP Calculation
229 TList *fRPList=new TList();
230 fRPList->SetName("Event Plane");
231 fRPList->SetOwner(kTRUE);
232 fOutputList->Add(fRPList);
233
234 hRPTPC=new TH2F("TPCAC" ,"TPC_AC" , fNCentralityBins,fCentralityBins, 180, 0, TMath::Pi());
235 hRPTPC->Sumw2();
236 fRPList->Add(hRPTPC);
237 hRPTPCEtaA=new TH2F("TPCEtaA" ,"TPC_EtaA" , fNCentralityBins,fCentralityBins, 180, 0, TMath::Pi());
238 hRPTPCEtaA->Sumw2();
239 fRPList->Add(hRPTPCEtaA);
240 hRPTPCEtaC=new TH2F("TPCEtaC" ,"TPC_EtaC" , fNCentralityBins,fCentralityBins, 180, 0, TMath::Pi());
241 hRPTPCEtaC->Sumw2();
242 fRPList->Add(hRPTPCEtaC);
243 hRPV0A=new TH2F("V0A" ,"VZERO_A" , fNCentralityBins,fCentralityBins, 180, 0, TMath::Pi());
244 hRPV0A->Sumw2();
245 fRPList->Add(hRPV0A);
246 hRPV0C=new TH2F("V0C" ,"VZERO_C" , fNCentralityBins,fCentralityBins, 180, 0, TMath::Pi());
247 hRPV0C->Sumw2();
248 fRPList->Add(hRPV0C);
249
250 hRPTPCAC=new TH2F("TPCA_TPCC" ,"TPCA_TPCC" , 180, 0, TMath::Pi(), 180, 0, TMath::Pi());
251 hRPTPCAC->Sumw2();
252 fRPList->Add(hRPTPCAC);
253 hRPV0ATPC=new TH2F("V0A_TPC" ,"V0A_TPC" , 180, 0, TMath::Pi(), 180, 0, TMath::Pi());
254 hRPV0ATPC->Sumw2();
255 fRPList->Add(hRPV0ATPC);
256 hRPV0CTPC=new TH2F("V0C_TPC" ,"V0C_TPC" , 180, 0, TMath::Pi(), 180, 0, TMath::Pi());
257 hRPV0CTPC->Sumw2();
258 fRPList->Add(hRPV0CTPC);
259 hRPV0AC=new TH2F("V0A_V0C" ,"V0A_V0C" , 180, 0, TMath::Pi(), 180, 0, TMath::Pi());
260 hRPV0AC->Sumw2();
261 fRPList->Add(hRPV0AC);
262 hRPTPCEtaAC=new TH2F("TPCEtaA_TPCEtaC" ,"TPCEtaA_TPCEtaC" , 180, 0, TMath::Pi(), 180, 0, TMath::Pi());
263 hRPTPCEtaAC->Sumw2();
264 fRPList->Add(hRPTPCEtaAC);
265
266 hCos2TPC=new TH2F("Cos2_TPCAC" ,"Cos2_TPCAC" ,fNCentralityBins,fCentralityBins,100,-1,1);
267 hCos2TPC->Sumw2();
268 fRPList->Add(hCos2TPC);
269 hCos2TPCEta=new TH2F("Cos2_TPCEtaAC" ,"Cos2_TPCEtaAC" ,fNCentralityBins,fCentralityBins,100,-1,1);
270 hCos2TPCEta->Sumw2();
271 fRPList->Add(hCos2TPCEta);
272 hCos2V0ATPC=new TH2F("Cos2_V0ATPC" ,"Cos2_V0ATPC" ,fNCentralityBins,fCentralityBins,100,-1,1);
273 hCos2V0ATPC->Sumw2();
274 fRPList->Add(hCos2V0ATPC);
275 hCos2V0CTPC=new TH2F("Cos2_V0CTPC" ,"Cos2_V0CTPC" ,fNCentralityBins,fCentralityBins,100,-1,1);
276 hCos2V0CTPC->Sumw2();
277 fRPList->Add(hCos2V0CTPC);
278 hCos2V0AC=new TH2F("Cos2_V0AC" ,"Cos2_V0AC" ,fNCentralityBins,fCentralityBins,100,-1,1);
279 hCos2V0AC->Sumw2();
280 fRPList->Add(hCos2V0AC);
281 }
282
283 TList *fPhotonQAList=new TList();
284 fPhotonQAList->SetName("Gamma_QA");
285 fPhotonQAList->SetOwner(kTRUE);
286 fOutputList->Add(fPhotonQAList);
287
288 if(fFillQA){
289 // Gamma QA
290 hGammaPhi=new TH3F*[fNCentralityBins];
291 for(Int_t m=0;m<fNCentralityBins;m++){
292 hGammaPhi[m]=new TH3F(Form("%d_GammaPhi",m),"GammaPhi",kGCnYBinsSpectra,kGCfirstYBinSpectra,kGClastYBinSpectra,fNRadialBins,-0.5,fNRadialBins-0.5,360,0,2*TMath::Pi());
293 hGammaPhi[m]->Sumw2();
294 fPhotonQAList->Add(hGammaPhi[m]);
295 }
296
297 hGammaMultCent=new TH2F("GammaMultvsCent","GammaMultvsCent",fNCentralityBins,fCentralityBins, 60,-0.5,59.5);
298 hGammaMultCent->Sumw2();
299 fPhotonQAList->Add(hGammaMultCent);
300
301 hMultChargedvsSPD=new TH2F("Mult_ChargedvsSPD","Mult_ChargedvsSPD",250,0,2500, 250,0,5000);
302 hMultChargedvsSPD->Sumw2();
303 fPhotonQAList->Add(hMultChargedvsSPD);
304 hMultChargedvsVZERO=new TH2F("Mult_ChargedvsVZERO","Mult_ChargedvsVZERO",250,0,2500, 200,0,20000);
305 hMultChargedvsVZERO->Sumw2();
306 fPhotonQAList->Add(hMultChargedvsVZERO);
307 hMultChargedvsNGamma=new TH2F("Mult_ChargedvsNGamma","Mult_ChargedvsNGamma",250,0,2500,60,-0.5,59.5);
308 hMultChargedvsNGamma->Sumw2();
309 fPhotonQAList->Add(hMultChargedvsNGamma);
310
311 Int_t nbinsgmult[knbinsGammaMult]={kGCnYBinsSpectra,fNRadialBins,400,fNCentralityBins};
312 Double_t mingmult[knbinsGammaMult]={kGCfirstYBinSpectra,-0.5,0,-0.5};
313 Double_t maxgmult[knbinsGammaMult]={kGClastYBinSpectra,fNRadialBins-0.5,8000,fNCentralityBins-0.5};
314 Double_t maxgmultdPhi[knbinsGammaMult]={kGClastYBinSpectra,fNRadialBins-0.5,2000,fNCentralityBins-0.5};
315 const char *bingammamult[knbinsGammaMult]={"pt","R","gmult","centr"};
316
317 if(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()){
318
319 hGammaMultdPhiTRUE=new THnSparseF("Gamma_MultdPhi_TRUE","Gamma_MultdPhi_TRUE",knbinsGammaMult,nbinsgmult,mingmult,maxgmultdPhi);
320 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMultdPhiTRUE->GetAxis(i)->SetName(bingammamult[i]);
321 hGammaMultdPhiTRUE->Sumw2();
322 fPhotonQAList->Add(hGammaMultdPhiTRUE);
323
324 hGammaMultTRUE=new THnSparseF("Gamma_Mult_TRUE","Gamma_Mult_TRUE",knbinsGammaMult,nbinsgmult,mingmult,maxgmult);
325 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMultTRUE->GetAxis(i)->SetName(bingammamult[i]);
326 hGammaMultTRUE->Sumw2();
327 fPhotonQAList->Add(hGammaMultTRUE);
328
329 hGammaMultdPhiRECOTRUE=new THnSparseF("Gamma_MultdPhi_RECOTRUE","Gamma_MultdPhi_RECOTRUE",knbinsGammaMult,nbinsgmult,mingmult,maxgmultdPhi);
330 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMultdPhiRECOTRUE->GetAxis(i)->SetName(bingammamult[i]);
331 hGammaMultdPhiRECOTRUE->Sumw2();
332 fPhotonQAList->Add(hGammaMultdPhiRECOTRUE);
333
334 hGammaMultRECOTRUE=new THnSparseF("Gamma_Mult_RECOTRUE","Gamma_Mult_RECOTRUE",knbinsGammaMult,nbinsgmult,mingmult,maxgmult);
335 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMultRECOTRUE->GetAxis(i)->SetName(bingammamult[i]);
336 hGammaMultRECOTRUE->Sumw2();
337 fPhotonQAList->Add(hGammaMultRECOTRUE);
338 }
339
340 hGammaMultdPhi=new THnSparseF*[fNEPMethods];
341 hGammaMult=new THnSparseF*[fNEPMethods];
342
343 hGammadNdPhi=new THnSparseF("Gamma_dNdPhi","Gamma_dNdPhi",knbinsGamma,nbinsg,ming,maxg);
344 for(Int_t i=0;i<knbinsGamma;i++) hGammadNdPhi->GetAxis(i)->SetName(bingamma[i]);
345 hGammadNdPhi->Sumw2();
346 fPhotonQAList->Add(hGammadNdPhi);
347
348 for(Int_t iEP=0;iEP<fNEPMethods;iEP++){
349 hGammaMultdPhi[iEP]=new THnSparseF(Form("Gamma_MultdPhi_%d",iEP),"Gamma_MultdPhi",knbinsGammaMult,nbinsgmult,mingmult,maxgmultdPhi);
350 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMultdPhi[iEP]->GetAxis(i)->SetName(bingammamult[i]);
351 hGammaMultdPhi[iEP]->Sumw2();
352 fPhotonQAList->Add(hGammaMultdPhi[iEP]);
353
354 hGammaMult[iEP]=new THnSparseF(Form("Gamma_Mult_%d",iEP),"Gamma_Mult",knbinsGammaMult,nbinsgmult,mingmult,maxgmult);
355 for(Int_t i=0;i<knbinsGammaMult;i++) hGammaMult[iEP]->GetAxis(i)->SetName(bingammamult[i]);
356 hGammaMult[iEP]->Sumw2();
357 fPhotonQAList->Add(hGammaMult[iEP]);
358 }
359
360 const Int_t knbinsCharged=3;
361 Int_t nbinscharged[knbinsCharged]={6,fNCentralityBins,fNEPMethods};
362 Double_t mincharged[knbinsCharged]={0,-0.5,-0.5};
363 Double_t maxcharged[knbinsCharged]={TMath::Pi()/2,fNCentralityBins-0.5,fNEPMethods-0.5};
364 hCharged=new THnSparseF("Charged","Charged",knbinsCharged,nbinscharged,mincharged,maxcharged);
365 hCharged->GetAxis(0)->SetName("dPhi");
366 hCharged->GetAxis(1)->SetName("centr");
367 hCharged->GetAxis(2)->SetName("EPm");
368 hCharged->Sumw2();
369 fPhotonQAList->Add(hCharged);
370 }
371 hNEvents=new TH1F("NEvents","NEvents",fNCentralityBins,fCentralityBins);
372 fPhotonQAList->Add(hNEvents);
373
374
375 // V0 Reader Cuts
376 TList *fV0ReaderCuts=fConversionCuts->GetCutHistograms();
377 fV0ReaderCuts->SetOwner(kTRUE);
378 fOutputList->Add(fV0ReaderCuts);
379
380 PostData(1, fOutputList);
381
382}
383
384//________________________________________________________________________
385void AliAnalysisTaskPi0v2::InitConversionSelection(){
386
387 fConversionSelection=new AliConversionSelection*[fNCuts];
388 for(Int_t iCut=0;iCut<fNCuts;iCut++){
389 AliConversionCuts *fPhotonv2Cuts=(AliConversionCuts*)fCutList->At(iCut);
390 if(iCut==0)fEtaMax=fPhotonv2Cuts->GetEtaCut();
391 fConversionSelection[iCut]=new AliConversionSelection(fPhotonv2Cuts);
392 fConversionSelection[iCut]->SetInvMassRange(fInvMassRange);
393 }
394}
395
396//________________________________________________________________________
397Bool_t AliAnalysisTaskPi0v2::InitEvent(){
398
399 if(!fV0Reader){AliError("Error: No V0 Reader and Pi0 Reconstructor");return kFALSE;}
400 if(!fV0Reader->IsEventSelected())return kFALSE;
401 fConversionGammas=fV0Reader->GetReconstructedGammas();
402
403 if(!fConversionSelection){
404 InitConversionSelection();
405 }
406
407 if(!SetCentrality()){return kFALSE;}
408
409 if(fConversionCuts->IsHeavyIon()){
410
411 if(fRunNumber!=fInputEvent->GetRunNumber()){
412 fRunNumber=fInputEvent->GetRunNumber();
413 fRunIndex=GetRunIndex(fRunNumber);
414 OpenInfoCalibration(fRunNumber); // Load Calibration of V0 Event Plane
415 }
416 if(fRunIndex<0)return kFALSE;
417
418 fEP=fInputEvent->GetEventplane();
419 if(!fEP)return kFALSE;
420
421 fRPTPCEtaA=GetTPCSubEPEta(fEtaGap/2,1);
422 fRPTPCEtaC=GetTPCSubEPEta(-1,-fEtaGap/2);
423
424 // GetV0 Event Plane
425 GetV0EP(fInputEvent);
426 }
427 return kTRUE;
428}
429
430//________________________________________________________________________
431void AliAnalysisTaskPi0v2::UserExec(Option_t *)
432{
433
434 if(!InitEvent())return;
435
436 // Process Cuts
437 for(Int_t iCut=0;iCut<fNCuts;iCut++){
438
439 if(fConversionSelection[iCut]->ProcessEvent(fConversionGammas,fInputEvent,fMCEvent)){ // May only be called once per event!
440
441 if(!fMCEvent){
442
443 // Process EP methods
444 for(Int_t iEP=0;iEP<fNEPMethods;iEP++){
445
446 ProcessPi0s(iCut,EEventPlaneMethod(iEP));
447
448 ProcessGammas(iCut,EEventPlaneMethod(iEP));
449 }
450 }
451 }
452 }
453
454 // Fill N Events
455 hNEvents->Fill(fCentrality);
456
457 // EventPlaneResolution
458 ProcessEventPlane();
459
460 // QA
461 if(fFillQA)ProcessQA();
462
463 PostData(1, fOutputList);
464
465}
466
467//________________________________________________________________________
468void AliAnalysisTaskPi0v2::ProcessPi0s(Int_t iCut,EEventPlaneMethod iEP){
469
470 if(!fConversionSelection[iCut])return;
471
472 Int_t photonMultiplicity=fConversionSelection[iCut]->GetMultiplicity(fInputEvent);
473 if(photonMultiplicity==0)return;
474
475 // Process Pi0s
476 Float_t weight=1; // Weight for Multiplicity
477 if(fWeightMultiplicity){
478 weight*=1/Float_t(photonMultiplicity);
479 }
480
481 for(Int_t ii=0;ii<fConversionSelection[iCut]->GetNumberOfPi0s();ii++){
482
483 AliAODConversionMother *pi0cand=fConversionSelection[iCut]->GetPi0(ii);
484
485 if(!pi0cand)continue;
486
487 Double_t val[knbinsPi0];
488 val[0]=pi0cand->Pt();
489 val[1]=pi0cand->M();
490 val[2]=GetPi0PhiwrtRP(pi0cand,iEP);
491 val[3]=fCentralityBin;
492 val[4]=Int_t(iEP);
493
494 hPi0[iCut]->Fill(val,weight);
495
496 }
497
498 // Pi0 BG
499 for(Int_t ii=0;ii<fConversionSelection[iCut]->GetNumberOfBGs();ii++){
500
501 AliAODConversionMother *pi0cand=fConversionSelection[iCut]->GetBG(ii);
502
503 if(!pi0cand)continue;
504
505 Double_t val[knbinsPi0];
506 val[0]=pi0cand->Pt();
507 val[1]=pi0cand->M();
508 val[2]=GetPi0PhiwrtRP(pi0cand,iEP);
509 val[3]=fCentralityBin;
510 val[4]=Int_t(iEP);
511
512 hPi0BG[iCut]->Fill(val,pi0cand->GetWeight()*weight);
513 }
514}
515
516//________________________________________________________________________
517void AliAnalysisTaskPi0v2::ProcessGammas(Int_t iCut,EEventPlaneMethod iEP){
518
519 if(!fConversionSelection[iCut])return;
520
521 Int_t photonMultiplicity=fConversionSelection[iCut]->GetMultiplicity(fInputEvent);
522 if(photonMultiplicity==0)return;
523
524 Float_t weight=1; // Weight for Multiplicity
525 if(fWeightMultiplicity){
526 weight*=1/Float_t(photonMultiplicity);
527 }
528
529 for(Int_t ii=0;ii<fConversionSelection[iCut]->GetNumberOfPhotons();ii++){
530
531 AliAODConversionPhoton *gamma=fConversionSelection[iCut]->GetPhoton(ii);
532
533 Double_t val[knbinsGamma];
534 val[0]=gamma->Pt();
535 val[1]=GetRadialBin(gamma->GetConversionRadius());
536 val[2]=GetPhotonPhiwrtRP(gamma,iEP);
537 val[3]=fCentralityBin;
538 val[4]=Int_t(iEP);
539
540 hGamma[iCut]->Fill(val,weight);
541
542 if(iCut==0&&fFillQA)hGammadNdPhi->Fill(val);
543 }
544}
545
546//________________________________________________________________________
547void AliAnalysisTaskPi0v2::ProcessQA(){
548
549 AliStack *fMCStack=NULL;
550 if(fMCEvent)fMCStack=fMCEvent->Stack();
551
552 // Multiplicity
553
554 Double_t multcharged=fConversionSelection[0]->GetNumberOfChargedTracks(fInputEvent);
555 Double_t multVZERO=fConversionSelection[0]->GetVZEROMult(fInputEvent);
556 Double_t multSPD=fConversionSelection[0]->GetSPDMult(fInputEvent);
557
558 hMultChargedvsNGamma->Fill(multcharged,fConversionSelection[0]->GetNumberOfPhotons());
559 hMultChargedvsVZERO->Fill(multcharged,multVZERO);
560 hMultChargedvsSPD->Fill(multcharged,multSPD);
561
562 // Efficiency Purity
563
564 Int_t photonMultiplicity=fConversionSelection[0]->GetMultiplicity(fInputEvent);
565 Float_t weight=1; // Weight for Multiplicity
566 if(fWeightMultiplicity){
567 weight*=1/Float_t(photonMultiplicity);
568 }
569
570 Double_t valdPhi[knbinsGammaMult];
571 Double_t val[knbinsGammaMult];
572
573 Int_t dNdPhi[fNBinsPhi];
574 Int_t ncharged;
575
576 for(Int_t iEP=0;iEP<fNEPMethods;iEP++){
577 GetChargeddNdPhi(&dNdPhi[0],ncharged,iEP);
578
579 // Reco
580 for(Int_t ii=0;ii<fConversionSelection[0]->GetNumberOfPhotons();ii++){
581 AliAODConversionPhoton *gamma=fConversionSelection[0]->GetPhoton(ii);
582 val[0]=gamma->Pt();
583 val[1]=GetRadialBin(gamma->GetConversionRadius());
584 val[2]=ncharged;
585 val[3]=fCentralityBin;
586
587 valdPhi[0]=gamma->Pt();
588 valdPhi[1]=GetRadialBin(gamma->GetConversionRadius());
589 valdPhi[2]=dNdPhi[GetPhotonPhiBin(gamma,iEP)];
590 valdPhi[3]=fCentralityBin;
591
592 hGammaMult[iEP]->Fill(val,weight);
593 hGammaMultdPhi[iEP]->Fill(valdPhi,weight);
594
595 // Gamma Phi
596 hGammaPhi[fCentralityBin]->Fill(gamma->Pt(),GetRadialBin(gamma->GetConversionRadius()),gamma->Phi());
597 hGammaMultCent->Fill(fCentrality,Float_t(fConversionSelection[0]->GetNumberOfPhotons()));
598
599 if(fMCStack){
600 if(gamma->IsTruePhoton(fMCStack)){
601 hGammaMultRECOTRUE->Fill(val,weight);
602 hGammaMultdPhiRECOTRUE->Fill(valdPhi,weight);
603 }
604 }
605 }
606
607 // MC Truth
608 if(fMCEvent){
609 for(Int_t i = 0; i < fMCStack->GetNprimary(); i++) {
610 TParticle* particle = (TParticle *)fMCStack->Particle(i);
611 if (!particle) continue;
612 if(fConversionCuts->PhotonIsSelectedMC(particle,fMCStack)){
613 TParticle *daughter=(TParticle *)fMCStack->Particle(particle->GetDaughter(0));
614 if(daughter){
615 val[0]=particle->Pt();
616 val[1]=GetRadialBin(daughter->R());
617 val[2]=ncharged;
618 val[3]=fCentralityBin;
619
620 valdPhi[0]=particle->Pt();
621 valdPhi[1]=GetRadialBin(daughter->R());
622 valdPhi[2]=dNdPhi[GetPhiBin(GetMCPhotonPhiwrtRP(particle,EEventPlaneMethod(iEP)))];
623 valdPhi[3]=fCentralityBin;
624
625 hGammaMultTRUE->Fill(val,weight);
626 hGammaMultdPhiTRUE->Fill(valdPhi,weight);
627 }
628 }
629 }
630 }
631 }
632}
633
634//________________________________________________________________________
635void AliAnalysisTaskPi0v2::GetPhotondNdPhi(Int_t *dNdPhi,Int_t iEP,Int_t iCut){
636
637 for(Int_t iPhi=0;iPhi<fNBinsPhi;iPhi++)dNdPhi[iPhi]=0;
638
639 for(Int_t ii=0;ii<fConversionSelection[iCut]->GetNumberOfPhotons();ii++){
640 AliAODConversionPhoton *gamma=fConversionSelection[iCut]->GetPhoton(ii);
641 Int_t phibin=GetPhotonPhiBin(gamma,iEP);
642 dNdPhi[phibin]++;
643 }
644}
645
646//________________________________________________________________________
647void AliAnalysisTaskPi0v2::GetChargeddNdPhi(Int_t *dNdPhi,Int_t &ntot,Int_t iEP){
648
649 for(Int_t iPhi=0;iPhi<fNBinsPhi;iPhi++)dNdPhi[iPhi]=0;
650 ntot=0;
651
652 Double_t val[3];
653 val[1]=fCentralityBin;
654 val[2]=Int_t(iEP);
655
656 for(Int_t iTracks = 0; iTracks < fInputEvent->GetNumberOfTracks(); iTracks++){
657 AliVTrack* currentTrack = (AliVTrack*)fInputEvent->GetTrack(iTracks);
658 if(!currentTrack) continue;
659 if(TMath::Abs(currentTrack->Eta())>fEtaMax)continue;
660
661 Double_t phiwrt=GetChargedPhiwrtRP(currentTrack,EEventPlaneMethod(iEP));
662 Int_t phibin=GetPhiBin(phiwrt);
663
664 val[0]=phiwrt;
665 hCharged->Fill(val);
666
667 dNdPhi[phibin]++;
668 ntot++;
669 }
670}
671
672//________________________________________________________________________
673Int_t AliAnalysisTaskPi0v2::GetPhiBin(Double_t phiwrt){
674 Double_t binrange=TMath::Pi()/(2*Double_t(fNBinsPhi));
675 for(Int_t iPhi=0;iPhi<fNBinsPhi;iPhi++){
676 if(phiwrt>=(binrange*iPhi)&&phiwrt<(binrange*(iPhi+1)))return iPhi;
677 }
678 return -1;
679}
680
681//________________________________________________________________________
682Int_t AliAnalysisTaskPi0v2::GetPhotonPhiBin(AliAODConversionPhoton *gamma,Int_t iEP){
683 Double_t phiwrt=GetPhotonPhiwrtRP(gamma,EEventPlaneMethod(iEP));
684 return GetPhiBin(phiwrt);
685}
686
687//________________________________________________________________________
688Double_t AliAnalysisTaskPi0v2::GetPi0PhiwrtRP(AliAODConversionMother *pi0,EEventPlaneMethod iEP){
689
690 AliAODConversionPhoton *gamma0=dynamic_cast<AliAODConversionPhoton*>(fConversionGammas->At(pi0->GetLabel1()));
691 AliAODConversionPhoton *gamma1=dynamic_cast<AliAODConversionPhoton*>(fConversionGammas->At(pi0->GetLabel2()));
692
693 Double_t EPAngle=GetEventPlaneAngle(iEP,pi0->Eta(),gamma0,gamma1);
694
695 return TMath::ACos(TMath::Abs(TMath::Cos(pi0->Phi()-EPAngle)));
696}
697
698//________________________________________________________________________
699Double_t AliAnalysisTaskPi0v2::GetPhotonPhiwrtRP(AliAODConversionPhoton *gamma,EEventPlaneMethod iEP){
700
701 Double_t EPAngle=GetEventPlaneAngle(iEP,gamma->Eta(),gamma,NULL);
702
703 return TMath::ACos(TMath::Abs(TMath::Cos(gamma->Phi()-EPAngle)));
704}
705
706//________________________________________________________________________
707Double_t AliAnalysisTaskPi0v2::GetMCPhotonPhiwrtRP(TParticle *gamma,EEventPlaneMethod iEP){
708
709 Double_t EPAngle=GetEventPlaneAngle(iEP,gamma->Eta(),NULL,NULL);
710
711 return TMath::ACos(TMath::Abs(TMath::Cos(gamma->Phi()-EPAngle)));
712}
713//________________________________________________________________________
714Double_t AliAnalysisTaskPi0v2::GetChargedPhiwrtRP(AliVTrack *track,EEventPlaneMethod iEP){
715
716 Double_t EPAngle=GetEventPlaneAngle(iEP,track->Eta(),NULL,NULL);
717
718 return TMath::ACos(TMath::Abs(TMath::Cos(track->Phi()-EPAngle)));
719}
720
721//________________________________________________________________________
722void AliAnalysisTaskPi0v2::Terminate(Option_t *)
723{
724
725}
726
727//________________________________________________________________________
728void AliAnalysisTaskPi0v2::ProcessEventPlane()
729{
730
731 if(!fMCEvent&&fConversionCuts->IsHeavyIon()){
732
733 // TPC EP
734 Double_t PsiRP1=fEP->GetQsub1()->Phi()/2;
735 Double_t PsiRP2=fEP->GetQsub2()->Phi()/2;
736 Double_t dPsi=fEP->GetQsubRes();
737 Double_t EPTPC=GetEventPlaneAngle(kTPC);
738
739 hRPTPC->Fill(fCentrality,EPTPC);
740 hRPTPCAC->Fill(PsiRP1,PsiRP2);
741 hCos2TPC->Fill(fCentrality,TMath::Cos(2.*dPsi));
742
743 // TPC Eta Gap
744 hRPTPCEtaA->Fill(fCentrality,fRPTPCEtaA);
745 hRPTPCEtaC->Fill(fCentrality,fRPTPCEtaC);
746 hRPTPCEtaAC->Fill(fRPTPCEtaA,fRPTPCEtaC);
747 dPsi=fRPTPCEtaA-fRPTPCEtaC;
748 hCos2TPCEta->Fill(fCentrality,TMath::Cos(2.*dPsi));
749
750 // V0
751 hCos2V0ATPC->Fill(fCentrality,TMath::Cos(2*(EPTPC-fRPV0A)));
752 hCos2V0CTPC->Fill(fCentrality,TMath::Cos(2*(EPTPC-fRPV0C)));
753 hCos2V0AC->Fill(fCentrality,TMath::Cos(2*(fRPV0C-fRPV0A)));
754
755 hRPV0A->Fill(fCentrality,fRPV0A);
756 hRPV0C->Fill(fCentrality,fRPV0C);
757
758 hRPV0ATPC->Fill(fRPV0A,EPTPC);
759 hRPV0CTPC->Fill(fRPV0C,EPTPC);
760 hRPV0AC->Fill(fRPV0A,fRPV0C);
761 }
762}
763
764//________________________________________________________________________
765TVector2 AliAnalysisTaskPi0v2::GetEPContribution(AliAODConversionPhoton *gamma){
766 TVector2 q;
767 for(Int_t ii=0;ii<2;ii++){
768 AliVTrack *fCurrentTrack=AliConversionCuts::GetTrack(fInputEvent,gamma->GetTrackLabel(ii));
769 TVector2 qtrack(fEP->GetQContributionX(fCurrentTrack),fEP->GetQContributionY(fCurrentTrack));
770 q+=qtrack;
771 }
772 return q;
773}
774
775//________________________________________________________________________
776Double_t AliAnalysisTaskPi0v2::GetCorrectedTPCEPAngle(AliAODConversionPhoton *gamma0,AliAODConversionPhoton *gamma1){
777 // Correct Event Plane for Dilepton Tracks
778 TVector2 q0(*fEP->GetQVector());
779 if(gamma0)q0-=GetEPContribution(gamma0);
780 if(gamma1)q0-=GetEPContribution(gamma1);
781 Double_t EPangle=q0.Phi()/2;
782 EPangle=ApplyFlatteningTPC(EPangle,fCentrality);
783
784 return EPangle;
785}
786
787//________________________________________________________________________
788Double_t AliAnalysisTaskPi0v2::GetTPCSubEPEta(Double_t etamin,Double_t etamax){
789 TVector2 q;
790 for(Int_t ii=0;ii<fInputEvent->GetNumberOfTracks();ii++){
791 AliVTrack *fCurrentTrack=dynamic_cast<AliVTrack*>(fInputEvent->GetTrack(ii));
792 if(!fCurrentTrack)continue;
793 if(fCurrentTrack->Eta()>=etamin&&fCurrentTrack->Eta()<=etamax){
794 TVector2 qtrack(fEP->GetQContributionX(fCurrentTrack),fEP->GetQContributionY(fCurrentTrack));
795 q+=qtrack;
796 }
797 }
798 return q.Phi()/2;
799}
800
801//________________________________________________________________________
802Double_t AliAnalysisTaskPi0v2::GetEventPlaneAngle(EEventPlaneMethod EPmethod,Double_t eta,AliAODConversionPhoton *gamma0,AliAODConversionPhoton *gamma1)
803{
804 // If arguments are not null, the contribution of these photons is subtracted from the TPC EP
805
806 if(fConversionCuts->IsHeavyIon()){
807
808 // For MC select random EP angle in order to avoid correlations due to azimuth dependent efficiencies (ITS holes)
809 if(fMCEvent){
810 return fRandomizer->Uniform(2*TMath::Pi());
811 }
812
813 switch(EPmethod){
814 case kTPC:
815 return GetCorrectedTPCEPAngle(gamma0,gamma1);
816 case kTPCEtaGap:
817 if(eta<0)return fRPTPCEtaA; // Use opposite EP
818 else return fRPTPCEtaC;
819 case kV0A:
820 return fRPV0A;
821 case kV0C:
822 return fRPV0C;
823 default:
824 return 0;
825 }
826 }
827
828 // NO EP in pp mode
829 return 0;
830}
831
832///________________________________________________________________________
833Bool_t AliAnalysisTaskPi0v2::SetCentrality(){
834
835 // Set centrality bin for current event
836
837 if(!fConversionCuts->IsHeavyIon()){
838 fCentrality=0;
839 fCentralityBin=0;
840 return kTRUE;
841 }
842
843 fCentrality=fConversionCuts->GetCentrality(fInputEvent);
844
845 if(fNCentralityBins>1){
846 for(fCentralityBin=0;fCentralityBin<fNCentralityBins;fCentralityBin++){
847 if(fCentrality>=fCentralityBins[fCentralityBin]&&fCentrality<fCentralityBins[fCentralityBin+1])return kTRUE;
848 }
849 return kFALSE;
850 }
851 fCentralityBin=0;
852 return kTRUE;
853}
854
855//________________________________________________________________________
856void AliAnalysisTaskPi0v2::SetCentralityBins(Double_t *bins,Int_t nbins)
857{
858 // Set Centrality bins for analysis
859
860 fNCentralityBins=nbins;
861
862 if(fCentralityBins){
863 delete[] fCentralityBins;
864 fCentralityBins=NULL;
865 }
866
867 fCentralityBins=new Double_t[fNCentralityBins+1];
868 for(Int_t ii=0;ii<=fNCentralityBins;ii++){
869 fCentralityBins[ii]=bins[ii];
870 }
871}
872
873//________________________________________________________________________
874void AliAnalysisTaskPi0v2::SetCuts(TString *cutarray,Int_t ncuts){
875
876 if(fCutList){
877 delete fCutList;
878 fCutList=NULL;
879 }
880
881 fNCuts=ncuts;
882
883 fCutList = new TList();
884 fCutList->SetOwner(kTRUE);
885 AliConversionCuts **analysisCuts = new AliConversionCuts*[fNCuts];
886 for(Int_t i = 0; i<fNCuts; i++){
887 analysisCuts[i] = new AliConversionCuts();
888 analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
889 fCutList->Add(analysisCuts[i]);
890 }
891}
892
893//________________________________________________________________________
894Int_t AliAnalysisTaskPi0v2::GetRadialBin(Double_t radius){
895 if(fNRadialBins<=1)return 0;
896 for(Int_t ii=0;ii<fNRadialBins;ii++){
897 if(radius>fRadialBins[ii]){
898 if(radius<fRadialBins[ii+1]){return ii;}
899 }
900 }
901 return -1;
902}
903
904
905
906//________________________________________________________________________
907void AliAnalysisTaskPi0v2::SetRadialBins(Float_t *bins,Int_t nbins)
908{
909 // Set Centrality bins for analysis
910
911 fNRadialBins=nbins;
912
913 if(fRadialBins){
914 delete[] fRadialBins;
915 fRadialBins=NULL;
916 }
917
918 fRadialBins=new Double_t[fNRadialBins+1];
919 for(Int_t ii=0;ii<=fNRadialBins;ii++){
920 fRadialBins[ii]=bins[ii];
921 }
922}
923
924//___________________________________________________________________________
925Int_t AliAnalysisTaskPi0v2::GetRunIndex(Int_t run){
926
927 switch(run){
928 case 139517 : return 137;
929 case 139514 : return 136;
930 case 139513 : return 135;
931 case 139511 : return 134;
932 case 139510 : return 133;
933 case 139507 : return 132;
934 case 139505 : return 131;
935 case 139504 : return 130;
936 case 139503 : return 129;
937 case 139470 : return 128;
938 case 139467 : return 127;
939 case 139466 : return 126;
940 case 139465 : return 125;
941 case 139440 : return 124;
942 case 139439 : return 123;
943 case 139438 : return 122;
944 case 139437 : return 121;
945 case 139360 : return 120;
946 case 139329 : return 119;
947 case 139328 : return 118;
948 case 139314 : return 117;
949 case 139311 : return 116;
950 case 139310 : return 115;
951 case 139309 : return 114;
952 case 139308 : return 113;
953 case 139173 : return 112;
954 case 139172 : return 111;
955 case 139110 : return 110;
956 case 139107 : return 109;
957 case 139105 : return 108;
958 case 139104 : return 107;
959 case 139042 : return 106;
960 case 139038 : return 105;
961 case 139037 : return 104;
962 case 139036 : return 103;
963 case 139029 : return 102;
964 case 139028 : return 101;
965 case 138983 : return 100;
966 case 138982 : return 99;
967 case 138980 : return 98;
968 case 138979 : return 97;
969 case 138978 : return 96;
970 case 138977 : return 95;
971 case 138976 : return 94;
972 case 138973 : return 93;
973 case 138972 : return 92;
974 case 138965 : return 91;
975 case 138924 : return 90;
976 case 138872 : return 89;
977 case 138871 : return 88;
978 case 138870 : return 87;
979 case 138837 : return 86;
980 case 138830 : return 85;
981 case 138828 : return 84;
982 case 138826 : return 83;
983 case 138796 : return 82;
984 case 138795 : return 81;
985 case 138742 : return 80;
986 case 138732 : return 79;
987 case 138730 : return 78;
988 case 138666 : return 77;
989 case 138662 : return 76;
990 case 138653 : return 75;
991 case 138652 : return 74;
992 case 138638 : return 73;
993 case 138624 : return 72;
994 case 138621 : return 71;
995 case 138583 : return 70;
996 case 138582 : return 69;
997 case 138579 : return 68;
998 case 138578 : return 67;
999 case 138534 : return 66;
1000 case 138469 : return 65;
1001 case 138442 : return 64;
1002 case 138439 : return 63;
1003 case 138438 : return 62;
1004 case 138396 : return 61;
1005 case 138364 : return 60;
1006 case 138359 : return 59;
1007 case 138275 : return 58;
1008 case 138225 : return 57;
1009 case 138201 : return 56;
1010 case 138200 : return 55;
1011 case 138197 : return 54;
1012 case 138192 : return 53;
1013 case 138190 : return 52;
1014 case 138154 : return 51;
1015 case 138153 : return 50;
1016 case 138151 : return 49;
1017 case 138150 : return 48;
1018 case 138126 : return 47;
1019 case 138125 : return 46;
1020 case 137848 : return 45;
1021 case 137847 : return 44;
1022 case 137844 : return 43;
1023 case 137843 : return 42;
1024 case 137752 : return 41;
1025 case 137751 : return 40;
1026 case 137748 : return 39;
1027 case 137724 : return 38;
1028 case 137722 : return 37;
1029 case 137718 : return 36;
1030 case 137704 : return 35;
1031 case 137693 : return 34;
1032 case 137692 : return 33;
1033 case 137691 : return 32;
1034 case 137689 : return 31;
1035 case 137686 : return 30;
1036 case 137685 : return 29;
1037 case 137639 : return 28;
1038 case 137638 : return 27;
1039 case 137608 : return 26;
1040 case 137595 : return 25;
1041 case 137549 : return 24;
1042 case 137546 : return 23;
1043 case 137544 : return 22;
1044 case 137541 : return 21;
1045 case 137539 : return 20;
1046 case 137531 : return 19;
1047 case 137530 : return 18;
1048 case 137443 : return 17;
1049 case 137441 : return 16;
1050 case 137440 : return 15;
1051 case 137439 : return 14;
1052 case 137434 : return 13;
1053 case 137432 : return 12;
1054 case 137431 : return 11;
1055 case 137430 : return 10;
1056 case 137366 : return 9;
1057 case 137243 : return 8;
1058 case 137236 : return 7;
1059 case 137235 : return 6;
1060 case 137232 : return 5;
1061 case 137231 : return 4;
1062 case 137165 : return 3;
1063 case 137162 : return 2;
1064 case 137161 : return 1;
1065 default : return -1;
1066 }
1067}
1068
1069//____________________________________________________________________________
1070void AliAnalysisTaskPi0v2::GetV0EP(AliVEvent * event){
1071
1072 // Corrected VZERO EP (from AliAnalysisTaskPi0Flow)
1073
1074 //VZERO data
1075 AliESDVZERO* esdV0 = (AliESDVZERO*)event->GetVZEROData();
1076
1077 //reset Q vector info
1078 Double_t Qxa2 = 0, Qya2 = 0;
1079 Double_t Qxc2 = 0, Qyc2 = 0;
1080
1081 for (Int_t iv0 = 0; iv0 < 64; iv0++) {
1082 Double_t phiV0 = TMath::PiOver4()*(0.5 + iv0 % 8);
1083 Float_t multv0 = esdV0->GetMultiplicity(iv0);
1084 if (iv0 < 32){ // V0C
1085 Qxc2 += TMath::Cos(2*phiV0) * multv0*fV0Cpol/fMultV0->GetBinContent(iv0+1);
1086 Qyc2 += TMath::Sin(2*phiV0) * multv0*fV0Cpol/fMultV0->GetBinContent(iv0+1);
1087 } else { // V0A
1088 Qxa2 += TMath::Cos(2*phiV0) * multv0*fV0Apol/fMultV0->GetBinContent(iv0+1);
1089 Qya2 += TMath::Sin(2*phiV0) * multv0*fV0Apol/fMultV0->GetBinContent(iv0+1);
1090 }
1091 }
1092
1093 Int_t iC = -1;
1094 // centrality bins
1095 if(fCentrality < 5) iC = 0;
1096 else if(fCentrality < 10) iC = 1;
1097 else if(fCentrality < 20) iC = 2;
1098 else if(fCentrality < 30) iC = 3;
1099 else if(fCentrality < 40) iC = 4;
1100 else if(fCentrality < 50) iC = 5;
1101 else if(fCentrality < 60) iC = 6;
1102 else if(fCentrality < 70) iC = 7;
1103 else iC = 8;
1104
1105 //grab for each centrality the proper histo with the Qx and Qy to do the recentering
1106 Double_t Qxamean2 = fMeanQ[iC][1][0];
1107 Double_t Qxarms2 = fWidthQ[iC][1][0];
1108 Double_t Qyamean2 = fMeanQ[iC][1][1];
1109 Double_t Qyarms2 = fWidthQ[iC][1][1];
1110
1111 Double_t Qxcmean2 = fMeanQ[iC][0][0];
1112 Double_t Qxcrms2 = fWidthQ[iC][0][0];
1113 Double_t Qycmean2 = fMeanQ[iC][0][1];
1114 Double_t Qycrms2 = fWidthQ[iC][0][1];
1115
1116 Double_t QxaCor2 = (Qxa2 - Qxamean2)/Qxarms2;
1117 Double_t QyaCor2 = (Qya2 - Qyamean2)/Qyarms2;
1118 Double_t QxcCor2 = (Qxc2 - Qxcmean2)/Qxcrms2;
1119 Double_t QycCor2 = (Qyc2 - Qycmean2)/Qycrms2;
1120
1121 fRPV0A = TMath::ATan2(QyaCor2, QxaCor2)/2.;
1122 fRPV0C = TMath::ATan2(QycCor2, QxcCor2)/2.;
1123
1124 while(fRPV0A<0)fRPV0A+=TMath::Pi() ;
1125 while(fRPV0A>TMath::Pi())fRPV0A-=TMath::Pi() ;
1126 fRPV0A=ApplyFlatteningV0A(fRPV0A,fCentrality) ;
1127 while(fRPV0A<0)fRPV0A+=TMath::Pi() ;
1128 while(fRPV0A>TMath::Pi())fRPV0A-=TMath::Pi() ;
1129
1130 while(fRPV0C<0)fRPV0C+=TMath::Pi() ;
1131 while(fRPV0C>TMath::Pi())fRPV0C-=TMath::Pi() ;
1132 fRPV0C=ApplyFlatteningV0C(fRPV0C,fCentrality) ;
1133 while(fRPV0C<0)fRPV0C+=TMath::Pi() ;
1134 while(fRPV0C>TMath::Pi())fRPV0C-=TMath::Pi() ;
1135
1136}
1137
1138
1139//_____________________________________________________________________________
1140void AliAnalysisTaskPi0v2::OpenInfoCalibration(Int_t run){
1141 TString oadbfilename = "$ALICE_ROOT/OADB/PWGCF/VZERO/VZEROcalibEP.root";
1142 TFile *foadb = TFile::Open(oadbfilename.Data());
1143
1144 if(!foadb){
1145 printf("OADB file %s cannot be opened\n",oadbfilename.Data());
1146 return;
1147 }
1148
1149 AliOADBContainer *cont = (AliOADBContainer*) foadb->Get("hMultV0BefCorr");
1150 if(!cont){
1151 printf("OADB object hMultV0BefCorr is not available in the file\n");
1152 return;
1153 }
1154
1155 if(!(cont->GetObject(run))){
1156 printf("OADB object hMultV0BefCorr is not available for run %i (used run 137366)\n",run);
1157 run = 137366;
1158 }
1159 printf("Setting V0 calibration \n") ;
1160 fMultV0 = ((TH2F *) cont->GetObject(run))->ProfileX();
1161
1162 TF1 *fpol0 = new TF1("fpol0","pol0");
1163 fMultV0->Fit(fpol0,"","",0,31);
1164 fV0Cpol = fpol0->GetParameter(0);
1165 fMultV0->Fit(fpol0,"","",32,64);
1166 fV0Apol = fpol0->GetParameter(0);
1167
1168 for(Int_t iside=0;iside<2;iside++){
1169 for(Int_t icoord=0;icoord<2;icoord++){
1170 for(Int_t i=0;i < nCentrBinV0;i++){
1171 char namecont[100];
1172 if(iside==0 && icoord==0)
1173 snprintf(namecont,100,"hQxc2_%i",i);
1174 else if(iside==1 && icoord==0)
1175 snprintf(namecont,100,"hQxa2_%i",i);
1176 else if(iside==0 && icoord==1)
1177 snprintf(namecont,100,"hQyc2_%i",i);
1178 else if(iside==1 && icoord==1)
1179 snprintf(namecont,100,"hQya2_%i",i);
1180
1181 cont = (AliOADBContainer*) foadb->Get(namecont);
1182 if(!cont){
1183 printf("OADB object %s is not available in the file\n",namecont);
1184 return;
1185 }
1186
1187 if(!(cont->GetObject(run))){
1188 printf("OADB object %s is not available for run %i (used run 137366)\n",namecont,run);
1189 run = 137366;
1190 }
1191 fMeanQ[i][iside][icoord] = ((TH1F *) cont->GetObject(run))->GetMean();
1192 fWidthQ[i][iside][icoord] = ((TH1F *) cont->GetObject(run))->GetRMS();
1193
1194 //for v3
1195 if(iside==0 && icoord==0)
1196 snprintf(namecont,100,"hQxc3_%i",i);
1197 else if(iside==1 && icoord==0)
1198 snprintf(namecont,100,"hQxa3_%i",i);
1199 else if(iside==0 && icoord==1)
1200 snprintf(namecont,100,"hQyc3_%i",i);
1201 else if(iside==1 && icoord==1)
1202 snprintf(namecont,100,"hQya3_%i",i);
1203
1204 cont = (AliOADBContainer*) foadb->Get(namecont);
1205 if(!cont){
1206 printf("OADB object %s is not available in the file\n",namecont);
1207 return;
1208 }
1209
1210 if(!(cont->GetObject(run))){
1211 printf("OADB object %s is not available for run %i (used run 137366)\n",namecont,run);
1212 run = 137366;
1213 }
1214 // fMeanQv3[i][iside][icoord] = ((TH1F *) cont->GetObject(run))->GetMean();
1215 // fWidthQv3[i][iside][icoord] = ((TH1F *) cont->GetObject(run))->GetRMS();
1216
1217 }
1218 }
1219 }
1220}
1221
1222//____________________________________________________________________________
1223Double_t AliAnalysisTaskPi0v2::ApplyFlatteningV0A(Double_t phi, Double_t c){
1224 //LHC10h
1225
1226 Double_t v2c =0.;
1227 Double_t v2s =0;
1228
1229 if(fRunIndex==1){
1230 v2c=0.00271301 + (0.00340046)*c + (-0.000376369)*c*c + (8.56362e-06)*c*c*c + (-5.70955e-08)*c*c*c*c ;
1231 v2s=0.024278 + (0.000586593)*c + (-1.07748e-06)*c*c + (-5.01567e-07)*c*c*c + (3.65956e-09)*c*c*c*c ;
1232 }
1233 if(fRunIndex==2){
1234 v2c=0.111374 + (-0.0111331)*c + (0.000212342)*c*c + (-9.6774e-07)*c*c*c + (-2.54199e-09)*c*c*c*c ;
1235 v2s=-0.0252636 + (0.0152193)*c + (-0.000830793)*c*c + (1.5099e-05)*c*c*c + (-9.04687e-08)*c*c*c*c ;
1236 }
1237 if(fRunIndex==4){
1238 v2c=-0.0323376 + (0.0058672)*c + (-0.000257135)*c*c + (3.54353e-06)*c*c*c + (-1.50431e-08)*c*c*c*c ;
1239 v2s=0.0596758 + (-0.00520203)*c + (0.000178495)*c*c + (-2.94143e-06)*c*c*c + (1.55424e-08)*c*c*c*c ;
1240 }
1241 if(fRunIndex==5){
1242 v2c=-0.0415369 + (0.0098273)*c + (-0.000486907)*c*c + (8.00188e-06)*c*c*c + (-4.29107e-08)*c*c*c*c ;
1243 v2s=-0.0146049 + (0.00693627)*c + (-0.000348663)*c*c + (5.79331e-06)*c*c*c + (-3.33179e-08)*c*c*c*c ;
1244 }
1245 if(fRunIndex==6){
1246 v2c=0.251902 + (-0.0415419)*c + (0.0018553)*c*c + (-3.15181e-05)*c*c*c + (1.80073e-07)*c*c*c*c ;
1247 v2s=0.0771839 + (-0.00962468)*c + (0.000394453)*c*c + (-7.22485e-06)*c*c*c + (4.53562e-08)*c*c*c*c ;
1248 }
1249 if(fRunIndex==7){
1250 v2c=0.157038 + (-0.0204377)*c + (0.000847477)*c*c + (-1.4232e-05)*c*c*c + (8.04965e-08)*c*c*c*c ;
1251 v2s=0.0256337 + (-0.00242333)*c + (0.000140455)*c*c + (-3.5149e-06)*c*c*c + (2.44854e-08)*c*c*c*c ;
1252 }
1253 if(fRunIndex==8){
1254 v2c=-0.0916108 + (0.0146996)*c + (-0.000624971)*c*c + (9.63331e-06)*c*c*c + (-4.8574e-08)*c*c*c*c ;
1255 v2s=0.0332113 + (-0.00332475)*c + (0.000151857)*c*c + (-2.86187e-06)*c*c*c + (1.54796e-08)*c*c*c*c ;
1256 }
1257 if(fRunIndex==9){
1258 v2c=0.00925348 + (1.58854e-05)*c + (-3.85178e-05)*c*c + (6.75167e-07)*c*c*c + (-3.1502e-09)*c*c*c*c ;
1259 v2s=0.0249253 + (-0.00300902)*c + (0.000152653)*c*c + (-3.20923e-06)*c*c*c + (2.08818e-08)*c*c*c*c ;
1260 }
1261 if(fRunIndex==10){
1262 v2c=0.040903 + (0.00142711)*c + (-0.000264675)*c*c + (6.50624e-06)*c*c*c + (-4.46628e-08)*c*c*c*c ;
1263 v2s=0.0914196 + (-0.00957901)*c + (0.000408998)*c*c + (-7.49412e-06)*c*c*c + (4.63915e-08)*c*c*c*c ;
1264 }
1265 if(fRunIndex==11){
1266 v2c=-0.0130133 + (0.00296474)*c + (-0.000183518)*c*c + (3.65755e-06)*c*c*c + (-2.24484e-08)*c*c*c*c ;
1267 v2s=-0.00449031 + (0.00365474)*c + (-0.000226455)*c*c + (4.49792e-06)*c*c*c + (-2.96018e-08)*c*c*c*c ;
1268 }
1269 if(fRunIndex==12){
1270 v2c=-0.021987 + (0.00506675)*c + (-0.000301768)*c*c + (5.66088e-06)*c*c*c + (-3.21447e-08)*c*c*c*c ;
1271 v2s=0.0830858 + (-0.00871637)*c + (0.000299212)*c*c + (-4.10705e-06)*c*c*c + (1.76888e-08)*c*c*c*c ;
1272 }
1273 if(fRunIndex==13){
1274 v2c=0.0502952 + (-0.00152478)*c + (-0.000106289)*c*c + (3.39398e-06)*c*c*c + (-2.45015e-08)*c*c*c*c ;
1275 v2s=0.0639936 + (-0.00775191)*c + (0.000347449)*c*c + (-6.46454e-06)*c*c*c + (4.00439e-08)*c*c*c*c ;
1276 }
1277 if(fRunIndex==14){
1278 v2c=-0.116922 + (0.0122105)*c + (-0.000373923)*c*c + (4.08464e-06)*c*c*c + (-1.05659e-08)*c*c*c*c ;
1279 v2s=0.128821 + (-0.0204399)*c + (0.000831009)*c*c + (-1.33147e-05)*c*c*c + (7.24876e-08)*c*c*c*c ;
1280 }
1281 if(fRunIndex==15){
1282 v2c=0.0853741 + (-0.00924528)*c + (0.00034257)*c*c + (-5.69551e-06)*c*c*c + (3.44998e-08)*c*c*c*c ;
1283 v2s=0.0205492 + (-0.00052015)*c + (-5.42051e-05)*c*c + (1.82754e-06)*c*c*c + (-1.58393e-08)*c*c*c*c ;
1284 }
1285 if(fRunIndex==16){
1286 v2c=-0.0233639 + (0.00703612)*c + (-0.000403314)*c*c + (7.45693e-06)*c*c*c + (-4.4522e-08)*c*c*c*c ;
1287 v2s=-0.0712605 + (0.0124469)*c + (-0.000543644)*c*c + (8.22618e-06)*c*c*c + (-4.1073e-08)*c*c*c*c ;
1288 }
1289 if(fRunIndex==17){
1290 v2c=0.0113215 + (-0.0019788)*c + (9.81834e-05)*c*c + (-1.30009e-06)*c*c*c + (5.03998e-09)*c*c*c*c ;
1291 v2s=0.0459851 + (-0.00170929)*c + (-5.96459e-05)*c*c + (2.90146e-06)*c*c*c + (-2.85325e-08)*c*c*c*c ;
1292 }
1293 if(fRunIndex==18){
1294 v2c=-0.0924521 + (0.0162203)*c + (-0.000702717)*c*c + (1.0269e-05)*c*c*c + (-4.632e-08)*c*c*c*c ;
1295 v2s=-0.0456605 + (0.00706138)*c + (-0.000243957)*c*c + (2.7335e-06)*c*c*c + (-8.83296e-09)*c*c*c*c ;
1296 }
1297 if(fRunIndex==19){
1298 v2c=0.00150974 + (0.0113356)*c + (-0.000923111)*c*c + (2.0452e-05)*c*c*c + (-1.36949e-07)*c*c*c*c ;
1299 v2s=0.0673845 + (-0.0187379)*c + (0.00112044)*c*c + (-2.19997e-05)*c*c*c + (1.33115e-07)*c*c*c*c ;
1300 }
1301 if(fRunIndex==20){
1302 v2c=0.0223891 + (-0.000921021)*c + (-3.56276e-05)*c*c + (9.95622e-07)*c*c*c + (-5.54486e-09)*c*c*c*c ;
1303 v2s=0.0657305 + (-0.00354234)*c + (2.90428e-05)*c*c + (5.82994e-07)*c*c*c + (-7.53656e-09)*c*c*c*c ;
1304 }
1305 if(fRunIndex==21){
1306 v2c=0.00969434 + (-0.00195313)*c + (3.56409e-05)*c*c + (2.15514e-07)*c*c*c + (-5.5161e-09)*c*c*c*c ;
1307 v2s=0.000236584 + (0.00458234)*c + (-0.000306408)*c*c + (6.04476e-06)*c*c*c + (-3.79354e-08)*c*c*c*c ;
1308 }
1309 if(fRunIndex==22){
1310 v2c=0.0973536 + (-0.00907174)*c + (0.000239041)*c*c + (-2.64399e-06)*c*c*c + (1.12038e-08)*c*c*c*c ;
1311 v2s=0.0506578 + (-0.00256521)*c + (1.53982e-05)*c*c + (4.18193e-07)*c*c*c + (-4.99622e-09)*c*c*c*c ;
1312 }
1313 if(fRunIndex==23){
1314 v2c=0.120714 + (-0.0190965)*c + (0.000887955)*c*c + (-1.57144e-05)*c*c*c + (9.22823e-08)*c*c*c*c ;
1315 v2s=-0.00178837 + (0.00377706)*c + (-0.000358619)*c*c + (9.54949e-06)*c*c*c + (-7.39141e-08)*c*c*c*c ;
1316 }
1317 if(fRunIndex==24){
1318 v2c=0.0593053 + (-0.00687397)*c + (0.000200607)*c*c + (-2.34483e-06)*c*c*c + (9.60642e-09)*c*c*c*c ;
1319 v2s=0.0235174 + (0.00219788)*c + (-0.000225065)*c*c + (4.9594e-06)*c*c*c + (-3.2695e-08)*c*c*c*c ;
1320 }
1321 if(fRunIndex==25){
1322 v2c=0.0551353 + (-0.00340874)*c + (2.26718e-05)*c*c + (4.1675e-07)*c*c*c + (-3.37457e-09)*c*c*c*c ;
1323 v2s=0.0110038 + (-0.000180543)*c + (3.28645e-05)*c*c + (-1.08119e-06)*c*c*c + (6.7157e-09)*c*c*c*c ;
1324 }
1325 if(fRunIndex==26){
1326 v2c=-0.0138673 + (0.00475006)*c + (-0.00030584)*c*c + (6.0294e-06)*c*c*c + (-3.69906e-08)*c*c*c*c ;
1327 v2s=0.0395065 + (-0.00331262)*c + (9.40592e-05)*c*c + (-1.1541e-06)*c*c*c + (4.38329e-09)*c*c*c*c ;
1328 }
1329 if(fRunIndex==27){
1330 v2c=0.0547504 + (-0.00669067)*c + (0.000269785)*c*c + (-4.75217e-06)*c*c*c + (2.96254e-08)*c*c*c*c ;
1331 v2s=0.101508 + (-0.0107211)*c + (0.000385539)*c*c + (-6.28713e-06)*c*c*c + (3.75932e-08)*c*c*c*c ;
1332 }
1333 if(fRunIndex==28){
1334 v2c=-0.102173 + (0.0158262)*c + (-0.000706394)*c*c + (1.14263e-05)*c*c*c + (-6.14646e-08)*c*c*c*c ;
1335 v2s=0.02013 + (9.82197e-05)*c + (-1.71572e-05)*c*c + (-7.69646e-08)*c*c*c + (1.45113e-09)*c*c*c*c ;
1336 }
1337 if(fRunIndex==29){
1338 v2c=0.100539 + (-0.0115012)*c + (0.000295367)*c*c + (-2.57489e-06)*c*c*c + (5.41526e-09)*c*c*c*c ;
1339 v2s=0.138027 + (-0.00342397)*c + (-0.000392569)*c*c + (1.34626e-05)*c*c*c + (-1.06683e-07)*c*c*c*c ;
1340 }
1341 if(fRunIndex==30){
1342 v2c=0.0695326 + (-0.00752755)*c + (0.000237872)*c*c + (-3.31513e-06)*c*c*c + (1.69653e-08)*c*c*c*c ;
1343 v2s=0.0290133 + (-0.00073999)*c + (-3.35653e-05)*c*c + (7.61029e-07)*c*c*c + (-2.95792e-09)*c*c*c*c ;
1344 }
1345 if(fRunIndex==31){
1346 v2c=-0.103171 + (0.0134058)*c + (-0.000631602)*c*c + (1.10298e-05)*c*c*c + (-6.41177e-08)*c*c*c*c ;
1347 v2s=-0.137648 + (0.0111259)*c + (-0.000400035)*c*c + (5.90535e-06)*c*c*c + (-3.24676e-08)*c*c*c*c ;
1348 }
1349 if(fRunIndex==32){
1350 v2c=0.0522957 + (-0.00718492)*c + (0.000288586)*c*c + (-5.14534e-06)*c*c*c + (3.31303e-08)*c*c*c*c ;
1351 v2s=0.0731685 + (-0.00469449)*c + (9.57641e-05)*c*c + (-1.11891e-06)*c*c*c + (6.13984e-09)*c*c*c*c ;
1352 }
1353 if(fRunIndex==33){
1354 v2c=0.0882828 + (-0.00916109)*c + (0.00027821)*c*c + (-3.61312e-06)*c*c*c + (1.76418e-08)*c*c*c*c ;
1355 v2s=0.123302 + (-0.0156303)*c + (0.000623468)*c*c + (-9.83474e-06)*c*c*c + (5.24687e-08)*c*c*c*c ;
1356 }
1357 if(fRunIndex==34){
1358 v2c=0.0848636 + (-0.0142808)*c + (0.000723319)*c*c + (-1.32627e-05)*c*c*c + (7.77064e-08)*c*c*c*c ;
1359 v2s=-0.138483 + (0.0268553)*c + (-0.0013585)*c*c + (2.56034e-05)*c*c*c + (-1.63065e-07)*c*c*c*c ;
1360 }
1361 if(fRunIndex==35){
1362 v2c=-0.0134575 + (0.00267433)*c + (-0.000168784)*c*c + (3.24012e-06)*c*c*c + (-1.92766e-08)*c*c*c*c ;
1363 v2s=0.0596795 + (-0.00350808)*c + (2.97282e-05)*c*c + (3.52115e-07)*c*c*c + (-4.17462e-09)*c*c*c*c ;
1364 }
1365 if(fRunIndex==36){
1366 v2c=0.0784813 + (-0.00718389)*c + (0.000145303)*c*c + (-1.40121e-06)*c*c*c + (7.52775e-09)*c*c*c*c ;
1367 v2s=-0.0439181 + (0.0101363)*c + (-0.000494075)*c*c + (8.43356e-06)*c*c*c + (-4.68549e-08)*c*c*c*c ;
1368 }
1369 if(fRunIndex==37){
1370 v2c=-0.0118549 + (0.00328751)*c + (-0.000207156)*c*c + (3.92678e-06)*c*c*c + (-2.30059e-08)*c*c*c*c ;
1371 v2s=0.0344878 + (-0.000363669)*c + (-7.62481e-05)*c*c + (1.67135e-06)*c*c*c + (-9.64004e-09)*c*c*c*c ;
1372 }
1373 if(fRunIndex==38){
1374 v2c=-0.0254171 + (0.00297494)*c + (-0.000122361)*c*c + (1.78861e-06)*c*c*c + (-7.87609e-09)*c*c*c*c ;
1375 v2s=-0.0321399 + (0.00705636)*c + (-0.000325906)*c*c + (5.17669e-06)*c*c*c + (-2.86143e-08)*c*c*c*c ;
1376 }
1377 if(fRunIndex==39){
1378 v2c=-0.0368721 + (0.00253133)*c + (-0.000136793)*c*c + (2.51874e-06)*c*c*c + (-1.46055e-08)*c*c*c*c ;
1379 v2s=-0.0757763 + (0.005728)*c + (-0.000217489)*c*c + (3.11936e-06)*c*c*c + (-1.56919e-08)*c*c*c*c ;
1380 }
1381 if(fRunIndex==40){
1382 v2c=0.0166884 + (0.000509256)*c + (-0.000117911)*c*c + (2.65839e-06)*c*c*c + (-1.60875e-08)*c*c*c*c ;
1383 v2s=0.0336629 + (-0.00123023)*c + (-4.90123e-05)*c*c + (2.17862e-06)*c*c*c + (-1.98445e-08)*c*c*c*c ;
1384 }
1385 if(fRunIndex==41){
1386 v2c=-0.00887514 + (0.00491344)*c + (-0.000325782)*c*c + (6.54782e-06)*c*c*c + (-4.14143e-08)*c*c*c*c ;
1387 v2s=0.0720999 + (-0.00650524)*c + (0.000214752)*c*c + (-3.19681e-06)*c*c*c + (1.64573e-08)*c*c*c*c ;
1388 }
1389 if(fRunIndex==42){
1390 v2c=-0.0495903 + (0.00428229)*c + (-0.000181618)*c*c + (2.3913e-06)*c*c*c + (-8.55469e-09)*c*c*c*c ;
1391 v2s=0.0127982 + (-0.00690065)*c + (0.000342707)*c*c + (-6.40515e-06)*c*c*c + (3.83973e-08)*c*c*c*c ;
1392 }
1393 if(fRunIndex==43){
1394 v2c=0.0196783 + (-0.000221524)*c + (-7.09556e-05)*c*c + (1.76517e-06)*c*c*c + (-1.19612e-08)*c*c*c*c ;
1395 v2s=0.0350359 + (-0.0020458)*c + (5.18727e-05)*c*c + (-6.58401e-07)*c*c*c + (1.10936e-09)*c*c*c*c ;
1396 }
1397 if(fRunIndex==44){
1398 v2c=-0.0196041 + (0.00257682)*c + (-0.000174854)*c*c + (3.2168e-06)*c*c*c + (-1.85569e-08)*c*c*c*c ;
1399 v2s=-0.0870278 + (0.00634528)*c + (-0.000273588)*c*c + (4.77944e-06)*c*c*c + (-2.98422e-08)*c*c*c*c ;
1400 }
1401 if(fRunIndex==45){
1402 v2c=-0.0105827 + (0.00271714)*c + (-0.000152357)*c*c + (2.71899e-06)*c*c*c + (-1.5435e-08)*c*c*c*c ;
1403 v2s=0.135579 + (-0.0135952)*c + (0.000424184)*c*c + (-5.12316e-06)*c*c*c + (1.84497e-08)*c*c*c*c ;
1404 }
1405 if(fRunIndex==46){
1406 v2c=-0.017755 + (0.00411137)*c + (-0.000231724)*c*c + (4.17641e-06)*c*c*c + (-2.39164e-08)*c*c*c*c ;
1407 v2s=-0.0198628 + (0.0010032)*c + (-4.9708e-05)*c*c + (7.12455e-07)*c*c*c + (-4.22231e-09)*c*c*c*c ;
1408 }
1409 if(fRunIndex==47){
1410 v2c=-0.0111387 + (0.00258591)*c + (-0.000194287)*c*c + (3.84062e-06)*c*c*c + (-2.37192e-08)*c*c*c*c ;
1411 v2s=-0.0478616 + (0.00365349)*c + (-0.000152037)*c*c + (2.24323e-06)*c*c*c + (-1.23308e-08)*c*c*c*c ;
1412 }
1413 if(fRunIndex==48){
1414 v2c=0.070344 + (-0.0137332)*c + (0.000638536)*c*c + (-1.188e-05)*c*c*c + (7.4801e-08)*c*c*c*c ;
1415 v2s=-0.00208823 + (-0.00314157)*c + (9.9553e-05)*c*c + (-8.73891e-07)*c*c*c + (-1.72944e-09)*c*c*c*c ;
1416 }
1417 if(fRunIndex==51){
1418 v2c=0.017789 + (-0.00598224)*c + (0.000328477)*c*c + (-6.80787e-06)*c*c*c + (4.46905e-08)*c*c*c*c ;
1419 v2s=-0.0894734 + (0.007351)*c + (-0.000291085)*c*c + (4.38654e-06)*c*c*c + (-2.32261e-08)*c*c*c*c ;
1420 }
1421 if(fRunIndex==52){
1422 v2c=0.0555614 + (-0.00430576)*c + (0.00011519)*c*c + (-1.81584e-06)*c*c*c + (1.19214e-08)*c*c*c*c ;
1423 v2s=0.0901479 + (-0.0091505)*c + (0.000343757)*c*c + (-5.71204e-06)*c*c*c + (3.31217e-08)*c*c*c*c ;
1424 }
1425 if(fRunIndex==53){
1426 v2c=0.0181761 + (-0.00186378)*c + (4.60531e-05)*c*c + (-7.59581e-07)*c*c*c + (5.52465e-09)*c*c*c*c ;
1427 v2s=0.0456713 + (-0.00526813)*c + (0.000232978)*c*c + (-3.98724e-06)*c*c*c + (2.10384e-08)*c*c*c*c ;
1428 }
1429 if(fRunIndex==54){
1430 v2c=0.0456225 + (-0.00611173)*c + (0.000228414)*c*c + (-3.70527e-06)*c*c*c + (2.10774e-08)*c*c*c*c ;
1431 v2s=-0.00553874 + (0.00501799)*c + (-0.000331323)*c*c + (6.96978e-06)*c*c*c + (-4.72472e-08)*c*c*c*c ;
1432 }
1433 if(fRunIndex==55){
1434 v2c=-0.0578879 + (0.00664986)*c + (-0.000335443)*c*c + (6.0086e-06)*c*c*c + (-3.54293e-08)*c*c*c*c ;
1435 v2s=-0.0675269 + (0.00407581)*c + (-0.000148063)*c*c + (1.922e-06)*c*c*c + (-8.8635e-09)*c*c*c*c ;
1436 }
1437 if(fRunIndex==56){
1438 v2c=-0.000898903 + (0.00199222)*c + (-0.000207257)*c*c + (4.46148e-06)*c*c*c + (-2.71361e-08)*c*c*c*c ;
1439 v2s=-0.103794 + (0.0198667)*c + (-0.000879042)*c*c + (1.43177e-05)*c*c*c + (-8.09672e-08)*c*c*c*c ;
1440 }
1441 if(fRunIndex==57){
1442 v2c=0.00953087 + (0.000947327)*c + (-3.45822e-05)*c*c + (-4.22294e-07)*c*c*c + (9.08727e-09)*c*c*c*c ;
1443 v2s=0.0662314 + (-0.0076136)*c + (0.000288482)*c*c + (-4.82931e-06)*c*c*c + (2.67751e-08)*c*c*c*c ;
1444 }
1445 if(fRunIndex==58){
1446 v2c=0.00682728 + (-0.000293407)*c + (-2.88037e-05)*c*c + (6.63629e-07)*c*c*c + (-3.93121e-09)*c*c*c*c ;
1447 v2s=0.0085365 + (0.00138919)*c + (-8.15068e-05)*c*c + (1.26437e-06)*c*c*c + (-7.50322e-09)*c*c*c*c ;
1448 }
1449 if(fRunIndex==59){
1450 v2c=-0.0185226 + (0.00377881)*c + (-0.000226383)*c*c + (4.22897e-06)*c*c*c + (-2.51383e-08)*c*c*c*c ;
1451 v2s=-0.222605 + (0.0100699)*c + (-0.00029256)*c*c + (3.3855e-06)*c*c*c + (-1.33648e-08)*c*c*c*c ;
1452 }
1453 if(fRunIndex==60){
1454 v2c=-0.00111322 + (0.00291259)*c + (-0.000215769)*c*c + (4.13355e-06)*c*c*c + (-2.36883e-08)*c*c*c*c ;
1455 v2s=0.033686 + (-0.00124534)*c + (1.52616e-05)*c*c + (-2.99302e-07)*c*c*c + (1.05198e-09)*c*c*c*c ;
1456 }
1457 if(fRunIndex==61){
1458 v2c=0.00041873 + (-0.000231092)*c + (1.19226e-05)*c*c + (-6.29604e-07)*c*c*c + (7.26582e-09)*c*c*c*c ;
1459 v2s=0.00435438 + (0.0025265)*c + (-0.000153028)*c*c + (2.71169e-06)*c*c*c + (-1.67127e-08)*c*c*c*c ;
1460 }
1461 if(fRunIndex==62){
1462 v2c=-0.0573899 + (0.00972389)*c + (-0.000495772)*c*c + (9.0103e-06)*c*c*c + (-5.45756e-08)*c*c*c*c ;
1463 v2s=0.00702691 + (0.00698504)*c + (-0.000541929)*c*c + (1.19971e-05)*c*c*c + (-8.26942e-08)*c*c*c*c ;
1464 }
1465 if(fRunIndex==63){
1466 v2c=0.047222 + (-0.00604037)*c + (0.0002092)*c*c + (-3.07144e-06)*c*c*c + (1.68579e-08)*c*c*c*c ;
1467 v2s=0.0568815 + (-0.00607571)*c + (0.000244243)*c*c + (-4.22237e-06)*c*c*c + (2.38928e-08)*c*c*c*c ;
1468 }
1469 if(fRunIndex==64){
1470 v2c=0.0464686 + (-0.00463482)*c + (0.000146379)*c*c + (-2.23008e-06)*c*c*c + (1.29549e-08)*c*c*c*c ;
1471 v2s=-0.00571995 + (0.00273485)*c + (-0.000131945)*c*c + (2.08986e-06)*c*c*c + (-1.26717e-08)*c*c*c*c ;
1472 }
1473 if(fRunIndex==65){
1474 v2c=0.0544794 + (-0.00616918)*c + (0.000173808)*c*c + (-2.01338e-06)*c*c*c + (8.78123e-09)*c*c*c*c ;
1475 v2s=0.0557002 + (-0.00701287)*c + (0.000340418)*c*c + (-6.12754e-06)*c*c*c + (3.3384e-08)*c*c*c*c ;
1476 }
1477 if(fRunIndex==66){
1478 v2c=0.0111176 + (0.000262254)*c + (-7.20349e-05)*c*c + (1.38755e-06)*c*c*c + (-6.65016e-09)*c*c*c*c ;
1479 v2s=0.0644184 + (-0.00518515)*c + (0.000174227)*c*c + (-2.7347e-06)*c*c*c + (1.29547e-08)*c*c*c*c ;
1480 }
1481 if(fRunIndex==67){
1482 v2c=-0.0481037 + (0.00823686)*c + (-0.000418053)*c*c + (7.40186e-06)*c*c*c + (-4.19609e-08)*c*c*c*c ;
1483 v2s=0.0560171 + (-0.002556)*c + (2.54529e-05)*c*c + (-5.22269e-09)*c*c*c + (-2.19983e-09)*c*c*c*c ;
1484 }
1485 if(fRunIndex==68){
1486 v2c=-0.0139809 + (0.00274305)*c + (-0.00017744)*c*c + (3.20286e-06)*c*c*c + (-1.74137e-08)*c*c*c*c ;
1487 v2s=-0.218673 + (0.00764806)*c + (-0.000222739)*c*c + (3.14656e-06)*c*c*c + (-1.73548e-08)*c*c*c*c ;
1488 }
1489 if(fRunIndex==69){
1490 v2c=0.0736938 + (-0.0122004)*c + (0.000498227)*c*c + (-8.08779e-06)*c*c*c + (4.52749e-08)*c*c*c*c ;
1491 v2s=-0.236254 + (0.00979816)*c + (-0.000280698)*c*c + (3.4103e-06)*c*c*c + (-1.47299e-08)*c*c*c*c ;
1492 }
1493 if(fRunIndex==70){
1494 v2c=-0.0366545 + (0.00399604)*c + (-0.000196876)*c*c + (3.49686e-06)*c*c*c + (-2.01142e-08)*c*c*c*c ;
1495 v2s=0.0220911 + (-0.000529459)*c + (9.10408e-06)*c*c + (-3.57334e-07)*c*c*c + (1.87909e-09)*c*c*c*c ;
1496 }
1497 if(fRunIndex==71){
1498 v2c=-0.039698 + (0.00654615)*c + (-0.000323865)*c*c + (5.50998e-06)*c*c*c + (-3.00866e-08)*c*c*c*c ;
1499 v2s=0.0650078 + (-0.00502636)*c + (0.000128459)*c*c + (-1.3877e-06)*c*c*c + (3.21917e-09)*c*c*c*c ;
1500 }
1501 if(fRunIndex==72){
1502 v2c=0.0594227 + (-0.00872513)*c + (0.000368684)*c*c + (-6.34986e-06)*c*c*c + (3.79523e-08)*c*c*c*c ;
1503 v2s=0.050348 + (-0.00512087)*c + (0.000228643)*c*c + (-4.40797e-06)*c*c*c + (2.62482e-08)*c*c*c*c ;
1504 }
1505 if(fRunIndex==73){
1506 v2c=0.131281 + (-0.0162697)*c + (0.000584063)*c*c + (-8.64349e-06)*c*c*c + (4.61009e-08)*c*c*c*c ;
1507 v2s=0.0271721 + (-0.000487039)*c + (1.60624e-05)*c*c + (-5.38029e-07)*c*c*c + (2.08589e-09)*c*c*c*c ;
1508 }
1509 if(fRunIndex==74){
1510 v2c=-0.0128651 + (0.00946426)*c + (-0.00055547)*c*c + (9.37245e-06)*c*c*c + (-4.78224e-08)*c*c*c*c ;
1511 v2s=0.0903476 + (-0.0164972)*c + (0.000927481)*c*c + (-1.88459e-05)*c*c*c + (1.20249e-07)*c*c*c*c ;
1512 }
1513 if(fRunIndex==75){
1514 v2c=0.0204147 + (-0.00139139)*c + (-2.47732e-05)*c*c + (1.1214e-06)*c*c*c + (-7.22533e-09)*c*c*c*c ;
1515 v2s=0.0214247 + (-0.000147929)*c + (-7.58241e-06)*c*c + (-7.05565e-08)*c*c*c + (-4.86528e-10)*c*c*c*c ;
1516 }
1517 if(fRunIndex==76){
1518 v2c=0.0469095 + (-0.00541793)*c + (0.000133238)*c*c + (-1.14334e-06)*c*c*c + (3.4826e-09)*c*c*c*c ;
1519 v2s=0.0304763 + (-0.00213957)*c + (8.27588e-05)*c*c + (-1.34495e-06)*c*c*c + (4.81134e-09)*c*c*c*c ;
1520 }
1521 if(fRunIndex==77){
1522 v2c=0.113822 + (-0.0154286)*c + (0.000574533)*c*c + (-8.67874e-06)*c*c*c + (4.76727e-08)*c*c*c*c ;
1523 v2s=0.010623 + (0.00513467)*c + (-0.000372733)*c*c + (7.40692e-06)*c*c*c + (-4.64093e-08)*c*c*c*c ;
1524 }
1525 if(fRunIndex==78){
1526 v2c=0.020658 + (-0.00274725)*c + (5.8872e-05)*c*c + (-3.43999e-07)*c*c*c + (-1.04273e-10)*c*c*c*c ;
1527 v2s=0.00219379 + (0.000972033)*c + (-1.41816e-05)*c*c + (-4.41149e-07)*c*c*c + (3.81037e-09)*c*c*c*c ;
1528 }
1529 if(fRunIndex==79){
1530 v2c=0.227324 + (-0.0264915)*c + (0.000915298)*c*c + (-1.3343e-05)*c*c*c + (7.1044e-08)*c*c*c*c ;
1531 v2s=0.0474072 + (-0.00939515)*c + (0.000730108)*c*c + (-1.80811e-05)*c*c*c + (1.3003e-07)*c*c*c*c ;
1532 }
1533 if(fRunIndex==81){
1534 v2c=-0.0368028 + (0.00368927)*c + (-0.000228321)*c*c + (4.48174e-06)*c*c*c + (-2.75298e-08)*c*c*c*c ;
1535 v2s=-0.238807 + (0.0108667)*c + (-0.000350649)*c*c + (4.95568e-06)*c*c*c + (-2.53022e-08)*c*c*c*c ;
1536 }
1537 if(fRunIndex==82){
1538 v2c=-0.0214072 + (0.00119258)*c + (-7.80215e-05)*c*c + (1.37682e-06)*c*c*c + (-7.31103e-09)*c*c*c*c ;
1539 v2s=-0.233825 + (0.00735109)*c + (-0.000142737)*c*c + (9.15013e-07)*c*c*c + (-1.98463e-10)*c*c*c*c ;
1540 }
1541 if(fRunIndex==83){
1542 v2c=-0.0716708 + (0.00855289)*c + (-0.00037249)*c*c + (6.41638e-06)*c*c*c + (-3.66318e-08)*c*c*c*c ;
1543 v2s=-0.329702 + (0.033231)*c + (-0.00131012)*c*c + (2.1101e-05)*c*c*c + (-1.18905e-07)*c*c*c*c ;
1544 }
1545 if(fRunIndex==84){
1546 v2c=-0.0259022 + (0.00711773)*c + (-0.000434778)*c*c + (8.94774e-06)*c*c*c + (-5.84408e-08)*c*c*c*c ;
1547 v2s=-0.0936287 + (-0.00455282)*c + (0.000390982)*c*c + (-7.95224e-06)*c*c*c + (4.95612e-08)*c*c*c*c ;
1548 }
1549 if(fRunIndex==85){
1550 v2c=0.0263223 + (-0.00475593)*c + (9.63236e-05)*c*c + (-4.95311e-07)*c*c*c + (-5.47428e-10)*c*c*c*c ;
1551 v2s=-0.220593 + (0.00823892)*c + (-0.000248477)*c*c + (3.28082e-06)*c*c*c + (-1.45098e-08)*c*c*c*c ;
1552 }
1553 if(fRunIndex==86){
1554 v2c=-0.0157962 + (0.00408922)*c + (-0.000282087)*c*c + (5.81961e-06)*c*c*c + (-3.6481e-08)*c*c*c*c ;
1555 v2s=0.0340979 + (-0.00139471)*c + (6.48561e-06)*c*c + (3.02184e-08)*c*c*c + (-1.05879e-09)*c*c*c*c ;
1556 }
1557 if(fRunIndex==87){
1558 v2c=0.0577366 + (-0.00471394)*c + (0.000107329)*c*c + (-7.28927e-07)*c*c*c + (-1.51199e-09)*c*c*c*c ;
1559 v2s=0.0350171 + (-0.00402687)*c + (0.00016046)*c*c + (-2.39625e-06)*c*c*c + (9.98146e-09)*c*c*c*c ;
1560 }
1561 if(fRunIndex==88){
1562 v2c=-0.00547971 + (0.00161905)*c + (-0.00013831)*c*c + (2.94662e-06)*c*c*c + (-1.92504e-08)*c*c*c*c ;
1563 v2s=-0.0722251 + (0.0127634)*c + (-0.000612284)*c*c + (1.07553e-05)*c*c*c + (-6.34054e-08)*c*c*c*c ;
1564 }
1565 if(fRunIndex==89){
1566 v2c=0.165574 + (-0.0241913)*c + (0.00101221)*c*c + (-1.62279e-05)*c*c*c + (8.77252e-08)*c*c*c*c ;
1567 v2s=0.0692954 + (-0.0147219)*c + (0.00068782)*c*c + (-1.11091e-05)*c*c*c + (5.37295e-08)*c*c*c*c ;
1568 }
1569 if(fRunIndex==91){
1570 v2c=-0.0334666 + (0.00583866)*c + (-0.000329883)*c*c + (6.08932e-06)*c*c*c + (-3.59853e-08)*c*c*c*c ;
1571 v2s=-0.229775 + (0.010061)*c + (-0.00028144)*c*c + (3.3231e-06)*c*c*c + (-1.40799e-08)*c*c*c*c ;
1572 }
1573 if(fRunIndex==95){
1574 v2c=-0.0993203 + (0.0150547)*c + (-0.000713531)*c*c + (1.22733e-05)*c*c*c + (-6.98669e-08)*c*c*c*c ;
1575 v2s=-0.217746 + (0.0068375)*c + (-0.000171962)*c*c + (2.05506e-06)*c*c*c + (-8.78703e-09)*c*c*c*c ;
1576 }
1577 if(fRunIndex==96){
1578 v2c=-0.0398801 + (0.00406028)*c + (-0.000214066)*c*c + (3.69056e-06)*c*c*c + (-1.99339e-08)*c*c*c*c ;
1579 v2s=-0.240195 + (0.00782806)*c + (-0.000131365)*c*c + (3.18193e-07)*c*c*c + (4.88009e-09)*c*c*c*c ;
1580 }
1581 if(fRunIndex==97){
1582 v2c=-0.0325538 + (0.00250986)*c + (-0.000136762)*c*c + (2.35572e-06)*c*c*c + (-1.23413e-08)*c*c*c*c ;
1583 v2s=-0.240415 + (0.010435)*c + (-0.000291848)*c*c + (3.35357e-06)*c*c*c + (-1.35833e-08)*c*c*c*c ;
1584 }
1585 if(fRunIndex==98){
1586 v2c=-0.00692338 + (0.00379335)*c + (-0.000332833)*c*c + (7.35052e-06)*c*c*c + (-4.89956e-08)*c*c*c*c ;
1587 v2s=-0.263042 + (0.0134291)*c + (-0.000465334)*c*c + (7.09595e-06)*c*c*c + (-3.90332e-08)*c*c*c*c ;
1588 }
1589 if(fRunIndex==101){
1590 v2c=-0.049409 + (0.0073433)*c + (-0.000376933)*c*c + (7.39945e-06)*c*c*c + (-4.967e-08)*c*c*c*c ;
1591 v2s=0.00564511 + (-0.00101958)*c + (0.000121077)*c*c + (-3.44889e-06)*c*c*c + (2.42733e-08)*c*c*c*c ;
1592 }
1593 if(fRunIndex==102){
1594 v2c=-0.0393185 + (0.00610331)*c + (-0.000286099)*c*c + (4.61693e-06)*c*c*c + (-2.44104e-08)*c*c*c*c ;
1595 v2s=-0.00898987 + (0.00463766)*c + (-0.000269122)*c*c + (5.46456e-06)*c*c*c + (-3.85677e-08)*c*c*c*c ;
1596 }
1597 if(fRunIndex==103){
1598 v2c=-0.0171815 + (0.00271616)*c + (-9.49063e-05)*c*c + (5.68888e-07)*c*c*c + (2.77777e-09)*c*c*c*c ;
1599 v2s=-0.0348362 + (0.00661647)*c + (-0.000196625)*c*c + (1.60914e-06)*c*c*c + (-4.34769e-09)*c*c*c*c ;
1600 }
1601 if(fRunIndex==104){
1602 v2c=0.0453893 + (-0.00220981)*c + (-6.34642e-05)*c*c + (2.5268e-06)*c*c*c + (-1.81319e-08)*c*c*c*c ;
1603 v2s=0.0131555 + (0.00115036)*c + (-5.70975e-05)*c*c + (5.20224e-07)*c*c*c + (-3.11332e-09)*c*c*c*c ;
1604 }
1605 if(fRunIndex==105){
1606 v2c=0.0052505 + (0.000251094)*c + (-3.95143e-05)*c*c + (5.6354e-07)*c*c*c + (-1.86603e-09)*c*c*c*c ;
1607 v2s=0.079882 + (-0.00841369)*c + (0.000342924)*c*c + (-5.64554e-06)*c*c*c + (2.89511e-08)*c*c*c*c ;
1608 }
1609 if(fRunIndex==106){
1610 v2c=0.0561693 + (-0.00693389)*c + (0.000243013)*c*c + (-3.7816e-06)*c*c*c + (2.21666e-08)*c*c*c*c ;
1611 v2s=0.115847 + (-0.0145063)*c + (0.000597908)*c*c + (-9.37418e-06)*c*c*c + (4.57406e-08)*c*c*c*c ;
1612 }
1613 if(fRunIndex==107){
1614 v2c=-0.160892 + (0.0259075)*c + (-0.00120084)*c*c + (1.99814e-05)*c*c*c + (-1.09708e-07)*c*c*c*c ;
1615 v2s=0.0922001 + (-0.00861465)*c + (0.000196846)*c*c + (-5.34939e-07)*c*c*c + (-1.20442e-08)*c*c*c*c ;
1616 }
1617 if(fRunIndex==108){
1618 v2c=0.0325605 + (-0.00795036)*c + (0.000436987)*c*c + (-9.37794e-06)*c*c*c + (6.59912e-08)*c*c*c*c ;
1619 v2s=0.0894763 + (-0.00694006)*c + (0.000214928)*c*c + (-3.05745e-06)*c*c*c + (1.33679e-08)*c*c*c*c ;
1620 }
1621 if(fRunIndex==109){
1622 v2c=-8.30505e-05 + (0.00304264)*c + (-0.000228365)*c*c + (4.43511e-06)*c*c*c + (-2.5788e-08)*c*c*c*c ;
1623 v2s=0.0276356 + (-0.000486726)*c + (-1.12152e-05)*c*c + (1.56258e-07)*c*c*c + (-2.17826e-09)*c*c*c*c ;
1624 }
1625 if(fRunIndex==112){
1626 v2c=0.0273919 + (-0.00243837)*c + (2.41667e-05)*c*c + (3.06367e-07)*c*c*c + (-3.9028e-09)*c*c*c*c ;
1627 v2s=0.0119428 + (0.00161447)*c + (-0.00010526)*c*c + (1.76239e-06)*c*c*c + (-1.12844e-08)*c*c*c*c ;
1628 }
1629 if(fRunIndex==114){
1630 v2c=0.114287 + (-0.016125)*c + (0.000673002)*c*c + (-1.11635e-05)*c*c*c + (6.39533e-08)*c*c*c*c ;
1631 v2s=0.0353709 + (-0.00169691)*c + (-3.93183e-05)*c*c + (2.08815e-06)*c*c*c + (-2.03159e-08)*c*c*c*c ;
1632 }
1633 if(fRunIndex==115){
1634 v2c=-0.0109276 + (0.00378322)*c + (-0.000252651)*c*c + (4.599e-06)*c*c*c + (-2.36568e-08)*c*c*c*c ;
1635 v2s=0.0185448 + (0.00102852)*c + (-6.04966e-05)*c*c + (5.89933e-07)*c*c*c + (-2.5921e-09)*c*c*c*c ;
1636 }
1637 if(fRunIndex==116){
1638 v2c=-0.117563 + (0.00405449)*c + (-0.000112806)*c*c + (4.01903e-06)*c*c*c + (-4.25358e-08)*c*c*c*c ;
1639 v2s=0.0280797 + (-0.00528592)*c + (0.000375784)*c*c + (-9.96825e-06)*c*c*c + (7.72147e-08)*c*c*c*c ;
1640 }
1641 if(fRunIndex==117){
1642 v2c=0.0301869 + (-0.000751512)*c + (-9.43748e-05)*c*c + (2.60697e-06)*c*c*c + (-1.73948e-08)*c*c*c*c ;
1643 v2s=0.0594442 + (-0.00499442)*c + (0.000170617)*c*c + (-2.80903e-06)*c*c*c + (1.51075e-08)*c*c*c*c ;
1644 }
1645 if(fRunIndex==118){
1646 v2c=0.0118942 + (-0.000669595)*c + (-3.6634e-05)*c*c + (1.40794e-06)*c*c*c + (-1.15253e-08)*c*c*c*c ;
1647 v2s=-0.0174395 + (0.00460355)*c + (-0.00023192)*c*c + (3.83941e-06)*c*c*c + (-2.30761e-08)*c*c*c*c ;
1648 }
1649 if(fRunIndex==119){
1650 v2c=0.018229 + (-0.00144698)*c + (-1.07295e-06)*c*c + (6.22528e-07)*c*c*c + (-4.90504e-09)*c*c*c*c ;
1651 v2s=0.0122761 + (0.000985351)*c + (-6.02936e-05)*c*c + (9.41867e-07)*c*c*c + (-6.59756e-09)*c*c*c*c ;
1652 }
1653 if(fRunIndex==120){
1654 v2c=-0.0968065 + (0.0131861)*c + (-0.000585224)*c*c + (9.82175e-06)*c*c*c + (-5.51313e-08)*c*c*c*c ;
1655 v2s=0.0719334 + (-0.00752912)*c + (0.000362464)*c*c + (-6.65114e-06)*c*c*c + (3.70295e-08)*c*c*c*c ;
1656 }
1657 if(fRunIndex==121){
1658 v2c=0.0296329 + (-0.00448004)*c + (0.00010239)*c*c + (-2.40647e-07)*c*c*c + (-5.4555e-09)*c*c*c*c ;
1659 v2s=0.0774925 + (-0.00635491)*c + (0.0001939)*c*c + (-3.04755e-06)*c*c*c + (1.66055e-08)*c*c*c*c ;
1660 }
1661 if(fRunIndex==122){
1662 v2c=0.00956939 + (0.000460793)*c + (-0.000105912)*c*c + (2.38286e-06)*c*c*c + (-1.446e-08)*c*c*c*c ;
1663 v2s=0.0152677 + (0.00345263)*c + (-0.000262399)*c*c + (5.5895e-06)*c*c*c + (-3.84974e-08)*c*c*c*c ;
1664 }
1665 if(fRunIndex==123){
1666 v2c=0.0495618 + (-0.00059293)*c + (-0.00025227)*c*c + (6.61088e-06)*c*c*c + (-4.24063e-08)*c*c*c*c ;
1667 v2s=0.0862227 + (-0.0190161)*c + (0.0010647)*c*c + (-1.94957e-05)*c*c*c + (1.08665e-07)*c*c*c*c ;
1668 }
1669 if(fRunIndex==124){
1670 v2c=-0.0408078 + (0.00565108)*c + (-0.00027528)*c*c + (5.2165e-06)*c*c*c + (-3.3244e-08)*c*c*c*c ;
1671 v2s=0.00169862 + (0.00371982)*c + (-0.000169757)*c*c + (2.07347e-06)*c*c*c + (-8.30142e-09)*c*c*c*c ;
1672 }
1673 if(fRunIndex==125){
1674 v2c=0.0285065 + (-0.00348583)*c + (0.000106189)*c*c + (-1.56726e-06)*c*c*c + (9.40552e-09)*c*c*c*c ;
1675 v2s=0.0363645 + (-0.00252206)*c + (0.000122553)*c*c + (-2.74066e-06)*c*c*c + (1.66856e-08)*c*c*c*c ;
1676 }
1677 if(fRunIndex==126){
1678 v2c=-0.00864424 + (0.00157555)*c + (-0.000127367)*c*c + (2.47447e-06)*c*c*c + (-1.42302e-08)*c*c*c*c ;
1679 v2s=-0.213003 + (0.00643156)*c + (-0.000121219)*c*c + (6.06838e-07)*c*c*c + (1.48381e-09)*c*c*c*c ;
1680 }
1681 if(fRunIndex==127){
1682 v2c=-0.0336884 + (0.00299236)*c + (-0.000147557)*c*c + (2.47908e-06)*c*c*c + (-1.33495e-08)*c*c*c*c ;
1683 v2s=-0.271925 + (0.0158755)*c + (-0.000583802)*c*c + (9.14376e-06)*c*c*c + (-5.05302e-08)*c*c*c*c ;
1684 }
1685 if(fRunIndex==128){
1686 v2c=-0.0639711 + (0.00822881)*c + (-0.000428416)*c*c + (7.98075e-06)*c*c*c + (-4.81995e-08)*c*c*c*c ;
1687 v2s=-0.22634 + (0.00908632)*c + (-0.00029)*c*c + (4.14095e-06)*c*c*c + (-2.15823e-08)*c*c*c*c ;
1688 }
1689 if(fRunIndex==129){
1690 v2c=0.0652753 + (-0.00529737)*c + (0.000133883)*c*c + (-1.63339e-06)*c*c*c + (7.88294e-09)*c*c*c*c ;
1691 v2s=0.00357003 + (0.00503303)*c + (-0.000316531)*c*c + (6.45793e-06)*c*c*c + (-4.44467e-08)*c*c*c*c ;
1692 }
1693 if(fRunIndex==131){
1694 v2c=0.117255 + (-0.0149009)*c + (0.000560178)*c*c + (-8.75754e-06)*c*c*c + (4.85024e-08)*c*c*c*c ;
1695 v2s=-0.00180786 + (0.00224258)*c + (-9.84056e-05)*c*c + (1.39525e-06)*c*c*c + (-7.63512e-09)*c*c*c*c ;
1696 }
1697 if(fRunIndex==132){
1698 v2c=0.0100287 + (-0.00051338)*c + (-2.7702e-05)*c*c + (5.0205e-07)*c*c*c + (-1.12682e-09)*c*c*c*c ;
1699 v2s=0.00556895 + (0.00211082)*c + (-0.000125448)*c*c + (2.30409e-06)*c*c*c + (-1.63703e-08)*c*c*c*c ;
1700 }
1701 if(fRunIndex==133){
1702 v2c=-0.00755913 + (0.00327793)*c + (-0.00024545)*c*c + (5.03417e-06)*c*c*c + (-3.11617e-08)*c*c*c*c ;
1703 v2s=0.0398477 + (-0.000728372)*c + (-2.63935e-05)*c*c + (7.85276e-07)*c*c*c + (-7.72703e-09)*c*c*c*c ;
1704 }
1705 if(fRunIndex==134){
1706 v2c=0.00344056 + (0.000171769)*c + (-8.00887e-05)*c*c + (1.76696e-06)*c*c*c + (-1.08652e-08)*c*c*c*c ;
1707 v2s=-0.200474 + (0.00736532)*c + (-0.000221809)*c*c + (2.91763e-06)*c*c*c + (-1.40827e-08)*c*c*c*c ;
1708 }
1709 if(fRunIndex==135){
1710 v2c=-0.0251554 + (0.00227996)*c + (-0.000133757)*c*c + (2.50604e-06)*c*c*c + (-1.48528e-08)*c*c*c*c ;
1711 v2s=-0.229223 + (0.0091122)*c + (-0.000245035)*c*c + (2.64826e-06)*c*c*c + (-9.57268e-09)*c*c*c*c ;
1712 }
1713 if(fRunIndex==136){
1714 v2c=-0.0308547 + (0.00209237)*c + (-0.000183034)*c*c + (4.54164e-06)*c*c*c + (-3.26557e-08)*c*c*c*c ;
1715 v2s=-0.241702 + (0.00988355)*c + (-0.000263014)*c*c + (2.73376e-06)*c*c*c + (-9.79971e-09)*c*c*c*c ;
1716 }
1717
1718/*
1719
1720 if(fRunIndex>=1 && fRunIndex<=45){ // 137161 -137848 //period 1
1721// v2c = 1.133676e-02-1.759454e-04*x-6.217033e-05*x*x+1.532507e-06*x*x*x-9.665265e-09*x*x*x*x;
1722// v2s = 2.173690e-02-8.410340e-04*x-1.391841e-05*x*x+5.003301e-07*x*x*x-4.175974e-09*x*x*x*x;
1723 v2c=2.295070e-02-1.050171e-03*x-3.432843e-05*x*x+1.156549e-06*x*x*x-7.781377e-09*x*x*x*x ;
1724 v2s=3.994886e-02-1.961957e-03*x+1.641958e-05*x*x+1.234025e-07*x*x*x-2.175698e-09*x*x*x*x ;
1725 }
1726 if(fRunIndex>=46 && fRunIndex<=58){ //138125-138275//period 2
1727 v2c =-1.394612e-03+1.333726e-03*x-1.115323e-04*x*x+2.103762e-06*x*x*x-1.186760e-08*x*x*x*x ;
1728 v2s =-7.788201e-03+7.178984e-04*x-4.553597e-05*x*x+6.941817e-07*x*x*x-4.465249e-09*x*x*x*x ;
1729 }
1730 if(fRunIndex>=59 && fRunIndex<=78){ //138359-138730//period 3
1731 v2c = 1.448809e-02+-7.894029e-04*x-2.693871e-05*x*x+7.944017e-07*x*x*x-4.272135e-09*x*x*x*x ;
1732 v2s = 3.442216e-03+-5.574560e-05*x-1.758343e-05*x*x+2.386547e-07*x*x*x-2.404210e-09*x*x*x*x ;
1733 }
1734 if(fRunIndex>=79){//period 4
1735 v2c =-1.553536e-03+7.761106e-04*x-9.529723e-05*x*x+2.026750e-06*x*x*x-1.217009e-08*x*x*x*x ;
1736 v2s =-7.787317e-02+3.843455e-03*x-1.221144e-04*x*x+1.478448e-06*x*x*x-7.442702e-09*x*x*x*x ;
1737 }
1738*/
1739
1740 return phi+v2c*TMath::Sin(2.*phi)-v2s*TMath::Cos(2.*phi) ;
1741
1742}
1743//____________________________________________________________________________
1744Double_t AliAnalysisTaskPi0v2::ApplyFlatteningV0C(Double_t phi, Double_t c){
1745 //LHC10h
1746
1747 Double_t v2c =0.;
1748 Double_t v2s =0;
1749
1750 if(fRunIndex==1){
1751 v2c=-0.038901 + (0.00539521)*c + (-0.000213271)*c*c + (3.01147e-06)*c*c*c + (-1.35977e-08)*c*c*c*c ;
1752 v2s=0.0693294 + (-0.00851319)*c + (0.000337408)*c*c + (-5.33097e-06)*c*c*c + (2.90883e-08)*c*c*c*c ;
1753 }
1754 if(fRunIndex==2){
1755 v2c=-0.0943695 + (0.0137829)*c + (-0.000569334)*c*c + (8.95522e-06)*c*c*c + (-4.78396e-08)*c*c*c*c ;
1756 v2s=0.00971525 + (-0.00181026)*c + (6.15798e-05)*c*c + (-8.51046e-07)*c*c*c + (4.57403e-09)*c*c*c*c ;
1757 }
1758 if(fRunIndex==4){
1759 v2c=-0.00585475 + (0.000307793)*c + (1.89109e-06)*c*c + (-3.14804e-07)*c*c*c + (3.54341e-09)*c*c*c*c ;
1760 v2s=0.00374568 + (-0.00123435)*c + (8.01941e-05)*c*c + (-1.70807e-06)*c*c*c + (1.14576e-08)*c*c*c*c ;
1761 }
1762 if(fRunIndex==5){
1763 v2c=-0.0301267 + (0.00206698)*c + (-4.82282e-05)*c*c + (4.39355e-07)*c*c*c + (-1.3747e-09)*c*c*c*c ;
1764 v2s=-0.0441947 + (0.00501774)*c + (-0.000181736)*c*c + (2.69235e-06)*c*c*c + (-1.37411e-08)*c*c*c*c ;
1765 }
1766 if(fRunIndex==6){
1767 v2c=0.190652 + (-0.0254912)*c + (0.00105707)*c*c + (-1.63872e-05)*c*c*c + (8.51113e-08)*c*c*c*c ;
1768 v2s=0.0803124 + (-0.0119724)*c + (0.000575997)*c*c + (-1.07892e-05)*c*c*c + (6.67214e-08)*c*c*c*c ;
1769 }
1770 if(fRunIndex==7){
1771 v2c=0.00875044 + (-0.00343147)*c + (0.0002198)*c*c + (-4.76292e-06)*c*c*c + (3.24844e-08)*c*c*c*c ;
1772 v2s=-0.081288 + (0.0146429)*c + (-0.000728327)*c*c + (1.32462e-05)*c*c*c + (-7.95362e-08)*c*c*c*c ;
1773 }
1774 if(fRunIndex==8){
1775 v2c=0.0381367 + (-0.0082048)*c + (0.000410914)*c*c + (-7.28851e-06)*c*c*c + (4.22008e-08)*c*c*c*c ;
1776 v2s=-0.023338 + (0.003381)*c + (-0.000167791)*c*c + (3.13998e-06)*c*c*c + (-1.97233e-08)*c*c*c*c ;
1777 }
1778 if(fRunIndex==9){
1779 v2c=-0.00734261 + (0.00121272)*c + (-5.31181e-05)*c*c + (7.88974e-07)*c*c*c + (-3.86356e-09)*c*c*c*c ;
1780 v2s=-0.0228333 + (0.00266052)*c + (-0.000112824)*c*c + (1.89515e-06)*c*c*c + (-1.07446e-08)*c*c*c*c ;
1781 }
1782 if(fRunIndex==10){
1783 v2c=-0.0456385 + (0.0102956)*c + (-0.000529689)*c*c + (9.62809e-06)*c*c*c + (-5.7301e-08)*c*c*c*c ;
1784 v2s=0.0142197 + (-0.0031677)*c + (0.000190389)*c*c + (-4.10911e-06)*c*c*c + (2.81791e-08)*c*c*c*c ;
1785 }
1786 if(fRunIndex==11){
1787 v2c=-0.00999249 + (0.00123939)*c + (-4.23895e-05)*c*c + (6.55415e-07)*c*c*c + (-3.98256e-09)*c*c*c*c ;
1788 v2s=-0.00809817 + (0.00199727)*c + (-0.000116718)*c*c + (2.33449e-06)*c*c*c + (-1.46407e-08)*c*c*c*c ;
1789 }
1790 if(fRunIndex==12){
1791 v2c=0.0183566 + (-0.0020805)*c + (6.79726e-05)*c*c + (-7.96708e-07)*c*c*c + (3.1214e-09)*c*c*c*c ;
1792 v2s=-0.00714715 + (0.000721418)*c + (-3.50355e-05)*c*c + (9.09996e-07)*c*c*c + (-7.92198e-09)*c*c*c*c ;
1793 }
1794 if(fRunIndex==13){
1795 v2c=0.0056915 + (0.000167211)*c + (-4.68018e-05)*c*c + (1.21891e-06)*c*c*c + (-8.70124e-09)*c*c*c*c ;
1796 v2s=0.00986223 + (-0.00247957)*c + (0.000133539)*c*c + (-2.51362e-06)*c*c*c + (1.5397e-08)*c*c*c*c ;
1797 }
1798 if(fRunIndex==14){
1799 v2c=0.0479623 + (-0.00948196)*c + (0.000551795)*c*c + (-1.05692e-05)*c*c*c + (6.38754e-08)*c*c*c*c ;
1800 v2s=-0.00438456 + (-0.00152976)*c + (2.00645e-05)*c*c + (9.91017e-07)*c*c*c + (-1.35776e-08)*c*c*c*c ;
1801 }
1802 if(fRunIndex==15){
1803 v2c=-0.00384088 + (0.00121719)*c + (-6.88665e-05)*c*c + (1.26138e-06)*c*c*c + (-7.13987e-09)*c*c*c*c ;
1804 v2s=0.0273017 + (-0.00314818)*c + (7.85031e-05)*c*c + (-9.07336e-08)*c*c*c + (-6.42421e-09)*c*c*c*c ;
1805 }
1806 if(fRunIndex==16){
1807 v2c=-0.0313175 + (0.00368025)*c + (-0.000126895)*c*c + (1.7443e-06)*c*c*c + (-9.12442e-09)*c*c*c*c ;
1808 v2s=-0.0812178 + (0.0117402)*c + (-0.000506287)*c*c + (8.36766e-06)*c*c*c + (-4.69866e-08)*c*c*c*c ;
1809 }
1810 if(fRunIndex==17){
1811 v2c=0.062893 + (-0.0128204)*c + (0.000591931)*c*c + (-1.00835e-05)*c*c*c + (5.72782e-08)*c*c*c*c ;
1812 v2s=-0.143105 + (0.0183921)*c + (-0.000787289)*c*c + (1.33237e-05)*c*c*c + (-7.71707e-08)*c*c*c*c ;
1813 }
1814 if(fRunIndex==18){
1815 v2c=-0.0029353 + (-0.0016097)*c + (0.000199618)*c*c + (-5.42008e-06)*c*c*c + (4.08897e-08)*c*c*c*c ;
1816 v2s=0.0482015 + (-0.00115233)*c + (3.96019e-05)*c*c + (-7.7193e-07)*c*c*c + (3.90004e-09)*c*c*c*c ;
1817 }
1818 if(fRunIndex==19){
1819 v2c=-0.0397728 + (0.0129154)*c + (-0.000703903)*c*c + (1.30972e-05)*c*c*c + (-7.93547e-08)*c*c*c*c ;
1820 v2s=0.181945 + (-0.0289807)*c + (0.00133685)*c*c + (-2.30866e-05)*c*c*c + (1.32746e-07)*c*c*c*c ;
1821 }
1822 if(fRunIndex==20){
1823 v2c=-0.00927858 + (0.0021033)*c + (-0.000113971)*c*c + (2.17508e-06)*c*c*c + (-1.36812e-08)*c*c*c*c ;
1824 v2s=0.0201645 + (-0.00331311)*c + (0.000164916)*c*c + (-3.08376e-06)*c*c*c + (1.92565e-08)*c*c*c*c ;
1825 }
1826 if(fRunIndex==21){
1827 v2c=-0.0254895 + (0.00388685)*c + (-0.000171429)*c*c + (2.87209e-06)*c*c*c + (-1.60625e-08)*c*c*c*c ;
1828 v2s=0.00237609 + (0.000664255)*c + (-4.93436e-05)*c*c + (8.86415e-07)*c*c*c + (-4.51913e-09)*c*c*c*c ;
1829 }
1830 if(fRunIndex==22){
1831 v2c=0.0308671 + (-0.00365824)*c + (0.000148483)*c*c + (-2.52946e-06)*c*c*c + (1.49859e-08)*c*c*c*c ;
1832 v2s=-0.00805023 + (0.00154452)*c + (-6.89678e-05)*c*c + (1.14317e-06)*c*c*c + (-5.98991e-09)*c*c*c*c ;
1833 }
1834 if(fRunIndex==23){
1835 v2c=0.122018 + (-0.0233973)*c + (0.00116385)*c*c + (-2.11232e-05)*c*c*c + (1.25836e-07)*c*c*c*c ;
1836 v2s=-0.0511879 + (0.00708815)*c + (-0.000227282)*c*c + (2.27542e-06)*c*c*c + (-4.83128e-09)*c*c*c*c ;
1837 }
1838 if(fRunIndex==24){
1839 v2c=0.00888562 + (-0.000839847)*c + (2.52839e-05)*c*c + (-2.28921e-07)*c*c*c + (1.01057e-11)*c*c*c*c ;
1840 v2s=0.00814357 + (-0.000788997)*c + (2.43948e-05)*c*c + (-2.99558e-07)*c*c*c + (1.21675e-09)*c*c*c*c ;
1841 }
1842 if(fRunIndex==25){
1843 v2c=0.00360833 + (0.000379977)*c + (-3.68193e-05)*c*c + (5.86861e-07)*c*c*c + (-2.01826e-09)*c*c*c*c ;
1844 v2s=-0.00642047 + (0.00060076)*c + (-2.42907e-05)*c*c + (4.12933e-07)*c*c*c + (-2.3557e-09)*c*c*c*c ;
1845 }
1846 if(fRunIndex==26){
1847 v2c=-0.0131971 + (0.00189761)*c + (-7.7551e-05)*c*c + (1.25884e-06)*c*c*c + (-7.38989e-09)*c*c*c*c ;
1848 v2s=0.00679891 + (-0.000614692)*c + (1.04559e-06)*c*c + (4.34463e-07)*c*c*c + (-4.75821e-09)*c*c*c*c ;
1849 }
1850 if(fRunIndex==27){
1851 v2c=-0.00887089 + (0.00126865)*c + (-7.15766e-05)*c*c + (1.43811e-06)*c*c*c + (-9.408e-09)*c*c*c*c ;
1852 v2s=-0.0302105 + (0.00349565)*c + (-0.000126108)*c*c + (1.835e-06)*c*c*c + (-9.27845e-09)*c*c*c*c ;
1853 }
1854 if(fRunIndex==28){
1855 v2c=-0.0187616 + (0.0021974)*c + (-5.11412e-05)*c*c + (2.94628e-07)*c*c*c + (7.4194e-10)*c*c*c*c ;
1856 v2s=-0.0119204 + (-0.00156879)*c + (0.000132746)*c*c + (-2.67371e-06)*c*c*c + (1.61097e-08)*c*c*c*c ;
1857 }
1858 if(fRunIndex==29){
1859 v2c=-0.121018 + (0.0228924)*c + (-0.00127093)*c*c + (2.48724e-05)*c*c*c + (-1.56121e-07)*c*c*c*c ;
1860 v2s=-0.0687197 + (0.015702)*c + (-0.000764074)*c*c + (1.38991e-05)*c*c*c + (-8.48812e-08)*c*c*c*c ;
1861 }
1862 if(fRunIndex==30){
1863 v2c=0.01582 + (-0.0012293)*c + (1.62768e-05)*c*c + (2.36172e-07)*c*c*c + (-3.67153e-09)*c*c*c*c ;
1864 v2s=-0.00875651 + (0.00135376)*c + (-5.16763e-05)*c*c + (6.17614e-07)*c*c*c + (-1.64099e-09)*c*c*c*c ;
1865 }
1866 if(fRunIndex==31){
1867 v2c=-0.044721 + (0.00756113)*c + (-0.000365759)*c*c + (6.15804e-06)*c*c*c + (-3.40075e-08)*c*c*c*c ;
1868 v2s=-0.0500828 + (-0.000695112)*c + (0.000130909)*c*c + (-3.38476e-06)*c*c*c + (2.41179e-08)*c*c*c*c ;
1869 }
1870 if(fRunIndex==32){
1871 v2c=0.0306236 + (-0.00443307)*c + (0.000196517)*c*c + (-3.48059e-06)*c*c*c + (2.10266e-08)*c*c*c*c ;
1872 v2s=0.00819214 + (-0.000560558)*c + (-5.03783e-06)*c*c + (4.37862e-07)*c*c*c + (-3.93834e-09)*c*c*c*c ;
1873 }
1874 if(fRunIndex==33){
1875 v2c=0.0120077 + (-0.00139976)*c + (5.85677e-05)*c*c + (-9.06209e-07)*c*c*c + (4.61403e-09)*c*c*c*c ;
1876 v2s=-0.0197677 + (0.00253589)*c + (-0.000101374)*c*c + (1.59347e-06)*c*c*c + (-8.65961e-09)*c*c*c*c ;
1877 }
1878 if(fRunIndex==34){
1879 v2c=0.0759546 + (-0.0131538)*c + (0.000577001)*c*c + (-9.39636e-06)*c*c*c + (5.16814e-08)*c*c*c*c ;
1880 v2s=-0.0727864 + (0.00956655)*c + (-0.000406191)*c*c + (7.50453e-06)*c*c*c + (-4.93819e-08)*c*c*c*c ;
1881 }
1882 if(fRunIndex==35){
1883 v2c=-0.00622564 + (0.00107647)*c + (-6.25644e-05)*c*c + (1.28553e-06)*c*c*c + (-8.46286e-09)*c*c*c*c ;
1884 v2s=0.0051704 + (-0.00165348)*c + (8.69531e-05)*c*c + (-1.69895e-06)*c*c*c + (1.09762e-08)*c*c*c*c ;
1885 }
1886 if(fRunIndex==36){
1887 v2c=0.0117872 + (-0.000807391)*c + (-1.70999e-05)*c*c + (1.19089e-06)*c*c*c + (-1.13278e-08)*c*c*c*c ;
1888 v2s=-0.000162495 + (1.04465e-05)*c + (-1.07156e-05)*c*c + (2.54308e-07)*c*c*c + (-7.04948e-10)*c*c*c*c ;
1889 }
1890 if(fRunIndex==37){
1891 v2c=-0.0119236 + (0.00290105)*c + (-0.000162086)*c*c + (3.22645e-06)*c*c*c + (-2.07075e-08)*c*c*c*c ;
1892 v2s=-0.0178535 + (0.00283024)*c + (-0.00012812)*c*c + (2.02404e-06)*c*c*c + (-1.0248e-08)*c*c*c*c ;
1893 }
1894 if(fRunIndex==38){
1895 v2c=0.023002 + (-0.00433112)*c + (0.000229879)*c*c + (-4.52792e-06)*c*c*c + (2.92715e-08)*c*c*c*c ;
1896 v2s=-0.0217933 + (0.00235291)*c + (-0.000108054)*c*c + (2.08937e-06)*c*c*c + (-1.3524e-08)*c*c*c*c ;
1897 }
1898 if(fRunIndex==39){
1899 v2c=0.00795585 + (-0.000372683)*c + (1.85243e-05)*c*c + (-4.28133e-07)*c*c*c + (3.02452e-09)*c*c*c*c ;
1900 v2s=-0.0849826 + (0.00529466)*c + (-0.000163148)*c*c + (2.1186e-06)*c*c*c + (-9.48208e-09)*c*c*c*c ;
1901 }
1902 if(fRunIndex==40){
1903 v2c=0.00653019 + (-0.00127651)*c + (7.59264e-05)*c*c + (-1.6363e-06)*c*c*c + (1.11375e-08)*c*c*c*c ;
1904 v2s=-0.0107686 + (0.000982587)*c + (-2.97061e-05)*c*c + (4.04795e-07)*c*c*c + (-1.91289e-09)*c*c*c*c ;
1905 }
1906 if(fRunIndex==41){
1907 v2c=0.012007 + (-0.00102703)*c + (3.48714e-05)*c*c + (-5.13011e-07)*c*c*c + (2.57827e-09)*c*c*c*c ;
1908 v2s=-0.00713083 + (0.00101952)*c + (-4.7675e-05)*c*c + (9.12723e-07)*c*c*c + (-5.99694e-09)*c*c*c*c ;
1909 }
1910 if(fRunIndex==42){
1911 v2c=-0.0486696 + (0.00723345)*c + (-0.000324232)*c*c + (5.25317e-06)*c*c*c + (-2.78981e-08)*c*c*c*c ;
1912 v2s=-0.0485877 + (-0.00064506)*c + (9.15858e-05)*c*c + (-1.97486e-06)*c*c*c + (1.28831e-08)*c*c*c*c ;
1913 }
1914 if(fRunIndex==43){
1915 v2c=0.0003541 + (-0.000236832)*c + (1.82448e-05)*c*c + (-4.20764e-07)*c*c*c + (2.96882e-09)*c*c*c*c ;
1916 v2s=0.00107616 + (-0.000462058)*c + (2.46761e-05)*c*c + (-4.33199e-07)*c*c*c + (2.47324e-09)*c*c*c*c ;
1917 }
1918 if(fRunIndex==44){
1919 v2c=-0.0778259 + (0.0128797)*c + (-0.000514371)*c*c + (7.70095e-06)*c*c*c + (-3.89026e-08)*c*c*c*c ;
1920 v2s=-0.106839 + (0.00670369)*c + (-0.000184721)*c*c + (1.64581e-06)*c*c*c + (-1.47049e-09)*c*c*c*c ;
1921 }
1922 if(fRunIndex==45){
1923 v2c=0.0513049 + (-0.00715356)*c + (0.000289106)*c*c + (-4.46938e-06)*c*c*c + (2.36599e-08)*c*c*c*c ;
1924 v2s=-0.0252464 + (0.00170124)*c + (-3.04376e-05)*c*c + (9.4935e-08)*c*c*c + (9.41807e-10)*c*c*c*c ;
1925 }
1926 if(fRunIndex==46){
1927 v2c=0.0317697 + (-0.00121243)*c + (5.94543e-05)*c*c + (-1.26341e-06)*c*c*c + (8.49329e-09)*c*c*c*c ;
1928 v2s=-0.092276 + (0.00299941)*c + (-5.24713e-05)*c*c + (3.07331e-07)*c*c*c + (7.44139e-10)*c*c*c*c ;
1929 }
1930 if(fRunIndex==47){
1931 v2c=0.0429796 + (-0.000861588)*c + (6.82123e-07)*c*c + (1.96388e-07)*c*c*c + (-1.56984e-09)*c*c*c*c ;
1932 v2s=-0.133047 + (0.00853153)*c + (-0.000287983)*c*c + (4.1155e-06)*c*c*c + (-2.00239e-08)*c*c*c*c ;
1933 }
1934 if(fRunIndex==48){
1935 v2c=0.0480979 + (-0.000354969)*c + (-6.45176e-05)*c*c + (2.34807e-06)*c*c*c + (-1.9069e-08)*c*c*c*c ;
1936 v2s=-0.138315 + (0.0126874)*c + (-0.000672052)*c*c + (1.34449e-05)*c*c*c + (-8.64756e-08)*c*c*c*c ;
1937 }
1938 if(fRunIndex==51){
1939 v2c=-0.0114422 + (0.00506157)*c + (-0.000227041)*c*c + (3.97457e-06)*c*c*c + (-2.35322e-08)*c*c*c*c ;
1940 v2s=-0.208135 + (0.0160488)*c + (-0.000536132)*c*c + (7.2714e-06)*c*c*c + (-3.28328e-08)*c*c*c*c ;
1941 }
1942 if(fRunIndex==52){
1943 v2c=0.00113321 + (0.000922798)*c + (-6.5971e-05)*c*c + (1.40416e-06)*c*c*c + (-9.29905e-09)*c*c*c*c ;
1944 v2s=0.0550392 + (-0.00725661)*c + (0.000298984)*c*c + (-4.85473e-06)*c*c*c + (2.72063e-08)*c*c*c*c ;
1945 }
1946 if(fRunIndex==53){
1947 v2c=-0.003508 + (6.54451e-06)*c + (1.28603e-05)*c*c + (-3.83674e-07)*c*c*c + (3.1053e-09)*c*c*c*c ;
1948 v2s=0.00714978 + (-0.00198147)*c + (0.000113909)*c*c + (-2.1949e-06)*c*c*c + (1.34933e-08)*c*c*c*c ;
1949 }
1950 if(fRunIndex==54){
1951 v2c=-0.0024598 + (5.87979e-05)*c + (9.36516e-06)*c*c + (-2.66045e-07)*c*c*c + (1.80016e-09)*c*c*c*c ;
1952 v2s=-0.00651135 + (0.00109883)*c + (-4.41501e-05)*c*c + (6.19398e-07)*c*c*c + (-2.61977e-09)*c*c*c*c ;
1953 }
1954 if(fRunIndex==55){
1955 v2c=0.0472936 + (-0.00156718)*c + (3.03713e-05)*c*c + (-1.97139e-07)*c*c*c + (-3.55205e-10)*c*c*c*c ;
1956 v2s=-0.155014 + (0.00717506)*c + (-0.000168137)*c*c + (1.55419e-06)*c*c*c + (-3.64579e-09)*c*c*c*c ;
1957 }
1958 if(fRunIndex==56){
1959 v2c=0.011328 + (0.000924574)*c + (-6.803e-05)*c*c + (9.76956e-07)*c*c*c + (-2.90494e-09)*c*c*c*c ;
1960 v2s=-0.0412661 + (0.0107218)*c + (-0.000533108)*c*c + (9.75563e-06)*c*c*c + (-5.95608e-08)*c*c*c*c ;
1961 }
1962 if(fRunIndex==57){
1963 v2c=0.00900452 + (0.000588343)*c + (-3.02208e-05)*c*c + (2.13239e-07)*c*c*c + (9.25727e-10)*c*c*c*c ;
1964 v2s=0.13086 + (-0.0204768)*c + (0.000928956)*c*c + (-1.63917e-05)*c*c*c + (9.8273e-08)*c*c*c*c ;
1965 }
1966 if(fRunIndex==58){
1967 v2c=-0.0143951 + (0.00247373)*c + (-0.000107909)*c*c + (1.69603e-06)*c*c*c + (-9.05555e-09)*c*c*c*c ;
1968 v2s=-0.0209397 + (0.00180446)*c + (-5.97034e-05)*c*c + (8.17388e-07)*c*c*c + (-3.77989e-09)*c*c*c*c ;
1969 }
1970 if(fRunIndex==59){
1971 v2c=-0.0619279 + (0.00393575)*c + (-0.00011867)*c*c + (1.45001e-06)*c*c*c + (-5.98792e-09)*c*c*c*c ;
1972 v2s=-0.248103 + (0.0158553)*c + (-0.000529631)*c*c + (7.46158e-06)*c*c*c + (-3.60091e-08)*c*c*c*c ;
1973 }
1974 if(fRunIndex==60){
1975 v2c=-0.00886115 + (0.00152741)*c + (-6.94062e-05)*c*c + (1.11727e-06)*c*c*c + (-5.81351e-09)*c*c*c*c ;
1976 v2s=-0.000911417 + (0.000122344)*c + (-1.24605e-06)*c*c + (-6.29249e-08)*c*c*c + (8.81971e-10)*c*c*c*c ;
1977 }
1978 if(fRunIndex==61){
1979 v2c=-0.00510166 + (0.00040596)*c + (-4.37341e-06)*c*c + (-1.15302e-07)*c*c*c + (1.53428e-09)*c*c*c*c ;
1980 v2s=-0.0170123 + (0.00220818)*c + (-9.50905e-05)*c*c + (1.59566e-06)*c*c*c + (-8.92053e-09)*c*c*c*c ;
1981 }
1982 if(fRunIndex==62){
1983 v2c=-0.00332885 + (-0.000430498)*c + (5.60341e-05)*c*c + (-1.47922e-06)*c*c*c + (1.09723e-08)*c*c*c*c ;
1984 v2s=-0.0519534 + (0.00892818)*c + (-0.000434243)*c*c + (7.8407e-06)*c*c*c + (-4.69363e-08)*c*c*c*c ;
1985 }
1986 if(fRunIndex==63){
1987 v2c=-0.00176529 + (0.00123097)*c + (-6.59677e-05)*c*c + (1.24591e-06)*c*c*c + (-7.90285e-09)*c*c*c*c ;
1988 v2s=-0.0181792 + (0.00222004)*c + (-7.96902e-05)*c*c + (1.13047e-06)*c*c*c + (-5.72272e-09)*c*c*c*c ;
1989 }
1990 if(fRunIndex==64){
1991 v2c=0.00546184 + (-0.000449463)*c + (1.50833e-05)*c*c + (-1.72501e-07)*c*c*c + (3.58774e-10)*c*c*c*c ;
1992 v2s=0.013284 + (-0.00260671)*c + (0.000134403)*c*c + (-2.461e-06)*c*c*c + (1.47694e-08)*c*c*c*c ;
1993 }
1994 if(fRunIndex==65){
1995 v2c=0.0108447 + (-0.0017699)*c + (4.62716e-05)*c*c + (-3.4647e-07)*c*c*c + (5.24536e-10)*c*c*c*c ;
1996 v2s=0.0302247 + (-0.00605136)*c + (0.000290676)*c*c + (-4.87706e-06)*c*c*c + (2.65104e-08)*c*c*c*c ;
1997 }
1998 if(fRunIndex==66){
1999 v2c=-0.0100092 + (0.00144457)*c + (-6.512e-05)*c*c + (1.05356e-06)*c*c*c + (-5.52231e-09)*c*c*c*c ;
2000 v2s=-0.00380228 + (0.000610003)*c + (-3.47943e-05)*c*c + (7.21188e-07)*c*c*c + (-4.76727e-09)*c*c*c*c ;
2001 }
2002 if(fRunIndex==67){
2003 v2c=-0.0229901 + (0.00400737)*c + (-0.000206024)*c*c + (3.93038e-06)*c*c*c + (-2.46388e-08)*c*c*c*c ;
2004 v2s=-0.00286041 + (0.000648607)*c + (-2.37397e-05)*c*c + (2.0372e-07)*c*c*c + (2.97679e-10)*c*c*c*c ;
2005 }
2006 if(fRunIndex==68){
2007 v2c=-0.0603647 + (0.0045685)*c + (-0.000171555)*c*c + (2.66377e-06)*c*c*c + (-1.45892e-08)*c*c*c*c ;
2008 v2s=-0.21447 + (0.00909524)*c + (-0.000265298)*c*c + (3.73262e-06)*c*c*c + (-1.84484e-08)*c*c*c*c ;
2009 }
2010 if(fRunIndex==69){
2011 v2c=-0.0207879 + (0.000139613)*c + (1.27802e-05)*c*c + (-6.40328e-07)*c*c*c + (6.35949e-09)*c*c*c*c ;
2012 v2s=-0.266459 + (0.0156802)*c + (-0.000481647)*c*c + (6.4597e-06)*c*c*c + (-3.047e-08)*c*c*c*c ;
2013 }
2014 if(fRunIndex==70){
2015 v2c=-0.0324843 + (0.00376852)*c + (-0.000136873)*c*c + (2.01389e-06)*c*c*c + (-1.03877e-08)*c*c*c*c ;
2016 v2s=-0.0178378 + (0.00145491)*c + (-3.55421e-05)*c*c + (2.25993e-07)*c*c*c + (7.06258e-10)*c*c*c*c ;
2017 }
2018 if(fRunIndex==71){
2019 v2c=-0.0138645 + (0.00167315)*c + (-6.68544e-05)*c*c + (1.07133e-06)*c*c*c + (-5.92659e-09)*c*c*c*c ;
2020 v2s=0.0359716 + (-0.00461076)*c + (0.000192868)*c*c + (-3.21409e-06)*c*c*c + (1.84468e-08)*c*c*c*c ;
2021 }
2022 if(fRunIndex==72){
2023 v2c=0.0111712 + (-0.00175546)*c + (6.87862e-05)*c*c + (-8.99029e-07)*c*c*c + (3.44492e-09)*c*c*c*c ;
2024 v2s=-0.00352823 + (0.000480673)*c + (-2.66091e-05)*c*c + (5.54922e-07)*c*c*c + (-3.84534e-09)*c*c*c*c ;
2025 }
2026 if(fRunIndex==73){
2027 v2c=-0.000939145 + (-0.000235945)*c + (2.51908e-05)*c*c + (-6.52435e-07)*c*c*c + (4.8218e-09)*c*c*c*c ;
2028 v2s=0.00666852 + (-0.000373105)*c + (-1.22838e-05)*c*c + (7.3621e-07)*c*c*c + (-7.2867e-09)*c*c*c*c ;
2029 }
2030 if(fRunIndex==74){
2031 v2c=0.0474235 + (-0.00890967)*c + (0.000427741)*c*c + (-7.38723e-06)*c*c*c + (4.17242e-08)*c*c*c*c ;
2032 v2s=0.0142517 + (-0.00326338)*c + (0.000142729)*c*c + (-2.27831e-06)*c*c*c + (1.14084e-08)*c*c*c*c ;
2033 }
2034 if(fRunIndex==75){
2035 v2c=0.0237808 + (-0.00310491)*c + (0.000134288)*c*c + (-2.27429e-06)*c*c*c + (1.32749e-08)*c*c*c*c ;
2036 v2s=-0.0258905 + (0.00417719)*c + (-0.000193901)*c*c + (3.41383e-06)*c*c*c + (-2.01695e-08)*c*c*c*c ;
2037 }
2038 if(fRunIndex==76){
2039 v2c=0.0149865 + (-0.00140739)*c + (4.24537e-05)*c*c + (-4.06273e-07)*c*c*c + (5.44496e-10)*c*c*c*c ;
2040 v2s=0.00820134 + (-0.000519432)*c + (6.82211e-06)*c*c + (1.06291e-07)*c*c*c + (-1.65344e-09)*c*c*c*c ;
2041 }
2042 if(fRunIndex==77){
2043 v2c=0.048123 + (-0.00608197)*c + (0.000246409)*c*c + (-4.05028e-06)*c*c*c + (2.33536e-08)*c*c*c*c ;
2044 v2s=-0.0214729 + (0.0029389)*c + (-0.000133826)*c*c + (2.29845e-06)*c*c*c + (-1.29253e-08)*c*c*c*c ;
2045 }
2046 if(fRunIndex==78){
2047 v2c=-0.0172924 + (0.00202216)*c + (-6.478e-05)*c*c + (8.08739e-07)*c*c*c + (-3.73211e-09)*c*c*c*c ;
2048 v2s=-0.0355638 + (0.00121518)*c + (3.97839e-05)*c*c + (-1.62346e-06)*c*c*c + (1.27105e-08)*c*c*c*c ;
2049 }
2050 if(fRunIndex==79){
2051 v2c=0.040762 + (-0.00311921)*c + (6.66868e-05)*c*c + (-1.18554e-06)*c*c*c + (1.07519e-08)*c*c*c*c ;
2052 v2s=-0.00558438 + (-0.00154821)*c + (0.000165732)*c*c + (-4.56377e-06)*c*c*c + (3.56988e-08)*c*c*c*c ;
2053 }
2054 if(fRunIndex==81){
2055 v2c=-0.0247995 + (0.00179444)*c + (-0.000102085)*c*c + (2.31093e-06)*c*c*c + (-1.61924e-08)*c*c*c*c ;
2056 v2s=-0.249009 + (0.0133297)*c + (-0.000415717)*c*c + (5.84494e-06)*c*c*c + (-2.88179e-08)*c*c*c*c ;
2057 }
2058 if(fRunIndex==82){
2059 v2c=-0.0249921 + (0.00199208)*c + (-7.86121e-05)*c*c + (1.40547e-06)*c*c*c + (-9.39651e-09)*c*c*c*c ;
2060 v2s=-0.254469 + (0.0150694)*c + (-0.000499848)*c*c + (7.08318e-06)*c*c*c + (-3.40806e-08)*c*c*c*c ;
2061 }
2062 if(fRunIndex==83){
2063 v2c=-0.0721 + (0.00951559)*c + (-0.000431762)*c*c + (8.4632e-06)*c*c*c + (-5.79485e-08)*c*c*c*c ;
2064 v2s=-0.00864625 + (-0.00706196)*c + (0.000313773)*c*c + (-5.10066e-06)*c*c*c + (2.9662e-08)*c*c*c*c ;
2065 }
2066 if(fRunIndex==84){
2067 v2c=-0.0198112 + (-0.0062788)*c + (0.000353589)*c*c + (-5.70525e-06)*c*c*c + (2.82617e-08)*c*c*c*c ;
2068 v2s=-0.0704647 + (0.0031015)*c + (-0.000244985)*c*c + (6.17258e-06)*c*c*c + (-4.43331e-08)*c*c*c*c ;
2069 }
2070 if(fRunIndex==85){
2071 v2c=-0.0161955 + (-0.0017361)*c + (0.000130612)*c*c + (-2.9548e-06)*c*c*c + (2.05684e-08)*c*c*c*c ;
2072 v2s=-0.204119 + (0.00714989)*c + (-0.000162549)*c*c + (1.87459e-06)*c*c*c + (-7.79143e-09)*c*c*c*c ;
2073 }
2074 if(fRunIndex==86){
2075 v2c=0.00995155 + (-0.00196822)*c + (0.00010617)*c*c + (-2.01071e-06)*c*c*c + (1.22728e-08)*c*c*c*c ;
2076 v2s=0.0148901 + (-0.002581)*c + (0.000132002)*c*c + (-2.44601e-06)*c*c*c + (1.48144e-08)*c*c*c*c ;
2077 }
2078 if(fRunIndex==87){
2079 v2c=-0.00411031 + (0.000510639)*c + (-2.27814e-05)*c*c + (6.41192e-07)*c*c*c + (-5.85467e-09)*c*c*c*c ;
2080 v2s=0.0121118 + (-0.00281)*c + (0.000143754)*c*c + (-2.62557e-06)*c*c*c + (1.61207e-08)*c*c*c*c ;
2081 }
2082 if(fRunIndex==88){
2083 v2c=-0.00570308 + (-0.000945116)*c + (7.42606e-05)*c*c + (-1.51212e-06)*c*c*c + (9.57993e-09)*c*c*c*c ;
2084 v2s=-0.0282666 + (0.00398475)*c + (-0.000165678)*c*c + (2.68542e-06)*c*c*c + (-1.4688e-08)*c*c*c*c ;
2085 }
2086 if(fRunIndex==89){
2087 v2c=0.0845785 + (-0.012611)*c + (0.000561197)*c*c + (-9.39333e-06)*c*c*c + (5.16139e-08)*c*c*c*c ;
2088 v2s=0.149689 + (-0.0193648)*c + (0.000771629)*c*c + (-1.21901e-05)*c*c*c + (6.6443e-08)*c*c*c*c ;
2089 }
2090 if(fRunIndex==91){
2091 v2c=-0.0289605 + (0.00127922)*c + (-4.47445e-06)*c*c + (-4.00293e-07)*c*c*c + (3.82966e-09)*c*c*c*c ;
2092 v2s=-0.241098 + (0.00906132)*c + (-0.000159552)*c*c + (6.47243e-07)*c*c*c + (4.61048e-09)*c*c*c*c ;
2093 }
2094 if(fRunIndex==95){
2095 v2c=0.0144072 + (-0.000966056)*c + (6.47254e-05)*c*c + (-1.84164e-06)*c*c*c + (1.51588e-08)*c*c*c*c ;
2096 v2s=-0.306678 + (0.0184746)*c + (-0.000540218)*c*c + (6.36681e-06)*c*c*c + (-2.37512e-08)*c*c*c*c ;
2097 }
2098 if(fRunIndex==96){
2099 v2c=0.00506182 + (-0.000632672)*c + (1.89752e-05)*c*c + (-2.14999e-07)*c*c*c + (3.12179e-10)*c*c*c*c ;
2100 v2s=-0.218693 + (0.00676786)*c + (-0.000115485)*c*c + (5.56737e-07)*c*c*c + (2.68891e-09)*c*c*c*c ;
2101 }
2102 if(fRunIndex==97){
2103 v2c=-0.0225062 + (0.00204445)*c + (-6.76028e-05)*c*c + (8.92976e-07)*c*c*c + (-4.00576e-09)*c*c*c*c ;
2104 v2s=-0.260286 + (0.0139325)*c + (-0.000425631)*c*c + (5.64038e-06)*c*c*c + (-2.51554e-08)*c*c*c*c ;
2105 }
2106 if(fRunIndex==98){
2107 v2c=-0.0261849 + (0.0063862)*c + (-0.000293076)*c*c + (4.37234e-06)*c*c*c + (-2.04899e-08)*c*c*c*c ;
2108 v2s=-0.17744 + (0.00236157)*c + (9.2782e-06)*c*c + (-5.63126e-07)*c*c*c + (4.9714e-09)*c*c*c*c ;
2109 }
2110 if(fRunIndex==101){
2111 v2c=0.101445 + (-0.015664)*c + (0.000747255)*c*c + (-1.36979e-05)*c*c*c + (8.36903e-08)*c*c*c*c ;
2112 v2s=0.0152975 + (-0.00329946)*c + (0.00018627)*c*c + (-3.88572e-06)*c*c*c + (2.64895e-08)*c*c*c*c ;
2113 }
2114 if(fRunIndex==102){
2115 v2c=-0.0344927 + (0.00470964)*c + (-0.000188941)*c*c + (3.0206e-06)*c*c*c + (-1.69735e-08)*c*c*c*c ;
2116 v2s=0.0142624 + (-0.00240623)*c + (9.16657e-05)*c*c + (-1.14663e-06)*c*c*c + (4.14195e-09)*c*c*c*c ;
2117 }
2118 if(fRunIndex==103){
2119 v2c=-0.0238316 + (0.00322675)*c + (-0.000123462)*c*c + (1.76428e-06)*c*c*c + (-8.66176e-09)*c*c*c*c ;
2120 v2s=0.0227184 + (-0.00299954)*c + (0.00014422)*c*c + (-2.81018e-06)*c*c*c + (1.82546e-08)*c*c*c*c ;
2121 }
2122 if(fRunIndex==104){
2123 v2c=0.0118014 + (-0.000998833)*c + (3.77271e-05)*c*c + (-6.19183e-07)*c*c*c + (3.49065e-09)*c*c*c*c ;
2124 v2s=0.000977264 + (0.00012542)*c + (1.40219e-08)*c*c + (-2.24472e-09)*c*c*c + (-3.39763e-10)*c*c*c*c ;
2125 }
2126 if(fRunIndex==105){
2127 v2c=0.0312332 + (-0.00348485)*c + (0.000108654)*c*c + (-1.28443e-06)*c*c*c + (4.80188e-09)*c*c*c*c ;
2128 v2s=-0.0344234 + (0.00443313)*c + (-0.000173309)*c*c + (2.7124e-06)*c*c*c + (-1.42518e-08)*c*c*c*c ;
2129 }
2130 if(fRunIndex==106){
2131 v2c=0.0369386 + (-0.00434063)*c + (0.00018894)*c*c + (-3.45688e-06)*c*c*c + (2.16617e-08)*c*c*c*c ;
2132 v2s=-0.00675251 + (0.000691284)*c + (-1.08033e-06)*c*c + (-5.57352e-07)*c*c*c + (6.28148e-09)*c*c*c*c ;
2133 }
2134 if(fRunIndex==107){
2135 v2c=-0.0441483 + (0.00803486)*c + (-0.00034871)*c*c + (4.76349e-06)*c*c*c + (-1.81393e-08)*c*c*c*c ;
2136 v2s=-0.141046 + (0.0165183)*c + (-0.000592546)*c*c + (8.956e-06)*c*c*c + (-4.86662e-08)*c*c*c*c ;
2137 }
2138 if(fRunIndex==108){
2139 v2c=-0.0531169 + (0.00637348)*c + (-0.000218334)*c*c + (2.66873e-06)*c*c*c + (-9.86859e-09)*c*c*c*c ;
2140 v2s=-0.00809622 + (-0.000552159)*c + (6.27984e-05)*c*c + (-1.33136e-06)*c*c*c + (7.45069e-09)*c*c*c*c ;
2141 }
2142 if(fRunIndex==109){
2143 v2c=-0.00189593 + (0.00153994)*c + (-9.92105e-05)*c*c + (1.98756e-06)*c*c*c + (-1.24677e-08)*c*c*c*c ;
2144 v2s=0.0252107 + (-0.00330651)*c + (0.000128362)*c*c + (-1.89951e-06)*c*c*c + (9.45477e-09)*c*c*c*c ;
2145 }
2146 if(fRunIndex==112){
2147 v2c=-0.00939028 + (0.00128226)*c + (-5.24125e-05)*c*c + (8.50382e-07)*c*c*c + (-4.65785e-09)*c*c*c*c ;
2148 v2s=0.00201633 + (-0.00116151)*c + (7.59727e-05)*c*c + (-1.59962e-06)*c*c*c + (1.04645e-08)*c*c*c*c ;
2149 }
2150 if(fRunIndex==114){
2151 v2c=-0.00783351 + (0.00122277)*c + (-3.37736e-05)*c*c + (3.12136e-07)*c*c*c + (-9.19359e-10)*c*c*c*c ;
2152 v2s=0.0103833 + (-0.00215968)*c + (9.96942e-05)*c*c + (-1.71974e-06)*c*c*c + (1.01511e-08)*c*c*c*c ;
2153 }
2154 if(fRunIndex==115){
2155 v2c=0.00588164 + (-0.00083683)*c + (5.13684e-05)*c*c + (-1.14177e-06)*c*c*c + (8.11899e-09)*c*c*c*c ;
2156 v2s=0.0037561 + (-0.000786163)*c + (3.88084e-05)*c*c + (-7.73678e-07)*c*c*c + (5.27396e-09)*c*c*c*c ;
2157 }
2158 if(fRunIndex==116){
2159 v2c=0.0610685 + (-0.00392239)*c + (6.23272e-05)*c*c + (9.84282e-07)*c*c*c + (-1.75418e-08)*c*c*c*c ;
2160 v2s=-0.106295 + (0.0124333)*c + (-0.000510299)*c*c + (8.51235e-06)*c*c*c + (-4.9787e-08)*c*c*c*c ;
2161 }
2162 if(fRunIndex==117){
2163 v2c=-0.0132264 + (0.0024143)*c + (-0.000106307)*c*c + (1.74632e-06)*c*c*c + (-9.79752e-09)*c*c*c*c ;
2164 v2s=-0.00788943 + (0.000443907)*c + (-1.83984e-05)*c*c + (3.81018e-07)*c*c*c + (-2.52604e-09)*c*c*c*c ;
2165 }
2166 if(fRunIndex==118){
2167 v2c=-0.045835 + (0.00849686)*c + (-0.000408675)*c*c + (7.3031e-06)*c*c*c + (-4.34655e-08)*c*c*c*c ;
2168 v2s=-0.0455199 + (0.00718851)*c + (-0.000340502)*c*c + (6.00457e-06)*c*c*c + (-3.51883e-08)*c*c*c*c ;
2169 }
2170 if(fRunIndex==119){
2171 v2c=-0.00627236 + (0.00146712)*c + (-9.05375e-05)*c*c + (1.96962e-06)*c*c*c + (-1.35867e-08)*c*c*c*c ;
2172 v2s=-0.00544096 + (0.000598445)*c + (-7.61032e-06)*c*c + (-2.09579e-07)*c*c*c + (3.00091e-09)*c*c*c*c ;
2173 }
2174 if(fRunIndex==120){
2175 v2c=-0.0376823 + (0.00376965)*c + (-8.01573e-05)*c*c + (-2.80792e-09)*c*c*c + (6.79508e-09)*c*c*c*c ;
2176 v2s=-0.0675803 + (0.00849852)*c + (-0.000298266)*c*c + (3.93178e-06)*c*c*c + (-1.68582e-08)*c*c*c*c ;
2177 }
2178 if(fRunIndex==121){
2179 v2c=0.0285153 + (-0.00340253)*c + (0.000130854)*c*c + (-2.02449e-06)*c*c*c + (1.09594e-08)*c*c*c*c ;
2180 v2s=0.0224991 + (-0.00524188)*c + (0.000269226)*c*c + (-4.93926e-06)*c*c*c + (2.97604e-08)*c*c*c*c ;
2181 }
2182 if(fRunIndex==122){
2183 v2c=0.00379321 + (0.000532446)*c + (-6.16649e-05)*c*c + (1.67786e-06)*c*c*c + (-1.29277e-08)*c*c*c*c ;
2184 v2s=-0.011158 + (0.00256124)*c + (-0.000150921)*c*c + (3.22488e-06)*c*c*c + (-2.16809e-08)*c*c*c*c ;
2185 }
2186 if(fRunIndex==123){
2187 v2c=-0.0607006 + (0.00831906)*c + (-0.000398705)*c*c + (6.66385e-06)*c*c*c + (-3.5222e-08)*c*c*c*c ;
2188 v2s=0.0389204 + (-0.00798361)*c + (0.000297135)*c*c + (-3.79751e-06)*c*c*c + (1.47914e-08)*c*c*c*c ;
2189 }
2190 if(fRunIndex==124){
2191 v2c=-0.0288952 + (0.00325728)*c + (-0.000104126)*c*c + (1.29062e-06)*c*c*c + (-6.03539e-09)*c*c*c*c ;
2192 v2s=0.0318884 + (-0.00301585)*c + (0.000131913)*c*c + (-2.5702e-06)*c*c*c + (1.71313e-08)*c*c*c*c ;
2193 }
2194 if(fRunIndex==125){
2195 v2c=-0.00413859 + (0.00024613)*c + (5.33902e-06)*c*c + (-3.86963e-07)*c*c*c + (3.89277e-09)*c*c*c*c ;
2196 v2s=-0.0492298 + (0.00624793)*c + (-0.000243844)*c*c + (3.69001e-06)*c*c*c + (-1.9156e-08)*c*c*c*c ;
2197 }
2198 if(fRunIndex==126){
2199 v2c=-0.018816 + (0.00134076)*c + (-4.50486e-05)*c*c + (5.84259e-07)*c*c*c + (-2.34235e-09)*c*c*c*c ;
2200 v2s=-0.295912 + (0.0159608)*c + (-0.000497569)*c*c + (6.92442e-06)*c*c*c + (-3.35159e-08)*c*c*c*c ;
2201 }
2202 if(fRunIndex==127){
2203 v2c=-0.039222 + (0.00400933)*c + (-0.0001341)*c*c + (1.65214e-06)*c*c*c + (-6.33829e-09)*c*c*c*c ;
2204 v2s=-0.290884 + (0.0147027)*c + (-0.000429021)*c*c + (5.64116e-06)*c*c*c + (-2.58784e-08)*c*c*c*c ;
2205 }
2206 if(fRunIndex==128){
2207 v2c=0.0168685 + (-0.00285775)*c + (0.000141099)*c*c + (-2.48846e-06)*c*c*c + (1.41845e-08)*c*c*c*c ;
2208 v2s=-0.278776 + (0.0113771)*c + (-0.000279052)*c*c + (3.16144e-06)*c*c*c + (-1.18545e-08)*c*c*c*c ;
2209 }
2210 if(fRunIndex==129){
2211 v2c=0.0490322 + (-0.00736827)*c + (0.000304529)*c*c + (-4.89472e-06)*c*c*c + (2.60746e-08)*c*c*c*c ;
2212 v2s=-0.0961198 + (0.0179264)*c + (-0.000914589)*c*c + (1.70145e-05)*c*c*c + (-1.03423e-07)*c*c*c*c ;
2213 }
2214 if(fRunIndex==131){
2215 v2c=0.075924 + (-0.0108106)*c + (0.000493471)*c*c + (-8.39074e-06)*c*c*c + (4.67109e-08)*c*c*c*c ;
2216 v2s=0.0700446 + (-0.0105027)*c + (0.00048267)*c*c + (-8.4422e-06)*c*c*c + (4.91292e-08)*c*c*c*c ;
2217 }
2218 if(fRunIndex==132){
2219 v2c=0.00752887 + (-0.00142672)*c + (7.37857e-05)*c*c + (-1.36894e-06)*c*c*c + (8.28054e-09)*c*c*c*c ;
2220 v2s=-0.0121719 + (0.00188256)*c + (-8.67104e-05)*c*c + (1.50514e-06)*c*c*c + (-8.79759e-09)*c*c*c*c ;
2221 }
2222 if(fRunIndex==133){
2223 v2c=0.00242068 + (-0.000274576)*c + (1.06631e-05)*c*c + (-1.45829e-07)*c*c*c + (7.18001e-10)*c*c*c*c ;
2224 v2s=0.0113182 + (-0.000529549)*c + (-7.75839e-06)*c*c + (4.86382e-07)*c*c*c + (-3.98634e-09)*c*c*c*c ;
2225 }
2226 if(fRunIndex==134){
2227 v2c=-0.0438721 + (0.00906564)*c + (-0.000407471)*c*c + (6.6254e-06)*c*c*c + (-3.61035e-08)*c*c*c*c ;
2228 v2s=-0.276264 + (0.0107322)*c + (-0.000281971)*c*c + (3.27634e-06)*c*c*c + (-1.17426e-08)*c*c*c*c ;
2229 }
2230 if(fRunIndex==135){
2231 v2c=0.0224469 + (-0.00324762)*c + (0.000165861)*c*c + (-3.16496e-06)*c*c*c + (2.011e-08)*c*c*c*c ;
2232 v2s=-0.308891 + (0.0136721)*c + (-0.00032343)*c*c + (2.93575e-06)*c*c*c + (-5.7637e-09)*c*c*c*c ;
2233 }
2234 if(fRunIndex==136){
2235 v2c=0.0279685 + (-0.00441)*c + (0.000294342)*c*c + (-6.56635e-06)*c*c*c + (4.53669e-08)*c*c*c*c ;
2236 v2s=-0.273287 + (0.00917778)*c + (-0.000144155)*c*c + (1.59881e-07)*c*c*c + (8.63874e-09)*c*c*c*c ;
2237 }
2238
2239 /*
2240 if(fRunIndex>=1 && fRunIndex<=45){ // 137161 -137848 //period 1
2241 // v2c =-3.523975e-03+8.132703e-04*x+-4.107765e-05*x*x+7.308763e-07*x*x*x+-4.285144e-09*x*x*x*x;
2242 // v2s =-1.555676e-02+1.075371e-03*x+-3.597818e-05*x*x+4.927188e-07*x*x*x+-2.217080e-09*x*x*x*x;
2243 v2c=-2.007559e-05+3.143290e-04*x-2.065995e-05*x*x+4.071671e-07*x*x*x-2.562114e-09*x*x*x*x ;
2244 v2s=-4.920033e-03+6.807606e-04*x-3.065278e-05*x*x+5.132261e-07*x*x*x-2.767644e-09*x*x*x*x ;
2245
2246 }
2247 if(fRunIndex>=46 && fRunIndex<=58){ //138125-138275//period 2
2248 v2c = 1.591471e-02+2.772445e-04*x-2.229627e-05*x*x+3.803341e-07*x*x*x-2.056551e-09*x*x*x*x;
2249 v2s =-5.800393e-02+2.007102e-03*x-3.292094e-05*x*x+7.266378e-08*x*x*x+1.828472e-09*x*x*x*x ;
2250 }
2251 if(fRunIndex>=59 && fRunIndex<=78){ //138359-138730//period 3
2252 v2c =-8.010228e-03+8.947184e-04*x-4.192609e-05*x*x+7.388210e-07*x*x*x-4.280541e-09*x*x*x*x ;
2253 v2s =-3.203888e-02+2.122650e-03*x-7.549389e-05*x*x+1.161380e-06*x*x*x-6.156621e-09*x*x*x*x ;
2254 }
2255 if(fRunIndex>=79){//period 4
2256 v2c =-1.025964e-02+1.128454e-03*x-4.603917e-05*x*x+7.386637e-07*x*x*x-4.035342e-09*x*x*x*x ;
2257 v2s =-9.746961e-02+4.672675e-03*x-1.311440e-04*x*x+1.650922e-06*x*x*x-7.083187e-09*x*x*x*x ;
2258 }
2259*/
2260 return phi+v2c*TMath::Sin(2.*phi)-v2s*TMath::Cos(2.*phi) ;
2261}
2262//____________________________________________________________________________
2263Double_t AliAnalysisTaskPi0v2::ApplyFlatteningTPC(Double_t phi, Double_t c){
2264
2265
2266 if(fRunIndex==39){ //137748
2267 Double_t v2c= 0.0222974-0.00132297*c+0.000204021*c*c-4.49827e-06*c*c*c+2.76986e-08*c*c*c*c ;
2268 Double_t v2s=0.00615401-7.76462e-05*c ;
2269 return phi+v2c*TMath::Sin(2.*phi)-v2s*TMath::Cos(2.*phi) ;
2270 }
2271
2272
2273 //Periods 1,2,3
2274 //fRunIndex - run index
2275 if(fRunIndex>=1 && fRunIndex<79){
2276 Double_t v2c=4.40516e-04+TMath::Exp(-4.71923-7.62089e-02*c) ;
2277 Double_t v2s=1.79859e-03+TMath::Exp(-4.99649-7.87523e-02*c) ;
2278
2279 return phi+v2c*TMath::Sin(2.*phi)-v2s*TMath::Cos(2.*phi) ;
2280 }
2281
2282 //period4
2283 //138826, 138828, 138830
2284 if(fRunIndex==83 ||fRunIndex==84 ||fRunIndex==85 ){
2285 Double_t v2c = -1.4e-03 ;
2286 Double_t v2s = -5.60117e-02-TMath::Exp( -1.62827e+00-9.96071e-02*c);
2287 return phi + v2c*TMath::Sin(2.*phi) - v2s*TMath::Cos(2.*phi) ;
2288 }
2289 //Runs 138871,138872
2290 if(fRunIndex==87 ||fRunIndex==88 ){
2291 Double_t v2c=-0.00518371 ;
2292 Double_t v2s= 0.00633324 ;
2293 return phi + v2c*TMath::Sin(2.*phi) - v2s*TMath::Cos(2.*phi) ;
2294 }
2295 //run 139029
2296 if(fRunIndex==102 ){
2297 Double_t v2c = -0.00354633 ;
2298 Double_t v2s = 0.00418512 ;
2299 return phi+v2c*TMath::Sin(2.*phi)-v2s*TMath::Cos(2.*phi) ;
2300 }
2301 //run 139110
2302/*
2303 if(fRunIndex==110 ){
2304 Double_t v2c=4.78327e-04+TMath::Exp(-3.21625 -6.74306e-02*c) ;
2305 Double_t v2s=3.50176e-02+TMath::Exp(-9.73541e-01-9.19214e-02*c) ;
2306 return phi + v2c*TMath::Sin(2.*phi) - v2s*TMath::Cos(2.*phi) ;
2307 }
2308*/
2309 if(fRunIndex>=79 && fRunIndex<=139){
2310 Double_t v2c= 4.78327e-04 + TMath::Exp(-6.70587-2.33120e-02*c) ;
2311 Double_t v2s=-2.57731e-03 + TMath::Exp(-2.75493-1.05166e-01*c) ;
2312
2313 return phi + v2c*TMath::Sin(2.*phi) - v2s*TMath::Cos(2.*phi) ;
2314 }
2315 return phi ;
2316
2317}