const Int_t numberOfCentralityBins = 14; TString centralityArray[numberOfCentralityBins] = {"0-80","10-20","20-30","30-40","40-50","50-60","60-70","70-80","0-100","0-1","1-2","2-3","92-8500","0-10000"}; const Int_t gRebin = 1; void drawCorrelationFunctionPsi(const char* filename = "AnalysisResultsPsi_train_06feb.root", Int_t gCentrality = 1, Int_t gBit = -1, const char* gCentralityEstimator = 0x0, Bool_t kShowShuffled = kFALSE, Bool_t kShowMixed = kTRUE, Double_t psiMin = -0.5, Double_t psiMax = 3.5, Double_t vertexZMin = -10., Double_t vertexZMax = 10., Double_t ptTriggerMin = -1., Double_t ptTriggerMax = -1., Double_t ptAssociatedMin = -1., Double_t ptAssociatedMax = -1., Bool_t normToTrig = kTRUE, Double_t normalizationRangePhi = TMath::Pi()/6., Bool_t kUseVzBinning = kTRUE, Int_t rebinEta = 1, Int_t rebinPhi = 1, TString eventClass = "EventPlane", //Can be "EventPlane", "Centrality", "Multiplicity" Bool_t bToy = kFALSE, TString listNameAdd = "") { //Macro that draws the correlation functions from the balance function //analysis vs the reaction plane //Author: Panos.Christakoglou@nikhef.nl //gROOT->LoadMacro("~/SetPlotStyle.C"); //SetPlotStyle(); gStyle->SetPalette(1,0); //Load the PWG2ebye library gSystem->Load("libCore.so"); gSystem->Load("libGeom.so"); gSystem->Load("libVMC.so"); gSystem->Load("libPhysics.so"); gSystem->Load("libTree.so"); gSystem->Load("libSTEERBase.so"); gSystem->Load("libESD.so"); gSystem->Load("libAOD.so"); gSystem->Load("libANALYSIS.so"); gSystem->Load("libANALYSISalice.so"); gSystem->Load("libEventMixing.so"); gSystem->Load("libCORRFW.so"); gSystem->Load("libPWGTools.so"); gSystem->Load("libPWGCFebye.so"); //Prepare the objects and return them TList *listQA = NULL; TList *list = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,0,bToy,listNameAdd); TList *listShuffled = NULL; if(kShowShuffled) listShuffled = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,1,bToy,listNameAdd); TList *listMixed = NULL; if(kShowMixed) listMixed = GetListOfObjects(filename,gCentrality,gBit,gCentralityEstimator,2,bToy,listNameAdd); //get the QA histograms (for convolution) //Open the file again TFile *f = TFile::Open(filename,"UPDATE"); if((!f)||(!f->IsOpen())) { Printf("The file %s is not found.",filename); } TDirectoryFile *dir = dynamic_cast<TDirectoryFile *>(f->Get("PWGCFEbyE.outputBalanceFunctionPsiAnalysis")); if(!dir) { Printf("The TDirectoryFile is not found.",filename); } TString listQAName = "listQAPsi_"; listQAName += centralityArray[gCentrality-1]; if(gBit > -1) { listQAName += "_Bit"; listQAName += gBit; } if(gCentralityEstimator) { listQAName += "_"; listQAName += gCentralityEstimator;} listQAName += listNameAdd; listQA = (TList*)dir->Get(Form("%s",listQAName.Data())); if(!listQA) { Printf("TList %s not found!!!",listQAName.Data()); } if(!list) { Printf("The TList object was not created"); return; } else draw(list,listShuffled,listMixed,listQA, gCentralityEstimator,gCentrality,psiMin,psiMax,vertexZMin,vertexZMax, ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,normToTrig, normalizationRangePhi, kUseVzBinning,rebinEta,rebinPhi,eventClass,bToy); } //______________________________________________________// TList *GetListOfObjects(const char* filename, Int_t gCentrality, Int_t gBit, const char *gCentralityEstimator, Int_t kData = 1, Bool_t bToy = kFALSE, TString listNameAdd = "") { //Get the TList objects (QA, bf, bf shuffled) TList *listBF = 0x0; //Open the file TFile *f = TFile::Open(filename,"UPDATE"); if((!f)||(!f->IsOpen())) { Printf("The file %s is not found. Aborting...",filename); return listBF; } //f->ls(); TDirectoryFile *dir = dynamic_cast<TDirectoryFile *>(f->Get("PWGCFEbyE.outputBalanceFunctionPsiAnalysis")); if(!dir) { Printf("The TDirectoryFile is not found. Aborting...",filename); return listBF; } //dir->ls(); TString listBFName; if(kData == 0) { //cout<<"no shuffling - no mixing"<<endl; listBFName = "listBFPsi"; } else if(kData == 1) { //cout<<"shuffling - no mixing"<<endl; listBFName = "listBFPsiShuffled"; } else if(kData == 2) { //cout<<"no shuffling - mixing"<<endl; listBFName = "listBFPsiMixed"; } // different list names in case of toy model if(!bToy){ listBFName += "_"; listBFName += centralityArray[gCentrality-1]; if(gBit > -1) { listBFName += "_Bit"; listBFName += gBit; } if(gCentralityEstimator) { listBFName += "_"; listBFName += gCentralityEstimator;} } else{ listBFName.ReplaceAll("Psi",""); } listBFName += listNameAdd; // histograms were already retrieved (in first iteration) if(dir->Get(Form("%s_histograms",listBFName.Data()))){ listBF = dynamic_cast<TList *>(dir->Get(Form("%s_histograms",listBFName.Data()))); } // histograms were not yet retrieved (this is the first iteration) else{ listBF = dynamic_cast<TList *>(dir->Get(listBFName.Data())); if(!listBF) dir->ls(); cout<<"======================================================="<<endl; cout<<"List name (control): "<<listBFName.Data()<<endl; cout<<"List name: "<<listBF->GetName()<<endl; //listBF->ls(); //Get the histograms TString histoName; if(kData == 0) histoName = "fHistP"; else if(kData == 1) histoName = "fHistP_shuffle"; else if(kData == 2) histoName = "fHistP"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistP = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistP) { Printf("fHistP %s not found!!!",histoName.Data()); break; } fHistP->FillParent(); fHistP->DeleteContainers(); if(kData == 0) histoName = "fHistN"; if(kData == 1) histoName = "fHistN_shuffle"; if(kData == 2) histoName = "fHistN"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistN = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistN) { Printf("fHistN %s not found!!!",histoName.Data()); break; } fHistN->FillParent(); fHistN->DeleteContainers(); if(kData == 0) histoName = "fHistPN"; if(kData == 1) histoName = "fHistPN_shuffle"; if(kData == 2) histoName = "fHistPN"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistPN = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistPN) { Printf("fHistPN %s not found!!!",histoName.Data()); break; } fHistPN->FillParent(); fHistPN->DeleteContainers(); if(kData == 0) histoName = "fHistNP"; if(kData == 1) histoName = "fHistNP_shuffle"; if(kData == 2) histoName = "fHistNP"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistNP = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistNP) { Printf("fHistNP %s not found!!!",histoName.Data()); break; } fHistNP->FillParent(); fHistNP->DeleteContainers(); if(kData == 0) histoName = "fHistPP"; if(kData == 1) histoName = "fHistPP_shuffle"; if(kData == 2) histoName = "fHistPP"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistPP = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistPP) { Printf("fHistPP %s not found!!!",histoName.Data()); break; } fHistPP->FillParent(); fHistPP->DeleteContainers(); if(kData == 0) histoName = "fHistNN"; if(kData == 1) histoName = "fHistNN_shuffle"; if(kData == 2) histoName = "fHistNN"; if(gCentralityEstimator) histoName += gCentralityEstimator; AliTHn *fHistNN = dynamic_cast<AliTHn *>(listBF->FindObject(histoName.Data())); if(!fHistNN) { Printf("fHistNN %s not found!!!",histoName.Data()); break; } fHistNN->FillParent(); fHistNN->DeleteContainers(); dir->cd(); listBF->Write(Form("%s_histograms",listBFName.Data()), TObject::kSingleKey); }// first iteration f->Close(); return listBF; } //______________________________________________________// void draw(TList *list, TList *listBFShuffled, TList *listBFMixed, TList *listQA, const char *gCentralityEstimator, Int_t gCentrality, Double_t psiMin, Double_t psiMax, Double_t vertexZMin,Double_t vertexZMax, Double_t ptTriggerMin, Double_t ptTriggerMax, Double_t ptAssociatedMin, Double_t ptAssociatedMax, Bool_t normToTrig, Double_t normalizationRangePhi, Bool_t kUseVzBinning, Int_t rebinEta, Int_t rebinPhi,TString eventClass,Bool_t bToy) { //Draws the correlation functions for every centrality bin //(+-), (-+), (++), (--) AliTHn *hP = NULL; AliTHn *hN = NULL; AliTHn *hPN = NULL; AliTHn *hNP = NULL; AliTHn *hPP = NULL; AliTHn *hNN = NULL; TString gHistPname = "fHistP"; if(gCentralityEstimator) gHistPname += gCentralityEstimator; TString gHistNname = "fHistN"; if(gCentralityEstimator) gHistNname += gCentralityEstimator; TString gHistPNname = "fHistPN"; if(gCentralityEstimator) gHistPNname += gCentralityEstimator; TString gHistNPname = "fHistNP"; if(gCentralityEstimator) gHistNPname += gCentralityEstimator; TString gHistPPname = "fHistPP"; if(gCentralityEstimator) gHistPPname += gCentralityEstimator; TString gHistNNname = "fHistNN"; if(gCentralityEstimator) gHistNNname += gCentralityEstimator; hP = (AliTHn*) list->FindObject(gHistPname.Data()); hN = (AliTHn*) list->FindObject(gHistNname.Data()); hPN = (AliTHn*) list->FindObject(gHistPNname.Data()); hNP = (AliTHn*) list->FindObject(gHistNPname.Data()); hPP = (AliTHn*) list->FindObject(gHistPPname.Data()); hNN = (AliTHn*) list->FindObject(gHistNNname.Data()); hNN->Print(); //Create the AliBalancePsi object and fill it with the AliTHn objects AliBalancePsi *b = new AliBalancePsi(); b->SetEventClass(eventClass); b->SetHistNp(hP); b->SetHistNn(hN); b->SetHistNpn(hPN); b->SetHistNnp(hNP); b->SetHistNpp(hPP); b->SetHistNnn(hNN); if(kUseVzBinning) b->SetVertexZBinning(kTRUE); //balance function shuffling AliTHn *hPShuffled = NULL; AliTHn *hNShuffled = NULL; AliTHn *hPNShuffled = NULL; AliTHn *hNPShuffled = NULL; AliTHn *hPPShuffled = NULL; AliTHn *hNNShuffled = NULL; if(listBFShuffled) { //listBFShuffled->ls(); gHistPname = "fHistP_shuffle"; if(gCentralityEstimator) gHistPname += gCentralityEstimator; gHistNname = "fHistN_shuffle"; if(gCentralityEstimator) gHistNname += gCentralityEstimator; gHistPNname = "fHistPN_shuffle"; if(gCentralityEstimator) gHistPNname += gCentralityEstimator; gHistNPname = "fHistNP_shuffle"; if(gCentralityEstimator) gHistNPname += gCentralityEstimator; gHistPPname = "fHistPP_shuffle"; if(gCentralityEstimator) gHistPPname += gCentralityEstimator; gHistNNname = "fHistNN_shuffle"; if(gCentralityEstimator) gHistNNname += gCentralityEstimator; hPShuffled = (AliTHn*) listBFShuffled->FindObject(gHistPname.Data()); hPShuffled->SetName("gHistPShuffled"); hNShuffled = (AliTHn*) listBFShuffled->FindObject(gHistNname.Data()); hNShuffled->SetName("gHistNShuffled"); hPNShuffled = (AliTHn*) listBFShuffled->FindObject(gHistPNname.Data()); hPNShuffled->SetName("gHistPNShuffled"); hNPShuffled = (AliTHn*) listBFShuffled->FindObject(gHistNPname.Data()); hNPShuffled->SetName("gHistNPShuffled"); hPPShuffled = (AliTHn*) listBFShuffled->FindObject(gHistPPname.Data()); hPPShuffled->SetName("gHistPPShuffled"); hNNShuffled = (AliTHn*) listBFShuffled->FindObject(gHistNNname.Data()); hNNShuffled->SetName("gHistNNShuffled"); AliBalancePsi *bShuffled = new AliBalancePsi(); bShuffled->SetEventClass(eventClass); bShuffled->SetHistNp(hPShuffled); bShuffled->SetHistNn(hNShuffled); bShuffled->SetHistNpn(hPNShuffled); bShuffled->SetHistNnp(hNPShuffled); bShuffled->SetHistNpp(hPPShuffled); bShuffled->SetHistNnn(hNNShuffled); if(kUseVzBinning) bShuffled->SetVertexZBinning(kTRUE); } //balance function mixing AliTHn *hPMixed = NULL; AliTHn *hNMixed = NULL; AliTHn *hPNMixed = NULL; AliTHn *hNPMixed = NULL; AliTHn *hPPMixed = NULL; AliTHn *hNNMixed = NULL; if(listBFMixed) { //listBFMixed->ls(); gHistPname = "fHistP"; if(gCentralityEstimator) gHistPname += gCentralityEstimator; gHistNname = "fHistN"; if(gCentralityEstimator) gHistNname += gCentralityEstimator; gHistPNname = "fHistPN"; if(gCentralityEstimator) gHistPNname += gCentralityEstimator; gHistNPname = "fHistNP"; if(gCentralityEstimator) gHistNPname += gCentralityEstimator; gHistPPname = "fHistPP"; if(gCentralityEstimator) gHistPPname += gCentralityEstimator; gHistNNname = "fHistNN"; if(gCentralityEstimator) gHistNNname += gCentralityEstimator; hPMixed = (AliTHn*) listBFMixed->FindObject(gHistPname.Data()); hPMixed->SetName("gHistPMixed"); hNMixed = (AliTHn*) listBFMixed->FindObject(gHistNname.Data()); hNMixed->SetName("gHistNMixed"); hPNMixed = (AliTHn*) listBFMixed->FindObject(gHistPNname.Data()); hPNMixed->SetName("gHistPNMixed"); hNPMixed = (AliTHn*) listBFMixed->FindObject(gHistNPname.Data()); hNPMixed->SetName("gHistNPMixed"); hPPMixed = (AliTHn*) listBFMixed->FindObject(gHistPPname.Data()); hPPMixed->SetName("gHistPPMixed"); hNNMixed = (AliTHn*) listBFMixed->FindObject(gHistNNname.Data()); hNNMixed->SetName("gHistNNMixed"); AliBalancePsi *bMixed = new AliBalancePsi(); bMixed->SetEventClass(eventClass); bMixed->SetHistNp(hPMixed); bMixed->SetHistNn(hNMixed); bMixed->SetHistNpn(hPNMixed); bMixed->SetHistNnp(hNPMixed); bMixed->SetHistNpp(hPPMixed); bMixed->SetHistNnn(hNNMixed); if(kUseVzBinning) bMixed->SetVertexZBinning(kTRUE); } TH2D *gHistPN[4]; TH2D *gHistNP[4]; TH2D *gHistPP[4]; TH2D *gHistNN[4]; TCanvas *cPN[4]; TCanvas *cNP[4]; TCanvas *cPP[4]; TCanvas *cNN[4]; TString histoTitle, pngName; // if no mixing then divide by convoluted histograms if(!listBFMixed && listQA){ if(!listQA->FindObject("fHistEtaPhiPos") || !listQA->FindObject("fHistEtaPhiNeg")){ Printf("fHistEtaPhiPos or fHistEtaPhiNeg not found! --> no convolution"); listQA = NULL; } else{ TH2D* fHistPos = (TH2D*)((TH3D*)listQA->FindObject("fHistEtaPhiPos"))->Project3D("xy"); fHistPos->GetYaxis()->SetRangeUser(-0.79,0.79); TH2D* fHistNeg = (TH2D*)((TH3D*)listQA->FindObject("fHistEtaPhiNeg"))->Project3D("xy"); fHistNeg->GetYaxis()->SetRangeUser(-0.79,0.79); gHistPN[2] = convolute2D(fHistPos, fHistNeg, "hConvPN"); gHistPN[2]->Scale(1./gHistPN[2]->GetBinContent(gHistPN[2]->FindBin(0,0))); gHistNP[2] = convolute2D(fHistNeg, fHistPos, "hConvNP"); gHistNP[2]->Scale(1./gHistNP[2]->GetBinContent(gHistNP[2]->FindBin(0,0))); gHistNN[2] = convolute2D(fHistNeg, fHistNeg, "hConvNN"); gHistNN[2]->Scale(1./gHistNN[2]->GetBinContent(gHistNN[2]->FindBin(0,0))); gHistPP[2] = convolute2D(fHistPos, fHistPos, "hConvPP"); gHistPP[2]->Scale(1./gHistPP[2]->GetBinContent(gHistPP[2]->FindBin(0,0))); } } //(+-) if(eventClass == "Centrality"){ histoTitle = "(+-) | Centrality: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " % "; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else if(eventClass == "Multiplicity"){ histoTitle = "(+-) | Multiplicity: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " tracks"; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else{ // "EventPlane" (default) histoTitle = "(+-) | Centrality: "; histoTitle += centralityArray[gCentrality-1]; histoTitle += "%"; if((psiMin == -0.5)&&(psiMax == 0.5)) histoTitle += " (-7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o})"; else if((psiMin == 0.5)&&(psiMax == 1.5)) histoTitle += " (37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o})"; else if((psiMin == 1.5)&&(psiMax == 2.5)) histoTitle += " (82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o})"; else histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } gHistPN[0] = b->GetCorrelationFunctionPN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPN[0]->Rebin2D(rebinEta,rebinPhi); gHistPN[0]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistPN[0]->GetYaxis()->SetTitleOffset(1.5); gHistPN[0]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPN[0]->SetTitle(histoTitle.Data()); cPN[0] = new TCanvas("cPN0","",0,0,600,500); cPN[0]->SetFillColor(10); cPN[0]->SetHighLightColor(10); gHistPN[0]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEta.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[0]->SaveAs(pngName.Data()); if(listBFShuffled) { histoTitle.ReplaceAll("(+-)","(+-) shuffled"); gHistPN[1] = bShuffled->GetCorrelationFunctionPN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPN[1]->Rebin2D(rebinEta,rebinPhi); gHistPN[1]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistPN[1]->GetYaxis()->SetTitleOffset(1.5); gHistPN[1]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPN[1]->SetTitle(histoTitle.Data()); cPN[1] = new TCanvas("cPN1","",0,100,600,500); cPN[1]->SetFillColor(10); cPN[1]->SetHighLightColor(10); gHistPN[1]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaShuffled.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[1]->SaveAs(pngName.Data()); } if(listBFMixed) { if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(+-) shuffled","(+-) mixed"); else histoTitle.ReplaceAll("(+-)","(+-) mixed"); // if normalization to trigger then do not divide Event mixing by number of trigger particles gHistPN[2] = bMixed->GetCorrelationFunctionPN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPN[2]->Rebin2D(rebinEta,rebinPhi); gHistPN[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistPN[2]->Integral(gHistPN[2]->GetXaxis()->FindBin(0-10e-5),gHistPN[2]->GetXaxis()->FindBin(0+10e-5),gHistPN[2]->GetNbinsY()/2 + 1,gHistPN[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistPN[2]->GetNbinsY()*(gHistPN[2]->GetXaxis()->FindBin(0.01) - gHistPN[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistPN[2]->GetXaxis()->GetBinWidth(gHistPN[2]->GetNbinsX()); Double_t maxEta = gHistPN[2]->GetXaxis()->GetBinUpEdge(gHistPN[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistPN[2]->Scale(1./mixedNorm); } gHistPN[2]->GetYaxis()->SetTitleOffset(1.5); gHistPN[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPN[2]->SetTitle(histoTitle.Data()); cPN[2] = new TCanvas("cPN2","",0,200,600,500); cPN[2]->SetFillColor(10); cPN[2]->SetHighLightColor(10); gHistPN[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[2]->SaveAs(pngName.Data()); //Correlation function (+-) gHistPN[3] = b->GetCorrelationFunction("PN",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi); //gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistPN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistPN[3]->GetZaxis()->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); //gHistPN[3]->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); cPN[3] = new TCanvas("cPN3","",0,300,600,500); cPN[3]->SetFillColor(10); cPN[3]->SetHighLightColor(10); gHistPN[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[3]->SaveAs(pngName.Data()); } // if no mixing then divide by convoluted histograms else if(listQA){ if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(+-) shuffled","(+-) convoluted"); else histoTitle.ReplaceAll("(+-)","(+-) convoluted"); if(rebinEta > 1 || rebinPhi > 1){ gHistPN[2]->Rebin2D(rebinEta,rebinPhi); gHistPN[0]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistPN[2]->Integral(gHistPN[2]->GetXaxis()->FindBin(0-10e-5),gHistPN[2]->GetXaxis()->FindBin(0+10e-5),gHistPN[2]->GetNbinsY()/2 + 1,gHistPN[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistPN[2]->GetNbinsY()*(gHistPN[2]->GetXaxis()->FindBin(0.01) - gHistPN[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistPN[2]->GetXaxis()->GetBinWidth(gHistPN[2]->GetNbinsX()); Double_t maxEta = gHistPN[2]->GetXaxis()->GetBinUpEdge(gHistPN[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistPN[2]->Scale(1./mixedNorm); } gHistPN[2]->GetYaxis()->SetTitleOffset(1.5); gHistPN[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPN[2]->SetTitle(histoTitle.Data()); cPN[2] = new TCanvas("cPN2","",0,200,600,500); cPN[2]->SetFillColor(10); cPN[2]->SetHighLightColor(10); gHistPN[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[2]->SaveAs(pngName.Data()); //Correlation function (+-) gHistPN[3] = dynamic_cast<TH2D *>(gHistPN[0]->Clone()); gHistPN[3]->Divide(gHistPN[2]); //gHistPN[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistPN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistPN[3]->GetZaxis()->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); //gHistPN[3]->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); cPN[3] = new TCanvas("cPN3","",0,300,600,500); cPN[3]->SetFillColor(10); cPN[3]->SetHighLightColor(10); gHistPN[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositiveNegative.png"; //cPN[3]->SaveAs(pngName.Data()); } //(-+) if(eventClass == "Centrality"){ histoTitle = "(-+) | Centrality: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " % "; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else if(eventClass == "Multiplicity"){ histoTitle = "(-+) | Multiplicity: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " tracks"; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else{ // "EventPlane" (default) histoTitle = "(-+) | Centrality: "; histoTitle += centralityArray[gCentrality-1]; histoTitle += "%"; if((psiMin == -0.5)&&(psiMax == 0.5)) histoTitle += " (-7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o})"; else if((psiMin == 0.5)&&(psiMax == 1.5)) histoTitle += " (37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o})"; else if((psiMin == 1.5)&&(psiMax == 2.5)) histoTitle += " (82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o})"; else histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } gHistNP[0] = b->GetCorrelationFunctionNP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNP[0]->Rebin2D(rebinEta,rebinPhi); gHistNP[0]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistNP[0]->GetYaxis()->SetTitleOffset(1.5); gHistNP[0]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNP[0]->SetTitle(histoTitle.Data()); cNP[0] = new TCanvas("cNP0","",100,0,600,500); cNP[0]->SetFillColor(10); cNP[0]->SetHighLightColor(10); gHistNP[0]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEta.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[0]->SaveAs(pngName.Data()); if(listBFShuffled) { histoTitle.ReplaceAll("(-+)","(-+) shuffled"); gHistNP[1] = bShuffled->GetCorrelationFunctionNP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNP[1]->Rebin2D(rebinEta,rebinPhi); gHistNP[1]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistNP[1]->GetYaxis()->SetTitleOffset(1.5); gHistNP[1]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNP[1]->SetTitle(histoTitle.Data()); cNP[1] = new TCanvas("cNP1","",100,100,600,500); cNP[1]->SetFillColor(10); cNP[1]->SetHighLightColor(10); gHistNP[1]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaShuffled.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[1]->SaveAs(pngName.Data()); } if(listBFMixed) { if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(-+) shuffled","(-+) mixed"); else histoTitle.ReplaceAll("(-+)","(-+) mixed"); // if normalization to trigger then do not divide Event mixing by number of trigger particles gHistNP[2] = bMixed->GetCorrelationFunctionNP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNP[2]->Rebin2D(rebinEta,rebinPhi); gHistNP[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistNP[2]->Integral(gHistNP[2]->GetXaxis()->FindBin(0-10e-5),gHistNP[2]->GetXaxis()->FindBin(0+10e-5),gHistNP[2]->GetNbinsY()/2 + 1,gHistNP[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistNP[2]->GetNbinsY()*(gHistNP[2]->GetXaxis()->FindBin(0.01) - gHistNP[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistNP[2]->GetXaxis()->GetBinWidth(gHistNP[2]->GetNbinsX()); Double_t maxEta = gHistNP[2]->GetXaxis()->GetBinUpEdge(gHistNP[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistNP[2]->Scale(1./mixedNorm); } gHistNP[2]->GetYaxis()->SetTitleOffset(1.5); gHistNP[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNP[2]->SetTitle(histoTitle.Data()); cNP[2] = new TCanvas("cNP2","",100,200,600,500); cNP[2]->SetFillColor(10); cNP[2]->SetHighLightColor(10); gHistNP[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[2]->SaveAs(pngName.Data()); //Correlation function (-+) gHistNP[3] = b->GetCorrelationFunction("NP",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi); //gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistNP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistNP[3]->GetZaxis()->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)"); //gHistNP[3]->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); cNP[3] = new TCanvas("cNP3","",100,300,600,500); cNP[3]->SetFillColor(10); cNP[3]->SetHighLightColor(10); gHistNP[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[3]->SaveAs(pngName.Data()); } // if no mixing then divide by convoluted histograms else if(listQA){ if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(-+) shuffled","(-+) convoluted"); else histoTitle.ReplaceAll("(-+)","(-+) convoluted"); if(rebinEta > 1 || rebinPhi > 1){ gHistNP[2]->Rebin2D(rebinEta,rebinPhi); gHistNP[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistNP[2]->Integral(gHistNP[2]->GetXaxis()->FindBin(0-10e-5),gHistNP[2]->GetXaxis()->FindBin(0+10e-5),gHistNP[2]->GetNbinsY()/2 + 1,gHistNP[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistNP[2]->GetNbinsY()*(gHistNP[2]->GetXaxis()->FindBin(0.01) - gHistNP[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistNP[2]->GetXaxis()->GetBinWidth(gHistNP[2]->GetNbinsX()); Double_t maxEta = gHistNP[2]->GetXaxis()->GetBinUpEdge(gHistNP[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistNP[2]->Scale(1./mixedNorm); } gHistNP[2]->GetYaxis()->SetTitleOffset(1.5); gHistNP[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNP[2]->SetTitle(histoTitle.Data()); cNP[2] = new TCanvas("cNP2","",100,200,600,500); cNP[2]->SetFillColor(10); cNP[2]->SetHighLightColor(10); gHistNP[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[2]->SaveAs(pngName.Data()); //Correlation function (-+) gHistNP[3] = dynamic_cast<TH2D *>(gHistNP[0]->Clone()); gHistNP[3]->Divide(gHistNP[2]); //gHistNP[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistNP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistNP[3]->GetZaxis()->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)"); //gHistNP[3]->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)"); cNP[3] = new TCanvas("cNP3","",100,300,600,500); cNP[3]->SetFillColor(10); cNP[3]->SetHighLightColor(10); gHistNP[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativePositive.png"; //cNP[3]->SaveAs(pngName.Data()); } //(++) if(eventClass == "Centrality"){ histoTitle = "(++) | Centrality: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " % "; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else if(eventClass == "Multiplicity"){ histoTitle = "(++) | Multiplicity: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " tracks"; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else{ // "EventPlane" (default) histoTitle = "(++) | Centrality: "; histoTitle += centralityArray[gCentrality-1]; histoTitle += "%"; if((psiMin == -0.5)&&(psiMax == 0.5)) histoTitle += " (-7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o})"; else if((psiMin == 0.5)&&(psiMax == 1.5)) histoTitle += " (37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o})"; else if((psiMin == 1.5)&&(psiMax == 2.5)) histoTitle += " (82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o})"; else histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } gHistPP[0] = b->GetCorrelationFunctionPP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPP[0]->Rebin2D(rebinEta,rebinPhi); gHistPP[0]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistPP[0]->GetYaxis()->SetTitleOffset(1.5); gHistPP[0]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPP[0]->SetTitle(histoTitle.Data()); cPP[0] = new TCanvas("cPP0","",200,0,600,500); cPP[0]->SetFillColor(10); cPP[0]->SetHighLightColor(10); gHistPP[0]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEta.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[0]->SaveAs(pngName.Data()); if(listBFShuffled) { histoTitle.ReplaceAll("(++)","(++) shuffled"); gHistPP[1] = bShuffled->GetCorrelationFunctionPP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPP[1]->Rebin2D(rebinEta,rebinPhi); gHistPP[1]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistPP[1]->GetYaxis()->SetTitleOffset(1.5); gHistPP[1]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPP[1]->SetTitle(histoTitle.Data()); cPP[1] = new TCanvas("cPP1","",200,100,600,500); cPP[1]->SetFillColor(10); cPP[1]->SetHighLightColor(10); gHistPP[1]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaShuffled.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[1]->SaveAs(pngName.Data()); } if(listBFMixed) { if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(++) shuffled","(++) mixed"); else histoTitle.ReplaceAll("(++)","(++) mixed"); // if normalization to trigger then do not divide Event mixing by number of trigger particles gHistPP[2] = bMixed->GetCorrelationFunctionPP(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistPP[2]->Rebin2D(rebinEta,rebinPhi); gHistPP[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistPP[2]->Integral(gHistPP[2]->GetXaxis()->FindBin(0-10e-5),gHistPP[2]->GetXaxis()->FindBin(0+10e-5),gHistPP[2]->GetNbinsY()/2 + 1,gHistPP[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistPP[2]->GetNbinsY()*(gHistPP[2]->GetXaxis()->FindBin(0.01) - gHistPP[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistPP[2]->GetXaxis()->GetBinWidth(gHistPP[2]->GetNbinsX()); Double_t maxEta = gHistPP[2]->GetXaxis()->GetBinUpEdge(gHistPP[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistPP[2]->Scale(1./mixedNorm); } gHistPP[2]->GetYaxis()->SetTitleOffset(1.5); gHistPP[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPP[2]->SetTitle(histoTitle.Data()); cPP[2] = new TCanvas("cPP2","",200,200,600,500); cPP[2]->SetFillColor(10); cPP[2]->SetHighLightColor(10); gHistPP[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[2]->SaveAs(pngName.Data()); //Correlation function (++) gHistPP[3] = b->GetCorrelationFunction("PP",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi); //gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistPP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistPP[3]->GetZaxis()->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); // gHistPP[3]->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); cPP[3] = new TCanvas("cPP3","",200,300,600,500); cPP[3]->SetFillColor(10); cPP[3]->SetHighLightColor(10); gHistPP[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[3]->SaveAs(pngName.Data()); } // if no mixing then divide by convoluted histograms else if(listQA){ if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(++) shuffled","(++) convoluted"); else histoTitle.ReplaceAll("(++)","(++) convoluted"); if(rebinEta > 1 || rebinPhi > 1){ gHistPP[2]->Rebin2D(rebinEta,rebinPhi); gHistPP[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistPP[2]->Integral(gHistPP[2]->GetXaxis()->FindBin(0-10e-5),gHistPP[2]->GetXaxis()->FindBin(0+10e-5),gHistPP[2]->GetNbinsY()/2 + 1,gHistPP[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistPP[2]->GetNbinsY()*(gHistPP[2]->GetXaxis()->FindBin(0.01) - gHistPP[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistPP[2]->GetXaxis()->GetBinWidth(gHistPP[2]->GetNbinsX()); Double_t maxEta = gHistPP[2]->GetXaxis()->GetBinUpEdge(gHistPP[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistPP[2]->Scale(1./mixedNorm); } gHistPP[2]->GetYaxis()->SetTitleOffset(1.5); gHistPP[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistPP[2]->SetTitle(histoTitle.Data()); cPP[2] = new TCanvas("cPP2","",200,200,600,500); cPP[2]->SetFillColor(10); cPP[2]->SetHighLightColor(10); gHistPP[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[2]->SaveAs(pngName.Data()); //Correlation function (++) gHistPP[3] = dynamic_cast<TH2D *>(gHistPP[0]->Clone()); gHistPP[3]->Divide(gHistPP[2]); //gHistPP[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistPP[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistPP[3]->GetZaxis()->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); //gHistPP[3]->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); cPP[3] = new TCanvas("cPP3","",200,300,600,500); cPP[3]->SetFillColor(10); cPP[3]->SetHighLightColor(10); gHistPP[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".PositivePositive.png"; //cPP[3]->SaveAs(pngName.Data()); } //(--) if(eventClass == "Centrality"){ histoTitle = "(--) | Centrality: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " % "; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else if(eventClass == "Multiplicity"){ histoTitle = "(--) | Multiplicity: "; histoTitle += psiMin; histoTitle += " - "; histoTitle += psiMax; histoTitle += " tracks"; histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } else{ // "EventPlane" (default) histoTitle = "(--) | Centrality: "; histoTitle += centralityArray[gCentrality-1]; histoTitle += "%"; if((psiMin == -0.5)&&(psiMax == 0.5)) histoTitle += " (-7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o})"; else if((psiMin == 0.5)&&(psiMax == 1.5)) histoTitle += " (37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o})"; else if((psiMin == 1.5)&&(psiMax == 2.5)) histoTitle += " (82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o})"; else histoTitle += " (0^{o} < #varphi^{t} - #Psi_{2} < 180^{o})"; } gHistNN[0] = b->GetCorrelationFunctionNN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNN[0]->Rebin2D(rebinEta,rebinPhi); gHistNN[0]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistNN[0]->GetYaxis()->SetTitleOffset(1.5); gHistNN[0]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNN[0]->SetTitle(histoTitle.Data()); cNN[0] = new TCanvas("cNN0","",300,0,600,500); cNN[0]->SetFillColor(10); cNN[0]->SetHighLightColor(10); gHistNN[0]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 //gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEta.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[0]->SaveAs(pngName.Data()); if(listBFShuffled) { histoTitle.ReplaceAll("(--)","(--) shuffled"); gHistNN[1] = bShuffled->GetCorrelationFunctionNN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNN[1]->Rebin2D(rebinEta,rebinPhi); gHistNN[1]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } gHistNN[1]->GetYaxis()->SetTitleOffset(1.5); gHistNN[1]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNN[1]->SetTitle(histoTitle.Data()); cNN[1] = new TCanvas("cNN1","",300,100,600,500); cNN[1]->SetFillColor(10); cNN[1]->SetHighLightColor(10); gHistNN[1]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaShuffled.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[1]->SaveAs(pngName.Data()); } if(listBFMixed) { if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(--) shuffled","(--) mixed"); else histoTitle.ReplaceAll("(--)","(--) mixed"); // if normalization to trigger then do not divide Event mixing by number of trigger particles gHistNN[2] = bMixed->GetCorrelationFunctionNN(psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax); if(rebinEta > 1 || rebinPhi > 1){ gHistNN[2]->Rebin2D(rebinEta,rebinPhi); gHistNN[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistNN[2]->Integral(gHistNN[2]->GetXaxis()->FindBin(0-10e-5),gHistNN[2]->GetXaxis()->FindBin(0+10e-5),gHistNN[2]->GetNbinsY()/2 + 1,gHistNN[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistNN[2]->GetNbinsY()*(gHistNN[2]->GetXaxis()->FindBin(0.01) - gHistNN[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistNN[2]->GetXaxis()->GetBinWidth(gHistNN[2]->GetNbinsX()); Double_t maxEta = gHistNN[2]->GetXaxis()->GetBinUpEdge(gHistNN[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistNN[2]->Scale(1./mixedNorm); } gHistNN[2]->GetYaxis()->SetTitleOffset(1.5); gHistNN[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNN[2]->SetTitle(histoTitle.Data()); cNN[2] = new TCanvas("cNN2","",300,200,600,500); cNN[2]->SetFillColor(10); cNN[2]->SetHighLightColor(10); gHistNN[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[2]->SaveAs(pngName.Data()); //Correlation function (--) gHistNN[3] = b->GetCorrelationFunction("NN",psiMin,psiMax,vertexZMin,vertexZMax,ptTriggerMin,ptTriggerMax,ptAssociatedMin,ptAssociatedMax,bMixed,normToTrig,normalizationRangePhi); //gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistNN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistNN[3]->GetZaxis()->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); // gHistNN[3]->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); cNN[3] = new TCanvas("cNN3","",300,300,600,500); cNN[3]->SetFillColor(10); cNN[3]->SetHighLightColor(10); gHistNN[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[3]->SaveAs(pngName.Data()); } // if no mixing then divide by convoluted histograms else if(listQA){ if(histoTitle.Contains("shuffled")) histoTitle.ReplaceAll("(--) shuffled","(--) convoluted"); else histoTitle.ReplaceAll("(--)","(--) convoluted"); if(rebinEta > 1 || rebinPhi > 1){ gHistNN[2]->Rebin2D(rebinEta,rebinPhi); gHistNN[2]->Scale(1./(Double_t)(rebinEta*rebinPhi)); } // normalization to 1 at (0,0) --> Jan Fietes method if(normToTrig){ Double_t mixedNorm = gHistNN[2]->Integral(gHistNN[2]->GetXaxis()->FindBin(0-10e-5),gHistNN[2]->GetXaxis()->FindBin(0+10e-5),gHistNN[2]->GetNbinsY()/2 + 1,gHistNN[2]->GetNbinsY()); mixedNorm /= 0.5 * gHistNN[2]->GetNbinsY()*(gHistNN[2]->GetXaxis()->FindBin(0.01) - gHistNN[2]->GetXaxis()->FindBin(-0.01) + 1); // finite bin correction Double_t binWidthEta = gHistNN[2]->GetXaxis()->GetBinWidth(gHistNN[2]->GetNbinsX()); Double_t maxEta = gHistNN[2]->GetXaxis()->GetBinUpEdge(gHistNN[2]->GetNbinsX()); Double_t finiteBinCorrection = -1.0 / (2*maxEta) * binWidthEta / 2 + 1; //Printf("Finite bin correction: %f", finiteBinCorrection); mixedNorm /= finiteBinCorrection; gHistNN[2]->Scale(1./mixedNorm); } gHistNN[2]->GetYaxis()->SetTitleOffset(1.5); gHistNN[2]->GetYaxis()->SetTitle("#Delta #varphi (rad)"); gHistNN[2]->SetTitle(histoTitle.Data()); cNN[2] = new TCanvas("cNN2","",300,200,600,500); cNN[2]->SetFillColor(10); cNN[2]->SetHighLightColor(10); gHistNN[2]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "DeltaPhiDeltaEtaMixed.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[2]->SaveAs(pngName.Data()); //Correlation function (--) gHistNN[3] = dynamic_cast<TH2D *>(gHistNN[0]->Clone()); gHistNN[3]->Divide(gHistNN[2]); //gHistNN[3]->GetXaxis()->SetRangeUser(-1.5,1.5); if(normToTrig) gHistNN[3]->GetZaxis()->SetTitle("#frac{1}{N_{trig}}#frac{d^{2}N_{assoc}}{d#Delta#eta#Delta#varphi} (rad^{-1})"); else gHistNN[3]->GetZaxis()->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); //gHistNN[3]->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); cNN[3] = new TCanvas("cNN3","",300,300,600,500); cNN[3]->SetFillColor(10); cNN[3]->SetHighLightColor(10); gHistNN[3]->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 //gPad->SetPhi(130); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; pngName += psiMin; pngName += "To"; pngName += psiMax; pngName += ".NegativeNegative.png"; //cNN[3]->SaveAs(pngName.Data()); } //Write to output file TString newFileName = "correlationFunction."; if(eventClass == "Centrality"){ newFileName += Form("Centrality%.1fTo%.1f",psiMin,psiMax); newFileName += ".PsiAll.PttFrom"; } else if(eventClass == "Multiplicity"){ newFileName += Form("Multiplicity%.0fTo%.0f",psiMin,psiMax); newFileName += ".PsiAll.PttFrom"; } else{ // "EventPlane" (default) newFileName += "Centrality"; newFileName += gCentrality; newFileName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) newFileName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) newFileName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) newFileName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) newFileName += "Rest.PttFrom"; else newFileName += "All.PttFrom"; } newFileName += Form("%.1f",ptTriggerMin); newFileName += "To"; newFileName += Form("%.1f",ptTriggerMax); newFileName += "PtaFrom"; newFileName += Form("%.1f",ptAssociatedMin); newFileName += "To"; newFileName += Form("%.1f",ptAssociatedMax); newFileName += "_"; newFileName += Form("%.1f",psiMin); newFileName += "-"; newFileName += Form("%.1f",psiMax); newFileName += ".root"; TFile *newFile = TFile::Open(newFileName.Data(),"recreate"); gHistPN[0]->SetName("gHistPNRaw"); gHistPN[0]->Write(); gHistNP[0]->SetName("gHistNPRaw"); gHistNP[0]->Write(); gHistPP[0]->SetName("gHistPPRaw"); gHistPP[0]->Write(); gHistNN[0]->SetName("gHistNNRaw"); gHistNN[0]->Write(); if(listBFShuffled) { gHistPN[1]->SetName("gHistPNShuffled"); gHistPN[1]->Write(); gHistNP[1]->SetName("gHistNPShuffled"); gHistNP[1]->Write(); gHistPP[1]->SetName("gHistPPShuffled"); gHistPP[1]->Write(); gHistNN[1]->SetName("gHistNNShuffled"); gHistNN[1]->Write(); } if(listBFMixed || (!listBFMixed&&listQA)) { gHistPN[2]->SetName("gHistPNMixed"); gHistPN[2]->Write(); gHistNP[2]->SetName("gHistNPMixed"); gHistNP[2]->Write(); gHistPP[2]->SetName("gHistPPMixed"); gHistPP[2]->Write(); gHistNN[2]->SetName("gHistNNMixed"); gHistNN[2]->Write(); gHistPN[3]->SetName("gHistPNCorrelationFunctions"); gHistPN[3]->Write(); gHistNP[3]->SetName("gHistNPCorrelationFunctions"); gHistNP[3]->Write(); gHistPP[3]->SetName("gHistPPCorrelationFunctions"); gHistPP[3]->Write(); gHistNN[3]->SetName("gHistNNCorrelationFunctions"); gHistNN[3]->Write(); } newFile->Close(); // some cleaning for(Int_t i = 0; i < 4; i++){ if(!listBFShuffled && i == 1) continue; if(!listBFMixed && (i == 2 || i == 3)) continue; if(gHistPP[i]) delete gHistPP[i]; if(gHistPN[i]) delete gHistPN[i]; if(gHistNP[i]) delete gHistNP[i]; if(gHistNN[i]) delete gHistNN[i]; if(cPN[i]) delete cPN[i]; if(cNP[i]) delete cNP[i]; if(cPP[i]) delete cPP[i]; if(cNN[i]) delete cNN[i]; } delete hP; delete hN; delete hPP; delete hPN; delete hNP; delete hNN; delete hPMixed; delete hNMixed; delete hPPMixed; delete hPNMixed; delete hNPMixed; delete hNNMixed; delete hPShuffled; delete hNShuffled; delete hPPShuffled; delete hPNShuffled; delete hNPShuffled; delete hNNShuffled; } //____________________________________________________________// void drawCorrelationFunctions(const char* lhcPeriod = "LHC10h", const char* gCentralityEstimator = "V0M", Int_t gBit = 128, const char* gEventPlaneEstimator = "VZERO", Int_t gCentrality = 1, Double_t psiMin = -0.5, Double_t psiMax = 3.5, Double_t vertexZMin = -10., Double_t vertexZMax = 10., Double_t ptTriggerMin = -1., Double_t ptTriggerMax = -1., Double_t ptAssociatedMin = -1., Double_t ptAssociatedMax = -1., Bool_t kFit = kFALSE) { //Macro that draws the charge dependent correlation functions //for each centrality bin for the different pT of trigger and //associated particles //Author: Panos.Christakoglou@nikhef.nl TGaxis::SetMaxDigits(3); //Get the input file /* TString filename = lhcPeriod; filename += "/Centrality"; filename += gCentralityEstimator; filename += "_Bit"; filename += gBit; filename += "_"; filename += gEventPlaneEstimator; filename +="/PttFrom"; filename += Form("%.1f",ptTriggerMin); filename += "To"; filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom"; filename += Form("%.1f",ptAssociatedMin); filename += "To"; filename += Form("%.1f",ptAssociatedMax); */ TString filename = "correlationFunction.Centrality"; filename += gCentrality; filename += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) filename += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) filename += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) filename += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) filename += "Rest.Ptt"; else filename += "All.PttFrom"; filename += Form("%.1f",ptTriggerMin); filename += "To"; filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom"; filename += Form("%.1f",ptAssociatedMin); filename += "To"; filename += Form("%.1f",ptAssociatedMax); filename += "_"; filename += Form("%.1f",psiMin); filename += "-"; filename += Form("%.1f",psiMax); filename += ".root"; //Open the file TFile *f = TFile::Open(filename.Data()); if((!f)||(!f->IsOpen())) { Printf("The file %s is not found. Aborting...",filename); return; } //f->ls(); //Latex TString centralityLatex = "Centrality: "; centralityLatex += centralityArray[gCentrality-1]; centralityLatex += "%"; TString psiLatex; if((psiMin == -0.5)&&(psiMax == 0.5)) psiLatex = " -7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o}"; else if((psiMin == 0.5)&&(psiMax == 1.5)) psiLatex = " 37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o}"; else if((psiMin == 1.5)&&(psiMax == 2.5)) psiLatex = " 82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o}"; else psiLatex = " 0^{o} < #varphi^{t} - #Psi_{2} < 180^{o}"; TString pttLatex = Form("%.1f",ptTriggerMin); pttLatex += " < p_{T,trig} < "; pttLatex += Form("%.1f",ptTriggerMax); pttLatex += " GeV/c"; TString ptaLatex = Form("%.1f",ptAssociatedMin); ptaLatex += " < p_{T,assoc} < "; ptaLatex += Form("%.1f",ptAssociatedMax); ptaLatex += " GeV/c"; TLatex *latexInfo1 = new TLatex(); latexInfo1->SetNDC(); latexInfo1->SetTextSize(0.045); latexInfo1->SetTextColor(1); TString pngName; //============================================================// //Get the +- correlation function TH2D *gHistPN = dynamic_cast<TH2D *>(f->Get("gHistPNCorrelationFunctions")); gHistPN->SetStats(kFALSE); gHistPN->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); //gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4); gHistPN->GetXaxis()->CenterTitle(); gHistPN->GetXaxis()->SetTitleOffset(1.2); gHistPN->GetYaxis()->CenterTitle(); gHistPN->GetYaxis()->SetTitleOffset(1.2); gHistPN->GetZaxis()->SetTitleOffset(1.5); TCanvas *cPN = new TCanvas("cPN","",0,0,600,500); cPN->SetFillColor(10); cPN->SetHighLightColor(10); cPN->SetLeftMargin(0.15); gHistPN->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) pngName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) pngName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) pngName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) pngName += "Rest.Ptt"; else pngName += "All.PttFrom"; pngName += Form("%.1f",ptTriggerMin); pngName += "To"; pngName += Form("%.1f",ptTriggerMax); pngName += "PtaFrom"; pngName += Form("%.1f",ptAssociatedMin); pngName += "To"; pngName += Form("%.1f",ptAssociatedMax); pngName += ".PositiveNegative.png"; cPN->SaveAs(pngName.Data()); if(kFit) fitCorrelationFunctions(gCentrality, psiMin, psiMax,vertexZMin, vertexZMax, ptTriggerMin,ptTriggerMax, ptAssociatedMin, ptAssociatedMax,gHistPN); //============================================================// //Get the -+ correlation function TH2D *gHistNP = dynamic_cast<TH2D *>(f->Get("gHistNPCorrelationFunctions")); gHistNP->SetStats(kFALSE); gHistNP->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)"); //gHistNP->GetXaxis()->SetRangeUser(-1.4,1.4); gHistNP->GetXaxis()->CenterTitle(); gHistNP->GetXaxis()->SetTitleOffset(1.2); gHistNP->GetYaxis()->CenterTitle(); gHistNP->GetYaxis()->SetTitleOffset(1.2); gHistNP->GetZaxis()->SetTitleOffset(1.5); TCanvas *cNP = new TCanvas("cNP","",50,50,600,500); cNP->SetFillColor(10); cNP->SetHighLightColor(10); cNP->SetLeftMargin(0.15); gHistNP->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) pngName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) pngName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) pngName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) pngName += "Rest.Ptt"; else pngName += "All.PttFrom"; pngName += Form("%.1f",ptTriggerMin); pngName += "To"; pngName += Form("%.1f",ptTriggerMax); pngName += "PtaFrom"; pngName += Form("%.1f",ptAssociatedMin); pngName += "To"; pngName += Form("%.1f",ptAssociatedMax); pngName += ".NegativePositive.png"; cNP->SaveAs(pngName.Data()); if(kFit) fitCorrelationFunctions(gCentrality, psiMin, psiMax,vertexZMin, vertexZMax, ptTriggerMin,ptTriggerMax, ptAssociatedMin, ptAssociatedMax,gHistNP); //============================================================// //Get the ++ correlation function TH2D *gHistPP = dynamic_cast<TH2D *>(f->Get("gHistPPCorrelationFunctions")); gHistPP->SetStats(kFALSE); gHistPP->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); //gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4); gHistPP->GetXaxis()->CenterTitle(); gHistPP->GetXaxis()->SetTitleOffset(1.2); gHistPP->GetYaxis()->CenterTitle(); gHistPP->GetYaxis()->SetTitleOffset(1.2); gHistPP->GetZaxis()->SetTitleOffset(1.5); TCanvas *cPP = new TCanvas("cPP","",100,100,600,500); cPP->SetFillColor(10); cPP->SetHighLightColor(10); cPP->SetLeftMargin(0.15); gHistPP->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) pngName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) pngName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) pngName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) pngName += "Rest.Ptt"; else pngName += "All.PttFrom"; pngName += Form("%.1f",ptTriggerMin); pngName += "To"; pngName += Form("%.1f",ptTriggerMax); pngName += "PtaFrom"; pngName += Form("%.1f",ptAssociatedMin); pngName += "To"; pngName += Form("%.1f",ptAssociatedMax); pngName += ".PositivePositive.png"; cPP->SaveAs(pngName.Data()); if(kFit) fitCorrelationFunctions(gCentrality, psiMin, psiMax,vertexZMin, vertexZMax, ptTriggerMin,ptTriggerMax, ptAssociatedMin, ptAssociatedMax,gHistPP); //============================================================// //Get the -- correlation function TH2D *gHistNN = dynamic_cast<TH2D *>(f->Get("gHistNNCorrelationFunctions")); gHistNN->SetStats(kFALSE); gHistNN->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); //gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4); gHistNN->GetXaxis()->CenterTitle(); gHistNN->GetXaxis()->SetTitleOffset(1.2); gHistNN->GetYaxis()->CenterTitle(); gHistNN->GetYaxis()->SetTitleOffset(1.2); gHistNN->GetZaxis()->SetTitleOffset(1.5); TCanvas *cNN = new TCanvas("cNN","",150,150,600,500); cNN->SetFillColor(10); cNN->SetHighLightColor(10); cNN->SetLeftMargin(0.15); gHistNN->DrawCopy("surf1fb"); gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = "CorrelationFunction.Centrality"; pngName += centralityArray[gCentrality-1]; pngName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) pngName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) pngName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) pngName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) pngName += "Rest.Ptt"; else pngName += "All.PttFrom"; pngName += Form("%.1f",ptTriggerMin); pngName += "To"; pngName += Form("%.1f",ptTriggerMax); pngName += "PtaFrom"; pngName += Form("%.1f",ptAssociatedMin); pngName += "To"; pngName += Form("%.1f",ptAssociatedMax); pngName += ".NegativeNegative.png"; cNN->SaveAs(pngName.Data()); if(kFit) fitCorrelationFunctions(gCentrality, psiMin, psiMax,vertexZMin, vertexZMax, ptTriggerMin,ptTriggerMax, ptAssociatedMin, ptAssociatedMax,gHistNN); } //____________________________________________________________// void drawProjections(Bool_t kProjectInEta = kFALSE, Int_t binMin = 1, Int_t binMax = 80, Int_t gCentrality = 1, Double_t psiMin = -0.5, Double_t psiMax = 3.5, Double_t ptTriggerMin = -1., Double_t ptTriggerMax = -1., Double_t ptAssociatedMin = -1., Double_t ptAssociatedMax = -1., Bool_t kUseZYAM = kFALSE, TString eventClass = "Centrality") { //Macro that draws the charge dependent correlation functions PROJECTIONS //for each centrality bin for the different pT of trigger and //associated particles TGaxis::SetMaxDigits(3); //Get the input file /*TString filename = lhcPeriod; filename += "/Centrality"; filename += gCentralityEstimator; filename += "_Bit"; filename += gBit; filename += "_"; filename += gEventPlaneEstimator; filename +="/PttFrom"; filename += Form("%.1f",ptTriggerMin); filename += "To"; filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom"; filename += Form("%.1f",ptAssociatedMin); filename += "To"; filename += Form("%.1f",ptAssociatedMax); */ TString filename = "correlationFunction."; if(eventClass == "Centrality"){ filename += Form("Centrality%.1fTo%.1f",psiMin,psiMax); filename += ".PsiAll.PttFrom"; } else if(eventClass == "Multiplicity"){ filename += Form("Multiplicity%.0fTo%.0f",psiMin,psiMax); filename += ".PsiAll.PttFrom"; } else{ // "EventPlane" (default) filename += "Centrality"; filename += gCentrality; filename += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) filename += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) filename += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) filename += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) filename += "Rest.PttFrom"; else filename += "All.PttFrom"; } filename += Form("%.1f",ptTriggerMin); filename += "To"; filename += Form("%.1f",ptTriggerMax); filename += "PtaFrom"; filename += Form("%.1f",ptAssociatedMin); filename += "To"; filename += Form("%.1f",ptAssociatedMax); //if(k2pMethod) filename += "_2pMethod"; filename += "_"; filename += Form("%.1f",psiMin); filename += "-"; filename += Form("%.1f",psiMax); filename += ".root"; //Open the file TFile *f = TFile::Open(filename.Data()); if((!f)||(!f->IsOpen())) { Printf("The file %s is not found. Aborting...",filename); return listBF; } //f->ls(); //Latex TString centralityLatex = "Centrality: "; if(eventClass == "Centrality") centralityLatex += Form("%.0f - %.0f",psiMin,psiMax); else centralityLatex += centralityArray[gCentrality-1]; centralityLatex += " %"; TString psiLatex; if((psiMin == -0.5)&&(psiMax == 0.5)) psiLatex = " -7.5^{o} < #varphi^{t} - #Psi_{2} < 7.5^{o}"; else if((psiMin == 0.5)&&(psiMax == 1.5)) psiLatex = " 37.5^{o} < #varphi^{t} - #Psi_{2} < 52.5^{o}"; else if((psiMin == 1.5)&&(psiMax == 2.5)) psiLatex = " 82.5^{o} < #varphi^{t} - #Psi_{2} < 97.5^{o}"; else psiLatex = " 0^{o} < #varphi^{t} - #Psi_{2} < 180^{o}"; TString pttLatex = Form("%.1f",ptTriggerMin); pttLatex += " < p_{T,trig} < "; pttLatex += Form("%.1f",ptTriggerMax); pttLatex += " GeV/c"; TString ptaLatex = Form("%.1f",ptAssociatedMin); ptaLatex += " < p_{T,assoc} < "; ptaLatex += Form("%.1f",ptAssociatedMax); ptaLatex += " GeV/c"; TLatex *latexInfo1 = new TLatex(); latexInfo1->SetNDC(); latexInfo1->SetTextSize(0.045); latexInfo1->SetTextColor(1); TString pngName; //============================================================// //Get the +- correlation function TH2D *gHistPN = dynamic_cast<TH2D *>(f->Get("gHistPNCorrelationFunctions")); gHistPN->SetStats(kFALSE); gHistPN->SetTitle("C_{+-}(#Delta#eta,#Delta#varphi)"); //gHistPN->GetXaxis()->SetRangeUser(-1.4,1.4); gHistPN->GetXaxis()->CenterTitle(); gHistPN->GetXaxis()->SetTitleOffset(1.2); gHistPN->GetYaxis()->CenterTitle(); gHistPN->GetYaxis()->SetTitleOffset(1.2); gHistPN->GetZaxis()->SetTitleOffset(1.5); TCanvas *cPN = new TCanvas("cPN","",0,0,600,500); cPN->SetFillColor(10); cPN->SetHighLightColor(10); cPN->SetLeftMargin(0.15); //================// TH1D* gHistPNprojection = 0x0; Double_t sum = 0.0; Double_t gError = 0.0; Int_t nCounter = 0; //projection in delta eta if(kProjectInEta) { gHistPNprojection = new TH1D("gHistPNprojection","",gHistPN->GetNbinsX(),gHistPN->GetXaxis()->GetXmin(),gHistPN->GetXaxis()->GetXmax()); for(Int_t iBinX = 1; iBinX <= gHistPN->GetNbinsX(); iBinX++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinY = binMin; iBinY <= binMax; iBinY++) { sum += gHistPN->GetBinContent(iBinX,iBinY); if(gHistPN->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistPN->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistPNprojection->SetBinContent(iBinX,sum); gHistPNprojection->SetBinError(iBinX,gError); } //gHistPNprojection->GetXaxis()->SetRangeUser(-1.4,1.4); if(kUseZYAM) gHistPNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}"); else gHistPNprojection->GetYaxis()->SetTitle("C_{+-}(#Delta#eta)"); gHistPNprojection->GetXaxis()->SetTitle("#Delta#eta"); }//projection in delta eta //projection in delta phi else { gHistPNprojection = new TH1D("gHistPNprojection","",gHistPN->GetNbinsY(),gHistPN->GetYaxis()->GetXmin(),gHistPN->GetYaxis()->GetXmax()); for(Int_t iBinY = 1; iBinY <= gHistPN->GetNbinsY(); iBinY++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinX = binMin; iBinX <= binMax; iBinX++) { //for(Int_t iBinX = 1; iBinX <= gHistPN->GetNbinsX(); iBinX++) { sum += gHistPN->GetBinContent(iBinX,iBinY); if(gHistPN->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistPN->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistPNprojection->SetBinContent(iBinY,sum); gHistPNprojection->SetBinError(iBinY,gError); } if(kUseZYAM) gHistPNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#varphi} - b_{ZYAM} (rad^{-1})"); else gHistPNprojection->GetYaxis()->SetTitle("C_{+-}(#Delta#varphi)"); gHistPNprojection->GetXaxis()->SetTitle("#Delta#varphi (rad)"); } //ZYAM if(kUseZYAM) { Double_t reference = gHistPNprojection->GetBinContent(gHistPNprojection->GetMinimumBin()); for(Int_t iBinX = 1; iBinX <= gHistPNprojection->GetNbinsX(); iBinX++) gHistPNprojection->SetBinContent(iBinX,gHistPNprojection->GetBinContent(iBinX) - reference); } gHistPNprojection->GetYaxis()->SetTitleOffset(1.5); gHistPNprojection->SetMarkerStyle(20); gHistPNprojection->SetStats(kFALSE); gHistPNprojection->DrawCopy("E"); //=================// gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = filename; if(kProjectInEta) pngName.ReplaceAll(".root","_DeltaEtaProjection"); else pngName.ReplaceAll(".root","_DeltaPhiProjection"); pngName += ".PositiveNegative.png"; cPN->SaveAs(pngName.Data()); //============================================================// //Get the -+ correlation function TH2D *gHistNP = dynamic_cast<TH2D *>(f->Get("gHistNPCorrelationFunctions")); gHistNP->SetStats(kFALSE); gHistNP->SetTitle("C_{-+}(#Delta#eta,#Delta#varphi)"); //gHistNP->GetXaxis()->SetRangeUser(-1.4,1.4); gHistNP->GetXaxis()->CenterTitle(); gHistNP->GetXaxis()->SetTitleOffset(1.2); gHistNP->GetYaxis()->CenterTitle(); gHistNP->GetYaxis()->SetTitleOffset(1.2); gHistNP->GetZaxis()->SetTitleOffset(1.5); TCanvas *cNP = new TCanvas("cNP","",50,50,600,500); cNP->SetFillColor(10); cNP->SetHighLightColor(10); cNP->SetLeftMargin(0.15); //================// TH1D* gHistNPprojection = 0x0; Double_t sum = 0.0; Double_t gError = 0.0; Int_t nCounter = 0; if(kProjectInEta) { gHistNPprojection = new TH1D("gHistNPprojection","",gHistNP->GetNbinsX(),gHistNP->GetXaxis()->GetXmin(),gHistNP->GetXaxis()->GetXmax()); for(Int_t iBinX = 1; iBinX <= gHistNP->GetNbinsX(); iBinX++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinY = binMin; iBinY <= binMax; iBinY++) { //for(Int_t iBinY = 1; iBinY <= gHistNP->GetNbinsY(); iBinY++) { sum += gHistNP->GetBinContent(iBinX,iBinY); if(gHistNP->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistNP->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistNPprojection->SetBinContent(iBinX,sum); gHistNPprojection->SetBinError(iBinX,gError); } //gHistNPprojection->GetXaxis()->SetRangeUser(-1.4,1.4); if(kUseZYAM) gHistNPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}"); else gHistNPprojection->GetYaxis()->SetTitle("C_{-+}(#Delta#eta)"); gHistNPprojection->GetXaxis()->SetTitle("#Delta#eta"); } else { gHistNPprojection = new TH1D("gHistNPprojection","",gHistNP->GetNbinsY(),gHistNP->GetYaxis()->GetXmin(),gHistNP->GetYaxis()->GetXmax()); for(Int_t iBinY = 1; iBinY <= gHistNP->GetNbinsY(); iBinY++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinX = binMin; iBinX <= binMax; iBinX++) { //for(Int_t iBinX = 1; iBinX <= gHistNP->GetNbinsX(); iBinX++) { sum += gHistNP->GetBinContent(iBinX,iBinY); if(gHistNP->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistNP->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistNPprojection->SetBinContent(iBinY,sum); gHistNPprojection->SetBinError(iBinY,gError); } if(kUseZYAM) gHistNPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#varphi} - b_{ZYAM} (rad^{-1})"); else gHistNPprojection->GetYaxis()->SetTitle("C_{-+}(#Delta#varphi)"); gHistNPprojection->GetXaxis()->SetTitle("#Delta#varphi (rad)"); } //ZYAM if(kUseZYAM) { Double_t reference = gHistNPprojection->GetBinContent(gHistNPprojection->GetMinimumBin()); for(Int_t iBinX = 1; iBinX <= gHistNPprojection->GetNbinsX(); iBinX++) gHistNPprojection->SetBinContent(iBinX,gHistNPprojection->GetBinContent(iBinX) - reference); } gHistNPprojection->GetYaxis()->SetTitleOffset(1.5); gHistNPprojection->SetMarkerStyle(20); gHistNPprojection->SetStats(kFALSE); gHistNPprojection->DrawCopy("E"); //================// gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = filename; if(kProjectInEta) pngName.ReplaceAll(".root","_DeltaEtaProjection"); else pngName.ReplaceAll(".root","_DeltaPhiProjection"); pngName += ".NegativePositive.png"; cNP->SaveAs(pngName.Data()); //============================================================// //Get the ++ correlation function TH2D *gHistPP = dynamic_cast<TH2D *>(f->Get("gHistPPCorrelationFunctions")); gHistPP->SetStats(kFALSE); gHistPP->SetTitle("C_{++}(#Delta#eta,#Delta#varphi)"); //gHistPP->GetXaxis()->SetRangeUser(-1.4,1.4); gHistPP->GetXaxis()->CenterTitle(); gHistPP->GetXaxis()->SetTitleOffset(1.2); gHistPP->GetYaxis()->CenterTitle(); gHistPP->GetYaxis()->SetTitleOffset(1.2); gHistPP->GetZaxis()->SetTitleOffset(1.5); TCanvas *cPP = new TCanvas("cPP","",100,100,600,500); cPP->SetFillColor(10); cPP->SetHighLightColor(10); cPP->SetLeftMargin(0.15); //=================// TH1D* gHistPPprojection = 0x0; Double_t sum = 0.0; Double_t gError = 0.0; Int_t nCounter = 0; if(kProjectInEta) { gHistPPprojection = new TH1D("gHistPPprojection","",gHistPP->GetNbinsX(),gHistPP->GetXaxis()->GetXmin(),gHistPP->GetXaxis()->GetXmax()); for(Int_t iBinX = 1; iBinX <= gHistPP->GetNbinsX(); iBinX++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinY = binMin; iBinY <= binMax; iBinY++) { //for(Int_t iBinY = 1; iBinY <= gHistPP->GetNbinsY(); iBinY++) { sum += gHistPP->GetBinContent(iBinX,iBinY); if(gHistPP->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistPP->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistPPprojection->SetBinContent(iBinX,sum); gHistPPprojection->SetBinError(iBinX,gError); } //gHistPPprojection->GetXaxis()->SetRangeUser(-1.4,1.4); if(kUseZYAM) gHistPPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}"); else gHistPPprojection->GetYaxis()->SetTitle("C_{++}(#Delta#eta)"); gHistPPprojection->GetXaxis()->SetTitle("#Delta#eta"); } else { gHistPPprojection = new TH1D("gHistPPprojection","",gHistPP->GetNbinsY(),gHistPP->GetYaxis()->GetXmin(),gHistPP->GetYaxis()->GetXmax()); for(Int_t iBinY = 1; iBinY <= gHistPP->GetNbinsY(); iBinY++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinX = binMin; iBinX <= binMax; iBinX++) { //for(Int_t iBinX = 1; iBinX <= gHistPP->GetNbinsX(); iBinX++) { sum += gHistPP->GetBinContent(iBinX,iBinY); if(gHistPP->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistPP->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistPPprojection->SetBinContent(iBinY,sum); gHistPPprojection->SetBinError(iBinY,gError); } if(kUseZYAM) gHistPPprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#varphi} - b_{ZYAM} (rad^{-1})"); else gHistPPprojection->GetYaxis()->SetTitle("C_{++}(#Delta#varphi)"); gHistPPprojection->GetXaxis()->SetTitle("#Delta#varphi (rad)"); } //ZYAM if(kUseZYAM) { Double_t reference = gHistPPprojection->GetBinContent(gHistPPprojection->GetMinimumBin()); for(Int_t iBinX = 1; iBinX <= gHistPPprojection->GetNbinsX(); iBinX++) gHistPPprojection->SetBinContent(iBinX,gHistPPprojection->GetBinContent(iBinX) - reference); } gHistPPprojection->GetYaxis()->SetTitleOffset(1.5); gHistPPprojection->SetMarkerStyle(20); gHistPPprojection->SetStats(kFALSE); gHistPPprojection->DrawCopy("E"); //================// gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = filename; if(kProjectInEta) pngName.ReplaceAll(".root","_DeltaEtaProjection"); else pngName.ReplaceAll(".root","_DeltaPhiProjection"); pngName += ".PositivePositive.png"; cPP->SaveAs(pngName.Data()); //============================================================// //Get the -- correlation function TH2D *gHistNN = dynamic_cast<TH2D *>(f->Get("gHistNNCorrelationFunctions")); gHistNN->SetStats(kFALSE); gHistNN->SetTitle("C_{--}(#Delta#eta,#Delta#varphi)"); //gHistNN->GetXaxis()->SetRangeUser(-1.4,1.4); gHistNN->GetXaxis()->CenterTitle(); gHistNN->GetXaxis()->SetTitleOffset(1.2); gHistNN->GetYaxis()->CenterTitle(); gHistNN->GetYaxis()->SetTitleOffset(1.2); gHistNN->GetZaxis()->SetTitleOffset(1.5); TCanvas *cNN = new TCanvas("cNN","",150,150,600,500); cNN->SetFillColor(10); cNN->SetHighLightColor(10); cNN->SetLeftMargin(0.15); //=================// TH1D* gHistNNprojection = 0x0; Double_t sum = 0.0; Double_t gError = 0.0; Int_t nCounter = 0; if(kProjectInEta) { gHistNNprojection = new TH1D("gHistNNprojection","",gHistNN->GetNbinsX(),gHistNN->GetXaxis()->GetXmin(),gHistNN->GetXaxis()->GetXmax()); for(Int_t iBinX = 1; iBinX <= gHistNN->GetNbinsX(); iBinX++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinY = binMin; iBinY <= binMax; iBinY++) { //for(Int_t iBinY = 1; iBinY <= gHistNN->GetNbinsY(); iBinY++) { sum += gHistNN->GetBinContent(iBinX,iBinY); if(gHistNN->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistNN->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistNNprojection->SetBinContent(iBinX,sum); gHistNNprojection->SetBinError(iBinX,gError); } //gHistNNprojection->GetXaxis()->SetRangeUser(-1.4,1.4); if(kUseZYAM) gHistNNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#eta} - b_{ZYAM}"); else gHistNNprojection->GetYaxis()->SetTitle("C_{--}(#Delta#eta)"); gHistNNprojection->GetXaxis()->SetTitle("#Delta#eta"); } else { gHistNNprojection = new TH1D("gHistNNprojection","",gHistNN->GetNbinsY(),gHistNN->GetYaxis()->GetXmin(),gHistNN->GetYaxis()->GetXmax()); for(Int_t iBinY = 1; iBinY <= gHistNN->GetNbinsY(); iBinY++) { sum = 0.; gError = 0.0; nCounter = 0; for(Int_t iBinX = binMin; iBinX <= binMax; iBinX++) { //for(Int_t iBinX = 1; iBinX <= gHistNN->GetNbinsX(); iBinX++) { sum += gHistNN->GetBinContent(iBinX,iBinY); if(gHistNN->GetBinContent(iBinX,iBinY) != 0.) nCounter += 1; Double_t exy = gHistNN->GetCellError(iBinX,iBinY); gError += exy*exy; } if(nCounter != 0) { sum /= nCounter; gError = TMath::Sqrt(gError)/nCounter; } gHistNNprojection->SetBinContent(iBinY,sum); gHistNNprojection->SetBinError(iBinY,gError); } if(kUseZYAM) gHistNNprojection->GetYaxis()->SetTitle("#frac{1}{N_{trig}}#frac{dN_{assoc}}{#Delta#varphi} - b_{ZYAM} (rad^{-1})"); else gHistNNprojection->GetYaxis()->SetTitle("C_{--}(#Delta#varphi)"); gHistNNprojection->GetXaxis()->SetTitle("#Delta#varphi (rad)"); } //ZYAM if(kUseZYAM) { Double_t reference = gHistNNprojection->GetBinContent(gHistNNprojection->GetMinimumBin()); for(Int_t iBinX = 1; iBinX <= gHistNNprojection->GetNbinsX(); iBinX++) gHistNNprojection->SetBinContent(iBinX,gHistNNprojection->GetBinContent(iBinX) - reference); } gHistNNprojection->GetYaxis()->SetTitleOffset(1.5); gHistNNprojection->SetMarkerStyle(20); gHistNNprojection->SetStats(kFALSE); gHistNNprojection->DrawCopy("E"); //=================// gPad->SetTheta(30); // default is 30 gPad->SetPhi(-60); // default is 30 gPad->Update(); latexInfo1->DrawLatex(0.6,0.95,centralityLatex.Data()); latexInfo1->DrawLatex(0.6,0.89,psiLatex.Data()); latexInfo1->DrawLatex(0.6,0.83,pttLatex.Data()); latexInfo1->DrawLatex(0.6,0.77,ptaLatex.Data()); pngName = filename; if(kProjectInEta) pngName.ReplaceAll(".root","_DeltaEtaProjection"); else pngName.ReplaceAll(".root","_DeltaPhiProjection"); pngName += ".NegativeNegative.png"; cNN->SaveAs(pngName.Data()); TString outFileName = filename; if(kProjectInEta) outFileName.ReplaceAll(".root","_DeltaEtaProjection.root"); else outFileName.ReplaceAll(".root","_DeltaPhiProjection.root"); TFile *fProjection = TFile::Open(outFileName.Data(),"recreate"); gHistNPprojection->Write(); gHistPNprojection->Write(); gHistNNprojection->Write(); gHistPPprojection->Write(); fProjection->Close(); } //____________________________________________________________// void fitCorrelationFunctions(Int_t gCentrality = 1, Double_t psiMin = -0.5, Double_t psiMax = 3.5, Double_t vertexZMin = -10.,Double_t vertexZMax = 10., Double_t ptTriggerMin = -1., Double_t ptTriggerMax = -1., Double_t ptAssociatedMin = -1., Double_t ptAssociatedMax = -1., TH2D *gHist) { cout<<"FITTING FUNCTION"<<endl; //near side peak: [1]*TMath::Exp(-TMath::Power((0.5*TMath::Power((x/[2]),2)+0.5*TMath::Power((y/[3]),2)),[4])) //away side ridge: [5]*TMath::Exp(-TMath::Power((0.5*TMath::Power(((y-TMath::Pi())/[6]),2)),[7])) //longitudinal ridge: [8]*TMath::Exp(-TMath::Power((0.5*TMath::Power((x/[9]),2)),[10])) //wing structures: [11]*TMath::Power(x,2) //flow contribution (v1 up to v4): 2.*([12]*TMath::Cos(y) + [13]*TMath::Cos(2.*y) + [14]*TMath::Cos(3.*y) + [15]*TMath::Cos(4.*y)) TF2 *gFitFunction = new TF2("gFitFunction","[0]+[1]*TMath::Exp(-TMath::Power((0.5*TMath::Power((x/[2]),2)+0.5*TMath::Power((y/[3]),2)),[4]))+[5]*TMath::Exp(-TMath::Power((0.5*TMath::Power(((y-TMath::Pi())/[6]),2)),[7]))+[8]*TMath::Exp(-TMath::Power((0.5*TMath::Power(((x+[17])/[9]),2)),[10]))+[8]*TMath::Exp(-TMath::Power((0.5*TMath::Power(((x-[17])/[9]),2)),[10]))+[11]*TMath::Power(x,2)+2.*[12]*([13]*TMath::Cos(y) + [14]*TMath::Cos(2.*y) + [15]*TMath::Cos(3.*y) + [16]*TMath::Cos(4.*y))",-2.0,2.0,-TMath::Pi()/2.,3.*TMath::Pi()/2.); gFitFunction->SetName("gFitFunction"); //Normalization gFitFunction->SetParName(0,"N1"); //near side peak gFitFunction->SetParName(1,"N_{near side}");gFitFunction->FixParameter(1,0.0); gFitFunction->SetParName(2,"Sigma_{near side}(delta eta)"); gFitFunction->FixParameter(2,0.0); gFitFunction->SetParName(3,"Sigma_{near side}(delta phi)"); gFitFunction->FixParameter(3,0.0); gFitFunction->SetParName(4,"Exponent_{near side}"); gFitFunction->FixParameter(4,1.0); //away side ridge gFitFunction->SetParName(5,"N_{away side}"); gFitFunction->FixParameter(5,0.0); gFitFunction->SetParName(6,"Sigma_{away side}(delta phi)"); gFitFunction->FixParameter(6,0.0); gFitFunction->SetParName(7,"Exponent_{away side}"); gFitFunction->FixParameter(7,1.0); //longitudinal ridge gFitFunction->SetParName(8,"N_{long. ridge}"); gFitFunction->SetParameter(8,0.05);// gFitFunction->FixParameter(8,0.0); gFitFunction->SetParName(9,"Sigma_{long. ridge}(delta eta)"); gFitFunction->FixParameter(9,0.0); gFitFunction->SetParName(10,"Exponent_{long. ridge}"); gFitFunction->FixParameter(10,1.0); //wing structures gFitFunction->SetParName(11,"N_{wing}"); gFitFunction->FixParameter(11,0.0); //flow contribution gFitFunction->SetParName(12,"N_{flow}"); gFitFunction->SetParameter(12,0.2);gFitFunction->SetParLimits(12,0.0,10); gFitFunction->SetParName(13,"V1"); gFitFunction->SetParameter(13,0.005);gFitFunction->SetParLimits(13,0.0,10); gFitFunction->SetParName(14,"V2"); gFitFunction->SetParameter(14,0.1);gFitFunction->SetParLimits(14,0.0,10); gFitFunction->SetParName(15,"V3"); gFitFunction->SetParameter(15,0.05);gFitFunction->SetParLimits(15,0.0,10); gFitFunction->SetParName(16,"V4"); gFitFunction->SetParameter(16,0.005);gFitFunction->SetParLimits(16,0.0,10); gFitFunction->SetParName(17,"Offset"); gFitFunction->FixParameter(17,0.0); // flow parameters Double_t fNV = 0.; Double_t fV1 = 0.; Double_t fV2 = 0.; Double_t fV3 = 0.; Double_t fV4 = 0.; //Fitting the correlation function (first the edges to extract flow) gHist->Fit("gFitFunction","nm","",1.0,1.6); fNV += gFitFunction->GetParameter(12); fV1 += gFitFunction->GetParameter(13); fV2 += gFitFunction->GetParameter(14); fV3 += gFitFunction->GetParameter(15); fV4 += gFitFunction->GetParameter(16); gHist->Fit("gFitFunction","nm","",-1.6,-1.0); fNV += gFitFunction->GetParameter(12); fV1 += gFitFunction->GetParameter(13); fV2 += gFitFunction->GetParameter(14); fV3 += gFitFunction->GetParameter(15); fV4 += gFitFunction->GetParameter(16); // Now fit the whole with fixed flow gFitFunction->FixParameter(12,fNV/2.); gFitFunction->FixParameter(13,fV1/2.); gFitFunction->FixParameter(14,fV2/2.); gFitFunction->FixParameter(15,fV3/2.); gFitFunction->FixParameter(16,fV4/2.); gFitFunction->ReleaseParameter(0);gFitFunction->SetParameter(0,1.0); gFitFunction->ReleaseParameter(1);gFitFunction->SetParameter(1,0.3); gFitFunction->ReleaseParameter(2);gFitFunction->SetParameter(2,0.3);gFitFunction->SetParLimits(2,0.05,0.7); gFitFunction->ReleaseParameter(3);gFitFunction->SetParameter(3,0.3);gFitFunction->SetParLimits(3,0.05,1.7); //gFitFunction->ReleaseParameter(4);gFitFunction->SetParameter(4,1.0);gFitFunction->SetParLimits(4,0.0,2.0); gFitFunction->ReleaseParameter(5);gFitFunction->SetParameter(5,1.0);//gFitFunction->SetParLimits(5,0.0,10); gFitFunction->ReleaseParameter(6);gFitFunction->SetParameter(6,0.5);//gFitFunction->SetParLimits(6,0.0,10); //gFitFunction->ReleaseParameter(7);gFitFunction->SetParameter(7,1.0);gFitFunction->SetParLimits(7,0.0,2.0); //gFitFunction->ReleaseParameter(8);gFitFunction->SetParameter(8,0.05); //gFitFunction->ReleaseParameter(9);gFitFunction->SetParameter(9,0.6);gFitFunction->SetParLimits(9,0.1,10.0); //gFitFunction->ReleaseParameter(10);gFitFunction->SetParameter(10,1.0);gFitFunction->SetParLimits(10,0.0,2.0); gFitFunction->ReleaseParameter(17);gFitFunction->SetParameter(17,0.7);gFitFunction->SetParLimits(17,0.0,0.9); gHist->Fit("gFitFunction","nm"); //Cloning the histogram TString histoName = gHist->GetName(); histoName += "Fit"; TH2D *gHistFit = new TH2D(histoName.Data(),";#Delta#eta;#Delta#varphi (rad);C(#Delta#eta,#Delta#varphi)",gHist->GetNbinsX(),gHist->GetXaxis()->GetXmin(),gHist->GetXaxis()->GetXmax(),gHist->GetNbinsY(),gHist->GetYaxis()->GetXmin(),gHist->GetYaxis()->GetXmax()); TH2D *gHistResidual = dynamic_cast<TH2D *>(gHist->Clone()); gHistResidual->SetName("gHistResidual"); gHistResidual->Sumw2(); for(Int_t iBinDeltaEta = 1; iBinDeltaEta <= gHist->GetNbinsX(); iBinDeltaEta++) { for(Int_t iBinDeltaPhi = 1; iBinDeltaPhi <= gHist->GetNbinsY(); iBinDeltaPhi++) { gHistFit->SetBinContent(iBinDeltaEta,iBinDeltaPhi,gFitFunction->Eval(gHist->GetXaxis()->GetBinCenter(iBinDeltaEta),gHist->GetYaxis()->GetBinCenter(iBinDeltaPhi))); } } gHistResidual->Add(gHistFit,-1); //Write to output file TString newFileName = "correlationFunctionFit"; if(histoName.Contains("PN")) newFileName += "PN"; else if(histoName.Contains("NP")) newFileName += "NP"; else if(histoName.Contains("PP")) newFileName += "PP"; else if(histoName.Contains("NN")) newFileName += "NN"; newFileName += ".Centrality"; newFileName += gCentrality; newFileName += ".Psi"; if((psiMin == -0.5)&&(psiMax == 0.5)) newFileName += "InPlane.Ptt"; else if((psiMin == 0.5)&&(psiMax == 1.5)) newFileName += "Intermediate.Ptt"; else if((psiMin == 1.5)&&(psiMax == 2.5)) newFileName += "OutOfPlane.Ptt"; else if((psiMin == 2.5)&&(psiMax == 3.5)) newFileName += "Rest.PttFrom"; else newFileName += "All.PttFrom"; newFileName += Form("%.1f",ptTriggerMin); newFileName += "To"; newFileName += Form("%.1f",ptTriggerMax); newFileName += "PtaFrom"; newFileName += Form("%.1f",ptAssociatedMin); newFileName += "To"; newFileName += Form("%.1f",ptAssociatedMax); newFileName += "_"; newFileName += Form("%.1f",psiMin); newFileName += "-"; newFileName += Form("%.1f",psiMax); newFileName += ".root"; TFile *newFile = TFile::Open(newFileName.Data(),"recreate"); gHist->Write(); gHistFit->Write(); gHistResidual->Write(); gFitFunction->Write(); newFile->Close(); } //____________________________________________________________// TH2D *convolute2D(TH2D* h1, TH2D* h2, TString hname) { // // this function does the convolution of 2 eta or phi "efficiencies" in a deta or dphi "efficiency" // and returns a new histogram which is normalized to the number of bin combinations // histogram definition TH2D *hConv = NULL; hConv = new TH2D(hname.Data(),hname.Data(), h1->GetNbinsY(),-2,2,h1->GetNbinsX(),-TMath::Pi()/2.,3*TMath::Pi()/2.); Double_t x1 = 0.; Double_t x2 = 0.; Double_t y1 = 0.; Double_t y2 = 0.; Double_t z1 = 0.; Double_t z2 = 0.; Double_t dx = 0.; Double_t dy = 0.; Double_t dz = 0.; // convolution for(Int_t i = 0; i < h1->GetNbinsX(); i ++){ cout<<"CONVOLUTING BIN = "<<i<<" of "<<h1->GetNbinsX()<<endl; for(Int_t j = 0; j < h1->GetNbinsY(); j ++){ for(Int_t k = 0; k < h2->GetNbinsX(); k ++){ for(Int_t l = 0; l < h2->GetNbinsY(); l ++){ x1 = (Double_t)h1->GetXaxis()->GetBinCenter(i+1); y1 = (Double_t)h1->GetYaxis()->GetBinCenter(j+1); x2 = (Double_t)h2->GetXaxis()->GetBinCenter(k+1); y2 = (Double_t)h2->GetYaxis()->GetBinCenter(l+1); z1 = (Double_t)h1->GetBinContent(i+1,j+1); z2 = (Double_t)h2->GetBinContent(k+1,l+1); // need the gymnastics to keep the same binning dx = x1 - x2 - (h1->GetXaxis()->GetBinWidth(1)/2.); if(gRandom->Gaus() > 0) dy = y1 - y2 + (h1->GetYaxis()->GetBinWidth(1)/2.); else dy = y1 - y2 - (h1->GetYaxis()->GetBinWidth(1)/2.); if(dx>3./2.*TMath::Pi()) dx = dx - 2.*TMath::Pi(); if(dx<-1./2.*TMath::Pi()) dx = 2*TMath::Pi() + dx; dz = z1*z2; hConv->Fill(dy,dx,dz); } } } } return hConv; }
drawCorrelationFunctionPsi.C:1
drawCorrelationFunctionPsi.C:2
drawCorrelationFunctionPsi.C:3
drawCorrelationFunctionPsi.C:4
drawCorrelationFunctionPsi.C:5
drawCorrelationFunctionPsi.C:6
drawCorrelationFunctionPsi.C:7
drawCorrelationFunctionPsi.C:8
drawCorrelationFunctionPsi.C:9
drawCorrelationFunctionPsi.C:10
drawCorrelationFunctionPsi.C:11
drawCorrelationFunctionPsi.C:12
drawCorrelationFunctionPsi.C:13
drawCorrelationFunctionPsi.C:14
drawCorrelationFunctionPsi.C:15
drawCorrelationFunctionPsi.C:16
drawCorrelationFunctionPsi.C:17
drawCorrelationFunctionPsi.C:18
drawCorrelationFunctionPsi.C:19
drawCorrelationFunctionPsi.C:20
drawCorrelationFunctionPsi.C:21
drawCorrelationFunctionPsi.C:22
drawCorrelationFunctionPsi.C:23
drawCorrelationFunctionPsi.C:24
drawCorrelationFunctionPsi.C:25
drawCorrelationFunctionPsi.C:26
drawCorrelationFunctionPsi.C:27
drawCorrelationFunctionPsi.C:28
drawCorrelationFunctionPsi.C:29
drawCorrelationFunctionPsi.C:30
drawCorrelationFunctionPsi.C:31
drawCorrelationFunctionPsi.C:32
drawCorrelationFunctionPsi.C:33
drawCorrelationFunctionPsi.C:34
drawCorrelationFunctionPsi.C:35
drawCorrelationFunctionPsi.C:36
drawCorrelationFunctionPsi.C:37
drawCorrelationFunctionPsi.C:38
drawCorrelationFunctionPsi.C:39
drawCorrelationFunctionPsi.C:40
drawCorrelationFunctionPsi.C:41
drawCorrelationFunctionPsi.C:42
drawCorrelationFunctionPsi.C:43
drawCorrelationFunctionPsi.C:44
drawCorrelationFunctionPsi.C:45
drawCorrelationFunctionPsi.C:46
drawCorrelationFunctionPsi.C:47
drawCorrelationFunctionPsi.C:48
drawCorrelationFunctionPsi.C:49
drawCorrelationFunctionPsi.C:50
drawCorrelationFunctionPsi.C:51
drawCorrelationFunctionPsi.C:52
drawCorrelationFunctionPsi.C:53
drawCorrelationFunctionPsi.C:54
drawCorrelationFunctionPsi.C:55
drawCorrelationFunctionPsi.C:56
drawCorrelationFunctionPsi.C:57
drawCorrelationFunctionPsi.C:58
drawCorrelationFunctionPsi.C:59
drawCorrelationFunctionPsi.C:60
drawCorrelationFunctionPsi.C:61
drawCorrelationFunctionPsi.C:62
drawCorrelationFunctionPsi.C:63
drawCorrelationFunctionPsi.C:64
drawCorrelationFunctionPsi.C:65
drawCorrelationFunctionPsi.C:66
drawCorrelationFunctionPsi.C:67
drawCorrelationFunctionPsi.C:68
drawCorrelationFunctionPsi.C:69
drawCorrelationFunctionPsi.C:70
drawCorrelationFunctionPsi.C:71
drawCorrelationFunctionPsi.C:72
drawCorrelationFunctionPsi.C:73
drawCorrelationFunctionPsi.C:74
drawCorrelationFunctionPsi.C:75
drawCorrelationFunctionPsi.C:76
drawCorrelationFunctionPsi.C:77
drawCorrelationFunctionPsi.C:78
drawCorrelationFunctionPsi.C:79
drawCorrelationFunctionPsi.C:80
drawCorrelationFunctionPsi.C:81
drawCorrelationFunctionPsi.C:82
drawCorrelationFunctionPsi.C:83
drawCorrelationFunctionPsi.C:84
drawCorrelationFunctionPsi.C:85
drawCorrelationFunctionPsi.C:86
drawCorrelationFunctionPsi.C:87
drawCorrelationFunctionPsi.C:88
drawCorrelationFunctionPsi.C:89
drawCorrelationFunctionPsi.C:90
drawCorrelationFunctionPsi.C:91
drawCorrelationFunctionPsi.C:92
drawCorrelationFunctionPsi.C:93
drawCorrelationFunctionPsi.C:94
drawCorrelationFunctionPsi.C:95
drawCorrelationFunctionPsi.C:96
drawCorrelationFunctionPsi.C:97
drawCorrelationFunctionPsi.C:98
drawCorrelationFunctionPsi.C:99
drawCorrelationFunctionPsi.C:100
drawCorrelationFunctionPsi.C:101
drawCorrelationFunctionPsi.C:102
drawCorrelationFunctionPsi.C:103
drawCorrelationFunctionPsi.C:104
drawCorrelationFunctionPsi.C:105
drawCorrelationFunctionPsi.C:106
drawCorrelationFunctionPsi.C:107
drawCorrelationFunctionPsi.C:108
drawCorrelationFunctionPsi.C:109
drawCorrelationFunctionPsi.C:110
drawCorrelationFunctionPsi.C:111
drawCorrelationFunctionPsi.C:112
drawCorrelationFunctionPsi.C:113
drawCorrelationFunctionPsi.C:114
drawCorrelationFunctionPsi.C:115
drawCorrelationFunctionPsi.C:116
drawCorrelationFunctionPsi.C:117
drawCorrelationFunctionPsi.C:118
drawCorrelationFunctionPsi.C:119
drawCorrelationFunctionPsi.C:120
drawCorrelationFunctionPsi.C:121
drawCorrelationFunctionPsi.C:122
drawCorrelationFunctionPsi.C:123
drawCorrelationFunctionPsi.C:124
drawCorrelationFunctionPsi.C:125
drawCorrelationFunctionPsi.C:126
drawCorrelationFunctionPsi.C:127
drawCorrelationFunctionPsi.C:128
drawCorrelationFunctionPsi.C:129
drawCorrelationFunctionPsi.C:130
drawCorrelationFunctionPsi.C:131
drawCorrelationFunctionPsi.C:132
drawCorrelationFunctionPsi.C:133
drawCorrelationFunctionPsi.C:134
drawCorrelationFunctionPsi.C:135
drawCorrelationFunctionPsi.C:136
drawCorrelationFunctionPsi.C:137
drawCorrelationFunctionPsi.C:138
drawCorrelationFunctionPsi.C:139
drawCorrelationFunctionPsi.C:140
drawCorrelationFunctionPsi.C:141
drawCorrelationFunctionPsi.C:142
drawCorrelationFunctionPsi.C:143
drawCorrelationFunctionPsi.C:144
drawCorrelationFunctionPsi.C:145
drawCorrelationFunctionPsi.C:146
drawCorrelationFunctionPsi.C:147
drawCorrelationFunctionPsi.C:148
drawCorrelationFunctionPsi.C:149
drawCorrelationFunctionPsi.C:150
drawCorrelationFunctionPsi.C:151
drawCorrelationFunctionPsi.C:152
drawCorrelationFunctionPsi.C:153
drawCorrelationFunctionPsi.C:154
drawCorrelationFunctionPsi.C:155
drawCorrelationFunctionPsi.C:156
drawCorrelationFunctionPsi.C:157
drawCorrelationFunctionPsi.C:158
drawCorrelationFunctionPsi.C:159
drawCorrelationFunctionPsi.C:160
drawCorrelationFunctionPsi.C:161
drawCorrelationFunctionPsi.C:162
drawCorrelationFunctionPsi.C:163
drawCorrelationFunctionPsi.C:164
drawCorrelationFunctionPsi.C:165
drawCorrelationFunctionPsi.C:166
drawCorrelationFunctionPsi.C:167
drawCorrelationFunctionPsi.C:168
drawCorrelationFunctionPsi.C:169
drawCorrelationFunctionPsi.C:170
drawCorrelationFunctionPsi.C:171
drawCorrelationFunctionPsi.C:172
drawCorrelationFunctionPsi.C:173
drawCorrelationFunctionPsi.C:174
drawCorrelationFunctionPsi.C:175
drawCorrelationFunctionPsi.C:176
drawCorrelationFunctionPsi.C:177
drawCorrelationFunctionPsi.C:178
drawCorrelationFunctionPsi.C:179
drawCorrelationFunctionPsi.C:180
drawCorrelationFunctionPsi.C:181
drawCorrelationFunctionPsi.C:182
drawCorrelationFunctionPsi.C:183
drawCorrelationFunctionPsi.C:184
drawCorrelationFunctionPsi.C:185
drawCorrelationFunctionPsi.C:186
drawCorrelationFunctionPsi.C:187
drawCorrelationFunctionPsi.C:188
drawCorrelationFunctionPsi.C:189
drawCorrelationFunctionPsi.C:190
drawCorrelationFunctionPsi.C:191
drawCorrelationFunctionPsi.C:192
drawCorrelationFunctionPsi.C:193
drawCorrelationFunctionPsi.C:194
drawCorrelationFunctionPsi.C:195
drawCorrelationFunctionPsi.C:196
drawCorrelationFunctionPsi.C:197
drawCorrelationFunctionPsi.C:198
drawCorrelationFunctionPsi.C:199
drawCorrelationFunctionPsi.C:200
drawCorrelationFunctionPsi.C:201
drawCorrelationFunctionPsi.C:202
drawCorrelationFunctionPsi.C:203
drawCorrelationFunctionPsi.C:204
drawCorrelationFunctionPsi.C:205
drawCorrelationFunctionPsi.C:206
drawCorrelationFunctionPsi.C:207
drawCorrelationFunctionPsi.C:208
drawCorrelationFunctionPsi.C:209
drawCorrelationFunctionPsi.C:210
drawCorrelationFunctionPsi.C:211
drawCorrelationFunctionPsi.C:212
drawCorrelationFunctionPsi.C:213
drawCorrelationFunctionPsi.C:214
drawCorrelationFunctionPsi.C:215
drawCorrelationFunctionPsi.C:216
drawCorrelationFunctionPsi.C:217
drawCorrelationFunctionPsi.C:218
drawCorrelationFunctionPsi.C:219
drawCorrelationFunctionPsi.C:220
drawCorrelationFunctionPsi.C:221
drawCorrelationFunctionPsi.C:222
drawCorrelationFunctionPsi.C:223
drawCorrelationFunctionPsi.C:224
drawCorrelationFunctionPsi.C:225
drawCorrelationFunctionPsi.C:226
drawCorrelationFunctionPsi.C:227
drawCorrelationFunctionPsi.C:228
drawCorrelationFunctionPsi.C:229
drawCorrelationFunctionPsi.C:230
drawCorrelationFunctionPsi.C:231
drawCorrelationFunctionPsi.C:232
drawCorrelationFunctionPsi.C:233
drawCorrelationFunctionPsi.C:234
drawCorrelationFunctionPsi.C:235
drawCorrelationFunctionPsi.C:236
drawCorrelationFunctionPsi.C:237
drawCorrelationFunctionPsi.C:238
drawCorrelationFunctionPsi.C:239
drawCorrelationFunctionPsi.C:240
drawCorrelationFunctionPsi.C:241
drawCorrelationFunctionPsi.C:242
drawCorrelationFunctionPsi.C:243
drawCorrelationFunctionPsi.C:244
drawCorrelationFunctionPsi.C:245
drawCorrelationFunctionPsi.C:246
drawCorrelationFunctionPsi.C:247
drawCorrelationFunctionPsi.C:248
drawCorrelationFunctionPsi.C:249
drawCorrelationFunctionPsi.C:250
drawCorrelationFunctionPsi.C:251
drawCorrelationFunctionPsi.C:252
drawCorrelationFunctionPsi.C:253
drawCorrelationFunctionPsi.C:254
drawCorrelationFunctionPsi.C:255
drawCorrelationFunctionPsi.C:256
drawCorrelationFunctionPsi.C:257
drawCorrelationFunctionPsi.C:258
drawCorrelationFunctionPsi.C:259
drawCorrelationFunctionPsi.C:260
drawCorrelationFunctionPsi.C:261
drawCorrelationFunctionPsi.C:262
drawCorrelationFunctionPsi.C:263
drawCorrelationFunctionPsi.C:264
drawCorrelationFunctionPsi.C:265
drawCorrelationFunctionPsi.C:266
drawCorrelationFunctionPsi.C:267
drawCorrelationFunctionPsi.C:268
drawCorrelationFunctionPsi.C:269
drawCorrelationFunctionPsi.C:270
drawCorrelationFunctionPsi.C:271
drawCorrelationFunctionPsi.C:272
drawCorrelationFunctionPsi.C:273
drawCorrelationFunctionPsi.C:274
drawCorrelationFunctionPsi.C:275
drawCorrelationFunctionPsi.C:276
drawCorrelationFunctionPsi.C:277
drawCorrelationFunctionPsi.C:278
drawCorrelationFunctionPsi.C:279
drawCorrelationFunctionPsi.C:280
drawCorrelationFunctionPsi.C:281
drawCorrelationFunctionPsi.C:282
drawCorrelationFunctionPsi.C:283
drawCorrelationFunctionPsi.C:284
drawCorrelationFunctionPsi.C:285
drawCorrelationFunctionPsi.C:286
drawCorrelationFunctionPsi.C:287
drawCorrelationFunctionPsi.C:288
drawCorrelationFunctionPsi.C:289
drawCorrelationFunctionPsi.C:290
drawCorrelationFunctionPsi.C:291
drawCorrelationFunctionPsi.C:292
drawCorrelationFunctionPsi.C:293
drawCorrelationFunctionPsi.C:294
drawCorrelationFunctionPsi.C:295
drawCorrelationFunctionPsi.C:296
drawCorrelationFunctionPsi.C:297
drawCorrelationFunctionPsi.C:298
drawCorrelationFunctionPsi.C:299
drawCorrelationFunctionPsi.C:300
drawCorrelationFunctionPsi.C:301
drawCorrelationFunctionPsi.C:302
drawCorrelationFunctionPsi.C:303
drawCorrelationFunctionPsi.C:304
drawCorrelationFunctionPsi.C:305
drawCorrelationFunctionPsi.C:306
drawCorrelationFunctionPsi.C:307
drawCorrelationFunctionPsi.C:308
drawCorrelationFunctionPsi.C:309
drawCorrelationFunctionPsi.C:310
drawCorrelationFunctionPsi.C:311
drawCorrelationFunctionPsi.C:312
drawCorrelationFunctionPsi.C:313
drawCorrelationFunctionPsi.C:314
drawCorrelationFunctionPsi.C:315
drawCorrelationFunctionPsi.C:316
drawCorrelationFunctionPsi.C:317
drawCorrelationFunctionPsi.C:318
drawCorrelationFunctionPsi.C:319
drawCorrelationFunctionPsi.C:320
drawCorrelationFunctionPsi.C:321
drawCorrelationFunctionPsi.C:322
drawCorrelationFunctionPsi.C:323
drawCorrelationFunctionPsi.C:324
drawCorrelationFunctionPsi.C:325
drawCorrelationFunctionPsi.C:326
drawCorrelationFunctionPsi.C:327
drawCorrelationFunctionPsi.C:328
drawCorrelationFunctionPsi.C:329
drawCorrelationFunctionPsi.C:330
drawCorrelationFunctionPsi.C:331
drawCorrelationFunctionPsi.C:332
drawCorrelationFunctionPsi.C:333
drawCorrelationFunctionPsi.C:334
drawCorrelationFunctionPsi.C:335
drawCorrelationFunctionPsi.C:336
drawCorrelationFunctionPsi.C:337
drawCorrelationFunctionPsi.C:338
drawCorrelationFunctionPsi.C:339
drawCorrelationFunctionPsi.C:340
drawCorrelationFunctionPsi.C:341
drawCorrelationFunctionPsi.C:342
drawCorrelationFunctionPsi.C:343
drawCorrelationFunctionPsi.C:344
drawCorrelationFunctionPsi.C:345
drawCorrelationFunctionPsi.C:346
drawCorrelationFunctionPsi.C:347
drawCorrelationFunctionPsi.C:348
drawCorrelationFunctionPsi.C:349
drawCorrelationFunctionPsi.C:350
drawCorrelationFunctionPsi.C:351
drawCorrelationFunctionPsi.C:352
drawCorrelationFunctionPsi.C:353
drawCorrelationFunctionPsi.C:354
drawCorrelationFunctionPsi.C:355
drawCorrelationFunctionPsi.C:356
drawCorrelationFunctionPsi.C:357
drawCorrelationFunctionPsi.C:358
drawCorrelationFunctionPsi.C:359
drawCorrelationFunctionPsi.C:360
drawCorrelationFunctionPsi.C:361
drawCorrelationFunctionPsi.C:362
drawCorrelationFunctionPsi.C:363
drawCorrelationFunctionPsi.C:364
drawCorrelationFunctionPsi.C:365
drawCorrelationFunctionPsi.C:366
drawCorrelationFunctionPsi.C:367
drawCorrelationFunctionPsi.C:368
drawCorrelationFunctionPsi.C:369
drawCorrelationFunctionPsi.C:370
drawCorrelationFunctionPsi.C:371
drawCorrelationFunctionPsi.C:372
drawCorrelationFunctionPsi.C:373
drawCorrelationFunctionPsi.C:374
drawCorrelationFunctionPsi.C:375
drawCorrelationFunctionPsi.C:376
drawCorrelationFunctionPsi.C:377
drawCorrelationFunctionPsi.C:378
drawCorrelationFunctionPsi.C:379
drawCorrelationFunctionPsi.C:380
drawCorrelationFunctionPsi.C:381
drawCorrelationFunctionPsi.C:382
drawCorrelationFunctionPsi.C:383
drawCorrelationFunctionPsi.C:384
drawCorrelationFunctionPsi.C:385
drawCorrelationFunctionPsi.C:386
drawCorrelationFunctionPsi.C:387
drawCorrelationFunctionPsi.C:388
drawCorrelationFunctionPsi.C:389
drawCorrelationFunctionPsi.C:390
drawCorrelationFunctionPsi.C:391
drawCorrelationFunctionPsi.C:392
drawCorrelationFunctionPsi.C:393
drawCorrelationFunctionPsi.C:394
drawCorrelationFunctionPsi.C:395
drawCorrelationFunctionPsi.C:396
drawCorrelationFunctionPsi.C:397
drawCorrelationFunctionPsi.C:398
drawCorrelationFunctionPsi.C:399
drawCorrelationFunctionPsi.C:400
drawCorrelationFunctionPsi.C:401
drawCorrelationFunctionPsi.C:402
drawCorrelationFunctionPsi.C:403
drawCorrelationFunctionPsi.C:404
drawCorrelationFunctionPsi.C:405
drawCorrelationFunctionPsi.C:406
drawCorrelationFunctionPsi.C:407
drawCorrelationFunctionPsi.C:408
drawCorrelationFunctionPsi.C:409
drawCorrelationFunctionPsi.C:410
drawCorrelationFunctionPsi.C:411
drawCorrelationFunctionPsi.C:412
drawCorrelationFunctionPsi.C:413
drawCorrelationFunctionPsi.C:414
drawCorrelationFunctionPsi.C:415
drawCorrelationFunctionPsi.C:416
drawCorrelationFunctionPsi.C:417
drawCorrelationFunctionPsi.C:418
drawCorrelationFunctionPsi.C:419
drawCorrelationFunctionPsi.C:420
drawCorrelationFunctionPsi.C:421
drawCorrelationFunctionPsi.C:422
drawCorrelationFunctionPsi.C:423
drawCorrelationFunctionPsi.C:424
drawCorrelationFunctionPsi.C:425
drawCorrelationFunctionPsi.C:426
drawCorrelationFunctionPsi.C:427
drawCorrelationFunctionPsi.C:428
drawCorrelationFunctionPsi.C:429
drawCorrelationFunctionPsi.C:430
drawCorrelationFunctionPsi.C:431
drawCorrelationFunctionPsi.C:432
drawCorrelationFunctionPsi.C:433
drawCorrelationFunctionPsi.C:434
drawCorrelationFunctionPsi.C:435
drawCorrelationFunctionPsi.C:436
drawCorrelationFunctionPsi.C:437
drawCorrelationFunctionPsi.C:438
drawCorrelationFunctionPsi.C:439
drawCorrelationFunctionPsi.C:440
drawCorrelationFunctionPsi.C:441
drawCorrelationFunctionPsi.C:442
drawCorrelationFunctionPsi.C:443
drawCorrelationFunctionPsi.C:444
drawCorrelationFunctionPsi.C:445
drawCorrelationFunctionPsi.C:446
drawCorrelationFunctionPsi.C:447
drawCorrelationFunctionPsi.C:448
drawCorrelationFunctionPsi.C:449
drawCorrelationFunctionPsi.C:450
drawCorrelationFunctionPsi.C:451
drawCorrelationFunctionPsi.C:452
drawCorrelationFunctionPsi.C:453
drawCorrelationFunctionPsi.C:454
drawCorrelationFunctionPsi.C:455
drawCorrelationFunctionPsi.C:456
drawCorrelationFunctionPsi.C:457
drawCorrelationFunctionPsi.C:458
drawCorrelationFunctionPsi.C:459
drawCorrelationFunctionPsi.C:460
drawCorrelationFunctionPsi.C:461
drawCorrelationFunctionPsi.C:462
drawCorrelationFunctionPsi.C:463
drawCorrelationFunctionPsi.C:464
drawCorrelationFunctionPsi.C:465
drawCorrelationFunctionPsi.C:466
drawCorrelationFunctionPsi.C:467
drawCorrelationFunctionPsi.C:468
drawCorrelationFunctionPsi.C:469
drawCorrelationFunctionPsi.C:470
drawCorrelationFunctionPsi.C:471
drawCorrelationFunctionPsi.C:472
drawCorrelationFunctionPsi.C:473
drawCorrelationFunctionPsi.C:474
drawCorrelationFunctionPsi.C:475
drawCorrelationFunctionPsi.C:476
drawCorrelationFunctionPsi.C:477
drawCorrelationFunctionPsi.C:478
drawCorrelationFunctionPsi.C:479
drawCorrelationFunctionPsi.C:480
drawCorrelationFunctionPsi.C:481
drawCorrelationFunctionPsi.C:482
drawCorrelationFunctionPsi.C:483
drawCorrelationFunctionPsi.C:484
drawCorrelationFunctionPsi.C:485
drawCorrelationFunctionPsi.C:486
drawCorrelationFunctionPsi.C:487
drawCorrelationFunctionPsi.C:488
drawCorrelationFunctionPsi.C:489
drawCorrelationFunctionPsi.C:490
drawCorrelationFunctionPsi.C:491
drawCorrelationFunctionPsi.C:492
drawCorrelationFunctionPsi.C:493
drawCorrelationFunctionPsi.C:494
drawCorrelationFunctionPsi.C:495
drawCorrelationFunctionPsi.C:496
drawCorrelationFunctionPsi.C:497
drawCorrelationFunctionPsi.C:498
drawCorrelationFunctionPsi.C:499
drawCorrelationFunctionPsi.C:500
drawCorrelationFunctionPsi.C:501
drawCorrelationFunctionPsi.C:502
drawCorrelationFunctionPsi.C:503
drawCorrelationFunctionPsi.C:504
drawCorrelationFunctionPsi.C:505
drawCorrelationFunctionPsi.C:506
drawCorrelationFunctionPsi.C:507
drawCorrelationFunctionPsi.C:508
drawCorrelationFunctionPsi.C:509
drawCorrelationFunctionPsi.C:510
drawCorrelationFunctionPsi.C:511
drawCorrelationFunctionPsi.C:512
drawCorrelationFunctionPsi.C:513
drawCorrelationFunctionPsi.C:514
drawCorrelationFunctionPsi.C:515
drawCorrelationFunctionPsi.C:516
drawCorrelationFunctionPsi.C:517
drawCorrelationFunctionPsi.C:518
drawCorrelationFunctionPsi.C:519
drawCorrelationFunctionPsi.C:520
drawCorrelationFunctionPsi.C:521
drawCorrelationFunctionPsi.C:522
drawCorrelationFunctionPsi.C:523
drawCorrelationFunctionPsi.C:524
drawCorrelationFunctionPsi.C:525
drawCorrelationFunctionPsi.C:526
drawCorrelationFunctionPsi.C:527
drawCorrelationFunctionPsi.C:528
drawCorrelationFunctionPsi.C:529
drawCorrelationFunctionPsi.C:530
drawCorrelationFunctionPsi.C:531
drawCorrelationFunctionPsi.C:532
drawCorrelationFunctionPsi.C:533
drawCorrelationFunctionPsi.C:534
drawCorrelationFunctionPsi.C:535
drawCorrelationFunctionPsi.C:536
drawCorrelationFunctionPsi.C:537
drawCorrelationFunctionPsi.C:538
drawCorrelationFunctionPsi.C:539
drawCorrelationFunctionPsi.C:540
drawCorrelationFunctionPsi.C:541
drawCorrelationFunctionPsi.C:542
drawCorrelationFunctionPsi.C:543
drawCorrelationFunctionPsi.C:544
drawCorrelationFunctionPsi.C:545
drawCorrelationFunctionPsi.C:546
drawCorrelationFunctionPsi.C:547
drawCorrelationFunctionPsi.C:548
drawCorrelationFunctionPsi.C:549
drawCorrelationFunctionPsi.C:550
drawCorrelationFunctionPsi.C:551
drawCorrelationFunctionPsi.C:552
drawCorrelationFunctionPsi.C:553
drawCorrelationFunctionPsi.C:554
drawCorrelationFunctionPsi.C:555
drawCorrelationFunctionPsi.C:556
drawCorrelationFunctionPsi.C:557
drawCorrelationFunctionPsi.C:558
drawCorrelationFunctionPsi.C:559
drawCorrelationFunctionPsi.C:560
drawCorrelationFunctionPsi.C:561
drawCorrelationFunctionPsi.C:562
drawCorrelationFunctionPsi.C:563
drawCorrelationFunctionPsi.C:564
drawCorrelationFunctionPsi.C:565
drawCorrelationFunctionPsi.C:566
drawCorrelationFunctionPsi.C:567
drawCorrelationFunctionPsi.C:568
drawCorrelationFunctionPsi.C:569
drawCorrelationFunctionPsi.C:570
drawCorrelationFunctionPsi.C:571
drawCorrelationFunctionPsi.C:572
drawCorrelationFunctionPsi.C:573
drawCorrelationFunctionPsi.C:574
drawCorrelationFunctionPsi.C:575
drawCorrelationFunctionPsi.C:576
drawCorrelationFunctionPsi.C:577
drawCorrelationFunctionPsi.C:578
drawCorrelationFunctionPsi.C:579
drawCorrelationFunctionPsi.C:580
drawCorrelationFunctionPsi.C:581
drawCorrelationFunctionPsi.C:582
drawCorrelationFunctionPsi.C:583
drawCorrelationFunctionPsi.C:584
drawCorrelationFunctionPsi.C:585
drawCorrelationFunctionPsi.C:586
drawCorrelationFunctionPsi.C:587
drawCorrelationFunctionPsi.C:588
drawCorrelationFunctionPsi.C:589
drawCorrelationFunctionPsi.C:590
drawCorrelationFunctionPsi.C:591
drawCorrelationFunctionPsi.C:592
drawCorrelationFunctionPsi.C:593
drawCorrelationFunctionPsi.C:594
drawCorrelationFunctionPsi.C:595
drawCorrelationFunctionPsi.C:596
drawCorrelationFunctionPsi.C:597
drawCorrelationFunctionPsi.C:598
drawCorrelationFunctionPsi.C:599
drawCorrelationFunctionPsi.C:600
drawCorrelationFunctionPsi.C:601
drawCorrelationFunctionPsi.C:602
drawCorrelationFunctionPsi.C:603
drawCorrelationFunctionPsi.C:604
drawCorrelationFunctionPsi.C:605
drawCorrelationFunctionPsi.C:606
drawCorrelationFunctionPsi.C:607
drawCorrelationFunctionPsi.C:608
drawCorrelationFunctionPsi.C:609
drawCorrelationFunctionPsi.C:610
drawCorrelationFunctionPsi.C:611
drawCorrelationFunctionPsi.C:612
drawCorrelationFunctionPsi.C:613
drawCorrelationFunctionPsi.C:614
drawCorrelationFunctionPsi.C:615
drawCorrelationFunctionPsi.C:616
drawCorrelationFunctionPsi.C:617
drawCorrelationFunctionPsi.C:618
drawCorrelationFunctionPsi.C:619
drawCorrelationFunctionPsi.C:620
drawCorrelationFunctionPsi.C:621
drawCorrelationFunctionPsi.C:622
drawCorrelationFunctionPsi.C:623
drawCorrelationFunctionPsi.C:624
drawCorrelationFunctionPsi.C:625
drawCorrelationFunctionPsi.C:626
drawCorrelationFunctionPsi.C:627
drawCorrelationFunctionPsi.C:628
drawCorrelationFunctionPsi.C:629
drawCorrelationFunctionPsi.C:630
drawCorrelationFunctionPsi.C:631
drawCorrelationFunctionPsi.C:632
drawCorrelationFunctionPsi.C:633
drawCorrelationFunctionPsi.C:634
drawCorrelationFunctionPsi.C:635
drawCorrelationFunctionPsi.C:636
drawCorrelationFunctionPsi.C:637
drawCorrelationFunctionPsi.C:638
drawCorrelationFunctionPsi.C:639
drawCorrelationFunctionPsi.C:640
drawCorrelationFunctionPsi.C:641
drawCorrelationFunctionPsi.C:642
drawCorrelationFunctionPsi.C:643
drawCorrelationFunctionPsi.C:644
drawCorrelationFunctionPsi.C:645
drawCorrelationFunctionPsi.C:646
drawCorrelationFunctionPsi.C:647
drawCorrelationFunctionPsi.C:648
drawCorrelationFunctionPsi.C:649
drawCorrelationFunctionPsi.C:650
drawCorrelationFunctionPsi.C:651
drawCorrelationFunctionPsi.C:652
drawCorrelationFunctionPsi.C:653
drawCorrelationFunctionPsi.C:654
drawCorrelationFunctionPsi.C:655
drawCorrelationFunctionPsi.C:656
drawCorrelationFunctionPsi.C:657
drawCorrelationFunctionPsi.C:658
drawCorrelationFunctionPsi.C:659
drawCorrelationFunctionPsi.C:660
drawCorrelationFunctionPsi.C:661
drawCorrelationFunctionPsi.C:662
drawCorrelationFunctionPsi.C:663
drawCorrelationFunctionPsi.C:664
drawCorrelationFunctionPsi.C:665
drawCorrelationFunctionPsi.C:666
drawCorrelationFunctionPsi.C:667
drawCorrelationFunctionPsi.C:668
drawCorrelationFunctionPsi.C:669
drawCorrelationFunctionPsi.C:670
drawCorrelationFunctionPsi.C:671
drawCorrelationFunctionPsi.C:672
drawCorrelationFunctionPsi.C:673
drawCorrelationFunctionPsi.C:674
drawCorrelationFunctionPsi.C:675
drawCorrelationFunctionPsi.C:676
drawCorrelationFunctionPsi.C:677
drawCorrelationFunctionPsi.C:678
drawCorrelationFunctionPsi.C:679
drawCorrelationFunctionPsi.C:680
drawCorrelationFunctionPsi.C:681
drawCorrelationFunctionPsi.C:682
drawCorrelationFunctionPsi.C:683
drawCorrelationFunctionPsi.C:684
drawCorrelationFunctionPsi.C:685
drawCorrelationFunctionPsi.C:686
drawCorrelationFunctionPsi.C:687
drawCorrelationFunctionPsi.C:688
drawCorrelationFunctionPsi.C:689
drawCorrelationFunctionPsi.C:690
drawCorrelationFunctionPsi.C:691
drawCorrelationFunctionPsi.C:692
drawCorrelationFunctionPsi.C:693
drawCorrelationFunctionPsi.C:694
drawCorrelationFunctionPsi.C:695
drawCorrelationFunctionPsi.C:696
drawCorrelationFunctionPsi.C:697
drawCorrelationFunctionPsi.C:698
drawCorrelationFunctionPsi.C:699
drawCorrelationFunctionPsi.C:700
drawCorrelationFunctionPsi.C:701
drawCorrelationFunctionPsi.C:702
drawCorrelationFunctionPsi.C:703
drawCorrelationFunctionPsi.C:704
drawCorrelationFunctionPsi.C:705
drawCorrelationFunctionPsi.C:706
drawCorrelationFunctionPsi.C:707
drawCorrelationFunctionPsi.C:708
drawCorrelationFunctionPsi.C:709
drawCorrelationFunctionPsi.C:710
drawCorrelationFunctionPsi.C:711
drawCorrelationFunctionPsi.C:712
drawCorrelationFunctionPsi.C:713
drawCorrelationFunctionPsi.C:714
drawCorrelationFunctionPsi.C:715
drawCorrelationFunctionPsi.C:716
drawCorrelationFunctionPsi.C:717
drawCorrelationFunctionPsi.C:718
drawCorrelationFunctionPsi.C:719
drawCorrelationFunctionPsi.C:720
drawCorrelationFunctionPsi.C:721
drawCorrelationFunctionPsi.C:722
drawCorrelationFunctionPsi.C:723
drawCorrelationFunctionPsi.C:724
drawCorrelationFunctionPsi.C:725
drawCorrelationFunctionPsi.C:726
drawCorrelationFunctionPsi.C:727
drawCorrelationFunctionPsi.C:728
drawCorrelationFunctionPsi.C:729
drawCorrelationFunctionPsi.C:730
drawCorrelationFunctionPsi.C:731
drawCorrelationFunctionPsi.C:732
drawCorrelationFunctionPsi.C:733
drawCorrelationFunctionPsi.C:734
drawCorrelationFunctionPsi.C:735
drawCorrelationFunctionPsi.C:736
drawCorrelationFunctionPsi.C:737
drawCorrelationFunctionPsi.C:738
drawCorrelationFunctionPsi.C:739
drawCorrelationFunctionPsi.C:740
drawCorrelationFunctionPsi.C:741
drawCorrelationFunctionPsi.C:742
drawCorrelationFunctionPsi.C:743
drawCorrelationFunctionPsi.C:744
drawCorrelationFunctionPsi.C:745
drawCorrelationFunctionPsi.C:746
drawCorrelationFunctionPsi.C:747
drawCorrelationFunctionPsi.C:748
drawCorrelationFunctionPsi.C:749
drawCorrelationFunctionPsi.C:750
drawCorrelationFunctionPsi.C:751
drawCorrelationFunctionPsi.C:752
drawCorrelationFunctionPsi.C:753
drawCorrelationFunctionPsi.C:754
drawCorrelationFunctionPsi.C:755
drawCorrelationFunctionPsi.C:756
drawCorrelationFunctionPsi.C:757
drawCorrelationFunctionPsi.C:758
drawCorrelationFunctionPsi.C:759
drawCorrelationFunctionPsi.C:760
drawCorrelationFunctionPsi.C:761
drawCorrelationFunctionPsi.C:762
drawCorrelationFunctionPsi.C:763
drawCorrelationFunctionPsi.C:764
drawCorrelationFunctionPsi.C:765
drawCorrelationFunctionPsi.C:766
drawCorrelationFunctionPsi.C:767
drawCorrelationFunctionPsi.C:768
drawCorrelationFunctionPsi.C:769
drawCorrelationFunctionPsi.C:770
drawCorrelationFunctionPsi.C:771
drawCorrelationFunctionPsi.C:772
drawCorrelationFunctionPsi.C:773
drawCorrelationFunctionPsi.C:774
drawCorrelationFunctionPsi.C:775
drawCorrelationFunctionPsi.C:776
drawCorrelationFunctionPsi.C:777
drawCorrelationFunctionPsi.C:778
drawCorrelationFunctionPsi.C:779
drawCorrelationFunctionPsi.C:780
drawCorrelationFunctionPsi.C:781
drawCorrelationFunctionPsi.C:782
drawCorrelationFunctionPsi.C:783
drawCorrelationFunctionPsi.C:784
drawCorrelationFunctionPsi.C:785
drawCorrelationFunctionPsi.C:786
drawCorrelationFunctionPsi.C:787
drawCorrelationFunctionPsi.C:788
drawCorrelationFunctionPsi.C:789
drawCorrelationFunctionPsi.C:790
drawCorrelationFunctionPsi.C:791
drawCorrelationFunctionPsi.C:792
drawCorrelationFunctionPsi.C:793
drawCorrelationFunctionPsi.C:794
drawCorrelationFunctionPsi.C:795
drawCorrelationFunctionPsi.C:796
drawCorrelationFunctionPsi.C:797
drawCorrelationFunctionPsi.C:798
drawCorrelationFunctionPsi.C:799
drawCorrelationFunctionPsi.C:800
drawCorrelationFunctionPsi.C:801
drawCorrelationFunctionPsi.C:802
drawCorrelationFunctionPsi.C:803
drawCorrelationFunctionPsi.C:804
drawCorrelationFunctionPsi.C:805
drawCorrelationFunctionPsi.C:806
drawCorrelationFunctionPsi.C:807
drawCorrelationFunctionPsi.C:808
drawCorrelationFunctionPsi.C:809
drawCorrelationFunctionPsi.C:810
drawCorrelationFunctionPsi.C:811
drawCorrelationFunctionPsi.C:812
drawCorrelationFunctionPsi.C:813
drawCorrelationFunctionPsi.C:814
drawCorrelationFunctionPsi.C:815
drawCorrelationFunctionPsi.C:816
drawCorrelationFunctionPsi.C:817
drawCorrelationFunctionPsi.C:818
drawCorrelationFunctionPsi.C:819
drawCorrelationFunctionPsi.C:820
drawCorrelationFunctionPsi.C:821
drawCorrelationFunctionPsi.C:822
drawCorrelationFunctionPsi.C:823
drawCorrelationFunctionPsi.C:824
drawCorrelationFunctionPsi.C:825
drawCorrelationFunctionPsi.C:826
drawCorrelationFunctionPsi.C:827
drawCorrelationFunctionPsi.C:828
drawCorrelationFunctionPsi.C:829
drawCorrelationFunctionPsi.C:830
drawCorrelationFunctionPsi.C:831
drawCorrelationFunctionPsi.C:832
drawCorrelationFunctionPsi.C:833
drawCorrelationFunctionPsi.C:834
drawCorrelationFunctionPsi.C:835
drawCorrelationFunctionPsi.C:836
drawCorrelationFunctionPsi.C:837
drawCorrelationFunctionPsi.C:838
drawCorrelationFunctionPsi.C:839
drawCorrelationFunctionPsi.C:840
drawCorrelationFunctionPsi.C:841
drawCorrelationFunctionPsi.C:842
drawCorrelationFunctionPsi.C:843
drawCorrelationFunctionPsi.C:844
drawCorrelationFunctionPsi.C:845
drawCorrelationFunctionPsi.C:846
drawCorrelationFunctionPsi.C:847
drawCorrelationFunctionPsi.C:848
drawCorrelationFunctionPsi.C:849
drawCorrelationFunctionPsi.C:850
drawCorrelationFunctionPsi.C:851
drawCorrelationFunctionPsi.C:852
drawCorrelationFunctionPsi.C:853
drawCorrelationFunctionPsi.C:854
drawCorrelationFunctionPsi.C:855
drawCorrelationFunctionPsi.C:856
drawCorrelationFunctionPsi.C:857
drawCorrelationFunctionPsi.C:858
drawCorrelationFunctionPsi.C:859
drawCorrelationFunctionPsi.C:860
drawCorrelationFunctionPsi.C:861
drawCorrelationFunctionPsi.C:862
drawCorrelationFunctionPsi.C:863
drawCorrelationFunctionPsi.C:864
drawCorrelationFunctionPsi.C:865
drawCorrelationFunctionPsi.C:866
drawCorrelationFunctionPsi.C:867
drawCorrelationFunctionPsi.C:868
drawCorrelationFunctionPsi.C:869
drawCorrelationFunctionPsi.C:870
drawCorrelationFunctionPsi.C:871
drawCorrelationFunctionPsi.C:872
drawCorrelationFunctionPsi.C:873
drawCorrelationFunctionPsi.C:874
drawCorrelationFunctionPsi.C:875
drawCorrelationFunctionPsi.C:876
drawCorrelationFunctionPsi.C:877
drawCorrelationFunctionPsi.C:878
drawCorrelationFunctionPsi.C:879
drawCorrelationFunctionPsi.C:880
drawCorrelationFunctionPsi.C:881
drawCorrelationFunctionPsi.C:882
drawCorrelationFunctionPsi.C:883
drawCorrelationFunctionPsi.C:884
drawCorrelationFunctionPsi.C:885
drawCorrelationFunctionPsi.C:886
drawCorrelationFunctionPsi.C:887
drawCorrelationFunctionPsi.C:888
drawCorrelationFunctionPsi.C:889
drawCorrelationFunctionPsi.C:890
drawCorrelationFunctionPsi.C:891
drawCorrelationFunctionPsi.C:892
drawCorrelationFunctionPsi.C:893
drawCorrelationFunctionPsi.C:894
drawCorrelationFunctionPsi.C:895
drawCorrelationFunctionPsi.C:896
drawCorrelationFunctionPsi.C:897
drawCorrelationFunctionPsi.C:898
drawCorrelationFunctionPsi.C:899
drawCorrelationFunctionPsi.C:900
drawCorrelationFunctionPsi.C:901
drawCorrelationFunctionPsi.C:902
drawCorrelationFunctionPsi.C:903
drawCorrelationFunctionPsi.C:904
drawCorrelationFunctionPsi.C:905
drawCorrelationFunctionPsi.C:906
drawCorrelationFunctionPsi.C:907
drawCorrelationFunctionPsi.C:908
drawCorrelationFunctionPsi.C:909
drawCorrelationFunctionPsi.C:910
drawCorrelationFunctionPsi.C:911
drawCorrelationFunctionPsi.C:912
drawCorrelationFunctionPsi.C:913
drawCorrelationFunctionPsi.C:914
drawCorrelationFunctionPsi.C:915
drawCorrelationFunctionPsi.C:916
drawCorrelationFunctionPsi.C:917
drawCorrelationFunctionPsi.C:918
drawCorrelationFunctionPsi.C:919
drawCorrelationFunctionPsi.C:920
drawCorrelationFunctionPsi.C:921
drawCorrelationFunctionPsi.C:922
drawCorrelationFunctionPsi.C:923
drawCorrelationFunctionPsi.C:924
drawCorrelationFunctionPsi.C:925
drawCorrelationFunctionPsi.C:926
drawCorrelationFunctionPsi.C:927
drawCorrelationFunctionPsi.C:928
drawCorrelationFunctionPsi.C:929
drawCorrelationFunctionPsi.C:930
drawCorrelationFunctionPsi.C:931
drawCorrelationFunctionPsi.C:932
drawCorrelationFunctionPsi.C:933
drawCorrelationFunctionPsi.C:934
drawCorrelationFunctionPsi.C:935
drawCorrelationFunctionPsi.C:936
drawCorrelationFunctionPsi.C:937
drawCorrelationFunctionPsi.C:938
drawCorrelationFunctionPsi.C:939
drawCorrelationFunctionPsi.C:940
drawCorrelationFunctionPsi.C:941
drawCorrelationFunctionPsi.C:942
drawCorrelationFunctionPsi.C:943
drawCorrelationFunctionPsi.C:944
drawCorrelationFunctionPsi.C:945
drawCorrelationFunctionPsi.C:946
drawCorrelationFunctionPsi.C:947
drawCorrelationFunctionPsi.C:948
drawCorrelationFunctionPsi.C:949
drawCorrelationFunctionPsi.C:950
drawCorrelationFunctionPsi.C:951
drawCorrelationFunctionPsi.C:952
drawCorrelationFunctionPsi.C:953
drawCorrelationFunctionPsi.C:954
drawCorrelationFunctionPsi.C:955
drawCorrelationFunctionPsi.C:956
drawCorrelationFunctionPsi.C:957
drawCorrelationFunctionPsi.C:958
drawCorrelationFunctionPsi.C:959
drawCorrelationFunctionPsi.C:960
drawCorrelationFunctionPsi.C:961
drawCorrelationFunctionPsi.C:962
drawCorrelationFunctionPsi.C:963
drawCorrelationFunctionPsi.C:964
drawCorrelationFunctionPsi.C:965
drawCorrelationFunctionPsi.C:966
drawCorrelationFunctionPsi.C:967
drawCorrelationFunctionPsi.C:968
drawCorrelationFunctionPsi.C:969
drawCorrelationFunctionPsi.C:970
drawCorrelationFunctionPsi.C:971
drawCorrelationFunctionPsi.C:972
drawCorrelationFunctionPsi.C:973
drawCorrelationFunctionPsi.C:974
drawCorrelationFunctionPsi.C:975
drawCorrelationFunctionPsi.C:976
drawCorrelationFunctionPsi.C:977
drawCorrelationFunctionPsi.C:978
drawCorrelationFunctionPsi.C:979
drawCorrelationFunctionPsi.C:980
drawCorrelationFunctionPsi.C:981
drawCorrelationFunctionPsi.C:982
drawCorrelationFunctionPsi.C:983
drawCorrelationFunctionPsi.C:984
drawCorrelationFunctionPsi.C:985
drawCorrelationFunctionPsi.C:986
drawCorrelationFunctionPsi.C:987
drawCorrelationFunctionPsi.C:988
drawCorrelationFunctionPsi.C:989
drawCorrelationFunctionPsi.C:990
drawCorrelationFunctionPsi.C:991
drawCorrelationFunctionPsi.C:992
drawCorrelationFunctionPsi.C:993
drawCorrelationFunctionPsi.C:994
drawCorrelationFunctionPsi.C:995
drawCorrelationFunctionPsi.C:996
drawCorrelationFunctionPsi.C:997
drawCorrelationFunctionPsi.C:998
drawCorrelationFunctionPsi.C:999
drawCorrelationFunctionPsi.C:1000
drawCorrelationFunctionPsi.C:1001
drawCorrelationFunctionPsi.C:1002
drawCorrelationFunctionPsi.C:1003
drawCorrelationFunctionPsi.C:1004
drawCorrelationFunctionPsi.C:1005
drawCorrelationFunctionPsi.C:1006
drawCorrelationFunctionPsi.C:1007
drawCorrelationFunctionPsi.C:1008
drawCorrelationFunctionPsi.C:1009
drawCorrelationFunctionPsi.C:1010
drawCorrelationFunctionPsi.C:1011
drawCorrelationFunctionPsi.C:1012
drawCorrelationFunctionPsi.C:1013
drawCorrelationFunctionPsi.C:1014
drawCorrelationFunctionPsi.C:1015
drawCorrelationFunctionPsi.C:1016
drawCorrelationFunctionPsi.C:1017
drawCorrelationFunctionPsi.C:1018
drawCorrelationFunctionPsi.C:1019
drawCorrelationFunctionPsi.C:1020
drawCorrelationFunctionPsi.C:1021
drawCorrelationFunctionPsi.C:1022
drawCorrelationFunctionPsi.C:1023
drawCorrelationFunctionPsi.C:1024
drawCorrelationFunctionPsi.C:1025
drawCorrelationFunctionPsi.C:1026
drawCorrelationFunctionPsi.C:1027
drawCorrelationFunctionPsi.C:1028
drawCorrelationFunctionPsi.C:1029
drawCorrelationFunctionPsi.C:1030
drawCorrelationFunctionPsi.C:1031
drawCorrelationFunctionPsi.C:1032
drawCorrelationFunctionPsi.C:1033
drawCorrelationFunctionPsi.C:1034
drawCorrelationFunctionPsi.C:1035
drawCorrelationFunctionPsi.C:1036
drawCorrelationFunctionPsi.C:1037
drawCorrelationFunctionPsi.C:1038
drawCorrelationFunctionPsi.C:1039
drawCorrelationFunctionPsi.C:1040
drawCorrelationFunctionPsi.C:1041
drawCorrelationFunctionPsi.C:1042
drawCorrelationFunctionPsi.C:1043
drawCorrelationFunctionPsi.C:1044
drawCorrelationFunctionPsi.C:1045
drawCorrelationFunctionPsi.C:1046
drawCorrelationFunctionPsi.C:1047
drawCorrelationFunctionPsi.C:1048
drawCorrelationFunctionPsi.C:1049
drawCorrelationFunctionPsi.C:1050
drawCorrelationFunctionPsi.C:1051
drawCorrelationFunctionPsi.C:1052
drawCorrelationFunctionPsi.C:1053
drawCorrelationFunctionPsi.C:1054
drawCorrelationFunctionPsi.C:1055
drawCorrelationFunctionPsi.C:1056
drawCorrelationFunctionPsi.C:1057
drawCorrelationFunctionPsi.C:1058
drawCorrelationFunctionPsi.C:1059
drawCorrelationFunctionPsi.C:1060
drawCorrelationFunctionPsi.C:1061
drawCorrelationFunctionPsi.C:1062
drawCorrelationFunctionPsi.C:1063
drawCorrelationFunctionPsi.C:1064
drawCorrelationFunctionPsi.C:1065
drawCorrelationFunctionPsi.C:1066
drawCorrelationFunctionPsi.C:1067
drawCorrelationFunctionPsi.C:1068
drawCorrelationFunctionPsi.C:1069
drawCorrelationFunctionPsi.C:1070
drawCorrelationFunctionPsi.C:1071
drawCorrelationFunctionPsi.C:1072
drawCorrelationFunctionPsi.C:1073
drawCorrelationFunctionPsi.C:1074
drawCorrelationFunctionPsi.C:1075
drawCorrelationFunctionPsi.C:1076
drawCorrelationFunctionPsi.C:1077
drawCorrelationFunctionPsi.C:1078
drawCorrelationFunctionPsi.C:1079
drawCorrelationFunctionPsi.C:1080
drawCorrelationFunctionPsi.C:1081
drawCorrelationFunctionPsi.C:1082
drawCorrelationFunctionPsi.C:1083
drawCorrelationFunctionPsi.C:1084
drawCorrelationFunctionPsi.C:1085
drawCorrelationFunctionPsi.C:1086
drawCorrelationFunctionPsi.C:1087
drawCorrelationFunctionPsi.C:1088
drawCorrelationFunctionPsi.C:1089
drawCorrelationFunctionPsi.C:1090
drawCorrelationFunctionPsi.C:1091
drawCorrelationFunctionPsi.C:1092
drawCorrelationFunctionPsi.C:1093
drawCorrelationFunctionPsi.C:1094
drawCorrelationFunctionPsi.C:1095
drawCorrelationFunctionPsi.C:1096
drawCorrelationFunctionPsi.C:1097
drawCorrelationFunctionPsi.C:1098
drawCorrelationFunctionPsi.C:1099
drawCorrelationFunctionPsi.C:1100
drawCorrelationFunctionPsi.C:1101
drawCorrelationFunctionPsi.C:1102
drawCorrelationFunctionPsi.C:1103
drawCorrelationFunctionPsi.C:1104
drawCorrelationFunctionPsi.C:1105
drawCorrelationFunctionPsi.C:1106
drawCorrelationFunctionPsi.C:1107
drawCorrelationFunctionPsi.C:1108
drawCorrelationFunctionPsi.C:1109
drawCorrelationFunctionPsi.C:1110
drawCorrelationFunctionPsi.C:1111
drawCorrelationFunctionPsi.C:1112
drawCorrelationFunctionPsi.C:1113
drawCorrelationFunctionPsi.C:1114
drawCorrelationFunctionPsi.C:1115
drawCorrelationFunctionPsi.C:1116
drawCorrelationFunctionPsi.C:1117
drawCorrelationFunctionPsi.C:1118
drawCorrelationFunctionPsi.C:1119
drawCorrelationFunctionPsi.C:1120
drawCorrelationFunctionPsi.C:1121
drawCorrelationFunctionPsi.C:1122
drawCorrelationFunctionPsi.C:1123
drawCorrelationFunctionPsi.C:1124
drawCorrelationFunctionPsi.C:1125
drawCorrelationFunctionPsi.C:1126
drawCorrelationFunctionPsi.C:1127
drawCorrelationFunctionPsi.C:1128
drawCorrelationFunctionPsi.C:1129
drawCorrelationFunctionPsi.C:1130
drawCorrelationFunctionPsi.C:1131
drawCorrelationFunctionPsi.C:1132
drawCorrelationFunctionPsi.C:1133
drawCorrelationFunctionPsi.C:1134
drawCorrelationFunctionPsi.C:1135
drawCorrelationFunctionPsi.C:1136
drawCorrelationFunctionPsi.C:1137
drawCorrelationFunctionPsi.C:1138
drawCorrelationFunctionPsi.C:1139
drawCorrelationFunctionPsi.C:1140
drawCorrelationFunctionPsi.C:1141
drawCorrelationFunctionPsi.C:1142
drawCorrelationFunctionPsi.C:1143
drawCorrelationFunctionPsi.C:1144
drawCorrelationFunctionPsi.C:1145
drawCorrelationFunctionPsi.C:1146
drawCorrelationFunctionPsi.C:1147
drawCorrelationFunctionPsi.C:1148
drawCorrelationFunctionPsi.C:1149
drawCorrelationFunctionPsi.C:1150
drawCorrelationFunctionPsi.C:1151
drawCorrelationFunctionPsi.C:1152
drawCorrelationFunctionPsi.C:1153
drawCorrelationFunctionPsi.C:1154
drawCorrelationFunctionPsi.C:1155
drawCorrelationFunctionPsi.C:1156
drawCorrelationFunctionPsi.C:1157
drawCorrelationFunctionPsi.C:1158
drawCorrelationFunctionPsi.C:1159
drawCorrelationFunctionPsi.C:1160
drawCorrelationFunctionPsi.C:1161
drawCorrelationFunctionPsi.C:1162
drawCorrelationFunctionPsi.C:1163
drawCorrelationFunctionPsi.C:1164
drawCorrelationFunctionPsi.C:1165
drawCorrelationFunctionPsi.C:1166
drawCorrelationFunctionPsi.C:1167
drawCorrelationFunctionPsi.C:1168
drawCorrelationFunctionPsi.C:1169
drawCorrelationFunctionPsi.C:1170
drawCorrelationFunctionPsi.C:1171
drawCorrelationFunctionPsi.C:1172
drawCorrelationFunctionPsi.C:1173
drawCorrelationFunctionPsi.C:1174
drawCorrelationFunctionPsi.C:1175
drawCorrelationFunctionPsi.C:1176
drawCorrelationFunctionPsi.C:1177
drawCorrelationFunctionPsi.C:1178
drawCorrelationFunctionPsi.C:1179
drawCorrelationFunctionPsi.C:1180
drawCorrelationFunctionPsi.C:1181
drawCorrelationFunctionPsi.C:1182
drawCorrelationFunctionPsi.C:1183
drawCorrelationFunctionPsi.C:1184
drawCorrelationFunctionPsi.C:1185
drawCorrelationFunctionPsi.C:1186
drawCorrelationFunctionPsi.C:1187
drawCorrelationFunctionPsi.C:1188
drawCorrelationFunctionPsi.C:1189
drawCorrelationFunctionPsi.C:1190
drawCorrelationFunctionPsi.C:1191
drawCorrelationFunctionPsi.C:1192
drawCorrelationFunctionPsi.C:1193
drawCorrelationFunctionPsi.C:1194
drawCorrelationFunctionPsi.C:1195
drawCorrelationFunctionPsi.C:1196
drawCorrelationFunctionPsi.C:1197
drawCorrelationFunctionPsi.C:1198
drawCorrelationFunctionPsi.C:1199
drawCorrelationFunctionPsi.C:1200
drawCorrelationFunctionPsi.C:1201
drawCorrelationFunctionPsi.C:1202
drawCorrelationFunctionPsi.C:1203
drawCorrelationFunctionPsi.C:1204
drawCorrelationFunctionPsi.C:1205
drawCorrelationFunctionPsi.C:1206
drawCorrelationFunctionPsi.C:1207
drawCorrelationFunctionPsi.C:1208
drawCorrelationFunctionPsi.C:1209
drawCorrelationFunctionPsi.C:1210
drawCorrelationFunctionPsi.C:1211
drawCorrelationFunctionPsi.C:1212
drawCorrelationFunctionPsi.C:1213
drawCorrelationFunctionPsi.C:1214
drawCorrelationFunctionPsi.C:1215
drawCorrelationFunctionPsi.C:1216
drawCorrelationFunctionPsi.C:1217
drawCorrelationFunctionPsi.C:1218
drawCorrelationFunctionPsi.C:1219
drawCorrelationFunctionPsi.C:1220
drawCorrelationFunctionPsi.C:1221
drawCorrelationFunctionPsi.C:1222
drawCorrelationFunctionPsi.C:1223
drawCorrelationFunctionPsi.C:1224
drawCorrelationFunctionPsi.C:1225
drawCorrelationFunctionPsi.C:1226
drawCorrelationFunctionPsi.C:1227
drawCorrelationFunctionPsi.C:1228
drawCorrelationFunctionPsi.C:1229
drawCorrelationFunctionPsi.C:1230
drawCorrelationFunctionPsi.C:1231
drawCorrelationFunctionPsi.C:1232
drawCorrelationFunctionPsi.C:1233
drawCorrelationFunctionPsi.C:1234
drawCorrelationFunctionPsi.C:1235
drawCorrelationFunctionPsi.C:1236
drawCorrelationFunctionPsi.C:1237
drawCorrelationFunctionPsi.C:1238
drawCorrelationFunctionPsi.C:1239
drawCorrelationFunctionPsi.C:1240
drawCorrelationFunctionPsi.C:1241
drawCorrelationFunctionPsi.C:1242
drawCorrelationFunctionPsi.C:1243
drawCorrelationFunctionPsi.C:1244
drawCorrelationFunctionPsi.C:1245
drawCorrelationFunctionPsi.C:1246
drawCorrelationFunctionPsi.C:1247
drawCorrelationFunctionPsi.C:1248
drawCorrelationFunctionPsi.C:1249
drawCorrelationFunctionPsi.C:1250
drawCorrelationFunctionPsi.C:1251
drawCorrelationFunctionPsi.C:1252
drawCorrelationFunctionPsi.C:1253
drawCorrelationFunctionPsi.C:1254
drawCorrelationFunctionPsi.C:1255
drawCorrelationFunctionPsi.C:1256
drawCorrelationFunctionPsi.C:1257
drawCorrelationFunctionPsi.C:1258
drawCorrelationFunctionPsi.C:1259
drawCorrelationFunctionPsi.C:1260
drawCorrelationFunctionPsi.C:1261
drawCorrelationFunctionPsi.C:1262
drawCorrelationFunctionPsi.C:1263
drawCorrelationFunctionPsi.C:1264
drawCorrelationFunctionPsi.C:1265
drawCorrelationFunctionPsi.C:1266
drawCorrelationFunctionPsi.C:1267
drawCorrelationFunctionPsi.C:1268
drawCorrelationFunctionPsi.C:1269
drawCorrelationFunctionPsi.C:1270
drawCorrelationFunctionPsi.C:1271
drawCorrelationFunctionPsi.C:1272
drawCorrelationFunctionPsi.C:1273
drawCorrelationFunctionPsi.C:1274
drawCorrelationFunctionPsi.C:1275
drawCorrelationFunctionPsi.C:1276
drawCorrelationFunctionPsi.C:1277
drawCorrelationFunctionPsi.C:1278
drawCorrelationFunctionPsi.C:1279
drawCorrelationFunctionPsi.C:1280
drawCorrelationFunctionPsi.C:1281
drawCorrelationFunctionPsi.C:1282
drawCorrelationFunctionPsi.C:1283
drawCorrelationFunctionPsi.C:1284
drawCorrelationFunctionPsi.C:1285
drawCorrelationFunctionPsi.C:1286
drawCorrelationFunctionPsi.C:1287
drawCorrelationFunctionPsi.C:1288
drawCorrelationFunctionPsi.C:1289
drawCorrelationFunctionPsi.C:1290
drawCorrelationFunctionPsi.C:1291
drawCorrelationFunctionPsi.C:1292
drawCorrelationFunctionPsi.C:1293
drawCorrelationFunctionPsi.C:1294
drawCorrelationFunctionPsi.C:1295
drawCorrelationFunctionPsi.C:1296
drawCorrelationFunctionPsi.C:1297
drawCorrelationFunctionPsi.C:1298
drawCorrelationFunctionPsi.C:1299
drawCorrelationFunctionPsi.C:1300
drawCorrelationFunctionPsi.C:1301
drawCorrelationFunctionPsi.C:1302
drawCorrelationFunctionPsi.C:1303
drawCorrelationFunctionPsi.C:1304
drawCorrelationFunctionPsi.C:1305
drawCorrelationFunctionPsi.C:1306
drawCorrelationFunctionPsi.C:1307
drawCorrelationFunctionPsi.C:1308
drawCorrelationFunctionPsi.C:1309
drawCorrelationFunctionPsi.C:1310
drawCorrelationFunctionPsi.C:1311
drawCorrelationFunctionPsi.C:1312
drawCorrelationFunctionPsi.C:1313
drawCorrelationFunctionPsi.C:1314
drawCorrelationFunctionPsi.C:1315
drawCorrelationFunctionPsi.C:1316
drawCorrelationFunctionPsi.C:1317
drawCorrelationFunctionPsi.C:1318
drawCorrelationFunctionPsi.C:1319
drawCorrelationFunctionPsi.C:1320
drawCorrelationFunctionPsi.C:1321
drawCorrelationFunctionPsi.C:1322
drawCorrelationFunctionPsi.C:1323
drawCorrelationFunctionPsi.C:1324
drawCorrelationFunctionPsi.C:1325
drawCorrelationFunctionPsi.C:1326
drawCorrelationFunctionPsi.C:1327
drawCorrelationFunctionPsi.C:1328
drawCorrelationFunctionPsi.C:1329
drawCorrelationFunctionPsi.C:1330
drawCorrelationFunctionPsi.C:1331
drawCorrelationFunctionPsi.C:1332
drawCorrelationFunctionPsi.C:1333
drawCorrelationFunctionPsi.C:1334
drawCorrelationFunctionPsi.C:1335
drawCorrelationFunctionPsi.C:1336
drawCorrelationFunctionPsi.C:1337
drawCorrelationFunctionPsi.C:1338
drawCorrelationFunctionPsi.C:1339
drawCorrelationFunctionPsi.C:1340
drawCorrelationFunctionPsi.C:1341
drawCorrelationFunctionPsi.C:1342
drawCorrelationFunctionPsi.C:1343
drawCorrelationFunctionPsi.C:1344
drawCorrelationFunctionPsi.C:1345
drawCorrelationFunctionPsi.C:1346
drawCorrelationFunctionPsi.C:1347
drawCorrelationFunctionPsi.C:1348
drawCorrelationFunctionPsi.C:1349
drawCorrelationFunctionPsi.C:1350
drawCorrelationFunctionPsi.C:1351
drawCorrelationFunctionPsi.C:1352
drawCorrelationFunctionPsi.C:1353
drawCorrelationFunctionPsi.C:1354
drawCorrelationFunctionPsi.C:1355
drawCorrelationFunctionPsi.C:1356
drawCorrelationFunctionPsi.C:1357
drawCorrelationFunctionPsi.C:1358
drawCorrelationFunctionPsi.C:1359
drawCorrelationFunctionPsi.C:1360
drawCorrelationFunctionPsi.C:1361
drawCorrelationFunctionPsi.C:1362
drawCorrelationFunctionPsi.C:1363
drawCorrelationFunctionPsi.C:1364
drawCorrelationFunctionPsi.C:1365
drawCorrelationFunctionPsi.C:1366
drawCorrelationFunctionPsi.C:1367
drawCorrelationFunctionPsi.C:1368
drawCorrelationFunctionPsi.C:1369
drawCorrelationFunctionPsi.C:1370
drawCorrelationFunctionPsi.C:1371
drawCorrelationFunctionPsi.C:1372
drawCorrelationFunctionPsi.C:1373
drawCorrelationFunctionPsi.C:1374
drawCorrelationFunctionPsi.C:1375
drawCorrelationFunctionPsi.C:1376
drawCorrelationFunctionPsi.C:1377
drawCorrelationFunctionPsi.C:1378
drawCorrelationFunctionPsi.C:1379
drawCorrelationFunctionPsi.C:1380
drawCorrelationFunctionPsi.C:1381
drawCorrelationFunctionPsi.C:1382
drawCorrelationFunctionPsi.C:1383
drawCorrelationFunctionPsi.C:1384
drawCorrelationFunctionPsi.C:1385
drawCorrelationFunctionPsi.C:1386
drawCorrelationFunctionPsi.C:1387
drawCorrelationFunctionPsi.C:1388
drawCorrelationFunctionPsi.C:1389
drawCorrelationFunctionPsi.C:1390
drawCorrelationFunctionPsi.C:1391
drawCorrelationFunctionPsi.C:1392
drawCorrelationFunctionPsi.C:1393
drawCorrelationFunctionPsi.C:1394
drawCorrelationFunctionPsi.C:1395
drawCorrelationFunctionPsi.C:1396
drawCorrelationFunctionPsi.C:1397
drawCorrelationFunctionPsi.C:1398
drawCorrelationFunctionPsi.C:1399
drawCorrelationFunctionPsi.C:1400
drawCorrelationFunctionPsi.C:1401
drawCorrelationFunctionPsi.C:1402
drawCorrelationFunctionPsi.C:1403
drawCorrelationFunctionPsi.C:1404
drawCorrelationFunctionPsi.C:1405
drawCorrelationFunctionPsi.C:1406
drawCorrelationFunctionPsi.C:1407
drawCorrelationFunctionPsi.C:1408
drawCorrelationFunctionPsi.C:1409
drawCorrelationFunctionPsi.C:1410
drawCorrelationFunctionPsi.C:1411
drawCorrelationFunctionPsi.C:1412
drawCorrelationFunctionPsi.C:1413
drawCorrelationFunctionPsi.C:1414
drawCorrelationFunctionPsi.C:1415
drawCorrelationFunctionPsi.C:1416
drawCorrelationFunctionPsi.C:1417
drawCorrelationFunctionPsi.C:1418
drawCorrelationFunctionPsi.C:1419
drawCorrelationFunctionPsi.C:1420
drawCorrelationFunctionPsi.C:1421
drawCorrelationFunctionPsi.C:1422
drawCorrelationFunctionPsi.C:1423
drawCorrelationFunctionPsi.C:1424
drawCorrelationFunctionPsi.C:1425
drawCorrelationFunctionPsi.C:1426
drawCorrelationFunctionPsi.C:1427
drawCorrelationFunctionPsi.C:1428
drawCorrelationFunctionPsi.C:1429
drawCorrelationFunctionPsi.C:1430
drawCorrelationFunctionPsi.C:1431
drawCorrelationFunctionPsi.C:1432
drawCorrelationFunctionPsi.C:1433
drawCorrelationFunctionPsi.C:1434
drawCorrelationFunctionPsi.C:1435
drawCorrelationFunctionPsi.C:1436
drawCorrelationFunctionPsi.C:1437
drawCorrelationFunctionPsi.C:1438
drawCorrelationFunctionPsi.C:1439
drawCorrelationFunctionPsi.C:1440
drawCorrelationFunctionPsi.C:1441
drawCorrelationFunctionPsi.C:1442
drawCorrelationFunctionPsi.C:1443
drawCorrelationFunctionPsi.C:1444
drawCorrelationFunctionPsi.C:1445
drawCorrelationFunctionPsi.C:1446
drawCorrelationFunctionPsi.C:1447
drawCorrelationFunctionPsi.C:1448
drawCorrelationFunctionPsi.C:1449
drawCorrelationFunctionPsi.C:1450
drawCorrelationFunctionPsi.C:1451
drawCorrelationFunctionPsi.C:1452
drawCorrelationFunctionPsi.C:1453
drawCorrelationFunctionPsi.C:1454
drawCorrelationFunctionPsi.C:1455
drawCorrelationFunctionPsi.C:1456
drawCorrelationFunctionPsi.C:1457
drawCorrelationFunctionPsi.C:1458
drawCorrelationFunctionPsi.C:1459
drawCorrelationFunctionPsi.C:1460
drawCorrelationFunctionPsi.C:1461
drawCorrelationFunctionPsi.C:1462
drawCorrelationFunctionPsi.C:1463
drawCorrelationFunctionPsi.C:1464
drawCorrelationFunctionPsi.C:1465
drawCorrelationFunctionPsi.C:1466
drawCorrelationFunctionPsi.C:1467
drawCorrelationFunctionPsi.C:1468
drawCorrelationFunctionPsi.C:1469
drawCorrelationFunctionPsi.C:1470
drawCorrelationFunctionPsi.C:1471
drawCorrelationFunctionPsi.C:1472
drawCorrelationFunctionPsi.C:1473
drawCorrelationFunctionPsi.C:1474
drawCorrelationFunctionPsi.C:1475
drawCorrelationFunctionPsi.C:1476
drawCorrelationFunctionPsi.C:1477
drawCorrelationFunctionPsi.C:1478
drawCorrelationFunctionPsi.C:1479
drawCorrelationFunctionPsi.C:1480
drawCorrelationFunctionPsi.C:1481
drawCorrelationFunctionPsi.C:1482
drawCorrelationFunctionPsi.C:1483
drawCorrelationFunctionPsi.C:1484
drawCorrelationFunctionPsi.C:1485
drawCorrelationFunctionPsi.C:1486
drawCorrelationFunctionPsi.C:1487
drawCorrelationFunctionPsi.C:1488
drawCorrelationFunctionPsi.C:1489
drawCorrelationFunctionPsi.C:1490
drawCorrelationFunctionPsi.C:1491
drawCorrelationFunctionPsi.C:1492
drawCorrelationFunctionPsi.C:1493
drawCorrelationFunctionPsi.C:1494
drawCorrelationFunctionPsi.C:1495
drawCorrelationFunctionPsi.C:1496
drawCorrelationFunctionPsi.C:1497
drawCorrelationFunctionPsi.C:1498
drawCorrelationFunctionPsi.C:1499
drawCorrelationFunctionPsi.C:1500
drawCorrelationFunctionPsi.C:1501
drawCorrelationFunctionPsi.C:1502
drawCorrelationFunctionPsi.C:1503
drawCorrelationFunctionPsi.C:1504
drawCorrelationFunctionPsi.C:1505
drawCorrelationFunctionPsi.C:1506
drawCorrelationFunctionPsi.C:1507
drawCorrelationFunctionPsi.C:1508
drawCorrelationFunctionPsi.C:1509
drawCorrelationFunctionPsi.C:1510
drawCorrelationFunctionPsi.C:1511
drawCorrelationFunctionPsi.C:1512
drawCorrelationFunctionPsi.C:1513
drawCorrelationFunctionPsi.C:1514
drawCorrelationFunctionPsi.C:1515
drawCorrelationFunctionPsi.C:1516
drawCorrelationFunctionPsi.C:1517
drawCorrelationFunctionPsi.C:1518
drawCorrelationFunctionPsi.C:1519
drawCorrelationFunctionPsi.C:1520
drawCorrelationFunctionPsi.C:1521
drawCorrelationFunctionPsi.C:1522
drawCorrelationFunctionPsi.C:1523
drawCorrelationFunctionPsi.C:1524
drawCorrelationFunctionPsi.C:1525
drawCorrelationFunctionPsi.C:1526
drawCorrelationFunctionPsi.C:1527
drawCorrelationFunctionPsi.C:1528
drawCorrelationFunctionPsi.C:1529
drawCorrelationFunctionPsi.C:1530
drawCorrelationFunctionPsi.C:1531
drawCorrelationFunctionPsi.C:1532
drawCorrelationFunctionPsi.C:1533
drawCorrelationFunctionPsi.C:1534
drawCorrelationFunctionPsi.C:1535
drawCorrelationFunctionPsi.C:1536
drawCorrelationFunctionPsi.C:1537
drawCorrelationFunctionPsi.C:1538
drawCorrelationFunctionPsi.C:1539
drawCorrelationFunctionPsi.C:1540
drawCorrelationFunctionPsi.C:1541
drawCorrelationFunctionPsi.C:1542
drawCorrelationFunctionPsi.C:1543
drawCorrelationFunctionPsi.C:1544
drawCorrelationFunctionPsi.C:1545
drawCorrelationFunctionPsi.C:1546
drawCorrelationFunctionPsi.C:1547
drawCorrelationFunctionPsi.C:1548
drawCorrelationFunctionPsi.C:1549
drawCorrelationFunctionPsi.C:1550
drawCorrelationFunctionPsi.C:1551
drawCorrelationFunctionPsi.C:1552
drawCorrelationFunctionPsi.C:1553
drawCorrelationFunctionPsi.C:1554
drawCorrelationFunctionPsi.C:1555
drawCorrelationFunctionPsi.C:1556
drawCorrelationFunctionPsi.C:1557
drawCorrelationFunctionPsi.C:1558
drawCorrelationFunctionPsi.C:1559
drawCorrelationFunctionPsi.C:1560
drawCorrelationFunctionPsi.C:1561
drawCorrelationFunctionPsi.C:1562
drawCorrelationFunctionPsi.C:1563
drawCorrelationFunctionPsi.C:1564
drawCorrelationFunctionPsi.C:1565
drawCorrelationFunctionPsi.C:1566
drawCorrelationFunctionPsi.C:1567
drawCorrelationFunctionPsi.C:1568
drawCorrelationFunctionPsi.C:1569
drawCorrelationFunctionPsi.C:1570
drawCorrelationFunctionPsi.C:1571
drawCorrelationFunctionPsi.C:1572
drawCorrelationFunctionPsi.C:1573
drawCorrelationFunctionPsi.C:1574
drawCorrelationFunctionPsi.C:1575
drawCorrelationFunctionPsi.C:1576
drawCorrelationFunctionPsi.C:1577
drawCorrelationFunctionPsi.C:1578
drawCorrelationFunctionPsi.C:1579
drawCorrelationFunctionPsi.C:1580
drawCorrelationFunctionPsi.C:1581
drawCorrelationFunctionPsi.C:1582
drawCorrelationFunctionPsi.C:1583
drawCorrelationFunctionPsi.C:1584
drawCorrelationFunctionPsi.C:1585
drawCorrelationFunctionPsi.C:1586
drawCorrelationFunctionPsi.C:1587
drawCorrelationFunctionPsi.C:1588
drawCorrelationFunctionPsi.C:1589
drawCorrelationFunctionPsi.C:1590
drawCorrelationFunctionPsi.C:1591
drawCorrelationFunctionPsi.C:1592
drawCorrelationFunctionPsi.C:1593
drawCorrelationFunctionPsi.C:1594
drawCorrelationFunctionPsi.C:1595
drawCorrelationFunctionPsi.C:1596
drawCorrelationFunctionPsi.C:1597
drawCorrelationFunctionPsi.C:1598
drawCorrelationFunctionPsi.C:1599
drawCorrelationFunctionPsi.C:1600
drawCorrelationFunctionPsi.C:1601
drawCorrelationFunctionPsi.C:1602
drawCorrelationFunctionPsi.C:1603
drawCorrelationFunctionPsi.C:1604
drawCorrelationFunctionPsi.C:1605
drawCorrelationFunctionPsi.C:1606
drawCorrelationFunctionPsi.C:1607
drawCorrelationFunctionPsi.C:1608
drawCorrelationFunctionPsi.C:1609
drawCorrelationFunctionPsi.C:1610
drawCorrelationFunctionPsi.C:1611
drawCorrelationFunctionPsi.C:1612
drawCorrelationFunctionPsi.C:1613
drawCorrelationFunctionPsi.C:1614
drawCorrelationFunctionPsi.C:1615
drawCorrelationFunctionPsi.C:1616
drawCorrelationFunctionPsi.C:1617
drawCorrelationFunctionPsi.C:1618
drawCorrelationFunctionPsi.C:1619
drawCorrelationFunctionPsi.C:1620
drawCorrelationFunctionPsi.C:1621
drawCorrelationFunctionPsi.C:1622
drawCorrelationFunctionPsi.C:1623
drawCorrelationFunctionPsi.C:1624
drawCorrelationFunctionPsi.C:1625
drawCorrelationFunctionPsi.C:1626
drawCorrelationFunctionPsi.C:1627
drawCorrelationFunctionPsi.C:1628
drawCorrelationFunctionPsi.C:1629
drawCorrelationFunctionPsi.C:1630
drawCorrelationFunctionPsi.C:1631
drawCorrelationFunctionPsi.C:1632
drawCorrelationFunctionPsi.C:1633
drawCorrelationFunctionPsi.C:1634
drawCorrelationFunctionPsi.C:1635
drawCorrelationFunctionPsi.C:1636
drawCorrelationFunctionPsi.C:1637
drawCorrelationFunctionPsi.C:1638
drawCorrelationFunctionPsi.C:1639
drawCorrelationFunctionPsi.C:1640
drawCorrelationFunctionPsi.C:1641
drawCorrelationFunctionPsi.C:1642
drawCorrelationFunctionPsi.C:1643
drawCorrelationFunctionPsi.C:1644
drawCorrelationFunctionPsi.C:1645
drawCorrelationFunctionPsi.C:1646
drawCorrelationFunctionPsi.C:1647
drawCorrelationFunctionPsi.C:1648
drawCorrelationFunctionPsi.C:1649
drawCorrelationFunctionPsi.C:1650
drawCorrelationFunctionPsi.C:1651
drawCorrelationFunctionPsi.C:1652
drawCorrelationFunctionPsi.C:1653
drawCorrelationFunctionPsi.C:1654
drawCorrelationFunctionPsi.C:1655
drawCorrelationFunctionPsi.C:1656
drawCorrelationFunctionPsi.C:1657
drawCorrelationFunctionPsi.C:1658
drawCorrelationFunctionPsi.C:1659
drawCorrelationFunctionPsi.C:1660
drawCorrelationFunctionPsi.C:1661
drawCorrelationFunctionPsi.C:1662
drawCorrelationFunctionPsi.C:1663
drawCorrelationFunctionPsi.C:1664
drawCorrelationFunctionPsi.C:1665
drawCorrelationFunctionPsi.C:1666
drawCorrelationFunctionPsi.C:1667
drawCorrelationFunctionPsi.C:1668
drawCorrelationFunctionPsi.C:1669
drawCorrelationFunctionPsi.C:1670
drawCorrelationFunctionPsi.C:1671
drawCorrelationFunctionPsi.C:1672
drawCorrelationFunctionPsi.C:1673
drawCorrelationFunctionPsi.C:1674
drawCorrelationFunctionPsi.C:1675
drawCorrelationFunctionPsi.C:1676
drawCorrelationFunctionPsi.C:1677
drawCorrelationFunctionPsi.C:1678
drawCorrelationFunctionPsi.C:1679
drawCorrelationFunctionPsi.C:1680
drawCorrelationFunctionPsi.C:1681
drawCorrelationFunctionPsi.C:1682
drawCorrelationFunctionPsi.C:1683
drawCorrelationFunctionPsi.C:1684
drawCorrelationFunctionPsi.C:1685
drawCorrelationFunctionPsi.C:1686
drawCorrelationFunctionPsi.C:1687
drawCorrelationFunctionPsi.C:1688
drawCorrelationFunctionPsi.C:1689
drawCorrelationFunctionPsi.C:1690
drawCorrelationFunctionPsi.C:1691
drawCorrelationFunctionPsi.C:1692
drawCorrelationFunctionPsi.C:1693
drawCorrelationFunctionPsi.C:1694
drawCorrelationFunctionPsi.C:1695
drawCorrelationFunctionPsi.C:1696
drawCorrelationFunctionPsi.C:1697
drawCorrelationFunctionPsi.C:1698
drawCorrelationFunctionPsi.C:1699
drawCorrelationFunctionPsi.C:1700
drawCorrelationFunctionPsi.C:1701
drawCorrelationFunctionPsi.C:1702
drawCorrelationFunctionPsi.C:1703
drawCorrelationFunctionPsi.C:1704
drawCorrelationFunctionPsi.C:1705
drawCorrelationFunctionPsi.C:1706
drawCorrelationFunctionPsi.C:1707
drawCorrelationFunctionPsi.C:1708
drawCorrelationFunctionPsi.C:1709
drawCorrelationFunctionPsi.C:1710
drawCorrelationFunctionPsi.C:1711
drawCorrelationFunctionPsi.C:1712
drawCorrelationFunctionPsi.C:1713
drawCorrelationFunctionPsi.C:1714
drawCorrelationFunctionPsi.C:1715
drawCorrelationFunctionPsi.C:1716
drawCorrelationFunctionPsi.C:1717
drawCorrelationFunctionPsi.C:1718
drawCorrelationFunctionPsi.C:1719
drawCorrelationFunctionPsi.C:1720
drawCorrelationFunctionPsi.C:1721
drawCorrelationFunctionPsi.C:1722
drawCorrelationFunctionPsi.C:1723
drawCorrelationFunctionPsi.C:1724
drawCorrelationFunctionPsi.C:1725
drawCorrelationFunctionPsi.C:1726
drawCorrelationFunctionPsi.C:1727
drawCorrelationFunctionPsi.C:1728
drawCorrelationFunctionPsi.C:1729
drawCorrelationFunctionPsi.C:1730
drawCorrelationFunctionPsi.C:1731
drawCorrelationFunctionPsi.C:1732
drawCorrelationFunctionPsi.C:1733
drawCorrelationFunctionPsi.C:1734
drawCorrelationFunctionPsi.C:1735
drawCorrelationFunctionPsi.C:1736
drawCorrelationFunctionPsi.C:1737
drawCorrelationFunctionPsi.C:1738
drawCorrelationFunctionPsi.C:1739
drawCorrelationFunctionPsi.C:1740
drawCorrelationFunctionPsi.C:1741
drawCorrelationFunctionPsi.C:1742
drawCorrelationFunctionPsi.C:1743
drawCorrelationFunctionPsi.C:1744
drawCorrelationFunctionPsi.C:1745
drawCorrelationFunctionPsi.C:1746
drawCorrelationFunctionPsi.C:1747
drawCorrelationFunctionPsi.C:1748
drawCorrelationFunctionPsi.C:1749
drawCorrelationFunctionPsi.C:1750
drawCorrelationFunctionPsi.C:1751
drawCorrelationFunctionPsi.C:1752
drawCorrelationFunctionPsi.C:1753
drawCorrelationFunctionPsi.C:1754
drawCorrelationFunctionPsi.C:1755
drawCorrelationFunctionPsi.C:1756
drawCorrelationFunctionPsi.C:1757
drawCorrelationFunctionPsi.C:1758
drawCorrelationFunctionPsi.C:1759
drawCorrelationFunctionPsi.C:1760
drawCorrelationFunctionPsi.C:1761
drawCorrelationFunctionPsi.C:1762
drawCorrelationFunctionPsi.C:1763
drawCorrelationFunctionPsi.C:1764
drawCorrelationFunctionPsi.C:1765
drawCorrelationFunctionPsi.C:1766
drawCorrelationFunctionPsi.C:1767
drawCorrelationFunctionPsi.C:1768
drawCorrelationFunctionPsi.C:1769
drawCorrelationFunctionPsi.C:1770
drawCorrelationFunctionPsi.C:1771
drawCorrelationFunctionPsi.C:1772
drawCorrelationFunctionPsi.C:1773
drawCorrelationFunctionPsi.C:1774
drawCorrelationFunctionPsi.C:1775
drawCorrelationFunctionPsi.C:1776
drawCorrelationFunctionPsi.C:1777
drawCorrelationFunctionPsi.C:1778
drawCorrelationFunctionPsi.C:1779
drawCorrelationFunctionPsi.C:1780
drawCorrelationFunctionPsi.C:1781
drawCorrelationFunctionPsi.C:1782
drawCorrelationFunctionPsi.C:1783
drawCorrelationFunctionPsi.C:1784
drawCorrelationFunctionPsi.C:1785
drawCorrelationFunctionPsi.C:1786
drawCorrelationFunctionPsi.C:1787
drawCorrelationFunctionPsi.C:1788
drawCorrelationFunctionPsi.C:1789
drawCorrelationFunctionPsi.C:1790
drawCorrelationFunctionPsi.C:1791
drawCorrelationFunctionPsi.C:1792
drawCorrelationFunctionPsi.C:1793
drawCorrelationFunctionPsi.C:1794
drawCorrelationFunctionPsi.C:1795
drawCorrelationFunctionPsi.C:1796
drawCorrelationFunctionPsi.C:1797
drawCorrelationFunctionPsi.C:1798
drawCorrelationFunctionPsi.C:1799
drawCorrelationFunctionPsi.C:1800
drawCorrelationFunctionPsi.C:1801
drawCorrelationFunctionPsi.C:1802
drawCorrelationFunctionPsi.C:1803
drawCorrelationFunctionPsi.C:1804
drawCorrelationFunctionPsi.C:1805
drawCorrelationFunctionPsi.C:1806
drawCorrelationFunctionPsi.C:1807
drawCorrelationFunctionPsi.C:1808
drawCorrelationFunctionPsi.C:1809
drawCorrelationFunctionPsi.C:1810
drawCorrelationFunctionPsi.C:1811
drawCorrelationFunctionPsi.C:1812
drawCorrelationFunctionPsi.C:1813
drawCorrelationFunctionPsi.C:1814
drawCorrelationFunctionPsi.C:1815
drawCorrelationFunctionPsi.C:1816
drawCorrelationFunctionPsi.C:1817
drawCorrelationFunctionPsi.C:1818
drawCorrelationFunctionPsi.C:1819
drawCorrelationFunctionPsi.C:1820
drawCorrelationFunctionPsi.C:1821
drawCorrelationFunctionPsi.C:1822
drawCorrelationFunctionPsi.C:1823
drawCorrelationFunctionPsi.C:1824
drawCorrelationFunctionPsi.C:1825
drawCorrelationFunctionPsi.C:1826
drawCorrelationFunctionPsi.C:1827
drawCorrelationFunctionPsi.C:1828
drawCorrelationFunctionPsi.C:1829
drawCorrelationFunctionPsi.C:1830
drawCorrelationFunctionPsi.C:1831
drawCorrelationFunctionPsi.C:1832
drawCorrelationFunctionPsi.C:1833
drawCorrelationFunctionPsi.C:1834
drawCorrelationFunctionPsi.C:1835
drawCorrelationFunctionPsi.C:1836
drawCorrelationFunctionPsi.C:1837
drawCorrelationFunctionPsi.C:1838
drawCorrelationFunctionPsi.C:1839
drawCorrelationFunctionPsi.C:1840
drawCorrelationFunctionPsi.C:1841
drawCorrelationFunctionPsi.C:1842
drawCorrelationFunctionPsi.C:1843
drawCorrelationFunctionPsi.C:1844
drawCorrelationFunctionPsi.C:1845
drawCorrelationFunctionPsi.C:1846
drawCorrelationFunctionPsi.C:1847
drawCorrelationFunctionPsi.C:1848
drawCorrelationFunctionPsi.C:1849
drawCorrelationFunctionPsi.C:1850
drawCorrelationFunctionPsi.C:1851
drawCorrelationFunctionPsi.C:1852
drawCorrelationFunctionPsi.C:1853
drawCorrelationFunctionPsi.C:1854
drawCorrelationFunctionPsi.C:1855
drawCorrelationFunctionPsi.C:1856
drawCorrelationFunctionPsi.C:1857
drawCorrelationFunctionPsi.C:1858
drawCorrelationFunctionPsi.C:1859
drawCorrelationFunctionPsi.C:1860
drawCorrelationFunctionPsi.C:1861
drawCorrelationFunctionPsi.C:1862
drawCorrelationFunctionPsi.C:1863
drawCorrelationFunctionPsi.C:1864
drawCorrelationFunctionPsi.C:1865
drawCorrelationFunctionPsi.C:1866
drawCorrelationFunctionPsi.C:1867
drawCorrelationFunctionPsi.C:1868
drawCorrelationFunctionPsi.C:1869
drawCorrelationFunctionPsi.C:1870
drawCorrelationFunctionPsi.C:1871
drawCorrelationFunctionPsi.C:1872
drawCorrelationFunctionPsi.C:1873
drawCorrelationFunctionPsi.C:1874
drawCorrelationFunctionPsi.C:1875
drawCorrelationFunctionPsi.C:1876
drawCorrelationFunctionPsi.C:1877
drawCorrelationFunctionPsi.C:1878
drawCorrelationFunctionPsi.C:1879
drawCorrelationFunctionPsi.C:1880
drawCorrelationFunctionPsi.C:1881
drawCorrelationFunctionPsi.C:1882
drawCorrelationFunctionPsi.C:1883
drawCorrelationFunctionPsi.C:1884
drawCorrelationFunctionPsi.C:1885
drawCorrelationFunctionPsi.C:1886
drawCorrelationFunctionPsi.C:1887
drawCorrelationFunctionPsi.C:1888
drawCorrelationFunctionPsi.C:1889
drawCorrelationFunctionPsi.C:1890
drawCorrelationFunctionPsi.C:1891
drawCorrelationFunctionPsi.C:1892
drawCorrelationFunctionPsi.C:1893
drawCorrelationFunctionPsi.C:1894
drawCorrelationFunctionPsi.C:1895
drawCorrelationFunctionPsi.C:1896
drawCorrelationFunctionPsi.C:1897
drawCorrelationFunctionPsi.C:1898
drawCorrelationFunctionPsi.C:1899
drawCorrelationFunctionPsi.C:1900
drawCorrelationFunctionPsi.C:1901
drawCorrelationFunctionPsi.C:1902
drawCorrelationFunctionPsi.C:1903
drawCorrelationFunctionPsi.C:1904
drawCorrelationFunctionPsi.C:1905
drawCorrelationFunctionPsi.C:1906
drawCorrelationFunctionPsi.C:1907
drawCorrelationFunctionPsi.C:1908
drawCorrelationFunctionPsi.C:1909
drawCorrelationFunctionPsi.C:1910
drawCorrelationFunctionPsi.C:1911
drawCorrelationFunctionPsi.C:1912
drawCorrelationFunctionPsi.C:1913
drawCorrelationFunctionPsi.C:1914
drawCorrelationFunctionPsi.C:1915
drawCorrelationFunctionPsi.C:1916
drawCorrelationFunctionPsi.C:1917
drawCorrelationFunctionPsi.C:1918
drawCorrelationFunctionPsi.C:1919
drawCorrelationFunctionPsi.C:1920
drawCorrelationFunctionPsi.C:1921
drawCorrelationFunctionPsi.C:1922
drawCorrelationFunctionPsi.C:1923
drawCorrelationFunctionPsi.C:1924
drawCorrelationFunctionPsi.C:1925
drawCorrelationFunctionPsi.C:1926
drawCorrelationFunctionPsi.C:1927
drawCorrelationFunctionPsi.C:1928
drawCorrelationFunctionPsi.C:1929
drawCorrelationFunctionPsi.C:1930
drawCorrelationFunctionPsi.C:1931
drawCorrelationFunctionPsi.C:1932
drawCorrelationFunctionPsi.C:1933
drawCorrelationFunctionPsi.C:1934
drawCorrelationFunctionPsi.C:1935
drawCorrelationFunctionPsi.C:1936
drawCorrelationFunctionPsi.C:1937
drawCorrelationFunctionPsi.C:1938
drawCorrelationFunctionPsi.C:1939
drawCorrelationFunctionPsi.C:1940
drawCorrelationFunctionPsi.C:1941
drawCorrelationFunctionPsi.C:1942
drawCorrelationFunctionPsi.C:1943
drawCorrelationFunctionPsi.C:1944
drawCorrelationFunctionPsi.C:1945
drawCorrelationFunctionPsi.C:1946
drawCorrelationFunctionPsi.C:1947
drawCorrelationFunctionPsi.C:1948
drawCorrelationFunctionPsi.C:1949
drawCorrelationFunctionPsi.C:1950
drawCorrelationFunctionPsi.C:1951
drawCorrelationFunctionPsi.C:1952
drawCorrelationFunctionPsi.C:1953
drawCorrelationFunctionPsi.C:1954
drawCorrelationFunctionPsi.C:1955
drawCorrelationFunctionPsi.C:1956
drawCorrelationFunctionPsi.C:1957
drawCorrelationFunctionPsi.C:1958
drawCorrelationFunctionPsi.C:1959
drawCorrelationFunctionPsi.C:1960
drawCorrelationFunctionPsi.C:1961
drawCorrelationFunctionPsi.C:1962
drawCorrelationFunctionPsi.C:1963
drawCorrelationFunctionPsi.C:1964
drawCorrelationFunctionPsi.C:1965
drawCorrelationFunctionPsi.C:1966
drawCorrelationFunctionPsi.C:1967
drawCorrelationFunctionPsi.C:1968
drawCorrelationFunctionPsi.C:1969
drawCorrelationFunctionPsi.C:1970
drawCorrelationFunctionPsi.C:1971
drawCorrelationFunctionPsi.C:1972
drawCorrelationFunctionPsi.C:1973
drawCorrelationFunctionPsi.C:1974
drawCorrelationFunctionPsi.C:1975
drawCorrelationFunctionPsi.C:1976
drawCorrelationFunctionPsi.C:1977
drawCorrelationFunctionPsi.C:1978
drawCorrelationFunctionPsi.C:1979
drawCorrelationFunctionPsi.C:1980
drawCorrelationFunctionPsi.C:1981
drawCorrelationFunctionPsi.C:1982
drawCorrelationFunctionPsi.C:1983
drawCorrelationFunctionPsi.C:1984
drawCorrelationFunctionPsi.C:1985
drawCorrelationFunctionPsi.C:1986
drawCorrelationFunctionPsi.C:1987
drawCorrelationFunctionPsi.C:1988
drawCorrelationFunctionPsi.C:1989
drawCorrelationFunctionPsi.C:1990
drawCorrelationFunctionPsi.C:1991
drawCorrelationFunctionPsi.C:1992
drawCorrelationFunctionPsi.C:1993
drawCorrelationFunctionPsi.C:1994
drawCorrelationFunctionPsi.C:1995
drawCorrelationFunctionPsi.C:1996
drawCorrelationFunctionPsi.C:1997
drawCorrelationFunctionPsi.C:1998
drawCorrelationFunctionPsi.C:1999
drawCorrelationFunctionPsi.C:2000
drawCorrelationFunctionPsi.C:2001
drawCorrelationFunctionPsi.C:2002
drawCorrelationFunctionPsi.C:2003
drawCorrelationFunctionPsi.C:2004
drawCorrelationFunctionPsi.C:2005
drawCorrelationFunctionPsi.C:2006
drawCorrelationFunctionPsi.C:2007
drawCorrelationFunctionPsi.C:2008
drawCorrelationFunctionPsi.C:2009
drawCorrelationFunctionPsi.C:2010
drawCorrelationFunctionPsi.C:2011
drawCorrelationFunctionPsi.C:2012
drawCorrelationFunctionPsi.C:2013
drawCorrelationFunctionPsi.C:2014
drawCorrelationFunctionPsi.C:2015
drawCorrelationFunctionPsi.C:2016
drawCorrelationFunctionPsi.C:2017
drawCorrelationFunctionPsi.C:2018
drawCorrelationFunctionPsi.C:2019
drawCorrelationFunctionPsi.C:2020
drawCorrelationFunctionPsi.C:2021
drawCorrelationFunctionPsi.C:2022
drawCorrelationFunctionPsi.C:2023
drawCorrelationFunctionPsi.C:2024
drawCorrelationFunctionPsi.C:2025
drawCorrelationFunctionPsi.C:2026
drawCorrelationFunctionPsi.C:2027
drawCorrelationFunctionPsi.C:2028
drawCorrelationFunctionPsi.C:2029
drawCorrelationFunctionPsi.C:2030
drawCorrelationFunctionPsi.C:2031
drawCorrelationFunctionPsi.C:2032
drawCorrelationFunctionPsi.C:2033
drawCorrelationFunctionPsi.C:2034
drawCorrelationFunctionPsi.C:2035
drawCorrelationFunctionPsi.C:2036
drawCorrelationFunctionPsi.C:2037
drawCorrelationFunctionPsi.C:2038
drawCorrelationFunctionPsi.C:2039
drawCorrelationFunctionPsi.C:2040
drawCorrelationFunctionPsi.C:2041
drawCorrelationFunctionPsi.C:2042
drawCorrelationFunctionPsi.C:2043
drawCorrelationFunctionPsi.C:2044
drawCorrelationFunctionPsi.C:2045
drawCorrelationFunctionPsi.C:2046
drawCorrelationFunctionPsi.C:2047
drawCorrelationFunctionPsi.C:2048
drawCorrelationFunctionPsi.C:2049
drawCorrelationFunctionPsi.C:2050
drawCorrelationFunctionPsi.C:2051
drawCorrelationFunctionPsi.C:2052
drawCorrelationFunctionPsi.C:2053
drawCorrelationFunctionPsi.C:2054
drawCorrelationFunctionPsi.C:2055
drawCorrelationFunctionPsi.C:2056
drawCorrelationFunctionPsi.C:2057
drawCorrelationFunctionPsi.C:2058
drawCorrelationFunctionPsi.C:2059
drawCorrelationFunctionPsi.C:2060
drawCorrelationFunctionPsi.C:2061
drawCorrelationFunctionPsi.C:2062
drawCorrelationFunctionPsi.C:2063
drawCorrelationFunctionPsi.C:2064
drawCorrelationFunctionPsi.C:2065
drawCorrelationFunctionPsi.C:2066
drawCorrelationFunctionPsi.C:2067
drawCorrelationFunctionPsi.C:2068
drawCorrelationFunctionPsi.C:2069
drawCorrelationFunctionPsi.C:2070
drawCorrelationFunctionPsi.C:2071
drawCorrelationFunctionPsi.C:2072
drawCorrelationFunctionPsi.C:2073
drawCorrelationFunctionPsi.C:2074
drawCorrelationFunctionPsi.C:2075
drawCorrelationFunctionPsi.C:2076
drawCorrelationFunctionPsi.C:2077
drawCorrelationFunctionPsi.C:2078
drawCorrelationFunctionPsi.C:2079
drawCorrelationFunctionPsi.C:2080
drawCorrelationFunctionPsi.C:2081
drawCorrelationFunctionPsi.C:2082
drawCorrelationFunctionPsi.C:2083
drawCorrelationFunctionPsi.C:2084
drawCorrelationFunctionPsi.C:2085
drawCorrelationFunctionPsi.C:2086
drawCorrelationFunctionPsi.C:2087
drawCorrelationFunctionPsi.C:2088
drawCorrelationFunctionPsi.C:2089
drawCorrelationFunctionPsi.C:2090
drawCorrelationFunctionPsi.C:2091
drawCorrelationFunctionPsi.C:2092
drawCorrelationFunctionPsi.C:2093
drawCorrelationFunctionPsi.C:2094
drawCorrelationFunctionPsi.C:2095
drawCorrelationFunctionPsi.C:2096
drawCorrelationFunctionPsi.C:2097
drawCorrelationFunctionPsi.C:2098
drawCorrelationFunctionPsi.C:2099
drawCorrelationFunctionPsi.C:2100
drawCorrelationFunctionPsi.C:2101
drawCorrelationFunctionPsi.C:2102
drawCorrelationFunctionPsi.C:2103
drawCorrelationFunctionPsi.C:2104
drawCorrelationFunctionPsi.C:2105
drawCorrelationFunctionPsi.C:2106
drawCorrelationFunctionPsi.C:2107
drawCorrelationFunctionPsi.C:2108
drawCorrelationFunctionPsi.C:2109
drawCorrelationFunctionPsi.C:2110
drawCorrelationFunctionPsi.C:2111
drawCorrelationFunctionPsi.C:2112
drawCorrelationFunctionPsi.C:2113
drawCorrelationFunctionPsi.C:2114
drawCorrelationFunctionPsi.C:2115
drawCorrelationFunctionPsi.C:2116
drawCorrelationFunctionPsi.C:2117
drawCorrelationFunctionPsi.C:2118
drawCorrelationFunctionPsi.C:2119
drawCorrelationFunctionPsi.C:2120
drawCorrelationFunctionPsi.C:2121
drawCorrelationFunctionPsi.C:2122
drawCorrelationFunctionPsi.C:2123
drawCorrelationFunctionPsi.C:2124
drawCorrelationFunctionPsi.C:2125
drawCorrelationFunctionPsi.C:2126
drawCorrelationFunctionPsi.C:2127
drawCorrelationFunctionPsi.C:2128
drawCorrelationFunctionPsi.C:2129
drawCorrelationFunctionPsi.C:2130
drawCorrelationFunctionPsi.C:2131
drawCorrelationFunctionPsi.C:2132
drawCorrelationFunctionPsi.C:2133
drawCorrelationFunctionPsi.C:2134
drawCorrelationFunctionPsi.C:2135
drawCorrelationFunctionPsi.C:2136
drawCorrelationFunctionPsi.C:2137
drawCorrelationFunctionPsi.C:2138
drawCorrelationFunctionPsi.C:2139
drawCorrelationFunctionPsi.C:2140
drawCorrelationFunctionPsi.C:2141
drawCorrelationFunctionPsi.C:2142
drawCorrelationFunctionPsi.C:2143
drawCorrelationFunctionPsi.C:2144
drawCorrelationFunctionPsi.C:2145
drawCorrelationFunctionPsi.C:2146
drawCorrelationFunctionPsi.C:2147
drawCorrelationFunctionPsi.C:2148
drawCorrelationFunctionPsi.C:2149
drawCorrelationFunctionPsi.C:2150
drawCorrelationFunctionPsi.C:2151
drawCorrelationFunctionPsi.C:2152
drawCorrelationFunctionPsi.C:2153
drawCorrelationFunctionPsi.C:2154
drawCorrelationFunctionPsi.C:2155
drawCorrelationFunctionPsi.C:2156
drawCorrelationFunctionPsi.C:2157
drawCorrelationFunctionPsi.C:2158
drawCorrelationFunctionPsi.C:2159
drawCorrelationFunctionPsi.C:2160
drawCorrelationFunctionPsi.C:2161
drawCorrelationFunctionPsi.C:2162
drawCorrelationFunctionPsi.C:2163
drawCorrelationFunctionPsi.C:2164
drawCorrelationFunctionPsi.C:2165
drawCorrelationFunctionPsi.C:2166
drawCorrelationFunctionPsi.C:2167
drawCorrelationFunctionPsi.C:2168
drawCorrelationFunctionPsi.C:2169
drawCorrelationFunctionPsi.C:2170
drawCorrelationFunctionPsi.C:2171
drawCorrelationFunctionPsi.C:2172
drawCorrelationFunctionPsi.C:2173
drawCorrelationFunctionPsi.C:2174
drawCorrelationFunctionPsi.C:2175
drawCorrelationFunctionPsi.C:2176
drawCorrelationFunctionPsi.C:2177
drawCorrelationFunctionPsi.C:2178
drawCorrelationFunctionPsi.C:2179
drawCorrelationFunctionPsi.C:2180
drawCorrelationFunctionPsi.C:2181
drawCorrelationFunctionPsi.C:2182
drawCorrelationFunctionPsi.C:2183
drawCorrelationFunctionPsi.C:2184
drawCorrelationFunctionPsi.C:2185
drawCorrelationFunctionPsi.C:2186
drawCorrelationFunctionPsi.C:2187
drawCorrelationFunctionPsi.C:2188
drawCorrelationFunctionPsi.C:2189
drawCorrelationFunctionPsi.C:2190
drawCorrelationFunctionPsi.C:2191
drawCorrelationFunctionPsi.C:2192
drawCorrelationFunctionPsi.C:2193
drawCorrelationFunctionPsi.C:2194
drawCorrelationFunctionPsi.C:2195
drawCorrelationFunctionPsi.C:2196
drawCorrelationFunctionPsi.C:2197
drawCorrelationFunctionPsi.C:2198
drawCorrelationFunctionPsi.C:2199
drawCorrelationFunctionPsi.C:2200
drawCorrelationFunctionPsi.C:2201
drawCorrelationFunctionPsi.C:2202
drawCorrelationFunctionPsi.C:2203
drawCorrelationFunctionPsi.C:2204
drawCorrelationFunctionPsi.C:2205
drawCorrelationFunctionPsi.C:2206
drawCorrelationFunctionPsi.C:2207
drawCorrelationFunctionPsi.C:2208
drawCorrelationFunctionPsi.C:2209
drawCorrelationFunctionPsi.C:2210
drawCorrelationFunctionPsi.C:2211
drawCorrelationFunctionPsi.C:2212
drawCorrelationFunctionPsi.C:2213
drawCorrelationFunctionPsi.C:2214
drawCorrelationFunctionPsi.C:2215
drawCorrelationFunctionPsi.C:2216
drawCorrelationFunctionPsi.C:2217
drawCorrelationFunctionPsi.C:2218
drawCorrelationFunctionPsi.C:2219
drawCorrelationFunctionPsi.C:2220
drawCorrelationFunctionPsi.C:2221
drawCorrelationFunctionPsi.C:2222
drawCorrelationFunctionPsi.C:2223
drawCorrelationFunctionPsi.C:2224
drawCorrelationFunctionPsi.C:2225
drawCorrelationFunctionPsi.C:2226
drawCorrelationFunctionPsi.C:2227
drawCorrelationFunctionPsi.C:2228
drawCorrelationFunctionPsi.C:2229
drawCorrelationFunctionPsi.C:2230
drawCorrelationFunctionPsi.C:2231
drawCorrelationFunctionPsi.C:2232
drawCorrelationFunctionPsi.C:2233
drawCorrelationFunctionPsi.C:2234
drawCorrelationFunctionPsi.C:2235
drawCorrelationFunctionPsi.C:2236
drawCorrelationFunctionPsi.C:2237
drawCorrelationFunctionPsi.C:2238
drawCorrelationFunctionPsi.C:2239
drawCorrelationFunctionPsi.C:2240
drawCorrelationFunctionPsi.C:2241
drawCorrelationFunctionPsi.C:2242
drawCorrelationFunctionPsi.C:2243
drawCorrelationFunctionPsi.C:2244
drawCorrelationFunctionPsi.C:2245
drawCorrelationFunctionPsi.C:2246
drawCorrelationFunctionPsi.C:2247
drawCorrelationFunctionPsi.C:2248
drawCorrelationFunctionPsi.C:2249
drawCorrelationFunctionPsi.C:2250
drawCorrelationFunctionPsi.C:2251
drawCorrelationFunctionPsi.C:2252
drawCorrelationFunctionPsi.C:2253
drawCorrelationFunctionPsi.C:2254
drawCorrelationFunctionPsi.C:2255
drawCorrelationFunctionPsi.C:2256
drawCorrelationFunctionPsi.C:2257
drawCorrelationFunctionPsi.C:2258
drawCorrelationFunctionPsi.C:2259
drawCorrelationFunctionPsi.C:2260
drawCorrelationFunctionPsi.C:2261
drawCorrelationFunctionPsi.C:2262
drawCorrelationFunctionPsi.C:2263
drawCorrelationFunctionPsi.C:2264
drawCorrelationFunctionPsi.C:2265
drawCorrelationFunctionPsi.C:2266
drawCorrelationFunctionPsi.C:2267
drawCorrelationFunctionPsi.C:2268
drawCorrelationFunctionPsi.C:2269
drawCorrelationFunctionPsi.C:2270
drawCorrelationFunctionPsi.C:2271
drawCorrelationFunctionPsi.C:2272
drawCorrelationFunctionPsi.C:2273
drawCorrelationFunctionPsi.C:2274
drawCorrelationFunctionPsi.C:2275
drawCorrelationFunctionPsi.C:2276
drawCorrelationFunctionPsi.C:2277
drawCorrelationFunctionPsi.C:2278
drawCorrelationFunctionPsi.C:2279
drawCorrelationFunctionPsi.C:2280
drawCorrelationFunctionPsi.C:2281
drawCorrelationFunctionPsi.C:2282
drawCorrelationFunctionPsi.C:2283
drawCorrelationFunctionPsi.C:2284
drawCorrelationFunctionPsi.C:2285
drawCorrelationFunctionPsi.C:2286
drawCorrelationFunctionPsi.C:2287
drawCorrelationFunctionPsi.C:2288
drawCorrelationFunctionPsi.C:2289
drawCorrelationFunctionPsi.C:2290
drawCorrelationFunctionPsi.C:2291
drawCorrelationFunctionPsi.C:2292
drawCorrelationFunctionPsi.C:2293
drawCorrelationFunctionPsi.C:2294
drawCorrelationFunctionPsi.C:2295
drawCorrelationFunctionPsi.C:2296
drawCorrelationFunctionPsi.C:2297
drawCorrelationFunctionPsi.C:2298
drawCorrelationFunctionPsi.C:2299
drawCorrelationFunctionPsi.C:2300
drawCorrelationFunctionPsi.C:2301
drawCorrelationFunctionPsi.C:2302
drawCorrelationFunctionPsi.C:2303
drawCorrelationFunctionPsi.C:2304
drawCorrelationFunctionPsi.C:2305
drawCorrelationFunctionPsi.C:2306
drawCorrelationFunctionPsi.C:2307
drawCorrelationFunctionPsi.C:2308
drawCorrelationFunctionPsi.C:2309
drawCorrelationFunctionPsi.C:2310
drawCorrelationFunctionPsi.C:2311
drawCorrelationFunctionPsi.C:2312
drawCorrelationFunctionPsi.C:2313
drawCorrelationFunctionPsi.C:2314
drawCorrelationFunctionPsi.C:2315
drawCorrelationFunctionPsi.C:2316
drawCorrelationFunctionPsi.C:2317
drawCorrelationFunctionPsi.C:2318
drawCorrelationFunctionPsi.C:2319
drawCorrelationFunctionPsi.C:2320
drawCorrelationFunctionPsi.C:2321
drawCorrelationFunctionPsi.C:2322
drawCorrelationFunctionPsi.C:2323
drawCorrelationFunctionPsi.C:2324
drawCorrelationFunctionPsi.C:2325
drawCorrelationFunctionPsi.C:2326
drawCorrelationFunctionPsi.C:2327
drawCorrelationFunctionPsi.C:2328
drawCorrelationFunctionPsi.C:2329
drawCorrelationFunctionPsi.C:2330
drawCorrelationFunctionPsi.C:2331
drawCorrelationFunctionPsi.C:2332
drawCorrelationFunctionPsi.C:2333
drawCorrelationFunctionPsi.C:2334
drawCorrelationFunctionPsi.C:2335
drawCorrelationFunctionPsi.C:2336
drawCorrelationFunctionPsi.C:2337
drawCorrelationFunctionPsi.C:2338
drawCorrelationFunctionPsi.C:2339
drawCorrelationFunctionPsi.C:2340
drawCorrelationFunctionPsi.C:2341
drawCorrelationFunctionPsi.C:2342
drawCorrelationFunctionPsi.C:2343
drawCorrelationFunctionPsi.C:2344
drawCorrelationFunctionPsi.C:2345
drawCorrelationFunctionPsi.C:2346
drawCorrelationFunctionPsi.C:2347
drawCorrelationFunctionPsi.C:2348
drawCorrelationFunctionPsi.C:2349
drawCorrelationFunctionPsi.C:2350
drawCorrelationFunctionPsi.C:2351
drawCorrelationFunctionPsi.C:2352
drawCorrelationFunctionPsi.C:2353
drawCorrelationFunctionPsi.C:2354
drawCorrelationFunctionPsi.C:2355
drawCorrelationFunctionPsi.C:2356
drawCorrelationFunctionPsi.C:2357
drawCorrelationFunctionPsi.C:2358
drawCorrelationFunctionPsi.C:2359
drawCorrelationFunctionPsi.C:2360
drawCorrelationFunctionPsi.C:2361
drawCorrelationFunctionPsi.C:2362
drawCorrelationFunctionPsi.C:2363
drawCorrelationFunctionPsi.C:2364
drawCorrelationFunctionPsi.C:2365
drawCorrelationFunctionPsi.C:2366
drawCorrelationFunctionPsi.C:2367
drawCorrelationFunctionPsi.C:2368
drawCorrelationFunctionPsi.C:2369
drawCorrelationFunctionPsi.C:2370
drawCorrelationFunctionPsi.C:2371
drawCorrelationFunctionPsi.C:2372
drawCorrelationFunctionPsi.C:2373
drawCorrelationFunctionPsi.C:2374
drawCorrelationFunctionPsi.C:2375
drawCorrelationFunctionPsi.C:2376
drawCorrelationFunctionPsi.C:2377
drawCorrelationFunctionPsi.C:2378
drawCorrelationFunctionPsi.C:2379
drawCorrelationFunctionPsi.C:2380
drawCorrelationFunctionPsi.C:2381