#include <math.h> #include <time.h> #include <stdio.h> #include <stdlib.h> #include <Riostream.h> #include <assert.h> #include "TVector2.h" #include "TFile.h" #include "TString.h" #include "TF1.h" #include "TF3.h" #include "TH1.h" #include "TH2.h" #include "TH3.h" #include "TProfile.h" #include "TProfile2D.h" #include "TMath.h" #include "TText.h" #include "TRandom3.h" #include "TArray.h" #include "TLegend.h" #include "TStyle.h" #include "TMinuit.h" #include "TCanvas.h" #include "TLatex.h" #include "TPaveStats.h" #include "TASImage.h" #include "TSpline.h" #define BohrR 1963.6885 // Bohr Radius for pions #define FmToGeV 0.19733 // conversion to fm #define PI 3.1415926 #define masspiC 0.1395702 // pi+ mass (GeV/c^2) double kappa3 = 0.1; // 0.1 (default), double kappa4 = 0.5; // 0.5 (default), using namespace std; int CollisionType_def=0;// 0=PbPb, 1=pPb, 2=pp bool MCcase_def=0;// MC data? int CHARGE_def=-1;// -1 or +1: + or - pions for same-charge case, --+ or ++- for mixed-charge case bool SameCharge_def=kTRUE;// 3-pion same-charge? bool AddCC=kTRUE; bool Gaussian=1;// Gaussian or Exponential? bool UseC2Bkg=1;// use Pythia/DPMJET bkg? // bool MuonCorrection=1;// apply Muon correction? bool IncludeExpansion=kTRUE;// Include EdgeWorth coefficients? bool FixExpansionAvg=1; int Mbin_def=2;// 0-19 (0=1050-2000 pions, 19=0-5 pions) int EDbin_def=0;// 0-2: Kt3 bin int Ktbin_def=1;// 1-6, 10=full range double MRCShift=1.0;// 1.0=full Momentum Resolution Correction. 1.1 for 10% systematic increase bool FitRangeShift=0;// 30% reduction in Fit Range bool ExtremeFitRangeC2=0;// 1.2 GeV/c fit range for C2? // // bool SaveToFile_def=kFALSE;// Save outputs to file? int SourceType=0;// 0=Therminator, 1=Therminator with 50fm cut (keep at 0 for default), 2=Gaussian bool ConstantFSI=kFALSE;// Constant FSI's for each kt bin? bool GofP=kFALSE;// Include momentum dependence of coherent fraction? bool ChargeConstraint=kFALSE;// Include Charge Constraint for coherent states? bool LinkN=kFALSE;// Make N(++)=N(+-)? bool GeneratedSignal=kFALSE;// Was the QS+FSI signal generated in MC? bool Tricubic=kFALSE;// Tricubic or Trilinear interpolation? Tricubic was found to be unstable. bool IncludeSS=kTRUE;// Include same-charge two-pion correlations in the fit? bool IncludeOS=kFALSE;// Include mixed-charge two-pion correlations in the fit? // // double ExpPower=2;// default Gaussian(2), Exponential(1) // // const int Sbins_2=1;// 2-particle Species bins. 1=Pion-Pion only. max=6 (pi-pi, pi-k, pi-p, k-p, k-k, p-p) const int Sbins_3=1;// 3-particle Species bins. 1=Pion-Pion-Pion only. max=10 const int fitlimitLowBin = 3; const int fitlimitHighBin = 14;// max = 20 (100MeV) bool LHC11h=kTRUE;// always kTRUE bool ExplicitLoop=kFALSE;// always kFALSE bool PairCut=kTRUE;// always kTRUE int FSIindex; int Ktlowbin; int Kthighbin; int MomResCentBin;// momentum resolution centrality bin float TwoFrac;// Lambda parameter float OneFrac;// Lambda^{1/2} float ThreeFrac;// Lambda^{3/2} float Q2Limit; float Q3Limit; double lambda_PM = 0.7; TH1D *CoulCorr2SS; TH1D *CoulCorr2OS; TH1D *CoulCorr2SS_2;// for interpolation in transition from Therminator to Gauss K factors TH1D *CoulCorr2OS_2;// for interpolation in transition from Therminator to Gauss K factors // void ReadCoulCorrections(int); void ReadMomResFile(int); void ReadMuonCorrections(int,int); double CoulCorr2(int, double); double CoulCorrGRS(bool, double, double, double); double Dfitfunction_c3(Double_t *x, Double_t *par); double Gamov(int, double); double C2ssFitFunction(Double_t *x, Double_t *par); double C2osFitFunction(Double_t *x, Double_t *par); double cubicInterpolate(double[4], double); double bicubicInterpolate(double[4][4], double, double); double tricubicInterpolate(double[4][4][4], double, double, double); void DrawALICELogo(Bool_t, Float_t, Float_t, Float_t, Float_t); // void fcnC2_Global(int&, double*, double&, double[], int); const int BINRANGE_C2global=400;// 40 double C2ssFitting[BINRANGE_C2global]; double C2osFitting[BINRANGE_C2global]; double C2ssFitting_e[BINRANGE_C2global]; double C2osFitting_e[BINRANGE_C2global]; double K2SS[BINRANGE_C2global]; double K2OS[BINRANGE_C2global]; double K2SS_e[BINRANGE_C2global]; double K2OS_e[BINRANGE_C2global]; double Chi2_C2global; double NFitPoints_C2global; double Chi2_C3global; double NFitPoints_C3global; const int BINS_OSL=40;// out,side,long bins double A[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double B[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double A_e[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double B_e[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double avg_q[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double K_OSL[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double K_OSL_e[BINS_OSL][BINS_OSL][BINS_OSL];// out,side,long double Chi2_OSL; double NFitPoints_OSL; const int BINRANGE_3=40;// q12,q13,q23 int BINLIMIT_3=20; double A_3[BINRANGE_3][BINRANGE_3][BINRANGE_3];// q12,q13,q23 double A_3_e[BINRANGE_3][BINRANGE_3][BINRANGE_3];// q12,q13,q23 double B_3[BINRANGE_3][BINRANGE_3][BINRANGE_3];// q12,q13,q23 double BinCenters[400]; double BinWidthQ2; double Chi2_c3; double NFitPoints_c3; void fcn_c3(int&, double*, double&, double[], int); double C3_base[BINRANGE_3][BINRANGE_3][BINRANGE_3]; double N_term1[BINRANGE_3][BINRANGE_3][BINRANGE_3]; double N_term2[BINRANGE_3][BINRANGE_3][BINRANGE_3]; double N_term3[BINRANGE_3][BINRANGE_3][BINRANGE_3]; double N_term4[BINRANGE_3][BINRANGE_3][BINRANGE_3]; double N_term5[BINRANGE_3][BINRANGE_3][BINRANGE_3]; TH3D *MomRes3d[2][5];// SC/MC, term# TH1D *MomRes1d[2][5];// SC/MC, term# TH1D *MomResC2[2];// SC/MC int MRC2index;// index for C2 MRC double AvgP[6][20];// kt bin, qinv bin TH1D *C2muonCorrection_SC; TH1D *C2muonCorrection_MC; TH1D *C3muonCorrection; TF1 *StrongSC;// same-charge pion strong FSI TF1 *MixedChargeSysFit;// mixed-charge 3-pion cumulant residue obtained from Plot_plotsTPR.C TF1 *BkgMCFit;// Pythia or DPMJET fit // void Plot_PDCumulants_TPR(bool SaveToFile=SaveToFile_def, int CollisionType=CollisionType_def, bool MCcase=MCcase_def, bool SameCharge=SameCharge_def, int EDbin=EDbin_def, int CHARGE=CHARGE_def, int Mbin=Mbin_def, int Ktbin=Ktbin_def){ EDbin_def=EDbin; Ktbin_def=Ktbin; CollisionType_def=CollisionType; SaveToFile_def=SaveToFile; MCcase_def=MCcase; CHARGE_def=CHARGE; SameCharge_def=SameCharge;// 3-pion same-charge Mbin_def=Mbin; // if(CollisionType==0){ Q3Limit = 0.1 + 0.1*Mbin/16.; }else{ Q3Limit = 0.3 + 0.2*fabs(Mbin-10)/9.; } if(FitRangeShift) Q3Limit -= 0.3*Q3Limit; Q2Limit = Q3Limit/sqrt(2.); if(ExtremeFitRangeC2) Q2Limit = 1.2; //Q2Limit = (0.3 + 0.2*fabs(Mbin-10)/9.) / sqrt(2.); //cout<<Q2Limit<<endl; // PbPb then pPb then pp //double Nch_means[3][20]={{1828,1407,1071,896,757,619,496,402,324,252,171,117,83,63,50,36, 36,36,36,36}, {71,71,71,71,71,71,71,71,71,71,71,71, 71,57,45,32,23,16,9.7,5}, {47,47,47,47,47,47,47,47,47,47,47,47,47, 47,37,27,20,15,8.6,4.6}}; //kappa3 = 0.1; //kappa4 = 0.5; //kappa3 = 0.05 + 0.25*(Mbin/19.);// linear dependence of kappa3 //double MainFitParams[8]={0}; // //double RefMainFitParams[8]={6.71377, 5.39974, 1.3822, 1.82065, 6.71377, 5.39974, 1.3822, 1.82065}; if(Gaussian) ExpPower=2.; else ExpPower=1.; if(CollisionType==0){ if(Mbin==0) FSIindex = 0;// 0-5% Therm else if(Mbin==1) FSIindex = 1;// 0-10% Therm else if(Mbin<=3) FSIindex = 2;// 10-20% Therm else if(Mbin<=5) FSIindex = 3;// 20-30% Therm else if(Mbin<=7) FSIindex = 4;// 30-40% Therm else if(Mbin<=9) FSIindex = 5;// 40-50% Therm else FSIindex = 6;// EW R=4.0 fm // next bin is 3.0 fm which is used for interpolation }else{ if(Mbin<=14) FSIindex = 8;// EW R=2.5 fm else if(Mbin<=16) FSIindex = 9;// EW R=2.0 fm else FSIindex = 10;// EW R=1.7 fm // last bin is 1.25 fm which is used for interpolation } // chooses lambda=0.7 for R=10,8,6,4,2 if(Mbin<=2) MRC2index=138; else if(Mbin<=5) MRC2index=106; else if(Mbin<=9) MRC2index=74; else if(Mbin<=13) MRC2index=42;// was 9 by mistake, now 13 (corrected on 22/10/2013) else MRC2index=10; TwoFrac = lambda_PM; OneFrac = sqrt(TwoFrac); ThreeFrac = pow(TwoFrac,3/2.); if(CollisionType==0 && Mbin<10) BINLIMIT_3=20; else BINLIMIT_3=40; Ktlowbin=(Ktbin)*2+3;// kt bins are 0.5 GeV/c wide (0-0.5, 0.5-1.0 ...) Kthighbin=(Ktbin)*2+4; if(Ktbin==2) {cout<<"Kthighbin increased"<<endl; Kthighbin = 20;}// to make <KT3> comparible to <kT> BkgMCFit=new TF1("BkgMCFit","[0]*([1] + [2]*exp(-pow([3]*x,2)))",0,2.0); BkgMCFit->SetParameter(0,1); BkgMCFit->SetParameter(1,1.001); BkgMCFit->SetParameter(2,0.04102); BkgMCFit->SetParameter(3,1.874); BkgMCFit->SetLineColor(1); if(!MCcase){ // kt bin(0.2<kt<0.3, 0.3<kt<1.0), Mult bin(low to high) float paramsPP_0[2][10]={{1.00014,1.00019,9.99733e-01,9.99379e-01,9.99488e-01, 9.99488e-01,9.99488e-01,9.99488e-01,9.99488e-01,9.99488e-01},{9.97267e-01,9.98918e-01,9.99117e-01,9.98878e-01,9.99457e-01, 9.99457e-01,9.99457e-01,9.99457e-01,9.99457e-01,9.99457e-01}}; float paramsPP_1[2][10]={{1.00198,1.00156,1.00109,1.00068,1.00079, 1.00079,1.00079,1.00079,1.00079,1.00079}, {9.98091e-01,9.98893e-01,9.99782e-01,9.99595e-01,1.00037, 1.00037,1.00037,1.00037,1.00037,1.00037}}; float paramsPP_2[2][10]={{4.51532e-02,4.30415e-02,4.10306e-02,3.42739e-02,2.95607e-02, 2.95607e-02,2.95607e-02,2.95607e-02,2.95607e-02,2.95607e-02}, {9.89882e-02,1.00657e-01,8.93560e-02,7.55225e-02,6.52689e-02, 6.52689e-02,6.52689e-02,6.52689e-02,6.52689e-02,6.52689e-02}}; float paramsPP_3[2][10]={{1.57172,1.91807,1.87353,1.88549,2.22286, 2.22286,2.22286,2.22286,2.22286,2.22286}, {1.49834,1.65562,1.66501,1.75097,1.78597, 1.78597,1.78597,1.78597,1.78597,1.78597}}; ////////////// float paramsPPb_0[2][10]={{9.99479e-01,9.98772e-01,9.99295e-01,9.99532e-01,9.99430e-01,9.99642e-01,9.99424e-01, 9.99424e-01,9.99424e-01,9.99424e-01}, {9.91941e-01,9.96569e-01,9.97183e-01,9.98106e-01,9.98591e-01,9.98717e-01,9.98517e-01, 9.98517e-01,9.98517e-01,9.98517e-01}}; float paramsPPb_1[2][10]={{1.00093,1.00023,1.00068,1.00082,1.00077,1.00087,1.00069, 1.00069,1.00069,1.00069}, {9.93727e-01,9.98118e-01,9.98630e-01,9.99584e-01,1.00005,1.00017,9.99972e-01, 9.99972e-01,9.99972e-01,9.99972e-01}}; float paramsPPb_2[2][10]={{3.05502e-02,2.21034e-02,1.39907e-02,9.82480e-03,6.48672e-03,5.37259e-03,4.54498e-03, 4.54498e-03,4.54498e-03,4.54498e-03}, {8.86336e-02,5.53244e-02,3.59877e-02,2.57145e-02,1.82072e-02,1.44253e-02,1.28923e-02, 1.28923e-02,1.28923e-02, 1.28923e-02}}; float paramsPPb_3[2][10]={{1.69217,1.51289,1.63134,1.71567,1.71918,1.78037,2.10736, 2.10736,2.10736,2.10736}, {1.21996,1.34570,1.35660,1.37801,1.39951,1.39093,1.33747, 1.33747,1.33747,1.33747}}; int ktindexBkg = Ktbin - 1; if(ktindexBkg > 1) ktindexBkg = 1; if(CollisionType==2){ BkgMCFit->FixParameter(0, paramsPP_0[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(1, paramsPP_1[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(2, paramsPP_2[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(3, paramsPP_3[ktindexBkg][19-Mbin]); }else if(CollisionType==1){ BkgMCFit->FixParameter(0, paramsPPb_0[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(1, paramsPPb_1[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(2, paramsPPb_2[ktindexBkg][19-Mbin]); BkgMCFit->FixParameter(3, paramsPPb_3[ktindexBkg][19-Mbin]); }else{// no Bkg for PbPb BkgMCFit->FixParameter(0, 1); BkgMCFit->FixParameter(1, 1); BkgMCFit->FixParameter(2, 0); BkgMCFit->FixParameter(3, 1); } } // Core-Halo modeling of single-particle and triple-particle core fraction //float AvgN[10]={1.99266, 3.97789, 6.4624, 8.94042, 11.4194, 13.8987, 16.385, 18.8756, 21.3691, 26.0742};// 13c (avg total pion mult)/2. 2.0sigma // // bin centers from QS+FSI double BinCenterPbPbCentral[40]={0.00206385, 0.00818515, 0.0129022, 0.0177584, 0.0226881, 0.027647, 0.032622, 0.0376015, 0.042588, 0.0475767, 0.0525692, 0.0575625, 0.0625569, 0.0675511, 0.0725471, 0.0775436, 0.0825399, 0.0875364, 0.0925339, 0.0975321, 0.102529, 0.107527, 0.112525, 0.117523, 0.122522, 0.12752, 0.132519, 0.137518, 0.142516, 0.147515, 0.152514, 0.157513, 0.162513, 0.167512, 0.172511, 0.177511, 0.18251, 0.187509, 0.192509, 0.197509}; double BinCenterPbPbPeripheral[40]={0.00206385, 0.00818515, 0.0129022, 0.0177584, 0.0226881, 0.027647, 0.032622, 0.0376015, 0.042588, 0.0475767, 0.0525692, 0.0575625, 0.0625569, 0.0675511, 0.0725471, 0.0775436, 0.0825399, 0.0875364, 0.0925339, 0.0975321, 0.102529, 0.107527, 0.112525, 0.117523, 0.122522, 0.12752, 0.132519, 0.137518, 0.142516, 0.147515, 0.152514, 0.157513, 0.162513, 0.167512, 0.172511, 0.177511, 0.18251, 0.187509, 0.192509, 0.197509}; double BinCenterpPbAndpp[40]={0.00359275, 0.016357, 0.0257109, 0.035445, 0.045297, 0.0552251, 0.0651888, 0.0751397, 0.0851088, 0.0951108, 0.105084, 0.115079, 0.12507, 0.135059, 0.145053, 0.155049, 0.16505, 0.175038, 0.185039, 0.195034, 0.205023, 0.215027, 0.225024, 0.235023, 0.245011, 0.255017, 0.265017, 0.275021, 0.285021, 0.295017, 0.305018, 0.315018, 0.325013, 0.335011, 0.345016, 0.355019, 0.365012, 0.375016, 0.385017, 0.395016}; if(CollisionType==0){ for(int i=0; i<40; i++){ if(Mbin<10) BinCenters[i] = BinCenterPbPbCentral[i]; else BinCenters[i] = BinCenterPbPbPeripheral[i]; } BinWidthQ2 = 0.005; }else{ for(int i=0; i<40; i++){ BinCenters[i] = BinCenterpPbAndpp[i]; } BinWidthQ2 = 0.01; } // extend BinCenters for high q for(int index=40; index<400; index++){ if(CollisionType==0) BinCenters[index] = (index+0.5)*(0.005); else BinCenters[index] = (index+0.5)*(0.010); } // Set 0's to 3-particle fit arrays for(int i=1; i<=BINLIMIT_3; i++){// bin number for(int j=1; j<=BINLIMIT_3; j++){// bin number for(int k=1; k<=BINLIMIT_3; k++){// bin number A_3[i-1][j-1][k-1]=0; A_3_e[i-1][j-1][k-1]=0; B_3[i-1][j-1][k-1]=0; } } } // same-charge pion strong FSI (obtained from ratio of CoulStrong to Coul at R=7 Gaussian from Lednicky's code) StrongSC=new TF1("StrongSC","[0]+[1]*exp(-pow([2]*x,2))",0,1000); StrongSC->FixParameter(0,9.99192e-01); StrongSC->FixParameter(1,-8.01113e-03); StrongSC->FixParameter(2,5.35912e-02); // mixed-charge 3-pion cumulant residue obtained from Plot_plotsTPR.C MixedChargeSysFit=new TF1("MixedChargeSysFit","[0]+[1]*exp(-pow([2]*x/0.19733,2))",0,.5); MixedChargeSysFit->FixParameter(0,1); MixedChargeSysFit->FixParameter(1,.06); MixedChargeSysFit->FixParameter(2,5.5 - 6*(Mbin/19.)); gStyle->SetOptStat(0); gStyle->SetOptDate(0); gStyle->SetOptFit(0111); TFile *_file0; if(CollisionType==0){// PbPb if(MCcase) { if(Mbin<=1){ _file0 = new TFile("Results/PDC_12a17a_R10_2.root","READ"); }else{ _file0 = new TFile("Results/PDC_12a17e_R10.root","READ"); } }else{ //_file0 = new TFile("Results/PDC_10h_11h_0to50_50to100.root","READ"); //_file0 = new TFile("Results/PDC_10h_11h_0to50_50to100_3Ktbins.root","READ"); //_file0 = new TFile("Results/PDC_10h_2Kt3bins.root","READ"); //_file0 = new TFile("Results/PDC_10h_NewNorm.root","READ"); //_file0 = new TFile("Results/PDC_10h_dEta0p03dPhi0p03.root","READ"); //_file0 = new TFile("Results/PDC_10h_noPID.root","READ"); //_file0 = new TFile("Results/PDC_10h_1percentCentral.root","READ"); //_file0 = new TFile("Results/PDC_10h_11h_2Kt3bins.root","READ");// v5 and before //_file0 = new TFile("Results/PDC_10h_dEta0p025dPhi0p055.root","READ"); //_file0 = new TFile("Results/PDC_11h_3Kt3bins_FB5and7overlap.root","READ"); //_file0 = new TFile("Results/PDC_10h_11h_V0binning.root","READ"); //_file0 = new TFile("Results/PDC_10h_11h_lowerMultBins.root","READ"); _file0 = new TFile("Results/PDC_10h_11h_NclsFix.root","READ");// standard } }else if(CollisionType==1){// pPb if(!MCcase){ //_file0 = new TFile("Results/PDC_13c_kMB.root","READ"); //_file0 = new TFile("Results/PDC_13bc.root","READ"); //_file0 = new TFile("Results/PDC_13bc_3Ktbins.root","READ"); //_file0 = new TFile("Results/PDC_13bc_kAnyINT_kINT7.root","READ"); //_file0 = new TFile("Results/PDC_13bc_2Kt3bins.root","READ"); //_file0 = new TFile("Results/PDC_13bc_kINT7.root","READ");// paper v1 //_file0 = new TFile("Results/PDC_13c_posEta.root","READ");// eta sub-range //_file0 = new TFile("Results/PDC_13c_negEta.root","READ");// eta sub-range //_file0 = new TFile("Results/PDC_13c_0p4to0p8eta.root","READ");// eta sub-range //_file0 = new TFile("Results/PDC_13c_neg0p4toneg0p8eta.root","READ");// eta sub-range //_file0 = new TFile("Results/PDC_13bc_dEta0p025dPhi0p055_kINT7_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_13bc_FB5and7overlap_V0binning.root","READ"); //_file0 = new TFile("Results/PDC_13bc_kAnyINT_V0binning.root","READ"); //_file0 = new TFile("Results/PDC_13bcd_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_13bcd_kINT7.root","READ"); //_file0 = new TFile("Results/PDC_13f_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_13e_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_13bcd_kINT7_plus_kHighMult.root","READ");// v5 and before //_file0 = new TFile("Results/PDC_13bcde_kINT7_NclsFix.root","READ"); //_file0 = new TFile("Results/PDC_13e_kHighMult_NclsFix.root","READ"); _file0 = new TFile("Results/PDC_13bcde_kINT7_kHighMult_NclsFix.root","READ");// standard }else{ _file0 = new TFile("Results/PDC_13b2_efix_p1234_R2_2Ktbins.root","READ");// standard (gen level) //_file0 = new TFile("Results/PDC_13b2_efix_p1234_R2_RecLevel.root","READ");// Reconstructed tracks //_file0 = new TFile("Results/PDC_13b2_p1_noTTC.root","READ"); } }else{// pp if(!MCcase){ //_file0 = new TFile("Results/PDC_10cde.root","READ"); //_file0 = new TFile("Results/PDC_10cde_3Ktbins.root","READ"); //_file0 = new TFile("Results/PDC_10e_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_10e_kAnyINT.root","READ"); //_file0 = new TFile("Results/PDC_10cde_kAnyINT.root","READ"); //_file0 = new TFile("Results/PDC_10cde_kAnyINT_Plus_kHighMult.root","READ"); //_file0 = new TFile("Results/K0sRange_10cde.root","READ");// for K0s peak removal //_file0 = new TFile("Results/PDC_10cde_2Kt3bins.root","READ"); //_file0 = new TFile("Results/PDC_10e_kHighMult.root","READ"); //_file0 = new TFile("Results/PDC_10cde_kMB.root","READ"); //_file0 = new TFile("Results/PDC_10cde_kMB_plus_kHighMult_2Kt3bins.root","READ");// v2 paper //_file0 = new TFile("Results/PDC_10cde_kMB_plus_kHighMult_AOD147.root","READ");// v5 and before //_file0 = new TFile("Results/PDC_10cde_FB5and7overlap.root","READ"); //_file0 = new TFile("Results/PDC_10d_noTTC.root","READ"); //_file0 = new TFile("Results/PDC_10cde_kMB.root","READ"); _file0 = new TFile("Results/PDC_10cde_kMB_kHighMult_NclsFix.root","READ");// standard }else{ _file0 = new TFile("Results/PDC_10f6a_R2_2Ktbins.root","READ");// standard //_file0 = new TFile("Results/PDC_10f6a_R2_AOD161.root","READ"); //_file0 = new TFile("Results/PDC_10f6a_R10_genLevel_2Ktbins.root","READ"); } } ReadCoulCorrections(FSIindex); ReadMomResFile(Mbin); ReadMuonCorrections(CollisionType,Mbin); // ///////////////////////////////////////////////////////// double NormQcutLow; double NormQcutHigh; if(CollisionType==0) {// PbPb if(Mbin<10){ NormQcutLow = 0.15;// was 0.15 NormQcutHigh = 0.175;// was 0.175 }else{ NormQcutLow = 0.3;// was 0.3 NormQcutHigh = 0.35;// was 0.35 } }else{// pPb and pp NormQcutLow = 1.0;// was 1.0 NormQcutHigh = 1.2;// was 1.2 } TList *MyList; TDirectoryFile *tdir; if(!MCcase) tdir = (TDirectoryFile*)_file0->Get("PWGCF.outputThreePionRadiiAnalysis.root"); if(CollisionType==0){ if(!MCcase){ if(Mbin<6) MyList=(TList*)tdir->Get("ThreePionRadiiOutput_1"); else MyList=(TList*)tdir->Get("ThreePionRadiiOutput_2"); }else{ MyList=(TList*)_file0->Get("MyList"); } }else { if(!MCcase) { if(CollisionType==2 && Mbin<15) MyList=(TList*)tdir->Get("ThreePionRadiiOutput_2"); else MyList=(TList*)tdir->Get("ThreePionRadiiOutput_1"); }else MyList=(TList*)_file0->Get("MyList"); } _file0->Close(); TH1D *Events = (TH1D*)MyList->FindObject("fEvents2"); // cout<<"#Events = "<<Events->Integral(Mbin+1,Mbin+1)<<endl; // below lines determine fractional cross-sections without correction for efficiency //TH1F *fMultDist3 = (TH1F*)MyList->FindObject("fMultDist3");// total event count before track cuts //for(int mm=19; mm>=10; mm--){ //cout<<"fracton of cross-section = "<<Events->GetBinContent(mm+1)/fMultDist3->GetEntries()<<endl; //} TF1 *Unity = new TF1("Unity","1",0,10); Unity->SetLineStyle(2); Unity->SetLineColor(1); // Explicit Loop Histos TH2D *Two_ex_2d[2][2][Sbins_2][2]; TH1D *Two_ex[2][2][Sbins_2][2]; // Normalizations double NormH_pc[5]={0}; double norm_ex_2[6][2]={{0}}; // 3d method histograms TH3D *Three_3d[2][2][2][Sbins_3][5]; TH1D *Three_1d[2][2][2][Sbins_3][5]; /////////////////////////////////// // Create Histograms // 2-particle for(int c1=0; c1<2; c1++){ for(int c2=0; c2<2; c2++){ for(int sc=0; sc<Sbins_2; sc++){ for(int term=0; term<2; term++){ TString *name2_ex = new TString("Explicit2_Charge1_"); *name2_ex += c1; name2_ex->Append("_Charge2_"); *name2_ex += c2; name2_ex->Append("_SC_"); *name2_ex += sc; name2_ex->Append("_M_"); *name2_ex += Mbin; name2_ex->Append("_ED_"); *name2_ex += 0;// EDbin name2_ex->Append("_Term_"); *name2_ex += term+1; if(sc==0 || sc==3 || sc==5){ if( (c1+c2)==1 ) {if(c1!=0) continue;}// skip degenerate histogram } Two_ex_2d[c1][c2][sc][term] = (TH2D*)MyList->FindObject(name2_ex->Data()); Two_ex_2d[c1][c2][sc][term]->Sumw2(); TString *proname = new TString(name2_ex->Data()); proname->Append("_pro"); if(Ktbin==10) {Ktlowbin=1; Kthighbin=Two_ex_2d[c1][c2][sc][term]->GetNbinsX();}//full kt interval Two_ex[c1][c2][sc][term] = (TH1D*)Two_ex_2d[c1][c2][sc][term]->ProjectionY(proname->Data(),Ktlowbin,Kthighbin);// bins:5-6 (kt:0.2-0.3) norm_ex_2[sc][term] = Two_ex[c1][c2][sc][term]->Integral(Two_ex[c1][c2][sc][term]->GetXaxis()->FindBin(NormQcutLow),Two_ex[c1][c2][sc][term]->GetXaxis()->FindBin(NormQcutHigh)); Two_ex[c1][c2][sc][term]->Scale(norm_ex_2[sc][0]/norm_ex_2[sc][term]); Two_ex[c1][c2][sc][term]->SetMarkerStyle(20); Two_ex[c1][c2][sc][term]->GetXaxis()->SetTitle("#font[12]{q}_{inv}"); Two_ex[c1][c2][sc][term]->GetYaxis()->SetTitle("#font[12]{C}_{2}"); Two_ex[c1][c2][sc][term]->SetTitle(""); }// term }// sc // 3-particle for(int c3=0; c3<2; c3++){ for(int sc=0; sc<Sbins_3; sc++){ for(int term=0; term<5; term++){ TString *name3_ex = new TString("Explicit3_Charge1_"); *name3_ex += c1; name3_ex->Append("_Charge2_"); *name3_ex += c2; name3_ex->Append("_Charge3_"); *name3_ex += c3; name3_ex->Append("_SC_"); *name3_ex += sc; name3_ex->Append("_M_"); *name3_ex += Mbin; name3_ex->Append("_ED_"); *name3_ex += EDbin; name3_ex->Append("_Term_"); *name3_ex += term+1; TString *name3_pc = new TString("PairCut3_Charge1_"); *name3_pc += c1; name3_pc->Append("_Charge2_"); *name3_pc += c2; name3_pc->Append("_Charge3_"); *name3_pc += c3; name3_pc->Append("_SC_"); *name3_pc += sc; name3_pc->Append("_M_"); *name3_pc += Mbin; name3_pc->Append("_ED_"); *name3_pc += EDbin; name3_pc->Append("_Term_"); *name3_pc += term+1; /////////////////////////////////////// // skip degenerate histograms if(sc==0 || sc==6 || sc==9){// Identical species if( (c1+c2+c3)==1) {if(c1!=0 || c2!=0 || c3!=1) continue;} if( (c1+c2+c3)==2) {if(c1!=0) continue;} }else if(sc!=5){ if( (c1+c2)==1) {if(c1!=0) continue;} }else {}// do nothing for pi-k-p case ///////////////////////////////////////// if(PairCut){ TString *nameNorm=new TString("PairCut3_Charge1_"); *nameNorm += c1; nameNorm->Append("_Charge2_"); *nameNorm += c2; nameNorm->Append("_Charge3_"); *nameNorm += c3; nameNorm->Append("_SC_"); *nameNorm += sc; nameNorm->Append("_M_"); *nameNorm += Mbin; nameNorm->Append("_ED_"); *nameNorm += 0;// EDbin nameNorm->Append("_Term_"); *nameNorm += term+1; nameNorm->Append("_Norm"); NormH_pc[term] = ((TH1D*)MyList->FindObject(nameNorm->Data()))->Integral(); if(NormH_pc[term] > 0){ if(sc<=2) { TString *name_3d = new TString(name3_pc->Data()); name_3d->Append("_3D"); Three_3d[c1][c2][c3][sc][term] = (TH3D*)MyList->FindObject(name_3d->Data()); Three_3d[c1][c2][c3][sc][term]->Sumw2(); Three_3d[c1][c2][c3][sc][term]->Scale(NormH_pc[0]/NormH_pc[term]); TString *name_Q3 = new TString(name3_pc->Data()); name_Q3->Append("_Q3"); Three_1d[c1][c2][c3][sc][term] = (TH1D*)MyList->FindObject(name_Q3->Data()); Three_1d[c1][c2][c3][sc][term]->Sumw2(); Three_1d[c1][c2][c3][sc][term]->Scale(NormH_pc[0]/NormH_pc[term]); //cout<<"Scale factor = "<<NormH_pc[0]/NormH_pc[term]<<endl; } }else{ cout<<"normalization = 0. Skipping this SC. SC="<<sc<<" c1="<<c1<<" c2="<<c2<<" c3="<<c3<<endl; } }// pair cut }// term }// SC }// c3 }// c2 }// c1 cout<<"Done getting Histograms"<<endl; TCanvas *can1 = new TCanvas("can1", "can1",11,53,800,800); can1->SetHighLightColor(2); //can1->Range(-0.7838115,-1.033258,0.7838115,1.033258); gStyle->SetOptFit(0111); can1->SetFillColor(10); can1->SetBorderMode(0); can1->SetBorderSize(2); can1->SetFrameFillColor(0); can1->SetFrameBorderMode(0); can1->SetFrameBorderMode(0); gPad->SetRightMargin(0.025); gPad->SetLeftMargin(0.1); gPad->SetTopMargin(0.02); gPad->SetGridx(0); gPad->SetGridy(0); TLegend *legend1 = new TLegend(.4,.67,.87,.87,NULL,"brNDC"); legend1->SetBorderSize(1); legend1->SetTextSize(.04);// small .03; large .036 legend1->SetFillColor(0); TLatex *Specif_System = new TLatex(0.23, 0.9,"ALICE p-Pb #sqrt{#font[12]{s}_{NN}}=5.02 TeV, #LT#font[12]{N}_{ch}#GT = 9.8 #pm 0.5"); TLatex *Specif_KT3 = new TLatex(0.23, 0.8,"0.16<#font[12]{K}_{T,3}<0.3 GeV/#font[12]{c}"); TLatex *Specif_kT = new TLatex(0.23, 0.8,"0.2<#font[12]{k}_{T}<0.3 GeV/#font[12]{c}"); TLatex *Specif_FSI = new TLatex(0.23, 0.7,"FSI uncorrected"); TLatex *Specif_Disclaimer = new TLatex(0.23, 0.6,"Statistical errors only"); Specif_FSI->SetTextFont(42); Specif_System->SetTextFont(42); Specif_KT3->SetTextFont(42); Specif_Disclaimer->SetTextFont(42); Specif_kT->SetTextFont(42); Specif_FSI->SetTextSize(0.04); Specif_System->SetTextSize(0.04); Specif_KT3->SetTextSize(0.04); Specif_Disclaimer->SetTextSize(0.04); Specif_kT->SetTextSize(0.04); Specif_FSI->SetNDC(); Specif_System->SetNDC(); Specif_KT3->SetNDC(); Specif_Disclaimer->SetNDC(); Specif_kT->SetNDC(); ///////////////////////// // 2-particle legend // 0 = pi-pi // 1 = pi-k // 2 = pi-p // 3 = k-k // 4 = k-p // 5 = p-p ///////////////////////// TH1D *Two_pipi_mp = (TH1D*)Two_ex[0][1][0][0]->Clone(); Two_pipi_mp->Divide(Two_ex[0][1][0][1]); // Normalization range counting. Just to evaluate required normalization width in qinv. //cout<<Two_ex[0][0][0][0]->Integral(Two_ex[0][0][0][0]->GetXaxis()->FindBin(0), Two_ex[0][0][0][0]->GetXaxis()->FindBin(0.1))<<endl; //cout<<Two_ex[0][0][0][0]->Integral(Two_ex[0][0][0][0]->GetXaxis()->FindBin(1.06), Two_ex[0][0][0][0]->GetXaxis()->FindBin(1.1))<<endl; //cout<<Two_ex[0][0][0][0]->Integral(Two_ex[0][0][0][0]->GetXaxis()->FindBin(0.15), Two_ex[0][0][0][0]->GetXaxis()->FindBin(0.175))<<endl; const int SCOI_2=0; TH1D *Two_ex_clone_mm=(TH1D*)Two_ex[0][0][SCOI_2][0]->Clone(); Two_ex_clone_mm->Divide(Two_ex[0][0][SCOI_2][1]); TH1D *Two_ex_clone_mp=(TH1D*)Two_ex[0][1][SCOI_2][0]->Clone(); Two_ex_clone_mp->Divide(Two_ex[0][1][SCOI_2][1]); TH1D *Two_ex_clone_pp=(TH1D*)Two_ex[1][1][SCOI_2][0]->Clone(); Two_ex_clone_pp->Divide(Two_ex[1][1][SCOI_2][1]); Two_ex_clone_mm->GetYaxis()->SetTitle("#font[12]{C}_{2}"); Two_ex_clone_mm->SetTitle(""); Two_ex_clone_mp->GetYaxis()->SetTitle("#font[12]{C}_{2}"); Two_ex_clone_mm->SetMarkerColor(2); Two_ex_clone_mm->SetLineColor(2); Two_ex_clone_mp->SetMarkerColor(1); Two_ex_clone_pp->SetMarkerColor(4); Two_ex_clone_pp->SetLineColor(4); Two_ex_clone_mm->GetXaxis()->SetRangeUser(0,0.6); Two_ex_clone_mm->SetMinimum(0.95); Two_ex_clone_mm->SetMaximum(2.4); if(MCcase){ Two_ex_clone_mp->SetMarkerColor(4); Two_ex_clone_mp->SetLineColor(4); Two_ex_clone_mm->Add(Two_ex_clone_pp); Two_ex_clone_mm->Scale(0.5); Two_ex_clone_mm->GetYaxis()->SetTitle("#font[12]{C}_{2}^{MC}"); Two_ex_clone_mm->GetYaxis()->SetTitleOffset(1.3); Two_ex_clone_mm->SetMinimum(0.95); Two_ex_clone_mm->SetMaximum(1.3); Two_ex_clone_mm->DrawCopy(); gStyle->SetOptFit(1111); Two_ex_clone_mm->Fit(BkgMCFit,"","",0,1.8); //BkgMCFit->Draw("same"); Two_ex_clone_mp->DrawCopy("same"); legend1->AddEntry(Two_ex_clone_mm,"same-charge","p"); //legend1->AddEntry(Two_ex_clone_pp,"++","p"); legend1->AddEntry(Two_ex_clone_mp,"mixed-charge","p"); legend1->Draw("same"); } //for(int i=0; i<100; i++){ //cout<<Two_ex_clone_mm->GetBinContent(i+1)<<", "; //} TF1 *C2osFit_CMS = new TF1("C2osFit_CMS","[0] + [1]*exp(-pow([2]*x,2))",0,2); C2osFit_CMS->FixParameter(0, 9.93558e-01); C2osFit_CMS->FixParameter(1, 5.19578e-02); C2osFit_CMS->FixParameter(2, 1.89425e+00); ///////////////////////////////////////////////////// // Global fitting C2os and C2ss double C2ssSys_e[BINRANGE_C2global]={0}; double C2osSys_e[BINRANGE_C2global]={0}; // const int npar=10;// 10 TMinuit MyMinuit(npar); MyMinuit.SetFCN(fcnC2_Global); double arglist_C2 = 0; int dummy=0; MyMinuit.mnexcm("SET NOWarnings",&arglist_C2,1, dummy); arglist_C2 = -1; MyMinuit.mnexcm("SET PRint",&arglist_C2,1, dummy); double OutputPar[npar]={0}; double OutputPar_e[npar]={0}; double par[npar]; // the start values double stepSize[npar]; // step sizes double minVal[npar]; // minimum bound on parameter double maxVal[npar]; // maximum bound on parameter string parName[npar]; TH1D *temp_mm=(TH1D*)Two_ex[0][0][SCOI_2][0]->Clone(); temp_mm->Divide(Two_ex[0][0][SCOI_2][1]); TH1D *temp_mp=(TH1D*)Two_ex[0][1][SCOI_2][0]->Clone(); temp_mp->Divide(Two_ex[0][1][SCOI_2][1]); TH1D *temp_pp=(TH1D*)Two_ex[1][1][SCOI_2][0]->Clone(); temp_pp->Divide(Two_ex[1][1][SCOI_2][1]); TH1D *C2ssRaw=(TH1D*)temp_mm->Clone();// MRC uncorrected TH1D *C2osRaw=(TH1D*)temp_mp->Clone();// MRC uncorrected C2ssRaw->SetMarkerStyle(24); C2osRaw->SetMarkerStyle(21);//21 C2ssRaw->SetMarkerColor(2); C2osRaw->SetMarkerColor(4); for(int i=0; i<BINRANGE_C2global; i++){ if(i >= temp_mm->GetNbinsX()) continue;// bin limit C2ssFitting[i] = (temp_mm->GetBinContent(i+1) + temp_pp->GetBinContent(i+1))/2.; double MRC_SC = 1.0; double MRC_MC = 1.0; if(MomResC2[0]->GetBinContent(i+1) > 0) MRC_SC = (MomResC2[0]->GetBinContent(i+1)-1)*MRCShift + 1; if(MomResC2[1]->GetBinContent(i+1) > 0) MRC_MC = (MomResC2[1]->GetBinContent(i+1)-1)*MRCShift + 1; if(!GeneratedSignal && !MCcase) C2ssFitting[i] *= MRC_SC; C2ssFitting_e[i] = pow(MRC_SC*sqrt(pow(temp_mm->GetBinError(i+1),2) + pow(temp_pp->GetBinError(i+1),2))/sqrt(2.),2); C2ssRaw->SetBinContent(i+1, (temp_mm->GetBinContent(i+1) + temp_pp->GetBinContent(i+1))/2.); C2ssRaw->SetBinError(i+1, pow(sqrt(pow(temp_mm->GetBinError(i+1),2) + pow(temp_pp->GetBinError(i+1),2))/sqrt(2.),2)); C2ssFitting_e[i] = sqrt(C2ssFitting_e[i]); C2osFitting[i] = temp_mp->GetBinContent(i+1); if(!GeneratedSignal && !MCcase) C2osFitting[i] *= MRC_MC; C2osFitting_e[i] = pow(MRC_MC*temp_mp->GetBinError(i+1),2); C2osRaw->SetBinContent(i+1, temp_mp->GetBinContent(i+1)); C2osRaw->SetBinError(i+1, temp_mp->GetBinError(i+1)); C2osFitting_e[i] += pow((MRC_MC-1)*0.1 * temp_mp->GetBinContent(i+1),2); C2osFitting_e[i] = sqrt(C2osFitting_e[i]); // C2ssSys_e[i] = pow((MRC_SC-1)*0.1 * (temp_mm->GetBinContent(i+1) + temp_pp->GetBinContent(i+1))/2.,2); C2ssSys_e[i] = sqrt(C2ssSys_e[i]); C2osSys_e[i] = pow((MRC_MC-1)*0.1 * temp_mp->GetBinContent(i+1),2); C2osSys_e[i] = sqrt(C2osSys_e[i]); // if(UseC2Bkg) C2ssFitting[i] /= BkgMCFit->Eval(BinCenters[i]); // // method with undilution for plotting purposes //C2ssFitting[i] = (C2ssFitting[i] - (1-lambda_PM)) / (CoulCorr2(+1, BinCenters[i]) * lambda_PM); // divide CMS background //C2ssFitting[i] /= C2osFit_CMS->Eval(BinCenters[i]); // K2SS[i] = CoulCorr2(+1, BinCenters[i]); K2OS[i] = CoulCorr2(-1, BinCenters[i]); //K2SS[i] = 1; //K2OS[i] = 1; // K2SS_e[i] = 0.0; K2OS_e[i] = 0.0; } C2ssFitting[0]=-1000; C2osFitting[0]=-1000; C2ssFitting_e[0]=1000; C2osFitting_e[0]=1000; K2SS_e[0]=1000; K2OS_e[0]=1000; int ierflg=0; double arglist[10]; //arglist[0]=2;// improve Minimization Strategy (1 is default) //MyMinuit.mnexcm("SET STR",arglist,1,ierflg); //arglist[0] = 0; //MyMinuit.mnexcm("SCAN", arglist,1,ierflg); arglist[0] = 5000;// 5000 MyMinuit.mnexcm("MIGRAD", arglist ,1,ierflg); TF1 *fitC2ss_Base = new TF1("fitC2ss_Base",C2ssFitFunction, 0.005,2, npar);//0.2 TF1 *fitC2ss_Expan = new TF1("fitC2ss_Expan",C2ssFitFunction, 0.005,2, npar);//0.2 fitC2ss_Base->SetLineStyle(2); TH1D *fitC2ss_h = new TH1D("fitC2ss_h","",Two_ex_clone_mm->GetNbinsX(),Two_ex_clone_mm->GetXaxis()->GetBinLowEdge(1), Two_ex_clone_mm->GetXaxis()->GetBinUpEdge(Two_ex_clone_mm->GetNbinsX())); for(int ft=0; ft<2; ft++){// Gaussian or EW if(ft==1 && !IncludeExpansion) continue; par[0] = 1.0; par[1]=0.7; par[2]=0.5; par[3]=7.2; par[4] = .1; par[5] = .0; par[6] = .0; par[7] = 0.; par[8] = 0.; par[9] = 0.;// was par[1]=0.5 stepSize[0] = 0.01; stepSize[1] = 0.01; stepSize[2] = 0.02; stepSize[3] = 0.2; stepSize[4] = 0.01; stepSize[5] = 0.001; stepSize[6] = 0.001; stepSize[7] = 0.001; stepSize[8]=0.001; stepSize[9]=0.01; minVal[0] = 0.955; minVal[1] = 0.5; minVal[2] = 0.; minVal[3] = 0.1; minVal[4] = 0.001; minVal[5] = -10.; minVal[6] = -10.; minVal[7] = -10.; minVal[8]=-10; minVal[9] = 0.995;// minVal[1] was 0.2 maxVal[0] = 1.1; maxVal[1] = 1.3; maxVal[2] = 0.99; maxVal[3] = 15.; maxVal[4] = 2.; maxVal[5] = 10.; maxVal[6] = 10.; maxVal[7] = 10.; maxVal[8]=10.; maxVal[9]=1.1;// maxVal[1] was 1.0 if(Gaussian==kFALSE) {maxVal[1]=2.0; maxVal[3]=20.0;} parName[0] = "Norm"; parName[1] = "#Lambda"; parName[2] = "G"; parName[3] = "Rch"; parName[4] = "Rcoh"; parName[5] = "coeff_3"; parName[6] = "coeff_4"; parName[7] = "coeff_5"; parName[8]="coeff_6"; parName[9]="Norm_2"; for (int i=0; i<npar; i++){ MyMinuit.DefineParameter(i, parName[i].c_str(), par[i], stepSize[i], minVal[i], maxVal[i]); } //MyMinuit.DefineParameter(1, parName[1].c_str(), 0.7, stepSize[1], minVal[1], maxVal[1]); MyMinuit.FixParameter(1);// lambda //MyMinuit.DefineParameter(0, parName[0].c_str(), .998, stepSize[0], minVal[0], maxVal[0]); MyMinuit.FixParameter(0);// N MyMinuit.DefineParameter(2, parName[2].c_str(), 0., stepSize[2], minVal[2], maxVal[2]); MyMinuit.FixParameter(2);// G MyMinuit.DefineParameter(4, parName[4].c_str(), 0., stepSize[4], minVal[4], maxVal[4]); MyMinuit.FixParameter(4);// Rcoh MyMinuit.DefineParameter(7, parName[7].c_str(), 0, stepSize[7], minVal[7], maxVal[7]); MyMinuit.FixParameter(7);// k5 MyMinuit.DefineParameter(8, parName[8].c_str(), 0, stepSize[8], minVal[8], maxVal[8]); MyMinuit.FixParameter(8);// k6 // if(ft==0){ MyMinuit.DefineParameter(5, parName[5].c_str(), 0, stepSize[5], minVal[5], maxVal[5]); MyMinuit.FixParameter(5);// k3 MyMinuit.DefineParameter(6, parName[6].c_str(), 0, stepSize[6], minVal[6], maxVal[6]); MyMinuit.FixParameter(6);// k4 }else{// IncludeExpansion if(FixExpansionAvg){ if(Gaussian){ MyMinuit.DefineParameter(5, parName[5].c_str(), kappa3, stepSize[5], minVal[5], maxVal[5]); MyMinuit.FixParameter(5);// k3 MyMinuit.DefineParameter(6, parName[6].c_str(), kappa4, stepSize[6], minVal[6], maxVal[6]); MyMinuit.FixParameter(6);// k4 }else{ MyMinuit.DefineParameter(5, parName[5].c_str(), 0, stepSize[5], minVal[5], maxVal[5]); MyMinuit.FixParameter(5);// k3 MyMinuit.DefineParameter(6, parName[6].c_str(), 0, stepSize[6], minVal[6], maxVal[6]); MyMinuit.FixParameter(6);// k4 } }else{ Int_t err=0; Double_t tmp[1]; tmp[0] = 5+1; MyMinuit.mnexcm( "RELease", tmp, 1, err );// kappa3 tmp[0] = 6+1; MyMinuit.mnexcm( "RELease", tmp, 1, err );// kappa4 // } } if(IncludeSS==kFALSE){ MyMinuit.DefineParameter(3, parName[3].c_str(), 9.1, stepSize[3], minVal[3], maxVal[3]); MyMinuit.FixParameter(3);// Rch MyMinuit.DefineParameter(0, parName[0].c_str(), .998, stepSize[0], minVal[0], maxVal[0]); MyMinuit.FixParameter(0);// N MyMinuit.DefineParameter(5, parName[5].c_str(), 0, stepSize[5], minVal[5], maxVal[5]); MyMinuit.FixParameter(5);// k3 MyMinuit.DefineParameter(6, parName[6].c_str(), 0, stepSize[6], minVal[6], maxVal[6]); MyMinuit.FixParameter(6);// k4 } if(IncludeOS==kFALSE){ MyMinuit.DefineParameter(9, parName[9].c_str(), 1.002, stepSize[9], minVal[9], maxVal[9]); MyMinuit.FixParameter(9);// N_2 } // Do the minimization! if(!MCcase){ cout<<"Start C2 Global fit"<<endl; MyMinuit.Migrad();// generally the best minimization algorithm cout<<"End C2 Global fit"<<endl; } for (int i=0; i<npar; i++){ MyMinuit.GetParameter(i,OutputPar[i],OutputPar_e[i]); } MyMinuit.mnexcm("SHOw PARameters", &arglist_C2, 1, ierflg); cout<<"C2 fit: Chi2/NDF = "<<Chi2_C2global/(NFitPoints_C2global - MyMinuit.GetNumFreePars())<<" Chi^2="<<Chi2_C2global<<endl; if(ft==0){// Gaussian or Exponential for(int i=0; i<npar; i++) { fitC2ss_Base->FixParameter(i,OutputPar[i]); fitC2ss_Base->SetParError(i,OutputPar_e[i]); } for(int bin=1; bin<=Two_ex_clone_mm->GetNbinsX(); bin++){ double qinv = Two_ex_clone_mm->GetXaxis()->GetBinCenter(bin); if(!MCcase) fitC2ss_h->SetBinContent(bin, fitC2ss_Base->Eval(qinv));// Base fit } }else{// Edgeworth or Laguerre for(int i=0; i<npar; i++) { fitC2ss_Expan->FixParameter(i,OutputPar[i]); fitC2ss_Expan->SetParError(i,OutputPar_e[i]); } for(int bin=1; bin<=Two_ex_clone_mm->GetNbinsX(); bin++){ //double qinv = Two_ex_clone_mm->GetXaxis()->GetBinCenter(bin); //if(!MCcase) fitC2ss_h->SetBinContent(bin, fitC2ss_Expan->Eval(qinv));// uncomment to show expansion } double lambdaStar=0, lambdaStar_e=0; if(Gaussian){ lambdaStar = OutputPar[1]*pow(1 + OutputPar[6]/8.,2); lambdaStar_e = pow(OutputPar_e[1]*(1 + OutputPar[6]/8.),2); lambdaStar_e = sqrt(lambdaStar_e); }else{ lambdaStar = OutputPar[1]*pow(1 - OutputPar[5] + OutputPar[6],2); lambdaStar_e = pow(OutputPar_e[1]*(1 - OutputPar[5] + OutputPar[6]),2); lambdaStar_e = sqrt(lambdaStar_e); } cout<<"lambda* = "<<lambdaStar<<" +- "<<lambdaStar_e<<endl; } //if(ft==0) {MainFitParams[0]=OutputPar[3]; MainFitParams[2]=OutputPar[1];} //else {MainFitParams[4]=OutputPar[3]; MainFitParams[6]=OutputPar[1];} }// ft TH1D *C2_ss=(TH1D*)Two_ex_clone_mm->Clone(); TH1D *C2_os=(TH1D*)Two_ex_clone_mp->Clone(); TH1D *C2_ss_Momsys=(TH1D*)Two_ex_clone_mm->Clone(); TH1D *C2_os_Momsys=(TH1D*)Two_ex_clone_mp->Clone(); TH1D *C2_ss_Ksys=(TH1D*)Two_ex_clone_mm->Clone(); TH1D *C2_os_Ksys=(TH1D*)Two_ex_clone_mp->Clone(); TH1D *K2_ss = (TH1D*)Two_ex_clone_mm->Clone(); TH1D *K2_os = (TH1D*)Two_ex_clone_mp->Clone(); for(int i=0; i<BINRANGE_C2global; i++){ C2_ss->SetBinContent(i+1, C2ssFitting[i]); C2_os->SetBinContent(i+1, C2osFitting[i]); C2_ss_Momsys->SetBinContent(i+1, C2ssFitting[i]); C2_os_Momsys->SetBinContent(i+1, C2osFitting[i]); C2_ss_Ksys->SetBinContent(i+1, C2ssFitting[i]); C2_os_Ksys->SetBinContent(i+1, C2osFitting[i]); double Toterror_ss = sqrt(fabs(pow(C2ssFitting_e[i],2)-pow(C2ssSys_e[i],2))); double Toterror_os = sqrt(fabs(pow(C2osFitting_e[i],2)-pow(C2osSys_e[i],2))); C2_ss->SetBinError(i+1, Toterror_ss); C2_os->SetBinError(i+1, Toterror_os); C2_ss_Momsys->SetBinError(i+1, C2ssSys_e[i]); C2_os_Momsys->SetBinError(i+1, C2osSys_e[i]); C2_ss_Ksys->SetBinError(i+1, OutputPar[1]*K2SS_e[i]); C2_os_Ksys->SetBinError(i+1, OutputPar[1]*K2OS_e[i]); // K2_ss->SetBinContent(i+1, K2SS[i]); K2_ss->SetBinError(i+1, 0); K2_os->SetBinContent(i+1, K2OS[i]); K2_os->SetBinError(i+1, 0); } C2_ss_Momsys->SetMarkerSize(0); C2_ss_Momsys->SetFillStyle(1000); C2_ss_Momsys->SetFillColor(kRed-10); C2_os_Momsys->SetMarkerSize(0); C2_os_Momsys->SetFillStyle(1000); C2_os_Momsys->SetFillColor(kBlue-10); C2_ss_Ksys->SetMarkerSize(0); C2_ss_Ksys->SetFillStyle(3004); C2_ss_Ksys->SetFillColor(kRed); C2_os_Ksys->SetMarkerSize(0); C2_os_Ksys->SetFillStyle(3004); C2_os_Ksys->SetFillColor(kBlue); C2_ss->GetXaxis()->SetRangeUser(0,Q2Limit);// 0,0.098 C2_ss->GetYaxis()->SetRangeUser(0.95,2.4);//0.98,1.3 C2_ss->GetXaxis()->SetTitleOffset(.8); C2_ss->GetYaxis()->SetTitleOffset(.85); C2_ss->GetXaxis()->SetTitle("#font[12]{q} (GeV/c)"); C2_ss->GetYaxis()->SetTitle("#font[12]{C}_{2}"); C2_ss->GetXaxis()->SetTitleSize(0.05); C2_ss->GetYaxis()->SetTitleSize(0.05); C2_os->GetXaxis()->SetRangeUser(0,0.6); C2_os->GetYaxis()->SetRangeUser(0.98,1.3); C2_os->GetXaxis()->SetTitleOffset(.8); C2_os->GetYaxis()->SetTitleOffset(.8); C2_os->GetXaxis()->SetTitle("#font[12]{q} (GeV/c)"); C2_os->GetXaxis()->SetTitleSize(0.05); C2_os->GetYaxis()->SetTitleSize(0.05); //C2_ss->SetMarkerSize(1.5); //C2_os->SetMarkerSize(1.5); C2_os->SetMarkerStyle(25); C2_os->SetMarkerColor(4); fitC2ss_h->SetLineWidth(2); fitC2ss_h->SetLineColor(2); TH1D *temp_hist_ss=new TH1D("temp_hist_ss","",100,0,0.5); TH1D *temp_hist_os=new TH1D("temp_hist_os","",100,0,0.5); TH1D *temp_fit=new TH1D("temp_fit","",100,0,0.5); // PbPb, M16 double temp_points_ss[100]={-1000, 1.31587, 1.40841, 1.4426, 1.36602, 1.37548, 1.34059, 1.30318, 1.27143, 1.28487, 1.26105, 1.21905, 1.2025, 1.19743, 1.15297, 1.16947, 1.14259, 1.13562, 1.13135, 1.10088, 1.10227, 1.08037, 1.0824, 1.08069, 1.0657, 1.06808, 1.05341, 1.0493, 1.05081, 1.04546, 1.04644, 1.03703, 1.03975, 1.04479, 1.03218, 1.02657, 1.02588, 1.02756, 1.02345, 1.02306, 1.02024, 1.01551, 1.01133, 1.01025, 1.00697, 1.01819, 1.00651, 1.01293, 1.00529, 0.997971, 1.00141, 1.00665, 1.01479, 1.01668, 1.00657, 1.01188, 1.0042, 0.995877, 1.01094, 1.0036, 1.00312, 0.994031, 1.0084, 1.00604, 1.00329, 1.00677, 0.997553, 0.992874, 0.993799, 0.995368, 0.999215, 1.00508, 0.993053, 0.996744, 0.990297, 0.990715, 1.00114, 0.996376, 0.984527, 0.99474, 0.993457, 1.00723, 0.999042, 0.996589, 0.988385, 0.9916, 0.995184, 0.988197, 1.00019, 0.998052, 0.985506, 0.99024, 0.996102, 0.995797, 0.991887, 0.99276, 0.992561, 0.991261, 0.996243}; double temp_points_ss_e[100]={1000, 0.206703, 0.0920405, 0.0637848, 0.0465513, 0.0380257, 0.0318382, 0.0270753, 0.0236593, 0.0216127, 0.0195971, 0.0176654, 0.0163346, 0.0153927, 0.0141902, 0.0136715, 0.0128887, 0.0124019, 0.0119795, 0.0113642, 0.0111615, 0.0106718, 0.0104775, 0.0102842, 0.0100296, 0.00990822, 0.00962737, 0.0094588, 0.00935572, 0.00923423, 0.0091266, 0.00897368, 0.00889453, 0.00883246, 0.00864815, 0.00853381, 0.00846384, 0.00839173, 0.00829825, 0.00822044, 0.00817065, 0.00805352, 0.00798793, 0.00791266, 0.00784862, 0.00788769, 0.00775376, 0.00772075, 0.00766719, 0.00757909, 0.00754696, 0.00755056, 0.00756756, 0.00754031, 0.00746571, 0.00745114, 0.00739226, 0.00729679, 0.00737659, 0.00730319, 0.00730321, 0.0072228, 0.00728131, 0.00726075, 0.0072364, 0.00723868, 0.00716992, 0.00712198, 0.00713157, 0.00713724, 0.00714391, 0.007175, 0.00708838, 0.00711573, 0.00709965, 0.00707786, 0.00713124, 0.00712164, 0.00703573, 0.00712383, 0.00709407, 0.00720098, 0.00713966, 0.00716466, 0.00710145, 0.00711087, 0.00714815, 0.00710751, 0.00717636, 0.00718936, 0.00712875, 0.00715855, 0.00718908, 0.00719846, 0.00719133, 0.00720147, 0.00721384, 0.00720763, 0.0072536}; for(int bin=1; bin<100; bin++){ //cout<<C2_ss->GetBinContent(bin)<<", "; //cout<<C2_ss->GetBinError(bin)<<", "; //cout<<fitC2ss_h->GetBinContent(bin)<<", "; temp_hist_ss->SetBinContent(bin, temp_points_ss[bin-1]); temp_hist_ss->SetBinError(bin, temp_points_ss_e[bin-1]); //temp_hist_os->SetBinContent(bin, temp_points_os[bin-1]); //temp_hist_os->SetBinError(bin, temp_points_os_e[bin-1]); //temp_fit->SetBinContent(bin, temp_points_fit[bin-1]); } //cout<<endl; temp_hist_ss->SetMarkerStyle(20); temp_hist_ss->SetMarkerColor(1); temp_hist_os->SetMarkerStyle(24); temp_hist_os->SetMarkerColor(1); temp_fit->SetLineColor(1); C2_os->SetMarkerStyle(24); C2_os->SetMarkerColor(2); if(!MCcase){ //C2_ss->SetMaximum(1.6); C2_ss->DrawCopy(); //legend1->AddEntry(C2_ss,"same-charge","p"); C2_os->DrawCopy("same"); //legend1->AddEntry(C2_os,"mixed-charge","p"); /* Specif_System->Draw("same"); Specif_kT->Draw("same"); Specif_FSI->Draw("same"); Specif_Disclaimer->Draw("same"); TLatex *Stats1 = new TLatex(.5,.55, "#lambda = 0.41 #pm 0.004"); Stats1->SetNDC(); Stats1->SetTextSize(0.06); //Stats1->Draw("same"); TLatex *Stats2 = new TLatex(.5,.45, "R_{inv} = 1.59 #pm 0.01 fm"); Stats2->SetNDC(); Stats2->SetTextSize(0.06); */ //Stats2->Draw("same"); /*TF1 *BkgFitC2 = new TF1("BkgFitC2","1+[0]*exp(-pow([1]*x/0.19733,2))",0,1); BkgFitC2->SetParameter(0,0.08); BkgFitC2->SetParameter(1,0.5); BkgFitC2->SetLineColor(1); C2_ss->Fit(BkgFitC2,"IME","",0.2,0.8); BkgFitC2->Draw("same");*/ //fitC2ss_h->SetLineColor(1); Unity->Draw("same"); fitC2ss_h->Draw("C same"); //temp_hist_ss->Draw("same"); //temp_hist_os->Draw("same"); //temp_fit->Draw("C same"); //fitC2ss_Base->Draw("C same"); //fitC2ss_Expan->Draw("C same"); //legend1->AddEntry(C2_ss,"p-Pb (++)","p"); //legend1->AddEntry(C2_os,"p-Pb (-+)","p"); //legend1->AddEntry(temp_hist_ss,"Pb-Pb (++), Mbin 16","p"); //legend1->AddEntry(temp_hist_os,"Pb-Pb (-+), <Nch>=36","p"); legend1->Draw("same"); } cout<<"============================================"<<endl; cout<<"Start 3-pion section"<<endl; TCanvas *can2 = new TCanvas("can2", "can2",800,0,800,800);//800,0,600,900 can2->SetHighLightColor(2); gStyle->SetOptFit(0111); can2->SetFillColor(10); can2->SetBorderMode(0); can2->SetBorderSize(2); can2->SetGridx(); can2->SetGridy(); can2->SetFrameFillColor(0); can2->SetFrameBorderMode(0); can2->SetFrameBorderMode(0); can2->cd(); gPad->SetRightMargin(0.02); TLegend *legend2 = new TLegend(.58,.55, .97,.65,NULL,"brNDC"); legend2->SetBorderSize(1); legend2->SetTextSize(.03);// small .03; large .06 legend2->SetFillColor(0); ///////////////////////////////////////////// TH3D *C3QS_3d = new TH3D("C3QS_3d","",BINRANGE_3,0,.4, BINRANGE_3,0,.4, BINRANGE_3,0,.4); TH3D *Combinatorics_3d = new TH3D("Combinatorics_3d","",BINRANGE_3,0,.4, BINRANGE_3,0,.4, BINRANGE_3,0,.4); // const float Q3HistoLimit=0.5; const int Q3BINS = int((Q3HistoLimit+0.0001)/0.01); TH1D *num_withRS = new TH1D("num_withRS","",Q3BINS,0,Q3HistoLimit); TH1D *num_withGRS = new TH1D("num_withGRS","",Q3BINS,0,Q3HistoLimit); TH1D *num_withTM = new TH1D("num_withTM","",Q3BINS,0,Q3HistoLimit); TH1D *Cterm1 = new TH1D("Cterm1","",Q3BINS,0,Q3HistoLimit); TH1D *Cterm1_noMRC = new TH1D("Cterm1_noMRC","",Q3BINS,0,Q3HistoLimit); TH1D *Cterm2 = new TH1D("Cterm2","",Q3BINS,0,Q3HistoLimit); TH1D *Cterm3 = new TH1D("Cterm3","",Q3BINS,0,Q3HistoLimit); TH1D *Cterm4 = new TH1D("Cterm4","",Q3BINS,0,Q3HistoLimit); TH1D *num_QS = new TH1D("num_QS","",Q3BINS,0,Q3HistoLimit); TH1D *Combinatorics_1d = new TH1D("Combinatorics_1d","",Q3BINS,0,Q3HistoLimit); TH1D *Combinatorics_1d_noMRC = new TH1D("Combinatorics_1d_noMRC","",Q3BINS,0,Q3HistoLimit); TH1D *Coul_Riverside = new TH1D("Coul_Riverside","",Q3BINS,0,Q3HistoLimit); TH1D *Coul_GRiverside = new TH1D("Coul_GRiverside","",Q3BINS,0,Q3HistoLimit); TH1D *c3_hist = new TH1D("c3_hist","",Q3BINS,0,Q3HistoLimit); if(SameCharge) {Cterm1_noMRC->Sumw2(); Combinatorics_1d_noMRC->Sumw2();} // double num_QS_e[Q3BINS]; double c3_e[Q3BINS]; for(int ii=0; ii<Q3BINS; ii++){ num_QS_e[ii]=0.; c3_e[ii]=0.; } // CB=Charge Bin. 0 means pi-, 1 means pi+ int CB1=0, CB2=0, CB3=0; int CP12=1, CP13=1, CP23=1; int CB1_2=0, CB2_2=0, CB3_2=0; if(CHARGE==-1) { if(SameCharge) {CB1=0; CB2=0; CB3=0;} else {CB1=0; CB2=0; CB3=1; CP12=+1, CP13=-1, CP23=-1;} }else { if(SameCharge) {CB1=1; CB2=1; CB3=1;} else {CB1=0; CB2=1; CB3=1; CP12=-1, CP13=-1, CP23=+1;} } if(AddCC){ if(CHARGE==-1) { if(SameCharge) {CB1_2=1; CB2_2=1; CB3_2=1;} else {CB1_2=0; CB2_2=1; CB3_2=1;}; }else { if(SameCharge) {CB1_2=0; CB2_2=0; CB3_2=0;} else {CB1_2=0; CB2_2=0; CB3_2=1;} } } // SC = species combination. Always 0 meaning pi-pi-pi. int SCBin=0; // TH1D *GenSignalExpected_num=new TH1D("GenSignalExpected_num","",Q3BINS,0,Q3HistoLimit); TH1D *GenSignalExpected_den=new TH1D("GenSignalExpected_den","",Q3BINS,0,Q3HistoLimit); // double value_num; double value_num_e; double N3_QS; double N3_QS_e; for(int i=2; i<=BINLIMIT_3; i++){// bin number double Q12 = BinCenters[i-1];// true center //int Q12bin = int(Q12/0.01)+1; // for(int j=2; j<=BINLIMIT_3; j++){// bin number double Q13 = BinCenters[j-1];// true center //int Q13bin = int(Q13/0.01)+1; // for(int k=2; k<=BINLIMIT_3; k++){// bin number double Q23 = BinCenters[k-1];// true center //int Q23bin = int(Q23/0.01)+1; // double Q3 = sqrt(pow(Q12,2) + pow(Q13,2) + pow(Q23,2)); int Q3bin = Cterm1->GetXaxis()->FindBin(Q3); if(Q12 < sqrt(pow(Q13,2)+pow(Q23,2) - 2*Q13*Q23)) continue;// not all configurations are possible if(Q12 > sqrt(pow(Q13,2)+pow(Q23,2) + 2*Q13*Q23)) continue;// not all configurations are possible // double G_12 = Gamov(CP12, Q12);// point-source Coulomb correlation double G_13 = Gamov(CP13, Q13); double G_23 = Gamov(CP23, Q23); double K2_12 = CoulCorr2(CP12, Q12);// finite-source Coulomb+Strong correlation from Therminator. double K2_13 = CoulCorr2(CP13, Q13); double K2_23 = CoulCorr2(CP23, Q23); double K3 = 1.;// 3-body Coulomb+Strong correlation if(SameCharge || CHARGE==-1) K3 = CoulCorrGRS(SameCharge, Q12, Q13, Q23); else K3 = CoulCorrGRS(SameCharge, Q23, Q12, Q13); if(MCcase && !GeneratedSignal) { K2_12=1.; K2_13=1.; K2_23=1.; K3=1.;} if(K3==0) continue; double TERM1=Three_3d[CB1][CB2][CB3][SCBin][0]->GetBinContent(i,j,k);// N3 (3-pion yield per q12,q13,q23 cell). 3-pions from same-event double TERM2=Three_3d[CB1][CB2][CB3][SCBin][1]->GetBinContent(i,j,k);// N2*N1. 1 and 2 from same-event double TERM3=Three_3d[CB1][CB2][CB3][SCBin][2]->GetBinContent(i,j,k);// N2*N1. 1 and 3 from same-event double TERM4=Three_3d[CB1][CB2][CB3][SCBin][3]->GetBinContent(i,j,k);// N2*N1. 2 and 3 from same-event double TERM5=Three_3d[CB1][CB2][CB3][SCBin][4]->GetBinContent(i,j,k);// N1*N1*N1. All from different events (pure combinatorics) if(AddCC){ if(CHARGE==-1){// base is (SC,MC,MC) TERM1 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][0]->GetBinContent(j,k,i); TERM2 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][1]->GetBinContent(j,k,i); TERM3 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][2]->GetBinContent(j,k,i); TERM4 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][3]->GetBinContent(j,k,i); TERM5 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][4]->GetBinContent(j,k,i); }else {// base is (MC,MC,SC) TERM1 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][0]->GetBinContent(k,i,j); TERM2 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][1]->GetBinContent(k,i,j); TERM3 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][2]->GetBinContent(k,i,j); TERM4 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][3]->GetBinContent(k,i,j); TERM5 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][4]->GetBinContent(k,i,j); } } if(TERM1==0 && TERM2==0 && TERM3==0 && TERM4==0 && TERM5==0) continue; if(TERM1==0 || TERM2==0 || TERM3==0 || TERM4==0 || TERM5==0) continue; double muonCorr_C3=1.0, muonCorr_C2_12=1.0, muonCorr_C2_13=1.0, muonCorr_C2_23=1.0; if(MuonCorrection && !MCcase){ if(SameCharge){ muonCorr_C3 = C3muonCorrection->GetBinContent(C3muonCorrection->GetXaxis()->FindBin(Q3)); muonCorr_C2_12 = C2muonCorrection_SC->GetBinContent(C2muonCorrection_SC->GetXaxis()->FindBin(Q12)); muonCorr_C2_13 = C2muonCorrection_SC->GetBinContent(C2muonCorrection_SC->GetXaxis()->FindBin(Q13)); muonCorr_C2_23 = C2muonCorrection_SC->GetBinContent(C2muonCorrection_SC->GetXaxis()->FindBin(Q23)); }else{ muonCorr_C3 = C3muonCorrection->GetBinContent(C3muonCorrection->GetXaxis()->FindBin(Q3)); if(CHARGE==-1){// pi-pi-pi+ muonCorr_C2_12 = C2muonCorrection_SC->GetBinContent(C2muonCorrection_SC->GetXaxis()->FindBin(Q12)); muonCorr_C2_13 = C2muonCorrection_MC->GetBinContent(C2muonCorrection_MC->GetXaxis()->FindBin(Q13)); muonCorr_C2_23 = C2muonCorrection_MC->GetBinContent(C2muonCorrection_MC->GetXaxis()->FindBin(Q23)); }else{// pi-pi+pi+ muonCorr_C2_12 = C2muonCorrection_MC->GetBinContent(C2muonCorrection_MC->GetXaxis()->FindBin(Q12)); muonCorr_C2_13 = C2muonCorrection_MC->GetBinContent(C2muonCorrection_MC->GetXaxis()->FindBin(Q13)); muonCorr_C2_23 = C2muonCorrection_SC->GetBinContent(C2muonCorrection_SC->GetXaxis()->FindBin(Q23)); } } } // apply momentum resolution correction if(!MCcase && !GeneratedSignal){ int MRC_Q3bin = int(Q3/0.01) + 1; if(SameCharge){ TERM1 *= (MomRes1d[0][0]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM2 *= (MomRes1d[0][1]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM3 *= (MomRes1d[0][2]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM4 *= (MomRes1d[0][3]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM5 *= (MomRes1d[0][4]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; }else{ // removed due to over-correction of 1-D MRC approximation TERM1 *= (MomRes1d[1][0]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM2 *= (MomRes1d[1][1]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM3 *= (MomRes1d[1][2]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM4 *= (MomRes1d[1][3]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; TERM5 *= (MomRes1d[1][4]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; } } // TwoFrac=lambda_PM;// 0.7 OneFrac=pow(TwoFrac,.5); ThreeFrac=pow(TwoFrac,1.5); // Purify. Isolate pure 3-pion QS correlations using Lambda and K3 (removes lower order correlations) N3_QS = TERM1; N3_QS -= ( pow(1-OneFrac,3) + 3*OneFrac*pow(1-OneFrac,2) )*TERM5; N3_QS -= (1-OneFrac)*(TERM2 + TERM3 + TERM4 - 3*(1-TwoFrac)*TERM5); N3_QS /= ThreeFrac; N3_QS /= K3; N3_QS *= muonCorr_C3; num_QS->Fill(Q3, N3_QS); // Isolate 3-pion cumulant value_num = N3_QS; value_num -= (TERM2 - (1-TwoFrac)*TERM5)/K2_12/TwoFrac * muonCorr_C2_12; value_num -= (TERM3 - (1-TwoFrac)*TERM5)/K2_13/TwoFrac * muonCorr_C2_13; value_num -= (TERM4 - (1-TwoFrac)*TERM5)/K2_23/TwoFrac * muonCorr_C2_23; value_num += 2*TERM5; // errors N3_QS_e = TERM1; N3_QS_e += pow(( pow(1-OneFrac,3) +3*OneFrac*pow(1-OneFrac,2) )*sqrt(TERM5),2); N3_QS_e += (pow((1-OneFrac),2)*(TERM2 + TERM3 + TERM4) + pow((1-OneFrac)*3*(1-TwoFrac)*sqrt(TERM5),2)); N3_QS_e /= pow(ThreeFrac,2); N3_QS_e /= pow(K3,2); num_QS_e[Q3bin-1] += N3_QS_e; // errors value_num_e = N3_QS_e; value_num_e += (pow(1/K2_12/TwoFrac*sqrt(TERM2),2) + pow((1-TwoFrac)/K2_12/TwoFrac*sqrt(TERM5),2)); value_num_e += (pow(1/K2_13/TwoFrac*sqrt(TERM3),2) + pow((1-TwoFrac)/K2_13/TwoFrac*sqrt(TERM5),2)); value_num_e += (pow(1/K2_23/TwoFrac*sqrt(TERM4),2) + pow((1-TwoFrac)/K2_23/TwoFrac*sqrt(TERM5),2)); value_num_e += pow(2*sqrt(TERM5),2); c3_e[Q3bin-1] += value_num_e + TERM5;// add baseline stat error // Fill histograms c3_hist->Fill(Q3, value_num + TERM5);// for cumulant correlation function C3QS_3d->SetBinContent(i,j,k, N3_QS); C3QS_3d->SetBinError(i,j,k, N3_QS_e); // Coul_Riverside->Fill(Q3, G_12*G_13*G_23*TERM5); Coul_GRiverside->Fill(Q3, TERM5*CoulCorrGRS(SameCharge, Q12, Q13, Q23)); num_withGRS->Fill(Q3, N3_QS); Cterm1->Fill(Q3, TERM1); Cterm2->Fill(Q3, TERM2); Cterm3->Fill(Q3, TERM3); Cterm4->Fill(Q3, TERM4); Combinatorics_1d->Fill(Q3, TERM5); Combinatorics_3d->Fill(Q12,Q13,Q23, TERM5); double Q3_m = sqrt(pow((i-0.5)*(0.005),2) + pow((j-0.5)*(0.005),2) + pow((k-0.5)*(0.005),2)); Cterm1_noMRC->Fill(Q3_m, Three_3d[CB1][CB2][CB3][SCBin][0]->GetBinContent(i,j,k)); Combinatorics_1d_noMRC->Fill(Q3_m, Three_3d[CB1][CB2][CB3][SCBin][4]->GetBinContent(i,j,k)); // A_3[i-1][j-1][k-1] = value_num + TERM5; B_3[i-1][j-1][k-1] = TERM5; A_3_e[i-1][j-1][k-1] = sqrt(value_num_e + TERM5); //if(i==j && i==k && i==4) cout<<A_3[i-1][j-1][k-1]<<" "<<B_3[i-1][j-1][k-1]<<" "<<A_3_e[i-1][j-1][k-1]<<endl; /////////////////////////////////////////////////////////// } } } //////////////////////////// // Intermediate steps num_withRS->Sumw2(); num_withGRS->Sumw2(); num_withTM->Sumw2(); Cterm1->Sumw2(); Cterm2->Sumw2(); Cterm3->Sumw2(); Cterm4->Sumw2(); Combinatorics_1d->Sumw2(); Combinatorics_3d->Sumw2(); for(int i=0; i<Q3BINS; i++) {c3_hist->SetBinError(i+1, sqrt(c3_e[i])); num_QS->SetBinError(i+1, sqrt(num_QS_e[i]));} num_withRS->Divide(Combinatorics_1d); num_withGRS->Divide(Combinatorics_1d); num_withTM->Divide(Combinatorics_1d); Cterm1->Divide(Combinatorics_1d); Cterm1_noMRC->Divide(Combinatorics_1d_noMRC); Cterm2->Divide(Combinatorics_1d); Cterm3->Divide(Combinatorics_1d); Cterm4->Divide(Combinatorics_1d); c3_hist->Divide(Combinatorics_1d); num_QS->Divide(Combinatorics_1d); /////////////////////////////////////////////////////////////////////////////////////////////////// Cterm1->SetMarkerStyle(20); Cterm1->SetMarkerColor(4); Cterm1->SetLineColor(4); Cterm1->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/#font[12]{c})"); Cterm1->GetYaxis()->SetTitle("#font[12]{C}_{3}"); //Cterm1->SetTitle("#pi^{-}#pi^{-}#pi^{-}"); Cterm1->SetMinimum(0.97); Cterm1->SetMaximum(5.7);// 6.1 for same-charge Cterm1->GetXaxis()->SetRangeUser(0,Q3Limit); Cterm1->GetYaxis()->SetTitleOffset(1.4); Cterm1->Draw(); legend2->AddEntry(Cterm1,"#font[12]{C}_{3} raw","p"); // Cterm2->SetMarkerStyle(20); Cterm2->SetMarkerColor(7); Cterm3->SetMarkerStyle(25); Cterm3->SetMarkerColor(8); Cterm4->SetMarkerStyle(26); Cterm4->SetMarkerColor(9); //Cterm2->Draw("same"); //Cterm3->Draw("same"); //Cterm4->Draw("same"); //legend2->AddEntry(Cterm1,"++-","p"); if(MCcase){ double C3points_HIJING_mmm[10]={0, 0.961834, 1.01827, 0.999387, 1.00202, 1.00081, 1.00082, 1.00037, 0.999074, 0.999099}; double C3points_HIJING_mmm_e[10]={0, 0.0833895, 0.015158, 0.0047978, 0.00235067, 0.00138155, 0.00087485, 0.000612203, 0.000450162, 0.000346943}; double C3points_HIJING_ppp[10]={0, 1.13015, 1.00623, 1.00536, 1.00293, 0.999964, 1.00116, 1.0007, 1.00037, 1.00105}; double C3points_HIJING_ppp_e[10]={0, 0.0977058, 0.0150694, 0.0048196, 0.00235518, 0.00138376, 0.000877562, 0.000614069, 0.000452051, 0.00034856}; TH1D *C3_HIJING_mmm=(TH1D*)Cterm1->Clone(); TH1D *C3_HIJING_ppp=(TH1D*)Cterm1->Clone(); for(int i=0; i<10; i++){ C3_HIJING_mmm->SetBinContent(i+1, C3points_HIJING_mmm[i]); C3_HIJING_mmm->SetBinError(i+1, C3points_HIJING_mmm_e[i]); C3_HIJING_ppp->SetBinContent(i+1, C3points_HIJING_ppp[i]); C3_HIJING_ppp->SetBinError(i+1, C3points_HIJING_ppp_e[i]); } C3_HIJING_mmm->SetMarkerColor(2); C3_HIJING_mmm->SetLineColor(2); C3_HIJING_ppp->SetMarkerColor(4); C3_HIJING_ppp->SetLineColor(4); //C3_HIJING_mmm->Draw("same"); //C3_HIJING_ppp->Draw("same"); //legend2->AddEntry(C3_HIJING_mmm,"---","p"); //legend2->AddEntry(C3_HIJING_ppp,"+++","p"); } num_QS->SetMarkerStyle(24); num_QS->SetMarkerColor(4); num_QS->SetLineColor(4); num_QS->GetXaxis()->SetTitle("#font[12]{Q}_{3}"); num_QS->GetYaxis()->SetTitle("#font[12]{C}_{3}^{QS}"); num_QS->GetXaxis()->SetRangeUser(0,.12); num_QS->SetMaximum(6); //num_QS->Draw("same"); //legend2->AddEntry(num_QS,"C_{3}^{QS}","p"); c3_hist->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/#font[12]{c})"); c3_hist->GetYaxis()->SetTitle("#font[12]{C}_{3} or #font[12]{#bf{c}}_{3}"); c3_hist->GetYaxis()->SetTitleOffset(1.4); c3_hist->GetXaxis()->SetRangeUser(0,Q3Limit); c3_hist->GetYaxis()->SetRangeUser(0.9,4); c3_hist->SetMarkerStyle(25); c3_hist->SetMarkerColor(2); c3_hist->SetLineColor(2); c3_hist->SetMaximum(3); c3_hist->SetMinimum(.7); if(!MCcase) c3_hist->Draw("same"); legend2->AddEntry(c3_hist,"#font[12]{#bf{c}}_{3} (cumulant correlation)","p"); c3_hist->Draw(); //for(int i=1; i<=15; i++) cout<<c3_hist->GetBinContent(i)<<", "; //cout<<endl; //for(int i=1; i<=15; i++) cout<<c3_hist->GetBinError(i)<<", "; //cout<<endl; const int npar_c3=5;// 10 TMinuit MyMinuit_c3(npar_c3); MyMinuit_c3.SetFCN(fcn_c3); int ierflg_c3=0; double arglist_c3 = 0; MyMinuit_c3.mnexcm("SET NOWarnings",&arglist_c3,1, ierflg_c3); arglist_c3 = -1; MyMinuit_c3.mnexcm("SET PRint",&arglist_c3,1, ierflg_c3); //arglist_c3=2;// improve Minimization Strategy (1 is default) //MyMinuit_c3.mnexcm("SET STR",&arglist_c3,1,ierflg_c3); //arglist_c3 = 0; //MyMinuit_c3.mnexcm("SCAN", &arglist_c3,1,ierflg_c3); arglist_c3 = 5000; MyMinuit_c3.mnexcm("MIGRAD", &arglist_c3 ,1,ierflg_c3); TF1 *c3Fit1D_Base=new TF1("c3Fit1D_Base","[0]*(1+[1]*exp(-pow([2]*x/0.19733,[3]) * [4]/2.))",0,1); c3Fit1D_Base->FixParameter(3,ExpPower); TF1 *c3Fit1D_Expan; if(Gaussian) { c3Fit1D_Base->FixParameter(4,1.0); c3Fit1D_Expan=new TF1("c3Fit1D_Expan","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.) * pow(1 + ([3]/(6.*pow(2,1.5))*(8.*pow([2]*x/sqrt(3.)/0.19733,3) - 12.*pow([2]*x/sqrt(3.)/0.19733,1))) + ([4]/(24.*pow(2,2))*(16.*pow([2]*x/sqrt(3.)/0.19733,4) -48.*pow([2]*x/sqrt(3.)/0.19733,2) + 12)),3))",0,1); }else { c3Fit1D_Base->FixParameter(4,sqrt(3.)); c3Fit1D_Expan=new TF1("c3Fit1D_Expan","[0]*(1+[1]*exp(-pow([2]*x/0.19733,1) * sqrt(3.)/2.) * pow(1 + [3]*([2]*x/sqrt(3.)/0.19733 - 1) + [4]/2*(pow([2]*x/sqrt(3.)/0.19733,2) - 4*[2]*x/sqrt(3.)/0.19733 + 2),3))",0,1); } double OutputPar_c3[npar_c3]={0}; double OutputPar_c3_e[npar_c3]={0}; double par_c3[npar_c3]; // the start values double stepSize_c3[npar_c3]; // step sizes double minVal_c3[npar_c3]; // minimum bound on parameter double maxVal_c3[npar_c3]; // maximum bound on parameter string parName_c3[npar_c3]; if(SameCharge && !MCcase){ for(int ft=0; ft<2; ft++){// Gaussian, Edgeworth if(ft==1 && !IncludeExpansion) continue; par_c3[0] = 1.0; par_c3[1] = .5; par_c3[2] = 5; par_c3[3] = kappa3; par_c3[4] = kappa4;// was par_c3[3] = 0.; par_c3[4] = 0.; stepSize_c3[0] = 0.01; stepSize_c3[1] = 0.1; stepSize_c3[2] = 0.2; stepSize_c3[3] = 0.01; stepSize_c3[4] = 0.01; minVal_c3[0] = 0.95; minVal_c3[1] = 0.2; minVal_c3[2] = 0.5; minVal_c3[3] = -1; minVal_c3[4] = -1;// was minVal_c3[3] = -1; minVal_c3[4] = -1; maxVal_c3[0] = 1.1; maxVal_c3[1] = 2.5; maxVal_c3[2] = 15.; maxVal_c3[3] = +1; maxVal_c3[4] = +1;// was minVal_c3[3] = +1; minVal_c3[4] = +1; parName_c3[0] = "N"; parName_c3[1] = "#lambda_{3}"; parName_c3[2] = "R_{inv}"; parName_c3[3] = "coeff_{3}"; parName_c3[4] = "coeff_{4}"; c3Fit1D_Base->SetParName(0,"N"); c3Fit1D_Base->SetParName(1,"#lambda_{3}"); c3Fit1D_Base->SetParName(2,"R_{inv}"); if(!Gaussian) {maxVal_c3[1] = 5.0; maxVal_c3[2] = 20.;} for (int i=0; i<npar_c3; i++){ MyMinuit_c3.DefineParameter(i, parName_c3[i].c_str(), par_c3[i], stepSize_c3[i], minVal_c3[i], maxVal_c3[i]); } // if(!FixExpansionAvg) { maxVal_c3[1] = 2.0; /*double RadiusFix; if(CollisionType==0) RadiusFix = 10 - 7*Mbin/15.; else RadiusFix = 2.5 - 1.25*(Mbin-12)/7.; MyMinuit_c3.DefineParameter(2, parName_c3[2].c_str(), RadiusFix, stepSize_c3[2], minVal_c3[2], maxVal_c3[2]); MyMinuit_c3.FixParameter(2); MyMinuit_c3.DefineParameter(1, parName_c3[1].c_str(), 2.0, stepSize_c3[1], minVal_c3[1], maxVal_c3[1]); MyMinuit_c3.FixParameter(1); */ } // kappa_3 and kappa_4 if(ft==0){// Gaussian MyMinuit_c3.DefineParameter(3, parName_c3[3].c_str(), 0., stepSize_c3[3], minVal_c3[3], maxVal_c3[3]); MyMinuit_c3.FixParameter(3); MyMinuit_c3.DefineParameter(4, parName_c3[4].c_str(), 0., stepSize_c3[4], minVal_c3[4], maxVal_c3[4]); MyMinuit_c3.FixParameter(4); }else{// Edgeworth if(FixExpansionAvg){ if(Gaussian){ MyMinuit_c3.DefineParameter(3, parName_c3[3].c_str(), kappa3, stepSize_c3[3], minVal_c3[3], maxVal_c3[3]); MyMinuit_c3.FixParameter(3); MyMinuit_c3.DefineParameter(4, parName_c3[4].c_str(), kappa4, stepSize_c3[4], minVal_c3[4], maxVal_c3[4]); MyMinuit_c3.FixParameter(4); }else{ MyMinuit_c3.DefineParameter(3, parName_c3[3].c_str(), 0, stepSize_c3[3], minVal_c3[3], maxVal_c3[3]); MyMinuit_c3.FixParameter(3); MyMinuit_c3.DefineParameter(4, parName_c3[4].c_str(), 0, stepSize_c3[4], minVal_c3[4], maxVal_c3[4]); MyMinuit_c3.FixParameter(4); } }else{ Int_t err=0; Double_t tmp[1]; tmp[0] = 3+1; MyMinuit_c3.mnexcm( "RELease", tmp, 1, err );// kappa3 tmp[0] = 4+1; MyMinuit_c3.mnexcm( "RELease", tmp, 1, err );// kappa4 } } ///////////////////////////////////////////////////////////// // Do the minimization! cout<<"Start Three-d fit"<<endl; MyMinuit_c3.Migrad();// Minuit's best minimization algorithm cout<<"End Three-d fit"<<endl; ///////////////////////////////////////////////////////////// MyMinuit_c3.mnexcm("SHOw PARameters", &arglist_c3, 1, ierflg); cout<<"c3 Fit: Chi2 = "<<Chi2_c3<<" NDF = "<<(NFitPoints_c3-MyMinuit_c3.GetNumFreePars())<<endl; for(int i=0; i<npar_c3; i++){ MyMinuit_c3.GetParameter(i,OutputPar_c3[i],OutputPar_c3_e[i]); } if(ft==0){ c3Fit1D_Base->SetParameter(0,OutputPar_c3[0]); c3Fit1D_Base->SetParError(0,OutputPar_c3_e[0]); c3Fit1D_Base->SetParameter(1,OutputPar_c3[1]); c3Fit1D_Base->SetParError(1,OutputPar_c3_e[1]); c3Fit1D_Base->SetParameter(2,OutputPar_c3[2]); c3Fit1D_Base->SetParError(2,OutputPar_c3_e[2]); c3Fit1D_Base->SetLineStyle(2); c3Fit1D_Base->Draw("same"); }else{ c3Fit1D_Expan->SetParameter(0,OutputPar_c3[0]); c3Fit1D_Expan->SetParError(0,OutputPar_c3_e[0]); c3Fit1D_Expan->SetParameter(1,OutputPar_c3[1]); c3Fit1D_Expan->SetParError(1,OutputPar_c3_e[1]); c3Fit1D_Expan->SetParameter(2,OutputPar_c3[2]); c3Fit1D_Expan->SetParError(2,OutputPar_c3_e[2]); c3Fit1D_Expan->SetParameter(3,OutputPar_c3[3]); c3Fit1D_Expan->SetParError(3,OutputPar_c3_e[3]); c3Fit1D_Expan->SetParameter(4,OutputPar_c3[4]); c3Fit1D_Expan->SetParError(4,OutputPar_c3_e[4]); c3Fit1D_Expan->SetLineStyle(1); //c3Fit1D_Expan->Draw("same"); double lambdaStar=0, lambdaStar_e=0; if(Gaussian){ lambdaStar = OutputPar_c3[1]*pow(1 + OutputPar_c3[4]/8.,3); lambdaStar_e = pow(OutputPar_c3_e[1]*pow(1 + OutputPar_c3[4]/8.,3),2); lambdaStar_e = sqrt(lambdaStar_e); }else{ lambdaStar = OutputPar_c3[1]*pow(1 - OutputPar_c3[3] + OutputPar_c3[4],3); lambdaStar_e = pow(OutputPar_c3_e[1]*pow(1 - OutputPar_c3[3] + OutputPar_c3[4],3),2); lambdaStar_e = sqrt(lambdaStar_e); } cout<<"lambda*_3 = "<<lambdaStar<<" +- "<<lambdaStar_e<<endl; } //if(ft==0) {MainFitParams[1]=OutputPar_c3[2]; MainFitParams[3]=OutputPar_c3[1];} //else {MainFitParams[5]=OutputPar_c3[2]; MainFitParams[7]=OutputPar_c3[1];} }// fit type }// SC and !MC // project 3D EW fit onto 1D histogram if(SameCharge && !MCcase){ for(int i=2; i<=BINLIMIT_3; i++){// bin number double Q12 = BinCenters[i-1];// true center for(int j=2; j<=BINLIMIT_3; j++){// bin number double Q13 = BinCenters[j-1];// true center for(int k=2; k<=BINLIMIT_3; k++){// bin number double Q23 = BinCenters[k-1];// true center // double Q3 = sqrt(pow(Q12,2) + pow(Q13,2) + pow(Q23,2)); if(Q12 < sqrt(pow(Q13,2)+pow(Q23,2) - 2*Q13*Q23)) continue;// not all configurations are possible if(Q12 > sqrt(pow(Q13,2)+pow(Q23,2) + 2*Q13*Q23)) continue;// not all configurations are possible double TERM5=Three_3d[CB1][CB2][CB3][SCBin][4]->GetBinContent(i,j,k);// N1*N1*N1. All from different events (pure combinatorics) if(TERM5==0) continue; // if(AddCC){ if(CHARGE==-1){// base is (SC,MC,MC) TERM5 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][4]->GetBinContent(j,k,i); }else {// base is (MC,MC,SC) TERM5 += Three_3d[CB1_2][CB2_2][CB3_2][SCBin][4]->GetBinContent(k,i,j); } } // int MRC_Q3bin = int(Q3/0.01) + 1; TERM5 *= (MomRes1d[0][4]->GetBinContent(MRC_Q3bin)-1)*MRCShift+1; // double radius_exp = c3Fit1D_Expan->GetParameter(2)/FmToGeV; double arg12 = Q12*radius_exp; double arg13 = Q13*radius_exp; double arg23 = Q23*radius_exp; Float_t EW12 = 1 + c3Fit1D_Expan->GetParameter(3)/(6.*pow(2.,1.5))*(8.*pow(arg12,3) - 12.*arg12); EW12 += c3Fit1D_Expan->GetParameter(4)/(24.*pow(2.,2))*(16.*pow(arg12,4) -48.*pow(arg12,2) + 12); Float_t EW13 = 1 + c3Fit1D_Expan->GetParameter(3)/(6.*pow(2.,1.5))*(8.*pow(arg13,3) - 12.*arg13); EW13 += c3Fit1D_Expan->GetParameter(4)/(24.*pow(2.,2))*(16.*pow(arg13,4) -48.*pow(arg13,2) + 12); Float_t EW23 = 1 + c3Fit1D_Expan->GetParameter(3)/(6.*pow(2.,1.5))*(8.*pow(arg23,3) - 12.*arg23); EW23 += c3Fit1D_Expan->GetParameter(4)/(24.*pow(2.,2))*(16.*pow(arg23,4) -48.*pow(arg23,2) + 12); // double cumulant_fitvalue=0; cumulant_fitvalue = c3Fit1D_Expan->GetParameter(0)*(c3Fit1D_Expan->GetParameter(1)*EW12*EW13*EW23*exp(-pow(radius_exp,ExpPower)/2. * (pow(Q12,ExpPower)+pow(Q13,ExpPower)+pow(Q23,ExpPower)))*TERM5 + TERM5); GenSignalExpected_num->Fill(Q3, cumulant_fitvalue); GenSignalExpected_den->Fill(Q3, TERM5); /////////////////////////////////////////////////////////// } } } } GenSignalExpected_num->Sumw2(); GenSignalExpected_num->Divide(GenSignalExpected_den); TSpline3 *c3Fit1D_ExpanSpline = new TSpline3(GenSignalExpected_num); if(CollisionType==0) c3Fit1D_ExpanSpline->SetLineColor(1); if(CollisionType==1) c3Fit1D_ExpanSpline->SetLineColor(2); if(CollisionType==2) c3Fit1D_ExpanSpline->SetLineColor(4); c3Fit1D_ExpanSpline->SetLineWidth(2); double xpoints[1000], ypoints[1000]; bool splineOnly=kFALSE; for(int iii=0; iii<1000; iii++){ xpoints[iii] = 0 + (iii+0.5)*0.001; if(!Gaussian && CollisionType!=0 && c3Fit1D_Expan->Eval(xpoints[iii]) < 2.2) splineOnly=kTRUE; if(!Gaussian && CollisionType==0 && c3Fit1D_Expan->Eval(xpoints[iii]) < 2.0) splineOnly=kTRUE;// 2.0 or 1.4 for Mbin=3 in Fig 2 if(!splineOnly){ ypoints[iii] = c3Fit1D_Expan->Eval(xpoints[iii]); if(c3Fit1D_Expan->Eval(xpoints[iii])<2 && fabs(c3Fit1D_Expan->Eval(xpoints[iii])-c3Fit1D_ExpanSpline->Eval(xpoints[iii])) < 0.01) splineOnly=kTRUE; } else {ypoints[iii] = c3Fit1D_ExpanSpline->Eval(xpoints[iii]); splineOnly=kTRUE;} } TGraph *gr_c3Spline = new TGraph(1000,xpoints,ypoints); gr_c3Spline->SetLineWidth(2); if(CollisionType==0) gr_c3Spline->SetLineColor(1); if(CollisionType==1) gr_c3Spline->SetLineColor(2); if(CollisionType==2) gr_c3Spline->SetLineColor(4); gr_c3Spline->Draw("c same"); double ChiSqSum_1D=0, SumNDF_1D=0; for(int bin=1; bin<=300; bin++){ double binCenter = c3_hist->GetXaxis()->GetBinCenter(bin); if(binCenter > Q3Limit) continue; if(c3_hist->GetBinError(bin)==0) continue; int grIndex=1; for(int gr=0; gr<999; gr++){ if(binCenter > xpoints[gr] && (binCenter < xpoints[gr+1])) {grIndex=gr; break;} } //ChiSqSum_1D += pow((c3_hist->GetBinContent(bin)-ypoints[grIndex]) / c3_hist->GetBinError(bin),2); ChiSqSum_1D += pow((c3_hist->GetBinContent(bin)-c3Fit1D_Base->Eval(binCenter)) / c3_hist->GetBinError(bin),2); //cout<<c3_hist->GetBinContent(bin)<<" "<<c3Fit1D_Base->Eval(binCenter)<<" "<<c3_hist->GetBinError(bin)<<endl; //cout<<c3_hist->GetBinContent(bin)<<" "<<ypoints[grIndex]<<" "<<c3_hist->GetBinError(bin)<<endl; SumNDF_1D++; } cout<<"1D Chi2/NDF = "<<ChiSqSum_1D / (SumNDF_1D-3.)<<endl; // uncomment to display fit box //c3_hist->GetListOfFunctions()->Add(c3Fit1D_Base); // Fit Range Systematics //for(int i=0; i<8; i++){cout<<MainFitParams[i]<<", ";} //cout<<endl; //for(int i=0; i<4; i++){cout<<int(100*(MainFitParams[i]-RefMainFitParams[i])/MainFitParams[i])<<"$\\%$ & ";}// Gaussian //cout<<endl; //for(int i=4; i<8; i++){cout<<int(100*(MainFitParams[i]-RefMainFitParams[i])/MainFitParams[i])<<"$\\%$ & ";}// Edgeworth //cout<<endl; //TString *lamName=new TString("#lambda_{3}="); //TString *RName=new TString("#R_{inv,3}="); //*lamName += int(100*c3Fit1D_Base->GetParameter(1))/100.; lamName->Append(" #pm "); *lamName += int(100*c3Fit1D_Base->GetParError(1))/100.; //*RName += round(100*c3Fit1D_Base->GetParameter(2))/100; RName->Append(" #pm "); *RName += int(100*c3Fit1D_Base->GetParError(2))/100.; RName->Append(" fm"); //TLatex *fitBox1=new TLatex(0.5,0.9,lamName->Data()); //fitBox1->SetNDC(kTRUE); //fitBox1->Draw(); //TLatex *fitBox2=new TLatex(0.5,0.8,"R_{inv,3} = 2.62 #pm 0.12"); //fitBox2->SetNDC(kTRUE); //fitBox2->Draw(); //TLatex *fitBox3=new TLatex(0.5,0.5,"<N_{rec,pions}>=103"); //fitBox3->SetNDC(kTRUE); //fitBox3->Draw(); //TPaveStats *c3Stats=(TPaveStats*)c3Fit1D_Base->FindObject("stats"); //c3Stats->SetX1NDC(0.15); //c3Stats->SetX2NDC(0.52); //c3Stats->SetY1NDC(.2); //c3Stats->SetY2NDC(.3); //c3Stats->Draw("same"); // The below 1D fit method has less well defined bin centers //TF1 *c3Fit1D=new TF1("c3Fit1D","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.))",0,1); //TF1 *c3Fit1D=new TF1("c3Fit1D","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.) * (1 + (-0.12/(6.*pow(2,1.5))*(8.*pow([2]*x/0.19733,3) - 12.*pow([2]*x/0.19733,1))) + (0.43/(24.*pow(2,2))*(16.*pow([2]*x/0.19733,4) -48.*pow([2]*x/0.19733,2) + 12))))",0,1); //TF1 *c3Fit1D=new TF1("c3Fit1D","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.) * (1 + ([3]/(6.*pow(2,1.5))*(8.*pow([2]*x/0.19733,3) - 12.*pow([2]*x/0.19733,1))) + ([4]/(24.*pow(2,2))*(16.*pow([2]*x/0.19733,4) -48.*pow([2]*x/0.19733,2) + 12))))",0,1); //c3Fit1D->SetLineColor(4); //c3Fit1D->SetParameter(0,1); c3Fit1D->SetParName(0,"N"); //c3Fit1D->SetParameter(1,0.5); c3Fit1D->SetParName(1,"#lambda_{3}"); //c3Fit1D->SetParameter(2,5); c3Fit1D->SetParName(2,"R_{inv}"); //c3Fit1D->SetParLimits(2,0.5,15); //c3Fit1D->SetLineColor(1); // //c3Fit1D->FixParameter(0,OutputPar_c3[0]); //c3Fit1D->FixParameter(1,OutputPar_c3[1]); //c3Fit1D->FixParameter(2,OutputPar_c3[2]); //c3Fit1D->SetParName(3,"coeff_{3}"); c3Fit1D->SetParName(4,"coeff_{4}"); //c3Fit1D->SetParameter(3,.24); c3Fit1D->SetParameter(3,.16); //c3_hist->Fit(c3Fit1D,"IMEN","",0.0,Q3Limit); //c3Fit1D->Draw("same"); //cout<<"1d fit: Chi2/NDF = "<<c3Fit1D->GetChisquare()/c3Fit1D->GetNDF()<<endl; //dentimesFit_c3->DrawCopy("l same"); // //c3Fit1D->FixParameter(0, fitcopy_c3->GetParameter(0)); //c3Fit1D->FixParameter(1, fitcopy_c3->GetParameter(1)); //c3Fit1D->FixParameter(2, fitcopy_c3->GetParameter(2)); //TF1 *c3Fit1D_2=new TF1("c3Fit1D_2","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.))",0,1); //c3Fit1D=new TF1("c3Fit1D","[0]*(1+[1]*exp(-pow([2]*x/0.19733,2)/2.) * pow(1 + ([3]/(6.*pow(2,1.5))*(8.*pow([2]*x/sqrt(3.)/0.19733,3) - 12.*pow([2]*x/sqrt(3.)/0.19733,1))) + ([4]/(24.*pow(2,2))*(16.*pow([2]*x/sqrt(3.)/0.19733,4) -48.*pow([2]*x/sqrt(3.)/0.19733,2) + 12)),3))",0,1); //c3Fit1D->FixParameter(0,OutputPar_c3[0]); //c3Fit1D->FixParameter(1,OutputPar_c3[1]); //c3Fit1D->FixParameter(2,OutputPar_c3[2]); //c3Fit1D->FixParameter(3,OutputPar_c3[3]); //c3Fit1D->FixParameter(4,OutputPar_c3[4]); //c3Fit1D->SetLineColor(4); //c3Fit1D->Draw("same"); //MomRes1d[1][0]->Draw(); //for(int i=1; i<50; i++) cout<<c3_hist->GetBinContent(i)<<", "; //cout<<endl; //for(int i=1; i<50; i++) cout<<c3_hist->GetBinError(i)<<", "; // pp M2, pi- Cterm1->Draw("same"); legend2->Draw("same"); /*if(SameCharge==kTRUE && MCcase==kFALSE){ TString *savename=new TString("C3c3_SC_"); if(CollisionType==0) savename->Append("PbPb_K"); if(CollisionType==1) savename->Append("pPb_K"); if(CollisionType==2) savename->Append("pp_K"); *savename += EDbin; savename->Append("_M"); *savename += Mbin; savename->Append(".eps"); can2->SaveAs(savename->Data()); }*/ //for(int ii=0; ii<10; ii++) cout<<c3_hist->GetBinContent(ii+1)<<", "; //Coul_GRiverside->Draw(); //Coul_Riverside->Draw("same"); /* TLegend *legend4 = new TLegend(.3,.8, .5,.95,NULL,"brNDC"); legend4->SetBorderSize(0); legend4->SetTextFont(42);//42 legend4->SetTextSize(.04);// small .03; large .06 legend4->SetFillColor(0); gPad->SetRightMargin(0.025); gPad->SetLeftMargin(0.1); gPad->SetTopMargin(0.02); TwoFrac=lambda_PM; OneFrac=pow(TwoFrac,.5), ThreeFrac=pow(TwoFrac,1.5); gPad->SetGridx(0); gPad->SetGridy(0); TH1D *c3_Extended = (TH1D*)Three_1d[CB1][CB2][CB3][SCBin][0]->Clone(); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][4], -( pow(1-OneFrac,3) + 3*OneFrac*pow(1-OneFrac,2) )); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][1], -(1-OneFrac)); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][2], -(1-OneFrac)); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][3], -(1-OneFrac)); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][4], (1-OneFrac)*3*(1-TwoFrac)); c3_Extended->Scale(1/ThreeFrac); // c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][1], -1/TwoFrac); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][2], -1/TwoFrac); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][3], -1/TwoFrac); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][4], 3*(1-TwoFrac)/TwoFrac); c3_Extended->Add(Three_1d[CB1][CB2][CB3][SCBin][4], 3); // c3_Extended->Divide(Three_1d[CB1][CB2][CB3][SCBin][4]); c3_Extended->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/#font[12]{c})"); c3_Extended->GetYaxis()->SetTitle("#font[12]{#bf{c}}_{3}^{#pm#pm#pm}"); c3_Extended->SetMinimum(0.9); c3_Extended->SetMaximum(2.0); c3_Extended->SetMarkerStyle(24); c3_Extended->SetMarkerColor(2); c3_Extended->SetLineColor(2); // TH1D *C3_Extended = (TH1D*)Three_1d[CB1][CB2][CB3][SCBin][0]->Clone(); C3_Extended->Divide(Three_1d[CB1][CB2][CB3][SCBin][4]); C3_Extended->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/c)"); C3_Extended->GetYaxis()->SetTitle("#font[12]{C}_{3}"); c3_Extended->GetXaxis()->SetTitleOffset(0.83); c3_Extended->GetYaxis()->SetTitleOffset(0.8); c3_Extended->GetYaxis()->SetTitleSize(0.05); c3_Extended->GetXaxis()->SetTitleSize(0.05); C3_Extended->SetMarkerStyle(20); C3_Extended->SetMarkerColor(4); c3_Extended->GetXaxis()->SetRangeUser(0,2); c3_Extended->SetMarkerStyle(20); c3_Extended->Draw(); //cout<<c3_Extended->GetBinError(40)<<" "<<Three_1d[CB1][CB2][CB3][SCBin][0]->GetBinError(40)/Three_1d[CB1][CB2][CB3][SCBin][4]->GetBinContent(40)<<" "<<Three_1d[CB1][CB2][CB3][SCBin][1]->GetBinError(40)/Three_1d[CB1][CB2][CB3][SCBin][4]->GetBinContent(40)<<endl; //legend4->AddEntry(c3_Extended,"c_{3}, #pi^{#pm}#pi^{#pm}#pi^{#pm}, Coulomb Uncorrected","p"); //legend4->AddEntry(C3_Extended,"C_{3}, #pi^{#pm}#pi^{#pm}#pi^{#pm}, Coulomb Uncorrected","p"); //legend4->Draw("same"); //Unity->Draw("same"); */ //Specif_System->Draw("same"); //Specif_KT3->Draw("same"); //Specif_FSI->Draw("same"); //Specif_Disclaimer->Draw("same"); /* gPad->SetGridx(0); gPad->SetGridy(0); int binLow=20; int binHigh=50; //TwoFrac=0.9; OneFrac=pow(TwoFrac,0.5); ThreeFrac=pow(TwoFrac,1.5); TH1D *K0sProjection = (TH1D*)(Three_3d[0][0][1][SCBin][0]->ProjectionY("K0sProjection",binLow,binHigh,binLow,binHigh)); TH1D *K0sProjection_term1 = (TH1D*)(Three_3d[0][0][1][SCBin][0]->ProjectionY("K0sProjection_term1",binLow,binHigh,binLow,binHigh)); TH1D *K0sProjection_term2 = (TH1D*)(Three_3d[0][0][1][SCBin][1]->ProjectionY("K0sProjection_term2",binLow,binHigh,binLow,binHigh)); TH1D *K0sProjection_term3 = (TH1D*)(Three_3d[0][0][1][SCBin][2]->ProjectionY("K0sProjection_term3",binLow,binHigh,binLow,binHigh)); TH1D *K0sProjection_term4 = (TH1D*)(Three_3d[0][0][1][SCBin][3]->ProjectionY("K0sProjection_term4",binLow,binHigh,binLow,binHigh)); TH1D *K0sProjection_term5 = (TH1D*)(Three_3d[0][0][1][SCBin][4]->ProjectionY("K0sProjection_term5",binLow,binHigh,binLow,binHigh)); K0sProjection_term1->Add(K0sProjection_term5, -( pow(1-OneFrac,3) + 3*OneFrac*pow(1-OneFrac,2) )); K0sProjection_term1->Add(K0sProjection_term2, -(1-OneFrac)); K0sProjection_term1->Add(K0sProjection_term3, -(1-OneFrac)); K0sProjection_term1->Add(K0sProjection_term4, -(1-OneFrac)); K0sProjection_term1->Add(K0sProjection_term5, (1-OneFrac)*3*(1-TwoFrac)); K0sProjection_term1->Scale(1/ThreeFrac); // K0sProjection_term1->Add(K0sProjection_term2, -1/TwoFrac); K0sProjection_term1->Add(K0sProjection_term3, -1/TwoFrac); K0sProjection_term1->Add(K0sProjection_term4, -1/TwoFrac); K0sProjection_term1->Add(K0sProjection_term5, 3*(1-TwoFrac)/TwoFrac); K0sProjection_term1->Add(K0sProjection_term5, 3); // K0sProjection->Divide(K0sProjection_term5); K0sProjection_term1->Divide(K0sProjection_term5); K0sProjection_term1->GetXaxis()->SetRangeUser(0,.5); K0sProjection_term1->SetMarkerStyle(24); K0sProjection_term1->SetMarkerColor(4); K0sProjection->SetMarkerStyle(20); K0sProjection->SetMarkerColor(4); K0sProjection->SetMinimum(0.92); K0sProjection->GetXaxis()->SetTitle("#font[12]{q_{13}^{#pm#mp}} (GeV/#font[12]{c})"); K0sProjection->GetYaxis()->SetTitle("#font[12]{C_{3}} or #font[12]{#bf{c}_{3}}"); K0sProjection->GetXaxis()->SetTitleOffset(1.2); K0sProjection->GetYaxis()->SetTitleOffset(1.3); K0sProjection->Draw(); // Sys errors TH1D *Sys_K0sProjection=(TH1D*)K0sProjection->Clone(); TH1D *Sys_K0sProjection_term1=(TH1D*)K0sProjection_term1->Clone(); Sys_K0sProjection->SetMarkerSize(0); Sys_K0sProjection_term1->SetMarkerSize(0); Sys_K0sProjection->SetFillColor(kBlue-10); Sys_K0sProjection_term1->SetFillColor(kBlue-10); Sys_K0sProjection->SetFillStyle(1000); Sys_K0sProjection_term1->SetFillStyle(1000); for(int i=0; i<Sys_K0sProjection->GetNbinsX(); i++) { Sys_K0sProjection->SetBinError(i+1, 0.01 * Sys_K0sProjection->GetBinContent(i+1)); Sys_K0sProjection_term1->SetBinError(i+1, 0.01 * Sys_K0sProjection_term1->GetBinContent(i+1)); } for(int i=0; i<Sys_K0sProjection->GetNbinsX(); i++) cout<<K0sProjection->GetBinContent(i+1)<<", "; cout<<endl; for(int i=0; i<Sys_K0sProjection->GetNbinsX(); i++) cout<<K0sProjection->GetBinError(i+1)<<", "; cout<<endl; for(int i=0; i<Sys_K0sProjection->GetNbinsX(); i++) cout<<K0sProjection_term1->GetBinContent(i+1)<<", "; cout<<endl; for(int i=0; i<Sys_K0sProjection->GetNbinsX(); i++) cout<<K0sProjection_term1->GetBinError(i+1)<<", "; cout<<endl; Sys_K0sProjection->Draw("E2 same"); Sys_K0sProjection_term1->Draw("E2 same"); K0sProjection->Draw("same"); K0sProjection_term1->Draw("same"); //legend4->AddEntry(K0sProjection,"#font[12]{C_{3}^{#pm#pm#mp}} projection","p"); //legend4->AddEntry(K0sProjection_term1,"#font[12]{#bf{c}_{3}^{#pm#pm#mp}} projection","p"); //legend4->Draw("same"); */ ////////////////////////////////////////////////////////////////////////////////////// if(SaveToFile){ TString *savefilename = new TString("OutFiles/OutFile"); if(!Gaussian) savefilename->Append("ExpFit"); if(CollisionType==0) savefilename->Append("PbPb"); else if(CollisionType==1) savefilename->Append("pPb"); else savefilename->Append("pp"); // if(MCcase) savefilename->Append("MonteCarlo"); // if(SameCharge) savefilename->Append("SC"); else savefilename->Append("MC"); // if(CHARGE==1) savefilename->Append("Pos"); else savefilename->Append("Neg"); // savefilename->Append("Kt3_"); *savefilename += EDbin+1; savefilename->Append("_M"); *savefilename += Mbin; savefilename->Append(".root"); TFile *savefile = new TFile(savefilename->Data(),"RECREATE"); can1->Write("can1"); C2_ss->Write("C2_ss"); C2_os->Write("C2_os"); fitC2ss_h->Write("fitC2ss_h"); fitC2ss_Base->Write("fitC2ss_Base"); fitC2ss_Expan->Write("fitC2ss_Expan"); Cterm1->Write("C3"); c3_hist->Write("c3"); Combinatorics_1d->Write("Combinatorics_1d"); c3Fit1D_Base->Write("c3Fit1D_Base"); c3Fit1D_Expan->Write("c3Fit1D_Expan"); c3Fit1D_ExpanSpline->Write("c3Fit1D_ExpanSpline"); gr_c3Spline->Write("gr_c3Spline"); MyMinuit_c3.Write("MyMinuit_c3"); // savefile->Close(); } } void ReadCoulCorrections(int index){ /////////////////////// TString *fname2 = new TString("KFile.root"); TFile *File=new TFile(fname2->Data(),"READ"); if(index>=12) cout<<"FSI index not acceptable in 2-particle Coulomb read"<<endl; TString *nameSS = new TString("K2ss_"); TString *nameOS = new TString("K2os_"); *nameSS += index; *nameOS += index; TString *nameSS_2 = new TString("K2ss_"); TString *nameOS_2 = new TString("K2os_"); *nameSS_2 += index+1; *nameOS_2 += index+1; /*if(index<9) {*nameSS_2 += index+1; *nameOS_2 += index+1;} else {*nameSS_2 += index; *nameOS_2 += index;}*/ TH1D *temp_ss = (TH1D*)File->Get(nameSS->Data()); TH1D *temp_os = (TH1D*)File->Get(nameOS->Data()); TH1D *temp_ss_2 = (TH1D*)File->Get(nameSS_2->Data()); TH1D *temp_os_2 = (TH1D*)File->Get(nameOS_2->Data()); CoulCorr2SS = (TH1D*)temp_ss->Clone(); CoulCorr2OS = (TH1D*)temp_os->Clone(); CoulCorr2SS_2 = (TH1D*)temp_ss_2->Clone(); CoulCorr2OS_2 = (TH1D*)temp_os_2->Clone(); // CoulCorr2SS->SetDirectory(0); CoulCorr2OS->SetDirectory(0); CoulCorr2SS_2->SetDirectory(0); CoulCorr2OS_2->SetDirectory(0); File->Close(); } double CoulCorr2(int chargeproduct, double Q2){// returns K2 int indexL=0; int indexH=0; double slope=0; double value1=1.0, value2=1.0; indexL = int(fabs(Q2 - CoulCorr2SS->GetXaxis()->GetBinWidth(1)/2.)/(CoulCorr2SS->GetXaxis()->GetBinWidth(1))); indexH = indexL+1; if(indexH >= CoulCorr2SS->GetNbinsX()) { return (Gamov(chargeproduct, Q2) * CoulCorr2SS->GetBinContent(CoulCorr2SS->GetNbinsX()-1) / Gamov(chargeproduct, CoulCorr2SS->GetXaxis()->GetBinCenter(CoulCorr2SS->GetNbinsX()-1))); } // Low index value if(chargeproduct==1){ slope = (CoulCorr2SS->GetBinContent(indexL+1) - CoulCorr2SS->GetBinContent(indexH+1)); slope /= (CoulCorr2SS->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2SS->GetXaxis()->GetBinCenter(indexH+1)); value1 = slope*(Q2 - CoulCorr2SS->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2SS->GetBinContent(indexL+1); }else { slope = (CoulCorr2OS->GetBinContent(indexL+1) - CoulCorr2OS->GetBinContent(indexH+1)); slope /= (CoulCorr2OS->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2OS->GetXaxis()->GetBinCenter(indexH+1)); value1 = slope*(Q2 - CoulCorr2OS->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2OS->GetBinContent(indexL+1); } // High index value if(chargeproduct==1){ slope = (CoulCorr2SS_2->GetBinContent(indexL+1) - CoulCorr2SS_2->GetBinContent(indexH+1)); slope /= (CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexH+1)); value2 = slope*(Q2 - CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2SS_2->GetBinContent(indexL+1); }else { slope = (CoulCorr2OS_2->GetBinContent(indexL+1) - CoulCorr2OS_2->GetBinContent(indexH+1)); slope /= (CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexH+1)); value2 = slope*(Q2 - CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2OS_2->GetBinContent(indexL+1); } if(value1>0 && value2>0){// interpolation only done between 9 and 12 Mbin. return (value1 + (value2-value1)/(2.));// always use the half-way point. }else if(value1>0){ return value1; }else if(value2>0){ return value2; }else return 1.0; // // old way with Gaussian sources for low mult systems /*if(indexH >= CoulCorr2SS->GetNbinsX()) { return (Gamov(chargeproduct, Q2) * CoulCorr2SS->GetBinContent(CoulCorr2SS->GetNbinsX()-1) / Gamov(chargeproduct, CoulCorr2SS->GetXaxis()->GetBinCenter(CoulCorr2SS->GetNbinsX()-1))); } if(chargeproduct==1){ slope = (CoulCorr2SS->GetBinContent(indexL+1) - CoulCorr2SS->GetBinContent(indexH+1)); slope /= (CoulCorr2SS->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2SS->GetXaxis()->GetBinCenter(indexH+1)); value1 = slope*(Q2 - CoulCorr2SS->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2SS->GetBinContent(indexL+1); }else { slope = (CoulCorr2OS->GetBinContent(indexL+1) - CoulCorr2OS->GetBinContent(indexH+1)); slope /= (CoulCorr2OS->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2OS->GetXaxis()->GetBinCenter(indexH+1)); value1 = slope*(Q2 - CoulCorr2OS->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2OS->GetBinContent(indexL+1); } if(Mbin_def<=9 || Mbin_def>12) return value1; else{// interpolate in K factor transition region if(chargeproduct==1){ slope = (CoulCorr2SS_2->GetBinContent(indexL+1) - CoulCorr2SS_2->GetBinContent(indexH+1)); slope /= (CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexH+1)); value2 = slope*(Q2 - CoulCorr2SS_2->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2SS_2->GetBinContent(indexL+1); }else { slope = (CoulCorr2OS_2->GetBinContent(indexL+1) - CoulCorr2OS_2->GetBinContent(indexH+1)); slope /= (CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexL+1) - CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexH+1)); value2 = slope*(Q2 - CoulCorr2OS_2->GetXaxis()->GetBinCenter(indexL+1)) + CoulCorr2OS_2->GetBinContent(indexL+1); } if(value1>0 && value2>0){// interpolation only done between 9 and 12 Mbin. return (value1 + (Mbin_def-9)*(value2-value1)/(3.)); }else if(value1>0){ return value1; }else if(value2>0){ return value2; }else return 1.0; } */ // } //---------------------------------------------------------------------- //________________________________________________________________________ void fcnC2_Global(int& npar, double* deriv, double& f, double par[], int flag){ double qinvSS=0; double Rch=par[3]/FmToGeV; double CSS=0, COS=0; double SumChi2=0; double Expan=0; NFitPoints_C2global=0; if(LinkN) par[9]=par[0];// Link N factors for(int i=1; i<BINRANGE_C2global; i++){ qinvSS = BinCenters[i]; if(qinvSS > Q2Limit) continue; // double arg=qinvSS*Rch; if(Gaussian){// Edgeworth expansion Expan = 1 + par[5]/(6.*pow(2,1.5)) * (8.*pow(arg,3) - 12.*pow(arg,1)); Expan += par[6]/(24.*pow(2,2)) * (16.*pow(arg,4) -48.*pow(arg,2) + 12); Expan += par[7]/(120.*pow(2,2.5)) * (32.*pow(arg,5) - 160.*pow(arg,3) + 120*arg); Expan += par[8]/(720.*pow(2,3)) * (64.*pow(arg,6) - 480.*pow(arg,4) + 720.*pow(arg,2) - 120); }else{ Expan = 1 + par[5] * (arg - 1); Expan += par[6]/2. * (pow(arg,2) - 4*arg + 2); Expan += par[7] * 0.; Expan += par[8] * 0.; } //Expan = 1.0; // // old way without undilution /*CSS = 1 + exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); CSS *= par[1]*K2SS[i]; CSS += 1-par[1]; CSS *= par[0];*/ // undilution method //CSS = 1 + par[1]*exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); //CSS *= par[0]; // undilution method with correct normalization location (was previously applied to C2^QS) CSS = 1 + par[1]*exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); CSS *= lambda_PM * K2SS[i]; CSS += 1-lambda_PM; CSS *= par[0]; // COS = 1; COS *= par[1]*K2OS[i]; COS += 1-par[1]; COS *= par[9];// different Norm // if(C2ssFitting[i] > 0){ if(IncludeSS) { //double errorSS = sqrt(pow( (CSS/par[0] - (1-par[1]))/K2SS[i] * K2SS_e[i] * par[0],2) + pow(C2ssFitting_e[i],2)); double errorSS = sqrt(pow(C2ssFitting_e[i],2));// new way with undilution already done in C2ssFitting[i] SumChi2 += pow((CSS - C2ssFitting[i])/errorSS,2); NFitPoints_C2global++; } } if(IncludeOS) { double errorOS = sqrt(pow( (COS/par[9] - (1-par[1]))/K2OS[i] * K2OS_e[i] * par[9],2) + pow(C2osFitting_e[i],2)); SumChi2 += pow((COS - C2osFitting[i])/errorOS,2); NFitPoints_C2global++; } } f = SumChi2; Chi2_C2global = f; } //________________________________________________________________________ double C2ssFitFunction(Double_t *x, Double_t *par) { double Rch=par[3]/FmToGeV; int qbin = int(fabs(x[0]/BinWidthQ2)); double qinvSS = BinCenters[qbin]; double arg=qinvSS*Rch; double Expan = 1; if(Gaussian){// Edgeworth expansion Expan += par[5]/(6.*pow(2,1.5))*(8.*pow(arg,3) - 12.*pow(arg,1)); Expan += par[6]/(24.*pow(2,2))*(16.*pow(arg,4) -48.*pow(arg,2) + 12); Expan += par[7]/(120.*pow(2,2.5))*(32.*pow(arg,5) - 160.*pow(arg,3) + 120*arg); Expan += par[8]/(720.*pow(2,3))*(64.*pow(arg,6) - 480.*pow(arg,4) + 720.*pow(arg,2) - 120); }else{// Laguerre Expan = 1 + par[5] * (arg - 1); Expan += par[6]/2. * (pow(arg,2) - 4*arg + 2); Expan += par[7] * 0.; Expan += par[8] * 0.; } // old way without undilution /*double CSS = 1 + exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); double K2=1.0; if(qbin < BINRANGE_C2global) K2=K2SS[qbin]; CSS *= par[1]*K2; CSS += 1-par[1]; CSS *= par[0];*/ // undilution method //double CSS = 1 + par[1]*exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); //CSS *= par[0]; // undilution method with correct normalization location (was previously applied to C2^QS) double CSS = 1 + par[1]*exp(-pow(Rch*qinvSS,ExpPower))*pow(Expan,2); CSS *= lambda_PM * K2SS[qbin]; CSS += 1-lambda_PM; CSS *= par[0]; return CSS; } //________________________________________________________________________ double C2osFitFunction(Double_t *x, Double_t *par) { if(LinkN) par[9]=par[0];// Link N factors int qbin = int(fabs(x[0]/BinWidthQ2)); //double qinvOS = BinCenters[qbin]; double COS = 1; double K2=1.0; if(qbin < BINRANGE_C2global) K2=K2OS[qbin]; COS *= par[1]*K2; COS += 1-par[1]; COS *= par[9]; return COS; } //__________________________________________________________________________ void fcn_c3(int& npar, double* deriv, double& f, double par[], int flag){ double q12=0, q13=0, q23=0; double Expan12=0, Expan13=0, Expan23=0; double C=0; double Rch=par[2]/FmToGeV; double SumChi2=0; //double lnL=0, term1=0, term2=0; NFitPoints_c3=0; //double SumChi2_test=0; for(int i=0; i<=BINLIMIT_3; i++){// q12 for(int j=0; j<=BINLIMIT_3; j++){// q13 for(int k=0; k<=BINLIMIT_3; k++){// q23 if(B_3[i][j][k] == 0) continue; if(A_3[i][j][k] == 0) continue; if(A_3_e[i][j][k] == 0) continue; q12 = BinCenters[i]; q13 = BinCenters[j]; q23 = BinCenters[k]; double q3 = sqrt(pow(q12,2)+pow(q13,2)+pow(q23,2)); if(q3 > Q3Limit) continue; // double arg12 = q12*Rch; double arg13 = q13*Rch; double arg23 = q23*Rch; if(Gaussian){// Edgeworth expansion Expan12 = 1 + par[3]/(6.*pow(2,1.5))*(8*pow(arg12,3) - 12*pow(arg12,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg12,4) -48*pow(arg12,2) + 12); Expan13 = 1 + par[3]/(6.*pow(2,1.5))*(8*pow(arg13,3) - 12*pow(arg13,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg13,4) -48*pow(arg13,2) + 12); Expan23 = 1 + par[3]/(6.*pow(2,1.5))*(8*pow(arg23,3) - 12*pow(arg23,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg23,4) -48*pow(arg23,2) + 12); }else{// Laguerre expansion Expan12 = 1 + par[3]*(arg12 - 1) + par[4]/2.*(pow(arg12,2) - 4*arg12 + 2); Expan13 = 1 + par[3]*(arg13 - 1) + par[4]/2.*(pow(arg13,2) - 4*arg13 + 2); Expan23 = 1 + par[3]*(arg23 - 1) + par[4]/2.*(pow(arg23,2) - 4*arg23 + 2); } // C = 1 + par[1]*exp(-(pow(arg12,ExpPower)+pow(arg13,ExpPower)+pow(arg23,ExpPower))/2.)*Expan12*Expan13*Expan23; C *= par[0];// Norm double error = pow(A_3_e[i][j][k]/B_3[i][j][k],2); error += pow(sqrt(B_3[i][j][k])*A_3[i][j][k]/pow(B_3[i][j][k],2),2); error = sqrt(error); SumChi2 += pow( (C - A_3[i][j][k]/B_3[i][j][k])/error,2); //if(q3<0.05) SumChi2_test += pow( (C - A_3[i][j][k]/B_3[i][j][k])/error,2); // /*if(A_3[i][j][k] >= 1) term1 = C*(A_3[i][j][k]+B_3[i][j][k])/(A_3[i][j][k]*(C+1)); else term1 = 0; term2 = (A_3[i][j][k]+B_3[i][j][k])/(B_3[i][j][k]*(C+1)); if(term1 > 0.0 && term2 > 0.0){ lnL += A_3[i][j][k]*log(term1) + B_3[i][j][k]*log(term2); }else if(term1==0 && term2 > 0.0){ lnL += B_3[i][j][k]*log(term2); }else {cout<<"WARNING -- term1 and term2 are negative"<<endl;} */ float DegenerateCount=1.; //if(i==j && i==k) DegenerateCount=1; //else if(i==j || i==k || j==k) DegenerateCount=3;//3 //else DegenerateCount=6;//6 NFitPoints_c3 += 1/DegenerateCount; } } } //f = -2.0*lnL;// log-liklihood minimization f = SumChi2;// Chi2 minimization Chi2_c3 = f; //Chi2_c3 = SumChi2_test; } //________________________________________________________________________ double Dfitfunction_c3(Double_t *x, Double_t *par) { double Rch = par[2]/FmToGeV; double arg12 = x[0]*Rch; double arg13 = x[1]*Rch; double arg23 = x[2]*Rch; double Expan12 = 1, Expan13 = 1, Expan23 = 1; if(Gaussian){// Edworth Expansion Expan12 += par[3]/(6.*pow(2,1.5))*(8*pow(arg12,3) - 12*pow(arg12,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg12,4) -48*pow(arg12,2) + 12); Expan13 += par[3]/(6.*pow(2,1.5))*(8*pow(arg13,3) - 12*pow(arg13,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg13,4) -48*pow(arg13,2) + 12); Expan23 += par[3]/(6.*pow(2,1.5))*(8*pow(arg23,3) - 12*pow(arg23,1)) + par[4]/(24.*pow(2,2))*(16*pow(arg23,4) -48*pow(arg23,2) + 12); }else{// Laguerre Expansion Expan12 = 1 + par[3]*(arg12 - 1) + par[4]/2.*(pow(arg12,2) - 4*arg12 + 2); Expan13 = 1 + par[3]*(arg13 - 1) + par[4]/2.*(pow(arg13,2) - 4*arg13 + 2); Expan23 = 1 + par[3]*(arg23 - 1) + par[4]/2.*(pow(arg23,2) - 4*arg23 + 2); } // double C = 1 + par[1]*exp(-(pow(arg12,ExpPower)+pow(arg13,ExpPower)+pow(arg23,ExpPower))/2.)*Expan12*Expan13*Expan23; C *= par[0];// Norm return C; } //________________________________________________________________________ double CoulCorrGRS(bool SC, double Q_12, double Q_13, double Q_23){ int index12L = int(fabs(Q_12 - CoulCorr2SS->GetXaxis()->GetBinWidth(1)/2.)/(CoulCorr2SS->GetXaxis()->GetBinWidth(1))); int index12H = index12L+1; int index13L = int(fabs(Q_13 - CoulCorr2SS->GetXaxis()->GetBinWidth(1)/2.)/(CoulCorr2SS->GetXaxis()->GetBinWidth(1))); int index13H = index13L+1; int index23L = int(fabs(Q_23 - CoulCorr2SS->GetXaxis()->GetBinWidth(1)/2.)/(CoulCorr2SS->GetXaxis()->GetBinWidth(1))); int index23H = index23L+1; if(Tricubic){ // Tricubic Interpolation double arr[4][4][4]={{{0}}}; for(int x=0; x<4; x++){ for(int y=0; y<4; y++){ for(int z=0; z<4; z++){ if(SC){ arr[x][y][z] = CoulCorr2SS->GetBinContent((index12L)+x)*CoulCorr2SS->GetBinContent((index23L)+y)*CoulCorr2SS->GetBinContent((index13L)+z); }else{ arr[x][y][z] = CoulCorr2SS->GetBinContent((index12L)+x)*CoulCorr2OS->GetBinContent((index23L)+y)*CoulCorr2OS->GetBinContent((index13L)+z); } } } } return tricubicInterpolate(arr, Q_12, Q_23, Q_13); }else{ // Trilinear Interpolation. See for instance: https://en.wikipedia.org/wiki/Trilinear_interpolation // double value1=1.0, value2=1.0; // double xd = (Q_12-CoulCorr2SS->GetXaxis()->GetBinCenter(index12L+1)); xd /= (CoulCorr2SS->GetXaxis()->GetBinCenter(index12H+1)-CoulCorr2SS->GetXaxis()->GetBinCenter(index12L+1)); double yd = (Q_13-CoulCorr2SS->GetXaxis()->GetBinCenter(index13L+1)); yd /= (CoulCorr2SS->GetXaxis()->GetBinCenter(index13H+1)-CoulCorr2SS->GetXaxis()->GetBinCenter(index13L+1)); double zd = (Q_23-CoulCorr2SS->GetXaxis()->GetBinCenter(index23L+1)); zd /= (CoulCorr2SS->GetXaxis()->GetBinCenter(index23H+1)-CoulCorr2SS->GetXaxis()->GetBinCenter(index23L+1)); // if(SC){ double c00 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2SS->GetBinContent(index13L+1)*CoulCorr2SS->GetBinContent(index23L+1)*(1-xd); c00 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2SS->GetBinContent(index13L+1)*CoulCorr2SS->GetBinContent(index23L+1)*xd; double c10 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2SS->GetBinContent(index13H+1)*CoulCorr2SS->GetBinContent(index23L+1)*(1-xd); c10 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2SS->GetBinContent(index13H+1)*CoulCorr2SS->GetBinContent(index23L+1)*xd; double c01 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2SS->GetBinContent(index13L+1)*CoulCorr2SS->GetBinContent(index23H+1)*(1-xd); c01 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2SS->GetBinContent(index13L+1)*CoulCorr2SS->GetBinContent(index23H+1)*xd; double c11 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2SS->GetBinContent(index13H+1)*CoulCorr2SS->GetBinContent(index23H+1)*(1-xd); c11 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2SS->GetBinContent(index13H+1)*CoulCorr2SS->GetBinContent(index23H+1)*xd; // double c0 = c00*(1-yd) + c10*yd; double c1 = c01*(1-yd) + c11*yd; value1 = (c0*(1-zd) + c1*zd); }else{ double c00 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2OS->GetBinContent(index13L+1)*CoulCorr2OS->GetBinContent(index23L+1)*(1-xd); c00 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2OS->GetBinContent(index13L+1)*CoulCorr2OS->GetBinContent(index23L+1)*xd; double c10 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2OS->GetBinContent(index13H+1)*CoulCorr2OS->GetBinContent(index23L+1)*(1-xd); c10 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2OS->GetBinContent(index13H+1)*CoulCorr2OS->GetBinContent(index23L+1)*xd; double c01 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2OS->GetBinContent(index13L+1)*CoulCorr2OS->GetBinContent(index23H+1)*(1-xd); c01 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2OS->GetBinContent(index13L+1)*CoulCorr2OS->GetBinContent(index23H+1)*xd; double c11 = CoulCorr2SS->GetBinContent(index12L+1)*CoulCorr2OS->GetBinContent(index13H+1)*CoulCorr2OS->GetBinContent(index23H+1)*(1-xd); c11 += CoulCorr2SS->GetBinContent(index12H+1)*CoulCorr2OS->GetBinContent(index13H+1)*CoulCorr2OS->GetBinContent(index23H+1)*xd; // double c0 = c00*(1-yd) + c10*yd; double c1 = c01*(1-yd) + c11*yd; value1 = (c0*(1-zd) + c1*zd); } //if(Mbin_def<=9 || Mbin_def>12) return value1;// old way for v5 // // interpolation for K factor tansition if(SC){ double c00 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2SS_2->GetBinContent(index13L+1)*CoulCorr2SS_2->GetBinContent(index23L+1)*(1-xd); c00 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2SS_2->GetBinContent(index13L+1)*CoulCorr2SS_2->GetBinContent(index23L+1)*xd; double c10 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2SS_2->GetBinContent(index13H+1)*CoulCorr2SS_2->GetBinContent(index23L+1)*(1-xd); c10 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2SS_2->GetBinContent(index13H+1)*CoulCorr2SS_2->GetBinContent(index23L+1)*xd; double c01 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2SS_2->GetBinContent(index13L+1)*CoulCorr2SS_2->GetBinContent(index23H+1)*(1-xd); c01 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2SS_2->GetBinContent(index13L+1)*CoulCorr2SS_2->GetBinContent(index23H+1)*xd; double c11 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2SS_2->GetBinContent(index13H+1)*CoulCorr2SS_2->GetBinContent(index23H+1)*(1-xd); c11 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2SS_2->GetBinContent(index13H+1)*CoulCorr2SS_2->GetBinContent(index23H+1)*xd; // double c0 = c00*(1-yd) + c10*yd; double c1 = c01*(1-yd) + c11*yd; value2 = (c0*(1-zd) + c1*zd); }else{ double c00 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2OS_2->GetBinContent(index13L+1)*CoulCorr2OS_2->GetBinContent(index23L+1)*(1-xd); c00 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2OS_2->GetBinContent(index13L+1)*CoulCorr2OS_2->GetBinContent(index23L+1)*xd; double c10 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2OS_2->GetBinContent(index13H+1)*CoulCorr2OS_2->GetBinContent(index23L+1)*(1-xd); c10 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2OS_2->GetBinContent(index13H+1)*CoulCorr2OS_2->GetBinContent(index23L+1)*xd; double c01 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2OS_2->GetBinContent(index13L+1)*CoulCorr2OS_2->GetBinContent(index23H+1)*(1-xd); c01 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2OS_2->GetBinContent(index13L+1)*CoulCorr2OS_2->GetBinContent(index23H+1)*xd; double c11 = CoulCorr2SS_2->GetBinContent(index12L+1)*CoulCorr2OS_2->GetBinContent(index13H+1)*CoulCorr2OS_2->GetBinContent(index23H+1)*(1-xd); c11 += CoulCorr2SS_2->GetBinContent(index12H+1)*CoulCorr2OS_2->GetBinContent(index13H+1)*CoulCorr2OS_2->GetBinContent(index23H+1)*xd; // double c0 = c00*(1-yd) + c10*yd; double c1 = c01*(1-yd) + c11*yd; value2 = (c0*(1-zd) + c1*zd); } if(value1>0 && value2>0){ return (value1 + (value2-value1)/(2.));// always take the half-way point //return (value1 + (Mbin_def-9)*(value2-value1)/(3.));// old way for v5 }else if(value1>0){ return value1; }else if(value2>0){ return value2; }else return 1.0; } } double Gamov(int chargeProduct, double qinv){ double arg = chargeProduct*2.*PI/(BohrR*qinv/2.); return arg/(exp(arg)-1); } void ReadMomResFile(int Mbin){ int MRCindex_L=0, MRCindex_H=0; float MRCindex_weight=0; if(Mbin<=2) {MRCindex_L=0; MRCindex_H=1; MRCindex_weight = Mbin/2.;} else if(Mbin<=6) {MRCindex_L=1; MRCindex_H=2; MRCindex_weight = fabs(Mbin-3)/3.;} else if(Mbin<=11) {MRCindex_L=2; MRCindex_H=3; MRCindex_weight = fabs(Mbin-7)/4.;} else {MRCindex_L=3; MRCindex_H=4; MRCindex_weight = fabs(Mbin-12)/7.;} TH1D *temp_L[2][5]; TH1D *temp_H[2][5]; TH1D *tempC2_L[2]; TH1D *tempC2_H[2]; TString *momresfilenameL = new TString("MomResFile_M"); *momresfilenameL += MRCindex_L; momresfilenameL->Append(".root"); TFile *MomResFileL = new TFile(momresfilenameL->Data(),"READ"); TString *names1D_L[2][5];// SC/MC, term# for(int ChProd=0; ChProd<2; ChProd++){ for(int term=0; term<5; term++){ // if(ChProd==0) {names1D_L[ChProd][term] = new TString("MRC3_SC_term");} else {names1D_L[ChProd][term] = new TString("MRC3_MC_term");} *names1D_L[ChProd][term] += term+1; temp_L[ChProd][term] = (TH1D*)MomResFileL->Get(names1D_L[ChProd][term]->Data()); temp_L[ChProd][term]->SetDirectory(0); } TString *C2MRCname=new TString("MomResHisto_"); if(ChProd==0) C2MRCname->Append("pp"); else C2MRCname->Append("mp"); tempC2_L[ChProd]=(TH1D*)(((TH2D*)(MomResFileL->Get(C2MRCname->Data())))->ProjectionY("C2MRCproL",MRC2index,MRC2index)); tempC2_L[ChProd]->SetDirectory(0); } // TString *momresfilenameH = new TString("MomResFile_M"); *momresfilenameH += MRCindex_H; momresfilenameH->Append(".root"); TFile *MomResFileH = new TFile(momresfilenameH->Data(),"READ"); TString *names1D_H[2][5];// SC/MC, term# for(int ChProd=0; ChProd<2; ChProd++){ for(int term=0; term<5; term++){ // if(ChProd==0) {names1D_H[ChProd][term] = new TString("MRC3_SC_term");} else {names1D_H[ChProd][term] = new TString("MRC3_MC_term");} *names1D_H[ChProd][term] += term+1; temp_H[ChProd][term] = (TH1D*)MomResFileH->Get(names1D_H[ChProd][term]->Data()); temp_H[ChProd][term]->SetDirectory(0); MomRes1d[ChProd][term] = (TH1D*)MomResFileH->Get(names1D_H[ChProd][term]->Data()); MomRes1d[ChProd][term]->SetDirectory(0); } TString *C2MRCname=new TString("MomResHisto_"); if(ChProd==0) C2MRCname->Append("pp"); else C2MRCname->Append("mp"); tempC2_H[ChProd]=(TH1D*)(((TH2D*)(MomResFileH->Get(C2MRCname->Data())))->ProjectionY("C2MRCproH",MRC2index,MRC2index)); tempC2_H[ChProd]->SetDirectory(0); TString *name=new TString("MomResC2_"); *name += ChProd; MomResC2[ChProd] = (TH1D*)(((TH2D*)(MomResFileH->Get(C2MRCname->Data())))->ProjectionY(name->Data(),MRC2index,MRC2index)); MomResC2[ChProd]->SetDirectory(0); } for(int ChProd=0; ChProd<2; ChProd++){ // C3 MRC for(int term=0; term<5; term++){ for(int bin=1; bin<=temp_H[ChProd][term]->GetNbinsX(); bin++){ double value=1; if(temp_L[ChProd][term]->GetBinContent(bin)>0 && temp_H[ChProd][term]->GetBinContent(bin)>0){// both have entries value = temp_L[ChProd][term]->GetBinContent(bin) + MRCindex_weight * temp_H[ChProd][term]->GetBinContent(bin); }else if(temp_L[ChProd][term]->GetBinContent(bin)>0){ value = temp_L[ChProd][term]->GetBinContent(bin); }else if(temp_H[ChProd][term]->GetBinContent(bin)>0){ value = temp_H[ChProd][term]->GetBinContent(bin); }else value=1.0; MomRes1d[ChProd][term]->SetBinContent(bin,value); } } // C2 MRC for(int bin=1; bin<=tempC2_H[ChProd]->GetNbinsX(); bin++){ double value=1; if(tempC2_L[ChProd]->GetBinContent(bin)>0 && tempC2_H[ChProd]->GetBinContent(bin)>0){// both have entries value = tempC2_L[ChProd]->GetBinContent(bin)*(1-MRCindex_weight) + MRCindex_weight * tempC2_H[ChProd]->GetBinContent(bin); }else if(tempC2_L[ChProd]->GetBinContent(bin)>0){ value = tempC2_L[ChProd]->GetBinContent(bin); }else if(tempC2_H[ChProd]->GetBinContent(bin)>0){ value = tempC2_H[ChProd]->GetBinContent(bin); }else value=1.0; MomResC2[ChProd]->SetBinContent(bin,value); } } for(int ChProd=0; ChProd<2; ChProd++){ for(int term=0; term<5; term++){ for(int i=0; i<MomRes1d[0][0]->GetNbinsX(); i++){ if(SourceType==0 && Mbin<10){ if(MomRes1d[ChProd][term]->GetXaxis()->GetBinCenter(i+1) > 0.1) MomRes1d[ChProd][term]->SetBinContent(i+1, 1.0); }else if(SourceType==0 && Mbin>=10){ if(MomRes1d[ChProd][term]->GetXaxis()->GetBinCenter(i+1) > 0.2) MomRes1d[ChProd][term]->SetBinContent(i+1, 1.0); }else{ if(MomRes1d[ChProd][term]->GetXaxis()->GetBinCenter(i+1) > 0.5) MomRes1d[ChProd][term]->SetBinContent(i+1, 1.0); } } } } MomResFileL->Close(); MomResFileH->Close(); } double cubicInterpolate (double p[4], double x) { return p[1] + 0.5 * x*(p[2] - p[0] + x*(2.0*p[0] - 5.0*p[1] + 4.0*p[2] - p[3] + x*(3.0*(p[1] - p[2]) + p[3] - p[0])));// Paulinternet } double bicubicInterpolate (double p[4][4], double x, double y) { double arr[4]; arr[0] = cubicInterpolate(p[0], y); arr[1] = cubicInterpolate(p[1], y); arr[2] = cubicInterpolate(p[2], y); arr[3] = cubicInterpolate(p[3], y); return cubicInterpolate(arr, x); } double tricubicInterpolate (double p[4][4][4], double x, double y, double z) { double arr[4]; arr[0] = bicubicInterpolate(p[0], y, z); arr[1] = bicubicInterpolate(p[1], y, z); arr[2] = bicubicInterpolate(p[2], y, z); arr[3] = bicubicInterpolate(p[3], y, z); return cubicInterpolate(arr, x); } //____________________________________________________________________________________________________ void ReadMuonCorrections(int ct, int mbin){ TString *name = new TString("MuonCorrection_"); if(ct==0){// PbPb if(mbin<=1) *name += 0; else if(mbin<=3) *name += 1; else if(mbin<=5) *name += 2; else if(mbin<=7) *name += 3; else if(mbin<10) *name += 4; else *name += 5; }else *name += 6; name->Append(".root"); TFile *muonfile=new TFile(name->Data(),"READ"); C2muonCorrection_SC = (TH1D*)muonfile->Get("C2muonCorrection_SC"); C2muonCorrection_MC = (TH1D*)muonfile->Get("C2muonCorrection_MC"); if(SameCharge_def) C3muonCorrection = (TH1D*)muonfile->Get("C3muonCorrection_SC"); else C3muonCorrection = (TH1D*)muonfile->Get("C3muonCorrection_MC"); // C2muonCorrection_SC->SetDirectory(0); C2muonCorrection_MC->SetDirectory(0); C3muonCorrection->SetDirectory(0); // muonfile->Close(); } //____________________________________________________________________________________________________ void DrawALICELogo(Bool_t prel, Float_t x1, Float_t y1, Float_t x2, Float_t y2) { // revision on July 24th or 25th 2012 // correct for aspect ratio of figure plus aspect ratio of pad (coordinates are NDC!) x2 = x1 + (y2 - y1) * (466. / 523) * gPad->GetWh() * gPad->GetHNDC() / (gPad->GetWNDC() * gPad->GetWw()); // Printf("%f %f %f %f", x1, x2, y1, y2); TPad *myPadLogo = new TPad("myPadLogo", "Pad for ALICE Logo", x1, y1, x2, y2); myPadLogo->SetLeftMargin(0); myPadLogo->SetTopMargin(0); myPadLogo->SetRightMargin(0); myPadLogo->SetBottomMargin(0); myPadLogo->Draw(); myPadLogo->cd(); TASImage *myAliceLogo = new TASImage((prel) ? "~/Pictures/2011-Nov-24-ALICE_PRELIMARY_logo_BLACK_small_usage_design.eps" : "~/Pictures/2011-Nov-24-ALICE_PERFORMANCE_logo_BLACK_small_usage_design.eps"); myAliceLogo->Draw(); }
Plot_PDCumulants_TPR.C:1
Plot_PDCumulants_TPR.C:2
Plot_PDCumulants_TPR.C:3
Plot_PDCumulants_TPR.C:4
Plot_PDCumulants_TPR.C:5
Plot_PDCumulants_TPR.C:6
Plot_PDCumulants_TPR.C:7
Plot_PDCumulants_TPR.C:8
Plot_PDCumulants_TPR.C:9
Plot_PDCumulants_TPR.C:10
Plot_PDCumulants_TPR.C:11
Plot_PDCumulants_TPR.C:12
Plot_PDCumulants_TPR.C:13
Plot_PDCumulants_TPR.C:14
Plot_PDCumulants_TPR.C:15
Plot_PDCumulants_TPR.C:16
Plot_PDCumulants_TPR.C:17
Plot_PDCumulants_TPR.C:18
Plot_PDCumulants_TPR.C:19
Plot_PDCumulants_TPR.C:20
Plot_PDCumulants_TPR.C:21
Plot_PDCumulants_TPR.C:22
Plot_PDCumulants_TPR.C:23
Plot_PDCumulants_TPR.C:24
Plot_PDCumulants_TPR.C:25
Plot_PDCumulants_TPR.C:26
Plot_PDCumulants_TPR.C:27
Plot_PDCumulants_TPR.C:28
Plot_PDCumulants_TPR.C:29
Plot_PDCumulants_TPR.C:30
Plot_PDCumulants_TPR.C:31
Plot_PDCumulants_TPR.C:32
Plot_PDCumulants_TPR.C:33
Plot_PDCumulants_TPR.C:34
Plot_PDCumulants_TPR.C:35
Plot_PDCumulants_TPR.C:36
Plot_PDCumulants_TPR.C:37
Plot_PDCumulants_TPR.C:38
Plot_PDCumulants_TPR.C:39
Plot_PDCumulants_TPR.C:40
Plot_PDCumulants_TPR.C:41
Plot_PDCumulants_TPR.C:42
Plot_PDCumulants_TPR.C:43
Plot_PDCumulants_TPR.C:44
Plot_PDCumulants_TPR.C:45
Plot_PDCumulants_TPR.C:46
Plot_PDCumulants_TPR.C:47
Plot_PDCumulants_TPR.C:48
Plot_PDCumulants_TPR.C:49
Plot_PDCumulants_TPR.C:50
Plot_PDCumulants_TPR.C:51
Plot_PDCumulants_TPR.C:52
Plot_PDCumulants_TPR.C:53
Plot_PDCumulants_TPR.C:54
Plot_PDCumulants_TPR.C:55
Plot_PDCumulants_TPR.C:56
Plot_PDCumulants_TPR.C:57
Plot_PDCumulants_TPR.C:58
Plot_PDCumulants_TPR.C:59
Plot_PDCumulants_TPR.C:60
Plot_PDCumulants_TPR.C:61
Plot_PDCumulants_TPR.C:62
Plot_PDCumulants_TPR.C:63
Plot_PDCumulants_TPR.C:64
Plot_PDCumulants_TPR.C:65
Plot_PDCumulants_TPR.C:66
Plot_PDCumulants_TPR.C:67
Plot_PDCumulants_TPR.C:68
Plot_PDCumulants_TPR.C:69
Plot_PDCumulants_TPR.C:70
Plot_PDCumulants_TPR.C:71
Plot_PDCumulants_TPR.C:72
Plot_PDCumulants_TPR.C:73
Plot_PDCumulants_TPR.C:74
Plot_PDCumulants_TPR.C:75
Plot_PDCumulants_TPR.C:76
Plot_PDCumulants_TPR.C:77
Plot_PDCumulants_TPR.C:78
Plot_PDCumulants_TPR.C:79
Plot_PDCumulants_TPR.C:80
Plot_PDCumulants_TPR.C:81
Plot_PDCumulants_TPR.C:82
Plot_PDCumulants_TPR.C:83
Plot_PDCumulants_TPR.C:84
Plot_PDCumulants_TPR.C:85
Plot_PDCumulants_TPR.C:86
Plot_PDCumulants_TPR.C:87
Plot_PDCumulants_TPR.C:88
Plot_PDCumulants_TPR.C:89
Plot_PDCumulants_TPR.C:90
Plot_PDCumulants_TPR.C:91
Plot_PDCumulants_TPR.C:92
Plot_PDCumulants_TPR.C:93
Plot_PDCumulants_TPR.C:94
Plot_PDCumulants_TPR.C:95
Plot_PDCumulants_TPR.C:96
Plot_PDCumulants_TPR.C:97
Plot_PDCumulants_TPR.C:98
Plot_PDCumulants_TPR.C:99
Plot_PDCumulants_TPR.C:100
Plot_PDCumulants_TPR.C:101
Plot_PDCumulants_TPR.C:102
Plot_PDCumulants_TPR.C:103
Plot_PDCumulants_TPR.C:104
Plot_PDCumulants_TPR.C:105
Plot_PDCumulants_TPR.C:106
Plot_PDCumulants_TPR.C:107
Plot_PDCumulants_TPR.C:108
Plot_PDCumulants_TPR.C:109
Plot_PDCumulants_TPR.C:110
Plot_PDCumulants_TPR.C:111
Plot_PDCumulants_TPR.C:112
Plot_PDCumulants_TPR.C:113
Plot_PDCumulants_TPR.C:114
Plot_PDCumulants_TPR.C:115
Plot_PDCumulants_TPR.C:116
Plot_PDCumulants_TPR.C:117
Plot_PDCumulants_TPR.C:118
Plot_PDCumulants_TPR.C:119
Plot_PDCumulants_TPR.C:120
Plot_PDCumulants_TPR.C:121
Plot_PDCumulants_TPR.C:122
Plot_PDCumulants_TPR.C:123
Plot_PDCumulants_TPR.C:124
Plot_PDCumulants_TPR.C:125
Plot_PDCumulants_TPR.C:126
Plot_PDCumulants_TPR.C:127
Plot_PDCumulants_TPR.C:128
Plot_PDCumulants_TPR.C:129
Plot_PDCumulants_TPR.C:130
Plot_PDCumulants_TPR.C:131
Plot_PDCumulants_TPR.C:132
Plot_PDCumulants_TPR.C:133
Plot_PDCumulants_TPR.C:134
Plot_PDCumulants_TPR.C:135
Plot_PDCumulants_TPR.C:136
Plot_PDCumulants_TPR.C:137
Plot_PDCumulants_TPR.C:138
Plot_PDCumulants_TPR.C:139
Plot_PDCumulants_TPR.C:140
Plot_PDCumulants_TPR.C:141
Plot_PDCumulants_TPR.C:142
Plot_PDCumulants_TPR.C:143
Plot_PDCumulants_TPR.C:144
Plot_PDCumulants_TPR.C:145
Plot_PDCumulants_TPR.C:146
Plot_PDCumulants_TPR.C:147
Plot_PDCumulants_TPR.C:148
Plot_PDCumulants_TPR.C:149
Plot_PDCumulants_TPR.C:150
Plot_PDCumulants_TPR.C:151
Plot_PDCumulants_TPR.C:152
Plot_PDCumulants_TPR.C:153
Plot_PDCumulants_TPR.C:154
Plot_PDCumulants_TPR.C:155
Plot_PDCumulants_TPR.C:156
Plot_PDCumulants_TPR.C:157
Plot_PDCumulants_TPR.C:158
Plot_PDCumulants_TPR.C:159
Plot_PDCumulants_TPR.C:160
Plot_PDCumulants_TPR.C:161
Plot_PDCumulants_TPR.C:162
Plot_PDCumulants_TPR.C:163
Plot_PDCumulants_TPR.C:164
Plot_PDCumulants_TPR.C:165
Plot_PDCumulants_TPR.C:166
Plot_PDCumulants_TPR.C:167
Plot_PDCumulants_TPR.C:168
Plot_PDCumulants_TPR.C:169
Plot_PDCumulants_TPR.C:170
Plot_PDCumulants_TPR.C:171
Plot_PDCumulants_TPR.C:172
Plot_PDCumulants_TPR.C:173
Plot_PDCumulants_TPR.C:174
Plot_PDCumulants_TPR.C:175
Plot_PDCumulants_TPR.C:176
Plot_PDCumulants_TPR.C:177
Plot_PDCumulants_TPR.C:178
Plot_PDCumulants_TPR.C:179
Plot_PDCumulants_TPR.C:180
Plot_PDCumulants_TPR.C:181
Plot_PDCumulants_TPR.C:182
Plot_PDCumulants_TPR.C:183
Plot_PDCumulants_TPR.C:184
Plot_PDCumulants_TPR.C:185
Plot_PDCumulants_TPR.C:186
Plot_PDCumulants_TPR.C:187
Plot_PDCumulants_TPR.C:188
Plot_PDCumulants_TPR.C:189
Plot_PDCumulants_TPR.C:190
Plot_PDCumulants_TPR.C:191
Plot_PDCumulants_TPR.C:192
Plot_PDCumulants_TPR.C:193
Plot_PDCumulants_TPR.C:194
Plot_PDCumulants_TPR.C:195
Plot_PDCumulants_TPR.C:196
Plot_PDCumulants_TPR.C:197
Plot_PDCumulants_TPR.C:198
Plot_PDCumulants_TPR.C:199
Plot_PDCumulants_TPR.C:200
Plot_PDCumulants_TPR.C:201
Plot_PDCumulants_TPR.C:202
Plot_PDCumulants_TPR.C:203
Plot_PDCumulants_TPR.C:204
Plot_PDCumulants_TPR.C:205
Plot_PDCumulants_TPR.C:206
Plot_PDCumulants_TPR.C:207
Plot_PDCumulants_TPR.C:208
Plot_PDCumulants_TPR.C:209
Plot_PDCumulants_TPR.C:210
Plot_PDCumulants_TPR.C:211
Plot_PDCumulants_TPR.C:212
Plot_PDCumulants_TPR.C:213
Plot_PDCumulants_TPR.C:214
Plot_PDCumulants_TPR.C:215
Plot_PDCumulants_TPR.C:216
Plot_PDCumulants_TPR.C:217
Plot_PDCumulants_TPR.C:218
Plot_PDCumulants_TPR.C:219
Plot_PDCumulants_TPR.C:220
Plot_PDCumulants_TPR.C:221
Plot_PDCumulants_TPR.C:222
Plot_PDCumulants_TPR.C:223
Plot_PDCumulants_TPR.C:224
Plot_PDCumulants_TPR.C:225
Plot_PDCumulants_TPR.C:226
Plot_PDCumulants_TPR.C:227
Plot_PDCumulants_TPR.C:228
Plot_PDCumulants_TPR.C:229
Plot_PDCumulants_TPR.C:230
Plot_PDCumulants_TPR.C:231
Plot_PDCumulants_TPR.C:232
Plot_PDCumulants_TPR.C:233
Plot_PDCumulants_TPR.C:234
Plot_PDCumulants_TPR.C:235
Plot_PDCumulants_TPR.C:236
Plot_PDCumulants_TPR.C:237
Plot_PDCumulants_TPR.C:238
Plot_PDCumulants_TPR.C:239
Plot_PDCumulants_TPR.C:240
Plot_PDCumulants_TPR.C:241
Plot_PDCumulants_TPR.C:242
Plot_PDCumulants_TPR.C:243
Plot_PDCumulants_TPR.C:244
Plot_PDCumulants_TPR.C:245
Plot_PDCumulants_TPR.C:246
Plot_PDCumulants_TPR.C:247
Plot_PDCumulants_TPR.C:248
Plot_PDCumulants_TPR.C:249
Plot_PDCumulants_TPR.C:250
Plot_PDCumulants_TPR.C:251
Plot_PDCumulants_TPR.C:252
Plot_PDCumulants_TPR.C:253
Plot_PDCumulants_TPR.C:254
Plot_PDCumulants_TPR.C:255
Plot_PDCumulants_TPR.C:256
Plot_PDCumulants_TPR.C:257
Plot_PDCumulants_TPR.C:258
Plot_PDCumulants_TPR.C:259
Plot_PDCumulants_TPR.C:260
Plot_PDCumulants_TPR.C:261
Plot_PDCumulants_TPR.C:262
Plot_PDCumulants_TPR.C:263
Plot_PDCumulants_TPR.C:264
Plot_PDCumulants_TPR.C:265
Plot_PDCumulants_TPR.C:266
Plot_PDCumulants_TPR.C:267
Plot_PDCumulants_TPR.C:268
Plot_PDCumulants_TPR.C:269
Plot_PDCumulants_TPR.C:270
Plot_PDCumulants_TPR.C:271
Plot_PDCumulants_TPR.C:272
Plot_PDCumulants_TPR.C:273
Plot_PDCumulants_TPR.C:274
Plot_PDCumulants_TPR.C:275
Plot_PDCumulants_TPR.C:276
Plot_PDCumulants_TPR.C:277
Plot_PDCumulants_TPR.C:278
Plot_PDCumulants_TPR.C:279
Plot_PDCumulants_TPR.C:280
Plot_PDCumulants_TPR.C:281
Plot_PDCumulants_TPR.C:282
Plot_PDCumulants_TPR.C:283
Plot_PDCumulants_TPR.C:284
Plot_PDCumulants_TPR.C:285
Plot_PDCumulants_TPR.C:286
Plot_PDCumulants_TPR.C:287
Plot_PDCumulants_TPR.C:288
Plot_PDCumulants_TPR.C:289
Plot_PDCumulants_TPR.C:290
Plot_PDCumulants_TPR.C:291
Plot_PDCumulants_TPR.C:292
Plot_PDCumulants_TPR.C:293
Plot_PDCumulants_TPR.C:294
Plot_PDCumulants_TPR.C:295
Plot_PDCumulants_TPR.C:296
Plot_PDCumulants_TPR.C:297
Plot_PDCumulants_TPR.C:298
Plot_PDCumulants_TPR.C:299
Plot_PDCumulants_TPR.C:300
Plot_PDCumulants_TPR.C:301
Plot_PDCumulants_TPR.C:302
Plot_PDCumulants_TPR.C:303
Plot_PDCumulants_TPR.C:304
Plot_PDCumulants_TPR.C:305
Plot_PDCumulants_TPR.C:306
Plot_PDCumulants_TPR.C:307
Plot_PDCumulants_TPR.C:308
Plot_PDCumulants_TPR.C:309
Plot_PDCumulants_TPR.C:310
Plot_PDCumulants_TPR.C:311
Plot_PDCumulants_TPR.C:312
Plot_PDCumulants_TPR.C:313
Plot_PDCumulants_TPR.C:314
Plot_PDCumulants_TPR.C:315
Plot_PDCumulants_TPR.C:316
Plot_PDCumulants_TPR.C:317
Plot_PDCumulants_TPR.C:318
Plot_PDCumulants_TPR.C:319
Plot_PDCumulants_TPR.C:320
Plot_PDCumulants_TPR.C:321
Plot_PDCumulants_TPR.C:322
Plot_PDCumulants_TPR.C:323
Plot_PDCumulants_TPR.C:324
Plot_PDCumulants_TPR.C:325
Plot_PDCumulants_TPR.C:326
Plot_PDCumulants_TPR.C:327
Plot_PDCumulants_TPR.C:328
Plot_PDCumulants_TPR.C:329
Plot_PDCumulants_TPR.C:330
Plot_PDCumulants_TPR.C:331
Plot_PDCumulants_TPR.C:332
Plot_PDCumulants_TPR.C:333
Plot_PDCumulants_TPR.C:334
Plot_PDCumulants_TPR.C:335
Plot_PDCumulants_TPR.C:336
Plot_PDCumulants_TPR.C:337
Plot_PDCumulants_TPR.C:338
Plot_PDCumulants_TPR.C:339
Plot_PDCumulants_TPR.C:340
Plot_PDCumulants_TPR.C:341
Plot_PDCumulants_TPR.C:342
Plot_PDCumulants_TPR.C:343
Plot_PDCumulants_TPR.C:344
Plot_PDCumulants_TPR.C:345
Plot_PDCumulants_TPR.C:346
Plot_PDCumulants_TPR.C:347
Plot_PDCumulants_TPR.C:348
Plot_PDCumulants_TPR.C:349
Plot_PDCumulants_TPR.C:350
Plot_PDCumulants_TPR.C:351
Plot_PDCumulants_TPR.C:352
Plot_PDCumulants_TPR.C:353
Plot_PDCumulants_TPR.C:354
Plot_PDCumulants_TPR.C:355
Plot_PDCumulants_TPR.C:356
Plot_PDCumulants_TPR.C:357
Plot_PDCumulants_TPR.C:358
Plot_PDCumulants_TPR.C:359
Plot_PDCumulants_TPR.C:360
Plot_PDCumulants_TPR.C:361
Plot_PDCumulants_TPR.C:362
Plot_PDCumulants_TPR.C:363
Plot_PDCumulants_TPR.C:364
Plot_PDCumulants_TPR.C:365
Plot_PDCumulants_TPR.C:366
Plot_PDCumulants_TPR.C:367
Plot_PDCumulants_TPR.C:368
Plot_PDCumulants_TPR.C:369
Plot_PDCumulants_TPR.C:370
Plot_PDCumulants_TPR.C:371
Plot_PDCumulants_TPR.C:372
Plot_PDCumulants_TPR.C:373
Plot_PDCumulants_TPR.C:374
Plot_PDCumulants_TPR.C:375
Plot_PDCumulants_TPR.C:376
Plot_PDCumulants_TPR.C:377
Plot_PDCumulants_TPR.C:378
Plot_PDCumulants_TPR.C:379
Plot_PDCumulants_TPR.C:380
Plot_PDCumulants_TPR.C:381
Plot_PDCumulants_TPR.C:382
Plot_PDCumulants_TPR.C:383
Plot_PDCumulants_TPR.C:384
Plot_PDCumulants_TPR.C:385
Plot_PDCumulants_TPR.C:386
Plot_PDCumulants_TPR.C:387
Plot_PDCumulants_TPR.C:388
Plot_PDCumulants_TPR.C:389
Plot_PDCumulants_TPR.C:390
Plot_PDCumulants_TPR.C:391
Plot_PDCumulants_TPR.C:392
Plot_PDCumulants_TPR.C:393
Plot_PDCumulants_TPR.C:394
Plot_PDCumulants_TPR.C:395
Plot_PDCumulants_TPR.C:396
Plot_PDCumulants_TPR.C:397
Plot_PDCumulants_TPR.C:398
Plot_PDCumulants_TPR.C:399
Plot_PDCumulants_TPR.C:400
Plot_PDCumulants_TPR.C:401
Plot_PDCumulants_TPR.C:402
Plot_PDCumulants_TPR.C:403
Plot_PDCumulants_TPR.C:404
Plot_PDCumulants_TPR.C:405
Plot_PDCumulants_TPR.C:406
Plot_PDCumulants_TPR.C:407
Plot_PDCumulants_TPR.C:408
Plot_PDCumulants_TPR.C:409
Plot_PDCumulants_TPR.C:410
Plot_PDCumulants_TPR.C:411
Plot_PDCumulants_TPR.C:412
Plot_PDCumulants_TPR.C:413
Plot_PDCumulants_TPR.C:414
Plot_PDCumulants_TPR.C:415
Plot_PDCumulants_TPR.C:416
Plot_PDCumulants_TPR.C:417
Plot_PDCumulants_TPR.C:418
Plot_PDCumulants_TPR.C:419
Plot_PDCumulants_TPR.C:420
Plot_PDCumulants_TPR.C:421
Plot_PDCumulants_TPR.C:422
Plot_PDCumulants_TPR.C:423
Plot_PDCumulants_TPR.C:424
Plot_PDCumulants_TPR.C:425
Plot_PDCumulants_TPR.C:426
Plot_PDCumulants_TPR.C:427
Plot_PDCumulants_TPR.C:428
Plot_PDCumulants_TPR.C:429
Plot_PDCumulants_TPR.C:430
Plot_PDCumulants_TPR.C:431
Plot_PDCumulants_TPR.C:432
Plot_PDCumulants_TPR.C:433
Plot_PDCumulants_TPR.C:434
Plot_PDCumulants_TPR.C:435
Plot_PDCumulants_TPR.C:436
Plot_PDCumulants_TPR.C:437
Plot_PDCumulants_TPR.C:438
Plot_PDCumulants_TPR.C:439
Plot_PDCumulants_TPR.C:440
Plot_PDCumulants_TPR.C:441
Plot_PDCumulants_TPR.C:442
Plot_PDCumulants_TPR.C:443
Plot_PDCumulants_TPR.C:444
Plot_PDCumulants_TPR.C:445
Plot_PDCumulants_TPR.C:446
Plot_PDCumulants_TPR.C:447
Plot_PDCumulants_TPR.C:448
Plot_PDCumulants_TPR.C:449
Plot_PDCumulants_TPR.C:450
Plot_PDCumulants_TPR.C:451
Plot_PDCumulants_TPR.C:452
Plot_PDCumulants_TPR.C:453
Plot_PDCumulants_TPR.C:454
Plot_PDCumulants_TPR.C:455
Plot_PDCumulants_TPR.C:456
Plot_PDCumulants_TPR.C:457
Plot_PDCumulants_TPR.C:458
Plot_PDCumulants_TPR.C:459
Plot_PDCumulants_TPR.C:460
Plot_PDCumulants_TPR.C:461
Plot_PDCumulants_TPR.C:462
Plot_PDCumulants_TPR.C:463
Plot_PDCumulants_TPR.C:464
Plot_PDCumulants_TPR.C:465
Plot_PDCumulants_TPR.C:466
Plot_PDCumulants_TPR.C:467
Plot_PDCumulants_TPR.C:468
Plot_PDCumulants_TPR.C:469
Plot_PDCumulants_TPR.C:470
Plot_PDCumulants_TPR.C:471
Plot_PDCumulants_TPR.C:472
Plot_PDCumulants_TPR.C:473
Plot_PDCumulants_TPR.C:474
Plot_PDCumulants_TPR.C:475
Plot_PDCumulants_TPR.C:476
Plot_PDCumulants_TPR.C:477
Plot_PDCumulants_TPR.C:478
Plot_PDCumulants_TPR.C:479
Plot_PDCumulants_TPR.C:480
Plot_PDCumulants_TPR.C:481
Plot_PDCumulants_TPR.C:482
Plot_PDCumulants_TPR.C:483
Plot_PDCumulants_TPR.C:484
Plot_PDCumulants_TPR.C:485
Plot_PDCumulants_TPR.C:486
Plot_PDCumulants_TPR.C:487
Plot_PDCumulants_TPR.C:488
Plot_PDCumulants_TPR.C:489
Plot_PDCumulants_TPR.C:490
Plot_PDCumulants_TPR.C:491
Plot_PDCumulants_TPR.C:492
Plot_PDCumulants_TPR.C:493
Plot_PDCumulants_TPR.C:494
Plot_PDCumulants_TPR.C:495
Plot_PDCumulants_TPR.C:496
Plot_PDCumulants_TPR.C:497
Plot_PDCumulants_TPR.C:498
Plot_PDCumulants_TPR.C:499
Plot_PDCumulants_TPR.C:500
Plot_PDCumulants_TPR.C:501
Plot_PDCumulants_TPR.C:502
Plot_PDCumulants_TPR.C:503
Plot_PDCumulants_TPR.C:504
Plot_PDCumulants_TPR.C:505
Plot_PDCumulants_TPR.C:506
Plot_PDCumulants_TPR.C:507
Plot_PDCumulants_TPR.C:508
Plot_PDCumulants_TPR.C:509
Plot_PDCumulants_TPR.C:510
Plot_PDCumulants_TPR.C:511
Plot_PDCumulants_TPR.C:512
Plot_PDCumulants_TPR.C:513
Plot_PDCumulants_TPR.C:514
Plot_PDCumulants_TPR.C:515
Plot_PDCumulants_TPR.C:516
Plot_PDCumulants_TPR.C:517
Plot_PDCumulants_TPR.C:518
Plot_PDCumulants_TPR.C:519
Plot_PDCumulants_TPR.C:520
Plot_PDCumulants_TPR.C:521
Plot_PDCumulants_TPR.C:522
Plot_PDCumulants_TPR.C:523
Plot_PDCumulants_TPR.C:524
Plot_PDCumulants_TPR.C:525
Plot_PDCumulants_TPR.C:526
Plot_PDCumulants_TPR.C:527
Plot_PDCumulants_TPR.C:528
Plot_PDCumulants_TPR.C:529
Plot_PDCumulants_TPR.C:530
Plot_PDCumulants_TPR.C:531
Plot_PDCumulants_TPR.C:532
Plot_PDCumulants_TPR.C:533
Plot_PDCumulants_TPR.C:534
Plot_PDCumulants_TPR.C:535
Plot_PDCumulants_TPR.C:536
Plot_PDCumulants_TPR.C:537
Plot_PDCumulants_TPR.C:538
Plot_PDCumulants_TPR.C:539
Plot_PDCumulants_TPR.C:540
Plot_PDCumulants_TPR.C:541
Plot_PDCumulants_TPR.C:542
Plot_PDCumulants_TPR.C:543
Plot_PDCumulants_TPR.C:544
Plot_PDCumulants_TPR.C:545
Plot_PDCumulants_TPR.C:546
Plot_PDCumulants_TPR.C:547
Plot_PDCumulants_TPR.C:548
Plot_PDCumulants_TPR.C:549
Plot_PDCumulants_TPR.C:550
Plot_PDCumulants_TPR.C:551
Plot_PDCumulants_TPR.C:552
Plot_PDCumulants_TPR.C:553
Plot_PDCumulants_TPR.C:554
Plot_PDCumulants_TPR.C:555
Plot_PDCumulants_TPR.C:556
Plot_PDCumulants_TPR.C:557
Plot_PDCumulants_TPR.C:558
Plot_PDCumulants_TPR.C:559
Plot_PDCumulants_TPR.C:560
Plot_PDCumulants_TPR.C:561
Plot_PDCumulants_TPR.C:562
Plot_PDCumulants_TPR.C:563
Plot_PDCumulants_TPR.C:564
Plot_PDCumulants_TPR.C:565
Plot_PDCumulants_TPR.C:566
Plot_PDCumulants_TPR.C:567
Plot_PDCumulants_TPR.C:568
Plot_PDCumulants_TPR.C:569
Plot_PDCumulants_TPR.C:570
Plot_PDCumulants_TPR.C:571
Plot_PDCumulants_TPR.C:572
Plot_PDCumulants_TPR.C:573
Plot_PDCumulants_TPR.C:574
Plot_PDCumulants_TPR.C:575
Plot_PDCumulants_TPR.C:576
Plot_PDCumulants_TPR.C:577
Plot_PDCumulants_TPR.C:578
Plot_PDCumulants_TPR.C:579
Plot_PDCumulants_TPR.C:580
Plot_PDCumulants_TPR.C:581
Plot_PDCumulants_TPR.C:582
Plot_PDCumulants_TPR.C:583
Plot_PDCumulants_TPR.C:584
Plot_PDCumulants_TPR.C:585
Plot_PDCumulants_TPR.C:586
Plot_PDCumulants_TPR.C:587
Plot_PDCumulants_TPR.C:588
Plot_PDCumulants_TPR.C:589
Plot_PDCumulants_TPR.C:590
Plot_PDCumulants_TPR.C:591
Plot_PDCumulants_TPR.C:592
Plot_PDCumulants_TPR.C:593
Plot_PDCumulants_TPR.C:594
Plot_PDCumulants_TPR.C:595
Plot_PDCumulants_TPR.C:596
Plot_PDCumulants_TPR.C:597
Plot_PDCumulants_TPR.C:598
Plot_PDCumulants_TPR.C:599
Plot_PDCumulants_TPR.C:600
Plot_PDCumulants_TPR.C:601
Plot_PDCumulants_TPR.C:602
Plot_PDCumulants_TPR.C:603
Plot_PDCumulants_TPR.C:604
Plot_PDCumulants_TPR.C:605
Plot_PDCumulants_TPR.C:606
Plot_PDCumulants_TPR.C:607
Plot_PDCumulants_TPR.C:608
Plot_PDCumulants_TPR.C:609
Plot_PDCumulants_TPR.C:610
Plot_PDCumulants_TPR.C:611
Plot_PDCumulants_TPR.C:612
Plot_PDCumulants_TPR.C:613
Plot_PDCumulants_TPR.C:614
Plot_PDCumulants_TPR.C:615
Plot_PDCumulants_TPR.C:616
Plot_PDCumulants_TPR.C:617
Plot_PDCumulants_TPR.C:618
Plot_PDCumulants_TPR.C:619
Plot_PDCumulants_TPR.C:620
Plot_PDCumulants_TPR.C:621
Plot_PDCumulants_TPR.C:622
Plot_PDCumulants_TPR.C:623
Plot_PDCumulants_TPR.C:624
Plot_PDCumulants_TPR.C:625
Plot_PDCumulants_TPR.C:626
Plot_PDCumulants_TPR.C:627
Plot_PDCumulants_TPR.C:628
Plot_PDCumulants_TPR.C:629
Plot_PDCumulants_TPR.C:630
Plot_PDCumulants_TPR.C:631
Plot_PDCumulants_TPR.C:632
Plot_PDCumulants_TPR.C:633
Plot_PDCumulants_TPR.C:634
Plot_PDCumulants_TPR.C:635
Plot_PDCumulants_TPR.C:636
Plot_PDCumulants_TPR.C:637
Plot_PDCumulants_TPR.C:638
Plot_PDCumulants_TPR.C:639
Plot_PDCumulants_TPR.C:640
Plot_PDCumulants_TPR.C:641
Plot_PDCumulants_TPR.C:642
Plot_PDCumulants_TPR.C:643
Plot_PDCumulants_TPR.C:644
Plot_PDCumulants_TPR.C:645
Plot_PDCumulants_TPR.C:646
Plot_PDCumulants_TPR.C:647
Plot_PDCumulants_TPR.C:648
Plot_PDCumulants_TPR.C:649
Plot_PDCumulants_TPR.C:650
Plot_PDCumulants_TPR.C:651
Plot_PDCumulants_TPR.C:652
Plot_PDCumulants_TPR.C:653
Plot_PDCumulants_TPR.C:654
Plot_PDCumulants_TPR.C:655
Plot_PDCumulants_TPR.C:656
Plot_PDCumulants_TPR.C:657
Plot_PDCumulants_TPR.C:658
Plot_PDCumulants_TPR.C:659
Plot_PDCumulants_TPR.C:660
Plot_PDCumulants_TPR.C:661
Plot_PDCumulants_TPR.C:662
Plot_PDCumulants_TPR.C:663
Plot_PDCumulants_TPR.C:664
Plot_PDCumulants_TPR.C:665
Plot_PDCumulants_TPR.C:666
Plot_PDCumulants_TPR.C:667
Plot_PDCumulants_TPR.C:668
Plot_PDCumulants_TPR.C:669
Plot_PDCumulants_TPR.C:670
Plot_PDCumulants_TPR.C:671
Plot_PDCumulants_TPR.C:672
Plot_PDCumulants_TPR.C:673
Plot_PDCumulants_TPR.C:674
Plot_PDCumulants_TPR.C:675
Plot_PDCumulants_TPR.C:676
Plot_PDCumulants_TPR.C:677
Plot_PDCumulants_TPR.C:678
Plot_PDCumulants_TPR.C:679
Plot_PDCumulants_TPR.C:680
Plot_PDCumulants_TPR.C:681
Plot_PDCumulants_TPR.C:682
Plot_PDCumulants_TPR.C:683
Plot_PDCumulants_TPR.C:684
Plot_PDCumulants_TPR.C:685
Plot_PDCumulants_TPR.C:686
Plot_PDCumulants_TPR.C:687
Plot_PDCumulants_TPR.C:688
Plot_PDCumulants_TPR.C:689
Plot_PDCumulants_TPR.C:690
Plot_PDCumulants_TPR.C:691
Plot_PDCumulants_TPR.C:692
Plot_PDCumulants_TPR.C:693
Plot_PDCumulants_TPR.C:694
Plot_PDCumulants_TPR.C:695
Plot_PDCumulants_TPR.C:696
Plot_PDCumulants_TPR.C:697
Plot_PDCumulants_TPR.C:698
Plot_PDCumulants_TPR.C:699
Plot_PDCumulants_TPR.C:700
Plot_PDCumulants_TPR.C:701
Plot_PDCumulants_TPR.C:702
Plot_PDCumulants_TPR.C:703
Plot_PDCumulants_TPR.C:704
Plot_PDCumulants_TPR.C:705
Plot_PDCumulants_TPR.C:706
Plot_PDCumulants_TPR.C:707
Plot_PDCumulants_TPR.C:708
Plot_PDCumulants_TPR.C:709
Plot_PDCumulants_TPR.C:710
Plot_PDCumulants_TPR.C:711
Plot_PDCumulants_TPR.C:712
Plot_PDCumulants_TPR.C:713
Plot_PDCumulants_TPR.C:714
Plot_PDCumulants_TPR.C:715
Plot_PDCumulants_TPR.C:716
Plot_PDCumulants_TPR.C:717
Plot_PDCumulants_TPR.C:718
Plot_PDCumulants_TPR.C:719
Plot_PDCumulants_TPR.C:720
Plot_PDCumulants_TPR.C:721
Plot_PDCumulants_TPR.C:722
Plot_PDCumulants_TPR.C:723
Plot_PDCumulants_TPR.C:724
Plot_PDCumulants_TPR.C:725
Plot_PDCumulants_TPR.C:726
Plot_PDCumulants_TPR.C:727
Plot_PDCumulants_TPR.C:728
Plot_PDCumulants_TPR.C:729
Plot_PDCumulants_TPR.C:730
Plot_PDCumulants_TPR.C:731
Plot_PDCumulants_TPR.C:732
Plot_PDCumulants_TPR.C:733
Plot_PDCumulants_TPR.C:734
Plot_PDCumulants_TPR.C:735
Plot_PDCumulants_TPR.C:736
Plot_PDCumulants_TPR.C:737
Plot_PDCumulants_TPR.C:738
Plot_PDCumulants_TPR.C:739
Plot_PDCumulants_TPR.C:740
Plot_PDCumulants_TPR.C:741
Plot_PDCumulants_TPR.C:742
Plot_PDCumulants_TPR.C:743
Plot_PDCumulants_TPR.C:744
Plot_PDCumulants_TPR.C:745
Plot_PDCumulants_TPR.C:746
Plot_PDCumulants_TPR.C:747
Plot_PDCumulants_TPR.C:748
Plot_PDCumulants_TPR.C:749
Plot_PDCumulants_TPR.C:750
Plot_PDCumulants_TPR.C:751
Plot_PDCumulants_TPR.C:752
Plot_PDCumulants_TPR.C:753
Plot_PDCumulants_TPR.C:754
Plot_PDCumulants_TPR.C:755
Plot_PDCumulants_TPR.C:756
Plot_PDCumulants_TPR.C:757
Plot_PDCumulants_TPR.C:758
Plot_PDCumulants_TPR.C:759
Plot_PDCumulants_TPR.C:760
Plot_PDCumulants_TPR.C:761
Plot_PDCumulants_TPR.C:762
Plot_PDCumulants_TPR.C:763
Plot_PDCumulants_TPR.C:764
Plot_PDCumulants_TPR.C:765
Plot_PDCumulants_TPR.C:766
Plot_PDCumulants_TPR.C:767
Plot_PDCumulants_TPR.C:768
Plot_PDCumulants_TPR.C:769
Plot_PDCumulants_TPR.C:770
Plot_PDCumulants_TPR.C:771
Plot_PDCumulants_TPR.C:772
Plot_PDCumulants_TPR.C:773
Plot_PDCumulants_TPR.C:774
Plot_PDCumulants_TPR.C:775
Plot_PDCumulants_TPR.C:776
Plot_PDCumulants_TPR.C:777
Plot_PDCumulants_TPR.C:778
Plot_PDCumulants_TPR.C:779
Plot_PDCumulants_TPR.C:780
Plot_PDCumulants_TPR.C:781
Plot_PDCumulants_TPR.C:782
Plot_PDCumulants_TPR.C:783
Plot_PDCumulants_TPR.C:784
Plot_PDCumulants_TPR.C:785
Plot_PDCumulants_TPR.C:786
Plot_PDCumulants_TPR.C:787
Plot_PDCumulants_TPR.C:788
Plot_PDCumulants_TPR.C:789
Plot_PDCumulants_TPR.C:790
Plot_PDCumulants_TPR.C:791
Plot_PDCumulants_TPR.C:792
Plot_PDCumulants_TPR.C:793
Plot_PDCumulants_TPR.C:794
Plot_PDCumulants_TPR.C:795
Plot_PDCumulants_TPR.C:796
Plot_PDCumulants_TPR.C:797
Plot_PDCumulants_TPR.C:798
Plot_PDCumulants_TPR.C:799
Plot_PDCumulants_TPR.C:800
Plot_PDCumulants_TPR.C:801
Plot_PDCumulants_TPR.C:802
Plot_PDCumulants_TPR.C:803
Plot_PDCumulants_TPR.C:804
Plot_PDCumulants_TPR.C:805
Plot_PDCumulants_TPR.C:806
Plot_PDCumulants_TPR.C:807
Plot_PDCumulants_TPR.C:808
Plot_PDCumulants_TPR.C:809
Plot_PDCumulants_TPR.C:810
Plot_PDCumulants_TPR.C:811
Plot_PDCumulants_TPR.C:812
Plot_PDCumulants_TPR.C:813
Plot_PDCumulants_TPR.C:814
Plot_PDCumulants_TPR.C:815
Plot_PDCumulants_TPR.C:816
Plot_PDCumulants_TPR.C:817
Plot_PDCumulants_TPR.C:818
Plot_PDCumulants_TPR.C:819
Plot_PDCumulants_TPR.C:820
Plot_PDCumulants_TPR.C:821
Plot_PDCumulants_TPR.C:822
Plot_PDCumulants_TPR.C:823
Plot_PDCumulants_TPR.C:824
Plot_PDCumulants_TPR.C:825
Plot_PDCumulants_TPR.C:826
Plot_PDCumulants_TPR.C:827
Plot_PDCumulants_TPR.C:828
Plot_PDCumulants_TPR.C:829
Plot_PDCumulants_TPR.C:830
Plot_PDCumulants_TPR.C:831
Plot_PDCumulants_TPR.C:832
Plot_PDCumulants_TPR.C:833
Plot_PDCumulants_TPR.C:834
Plot_PDCumulants_TPR.C:835
Plot_PDCumulants_TPR.C:836
Plot_PDCumulants_TPR.C:837
Plot_PDCumulants_TPR.C:838
Plot_PDCumulants_TPR.C:839
Plot_PDCumulants_TPR.C:840
Plot_PDCumulants_TPR.C:841
Plot_PDCumulants_TPR.C:842
Plot_PDCumulants_TPR.C:843
Plot_PDCumulants_TPR.C:844
Plot_PDCumulants_TPR.C:845
Plot_PDCumulants_TPR.C:846
Plot_PDCumulants_TPR.C:847
Plot_PDCumulants_TPR.C:848
Plot_PDCumulants_TPR.C:849
Plot_PDCumulants_TPR.C:850
Plot_PDCumulants_TPR.C:851
Plot_PDCumulants_TPR.C:852
Plot_PDCumulants_TPR.C:853
Plot_PDCumulants_TPR.C:854
Plot_PDCumulants_TPR.C:855
Plot_PDCumulants_TPR.C:856
Plot_PDCumulants_TPR.C:857
Plot_PDCumulants_TPR.C:858
Plot_PDCumulants_TPR.C:859
Plot_PDCumulants_TPR.C:860
Plot_PDCumulants_TPR.C:861
Plot_PDCumulants_TPR.C:862
Plot_PDCumulants_TPR.C:863
Plot_PDCumulants_TPR.C:864
Plot_PDCumulants_TPR.C:865
Plot_PDCumulants_TPR.C:866
Plot_PDCumulants_TPR.C:867
Plot_PDCumulants_TPR.C:868
Plot_PDCumulants_TPR.C:869
Plot_PDCumulants_TPR.C:870
Plot_PDCumulants_TPR.C:871
Plot_PDCumulants_TPR.C:872
Plot_PDCumulants_TPR.C:873
Plot_PDCumulants_TPR.C:874
Plot_PDCumulants_TPR.C:875
Plot_PDCumulants_TPR.C:876
Plot_PDCumulants_TPR.C:877
Plot_PDCumulants_TPR.C:878
Plot_PDCumulants_TPR.C:879
Plot_PDCumulants_TPR.C:880
Plot_PDCumulants_TPR.C:881
Plot_PDCumulants_TPR.C:882
Plot_PDCumulants_TPR.C:883
Plot_PDCumulants_TPR.C:884
Plot_PDCumulants_TPR.C:885
Plot_PDCumulants_TPR.C:886
Plot_PDCumulants_TPR.C:887
Plot_PDCumulants_TPR.C:888
Plot_PDCumulants_TPR.C:889
Plot_PDCumulants_TPR.C:890
Plot_PDCumulants_TPR.C:891
Plot_PDCumulants_TPR.C:892
Plot_PDCumulants_TPR.C:893
Plot_PDCumulants_TPR.C:894
Plot_PDCumulants_TPR.C:895
Plot_PDCumulants_TPR.C:896
Plot_PDCumulants_TPR.C:897
Plot_PDCumulants_TPR.C:898
Plot_PDCumulants_TPR.C:899
Plot_PDCumulants_TPR.C:900
Plot_PDCumulants_TPR.C:901
Plot_PDCumulants_TPR.C:902
Plot_PDCumulants_TPR.C:903
Plot_PDCumulants_TPR.C:904
Plot_PDCumulants_TPR.C:905
Plot_PDCumulants_TPR.C:906
Plot_PDCumulants_TPR.C:907
Plot_PDCumulants_TPR.C:908
Plot_PDCumulants_TPR.C:909
Plot_PDCumulants_TPR.C:910
Plot_PDCumulants_TPR.C:911
Plot_PDCumulants_TPR.C:912
Plot_PDCumulants_TPR.C:913
Plot_PDCumulants_TPR.C:914
Plot_PDCumulants_TPR.C:915
Plot_PDCumulants_TPR.C:916
Plot_PDCumulants_TPR.C:917
Plot_PDCumulants_TPR.C:918
Plot_PDCumulants_TPR.C:919
Plot_PDCumulants_TPR.C:920
Plot_PDCumulants_TPR.C:921
Plot_PDCumulants_TPR.C:922
Plot_PDCumulants_TPR.C:923
Plot_PDCumulants_TPR.C:924
Plot_PDCumulants_TPR.C:925
Plot_PDCumulants_TPR.C:926
Plot_PDCumulants_TPR.C:927
Plot_PDCumulants_TPR.C:928
Plot_PDCumulants_TPR.C:929
Plot_PDCumulants_TPR.C:930
Plot_PDCumulants_TPR.C:931
Plot_PDCumulants_TPR.C:932
Plot_PDCumulants_TPR.C:933
Plot_PDCumulants_TPR.C:934
Plot_PDCumulants_TPR.C:935
Plot_PDCumulants_TPR.C:936
Plot_PDCumulants_TPR.C:937
Plot_PDCumulants_TPR.C:938
Plot_PDCumulants_TPR.C:939
Plot_PDCumulants_TPR.C:940
Plot_PDCumulants_TPR.C:941
Plot_PDCumulants_TPR.C:942
Plot_PDCumulants_TPR.C:943
Plot_PDCumulants_TPR.C:944
Plot_PDCumulants_TPR.C:945
Plot_PDCumulants_TPR.C:946
Plot_PDCumulants_TPR.C:947
Plot_PDCumulants_TPR.C:948
Plot_PDCumulants_TPR.C:949
Plot_PDCumulants_TPR.C:950
Plot_PDCumulants_TPR.C:951
Plot_PDCumulants_TPR.C:952
Plot_PDCumulants_TPR.C:953
Plot_PDCumulants_TPR.C:954
Plot_PDCumulants_TPR.C:955
Plot_PDCumulants_TPR.C:956
Plot_PDCumulants_TPR.C:957
Plot_PDCumulants_TPR.C:958
Plot_PDCumulants_TPR.C:959
Plot_PDCumulants_TPR.C:960
Plot_PDCumulants_TPR.C:961
Plot_PDCumulants_TPR.C:962
Plot_PDCumulants_TPR.C:963
Plot_PDCumulants_TPR.C:964
Plot_PDCumulants_TPR.C:965
Plot_PDCumulants_TPR.C:966
Plot_PDCumulants_TPR.C:967
Plot_PDCumulants_TPR.C:968
Plot_PDCumulants_TPR.C:969
Plot_PDCumulants_TPR.C:970
Plot_PDCumulants_TPR.C:971
Plot_PDCumulants_TPR.C:972
Plot_PDCumulants_TPR.C:973
Plot_PDCumulants_TPR.C:974
Plot_PDCumulants_TPR.C:975
Plot_PDCumulants_TPR.C:976
Plot_PDCumulants_TPR.C:977
Plot_PDCumulants_TPR.C:978
Plot_PDCumulants_TPR.C:979
Plot_PDCumulants_TPR.C:980
Plot_PDCumulants_TPR.C:981
Plot_PDCumulants_TPR.C:982
Plot_PDCumulants_TPR.C:983
Plot_PDCumulants_TPR.C:984
Plot_PDCumulants_TPR.C:985
Plot_PDCumulants_TPR.C:986
Plot_PDCumulants_TPR.C:987
Plot_PDCumulants_TPR.C:988
Plot_PDCumulants_TPR.C:989
Plot_PDCumulants_TPR.C:990
Plot_PDCumulants_TPR.C:991
Plot_PDCumulants_TPR.C:992
Plot_PDCumulants_TPR.C:993
Plot_PDCumulants_TPR.C:994
Plot_PDCumulants_TPR.C:995
Plot_PDCumulants_TPR.C:996
Plot_PDCumulants_TPR.C:997
Plot_PDCumulants_TPR.C:998
Plot_PDCumulants_TPR.C:999
Plot_PDCumulants_TPR.C:1000
Plot_PDCumulants_TPR.C:1001
Plot_PDCumulants_TPR.C:1002
Plot_PDCumulants_TPR.C:1003
Plot_PDCumulants_TPR.C:1004
Plot_PDCumulants_TPR.C:1005
Plot_PDCumulants_TPR.C:1006
Plot_PDCumulants_TPR.C:1007
Plot_PDCumulants_TPR.C:1008
Plot_PDCumulants_TPR.C:1009
Plot_PDCumulants_TPR.C:1010
Plot_PDCumulants_TPR.C:1011
Plot_PDCumulants_TPR.C:1012
Plot_PDCumulants_TPR.C:1013
Plot_PDCumulants_TPR.C:1014
Plot_PDCumulants_TPR.C:1015
Plot_PDCumulants_TPR.C:1016
Plot_PDCumulants_TPR.C:1017
Plot_PDCumulants_TPR.C:1018
Plot_PDCumulants_TPR.C:1019
Plot_PDCumulants_TPR.C:1020
Plot_PDCumulants_TPR.C:1021
Plot_PDCumulants_TPR.C:1022
Plot_PDCumulants_TPR.C:1023
Plot_PDCumulants_TPR.C:1024
Plot_PDCumulants_TPR.C:1025
Plot_PDCumulants_TPR.C:1026
Plot_PDCumulants_TPR.C:1027
Plot_PDCumulants_TPR.C:1028
Plot_PDCumulants_TPR.C:1029
Plot_PDCumulants_TPR.C:1030
Plot_PDCumulants_TPR.C:1031
Plot_PDCumulants_TPR.C:1032
Plot_PDCumulants_TPR.C:1033
Plot_PDCumulants_TPR.C:1034
Plot_PDCumulants_TPR.C:1035
Plot_PDCumulants_TPR.C:1036
Plot_PDCumulants_TPR.C:1037
Plot_PDCumulants_TPR.C:1038
Plot_PDCumulants_TPR.C:1039
Plot_PDCumulants_TPR.C:1040
Plot_PDCumulants_TPR.C:1041
Plot_PDCumulants_TPR.C:1042
Plot_PDCumulants_TPR.C:1043
Plot_PDCumulants_TPR.C:1044
Plot_PDCumulants_TPR.C:1045
Plot_PDCumulants_TPR.C:1046
Plot_PDCumulants_TPR.C:1047
Plot_PDCumulants_TPR.C:1048
Plot_PDCumulants_TPR.C:1049
Plot_PDCumulants_TPR.C:1050
Plot_PDCumulants_TPR.C:1051
Plot_PDCumulants_TPR.C:1052
Plot_PDCumulants_TPR.C:1053
Plot_PDCumulants_TPR.C:1054
Plot_PDCumulants_TPR.C:1055
Plot_PDCumulants_TPR.C:1056
Plot_PDCumulants_TPR.C:1057
Plot_PDCumulants_TPR.C:1058
Plot_PDCumulants_TPR.C:1059
Plot_PDCumulants_TPR.C:1060
Plot_PDCumulants_TPR.C:1061
Plot_PDCumulants_TPR.C:1062
Plot_PDCumulants_TPR.C:1063
Plot_PDCumulants_TPR.C:1064
Plot_PDCumulants_TPR.C:1065
Plot_PDCumulants_TPR.C:1066
Plot_PDCumulants_TPR.C:1067
Plot_PDCumulants_TPR.C:1068
Plot_PDCumulants_TPR.C:1069
Plot_PDCumulants_TPR.C:1070
Plot_PDCumulants_TPR.C:1071
Plot_PDCumulants_TPR.C:1072
Plot_PDCumulants_TPR.C:1073
Plot_PDCumulants_TPR.C:1074
Plot_PDCumulants_TPR.C:1075
Plot_PDCumulants_TPR.C:1076
Plot_PDCumulants_TPR.C:1077
Plot_PDCumulants_TPR.C:1078
Plot_PDCumulants_TPR.C:1079
Plot_PDCumulants_TPR.C:1080
Plot_PDCumulants_TPR.C:1081
Plot_PDCumulants_TPR.C:1082
Plot_PDCumulants_TPR.C:1083
Plot_PDCumulants_TPR.C:1084
Plot_PDCumulants_TPR.C:1085
Plot_PDCumulants_TPR.C:1086
Plot_PDCumulants_TPR.C:1087
Plot_PDCumulants_TPR.C:1088
Plot_PDCumulants_TPR.C:1089
Plot_PDCumulants_TPR.C:1090
Plot_PDCumulants_TPR.C:1091
Plot_PDCumulants_TPR.C:1092
Plot_PDCumulants_TPR.C:1093
Plot_PDCumulants_TPR.C:1094
Plot_PDCumulants_TPR.C:1095
Plot_PDCumulants_TPR.C:1096
Plot_PDCumulants_TPR.C:1097
Plot_PDCumulants_TPR.C:1098
Plot_PDCumulants_TPR.C:1099
Plot_PDCumulants_TPR.C:1100
Plot_PDCumulants_TPR.C:1101
Plot_PDCumulants_TPR.C:1102
Plot_PDCumulants_TPR.C:1103
Plot_PDCumulants_TPR.C:1104
Plot_PDCumulants_TPR.C:1105
Plot_PDCumulants_TPR.C:1106
Plot_PDCumulants_TPR.C:1107
Plot_PDCumulants_TPR.C:1108
Plot_PDCumulants_TPR.C:1109
Plot_PDCumulants_TPR.C:1110
Plot_PDCumulants_TPR.C:1111
Plot_PDCumulants_TPR.C:1112
Plot_PDCumulants_TPR.C:1113
Plot_PDCumulants_TPR.C:1114
Plot_PDCumulants_TPR.C:1115
Plot_PDCumulants_TPR.C:1116
Plot_PDCumulants_TPR.C:1117
Plot_PDCumulants_TPR.C:1118
Plot_PDCumulants_TPR.C:1119
Plot_PDCumulants_TPR.C:1120
Plot_PDCumulants_TPR.C:1121
Plot_PDCumulants_TPR.C:1122
Plot_PDCumulants_TPR.C:1123
Plot_PDCumulants_TPR.C:1124
Plot_PDCumulants_TPR.C:1125
Plot_PDCumulants_TPR.C:1126
Plot_PDCumulants_TPR.C:1127
Plot_PDCumulants_TPR.C:1128
Plot_PDCumulants_TPR.C:1129
Plot_PDCumulants_TPR.C:1130
Plot_PDCumulants_TPR.C:1131
Plot_PDCumulants_TPR.C:1132
Plot_PDCumulants_TPR.C:1133
Plot_PDCumulants_TPR.C:1134
Plot_PDCumulants_TPR.C:1135
Plot_PDCumulants_TPR.C:1136
Plot_PDCumulants_TPR.C:1137
Plot_PDCumulants_TPR.C:1138
Plot_PDCumulants_TPR.C:1139
Plot_PDCumulants_TPR.C:1140
Plot_PDCumulants_TPR.C:1141
Plot_PDCumulants_TPR.C:1142
Plot_PDCumulants_TPR.C:1143
Plot_PDCumulants_TPR.C:1144
Plot_PDCumulants_TPR.C:1145
Plot_PDCumulants_TPR.C:1146
Plot_PDCumulants_TPR.C:1147
Plot_PDCumulants_TPR.C:1148
Plot_PDCumulants_TPR.C:1149
Plot_PDCumulants_TPR.C:1150
Plot_PDCumulants_TPR.C:1151
Plot_PDCumulants_TPR.C:1152
Plot_PDCumulants_TPR.C:1153
Plot_PDCumulants_TPR.C:1154
Plot_PDCumulants_TPR.C:1155
Plot_PDCumulants_TPR.C:1156
Plot_PDCumulants_TPR.C:1157
Plot_PDCumulants_TPR.C:1158
Plot_PDCumulants_TPR.C:1159
Plot_PDCumulants_TPR.C:1160
Plot_PDCumulants_TPR.C:1161
Plot_PDCumulants_TPR.C:1162
Plot_PDCumulants_TPR.C:1163
Plot_PDCumulants_TPR.C:1164
Plot_PDCumulants_TPR.C:1165
Plot_PDCumulants_TPR.C:1166
Plot_PDCumulants_TPR.C:1167
Plot_PDCumulants_TPR.C:1168
Plot_PDCumulants_TPR.C:1169
Plot_PDCumulants_TPR.C:1170
Plot_PDCumulants_TPR.C:1171
Plot_PDCumulants_TPR.C:1172
Plot_PDCumulants_TPR.C:1173
Plot_PDCumulants_TPR.C:1174
Plot_PDCumulants_TPR.C:1175
Plot_PDCumulants_TPR.C:1176
Plot_PDCumulants_TPR.C:1177
Plot_PDCumulants_TPR.C:1178
Plot_PDCumulants_TPR.C:1179
Plot_PDCumulants_TPR.C:1180
Plot_PDCumulants_TPR.C:1181
Plot_PDCumulants_TPR.C:1182
Plot_PDCumulants_TPR.C:1183
Plot_PDCumulants_TPR.C:1184
Plot_PDCumulants_TPR.C:1185
Plot_PDCumulants_TPR.C:1186
Plot_PDCumulants_TPR.C:1187
Plot_PDCumulants_TPR.C:1188
Plot_PDCumulants_TPR.C:1189
Plot_PDCumulants_TPR.C:1190
Plot_PDCumulants_TPR.C:1191
Plot_PDCumulants_TPR.C:1192
Plot_PDCumulants_TPR.C:1193
Plot_PDCumulants_TPR.C:1194
Plot_PDCumulants_TPR.C:1195
Plot_PDCumulants_TPR.C:1196
Plot_PDCumulants_TPR.C:1197
Plot_PDCumulants_TPR.C:1198
Plot_PDCumulants_TPR.C:1199
Plot_PDCumulants_TPR.C:1200
Plot_PDCumulants_TPR.C:1201
Plot_PDCumulants_TPR.C:1202
Plot_PDCumulants_TPR.C:1203
Plot_PDCumulants_TPR.C:1204
Plot_PDCumulants_TPR.C:1205
Plot_PDCumulants_TPR.C:1206
Plot_PDCumulants_TPR.C:1207
Plot_PDCumulants_TPR.C:1208
Plot_PDCumulants_TPR.C:1209
Plot_PDCumulants_TPR.C:1210
Plot_PDCumulants_TPR.C:1211
Plot_PDCumulants_TPR.C:1212
Plot_PDCumulants_TPR.C:1213
Plot_PDCumulants_TPR.C:1214
Plot_PDCumulants_TPR.C:1215
Plot_PDCumulants_TPR.C:1216
Plot_PDCumulants_TPR.C:1217
Plot_PDCumulants_TPR.C:1218
Plot_PDCumulants_TPR.C:1219
Plot_PDCumulants_TPR.C:1220
Plot_PDCumulants_TPR.C:1221
Plot_PDCumulants_TPR.C:1222
Plot_PDCumulants_TPR.C:1223
Plot_PDCumulants_TPR.C:1224
Plot_PDCumulants_TPR.C:1225
Plot_PDCumulants_TPR.C:1226
Plot_PDCumulants_TPR.C:1227
Plot_PDCumulants_TPR.C:1228
Plot_PDCumulants_TPR.C:1229
Plot_PDCumulants_TPR.C:1230
Plot_PDCumulants_TPR.C:1231
Plot_PDCumulants_TPR.C:1232
Plot_PDCumulants_TPR.C:1233
Plot_PDCumulants_TPR.C:1234
Plot_PDCumulants_TPR.C:1235
Plot_PDCumulants_TPR.C:1236
Plot_PDCumulants_TPR.C:1237
Plot_PDCumulants_TPR.C:1238
Plot_PDCumulants_TPR.C:1239
Plot_PDCumulants_TPR.C:1240
Plot_PDCumulants_TPR.C:1241
Plot_PDCumulants_TPR.C:1242
Plot_PDCumulants_TPR.C:1243
Plot_PDCumulants_TPR.C:1244
Plot_PDCumulants_TPR.C:1245
Plot_PDCumulants_TPR.C:1246
Plot_PDCumulants_TPR.C:1247
Plot_PDCumulants_TPR.C:1248
Plot_PDCumulants_TPR.C:1249
Plot_PDCumulants_TPR.C:1250
Plot_PDCumulants_TPR.C:1251
Plot_PDCumulants_TPR.C:1252
Plot_PDCumulants_TPR.C:1253
Plot_PDCumulants_TPR.C:1254
Plot_PDCumulants_TPR.C:1255
Plot_PDCumulants_TPR.C:1256
Plot_PDCumulants_TPR.C:1257
Plot_PDCumulants_TPR.C:1258
Plot_PDCumulants_TPR.C:1259
Plot_PDCumulants_TPR.C:1260
Plot_PDCumulants_TPR.C:1261
Plot_PDCumulants_TPR.C:1262
Plot_PDCumulants_TPR.C:1263
Plot_PDCumulants_TPR.C:1264
Plot_PDCumulants_TPR.C:1265
Plot_PDCumulants_TPR.C:1266
Plot_PDCumulants_TPR.C:1267
Plot_PDCumulants_TPR.C:1268
Plot_PDCumulants_TPR.C:1269
Plot_PDCumulants_TPR.C:1270
Plot_PDCumulants_TPR.C:1271
Plot_PDCumulants_TPR.C:1272
Plot_PDCumulants_TPR.C:1273
Plot_PDCumulants_TPR.C:1274
Plot_PDCumulants_TPR.C:1275
Plot_PDCumulants_TPR.C:1276
Plot_PDCumulants_TPR.C:1277
Plot_PDCumulants_TPR.C:1278
Plot_PDCumulants_TPR.C:1279
Plot_PDCumulants_TPR.C:1280
Plot_PDCumulants_TPR.C:1281
Plot_PDCumulants_TPR.C:1282
Plot_PDCumulants_TPR.C:1283
Plot_PDCumulants_TPR.C:1284
Plot_PDCumulants_TPR.C:1285
Plot_PDCumulants_TPR.C:1286
Plot_PDCumulants_TPR.C:1287
Plot_PDCumulants_TPR.C:1288
Plot_PDCumulants_TPR.C:1289
Plot_PDCumulants_TPR.C:1290
Plot_PDCumulants_TPR.C:1291
Plot_PDCumulants_TPR.C:1292
Plot_PDCumulants_TPR.C:1293
Plot_PDCumulants_TPR.C:1294
Plot_PDCumulants_TPR.C:1295
Plot_PDCumulants_TPR.C:1296
Plot_PDCumulants_TPR.C:1297
Plot_PDCumulants_TPR.C:1298
Plot_PDCumulants_TPR.C:1299
Plot_PDCumulants_TPR.C:1300
Plot_PDCumulants_TPR.C:1301
Plot_PDCumulants_TPR.C:1302
Plot_PDCumulants_TPR.C:1303
Plot_PDCumulants_TPR.C:1304
Plot_PDCumulants_TPR.C:1305
Plot_PDCumulants_TPR.C:1306
Plot_PDCumulants_TPR.C:1307
Plot_PDCumulants_TPR.C:1308
Plot_PDCumulants_TPR.C:1309
Plot_PDCumulants_TPR.C:1310
Plot_PDCumulants_TPR.C:1311
Plot_PDCumulants_TPR.C:1312
Plot_PDCumulants_TPR.C:1313
Plot_PDCumulants_TPR.C:1314
Plot_PDCumulants_TPR.C:1315
Plot_PDCumulants_TPR.C:1316
Plot_PDCumulants_TPR.C:1317
Plot_PDCumulants_TPR.C:1318
Plot_PDCumulants_TPR.C:1319
Plot_PDCumulants_TPR.C:1320
Plot_PDCumulants_TPR.C:1321
Plot_PDCumulants_TPR.C:1322
Plot_PDCumulants_TPR.C:1323
Plot_PDCumulants_TPR.C:1324
Plot_PDCumulants_TPR.C:1325
Plot_PDCumulants_TPR.C:1326
Plot_PDCumulants_TPR.C:1327
Plot_PDCumulants_TPR.C:1328
Plot_PDCumulants_TPR.C:1329
Plot_PDCumulants_TPR.C:1330
Plot_PDCumulants_TPR.C:1331
Plot_PDCumulants_TPR.C:1332
Plot_PDCumulants_TPR.C:1333
Plot_PDCumulants_TPR.C:1334
Plot_PDCumulants_TPR.C:1335
Plot_PDCumulants_TPR.C:1336
Plot_PDCumulants_TPR.C:1337
Plot_PDCumulants_TPR.C:1338
Plot_PDCumulants_TPR.C:1339
Plot_PDCumulants_TPR.C:1340
Plot_PDCumulants_TPR.C:1341
Plot_PDCumulants_TPR.C:1342
Plot_PDCumulants_TPR.C:1343
Plot_PDCumulants_TPR.C:1344
Plot_PDCumulants_TPR.C:1345
Plot_PDCumulants_TPR.C:1346
Plot_PDCumulants_TPR.C:1347
Plot_PDCumulants_TPR.C:1348
Plot_PDCumulants_TPR.C:1349
Plot_PDCumulants_TPR.C:1350
Plot_PDCumulants_TPR.C:1351
Plot_PDCumulants_TPR.C:1352
Plot_PDCumulants_TPR.C:1353
Plot_PDCumulants_TPR.C:1354
Plot_PDCumulants_TPR.C:1355
Plot_PDCumulants_TPR.C:1356
Plot_PDCumulants_TPR.C:1357
Plot_PDCumulants_TPR.C:1358
Plot_PDCumulants_TPR.C:1359
Plot_PDCumulants_TPR.C:1360
Plot_PDCumulants_TPR.C:1361
Plot_PDCumulants_TPR.C:1362
Plot_PDCumulants_TPR.C:1363
Plot_PDCumulants_TPR.C:1364
Plot_PDCumulants_TPR.C:1365
Plot_PDCumulants_TPR.C:1366
Plot_PDCumulants_TPR.C:1367
Plot_PDCumulants_TPR.C:1368
Plot_PDCumulants_TPR.C:1369
Plot_PDCumulants_TPR.C:1370
Plot_PDCumulants_TPR.C:1371
Plot_PDCumulants_TPR.C:1372
Plot_PDCumulants_TPR.C:1373
Plot_PDCumulants_TPR.C:1374
Plot_PDCumulants_TPR.C:1375
Plot_PDCumulants_TPR.C:1376
Plot_PDCumulants_TPR.C:1377
Plot_PDCumulants_TPR.C:1378
Plot_PDCumulants_TPR.C:1379
Plot_PDCumulants_TPR.C:1380
Plot_PDCumulants_TPR.C:1381
Plot_PDCumulants_TPR.C:1382
Plot_PDCumulants_TPR.C:1383
Plot_PDCumulants_TPR.C:1384
Plot_PDCumulants_TPR.C:1385
Plot_PDCumulants_TPR.C:1386
Plot_PDCumulants_TPR.C:1387
Plot_PDCumulants_TPR.C:1388
Plot_PDCumulants_TPR.C:1389
Plot_PDCumulants_TPR.C:1390
Plot_PDCumulants_TPR.C:1391
Plot_PDCumulants_TPR.C:1392
Plot_PDCumulants_TPR.C:1393
Plot_PDCumulants_TPR.C:1394
Plot_PDCumulants_TPR.C:1395
Plot_PDCumulants_TPR.C:1396
Plot_PDCumulants_TPR.C:1397
Plot_PDCumulants_TPR.C:1398
Plot_PDCumulants_TPR.C:1399
Plot_PDCumulants_TPR.C:1400
Plot_PDCumulants_TPR.C:1401
Plot_PDCumulants_TPR.C:1402
Plot_PDCumulants_TPR.C:1403
Plot_PDCumulants_TPR.C:1404
Plot_PDCumulants_TPR.C:1405
Plot_PDCumulants_TPR.C:1406
Plot_PDCumulants_TPR.C:1407
Plot_PDCumulants_TPR.C:1408
Plot_PDCumulants_TPR.C:1409
Plot_PDCumulants_TPR.C:1410
Plot_PDCumulants_TPR.C:1411
Plot_PDCumulants_TPR.C:1412
Plot_PDCumulants_TPR.C:1413
Plot_PDCumulants_TPR.C:1414
Plot_PDCumulants_TPR.C:1415
Plot_PDCumulants_TPR.C:1416
Plot_PDCumulants_TPR.C:1417
Plot_PDCumulants_TPR.C:1418
Plot_PDCumulants_TPR.C:1419
Plot_PDCumulants_TPR.C:1420
Plot_PDCumulants_TPR.C:1421
Plot_PDCumulants_TPR.C:1422
Plot_PDCumulants_TPR.C:1423
Plot_PDCumulants_TPR.C:1424
Plot_PDCumulants_TPR.C:1425
Plot_PDCumulants_TPR.C:1426
Plot_PDCumulants_TPR.C:1427
Plot_PDCumulants_TPR.C:1428
Plot_PDCumulants_TPR.C:1429
Plot_PDCumulants_TPR.C:1430
Plot_PDCumulants_TPR.C:1431
Plot_PDCumulants_TPR.C:1432
Plot_PDCumulants_TPR.C:1433
Plot_PDCumulants_TPR.C:1434
Plot_PDCumulants_TPR.C:1435
Plot_PDCumulants_TPR.C:1436
Plot_PDCumulants_TPR.C:1437
Plot_PDCumulants_TPR.C:1438
Plot_PDCumulants_TPR.C:1439
Plot_PDCumulants_TPR.C:1440
Plot_PDCumulants_TPR.C:1441
Plot_PDCumulants_TPR.C:1442
Plot_PDCumulants_TPR.C:1443
Plot_PDCumulants_TPR.C:1444
Plot_PDCumulants_TPR.C:1445
Plot_PDCumulants_TPR.C:1446
Plot_PDCumulants_TPR.C:1447
Plot_PDCumulants_TPR.C:1448
Plot_PDCumulants_TPR.C:1449
Plot_PDCumulants_TPR.C:1450
Plot_PDCumulants_TPR.C:1451
Plot_PDCumulants_TPR.C:1452
Plot_PDCumulants_TPR.C:1453
Plot_PDCumulants_TPR.C:1454
Plot_PDCumulants_TPR.C:1455
Plot_PDCumulants_TPR.C:1456
Plot_PDCumulants_TPR.C:1457
Plot_PDCumulants_TPR.C:1458
Plot_PDCumulants_TPR.C:1459
Plot_PDCumulants_TPR.C:1460
Plot_PDCumulants_TPR.C:1461
Plot_PDCumulants_TPR.C:1462
Plot_PDCumulants_TPR.C:1463
Plot_PDCumulants_TPR.C:1464
Plot_PDCumulants_TPR.C:1465
Plot_PDCumulants_TPR.C:1466
Plot_PDCumulants_TPR.C:1467
Plot_PDCumulants_TPR.C:1468
Plot_PDCumulants_TPR.C:1469
Plot_PDCumulants_TPR.C:1470
Plot_PDCumulants_TPR.C:1471
Plot_PDCumulants_TPR.C:1472
Plot_PDCumulants_TPR.C:1473
Plot_PDCumulants_TPR.C:1474
Plot_PDCumulants_TPR.C:1475
Plot_PDCumulants_TPR.C:1476
Plot_PDCumulants_TPR.C:1477
Plot_PDCumulants_TPR.C:1478
Plot_PDCumulants_TPR.C:1479
Plot_PDCumulants_TPR.C:1480
Plot_PDCumulants_TPR.C:1481
Plot_PDCumulants_TPR.C:1482
Plot_PDCumulants_TPR.C:1483
Plot_PDCumulants_TPR.C:1484
Plot_PDCumulants_TPR.C:1485
Plot_PDCumulants_TPR.C:1486
Plot_PDCumulants_TPR.C:1487
Plot_PDCumulants_TPR.C:1488
Plot_PDCumulants_TPR.C:1489
Plot_PDCumulants_TPR.C:1490
Plot_PDCumulants_TPR.C:1491
Plot_PDCumulants_TPR.C:1492
Plot_PDCumulants_TPR.C:1493
Plot_PDCumulants_TPR.C:1494
Plot_PDCumulants_TPR.C:1495
Plot_PDCumulants_TPR.C:1496
Plot_PDCumulants_TPR.C:1497
Plot_PDCumulants_TPR.C:1498
Plot_PDCumulants_TPR.C:1499
Plot_PDCumulants_TPR.C:1500
Plot_PDCumulants_TPR.C:1501
Plot_PDCumulants_TPR.C:1502
Plot_PDCumulants_TPR.C:1503
Plot_PDCumulants_TPR.C:1504
Plot_PDCumulants_TPR.C:1505
Plot_PDCumulants_TPR.C:1506
Plot_PDCumulants_TPR.C:1507
Plot_PDCumulants_TPR.C:1508
Plot_PDCumulants_TPR.C:1509
Plot_PDCumulants_TPR.C:1510
Plot_PDCumulants_TPR.C:1511
Plot_PDCumulants_TPR.C:1512
Plot_PDCumulants_TPR.C:1513
Plot_PDCumulants_TPR.C:1514
Plot_PDCumulants_TPR.C:1515
Plot_PDCumulants_TPR.C:1516
Plot_PDCumulants_TPR.C:1517
Plot_PDCumulants_TPR.C:1518
Plot_PDCumulants_TPR.C:1519
Plot_PDCumulants_TPR.C:1520
Plot_PDCumulants_TPR.C:1521
Plot_PDCumulants_TPR.C:1522
Plot_PDCumulants_TPR.C:1523
Plot_PDCumulants_TPR.C:1524
Plot_PDCumulants_TPR.C:1525
Plot_PDCumulants_TPR.C:1526
Plot_PDCumulants_TPR.C:1527
Plot_PDCumulants_TPR.C:1528
Plot_PDCumulants_TPR.C:1529
Plot_PDCumulants_TPR.C:1530
Plot_PDCumulants_TPR.C:1531
Plot_PDCumulants_TPR.C:1532
Plot_PDCumulants_TPR.C:1533
Plot_PDCumulants_TPR.C:1534
Plot_PDCumulants_TPR.C:1535
Plot_PDCumulants_TPR.C:1536
Plot_PDCumulants_TPR.C:1537
Plot_PDCumulants_TPR.C:1538
Plot_PDCumulants_TPR.C:1539
Plot_PDCumulants_TPR.C:1540
Plot_PDCumulants_TPR.C:1541
Plot_PDCumulants_TPR.C:1542
Plot_PDCumulants_TPR.C:1543
Plot_PDCumulants_TPR.C:1544
Plot_PDCumulants_TPR.C:1545
Plot_PDCumulants_TPR.C:1546
Plot_PDCumulants_TPR.C:1547
Plot_PDCumulants_TPR.C:1548
Plot_PDCumulants_TPR.C:1549
Plot_PDCumulants_TPR.C:1550
Plot_PDCumulants_TPR.C:1551
Plot_PDCumulants_TPR.C:1552
Plot_PDCumulants_TPR.C:1553
Plot_PDCumulants_TPR.C:1554
Plot_PDCumulants_TPR.C:1555
Plot_PDCumulants_TPR.C:1556
Plot_PDCumulants_TPR.C:1557
Plot_PDCumulants_TPR.C:1558
Plot_PDCumulants_TPR.C:1559
Plot_PDCumulants_TPR.C:1560
Plot_PDCumulants_TPR.C:1561
Plot_PDCumulants_TPR.C:1562
Plot_PDCumulants_TPR.C:1563
Plot_PDCumulants_TPR.C:1564
Plot_PDCumulants_TPR.C:1565
Plot_PDCumulants_TPR.C:1566
Plot_PDCumulants_TPR.C:1567
Plot_PDCumulants_TPR.C:1568
Plot_PDCumulants_TPR.C:1569
Plot_PDCumulants_TPR.C:1570
Plot_PDCumulants_TPR.C:1571
Plot_PDCumulants_TPR.C:1572
Plot_PDCumulants_TPR.C:1573
Plot_PDCumulants_TPR.C:1574
Plot_PDCumulants_TPR.C:1575
Plot_PDCumulants_TPR.C:1576
Plot_PDCumulants_TPR.C:1577
Plot_PDCumulants_TPR.C:1578
Plot_PDCumulants_TPR.C:1579
Plot_PDCumulants_TPR.C:1580
Plot_PDCumulants_TPR.C:1581
Plot_PDCumulants_TPR.C:1582
Plot_PDCumulants_TPR.C:1583
Plot_PDCumulants_TPR.C:1584
Plot_PDCumulants_TPR.C:1585
Plot_PDCumulants_TPR.C:1586
Plot_PDCumulants_TPR.C:1587
Plot_PDCumulants_TPR.C:1588
Plot_PDCumulants_TPR.C:1589
Plot_PDCumulants_TPR.C:1590
Plot_PDCumulants_TPR.C:1591
Plot_PDCumulants_TPR.C:1592
Plot_PDCumulants_TPR.C:1593
Plot_PDCumulants_TPR.C:1594
Plot_PDCumulants_TPR.C:1595
Plot_PDCumulants_TPR.C:1596
Plot_PDCumulants_TPR.C:1597
Plot_PDCumulants_TPR.C:1598
Plot_PDCumulants_TPR.C:1599
Plot_PDCumulants_TPR.C:1600
Plot_PDCumulants_TPR.C:1601
Plot_PDCumulants_TPR.C:1602
Plot_PDCumulants_TPR.C:1603
Plot_PDCumulants_TPR.C:1604
Plot_PDCumulants_TPR.C:1605
Plot_PDCumulants_TPR.C:1606
Plot_PDCumulants_TPR.C:1607
Plot_PDCumulants_TPR.C:1608
Plot_PDCumulants_TPR.C:1609
Plot_PDCumulants_TPR.C:1610
Plot_PDCumulants_TPR.C:1611
Plot_PDCumulants_TPR.C:1612
Plot_PDCumulants_TPR.C:1613
Plot_PDCumulants_TPR.C:1614
Plot_PDCumulants_TPR.C:1615
Plot_PDCumulants_TPR.C:1616
Plot_PDCumulants_TPR.C:1617
Plot_PDCumulants_TPR.C:1618
Plot_PDCumulants_TPR.C:1619
Plot_PDCumulants_TPR.C:1620
Plot_PDCumulants_TPR.C:1621
Plot_PDCumulants_TPR.C:1622
Plot_PDCumulants_TPR.C:1623
Plot_PDCumulants_TPR.C:1624
Plot_PDCumulants_TPR.C:1625
Plot_PDCumulants_TPR.C:1626
Plot_PDCumulants_TPR.C:1627
Plot_PDCumulants_TPR.C:1628
Plot_PDCumulants_TPR.C:1629
Plot_PDCumulants_TPR.C:1630
Plot_PDCumulants_TPR.C:1631
Plot_PDCumulants_TPR.C:1632
Plot_PDCumulants_TPR.C:1633
Plot_PDCumulants_TPR.C:1634
Plot_PDCumulants_TPR.C:1635
Plot_PDCumulants_TPR.C:1636
Plot_PDCumulants_TPR.C:1637
Plot_PDCumulants_TPR.C:1638
Plot_PDCumulants_TPR.C:1639
Plot_PDCumulants_TPR.C:1640
Plot_PDCumulants_TPR.C:1641
Plot_PDCumulants_TPR.C:1642
Plot_PDCumulants_TPR.C:1643
Plot_PDCumulants_TPR.C:1644
Plot_PDCumulants_TPR.C:1645
Plot_PDCumulants_TPR.C:1646
Plot_PDCumulants_TPR.C:1647
Plot_PDCumulants_TPR.C:1648
Plot_PDCumulants_TPR.C:1649
Plot_PDCumulants_TPR.C:1650
Plot_PDCumulants_TPR.C:1651
Plot_PDCumulants_TPR.C:1652
Plot_PDCumulants_TPR.C:1653
Plot_PDCumulants_TPR.C:1654
Plot_PDCumulants_TPR.C:1655
Plot_PDCumulants_TPR.C:1656
Plot_PDCumulants_TPR.C:1657
Plot_PDCumulants_TPR.C:1658
Plot_PDCumulants_TPR.C:1659
Plot_PDCumulants_TPR.C:1660
Plot_PDCumulants_TPR.C:1661
Plot_PDCumulants_TPR.C:1662
Plot_PDCumulants_TPR.C:1663
Plot_PDCumulants_TPR.C:1664
Plot_PDCumulants_TPR.C:1665
Plot_PDCumulants_TPR.C:1666
Plot_PDCumulants_TPR.C:1667
Plot_PDCumulants_TPR.C:1668
Plot_PDCumulants_TPR.C:1669
Plot_PDCumulants_TPR.C:1670
Plot_PDCumulants_TPR.C:1671
Plot_PDCumulants_TPR.C:1672
Plot_PDCumulants_TPR.C:1673
Plot_PDCumulants_TPR.C:1674
Plot_PDCumulants_TPR.C:1675
Plot_PDCumulants_TPR.C:1676
Plot_PDCumulants_TPR.C:1677
Plot_PDCumulants_TPR.C:1678
Plot_PDCumulants_TPR.C:1679
Plot_PDCumulants_TPR.C:1680
Plot_PDCumulants_TPR.C:1681
Plot_PDCumulants_TPR.C:1682
Plot_PDCumulants_TPR.C:1683
Plot_PDCumulants_TPR.C:1684
Plot_PDCumulants_TPR.C:1685
Plot_PDCumulants_TPR.C:1686
Plot_PDCumulants_TPR.C:1687
Plot_PDCumulants_TPR.C:1688
Plot_PDCumulants_TPR.C:1689
Plot_PDCumulants_TPR.C:1690
Plot_PDCumulants_TPR.C:1691
Plot_PDCumulants_TPR.C:1692
Plot_PDCumulants_TPR.C:1693
Plot_PDCumulants_TPR.C:1694
Plot_PDCumulants_TPR.C:1695
Plot_PDCumulants_TPR.C:1696
Plot_PDCumulants_TPR.C:1697
Plot_PDCumulants_TPR.C:1698
Plot_PDCumulants_TPR.C:1699
Plot_PDCumulants_TPR.C:1700
Plot_PDCumulants_TPR.C:1701
Plot_PDCumulants_TPR.C:1702
Plot_PDCumulants_TPR.C:1703
Plot_PDCumulants_TPR.C:1704
Plot_PDCumulants_TPR.C:1705
Plot_PDCumulants_TPR.C:1706
Plot_PDCumulants_TPR.C:1707
Plot_PDCumulants_TPR.C:1708
Plot_PDCumulants_TPR.C:1709
Plot_PDCumulants_TPR.C:1710
Plot_PDCumulants_TPR.C:1711
Plot_PDCumulants_TPR.C:1712
Plot_PDCumulants_TPR.C:1713
Plot_PDCumulants_TPR.C:1714
Plot_PDCumulants_TPR.C:1715
Plot_PDCumulants_TPR.C:1716
Plot_PDCumulants_TPR.C:1717
Plot_PDCumulants_TPR.C:1718
Plot_PDCumulants_TPR.C:1719
Plot_PDCumulants_TPR.C:1720
Plot_PDCumulants_TPR.C:1721
Plot_PDCumulants_TPR.C:1722
Plot_PDCumulants_TPR.C:1723
Plot_PDCumulants_TPR.C:1724
Plot_PDCumulants_TPR.C:1725
Plot_PDCumulants_TPR.C:1726
Plot_PDCumulants_TPR.C:1727
Plot_PDCumulants_TPR.C:1728
Plot_PDCumulants_TPR.C:1729
Plot_PDCumulants_TPR.C:1730
Plot_PDCumulants_TPR.C:1731
Plot_PDCumulants_TPR.C:1732
Plot_PDCumulants_TPR.C:1733
Plot_PDCumulants_TPR.C:1734
Plot_PDCumulants_TPR.C:1735
Plot_PDCumulants_TPR.C:1736
Plot_PDCumulants_TPR.C:1737
Plot_PDCumulants_TPR.C:1738
Plot_PDCumulants_TPR.C:1739
Plot_PDCumulants_TPR.C:1740
Plot_PDCumulants_TPR.C:1741
Plot_PDCumulants_TPR.C:1742
Plot_PDCumulants_TPR.C:1743
Plot_PDCumulants_TPR.C:1744
Plot_PDCumulants_TPR.C:1745
Plot_PDCumulants_TPR.C:1746
Plot_PDCumulants_TPR.C:1747
Plot_PDCumulants_TPR.C:1748
Plot_PDCumulants_TPR.C:1749
Plot_PDCumulants_TPR.C:1750
Plot_PDCumulants_TPR.C:1751
Plot_PDCumulants_TPR.C:1752
Plot_PDCumulants_TPR.C:1753
Plot_PDCumulants_TPR.C:1754
Plot_PDCumulants_TPR.C:1755
Plot_PDCumulants_TPR.C:1756
Plot_PDCumulants_TPR.C:1757
Plot_PDCumulants_TPR.C:1758
Plot_PDCumulants_TPR.C:1759
Plot_PDCumulants_TPR.C:1760
Plot_PDCumulants_TPR.C:1761
Plot_PDCumulants_TPR.C:1762
Plot_PDCumulants_TPR.C:1763
Plot_PDCumulants_TPR.C:1764
Plot_PDCumulants_TPR.C:1765
Plot_PDCumulants_TPR.C:1766
Plot_PDCumulants_TPR.C:1767
Plot_PDCumulants_TPR.C:1768
Plot_PDCumulants_TPR.C:1769
Plot_PDCumulants_TPR.C:1770
Plot_PDCumulants_TPR.C:1771
Plot_PDCumulants_TPR.C:1772
Plot_PDCumulants_TPR.C:1773
Plot_PDCumulants_TPR.C:1774
Plot_PDCumulants_TPR.C:1775
Plot_PDCumulants_TPR.C:1776
Plot_PDCumulants_TPR.C:1777
Plot_PDCumulants_TPR.C:1778
Plot_PDCumulants_TPR.C:1779
Plot_PDCumulants_TPR.C:1780
Plot_PDCumulants_TPR.C:1781
Plot_PDCumulants_TPR.C:1782
Plot_PDCumulants_TPR.C:1783
Plot_PDCumulants_TPR.C:1784
Plot_PDCumulants_TPR.C:1785
Plot_PDCumulants_TPR.C:1786
Plot_PDCumulants_TPR.C:1787
Plot_PDCumulants_TPR.C:1788
Plot_PDCumulants_TPR.C:1789
Plot_PDCumulants_TPR.C:1790
Plot_PDCumulants_TPR.C:1791
Plot_PDCumulants_TPR.C:1792
Plot_PDCumulants_TPR.C:1793
Plot_PDCumulants_TPR.C:1794
Plot_PDCumulants_TPR.C:1795
Plot_PDCumulants_TPR.C:1796
Plot_PDCumulants_TPR.C:1797
Plot_PDCumulants_TPR.C:1798
Plot_PDCumulants_TPR.C:1799
Plot_PDCumulants_TPR.C:1800
Plot_PDCumulants_TPR.C:1801
Plot_PDCumulants_TPR.C:1802
Plot_PDCumulants_TPR.C:1803
Plot_PDCumulants_TPR.C:1804
Plot_PDCumulants_TPR.C:1805
Plot_PDCumulants_TPR.C:1806
Plot_PDCumulants_TPR.C:1807
Plot_PDCumulants_TPR.C:1808
Plot_PDCumulants_TPR.C:1809
Plot_PDCumulants_TPR.C:1810
Plot_PDCumulants_TPR.C:1811
Plot_PDCumulants_TPR.C:1812
Plot_PDCumulants_TPR.C:1813
Plot_PDCumulants_TPR.C:1814
Plot_PDCumulants_TPR.C:1815
Plot_PDCumulants_TPR.C:1816
Plot_PDCumulants_TPR.C:1817
Plot_PDCumulants_TPR.C:1818
Plot_PDCumulants_TPR.C:1819
Plot_PDCumulants_TPR.C:1820
Plot_PDCumulants_TPR.C:1821
Plot_PDCumulants_TPR.C:1822
Plot_PDCumulants_TPR.C:1823
Plot_PDCumulants_TPR.C:1824
Plot_PDCumulants_TPR.C:1825
Plot_PDCumulants_TPR.C:1826
Plot_PDCumulants_TPR.C:1827
Plot_PDCumulants_TPR.C:1828
Plot_PDCumulants_TPR.C:1829
Plot_PDCumulants_TPR.C:1830
Plot_PDCumulants_TPR.C:1831
Plot_PDCumulants_TPR.C:1832
Plot_PDCumulants_TPR.C:1833
Plot_PDCumulants_TPR.C:1834
Plot_PDCumulants_TPR.C:1835
Plot_PDCumulants_TPR.C:1836
Plot_PDCumulants_TPR.C:1837
Plot_PDCumulants_TPR.C:1838
Plot_PDCumulants_TPR.C:1839
Plot_PDCumulants_TPR.C:1840
Plot_PDCumulants_TPR.C:1841
Plot_PDCumulants_TPR.C:1842
Plot_PDCumulants_TPR.C:1843
Plot_PDCumulants_TPR.C:1844
Plot_PDCumulants_TPR.C:1845
Plot_PDCumulants_TPR.C:1846
Plot_PDCumulants_TPR.C:1847
Plot_PDCumulants_TPR.C:1848
Plot_PDCumulants_TPR.C:1849
Plot_PDCumulants_TPR.C:1850
Plot_PDCumulants_TPR.C:1851
Plot_PDCumulants_TPR.C:1852
Plot_PDCumulants_TPR.C:1853
Plot_PDCumulants_TPR.C:1854
Plot_PDCumulants_TPR.C:1855
Plot_PDCumulants_TPR.C:1856
Plot_PDCumulants_TPR.C:1857
Plot_PDCumulants_TPR.C:1858
Plot_PDCumulants_TPR.C:1859
Plot_PDCumulants_TPR.C:1860
Plot_PDCumulants_TPR.C:1861
Plot_PDCumulants_TPR.C:1862
Plot_PDCumulants_TPR.C:1863
Plot_PDCumulants_TPR.C:1864
Plot_PDCumulants_TPR.C:1865
Plot_PDCumulants_TPR.C:1866
Plot_PDCumulants_TPR.C:1867
Plot_PDCumulants_TPR.C:1868
Plot_PDCumulants_TPR.C:1869
Plot_PDCumulants_TPR.C:1870
Plot_PDCumulants_TPR.C:1871
Plot_PDCumulants_TPR.C:1872
Plot_PDCumulants_TPR.C:1873
Plot_PDCumulants_TPR.C:1874
Plot_PDCumulants_TPR.C:1875
Plot_PDCumulants_TPR.C:1876
Plot_PDCumulants_TPR.C:1877
Plot_PDCumulants_TPR.C:1878
Plot_PDCumulants_TPR.C:1879
Plot_PDCumulants_TPR.C:1880
Plot_PDCumulants_TPR.C:1881
Plot_PDCumulants_TPR.C:1882
Plot_PDCumulants_TPR.C:1883
Plot_PDCumulants_TPR.C:1884
Plot_PDCumulants_TPR.C:1885
Plot_PDCumulants_TPR.C:1886
Plot_PDCumulants_TPR.C:1887
Plot_PDCumulants_TPR.C:1888
Plot_PDCumulants_TPR.C:1889
Plot_PDCumulants_TPR.C:1890
Plot_PDCumulants_TPR.C:1891
Plot_PDCumulants_TPR.C:1892
Plot_PDCumulants_TPR.C:1893
Plot_PDCumulants_TPR.C:1894
Plot_PDCumulants_TPR.C:1895
Plot_PDCumulants_TPR.C:1896
Plot_PDCumulants_TPR.C:1897
Plot_PDCumulants_TPR.C:1898
Plot_PDCumulants_TPR.C:1899
Plot_PDCumulants_TPR.C:1900
Plot_PDCumulants_TPR.C:1901
Plot_PDCumulants_TPR.C:1902
Plot_PDCumulants_TPR.C:1903
Plot_PDCumulants_TPR.C:1904
Plot_PDCumulants_TPR.C:1905
Plot_PDCumulants_TPR.C:1906
Plot_PDCumulants_TPR.C:1907
Plot_PDCumulants_TPR.C:1908
Plot_PDCumulants_TPR.C:1909
Plot_PDCumulants_TPR.C:1910
Plot_PDCumulants_TPR.C:1911
Plot_PDCumulants_TPR.C:1912
Plot_PDCumulants_TPR.C:1913
Plot_PDCumulants_TPR.C:1914
Plot_PDCumulants_TPR.C:1915
Plot_PDCumulants_TPR.C:1916
Plot_PDCumulants_TPR.C:1917
Plot_PDCumulants_TPR.C:1918
Plot_PDCumulants_TPR.C:1919
Plot_PDCumulants_TPR.C:1920
Plot_PDCumulants_TPR.C:1921
Plot_PDCumulants_TPR.C:1922
Plot_PDCumulants_TPR.C:1923
Plot_PDCumulants_TPR.C:1924
Plot_PDCumulants_TPR.C:1925
Plot_PDCumulants_TPR.C:1926
Plot_PDCumulants_TPR.C:1927
Plot_PDCumulants_TPR.C:1928
Plot_PDCumulants_TPR.C:1929
Plot_PDCumulants_TPR.C:1930
Plot_PDCumulants_TPR.C:1931
Plot_PDCumulants_TPR.C:1932
Plot_PDCumulants_TPR.C:1933
Plot_PDCumulants_TPR.C:1934
Plot_PDCumulants_TPR.C:1935
Plot_PDCumulants_TPR.C:1936
Plot_PDCumulants_TPR.C:1937
Plot_PDCumulants_TPR.C:1938
Plot_PDCumulants_TPR.C:1939
Plot_PDCumulants_TPR.C:1940
Plot_PDCumulants_TPR.C:1941
Plot_PDCumulants_TPR.C:1942
Plot_PDCumulants_TPR.C:1943
Plot_PDCumulants_TPR.C:1944
Plot_PDCumulants_TPR.C:1945
Plot_PDCumulants_TPR.C:1946
Plot_PDCumulants_TPR.C:1947
Plot_PDCumulants_TPR.C:1948
Plot_PDCumulants_TPR.C:1949
Plot_PDCumulants_TPR.C:1950
Plot_PDCumulants_TPR.C:1951
Plot_PDCumulants_TPR.C:1952
Plot_PDCumulants_TPR.C:1953
Plot_PDCumulants_TPR.C:1954
Plot_PDCumulants_TPR.C:1955
Plot_PDCumulants_TPR.C:1956
Plot_PDCumulants_TPR.C:1957
Plot_PDCumulants_TPR.C:1958
Plot_PDCumulants_TPR.C:1959
Plot_PDCumulants_TPR.C:1960
Plot_PDCumulants_TPR.C:1961
Plot_PDCumulants_TPR.C:1962
Plot_PDCumulants_TPR.C:1963
Plot_PDCumulants_TPR.C:1964
Plot_PDCumulants_TPR.C:1965
Plot_PDCumulants_TPR.C:1966
Plot_PDCumulants_TPR.C:1967
Plot_PDCumulants_TPR.C:1968
Plot_PDCumulants_TPR.C:1969
Plot_PDCumulants_TPR.C:1970
Plot_PDCumulants_TPR.C:1971
Plot_PDCumulants_TPR.C:1972
Plot_PDCumulants_TPR.C:1973
Plot_PDCumulants_TPR.C:1974
Plot_PDCumulants_TPR.C:1975
Plot_PDCumulants_TPR.C:1976
Plot_PDCumulants_TPR.C:1977
Plot_PDCumulants_TPR.C:1978
Plot_PDCumulants_TPR.C:1979
Plot_PDCumulants_TPR.C:1980
Plot_PDCumulants_TPR.C:1981
Plot_PDCumulants_TPR.C:1982
Plot_PDCumulants_TPR.C:1983
Plot_PDCumulants_TPR.C:1984
Plot_PDCumulants_TPR.C:1985
Plot_PDCumulants_TPR.C:1986
Plot_PDCumulants_TPR.C:1987
Plot_PDCumulants_TPR.C:1988
Plot_PDCumulants_TPR.C:1989
Plot_PDCumulants_TPR.C:1990
Plot_PDCumulants_TPR.C:1991
Plot_PDCumulants_TPR.C:1992
Plot_PDCumulants_TPR.C:1993
Plot_PDCumulants_TPR.C:1994
Plot_PDCumulants_TPR.C:1995
Plot_PDCumulants_TPR.C:1996
Plot_PDCumulants_TPR.C:1997
Plot_PDCumulants_TPR.C:1998
Plot_PDCumulants_TPR.C:1999
Plot_PDCumulants_TPR.C:2000
Plot_PDCumulants_TPR.C:2001
Plot_PDCumulants_TPR.C:2002
Plot_PDCumulants_TPR.C:2003
Plot_PDCumulants_TPR.C:2004
Plot_PDCumulants_TPR.C:2005
Plot_PDCumulants_TPR.C:2006
Plot_PDCumulants_TPR.C:2007
Plot_PDCumulants_TPR.C:2008
Plot_PDCumulants_TPR.C:2009
Plot_PDCumulants_TPR.C:2010
Plot_PDCumulants_TPR.C:2011
Plot_PDCumulants_TPR.C:2012
Plot_PDCumulants_TPR.C:2013
Plot_PDCumulants_TPR.C:2014
Plot_PDCumulants_TPR.C:2015
Plot_PDCumulants_TPR.C:2016
Plot_PDCumulants_TPR.C:2017
Plot_PDCumulants_TPR.C:2018
Plot_PDCumulants_TPR.C:2019
Plot_PDCumulants_TPR.C:2020
Plot_PDCumulants_TPR.C:2021
Plot_PDCumulants_TPR.C:2022
Plot_PDCumulants_TPR.C:2023
Plot_PDCumulants_TPR.C:2024
Plot_PDCumulants_TPR.C:2025
Plot_PDCumulants_TPR.C:2026
Plot_PDCumulants_TPR.C:2027
Plot_PDCumulants_TPR.C:2028
Plot_PDCumulants_TPR.C:2029
Plot_PDCumulants_TPR.C:2030
Plot_PDCumulants_TPR.C:2031
Plot_PDCumulants_TPR.C:2032
Plot_PDCumulants_TPR.C:2033
Plot_PDCumulants_TPR.C:2034
Plot_PDCumulants_TPR.C:2035
Plot_PDCumulants_TPR.C:2036
Plot_PDCumulants_TPR.C:2037
Plot_PDCumulants_TPR.C:2038
Plot_PDCumulants_TPR.C:2039
Plot_PDCumulants_TPR.C:2040
Plot_PDCumulants_TPR.C:2041
Plot_PDCumulants_TPR.C:2042
Plot_PDCumulants_TPR.C:2043
Plot_PDCumulants_TPR.C:2044
Plot_PDCumulants_TPR.C:2045
Plot_PDCumulants_TPR.C:2046
Plot_PDCumulants_TPR.C:2047
Plot_PDCumulants_TPR.C:2048
Plot_PDCumulants_TPR.C:2049
Plot_PDCumulants_TPR.C:2050
Plot_PDCumulants_TPR.C:2051
Plot_PDCumulants_TPR.C:2052
Plot_PDCumulants_TPR.C:2053
Plot_PDCumulants_TPR.C:2054
Plot_PDCumulants_TPR.C:2055
Plot_PDCumulants_TPR.C:2056
Plot_PDCumulants_TPR.C:2057
Plot_PDCumulants_TPR.C:2058
Plot_PDCumulants_TPR.C:2059
Plot_PDCumulants_TPR.C:2060
Plot_PDCumulants_TPR.C:2061
Plot_PDCumulants_TPR.C:2062
Plot_PDCumulants_TPR.C:2063
Plot_PDCumulants_TPR.C:2064
Plot_PDCumulants_TPR.C:2065
Plot_PDCumulants_TPR.C:2066
Plot_PDCumulants_TPR.C:2067
Plot_PDCumulants_TPR.C:2068
Plot_PDCumulants_TPR.C:2069
Plot_PDCumulants_TPR.C:2070
Plot_PDCumulants_TPR.C:2071
Plot_PDCumulants_TPR.C:2072
Plot_PDCumulants_TPR.C:2073
Plot_PDCumulants_TPR.C:2074
Plot_PDCumulants_TPR.C:2075
Plot_PDCumulants_TPR.C:2076
Plot_PDCumulants_TPR.C:2077
Plot_PDCumulants_TPR.C:2078
Plot_PDCumulants_TPR.C:2079
Plot_PDCumulants_TPR.C:2080
Plot_PDCumulants_TPR.C:2081
Plot_PDCumulants_TPR.C:2082
Plot_PDCumulants_TPR.C:2083
Plot_PDCumulants_TPR.C:2084
Plot_PDCumulants_TPR.C:2085
Plot_PDCumulants_TPR.C:2086
Plot_PDCumulants_TPR.C:2087
Plot_PDCumulants_TPR.C:2088
Plot_PDCumulants_TPR.C:2089
Plot_PDCumulants_TPR.C:2090
Plot_PDCumulants_TPR.C:2091
Plot_PDCumulants_TPR.C:2092
Plot_PDCumulants_TPR.C:2093
Plot_PDCumulants_TPR.C:2094
Plot_PDCumulants_TPR.C:2095
Plot_PDCumulants_TPR.C:2096
Plot_PDCumulants_TPR.C:2097
Plot_PDCumulants_TPR.C:2098
Plot_PDCumulants_TPR.C:2099
Plot_PDCumulants_TPR.C:2100
Plot_PDCumulants_TPR.C:2101
Plot_PDCumulants_TPR.C:2102
Plot_PDCumulants_TPR.C:2103
Plot_PDCumulants_TPR.C:2104
Plot_PDCumulants_TPR.C:2105
Plot_PDCumulants_TPR.C:2106
Plot_PDCumulants_TPR.C:2107
Plot_PDCumulants_TPR.C:2108
Plot_PDCumulants_TPR.C:2109
Plot_PDCumulants_TPR.C:2110
Plot_PDCumulants_TPR.C:2111
Plot_PDCumulants_TPR.C:2112
Plot_PDCumulants_TPR.C:2113
Plot_PDCumulants_TPR.C:2114
Plot_PDCumulants_TPR.C:2115
Plot_PDCumulants_TPR.C:2116
Plot_PDCumulants_TPR.C:2117
Plot_PDCumulants_TPR.C:2118
Plot_PDCumulants_TPR.C:2119
Plot_PDCumulants_TPR.C:2120
Plot_PDCumulants_TPR.C:2121
Plot_PDCumulants_TPR.C:2122
Plot_PDCumulants_TPR.C:2123
Plot_PDCumulants_TPR.C:2124
Plot_PDCumulants_TPR.C:2125
Plot_PDCumulants_TPR.C:2126
Plot_PDCumulants_TPR.C:2127
Plot_PDCumulants_TPR.C:2128
Plot_PDCumulants_TPR.C:2129
Plot_PDCumulants_TPR.C:2130
Plot_PDCumulants_TPR.C:2131
Plot_PDCumulants_TPR.C:2132
Plot_PDCumulants_TPR.C:2133
Plot_PDCumulants_TPR.C:2134
Plot_PDCumulants_TPR.C:2135
Plot_PDCumulants_TPR.C:2136
Plot_PDCumulants_TPR.C:2137
Plot_PDCumulants_TPR.C:2138
Plot_PDCumulants_TPR.C:2139
Plot_PDCumulants_TPR.C:2140
Plot_PDCumulants_TPR.C:2141
Plot_PDCumulants_TPR.C:2142
Plot_PDCumulants_TPR.C:2143
Plot_PDCumulants_TPR.C:2144
Plot_PDCumulants_TPR.C:2145
Plot_PDCumulants_TPR.C:2146
Plot_PDCumulants_TPR.C:2147
Plot_PDCumulants_TPR.C:2148
Plot_PDCumulants_TPR.C:2149
Plot_PDCumulants_TPR.C:2150
Plot_PDCumulants_TPR.C:2151
Plot_PDCumulants_TPR.C:2152
Plot_PDCumulants_TPR.C:2153
Plot_PDCumulants_TPR.C:2154
Plot_PDCumulants_TPR.C:2155
Plot_PDCumulants_TPR.C:2156
Plot_PDCumulants_TPR.C:2157
Plot_PDCumulants_TPR.C:2158
Plot_PDCumulants_TPR.C:2159
Plot_PDCumulants_TPR.C:2160
Plot_PDCumulants_TPR.C:2161
Plot_PDCumulants_TPR.C:2162
Plot_PDCumulants_TPR.C:2163
Plot_PDCumulants_TPR.C:2164
Plot_PDCumulants_TPR.C:2165
Plot_PDCumulants_TPR.C:2166
Plot_PDCumulants_TPR.C:2167
Plot_PDCumulants_TPR.C:2168
Plot_PDCumulants_TPR.C:2169
Plot_PDCumulants_TPR.C:2170
Plot_PDCumulants_TPR.C:2171
Plot_PDCumulants_TPR.C:2172
Plot_PDCumulants_TPR.C:2173
Plot_PDCumulants_TPR.C:2174
Plot_PDCumulants_TPR.C:2175
Plot_PDCumulants_TPR.C:2176
Plot_PDCumulants_TPR.C:2177
Plot_PDCumulants_TPR.C:2178
Plot_PDCumulants_TPR.C:2179
Plot_PDCumulants_TPR.C:2180
Plot_PDCumulants_TPR.C:2181
Plot_PDCumulants_TPR.C:2182
Plot_PDCumulants_TPR.C:2183
Plot_PDCumulants_TPR.C:2184
Plot_PDCumulants_TPR.C:2185
Plot_PDCumulants_TPR.C:2186
Plot_PDCumulants_TPR.C:2187
Plot_PDCumulants_TPR.C:2188
Plot_PDCumulants_TPR.C:2189
Plot_PDCumulants_TPR.C:2190
Plot_PDCumulants_TPR.C:2191
Plot_PDCumulants_TPR.C:2192
Plot_PDCumulants_TPR.C:2193
Plot_PDCumulants_TPR.C:2194
Plot_PDCumulants_TPR.C:2195
Plot_PDCumulants_TPR.C:2196
Plot_PDCumulants_TPR.C:2197
Plot_PDCumulants_TPR.C:2198
Plot_PDCumulants_TPR.C:2199
Plot_PDCumulants_TPR.C:2200
Plot_PDCumulants_TPR.C:2201
Plot_PDCumulants_TPR.C:2202
Plot_PDCumulants_TPR.C:2203
Plot_PDCumulants_TPR.C:2204
Plot_PDCumulants_TPR.C:2205
Plot_PDCumulants_TPR.C:2206
Plot_PDCumulants_TPR.C:2207
Plot_PDCumulants_TPR.C:2208
Plot_PDCumulants_TPR.C:2209
Plot_PDCumulants_TPR.C:2210
Plot_PDCumulants_TPR.C:2211
Plot_PDCumulants_TPR.C:2212
Plot_PDCumulants_TPR.C:2213
Plot_PDCumulants_TPR.C:2214
Plot_PDCumulants_TPR.C:2215
Plot_PDCumulants_TPR.C:2216
Plot_PDCumulants_TPR.C:2217
Plot_PDCumulants_TPR.C:2218
Plot_PDCumulants_TPR.C:2219
Plot_PDCumulants_TPR.C:2220
Plot_PDCumulants_TPR.C:2221
Plot_PDCumulants_TPR.C:2222
Plot_PDCumulants_TPR.C:2223
Plot_PDCumulants_TPR.C:2224
Plot_PDCumulants_TPR.C:2225
Plot_PDCumulants_TPR.C:2226
Plot_PDCumulants_TPR.C:2227
Plot_PDCumulants_TPR.C:2228
Plot_PDCumulants_TPR.C:2229
Plot_PDCumulants_TPR.C:2230
Plot_PDCumulants_TPR.C:2231
Plot_PDCumulants_TPR.C:2232
Plot_PDCumulants_TPR.C:2233
Plot_PDCumulants_TPR.C:2234
Plot_PDCumulants_TPR.C:2235
Plot_PDCumulants_TPR.C:2236
Plot_PDCumulants_TPR.C:2237
Plot_PDCumulants_TPR.C:2238
Plot_PDCumulants_TPR.C:2239
Plot_PDCumulants_TPR.C:2240
Plot_PDCumulants_TPR.C:2241
Plot_PDCumulants_TPR.C:2242
Plot_PDCumulants_TPR.C:2243
Plot_PDCumulants_TPR.C:2244
Plot_PDCumulants_TPR.C:2245
Plot_PDCumulants_TPR.C:2246
Plot_PDCumulants_TPR.C:2247
Plot_PDCumulants_TPR.C:2248
Plot_PDCumulants_TPR.C:2249
Plot_PDCumulants_TPR.C:2250
Plot_PDCumulants_TPR.C:2251
Plot_PDCumulants_TPR.C:2252
Plot_PDCumulants_TPR.C:2253
Plot_PDCumulants_TPR.C:2254
Plot_PDCumulants_TPR.C:2255
Plot_PDCumulants_TPR.C:2256
Plot_PDCumulants_TPR.C:2257
Plot_PDCumulants_TPR.C:2258
Plot_PDCumulants_TPR.C:2259
Plot_PDCumulants_TPR.C:2260
Plot_PDCumulants_TPR.C:2261
Plot_PDCumulants_TPR.C:2262
Plot_PDCumulants_TPR.C:2263
Plot_PDCumulants_TPR.C:2264
Plot_PDCumulants_TPR.C:2265
Plot_PDCumulants_TPR.C:2266
Plot_PDCumulants_TPR.C:2267
Plot_PDCumulants_TPR.C:2268
Plot_PDCumulants_TPR.C:2269
Plot_PDCumulants_TPR.C:2270
Plot_PDCumulants_TPR.C:2271
Plot_PDCumulants_TPR.C:2272
Plot_PDCumulants_TPR.C:2273
Plot_PDCumulants_TPR.C:2274
Plot_PDCumulants_TPR.C:2275
Plot_PDCumulants_TPR.C:2276
Plot_PDCumulants_TPR.C:2277
Plot_PDCumulants_TPR.C:2278
Plot_PDCumulants_TPR.C:2279
Plot_PDCumulants_TPR.C:2280
Plot_PDCumulants_TPR.C:2281
Plot_PDCumulants_TPR.C:2282
Plot_PDCumulants_TPR.C:2283
Plot_PDCumulants_TPR.C:2284
Plot_PDCumulants_TPR.C:2285
Plot_PDCumulants_TPR.C:2286
Plot_PDCumulants_TPR.C:2287
Plot_PDCumulants_TPR.C:2288
Plot_PDCumulants_TPR.C:2289
Plot_PDCumulants_TPR.C:2290
Plot_PDCumulants_TPR.C:2291
Plot_PDCumulants_TPR.C:2292
Plot_PDCumulants_TPR.C:2293
Plot_PDCumulants_TPR.C:2294
Plot_PDCumulants_TPR.C:2295
Plot_PDCumulants_TPR.C:2296
Plot_PDCumulants_TPR.C:2297
Plot_PDCumulants_TPR.C:2298
Plot_PDCumulants_TPR.C:2299
Plot_PDCumulants_TPR.C:2300
Plot_PDCumulants_TPR.C:2301
Plot_PDCumulants_TPR.C:2302
Plot_PDCumulants_TPR.C:2303
Plot_PDCumulants_TPR.C:2304
Plot_PDCumulants_TPR.C:2305
Plot_PDCumulants_TPR.C:2306
Plot_PDCumulants_TPR.C:2307
Plot_PDCumulants_TPR.C:2308
Plot_PDCumulants_TPR.C:2309
Plot_PDCumulants_TPR.C:2310
Plot_PDCumulants_TPR.C:2311
Plot_PDCumulants_TPR.C:2312
Plot_PDCumulants_TPR.C:2313
Plot_PDCumulants_TPR.C:2314
Plot_PDCumulants_TPR.C:2315
Plot_PDCumulants_TPR.C:2316
Plot_PDCumulants_TPR.C:2317
Plot_PDCumulants_TPR.C:2318
Plot_PDCumulants_TPR.C:2319
Plot_PDCumulants_TPR.C:2320
Plot_PDCumulants_TPR.C:2321
Plot_PDCumulants_TPR.C:2322
Plot_PDCumulants_TPR.C:2323
Plot_PDCumulants_TPR.C:2324
Plot_PDCumulants_TPR.C:2325
Plot_PDCumulants_TPR.C:2326
Plot_PDCumulants_TPR.C:2327
Plot_PDCumulants_TPR.C:2328
Plot_PDCumulants_TPR.C:2329
Plot_PDCumulants_TPR.C:2330
Plot_PDCumulants_TPR.C:2331
Plot_PDCumulants_TPR.C:2332
Plot_PDCumulants_TPR.C:2333
Plot_PDCumulants_TPR.C:2334
Plot_PDCumulants_TPR.C:2335
Plot_PDCumulants_TPR.C:2336
Plot_PDCumulants_TPR.C:2337
Plot_PDCumulants_TPR.C:2338
Plot_PDCumulants_TPR.C:2339
Plot_PDCumulants_TPR.C:2340
Plot_PDCumulants_TPR.C:2341
Plot_PDCumulants_TPR.C:2342
Plot_PDCumulants_TPR.C:2343
Plot_PDCumulants_TPR.C:2344
Plot_PDCumulants_TPR.C:2345
Plot_PDCumulants_TPR.C:2346
Plot_PDCumulants_TPR.C:2347
Plot_PDCumulants_TPR.C:2348
Plot_PDCumulants_TPR.C:2349
Plot_PDCumulants_TPR.C:2350
Plot_PDCumulants_TPR.C:2351
Plot_PDCumulants_TPR.C:2352
Plot_PDCumulants_TPR.C:2353
Plot_PDCumulants_TPR.C:2354
Plot_PDCumulants_TPR.C:2355
Plot_PDCumulants_TPR.C:2356
Plot_PDCumulants_TPR.C:2357
Plot_PDCumulants_TPR.C:2358
Plot_PDCumulants_TPR.C:2359
Plot_PDCumulants_TPR.C:2360
Plot_PDCumulants_TPR.C:2361
Plot_PDCumulants_TPR.C:2362
Plot_PDCumulants_TPR.C:2363
Plot_PDCumulants_TPR.C:2364
Plot_PDCumulants_TPR.C:2365
Plot_PDCumulants_TPR.C:2366
Plot_PDCumulants_TPR.C:2367
Plot_PDCumulants_TPR.C:2368
Plot_PDCumulants_TPR.C:2369
Plot_PDCumulants_TPR.C:2370
Plot_PDCumulants_TPR.C:2371
Plot_PDCumulants_TPR.C:2372
Plot_PDCumulants_TPR.C:2373
Plot_PDCumulants_TPR.C:2374
Plot_PDCumulants_TPR.C:2375
Plot_PDCumulants_TPR.C:2376
Plot_PDCumulants_TPR.C:2377
Plot_PDCumulants_TPR.C:2378
Plot_PDCumulants_TPR.C:2379
Plot_PDCumulants_TPR.C:2380
Plot_PDCumulants_TPR.C:2381
Plot_PDCumulants_TPR.C:2382
Plot_PDCumulants_TPR.C:2383
Plot_PDCumulants_TPR.C:2384
Plot_PDCumulants_TPR.C:2385
Plot_PDCumulants_TPR.C:2386
Plot_PDCumulants_TPR.C:2387
Plot_PDCumulants_TPR.C:2388
Plot_PDCumulants_TPR.C:2389
Plot_PDCumulants_TPR.C:2390
Plot_PDCumulants_TPR.C:2391
Plot_PDCumulants_TPR.C:2392
Plot_PDCumulants_TPR.C:2393
Plot_PDCumulants_TPR.C:2394
Plot_PDCumulants_TPR.C:2395
Plot_PDCumulants_TPR.C:2396
Plot_PDCumulants_TPR.C:2397
Plot_PDCumulants_TPR.C:2398
Plot_PDCumulants_TPR.C:2399
Plot_PDCumulants_TPR.C:2400
Plot_PDCumulants_TPR.C:2401
Plot_PDCumulants_TPR.C:2402
Plot_PDCumulants_TPR.C:2403
Plot_PDCumulants_TPR.C:2404
Plot_PDCumulants_TPR.C:2405
Plot_PDCumulants_TPR.C:2406
Plot_PDCumulants_TPR.C:2407
Plot_PDCumulants_TPR.C:2408
Plot_PDCumulants_TPR.C:2409
Plot_PDCumulants_TPR.C:2410
Plot_PDCumulants_TPR.C:2411
Plot_PDCumulants_TPR.C:2412
Plot_PDCumulants_TPR.C:2413
Plot_PDCumulants_TPR.C:2414
Plot_PDCumulants_TPR.C:2415
Plot_PDCumulants_TPR.C:2416
Plot_PDCumulants_TPR.C:2417
Plot_PDCumulants_TPR.C:2418
Plot_PDCumulants_TPR.C:2419
Plot_PDCumulants_TPR.C:2420
Plot_PDCumulants_TPR.C:2421
Plot_PDCumulants_TPR.C:2422
Plot_PDCumulants_TPR.C:2423
Plot_PDCumulants_TPR.C:2424
Plot_PDCumulants_TPR.C:2425
Plot_PDCumulants_TPR.C:2426
Plot_PDCumulants_TPR.C:2427
Plot_PDCumulants_TPR.C:2428
Plot_PDCumulants_TPR.C:2429
Plot_PDCumulants_TPR.C:2430
Plot_PDCumulants_TPR.C:2431
Plot_PDCumulants_TPR.C:2432
Plot_PDCumulants_TPR.C:2433
Plot_PDCumulants_TPR.C:2434
Plot_PDCumulants_TPR.C:2435
Plot_PDCumulants_TPR.C:2436
Plot_PDCumulants_TPR.C:2437
Plot_PDCumulants_TPR.C:2438
Plot_PDCumulants_TPR.C:2439
Plot_PDCumulants_TPR.C:2440
Plot_PDCumulants_TPR.C:2441
Plot_PDCumulants_TPR.C:2442
Plot_PDCumulants_TPR.C:2443
Plot_PDCumulants_TPR.C:2444
Plot_PDCumulants_TPR.C:2445
Plot_PDCumulants_TPR.C:2446
Plot_PDCumulants_TPR.C:2447
Plot_PDCumulants_TPR.C:2448
Plot_PDCumulants_TPR.C:2449
Plot_PDCumulants_TPR.C:2450
Plot_PDCumulants_TPR.C:2451
Plot_PDCumulants_TPR.C:2452
Plot_PDCumulants_TPR.C:2453
Plot_PDCumulants_TPR.C:2454
Plot_PDCumulants_TPR.C:2455
Plot_PDCumulants_TPR.C:2456
Plot_PDCumulants_TPR.C:2457
Plot_PDCumulants_TPR.C:2458
Plot_PDCumulants_TPR.C:2459
Plot_PDCumulants_TPR.C:2460
Plot_PDCumulants_TPR.C:2461
Plot_PDCumulants_TPR.C:2462
Plot_PDCumulants_TPR.C:2463
Plot_PDCumulants_TPR.C:2464
Plot_PDCumulants_TPR.C:2465
Plot_PDCumulants_TPR.C:2466
Plot_PDCumulants_TPR.C:2467
Plot_PDCumulants_TPR.C:2468
Plot_PDCumulants_TPR.C:2469
Plot_PDCumulants_TPR.C:2470
Plot_PDCumulants_TPR.C:2471
Plot_PDCumulants_TPR.C:2472
Plot_PDCumulants_TPR.C:2473
Plot_PDCumulants_TPR.C:2474
Plot_PDCumulants_TPR.C:2475
Plot_PDCumulants_TPR.C:2476
Plot_PDCumulants_TPR.C:2477
Plot_PDCumulants_TPR.C:2478
Plot_PDCumulants_TPR.C:2479
Plot_PDCumulants_TPR.C:2480
Plot_PDCumulants_TPR.C:2481
Plot_PDCumulants_TPR.C:2482
Plot_PDCumulants_TPR.C:2483
Plot_PDCumulants_TPR.C:2484
Plot_PDCumulants_TPR.C:2485
Plot_PDCumulants_TPR.C:2486
Plot_PDCumulants_TPR.C:2487
Plot_PDCumulants_TPR.C:2488
Plot_PDCumulants_TPR.C:2489
Plot_PDCumulants_TPR.C:2490
Plot_PDCumulants_TPR.C:2491
Plot_PDCumulants_TPR.C:2492
Plot_PDCumulants_TPR.C:2493
Plot_PDCumulants_TPR.C:2494
Plot_PDCumulants_TPR.C:2495
Plot_PDCumulants_TPR.C:2496
Plot_PDCumulants_TPR.C:2497
Plot_PDCumulants_TPR.C:2498
Plot_PDCumulants_TPR.C:2499
Plot_PDCumulants_TPR.C:2500
Plot_PDCumulants_TPR.C:2501
Plot_PDCumulants_TPR.C:2502
Plot_PDCumulants_TPR.C:2503
Plot_PDCumulants_TPR.C:2504
Plot_PDCumulants_TPR.C:2505
Plot_PDCumulants_TPR.C:2506
Plot_PDCumulants_TPR.C:2507
Plot_PDCumulants_TPR.C:2508
Plot_PDCumulants_TPR.C:2509
Plot_PDCumulants_TPR.C:2510
Plot_PDCumulants_TPR.C:2511
Plot_PDCumulants_TPR.C:2512
Plot_PDCumulants_TPR.C:2513
Plot_PDCumulants_TPR.C:2514
Plot_PDCumulants_TPR.C:2515
Plot_PDCumulants_TPR.C:2516
Plot_PDCumulants_TPR.C:2517
Plot_PDCumulants_TPR.C:2518
Plot_PDCumulants_TPR.C:2519
Plot_PDCumulants_TPR.C:2520
Plot_PDCumulants_TPR.C:2521
Plot_PDCumulants_TPR.C:2522
Plot_PDCumulants_TPR.C:2523
Plot_PDCumulants_TPR.C:2524
Plot_PDCumulants_TPR.C:2525
Plot_PDCumulants_TPR.C:2526
Plot_PDCumulants_TPR.C:2527
Plot_PDCumulants_TPR.C:2528
Plot_PDCumulants_TPR.C:2529
Plot_PDCumulants_TPR.C:2530
Plot_PDCumulants_TPR.C:2531
Plot_PDCumulants_TPR.C:2532
Plot_PDCumulants_TPR.C:2533
Plot_PDCumulants_TPR.C:2534
Plot_PDCumulants_TPR.C:2535
Plot_PDCumulants_TPR.C:2536
Plot_PDCumulants_TPR.C:2537
Plot_PDCumulants_TPR.C:2538
Plot_PDCumulants_TPR.C:2539
Plot_PDCumulants_TPR.C:2540
Plot_PDCumulants_TPR.C:2541
Plot_PDCumulants_TPR.C:2542
Plot_PDCumulants_TPR.C:2543
Plot_PDCumulants_TPR.C:2544
Plot_PDCumulants_TPR.C:2545
Plot_PDCumulants_TPR.C:2546
Plot_PDCumulants_TPR.C:2547
Plot_PDCumulants_TPR.C:2548
Plot_PDCumulants_TPR.C:2549
Plot_PDCumulants_TPR.C:2550
Plot_PDCumulants_TPR.C:2551
Plot_PDCumulants_TPR.C:2552
Plot_PDCumulants_TPR.C:2553
Plot_PDCumulants_TPR.C:2554
Plot_PDCumulants_TPR.C:2555
Plot_PDCumulants_TPR.C:2556
Plot_PDCumulants_TPR.C:2557
Plot_PDCumulants_TPR.C:2558
Plot_PDCumulants_TPR.C:2559
Plot_PDCumulants_TPR.C:2560
Plot_PDCumulants_TPR.C:2561
Plot_PDCumulants_TPR.C:2562
Plot_PDCumulants_TPR.C:2563
Plot_PDCumulants_TPR.C:2564
Plot_PDCumulants_TPR.C:2565
Plot_PDCumulants_TPR.C:2566
Plot_PDCumulants_TPR.C:2567
Plot_PDCumulants_TPR.C:2568
Plot_PDCumulants_TPR.C:2569
Plot_PDCumulants_TPR.C:2570
Plot_PDCumulants_TPR.C:2571
Plot_PDCumulants_TPR.C:2572
Plot_PDCumulants_TPR.C:2573
Plot_PDCumulants_TPR.C:2574
Plot_PDCumulants_TPR.C:2575
Plot_PDCumulants_TPR.C:2576
Plot_PDCumulants_TPR.C:2577
Plot_PDCumulants_TPR.C:2578
Plot_PDCumulants_TPR.C:2579
Plot_PDCumulants_TPR.C:2580
Plot_PDCumulants_TPR.C:2581
Plot_PDCumulants_TPR.C:2582
Plot_PDCumulants_TPR.C:2583
Plot_PDCumulants_TPR.C:2584
Plot_PDCumulants_TPR.C:2585
Plot_PDCumulants_TPR.C:2586
Plot_PDCumulants_TPR.C:2587
Plot_PDCumulants_TPR.C:2588
Plot_PDCumulants_TPR.C:2589
Plot_PDCumulants_TPR.C:2590
Plot_PDCumulants_TPR.C:2591
Plot_PDCumulants_TPR.C:2592
Plot_PDCumulants_TPR.C:2593
Plot_PDCumulants_TPR.C:2594
Plot_PDCumulants_TPR.C:2595
Plot_PDCumulants_TPR.C:2596
Plot_PDCumulants_TPR.C:2597
Plot_PDCumulants_TPR.C:2598
Plot_PDCumulants_TPR.C:2599
Plot_PDCumulants_TPR.C:2600