From 0bf9d3e48dabef4d01aa5aa9d9e8b07a61918155 Mon Sep 17 00:00:00 2001 From: mfloris Date: Thu, 8 May 2014 15:10:51 +0200 Subject: [PATCH] Added interpolation for the 20-40% bin --- PWGLF/ThermalFits/FitNPartDependence.C | 40 +++++++++---------- .../ThermalFits/InterpolateRatiosAndYields.C | 28 +++++++++++-- 2 files changed, 44 insertions(+), 24 deletions(-) diff --git a/PWGLF/ThermalFits/FitNPartDependence.C b/PWGLF/ThermalFits/FitNPartDependence.C index 8e7d79b3b31..5e6089337f1 100644 --- a/PWGLF/ThermalFits/FitNPartDependence.C +++ b/PWGLF/ThermalFits/FitNPartDependence.C @@ -88,28 +88,28 @@ void FitNPartDependence() { // const char * centrToExtrapolate = "V0M0010"; // Int_t pdg = 1000010020; // TClonesArray * arr = AliParticleYield::ReadFromASCIIFile("PbPb_2760_DeuHelium3.txt"); - // Deuteron pPb - // centrFile = "dndeta_pPb.txt"; - // const char * centralityToPlot[] = { "V0A0010", "V0A1020", "V0A2040", "V0A4060", "V0A6000" ,0}; - // // const char * centrToExtrapolate = "V0A0005"; - // const char * centrToExtrapolate = "V0A6080"; - // Int_t pdg = 1000010020; - // TClonesArray * arr = AliParticleYield::ReadFromASCIIFile("pPb_5020_deuteron.txt"); - // maxy = 0.01; - // systemAndEnergy = "p-Pb #sqrt{#it{s}}_{NN} = 5.02 TeV"; - // energy = 5020; - // collSystem = 1; - // K* pPb - centrFile = "dndeta_pPb.txt"; - // centrFile = "dndeta_PbPb.txt"; - maxy = 1; - systemAndEnergy = "p-Pb #sqrt{#it{s}}_{NN} = 5.02 TeV"; - const char * centralityToPlot[] = { "V0A0020" , "V0A2040" , "V0A4060" , "V0A6080", "V0A8000" , 0}; + // Deuteron pPb + centrFile = "dndeta_pPb.txt"; + const char * centralityToPlot[] = { "V0A0010", "V0A1020", "V0A2040", "V0A4060", "V0A6000" ,0}; const char * centrToExtrapolate = "V0A0005"; - Int_t pdg = 313; - TClonesArray * arr = AliParticleYield::ReadFromASCIIFile("./pPb_5020_Kstar.txt"); + //const char * centrToExtrapolate = "V0A6080"; + Int_t pdg = -1000010020; + TClonesArray * arr = AliParticleYield::ReadFromASCIIFile("pPb_5020_deuteron.txt"); + maxy = 0.01; + systemAndEnergy = "p-Pb #sqrt{#it{s}}_{NN} = 5.02 TeV"; energy = 5020; collSystem = 1; + // K* pPb + // centrFile = "dndeta_pPb.txt"; + // // centrFile = "dndeta_PbPb.txt"; + // maxy = 1; + // systemAndEnergy = "p-Pb #sqrt{#it{s}}_{NN} = 5.02 TeV"; + // const char * centralityToPlot[] = { "V0A0020" , "V0A2040" , "V0A4060" , "V0A6080", "V0A8000" , 0}; + // const char * centrToExtrapolate = "V0A0005"; + // Int_t pdg = 313; + // TClonesArray * arr = AliParticleYield::ReadFromASCIIFile("./pPb_5020_Kstar.txt"); + // energy = 5020; + // collSystem = 1; // Helium3 @@ -181,8 +181,6 @@ void FitNPartDependence() { Double_t errorSystPlus = TMath::Abs(FitShiftedGraphAndExtrapolate(grSyst, kShiftHarder, f1, centrToExtrapolate, kRed)-yield); Double_t errorSystMinus = TMath::Abs(FitShiftedGraphAndExtrapolate(grSyst, kShiftSofter, f1, centrToExtrapolate, kRed)-yield); - // Double_t errorStatPlus = FitShiftedGraphAndExtrapolate(grStat, kShiftUp , f1, centrToExtrapolate, kBlue) -yield; - // Double_t errorStatMinus = FitShiftedGraphAndExtrapolate(grStat, kShiftDown, f1, centrToExtrapolate, kBlue) -yield; Double_t errorStat = fError->Eval(npartPbPb[centrToExtrapolate]); diff --git a/PWGLF/ThermalFits/InterpolateRatiosAndYields.C b/PWGLF/ThermalFits/InterpolateRatiosAndYields.C index ffe0922daa6..324ce73e158 100644 --- a/PWGLF/ThermalFits/InterpolateRatiosAndYields.C +++ b/PWGLF/ThermalFits/InterpolateRatiosAndYields.C @@ -8,6 +8,7 @@ TClonesArray * arr =0; void InterpolateRatios(Int_t pdg1, Int_t pdg2, TString centr1="V0M0005", TString centr2="V0M0510", TString centrfinal="V0M0010") ; void Interpolate0010(Int_t pdg) ; +void Interpolate2040(Int_t pdg) ; void Interpolate6080(Int_t pdg) ; void ExtrapolateWithConstantRatioToPions(Int_t pdg, TString centrOrigin, TString centrDest); @@ -38,14 +39,21 @@ void InterpolateRatiosAndYields() { // InterpolateRatios(2212,211, "V0M6070", "V0M7080", "V0M6080"); // InterpolateRatios(321,211, "V0M6070", "V0M7080", "V0M6080"); + Interpolate2040(211); + Interpolate2040(-211); + Interpolate2040(2212); + Interpolate2040(-2212); + Interpolate2040(321); + Interpolate2040(-321); + // *************** Lambda and K0 ***************** // arr= AliParticleYield::ReadFromASCIIFile("PbPb_2760_LambdaK0.txt"); // Interpolate0010(3122); // Interpolate0010(310); // *************** Helium 3 ***************** - arr = AliParticleYield::ReadFromASCIIFile("PbPb_2760_DeuHelium3.txt"); - arr->AbsorbObjects(AliParticleYield::ReadFromASCIIFile("./PbPb_2760_AveragedNumbers.txt")); - ExtrapolateWithConstantRatioToPions(1000020030, "V0M0020", "V0M0010"); + // arr = AliParticleYield::ReadFromASCIIFile("PbPb_2760_DeuHelium3.txt"); + // arr->AbsorbObjects(AliParticleYield::ReadFromASCIIFile("./PbPb_2760_AveragedNumbers.txt")); + // ExtrapolateWithConstantRatioToPions(1000020030, "V0M0020", "V0M0010"); // *************** Kstar ***************** // arr = AliParticleYield::ReadFromASCIIFile("PbPb_2760_Kstar892.txt"); // arr->AbsorbObjects(AliParticleYield::ReadFromASCIIFile("./PbPb_2760_AveragedNumbers.txt")); @@ -104,6 +112,20 @@ void Interpolate6080(Int_t pdg) { std::cout << pav->GetYield() << ", " << pav->GetStatError() << ", "<< pav->GetSystError() << std::endl; +} + +void Interpolate2040(Int_t pdg) { + + TString centrPrefix = collSystem == 2 ? "V0M" : "V0A"; + + AliParticleYield * p0 = AliParticleYield::FindParticle(arr, pdg, collSystem, energy, centrPrefix+"2030"); + AliParticleYield * p1 = AliParticleYield::FindParticle(arr, pdg, collSystem, energy, centrPrefix+"3040"); + p0->Scale(0.5); + p1->Scale(0.5); + AliParticleYield * pav = AliParticleYield::Add(p0,p1); + std::cout << pav->GetYield() << ", " << pav->GetStatError() << ", "<< pav->GetSystError() << std::endl; + + } void InterpolateRatios(Int_t pdg1, Int_t pdg2, TString centr1, TString centr2, TString centrfinal) { -- 2.43.0