#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" #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) using namespace std; // 0(standard), 1(fcSq=0.65), 2(fcSq=0.75), 3(B minus), 4(B plus), // 5(Linear instead of Cubic,also fc^2=0.75), 6(10h), 7(MRC 10% increase), 8(MuonCorr 92% pion-pair purity) int FileSetting=0; // bool MCcase_def=0;// MC data? int CollisionType=0;// PbPb, pPb, pp // int Mbin_def=0;// 0-9: centrality bin in widths of 5% bool SameCharge_def=1;// same-charge? int CHARGE_def=1;// -1 or +1: + or - pions for same-charge case, --+ or -++, ---+ or -+++ int MixedCharge4pionType_def = 2;// 1(---+) or 2(--++) // int EDbin_def=0;// 0 or 1: Kt3 bin int TPNbin=0;// TPN bin for r3 and r4 int Gbin = int( (0) /2. ) + 5;// +5 (Rcoh=0), +55 (Rcoh=Rch) int c3FitType = 1;// EW(1), LG(2) int Ktbin_def=1;// 1(0.2-0.3),..., 6(0.7-0.8), 10 = Full Range // bool MRC=1;// Momentum Resolution Corrections? bool MuonCorrection=1;// correct for Muons? bool FSICorrection=1;// correct For Final-State-Interactions bool InterpCorrection=1;// correct for finite bin interpolation // int f_choice=0;// 0(Core/Halo), 1(40fm), 2(70fm), 3(100fm) // // bool SaveToFile_def=kFALSE;// Save outputs to file? bool GeneratedSignal=kFALSE;// Was the QS+FSI signal generated in MC? // // // // int fFSIindex=0; int Ktlowbin; int Kthighbin; float TwoFrac;// Lambda parameter float OneFrac;// Lambda^{1/2} float ThreeFrac;// Lambda^{3/2} float FourFrac;// lambda^{4/2} double Qcut_pp = 0.6;// 0.6 double Qcut_PbPb = 0.1; double NormQcutLow_pp = 1.0; double NormQcutHigh_pp = 1.5; double NormQcutLow_PbPb = .15; double NormQcutHigh_PbPb = .2;// was .175 int ThreeParticleRebin; int FourParticleRebin; int RbinMRC; TH1D *fFSIss[12]; TH1D *fFSIos[12]; void SetFSICorrelations(); void SetFSIindex(Float_t); Float_t FSICorrelation(Int_t, Int_t, Float_t); void SetMuonCorrections(); void SetMomResCorrections(); double Gamov(int, double); double cubicInterpolate(double[4], double); double bicubicInterpolate(double[4][4], double, double); double tricubicInterpolate(double[4][4][4], double, double, double); // float fact(float); TH1D *MRC_SC_2[2]; TH1D *MRC_MC_2[2]; TH1D *MRC_SC_3[3]; TH1D *MRC_MC_3[4]; TH1D *MRC_SC_4[5]; TH1D *MRC_MC1_4[6]; TH1D *MRC_MC2_4[6]; double AvgQinvSS[30]; double AvgQinvOS[30]; double BinCentersQ4[20]; // TH1D *C2muonCorrectionSC; TH1D *C2muonCorrectionMC; TH1D *WeightmuonCorrection; TH1D *C3muonCorrectionSC[2]; TH1D *C3muonCorrectionMC[3]; TH1D *C4muonCorrectionSC[4]; TH1D *C4muonCorrectionMC1[5]; TH1D *C4muonCorrectionMC2[5]; void Plot_FourPion(bool SaveToFile=SaveToFile_def, bool MCcase=MCcase_def, bool SameCharge=SameCharge_def, int MixedCharge4pionType=MixedCharge4pionType_def, int EDbin=EDbin_def, int CHARGE=CHARGE_def, int Mbin=Mbin_def, int Ktbin=Ktbin_def){ EDbin_def=EDbin; SaveToFile_def=SaveToFile; MCcase_def=MCcase; CHARGE_def=CHARGE; SameCharge_def=SameCharge;// 3-pion same-charge MixedCharge4pionType_def=MixedCharge4pionType; Mbin_def=Mbin; Ktbin_def=Ktbin; // Ktlowbin=(Ktbin)*2+3;// kt bins are 0.5 GeV/c wide (0-0.5, 0.5-1.0 ...) Kthighbin=(Ktbin)*2+4; // if(FileSetting==1) TwoFrac=0.65; else if(FileSetting==2 || FileSetting==5) TwoFrac=0.75; else TwoFrac=0.7; OneFrac = sqrt(TwoFrac); ThreeFrac = pow(TwoFrac, 1.5); FourFrac = pow(TwoFrac, 2.); //// Core/Halo, 40fm, 70fm, 100fm float ThermShift_f33[4]={0., 0.06933, 0.01637, 0.006326}; float ThermShift_f32[4]={0., -0.0185, -0.005301, -0.002286}; float ThermShift_f31[4]={0., -0.01382, -0.0004682, 0.0005337}; float f33prime = ThreeFrac; float f32prime = TwoFrac*(1-OneFrac); float f31prime = pow(1-OneFrac,3) + 3*OneFrac*pow(1-OneFrac,2); f33prime += ThermShift_f33[f_choice]; f32prime += ThermShift_f32[f_choice]; f31prime += ThermShift_f31[f_choice]; float f33 = f33prime; float f32 = f32prime/TwoFrac; float f31 = f31prime - 3*((1-TwoFrac)*(1-OneFrac) + ThermShift_f32[f_choice]*(1-TwoFrac)/TwoFrac); //cout<<f33 + 3*f32 + f31<<endl; //// Core/Halo, 40fm, 70fm, 100fm float ThermShift_f44[4]={0., 0.08741, 0.005284, -0.01064}; float ThermShift_f43[4]={0., -0.01126, -0.001486, 0.001686}; float ThermShift_f42[4]={0., -0.006466, -7.683e-05, 0.0004572}; float ThermShift_f41[4]={0., -0.003556, 0.00112, 0.00115}; float f44prime = FourFrac; float f43prime = ThreeFrac*(1-OneFrac); float f42prime = TwoFrac*pow(1-OneFrac,2); float f41prime = pow(1-OneFrac,4) + 4*OneFrac*pow(1-OneFrac,3); f44prime += ThermShift_f44[f_choice]; f43prime += ThermShift_f43[f_choice]; f42prime += ThermShift_f42[f_choice]; f41prime += ThermShift_f41[f_choice]; float f44 = f44prime; float f43 = f43prime/f33prime; float f42 = f42prime/TwoFrac; f42 -= 2*f43prime*f32prime/f33prime/TwoFrac; float f41 = f41prime; f41 -= 4*f43prime*f31prime/f33prime; f41 -= 6*f42prime*(1-TwoFrac)/TwoFrac; f41 += 12*f43prime/f33prime*f32prime/TwoFrac*(1-TwoFrac); //cout<<f44 + 4*f43 + 6*f42 + f41<<endl; //cout<<f44<<" "<<f43<<" "<<f42<<" "<<f41<<endl; // // Core/Halo, 40fm, 70fm, 100fm //float TherminatorMod_f33[4]={1., 1.123, 1.022, 1.008};// 1.008 //float TherminatorMod_f32[4]={1., 0.844, 0.933, 0.967};// .967 //float TherminatorMod_f31[4]={1., 0.828, 0.982, 1.028};// 1.028 /*float TherminatorMod_f44[4]={1., 1.188, 1.008, 0.985};// .985 float TherminatorMod_f43[4]={1., 0.885, 0.979, 1.026};// 1.026 float TherminatorMod_f42[4]={1., 0.687, 0.99, 1.08};// 1.08 float TherminatorMod_f41[4]={1., 0.806, 1.33, 1.548};//1.548 float f44 = FourFrac; float f43 = (1-OneFrac); float f42 = -pow(1-OneFrac,2); float f41 = -3*pow(1-OneFrac,4) - 8*OneFrac*pow(1-OneFrac,3) + 6*(1-TwoFrac)*pow(1-OneFrac,2); f44 *= TherminatorMod_f44[f_choice]; f43 *= TherminatorMod_f43[f_choice]; f42 *= TherminatorMod_f42[f_choice]; f41 *= TherminatorMod_f41[f_choice];*/ //f44prime *= TherminatorMod_f44[f_choice]; //f43prime *= TherminatorMod_f43[f_choice]; //f42prime *= TherminatorMod_f42[f_choice]; //f41prime *= TherminatorMod_f41[f_choice]; //float f44 = f44prime; //float f43 = f43prime/ThreeFrac; //float f42 = f42prime/TwoFrac - 2*pow(1-OneFrac,2); //float f41 = f41prime - 4*pow(1-OneFrac,4) - 12*OneFrac*pow(1-OneFrac,3) + 6*(1-TwoFrac)*pow(1-OneFrac,2); cout<<"Mbin = "<<Mbin<<" KT3 = "<<EDbin<<" SameCharge = "<<SameCharge<<endl; if(!SameCharge) cout<<"4-pion MixedCharge type = "<<MixedCharge4pionType<<endl; if(CollisionType==0){ if(Mbin==0) {RbinMRC=10;} else if(Mbin==1) {RbinMRC=9;} else if(Mbin<=3) {RbinMRC=8;} else if(Mbin<=5) {RbinMRC=7;} else {RbinMRC=6;} }else{ RbinMRC=2; } if(MCcase) MuonCorrection=kFALSE; if(CollisionType==0){ ThreeParticleRebin=2; FourParticleRebin=3; }else{ ThreeParticleRebin=2; FourParticleRebin=6; } float Cutoff_FullWeight_Q3[10]={0.065, 0.065, 0.075, 0.075, 0.075, 0.075, 0.085, 0.085, 0.085, 0.085}; float Cutoff_FullWeight_Q4[10]={0.115, 0.115, 0.115, 0.130, 0.130, 0.130, 0.145, 0.145, 0.145, 0.145}; // // avg Qinv within the 5 MeV bins (0-5, 5-10,...) for true bin mean values. From unsmeared HIJING 0-5% with input signal double AvgQinvSS_temp[30]={0.00421164, 0.00796583, 0.0128473, 0.0178262, 0.0228416, 0.0276507, 0.0325368, 0.0376114, 0.0424707, 0.0476274, 0.0526015, 0.0575645, 0.0625478, 0.0675416, 0.0725359, 0.0775219, 0.0825521, 0.0875211, 0.0925303, 0.0975319, 0.102544, 0.10753, 0.112499, 0.117545, 0.122522, 0.127522, 0.132499, 0.137514, 0.142495, 0.147521}; double AvgQinvOS_temp[30]={0.00352789, 0.00797596, 0.0128895, 0.0177198, 0.0226397, 0.0276331, 0.0326309, 0.0376471, 0.0426217, 0.047567, 0.0525659, 0.0575472, 0.0625886, 0.0675261, 0.0725543, 0.077564, 0.0825617, 0.0875465, 0.092539, 0.0975348, 0.102529, 0.107527, 0.112506, 0.117531, 0.122536, 0.1275, 0.132508, 0.137508, 0.14251, 0.147511}; for(int i=0; i<30; i++){ AvgQinvSS[i] = AvgQinvSS_temp[i]; AvgQinvOS[i] = AvgQinvOS_temp[i]; } TH2D *TwoParticle_2d[2][2][2];// ch1,ch2,term TH1D *TwoParticle[2][2][2];// ch1,ch2,term TH2D *UnitMult_2d[2][2][2];// ch1,ch2,term TH1D *UnitMult[2][2][2];// ch1,ch2,term double norm_2[2]={0}; double norm_2_UM[2]={0}; // TH1D *ThreeParticle[2][2][2][5];// ch1,ch2,ch3,term TProfile *K3avg[2][2][2][4]; TH2D *TPFullWeight_ThreeParticle_2D[2];// charge TH2D *TPNegFullWeight_ThreeParticle_2D[2];// charge TH1D *TPN_ThreeParticle[2];// charge TH1D *TPFullWeight_ThreeParticle[2];// charge TH1D *TPNegFullWeight_ThreeParticle[2];// charge double norm_3[5]={0}; // TH1D *FourParticle[2][2][2][2][13];// ch1,ch2,ch3,ch4,term TProfile *K4avg[2][2][2][2][12]; TH2D *TPFullWeight_FourParticle_2D[2];// charge TH2D *TPNegFullWeight_FourParticle_2D[2];// charge TH1D *TPFullWeight_FourParticle[2];// charge TH1D *TPNegFullWeight_FourParticle[2];// charge TH1D *TPN_FourParticle[2];// charge TH3D *FullBuildFromFits_3D;// charge TH1D *FullBuildFromFits;// charge TH3D *PartialBuildFromFits_3D;// charge TH1D *PartialBuildFromFits;// charge double norm_4[13]={0}; 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_GeneratedSignal.root","READ"); //_file0 = new TFile("Results/PDC_12a17a_pTSpectrumWeight.root","READ"); _file0 = new TFile("Results/PDC_12a17a_Qweights.root","READ"); } }else{ //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_pT_0p2to1p0_FullRunWrongWeightsNoPadRowTTCandMCTTC_RawWeightFile.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_pT_0p16to0p25_0p25to1p0.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_pT_0p2to1p0.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_pT0p16to0p25_KT0p35.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_1percentCentEM.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_0p02eta0p045phi_0p03eta0p067phi.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_c3FitBuild.root","READ"); if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_LowNorm_HighNorm.root","READ"); //if(FileSetting==0) _file0 = new TFile("Results/PDC_11h_extendedGweights.root","READ");// Preliminary results if(FileSetting==5) _file0 = new TFile("Results/PDC_11h_Cubic_Linear.root","READ"); if(FileSetting==1 || FileSetting==2) _file0 = new TFile("Results/PDC_11h_Lam0p65_Lam0p75.root","READ"); if(FileSetting==3) _file0 = new TFile("Results/PDC_11h_Cubic_Linear_Bminus.root","READ"); if(FileSetting==4) _file0 = new TFile("Results/PDC_11h_Cubic_Linear_Bplus.root","READ"); if(FileSetting==6) _file0 = new TFile("Results/PDC_10h_Cubic_Linear.root","READ"); if(FileSetting==7 || FileSetting==8) _file0 = new TFile("Results/PDC_11h_MRC10percIncrease_Muon92percent.root","READ"); } }else if(CollisionType==1){// pPb _file0 = new TFile("Results/PDC_13bc_kINT7_LowNorm_HighNorm.root","READ"); }else{// pp _file0 = new TFile("Results/PDC_10bcde_kMB_LowNorm_HighNorm.root","READ"); } SetFSIindex(10.); SetFSICorrelations(); SetMomResCorrections(); SetMuonCorrections(); // ///////////////////////////////////////////////////////// double NormQcutLow; double NormQcutHigh; if(CollisionType==0) { NormQcutLow = NormQcutLow_PbPb; NormQcutHigh = NormQcutHigh_PbPb; }else { NormQcutLow = NormQcutLow_pp; NormQcutHigh = NormQcutHigh_pp; } TList *MyList; if(!MCcase){ TDirectoryFile *tdir = (TDirectoryFile*)_file0->Get("PWGCF.outputFourPionAnalysis.root"); if(FileSetting==2 || FileSetting==5 || FileSetting==8 || CollisionType!=0) MyList=(TList*)tdir->Get("FourPionOutput_2"); else MyList=(TList*)tdir->Get("FourPionOutput_1"); //MyList=(TList*)_file0->Get("MyList"); }else{ MyList=(TList*)_file0->Get("MyList"); } _file0->Close(); TH1D *Events = (TH1D*)MyList->FindObject("fEvents2"); // cout<<"#Events = "<<Events->Integral(Mbin+1,Mbin+1)<<endl; /////////////////////////////////// // Get Histograms // 2-particle for(int c1=0; c1<2; c1++){ for(int c2=0; c2<2; c2++){ for(int term=0; term<2; term++){ TString *name2 = new TString("TwoParticle_Charge1_"); *name2 += c1; name2->Append("_Charge2_"); *name2 += c2; name2->Append("_M_"); *name2 += Mbin; name2->Append("_ED_"); *name2 += 0;// EDbin name2->Append("_Term_"); *name2 += term+1; if( (c1+c2)==1 ) {if(c1!=0) continue;}// skip degenerate histogram TwoParticle_2d[c1][c2][term] = (TH2D*)MyList->FindObject(name2->Data()); TwoParticle_2d[c1][c2][term]->Sumw2(); TString *proname = new TString(name2->Data()); proname->Append("_pro"); if(Ktbin==10) {Ktlowbin=1; Kthighbin=TwoParticle_2d[c1][c2][term]->GetNbinsX();} TwoParticle[c1][c2][term] = (TH1D*)TwoParticle_2d[c1][c2][term]->ProjectionY(proname->Data(),Ktlowbin,Kthighbin);// bins:5-6 (kt:0.2-0.3) norm_2[term] = TwoParticle[c1][c2][term]->Integral(TwoParticle[c1][c2][term]->GetXaxis()->FindBin(NormQcutLow),TwoParticle[c1][c2][term]->GetXaxis()->FindBin(NormQcutHigh)); //cout<<"2-pion norms "<<norm_2[term]<<endl; TwoParticle[c1][c2][term]->Scale(norm_2[0]/norm_2[term]); TwoParticle[c1][c2][term]->SetMarkerStyle(20); TwoParticle[c1][c2][term]->GetXaxis()->SetTitle("q_{inv} (GeV/c)"); TwoParticle[c1][c2][term]->GetYaxis()->SetTitle("C_{2}"); TwoParticle[c1][c2][term]->SetTitle(""); // TString *nameUM=new TString(name2->Data()); nameUM->Append("_UnitMult"); UnitMult_2d[c1][c2][term] = (TH2D*)MyList->FindObject(nameUM->Data()); UnitMult_2d[c1][c2][term]->Sumw2(); TString *pronameUM = new TString(nameUM->Data()); pronameUM->Append("_pro"); UnitMult[c1][c2][term] = (TH1D*)UnitMult_2d[c1][c2][term]->ProjectionY(pronameUM->Data(),13,13);// 11 means 1000 pions norm_2_UM[term] = UnitMult[c1][c2][term]->Integral(UnitMult[c1][c2][term]->GetXaxis()->FindBin(1.0),UnitMult[c1][c2][term]->GetXaxis()->FindBin(1.2)); UnitMult[c1][c2][term]->Scale(norm_2_UM[0]/norm_2_UM[term]); // UnitMult[c1][c2][term]->SetMarkerStyle(20); UnitMult[c1][c2][term]->GetXaxis()->SetTitle("q_{inv} (GeV/c)"); UnitMult[c1][c2][term]->GetYaxis()->SetTitle("C_{2}"); UnitMult[c1][c2][term]->SetTitle(""); }// term // 3-particle for(int c3=0; c3<2; c3++){ for(int term=0; term<5; term++){ TString *name3 = new TString("ThreeParticle_Charge1_"); *name3 += c1; name3->Append("_Charge2_"); *name3 += c2; name3->Append("_Charge3_"); *name3 += c3; name3->Append("_M_"); *name3 += Mbin; name3->Append("_ED_"); *name3 += EDbin; name3->Append("_Term_"); *name3 += term+1; TString *nameNorm3=new TString(name3->Data()); nameNorm3->Append("_Norm"); // TString *nameK3=new TString(name3->Data()); nameK3->Append("_Kfactor"); // TString *nameTPN3=new TString(name3->Data()); nameTPN3->Append("_TwoPartNorm"); // TString *nameNegTPN3=new TString(name3->Data()); nameNegTPN3->Append("_TwoPartNegNorm"); // name3->Append("_1D"); /////////////////////////////////////// // skip degenerate histograms if( (c1+c2+c3)==1) {if(c1!=0 || c2!=0 || c3!=1) continue;} if( (c1+c2+c3)==2) {if(c1!=0) continue;} //////////////////////////////////////// norm_3[term] = ((TH1D*)MyList->FindObject(nameNorm3->Data()))->Integral(); ThreeParticle[c1][c2][c3][term] = (TH1D*)MyList->FindObject(name3->Data()); ThreeParticle[c1][c2][c3][term]->Sumw2(); //if(c1==c2 && c1==c3) cout<<"3-pion norms "<<norm_3[term]<<endl; ThreeParticle[c1][c2][c3][term]->Scale(norm_3[0]/norm_3[term]); ThreeParticle[c1][c2][c3][term]->GetXaxis()->SetTitle("Q_{3} (GeV/c)"); ThreeParticle[c1][c2][c3][term]->GetYaxis()->SetTitle("C_{3}"); ThreeParticle[c1][c2][c3][term]->SetMarkerStyle(20); ThreeParticle[c1][c2][c3][term]->SetTitle(""); // ThreeParticle[c1][c2][c3][term]->Rebin(ThreeParticleRebin); // if(term<4){ K3avg[c1][c2][c3][term] = (TProfile*)MyList->FindObject(nameK3->Data()); K3avg[c1][c2][c3][term]->Rebin(ThreeParticleRebin); if(MCcase || !FSICorrection) { K3avg[c1][c2][c3][term]->Reset(); for(int ii=1; ii<=K3avg[c1][c2][c3][term]->GetNbinsX(); ii++) { K3avg[c1][c2][c3][term]->Fill(K3avg[c1][c2][c3][term]->GetXaxis()->GetBinCenter(ii), 1); } } } // if(term==4 && c1==c2 && c1==c3){ TPFullWeight_ThreeParticle_2D[c1] = (TH2D*)MyList->FindObject(nameTPN3->Data()); TPFullWeight_ThreeParticle_2D[c1]->Scale(norm_3[0]/norm_3[term]); TPFullWeight_ThreeParticle_2D[c1]->RebinY(ThreeParticleRebin); // TPNegFullWeight_ThreeParticle_2D[c1] = (TH2D*)MyList->FindObject(nameNegTPN3->Data()); TPNegFullWeight_ThreeParticle_2D[c1]->Scale(norm_3[0]/norm_3[term]); TPNegFullWeight_ThreeParticle_2D[c1]->RebinY(ThreeParticleRebin); // for(int binx=1; binx<=TPFullWeight_ThreeParticle_2D[c1]->GetNbinsX(); binx++) { for(int biny=1; biny<=TPFullWeight_ThreeParticle_2D[c1]->GetNbinsY(); biny++) { TPFullWeight_ThreeParticle_2D[c1]->SetBinError(binx, biny, 0); if(binx!=4){ TPFullWeight_ThreeParticle_2D[c1]->SetBinContent(binx, biny, TPFullWeight_ThreeParticle_2D[c1]->GetBinContent(binx, biny) + TPFullWeight_ThreeParticle_2D[c1]->GetBinContent(4, biny)); TPNegFullWeight_ThreeParticle_2D[c1]->SetBinContent(binx, biny, TPNegFullWeight_ThreeParticle_2D[c1]->GetBinContent(binx, biny) + TPNegFullWeight_ThreeParticle_2D[c1]->GetBinContent(4, biny)); if(InterpCorrection){ double q3 = TPFullWeight_ThreeParticle_2D[c1]->GetYaxis()->GetBinCenter(biny); double InterCorr = 1.024 - 0.2765*q3; //if(q3<0.1) cout<<binx<<" "<<biny<<" "<<q3<<" "<<InterCorr<<endl; TPFullWeight_ThreeParticle_2D[c1]->SetBinContent(binx, biny, InterCorr * TPFullWeight_ThreeParticle_2D[c1]->GetBinContent(binx, biny)); TPNegFullWeight_ThreeParticle_2D[c1]->SetBinContent(binx, biny, InterCorr * TPNegFullWeight_ThreeParticle_2D[c1]->GetBinContent(binx, biny)); } } } } TString *proName=new TString(nameTPN3->Data()); TString *proNameNeg=new TString(nameNegTPN3->Data()); proName->Append("_pro"); proNameNeg->Append("_pro"); TPN_ThreeParticle[c1] = (TH1D*)TPFullWeight_ThreeParticle_2D[c1]->ProjectionY(proName->Data(), TPNbin, TPNbin); // proName->Append("_FullWeight"); proNameNeg->Append("_FullWeight"); TPFullWeight_ThreeParticle[c1] = (TH1D*)TPFullWeight_ThreeParticle_2D[c1]->ProjectionY(proName->Data(), Gbin, Gbin); TPNegFullWeight_ThreeParticle[c1] = (TH1D*)TPNegFullWeight_ThreeParticle_2D[c1]->ProjectionY(proNameNeg->Data(), Gbin, Gbin); proName->Append("_FullWeightDen"); proNameNeg->Append("_FullWeightDen"); TH1D *tempDen = (TH1D*)TPFullWeight_ThreeParticle_2D[c1]->ProjectionY(proName->Data(), 4, 4); TH1D *tempDenNeg = (TH1D*)TPNegFullWeight_ThreeParticle_2D[c1]->ProjectionY(proNameNeg->Data(), 4, 4); // Add Pos with Neg weights tempDen->Add(tempDenNeg); TPFullWeight_ThreeParticle[c1]->Add(TPNegFullWeight_ThreeParticle[c1]); // //TPFullWeight_ThreeParticle[c1]->Add(tempDen);// now added above in Interp section TPFullWeight_ThreeParticle[c1]->Divide(tempDen); TPFullWeight_ThreeParticle[c1]->SetLineColor(1); // } }// term 3-particle // 4-particle for(int c4=0; c4<2; c4++){ for(int term=0; term<13; term++){ TString *name4 = new TString("FourParticle_Charge1_"); *name4 += c1; name4->Append("_Charge2_"); *name4 += c2; name4->Append("_Charge3_"); *name4 += c3; name4->Append("_Charge4_"); *name4 += c4; name4->Append("_M_"); *name4 += Mbin; name4->Append("_ED_"); *name4 += EDbin; name4->Append("_Term_"); *name4 += term+1; TString *nameNorm4=new TString(name4->Data()); nameNorm4->Append("_Norm"); // TString *nameK4=new TString(name4->Data()); nameK4->Append("_Kfactor");// or Weighted // TString *nameTPN4=new TString(name4->Data()); nameTPN4->Append("_TwoPartNorm"); // TString *nameNegTPN4=new TString(name4->Data()); nameNegTPN4->Append("_TwoPartNegNorm"); // TString *nameFitBuild=new TString(name4->Data()); nameFitBuild->Append("_FullBuildFromFits"); // TString *namePartialFitBuild=new TString(name4->Data()); namePartialFitBuild->Append("_PartialBuildFromFits"); // name4->Append("_1D"); /////////////////////////////////////// // skip degenerate histograms if( (c1+c2+c3+c4)==1) {if(c4!=1) continue;} if( (c1+c2+c3+c4)==2) {if(c3+c4!=2) continue;} if( (c1+c2+c3+c4)==3) {if(c1!=0) continue;} ///////////////////////////////////////// norm_4[term] = ((TH1D*)MyList->FindObject(nameNorm4->Data()))->Integral(); //if( (c1+c2+c3+c4)==4) cout<<"4-pion norms "<<norm_4[term]<<endl; if(norm_4[term] > 0){ //if(c1==c2 && c1==c3 && c1==c4) cout<<term<<" "<<norm_4[0]/norm_4[term]<<endl; FourParticle[c1][c2][c3][c4][term] = (TH1D*)MyList->FindObject(name4->Data()); FourParticle[c1][c2][c3][c4][term]->Sumw2(); FourParticle[c1][c2][c3][c4][term]->Scale(norm_4[0]/norm_4[term]); FourParticle[c1][c2][c3][c4][term]->GetXaxis()->SetTitle("Q_{4} (GeV/c)"); FourParticle[c1][c2][c3][c4][term]->GetYaxis()->SetTitle("C_{4}"); FourParticle[c1][c2][c3][c4][term]->SetMarkerStyle(20); FourParticle[c1][c2][c3][c4][term]->SetTitle(""); // FourParticle[c1][c2][c3][c4][term]->Rebin(FourParticleRebin); }else{ cout<<"4-particle normalization = 0."<<endl; } if(term<12){ K4avg[c1][c2][c3][c4][term] = (TProfile*)MyList->FindObject(nameK4->Data()); K4avg[c1][c2][c3][c4][term]->Rebin(FourParticleRebin); if(MCcase || !FSICorrection) { K4avg[c1][c2][c3][c4][term]->Reset(); for(int ii=1; ii<=K4avg[c1][c2][c3][c4][term]->GetNbinsX(); ii++) { K4avg[c1][c2][c3][c4][term]->Fill(K4avg[c1][c2][c3][c4][term]->GetXaxis()->GetBinCenter(ii), 1); } } } if(term==12 && c1==c2 && c1==c3 && c1==c4){ TPFullWeight_FourParticle_2D[c1] = (TH2D*)MyList->FindObject(nameTPN4->Data()); TPFullWeight_FourParticle_2D[c1]->Scale(norm_4[0]/norm_4[term]); TPFullWeight_FourParticle_2D[c1]->RebinY(FourParticleRebin); // TPNegFullWeight_FourParticle_2D[c1] = (TH2D*)MyList->FindObject(nameNegTPN4->Data()); TPNegFullWeight_FourParticle_2D[c1]->Scale(norm_4[0]/norm_4[term]); TPNegFullWeight_FourParticle_2D[c1]->RebinY(FourParticleRebin); // if(c1==0 && !MCcase){ FullBuildFromFits_3D = (TH3D*)MyList->FindObject(nameFitBuild->Data()); FullBuildFromFits_3D->Scale(norm_4[0]/norm_4[term]); FullBuildFromFits_3D->RebinZ(FourParticleRebin); // PartialBuildFromFits_3D = (TH3D*)MyList->FindObject(namePartialFitBuild->Data()); PartialBuildFromFits_3D->Scale(norm_4[0]/norm_4[term]); PartialBuildFromFits_3D->RebinZ(FourParticleRebin); } // for(int binx=1; binx<=TPFullWeight_FourParticle_2D[c1]->GetNbinsX(); binx++) { for(int biny=1; biny<=TPFullWeight_FourParticle_2D[c1]->GetNbinsY(); biny++) { TPFullWeight_FourParticle_2D[c1]->SetBinError(binx, biny, 0); if(binx!=4){ TPFullWeight_FourParticle_2D[c1]->SetBinContent(binx, biny, TPFullWeight_FourParticle_2D[c1]->GetBinContent(binx, biny) + TPFullWeight_FourParticle_2D[c1]->GetBinContent(4, biny)); TPNegFullWeight_FourParticle_2D[c1]->SetBinContent(binx, biny, TPNegFullWeight_FourParticle_2D[c1]->GetBinContent(binx, biny) + TPNegFullWeight_FourParticle_2D[c1]->GetBinContent(4, biny)); if(InterpCorrection){// Interpolator correction double q4 = TPFullWeight_FourParticle_2D[c1]->GetYaxis()->GetBinCenter(biny); double InterCorr = 1.032 - 0.2452*q4; TPFullWeight_FourParticle_2D[c1]->SetBinContent(binx, biny, InterCorr * TPFullWeight_FourParticle_2D[c1]->GetBinContent(binx, biny)); TPNegFullWeight_FourParticle_2D[c1]->SetBinContent(binx, biny, InterCorr * TPNegFullWeight_FourParticle_2D[c1]->GetBinContent(binx, biny)); } } } } TString *proName=new TString(nameTPN4->Data()); TString *proNegName=new TString(nameNegTPN4->Data()); TString *proNameFitBuild=new TString(nameFitBuild->Data()); TString *proNamePartialFitBuild=new TString(namePartialFitBuild->Data()); proName->Append("_pro"); proNegName->Append("_pro"); proNameFitBuild->Append("_pro"); proNamePartialFitBuild->Append("_pro"); TPN_FourParticle[c1] = (TH1D*)TPFullWeight_FourParticle_2D[c1]->ProjectionY(proName->Data(), TPNbin, TPNbin); // proName->Append("_FullWeight"); proNegName->Append("_FullWeight"); TPFullWeight_FourParticle[c1] = (TH1D*)TPFullWeight_FourParticle_2D[c1]->ProjectionY(proName->Data(), Gbin, Gbin); TPNegFullWeight_FourParticle[c1] = (TH1D*)TPNegFullWeight_FourParticle_2D[c1]->ProjectionY(proNegName->Data(), Gbin, Gbin); proName->Append("_FullWeightDen"); proNegName->Append("_FullWeightDen"); TH1D *tempDen = (TH1D*)TPFullWeight_FourParticle_2D[c1]->ProjectionY(proName->Data(), 4, 4); TH1D *tempDenNeg = (TH1D*)TPNegFullWeight_FourParticle_2D[c1]->ProjectionY(proNegName->Data(), 4, 4); // // Add Pos and Neg Weights tempDen->Add(tempDenNeg); TPFullWeight_FourParticle[c1]->Add(TPNegFullWeight_FourParticle[c1]); // //TPFullWeight_FourParticle[c1]->Add(tempDen);// now added above in Interp section TPFullWeight_FourParticle[c1]->Divide(tempDen); TPFullWeight_FourParticle[c1]->SetLineColor(1); /*TString *ErrName=new TString(nameTPN4->Data()); ErrName->Append("Err"); TH2D *temperr2D = (TH2D*)MyList->FindObject(ErrName->Data()); TH1D *temperr = (TH1D*)temperr2D->ProjectionY("tesst",4,4); temperr->Rebin(FourParticleRebin); cout.precision(8); cout<<temperr->GetBinContent(3)<<endl; cout<<(temperr->GetBinContent(5) / tempDen->GetBinContent(5))<<" "<<TPFullWeight_FourParticle[c1]->GetBinContent(5)<<endl; */ if(c1==0 && !MCcase){ FullBuildFromFits = (TH1D*)FullBuildFromFits_3D->ProjectionZ(proNameFitBuild->Data(), c3FitType, c3FitType, Gbin, Gbin); TH1D *tempDen2 = (TH1D*)FullBuildFromFits_3D->ProjectionZ("tempDen2", c3FitType, c3FitType, 4, 4); tempDen2->Scale(1/100.);// It was filled 100 times with the same value FullBuildFromFits->Add(tempDen2); FullBuildFromFits->Divide(tempDen2); FullBuildFromFits->SetLineColor(4); // PartialBuildFromFits = (TH1D*)PartialBuildFromFits_3D->ProjectionZ(proNamePartialFitBuild->Data(), c3FitType, c3FitType, Gbin, Gbin); PartialBuildFromFits->Add(tempDen2); PartialBuildFromFits->Divide(tempDen2); PartialBuildFromFits->SetLineColor(2); } } }// term 4-particle }// c4 }// c3 }// c2 }// c1 TH1D *fTPNRejects3pion = (TH1D*)MyList->FindObject("fTPNRejects3pion2"); TH1D *fTPNRejects4pion1 = (TH1D*)MyList->FindObject("fTPNRejects4pion1"); //TH2D *c4QSFitNum2D = (TH2D*)MyList->FindObject("fc4QSFitNum"); //TH2D *c4QSFitDen2D = (TH2D*)MyList->FindObject("fc4QSFitDen"); //c4QSFitNum2D->RebinY(3); //c4QSFitDen2D->RebinY(3); cout<<"Done getting Histograms"<<endl; TF1 *Unity=new TF1("Unity","1",0,100); Unity->SetLineStyle(2); int ch1_2=0,ch2_2=0; int ch1_3=0,ch2_3=0,ch3_3=0; int ch1_4=0,ch2_4=0,ch3_4=0,ch4_4=0; if(SameCharge && CHARGE==+1) {ch1_2=1; ch2_2=1; ch1_3=1; ch2_3=1; ch3_3=1; ch1_4=1; ch2_4=1; ch3_4=1; ch4_4=1;} if(!SameCharge){ ch1_2=0; ch2_2=1; // if(CHARGE==-1) { ch1_3=0; ch2_3=0; ch3_3=1; } else { ch1_3=0; ch2_3=1; ch3_3=1; } // if(CHARGE==-1 && MixedCharge4pionType==1) {ch1_4=0; ch2_4=0; ch3_4=0; ch4_4=1;} if(CHARGE==+1 && MixedCharge4pionType==1) {ch1_4=0; ch2_4=1; ch3_4=1; ch4_4=1;} if(MixedCharge4pionType==2) {ch1_4=0; ch2_4=0; ch3_4=1; ch4_4=1;} } /////////////////////////////////////////////////////////// // 2-pion section cout<<"2-pion section"<<endl; TCanvas *can1 = new TCanvas("can1", "can1",11,53,700,500); 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->SetGridx(); can1->SetGridy(); can1->SetFrameFillColor(0); can1->SetFrameBorderMode(0); can1->SetFrameBorderMode(0); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02); TLegend *legend1 = new TLegend(.6,.67,.87,.87,NULL,"brNDC"); legend1->SetBorderSize(1); legend1->SetTextSize(.04);// small .03; large .036 legend1->SetFillColor(0); TH1D *TERM1_2=(TH1D*)TwoParticle[ch1_2][ch2_2][0]->Clone(); TH1D *TERM2_2=(TH1D*)TwoParticle[ch1_2][ch2_2][1]->Clone(); // if(SameCharge){ TERM1_2->Multiply(MRC_SC_2[0]); TERM2_2->Multiply(MRC_SC_2[1]); }else { TERM1_2->Multiply(MRC_MC_2[0]); TERM2_2->Multiply(MRC_MC_2[1]); } TH1D *C2=(TH1D*)TERM1_2->Clone(); C2->Divide(TERM2_2); C2->GetXaxis()->SetRangeUser(0,0.1); C2->SetMinimum(0.98); C2->SetMaximum(2.5); /*C2->GetXaxis()->SetTitleSize(0.06); C2->GetXaxis()->SetLabelSize(0.06); C2->GetYaxis()->SetTitleSize(0.06); C2->GetYaxis()->SetLabelSize(0.06); C2->GetXaxis()->SetTitleOffset(1.05); C2->GetYaxis()->SetTitleOffset(1.1); C2->GetXaxis()->SetNdivisions(606); C2->GetYaxis()->SetNdivisions(505);*/ C2->Draw(); TH1D *C2QS=(TH1D*)TERM1_2->Clone(); C2QS->Add(TERM2_2, -(1-TwoFrac)); C2QS->Scale(1/TwoFrac); for(int bin=1; bin<=C2QS->GetNbinsX(); bin++) { Float_t K2 = 1.0; if(FSICorrection) K2 = FSICorrelation(ch1_2, ch2_2, C2QS->GetXaxis()->GetBinCenter(bin)); C2QS->SetBinContent(bin, C2QS->GetBinContent(bin)/K2); C2QS->SetBinError(bin, C2QS->GetBinError(bin)/K2); } C2QS->Divide(TERM2_2); if(SameCharge && MuonCorrection) C2QS->Multiply(C2muonCorrectionSC); if(!SameCharge && MuonCorrection) C2QS->Multiply(C2muonCorrectionMC); C2QS->SetMarkerColor(2); C2QS->SetLineColor(2); if(!MCcase) C2QS->Draw("same"); // To visualize the Qcut and Norm Regions //TH1D *QcutRegion = new TH1D("QcutRegion","",400,0,2); //TH1D *NormRegion1 = new TH1D("NormRegion1","",400,0,2); //TH1D *NormRegion2 = new TH1D("NormRegion2","",400,0,2); //TERM1_2->Scale(1/TERM1_2->Integral()); //for(int bin=1; bin<=16; bin++) QcutRegion->SetBinContent(bin,TERM1_2->GetBinContent(bin)); //for(int bin=213; bin<=220; bin++) NormRegion1->SetBinContent(bin,Two_ex[0][0][0][0]->GetBinContent(bin)); //for(int bin=31; bin<=35; bin++) NormRegion2->SetBinContent(bin,Two_ex[0][0][0][0]->GetBinContent(bin)); //QcutRegion->SetFillColor(4); //NormRegion1->SetFillColor(2); //NormRegion2->SetFillColor(3); //TERM1_2->GetYaxis()->SetTitle("Pair probability"); //TERM1_2->GetYaxis()->SetTitleOffset(1.3); //TERM1_2->GetXaxis()->SetRangeUser(0,1.6); //TERM1_2->Draw(); //QcutRegion->Draw("same"); //cout<<TERM1_2->Integral(1,16) / TERM1_2->Integral()<<endl; // Print 2-pion values /*for(int bin=1; bin<=20; bin++){ cout<<C2QS->GetBinContent(bin)<<", "; } cout<<endl; for(int bin=1; bin<=20; bin++){ cout<<C2QS->GetBinError(bin)<<", "; } cout<<endl;*/ double C2refPoints_0p65[20]={-0.218603, 1.00087, 1.00215, 1.00209, 1.00202, 1.00215, 1.00162, 1.00129, 1.00118, 1.00076, 1.00049, 1.00021, 0.999996, 0.999733, 0.999622, 0.999554, 0.999764, 0.999679, 0.999677, 0.999641}; double C2refPoints_e_0p65[20]={0.309152, 0.000473757, 0.000305599, 0.00022603, 0.000179596, 0.000149295, 0.000128009, 0.000112371, 0.00010047, 9.11363e-05, 8.37061e-05, 7.76871e-05, 7.27445e-05, 6.86217e-05, 6.5155e-05, 6.22148e-05, 5.9718e-05, 5.75613e-05, 5.57066e-05, 5.41028e-05}; double C2refPoints_0p75[20]={-0.135326, 0.968257, 0.984681, 0.991465, 0.99503, 0.997249, 0.998093, 0.998645, 0.999135, 0.999172, 0.999237, 0.999213, 0.999197, 0.999095, 0.999105, 0.999123, 0.999385, 0.999366, 0.999408, 0.999411}; double C2refPoints_e_0p75[20]={0.19138, 0.000422359, 0.000272193, 0.000201183, 0.000159788, 0.000132797, 0.000113845, 9.99271e-05, 8.93377e-05, 8.10334e-05, 7.44239e-05, 6.90702e-05, 6.46743e-05, 6.10077e-05, 5.79248e-05, 5.53103e-05, 5.30903e-05, 5.11725e-05, 4.95234e-05, 4.80974e-05}; TH1D *C2ref_0p65=(TH1D*)C2QS->Clone(); TH1D *C2ref_0p75=(TH1D*)C2QS->Clone(); for(int bin=1; bin<=20; bin++){ C2ref_0p65->SetBinContent(bin, C2refPoints_0p65[bin-1]); C2ref_0p65->SetBinError(bin, C2refPoints_e_0p65[bin-1]); C2ref_0p75->SetBinContent(bin, C2refPoints_0p75[bin-1]); C2ref_0p75->SetBinError(bin, C2refPoints_e_0p75[bin-1]); } C2ref_0p65->SetMarkerColor(4); C2ref_0p65->SetLineColor(4); C2ref_0p75->SetMarkerColor(6); C2ref_0p75->SetLineColor(6); //C2ref_0p65->Draw("same"); //C2ref_0p75->Draw("same"); TH1D *UnitMultC2=(TH1D*)UnitMult[0][0][0]->Clone(); TH1D *UnitMultC2den=(TH1D*)UnitMult[0][0][1]->Clone(); UnitMultC2->Divide(UnitMultC2den); UnitMultC2->GetXaxis()->SetRangeUser(0,0.3); UnitMultC2->SetMinimum(0.97); UnitMultC2->SetMaximum(1.25); UnitMultC2->Draw(); TH1D *UnitMultC2QS=(TH1D*)UnitMult[0][0][0]->Clone(); UnitMultC2QS->Add(UnitMult[0][0][1], -(1-TwoFrac)); UnitMultC2QS->Scale(1/TwoFrac); for(int bin=1; bin<=UnitMultC2QS->GetNbinsX(); bin++) { Float_t K2 = 1.0; if(FSICorrection) K2 = FSICorrelation(ch1_2, ch2_2, UnitMultC2QS->GetXaxis()->GetBinCenter(bin)); UnitMultC2QS->SetBinContent(bin, UnitMultC2QS->GetBinContent(bin)/K2); UnitMultC2QS->SetBinError(bin, UnitMultC2QS->GetBinError(bin)/K2); } UnitMultC2QS->Divide(UnitMultC2den); UnitMultC2QS->SetMarkerColor(2); UnitMultC2QS->SetLineColor(2); UnitMultC2QS->Draw("same"); /////////////////////////////////////////////////////////// // 3-pion cout<<"3-pion section"<<endl; TCanvas *can2 = new TCanvas("can2", "can2",600,53,700,500); can2->SetHighLightColor(2); can2->Range(-0.7838115,-1.033258,0.7838115,1.033258); 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); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02); TLegend *legend2 = (TLegend*)legend1->Clone(); legend2->SetX1(0.45); legend2->SetX2(0.95); legend2->SetY1(0.6); legend2->SetY2(0.95); TH1D *TERM1_3=(TH1D*)ThreeParticle[ch1_3][ch2_3][ch3_3][0]->Clone(); TH1D *TERM2_3=(TH1D*)ThreeParticle[ch1_3][ch2_3][ch3_3][1]->Clone(); TH1D *TERM3_3=(TH1D*)ThreeParticle[ch1_3][ch2_3][ch3_3][2]->Clone(); TH1D *TERM4_3=(TH1D*)ThreeParticle[ch1_3][ch2_3][ch3_3][3]->Clone(); TH1D *TERM5_3=(TH1D*)ThreeParticle[ch1_3][ch2_3][ch3_3][4]->Clone(); // if(SameCharge){ TERM1_3->Multiply(MRC_SC_3[0]); TERM2_3->Multiply(MRC_SC_3[1]); TERM3_3->Multiply(MRC_SC_3[1]); TERM4_3->Multiply(MRC_SC_3[1]); TERM5_3->Multiply(MRC_SC_3[2]); }else{ if(CHARGE==-1){// --+ but MRC is stored as -++ TERM1_3->Multiply(MRC_MC_3[0]); TERM2_3->Multiply(MRC_MC_3[2]); TERM3_3->Multiply(MRC_MC_3[1]); TERM4_3->Multiply(MRC_MC_3[1]); TERM5_3->Multiply(MRC_MC_3[3]); }else{ TERM1_3->Multiply(MRC_MC_3[0]); TERM2_3->Multiply(MRC_MC_3[1]); TERM3_3->Multiply(MRC_MC_3[1]); TERM4_3->Multiply(MRC_MC_3[2]); TERM5_3->Multiply(MRC_MC_3[3]); } } TH1D *N3QS = (TH1D*)TERM1_3->Clone(); N3QS->Add(TERM5_3, -f31); N3QS->Add(TERM2_3, -f32); N3QS->Add(TERM3_3, -f32); N3QS->Add(TERM4_3, -f32); N3QS->Scale(1/f33); N3QS->Multiply(K3avg[ch1_3][ch2_3][ch3_3][0]); TH1D *C3QS = (TH1D*)N3QS->Clone(); C3QS->Divide(TERM5_3); if(SameCharge) C3QS->Multiply(C3muonCorrectionSC[0]); else C3QS->Multiply(C3muonCorrectionMC[0]); C3QS->SetMarkerStyle(20); C3QS->SetMarkerColor(4); C3QS->GetXaxis()->SetRangeUser(0,0.1); C3QS->SetMinimum(0.9); C3QS->SetMaximum(5.0); TH1D *n3QS = (TH1D*)N3QS->Clone(); TH1D *c3QS = (TH1D*)N3QS->Clone(); for(int bin=1; bin<=n3QS->GetNbinsX(); bin++){ if(n3QS->GetBinContent(bin)==0) continue; double MuonCorr1=1.0, MuonCorr2=1.0, MuonCorr3=1.0, MuonCorr4=1.0; if(SameCharge){ MuonCorr1 = C3muonCorrectionSC[0]->GetBinContent(bin); MuonCorr2 = C3muonCorrectionSC[1]->GetBinContent(bin); MuonCorr3 = MuonCorr2; MuonCorr4 = MuonCorr2; }else if(CHARGE==-1){ MuonCorr1 = C3muonCorrectionMC[0]->GetBinContent(bin); MuonCorr2 = C3muonCorrectionMC[2]->GetBinContent(bin);// -- MuonCorr3 = C3muonCorrectionMC[1]->GetBinContent(bin);// -+ MuonCorr4 = MuonCorr3;// -+ }else{ MuonCorr1 = C3muonCorrectionMC[0]->GetBinContent(bin); MuonCorr2 = C3muonCorrectionMC[1]->GetBinContent(bin);// -+ MuonCorr3 = MuonCorr2; MuonCorr4 = C3muonCorrectionMC[2]->GetBinContent(bin); } double value = n3QS->GetBinContent(bin) * MuonCorr1; value -= ( TERM2_3->GetBinContent(bin) - (1-TwoFrac)*TERM5_3->GetBinContent(bin)) * K3avg[ch1_3][ch2_3][ch3_3][1]->GetBinContent(bin) / TwoFrac * MuonCorr2; value -= ( TERM3_3->GetBinContent(bin) - (1-TwoFrac)*TERM5_3->GetBinContent(bin)) * K3avg[ch1_3][ch2_3][ch3_3][2]->GetBinContent(bin) / TwoFrac * MuonCorr3; value -= ( TERM4_3->GetBinContent(bin) - (1-TwoFrac)*TERM5_3->GetBinContent(bin)) * K3avg[ch1_3][ch2_3][ch3_3][3]->GetBinContent(bin) / TwoFrac * MuonCorr4; //value -= TERM5_3->GetBinContent(bin)*(MuonCorr1 - MuonCorr2 - MuonCorr3 - MuonCorr4); value += 2*TERM5_3->GetBinContent(bin); // n3QS->SetBinContent(bin, value); c3QS->SetBinContent(bin, value + TERM5_3->GetBinContent(bin)); } c3QS->Divide(TERM5_3); C3QS->Draw(); c3QS->Draw("same"); //int lowBin = C3QS->GetXaxis()->FindBin(Cutoff_FullWeight_Q3[Mbin]); //int highBin = C3QS->GetXaxis()->FindBin(Cutoff_FullWeight_Q3[Mbin]); //double SF=C3QS->Integral(lowBin, highBin); SF /= TPFullWeight_ThreeParticle[ch1_3]->Integral(lowBin, highBin); //TPFullWeight_ThreeParticle[ch1_3]->Scale(SF); // if(SameCharge) TPFullWeight_ThreeParticle[ch1_3]->Draw("same"); //cout<<TPFullWeight_ThreeParticle[ch1_3]->GetBinContent(2)<<endl; //TH1D *C3raw = (TH1D*)TERM1_3->Clone(); //C3raw->Divide(TERM5_3); //C3raw->SetMarkerColor(4); //C3raw->Draw("same"); legend2->AddEntry(C3QS,"C_{3}^{QS}","p"); legend2->AddEntry(c3QS,"c_{3}^{QS}{2-pion removal}","p"); if(SameCharge) legend2->AddEntry(TPFullWeight_ThreeParticle[ch1_3],"C_{3}^{QS} built","l"); legend2->Draw("same"); /////////////////////////////////////// // C3QS built ratio TCanvas *can2_2 = new TCanvas("can2_2", "can2_2",1200,53,700,500); can2_2->SetHighLightColor(2); can2_2->Range(-0.7838115,-1.033258,0.7838115,1.033258); gStyle->SetOptFit(0111); can2_2->SetFillColor(10); can2_2->SetBorderMode(0); can2_2->SetBorderSize(2); can2_2->SetGridx(); can2_2->SetGridy(); can2_2->SetFrameFillColor(0); can2_2->SetFrameBorderMode(0); can2_2->SetFrameBorderMode(0); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02); TLegend *legend2_2 = (TLegend*)legend1->Clone(); legend2_2->SetX1(0.25); legend2_2->SetX2(0.6); legend2_2->SetY1(0.8); legend2_2->SetY2(0.95); TH1D *C3QSratio = (TH1D*)C3QS->Clone(); //TH1D *Chi2NDF_PointSize_3 = new TH1D("Chi2NDF_PointSize_3","",40,-0.5,39.5); //TH1D *Chi2NDF_FullSize_3 = new TH1D("Chi2NDF_FullSize_3","",40,-0.5,39.5); TH1D *Chi2NDF_PointSize_3 = new TH1D("Chi2NDF_PointSize_3","",100,-0.5,99.5); TH1D *Chi2NDF_FullSize_3 = new TH1D("Chi2NDF_FullSize_3","",100,-0.5,99.5); Chi2NDF_PointSize_3->SetLineColor(4); Chi2NDF_FullSize_3->SetLineColor(2); Chi2NDF_PointSize_3->SetMarkerColor(4); Chi2NDF_FullSize_3->SetMarkerColor(2); Chi2NDF_PointSize_3->GetXaxis()->SetTitle("Coherent fraction (%)"); Chi2NDF_PointSize_3->GetYaxis()->SetTitle("#sqrt{#chi^{2}}"); if(SameCharge && CollisionType==0){ TH1D *tempDen = (TH1D*)TPFullWeight_ThreeParticle_2D[ch1_3]->ProjectionY("TPFullWeight3_Den", 4, 4); TH1D *tempDenNeg = (TH1D*)TPNegFullWeight_ThreeParticle_2D[ch1_3]->ProjectionY("TPNegFullWeight3_Den", 4, 4); tempDen->Add(tempDenNeg);// Add Pos and Neg Den for(int binG=5; binG<=104; binG++){// was 44 TString *proName=new TString("TPFullWeight3_"); *proName += binG; TH1D *tempNum = (TH1D*)TPFullWeight_ThreeParticle_2D[ch1_3]->ProjectionY(proName->Data(), binG, binG); proName->Append("_Neg"); TH1D *tempNumNeg = (TH1D*)TPNegFullWeight_ThreeParticle_2D[ch1_3]->ProjectionY(proName->Data(), binG, binG); // // Add Pos and Neg Num tempNum->Add(tempNumNeg); // //tempNum->Add(tempDen);// now added in histogram read section tempNum->Divide(tempDen); //lowBin = C3QS->GetXaxis()->FindBin(Cutoff_FullWeight_Q3[Mbin]); //highBin = C3QS->GetXaxis()->FindBin(Cutoff_FullWeight_Q3[Mbin]); //SF=C3QS->Integral(lowBin, highBin); //SF /= tempNum->Integral(lowBin, highBin); //tempNum->Scale(SF); double Chi2=0; double NDF=0; for(int binQ3=3; binQ3<=3; binQ3++){ if(tempNum->GetBinContent(binQ3) <=0) continue; double value = C3QS->GetBinContent(binQ3) - tempNum->GetBinContent(binQ3); double err = C3QS->GetBinError(binQ3); if(err <=0) continue; Chi2 += pow(value / err,2); // NDF += 1; } //if(binG<25) {Chi2NDF_PointSize_3->SetBinContent(1 + 2*(binG-5), sqrt(Chi2)/NDF); Chi2NDF_PointSize_3->SetBinError(1 + 2*(binG-5), 0.001);} //else {Chi2NDF_FullSize_3->SetBinContent(1 + 2*(binG-25), sqrt(Chi2)/NDF); Chi2NDF_FullSize_3->SetBinError(1 + 2*(binG-25), 0.001);} if(binG<55) {Chi2NDF_PointSize_3->SetBinContent(1 + 2*(binG-5), sqrt(Chi2)/NDF); Chi2NDF_PointSize_3->SetBinError(1 + 2*(binG-5), 0.001);} else {Chi2NDF_FullSize_3->SetBinContent(1 + 2*(binG-55), sqrt(Chi2)/NDF); Chi2NDF_FullSize_3->SetBinError(1 + 2*(binG-55), 0.001);} //if(NDF>0) cout<<binG<<" "<<sqrt(Chi2)/NDF<<endl; } Chi2NDF_PointSize_3->SetMarkerStyle(20); Chi2NDF_FullSize_3->SetMarkerStyle(20); C3QSratio->Divide(TPFullWeight_ThreeParticle[ch1_3]); // C3QSratio->SetMinimum(0.94); C3QSratio->SetMaximum(1.02); C3QSratio->GetYaxis()->SetTitleOffset(1.2); C3QSratio->GetYaxis()->SetTitle("C_{3}^{QS} / C_{3}(built)"); C3QSratio->Draw(); Unity->Draw("same"); Chi2NDF_PointSize_3->Draw(); Chi2NDF_FullSize_3->Draw("same"); legend2_2->AddEntry(Chi2NDF_PointSize_3,"R_{coh}=0 (Point Source)","p"); legend2_2->AddEntry(Chi2NDF_FullSize_3,"R_{coh}=R_{ch}","p"); legend2_2->Draw("same"); } // r3 TH1D *r3; if(SameCharge && CollisionType==0){ r3 = (TH1D*)n3QS->Clone(); TPN_ThreeParticle[ch1_3]->Multiply(MRC_SC_3[2]); r3->Divide(TPN_ThreeParticle[ch1_3]); r3->GetXaxis()->SetRangeUser(0,0.08); r3->SetMinimum(0.5); r3->SetMaximum(2.5); r3->GetYaxis()->SetTitle("r_{3}"); // r3->Draw(); //TPN_ThreeParticle[ch1_3]->Draw(); //fTPNRejects3pion->SetLineColor(2); //fTPNRejects3pion->Draw("same"); } // Print 3-pion points for(int bin=1; bin<=10; bin++){ //cout<<C3QS->GetBinContent(bin)<<", "; //cout<<c3QS->GetBinContent(bin)<<", "; //cout<<TPFullWeight_ThreeParticle[ch1_3]->GetBinContent(bin)<<", "; } //cout<<endl; for(int bin=1; bin<=10; bin++){ //cout<<C3QS->GetBinError(bin)<<", "; //cout<<c3QS->GetBinError(bin)<<", "; } //cout<<endl; //////////////////////////////////////////////////////////////// // 4-pion cout<<"4-pion section"<<endl; TCanvas *can3 = new TCanvas("can3", "can3",11,600,700,500);// 60 was 600 can3->SetHighLightColor(2); can3->Range(-0.7838115,-1.033258,0.7838115,1.033258); gStyle->SetOptFit(0111); can3->SetFillColor(10); can3->SetBorderMode(0); can3->SetBorderSize(2); can3->SetGridx(); can3->SetGridy(); can3->SetFrameFillColor(0); can3->SetFrameBorderMode(0); can3->SetFrameBorderMode(0); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02); TLegend *legend3 = (TLegend*)legend1->Clone(); legend3->SetX1(0.47); legend3->SetX2(0.98); legend3->SetY1(0.7); legend3->SetY2(0.98); TH1D *TERMS_4[13]; for(int index=0; index<13; index++){ TERMS_4[index] = (TH1D*)FourParticle[ch1_4][ch2_4][ch3_4][ch4_4][index]->Clone(); // TH1D *MRC_temp; if(SameCharge){ if(index==0) MRC_temp = (TH1D*)MRC_SC_4[0]->Clone(); else if(index<=4) MRC_temp = (TH1D*)MRC_SC_4[1]->Clone(); else if(index<=10) MRC_temp = (TH1D*)MRC_SC_4[2]->Clone(); else if(index==11) MRC_temp = (TH1D*)MRC_SC_4[3]->Clone(); else MRC_temp = (TH1D*)MRC_SC_4[4]->Clone(); }else if(CHARGE==-1 && MixedCharge4pionType==1){ if(index==0) MRC_temp = (TH1D*)MRC_MC1_4[0]->Clone();//0 else if(index!=1 && index<=4) MRC_temp = (TH1D*)MRC_MC1_4[1]->Clone();//1 else if(index==1) MRC_temp = (TH1D*)MRC_MC1_4[2]->Clone();//2 else if(index==7 || index==9 || index==10) MRC_temp = (TH1D*)MRC_MC1_4[3]->Clone();//3 else if(index!=12) MRC_temp = (TH1D*)MRC_MC1_4[4]->Clone();//4 else MRC_temp = (TH1D*)MRC_MC1_4[5]->Clone();//5 }else if(CHARGE==+1 && MixedCharge4pionType==1){ if(index==0) MRC_temp = (TH1D*)MRC_MC1_4[0]->Clone(); else if(index<=3) MRC_temp = (TH1D*)MRC_MC1_4[1]->Clone(); else if(index==4) MRC_temp = (TH1D*)MRC_MC1_4[2]->Clone(); else if(index==5 || index==6 || index==7) MRC_temp = (TH1D*)MRC_MC1_4[3]->Clone(); else if(index!=12) MRC_temp = (TH1D*)MRC_MC1_4[4]->Clone(); else MRC_temp = (TH1D*)MRC_MC1_4[5]->Clone(); }else{// --++ if(index==0) MRC_temp = (TH1D*)MRC_MC2_4[0]->Clone(); else if(index<=4) MRC_temp = (TH1D*)MRC_MC2_4[1]->Clone(); else if(index==5 || index==10) MRC_temp = (TH1D*)MRC_MC2_4[2]->Clone();//2 else if(index<=9) MRC_temp = (TH1D*)MRC_MC2_4[3]->Clone();// 3 else if(index==11) MRC_temp = (TH1D*)MRC_MC2_4[4]->Clone();// 4 else MRC_temp = (TH1D*)MRC_MC2_4[5]->Clone(); } // TERMS_4[index]->Multiply(MRC_temp); } TH1D *C4raw=(TH1D*)TERMS_4[2]->Clone(); C4raw->Divide(TERMS_4[12]); C4raw->GetXaxis()->SetRangeUser(0,0.2); TH1D *N4QS=(TH1D*)TERMS_4[0]->Clone(); // N4QS->Add(TERMS_4[1], -f43); N4QS->Add(TERMS_4[2], -f43); N4QS->Add(TERMS_4[3], -f43); N4QS->Add(TERMS_4[4], -f43); // N4QS->Add(TERMS_4[5], -f42); N4QS->Add(TERMS_4[6], -f42); N4QS->Add(TERMS_4[7], -f42); N4QS->Add(TERMS_4[8], -f42); N4QS->Add(TERMS_4[9], -f42); N4QS->Add(TERMS_4[10], -f42); // N4QS->Add(TERMS_4[12], -f41); // N4QS->Scale(1/f44); //K4avg[ch1_4][ch2_4][ch3_4][ch4_4][0]->Scale(1.005);// K scale variation N4QS->Multiply(K4avg[ch1_4][ch2_4][ch3_4][ch4_4][0]); //N4QS->Divide(K4avg[ch1_4][ch2_4][ch3_4][ch4_4][1]);// K factorization test (MC1) //N4QS->Divide(K4avg[ch1_4][ch2_4][ch3_4][ch4_4][11]);// K factorization test (MC2) TH1D *C4QS=(TH1D*)N4QS->Clone(); C4QS->GetXaxis()->SetRangeUser(0,0.179); C4QS->SetMinimum(0.5); C4QS->SetMarkerColor(4); // /*TH1D *C4QS_basic=(TH1D*)TERMS_4[0]->Clone(); for(int ii=1; ii<=C4QS_basic->GetNbinsX(); ii++){ double value = C4QS_basic->GetBinContent(ii) - TERMS_4[12]->GetBinContent(ii); value /= f44; value += TERMS_4[12]->GetBinContent(ii); C4QS_basic->SetBinContent(ii, value); } C4QS_basic->Divide(TERMS_4[12]); //C4QS_basic->Multiply(K4avg[ch1_4][ch2_4][ch3_4][ch4_4][0]); C4QS_basic->SetMarkerColor(1); //C4QS_basic->Draw("same"); */ // TH1D *C4_pieces[12]; for(int i=0; i<12; i++){ C4_pieces[i]=(TH1D*)TERMS_4[i]->Clone(); C4_pieces[i]->Divide(TERMS_4[12]); if(i==0) C4_pieces[i]->SetMarkerColor(1); else if(i<5) C4_pieces[i]->SetMarkerColor(2); else if(i<10) C4_pieces[i]->SetMarkerColor(3); else C4_pieces[i]->SetMarkerColor(4); } C4_pieces[0]->GetXaxis()->SetRangeUser(0,0.2); C4_pieces[0]->SetMinimum(0.5); C4_pieces[0]->SetMaximum(3); //C4_pieces[0]->Draw(); //C4_pieces[1]->Draw("same"); //C4_pieces[5]->Draw("same"); //C4_pieces[11]->Draw("same"); // TH1D *C4_2pairs_FromSquares = (TH1D*)TERMS_4[5]->Clone(); C4_2pairs_FromSquares->Divide(TERMS_4[12]); for(int bin=1; bin<=C4_2pairs_FromSquares->GetNbinsX(); bin++){ double value = C4_2pairs_FromSquares->GetBinContent(bin); value -= 1; value *= value; value += 1; C4_2pairs_FromSquares->SetBinContent(bin, value); } C4_2pairs_FromSquares->SetMarkerColor(6); //C4_2pairs_FromSquares->Draw("same"); // TH1D *C4QS_2pairs=(TH1D*)TERMS_4[11]->Clone(); C4QS_2pairs->Add(TERMS_4[12], -pow(1-TwoFrac,2));// N1^4 C4QS_2pairs->Add(TERMS_4[5], -2*TwoFrac*(1-TwoFrac));// N2 * N1^2 C4QS_2pairs->Scale(1/pow(TwoFrac,2)); C4QS_2pairs->Multiply(K4avg[ch1_4][ch2_4][ch3_4][ch4_4][11]); C4QS_2pairs->Divide(TERMS_4[12]); C4QS_2pairs->SetMarkerColor(6); //C4QS_2pairs->Draw("same"); // TH1D *n4QS = (TH1D*)N4QS->Clone(); TH1D *c4QS = (TH1D*)N4QS->Clone(); TH1D *c4QS_RemovalStage1 = (TH1D*)N4QS->Clone(); TH1D *c4QS_RemovalStage2 = (TH1D*)N4QS->Clone(); TH1D *c4QS_RemovalStage3 = (TH1D*)N4QS->Clone(); for(int bin=1; bin<=n4QS->GetNbinsX(); bin++){ if(n4QS->GetBinContent(bin)==0) continue; bool exitCode=kFALSE; for(int ii=0; ii<13; ii++) {if(TERMS_4[ii]->GetBinContent(bin) < 1) exitCode=kTRUE;} if(exitCode) continue; //cout<<TERMS_4[0]->GetBinContent(bin)<<endl; // double N4QSvalue = N4QS->GetBinContent(bin); double SubtractionTerm[11] = {0}; double ErrorTerm[12] = {0}; ErrorTerm[0] = n4QS->GetBinError(bin); // // 3-pion terms if(SameCharge || MixedCharge4pionType==1){ SubtractionTerm[0] = (TERMS_4[1]->GetBinContent(bin) - f32*TERMS_4[5]->GetBinContent(bin) - f32*TERMS_4[6]->GetBinContent(bin) - f32*TERMS_4[8]->GetBinContent(bin) - f31*TERMS_4[12]->GetBinContent(bin)) / f33;// 5 6 8 SubtractionTerm[0] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][1]->GetBinContent(bin); ErrorTerm[1] = sqrt(TERMS_4[1]->GetBinContent(bin) + f32*f32*TERMS_4[5]->GetBinContent(bin) + f32*f32*TERMS_4[6]->GetBinContent(bin) + f32*f32*TERMS_4[8]->GetBinContent(bin) + f31*f31*TERMS_4[12]->GetBinContent(bin)) / f33; ErrorTerm[1] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][1]->GetBinContent(bin); // SubtractionTerm[1] = (TERMS_4[2]->GetBinContent(bin) - f32*TERMS_4[5]->GetBinContent(bin) - f32*TERMS_4[7]->GetBinContent(bin) - f32*TERMS_4[9]->GetBinContent(bin) - f31*TERMS_4[12]->GetBinContent(bin)) / f33; SubtractionTerm[1] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][2]->GetBinContent(bin); ErrorTerm[2] = sqrt(TERMS_4[2]->GetBinContent(bin) + f32*f32*TERMS_4[5]->GetBinContent(bin) + f32*f32*TERMS_4[7]->GetBinContent(bin) + f32*f32*TERMS_4[9]->GetBinContent(bin) + f31*f31*TERMS_4[12]->GetBinContent(bin)) / f33; ErrorTerm[2] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][2]->GetBinContent(bin); // SubtractionTerm[2] = (TERMS_4[3]->GetBinContent(bin) - f32*TERMS_4[6]->GetBinContent(bin) - f32*TERMS_4[7]->GetBinContent(bin) - f32*TERMS_4[10]->GetBinContent(bin) - f31*TERMS_4[12]->GetBinContent(bin)) / f33; SubtractionTerm[2] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][3]->GetBinContent(bin); ErrorTerm[3] = sqrt(TERMS_4[3]->GetBinContent(bin) + f32*f32*TERMS_4[6]->GetBinContent(bin) + f32*f32*TERMS_4[7]->GetBinContent(bin) + f32*f32*TERMS_4[10]->GetBinContent(bin) + f31*f31*TERMS_4[12]->GetBinContent(bin)) / f33; ErrorTerm[3] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][3]->GetBinContent(bin); // SubtractionTerm[3] = (TERMS_4[4]->GetBinContent(bin) - f32*TERMS_4[8]->GetBinContent(bin) - f32*TERMS_4[9]->GetBinContent(bin) - f32*TERMS_4[10]->GetBinContent(bin) - f31*TERMS_4[12]->GetBinContent(bin)) / f33; SubtractionTerm[3] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][4]->GetBinContent(bin); ErrorTerm[4] = sqrt(TERMS_4[4]->GetBinContent(bin) + f32*f32*TERMS_4[8]->GetBinContent(bin) + f32*f32*TERMS_4[9]->GetBinContent(bin) + f32*f32*TERMS_4[10]->GetBinContent(bin) + f31*f31*TERMS_4[12]->GetBinContent(bin)) / f33; ErrorTerm[4] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][4]->GetBinContent(bin); // } // 2-pion terms SubtractionTerm[4] = ( TERMS_4[5]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) / TwoFrac * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][5]->GetBinContent(bin); ErrorTerm[5] = sqrt( TERMS_4[5]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][5]->GetBinContent(bin) / TwoFrac; // SubtractionTerm[5] = ( TERMS_4[6]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][6]->GetBinContent(bin) / TwoFrac; ErrorTerm[6] = sqrt( TERMS_4[6]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][6]->GetBinContent(bin) / TwoFrac; // SubtractionTerm[6] = ( TERMS_4[7]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][7]->GetBinContent(bin) / TwoFrac; ErrorTerm[7] = sqrt( TERMS_4[7]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][7]->GetBinContent(bin) / TwoFrac; // SubtractionTerm[7] = ( TERMS_4[8]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][8]->GetBinContent(bin) / TwoFrac; ErrorTerm[8] = sqrt( TERMS_4[8]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][8]->GetBinContent(bin) / TwoFrac; // SubtractionTerm[8] = ( TERMS_4[9]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][9]->GetBinContent(bin) / TwoFrac; ErrorTerm[9] = sqrt( TERMS_4[9]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][9]->GetBinContent(bin) / TwoFrac; // SubtractionTerm[9] = ( TERMS_4[10]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][10]->GetBinContent(bin) / TwoFrac; ErrorTerm[10] = sqrt( TERMS_4[10]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][10]->GetBinContent(bin) / TwoFrac; // // // 2 2-pion terms: cos(q12*x12)*cos(q34*x34) SubtractionTerm[10] = TERMS_4[11]->GetBinContent(bin);// N22 SubtractionTerm[10] -= 2*(1-TwoFrac) * TERMS_4[5]->GetBinContent(bin); SubtractionTerm[10] += pow(1-TwoFrac,2) * TERMS_4[12]->GetBinContent(bin); SubtractionTerm[10] /= pow(TwoFrac,2); SubtractionTerm[10] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][11]->GetBinContent(bin); SubtractionTerm[10] *= C4muonCorrectionSC[3]->GetBinContent(bin); double intermed = ( TERMS_4[5]->GetBinContent(bin) - (1-TwoFrac)*TERMS_4[12]->GetBinContent(bin)) / TwoFrac; intermed *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][5]->GetBinContent(bin); intermed = intermed * C4muonCorrectionSC[2]->GetBinContent(bin); SubtractionTerm[10] -= 2*intermed; SubtractionTerm[10] += 2*TERMS_4[12]->GetBinContent(bin); // ErrorTerm[11] = TERMS_4[11]->GetBinContent(bin); ErrorTerm[11] += pow(1-TwoFrac,4) * TERMS_4[12]->GetBinContent(bin); ErrorTerm[11] += pow(2*(1-TwoFrac)*TwoFrac,2) * TERMS_4[5]->GetBinContent(bin); ErrorTerm[11] = sqrt(ErrorTerm[11]); ErrorTerm[11] /= pow(TwoFrac,2); ErrorTerm[11] *= K4avg[ch1_4][ch2_4][ch3_4][ch4_4][11]->GetBinContent(bin); ErrorTerm[11] = pow(ErrorTerm[11],2); ErrorTerm[11] += pow( 2*sqrt(TERMS_4[5]->GetBinContent(bin) + pow(1-TwoFrac,2)*TERMS_4[12]->GetBinContent(bin)) * K4avg[ch1_4][ch2_4][ch3_4][ch4_4][5]->GetBinContent(bin) / TwoFrac, 2); ErrorTerm[11] += 4*TERMS_4[12]->GetBinContent(bin); ErrorTerm[11] = sqrt(ErrorTerm[11]); // double FinalValue[4]={0}; double FinalValue_e[4]={0}; if(SameCharge) { //N4QSvalue = (N4QSvalue - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionSC[0]->GetBinContent(bin) + TERMS_4[12]->GetBinContent(bin); N4QSvalue *= C4muonCorrectionSC[0]->GetBinContent(bin); ErrorTerm[0] *= C4muonCorrectionSC[0]->GetBinContent(bin); }else if(MixedCharge4pionType==1) { //N4QSvalue = (N4QSvalue - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC1[0]->GetBinContent(bin) + TERMS_4[12]->GetBinContent(bin); N4QSvalue *= C4muonCorrectionMC1[0]->GetBinContent(bin); ErrorTerm[0] *= C4muonCorrectionMC1[0]->GetBinContent(bin); }else { //N4QSvalue = (N4QSvalue - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC2[0]->GetBinContent(bin) + TERMS_4[12]->GetBinContent(bin); N4QSvalue *= C4muonCorrectionMC2[0]->GetBinContent(bin); ErrorTerm[0] *= C4muonCorrectionMC2[0]->GetBinContent(bin); } for(int ii=0; ii<4; ii++) { FinalValue[ii] = N4QSvalue; FinalValue_e[ii] = pow(ErrorTerm[0],2); } if(SameCharge) { //FinalValue[0] -= 4*(SubtractionTerm[0] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionSC[1]->GetBinContent(bin); //FinalValue[0] += 6*(SubtractionTerm[4] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionSC[2]->GetBinContent(bin); FinalValue[0] -= 4*(SubtractionTerm[0] * C4muonCorrectionSC[1]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); FinalValue[0] += 6*(SubtractionTerm[4] * C4muonCorrectionSC[2]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); FinalValue[0] -= 3*(SubtractionTerm[10] - TERMS_4[12]->GetBinContent(bin)); FinalValue_e[0] += 4*pow(ErrorTerm[1]*C4muonCorrectionSC[1]->GetBinContent(bin),2); FinalValue_e[0] += 6*pow(ErrorTerm[5]*C4muonCorrectionSC[2]->GetBinContent(bin),2); FinalValue_e[0] += 3*pow(ErrorTerm[11]*C4muonCorrectionSC[3]->GetBinContent(bin),2); FinalValue_e[0] += (4*C4muonCorrectionSC[1]->GetBinContent(bin) + 6*C4muonCorrectionSC[2]->GetBinContent(bin) + 3*C4muonCorrectionSC[3]->GetBinContent(bin)) * TERMS_4[12]->GetBinContent(bin); // FinalValue[1] -= 6*SubtractionTerm[4] - 6*TERMS_4[12]->GetBinContent(bin);// 2-pion removal FinalValue_e[1] += 6*pow(ErrorTerm[5],2) + 6*TERMS_4[12]->GetBinContent(bin); FinalValue[2] = FinalValue[1] - 3*SubtractionTerm[10] + 3*TERMS_4[12]->GetBinContent(bin);// further 2-pair removal FinalValue_e[2] += FinalValue_e[1] + 3*pow(ErrorTerm[11],2) + 3*TERMS_4[12]->GetBinContent(bin); FinalValue[3] = SubtractionTerm[10]; }else{ if(MixedCharge4pionType==1 && CHARGE==-1){ //FinalValue[0] -= (SubtractionTerm[0] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC1[2]->GetBinContent(bin);// [2] is +++ MuonCorr FinalValue[0] -= SubtractionTerm[0] * C4muonCorrectionMC1[2]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin);// [2] is +++ MuonCorr FinalValue[0] -= 3*(SubtractionTerm[6] * C4muonCorrectionMC1[3]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); // FinalValue_e[0] += pow(ErrorTerm[1]*C4muonCorrectionMC1[2]->GetBinContent(bin),2); FinalValue[1] -= 3*SubtractionTerm[4] - 3*TERMS_4[12]->GetBinContent(bin); FinalValue_e[1] += 3*pow(ErrorTerm[5],2) + 3*TERMS_4[12]->GetBinContent(bin); }else if(MixedCharge4pionType==1 && CHARGE==+1){ //FinalValue[0] -= (SubtractionTerm[3] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC1[2]->GetBinContent(bin); FinalValue[0] -= SubtractionTerm[3] * C4muonCorrectionMC1[2]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin); FinalValue[0] -= 3*(SubtractionTerm[6] * C4muonCorrectionMC1[3]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); // FinalValue_e[0] += pow(ErrorTerm[4]*C4muonCorrectionMC1[2]->GetBinContent(bin),2); FinalValue[1] -= 3*SubtractionTerm[9] - 3*TERMS_4[12]->GetBinContent(bin); FinalValue_e[1] += 3*pow(ErrorTerm[10],2) + 3*TERMS_4[12]->GetBinContent(bin); }else{// --++ case //FinalValue[0] -= 2*(SubtractionTerm[4] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC2[2]->GetBinContent(bin);// old way FinalValue[0] -= 2*(SubtractionTerm[4] * C4muonCorrectionMC2[2]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); FinalValue[0] -= (SubtractionTerm[10] - TERMS_4[12]->GetBinContent(bin)); FinalValue[0] -= 4*(SubtractionTerm[6] * C4muonCorrectionMC2[3]->GetBinContent(bin) - TERMS_4[12]->GetBinContent(bin)); FinalValue_e[0] += 2*pow(ErrorTerm[5],2) + pow(ErrorTerm[11],2) + 2*TERMS_4[12]->GetBinContent(bin); // FinalValue[1] -= 2*(SubtractionTerm[4] - TERMS_4[12]->GetBinContent(bin)) * C4muonCorrectionMC2[2]->GetBinContent(bin); FinalValue_e[1] += 2*pow(ErrorTerm[5],2) + 2*TERMS_4[12]->GetBinContent(bin); } } for(int ii=0; ii<4; ii++) FinalValue_e[ii] = sqrt(FinalValue_e[ii]); // n4QS->SetBinContent(bin, FinalValue[0] - TERMS_4[12]->GetBinContent(bin)); n4QS->SetBinError(bin, sqrt( pow(FinalValue_e[0],2) + TERMS_4[12]->GetBinContent(bin))); c4QS->SetBinContent(bin, FinalValue[0]); c4QS->SetBinError(bin, FinalValue_e[0]); C4QS->SetBinContent(bin, N4QSvalue); C4QS->SetBinError(bin, ErrorTerm[0]); // c4QS_RemovalStage1->SetBinContent(bin, FinalValue[1]); c4QS_RemovalStage1->SetBinError(bin, FinalValue_e[1]); // c4QS_RemovalStage2->SetBinContent(bin, FinalValue[2]); c4QS_RemovalStage2->SetBinError(bin, FinalValue_e[2]); // c4QS_RemovalStage3->SetBinContent(bin, FinalValue[3]); c4QS_RemovalStage3->SetBinError(bin, FinalValue_e[3]); } C4QS->Divide(TERMS_4[12]); c4QS->Divide(TERMS_4[12]); c4QS->SetMarkerColor(1); c4QS->SetLineColor(1); // c4QS_RemovalStage1->Divide(TERMS_4[12]); c4QS_RemovalStage1->SetMarkerColor(2); c4QS_RemovalStage1->SetLineColor(2); c4QS_RemovalStage2->Divide(TERMS_4[12]); c4QS_RemovalStage2->SetMarkerColor(6); c4QS_RemovalStage2->SetLineColor(6); c4QS_RemovalStage3->Divide(TERMS_4[12]); c4QS_RemovalStage3->SetMarkerColor(7); c4QS_RemovalStage3->SetLineColor(7); // // if(SameCharge) { if(CollisionType==0) C4QS->SetMaximum(7); else C4QS->SetMaximum(12); }else if(MixedCharge4pionType==1) C4QS->SetMaximum(4); else C4QS->SetMaximum(3); if(CollisionType!=0) C4QS->GetXaxis()->SetRangeUser(0,0.5); C4QS->Draw(); c4QS_RemovalStage1->Draw("same"); if(SameCharge) c4QS_RemovalStage2->Draw("same"); //c4QS_RemovalStage3->Draw("same"); c4QS->Draw("same"); //cout<<TPFullWeight_FourParticle[ch1_4]->GetBinContent(9)<<endl; if(SameCharge && !MCcase) { //TPFullWeight_FourParticle[ch1_4]->Draw("same"); FullBuildFromFits->Draw("same"); PartialBuildFromFits->Draw("same"); } legend3->AddEntry(C4QS,"C_{4}^{QS}","p"); legend3->AddEntry(c4QS_RemovalStage1,"c_{4}^{QS}{2-pion removal}","p"); if(SameCharge) legend3->AddEntry(c4QS_RemovalStage2,"c_{4}^{QS}{2-pion+2-pair removal}","p"); if(SameCharge) legend3->AddEntry(c4QS,"c_{4}^{QS}{2-pion+2-pair+3-pion removal}","p"); if(!SameCharge && MixedCharge4pionType==1) legend3->AddEntry(c4QS,"c_{4}^{QS}{2-pion+3-pion removal}","p"); if(!SameCharge && MixedCharge4pionType==2) legend3->AddEntry(c4QS,"c_{4}^{QS}{2-pion+2-pair removal}","p"); //legend3->AddEntry(c4QS,"c_{4}^{QS}{2-pion+2-pair removal}","p"); if(SameCharge) legend3->AddEntry(TPFullWeight_FourParticle[ch1_4],"C_{4}^{QS} built","l"); legend3->Draw("same"); //K4avg[ch1_4][ch2_4][ch3_4][ch4_4][0]->Draw(); //TH1D *c4QSFitNum = (TH1D*)c4QSFitNum2D->ProjectionY("c4QSFitNum",5,5); //TH1D *c4QSFitDen = (TH1D*)c4QSFitDen2D->ProjectionY("c4QSFitDen",5,5); //c4QSFitNum->Divide(c4QSFitDen); //for(int bin=1; bin<20; bin++){ //double fc4 = //c4QSFitNum->SetBinContent(bin, (c4QSFitNum->GetBinContent(bin) - 1.0)* //} //c4QSFitNum->Draw("same"); //C4QS_Syst->Draw("E2 same"); //C4QS->Draw("same"); //C4raw->Draw(); //TH1D *testTerm = (TH1D*)TERMS_4[11]->Clone(); //testTerm->Divide(TERMS_4[12]); //testTerm->SetMinimum(1); testTerm->SetMaximum(1.4); testTerm->GetXaxis()->SetRangeUser(0,0.14); //testTerm->Draw(); //////////////////////////////////////////////////////////////// if(SameCharge && CollisionType==0){ TCanvas *can4 = new TCanvas("can4", "can4",600,600,700,500); can4->SetHighLightColor(2); can4->Range(-0.7838115,-1.033258,0.7838115,1.033258); gStyle->SetOptFit(0111); can4->SetFillColor(10); can4->SetBorderMode(0); can4->SetBorderSize(2); can4->SetGridx(); can4->SetGridy(); can4->SetFrameFillColor(0); can4->SetFrameBorderMode(0); can4->SetFrameBorderMode(0); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02); TLegend *legend3_2 = (TLegend*)legend1->Clone(); legend3_2->SetX1(0.45); legend3_2->SetX2(0.98); legend3_2->SetY1(0.6); legend3_2->SetY2(0.95); // //TH1D *Chi2NDF_PointSize_4 = new TH1D("Chi2NDF_PointSize_4","",40,-0.5,39.5); //TH1D *Chi2NDF_FullSize_4 = new TH1D("Chi2NDF_FullSize_4","",40,-0.5,39.5); TH1D *Chi2NDF_PointSize_4 = new TH1D("Chi2NDF_PointSize_4","",100,-0.5,99.5); TH1D *Chi2NDF_FullSize_4 = new TH1D("Chi2NDF_FullSize_4","",100,-0.5,99.5); Chi2NDF_PointSize_4->SetLineColor(4); Chi2NDF_FullSize_4->SetLineColor(2); Chi2NDF_PointSize_4->SetMarkerColor(4); Chi2NDF_FullSize_4->SetMarkerColor(2); Chi2NDF_PointSize_4->GetXaxis()->SetTitle("Coherent fraction (%)"); Chi2NDF_PointSize_4->GetYaxis()->SetTitle("#sqrt{#chi^{2}}"); TH1D *tempDen = (TH1D*)TPFullWeight_FourParticle_2D[ch1_4]->ProjectionY("TPFullWeight4_Den", 4, 4); TH1D *tempDenNeg = (TH1D*)TPNegFullWeight_FourParticle_2D[ch1_4]->ProjectionY("TPNegFullWeight4_Den", 4, 4); tempDen->Add(tempDenNeg);// Add Pos and Neg Weight for(int binG=5; binG<=104; binG++){// 44 TString *proName=new TString("TPFullWeight4_"); *proName += binG; TH1D *tempNum = (TH1D*)TPFullWeight_FourParticle_2D[ch1_4]->ProjectionY(proName->Data(), binG, binG); proName->Append("_Neg"); TH1D *tempNumNeg = (TH1D*)TPNegFullWeight_FourParticle_2D[ch1_4]->ProjectionY(proName->Data(), binG, binG); // // Add Pos and Neg Weights tempNum->Add(tempNumNeg); // //tempNum->Add(tempDen);// now added in InterpCorr section tempNum->Divide(tempDen); double Chi2=0; double NDF=0; for(int binQ4=4; binQ4<=4; binQ4++){ if(tempNum->GetBinContent(binQ4) <=0) continue; double value = C4QS->GetBinContent(binQ4) - tempNum->GetBinContent(binQ4); double err = pow(C4QS->GetBinError(binQ4),2); err = sqrt(err); if(err<=0) continue; Chi2 += pow(value / err,2); // NDF += 1; } //if(binG<25) {Chi2NDF_PointSize_4->SetBinContent(1 + 2*(binG-5), sqrt(fabs(Chi2))/NDF); Chi2NDF_PointSize_4->SetBinError(1 + 2*(binG-5), 0.001);} //else {Chi2NDF_FullSize_4->SetBinContent(1 + 2*(binG-25), sqrt(fabs(Chi2))/NDF); Chi2NDF_FullSize_4->SetBinError(1 + 2*(binG-25), 0.001);} if(binG<55) {Chi2NDF_PointSize_4->SetBinContent(1 + 2*(binG-5), sqrt(fabs(Chi2))/NDF); Chi2NDF_PointSize_4->SetBinError(1 + 2*(binG-5), 0.001);} else {Chi2NDF_FullSize_4->SetBinContent(1 + 2*(binG-55), sqrt(fabs(Chi2))/NDF); Chi2NDF_FullSize_4->SetBinError(1 + 2*(binG-55), 0.001);} } Chi2NDF_PointSize_4->SetMarkerStyle(20); Chi2NDF_FullSize_4->SetMarkerStyle(20); Chi2NDF_PointSize_4->Draw(); Chi2NDF_FullSize_4->Draw("same"); legend3_2->AddEntry(Chi2NDF_PointSize_4,"R_{coh}=0 (Point Source)","p"); legend3_2->AddEntry(Chi2NDF_FullSize_4,"R_{coh}=R_{ch}","p"); legend3_2->Draw("same"); } // Print 4-pion points for(int bin=1; bin<=12; bin++){ //cout<<C4QS->GetBinContent(bin)<<", "; //cout<<c4QS->GetBinContent(bin)<<", "; //cout<<TPFullWeight_FourParticle[ch1_4]->GetBinContent(bin)<<", "; //cout<<C4raw->GetBinContent(bin)<<", "; //cout<<K4avg[ch1_4][ch2_4][ch3_4][ch4_4][0]->GetBinContent(bin)<<", "; } cout<<endl; for(int bin=1; bin<=12; bin++){ //cout<<c4QS->GetBinContent(bin)<<", "; //cout<<C4QS->GetBinError(bin)<<", "; //cout<<c4QS->GetBinError(bin)<<", "; //cout<<C4raw->GetBinError(bin)<<", "; } //cout<<endl; //////////////////////////////////////////////////////////////// // r4 TF1 *ChaoticLimit_r42 = new TF1("ChaoticLimit_r42","6",0,10); ChaoticLimit_r42->SetLineStyle(2); TH1D *r42; if(SameCharge && CollisionType==0){ /*TCanvas *can5 = new TCanvas("can5", "can5",1200,600,700,500); can5->SetHighLightColor(2); can5->Range(-0.7838115,-1.033258,0.7838115,1.033258); gStyle->SetOptFit(0111); can5->SetFillColor(10); can5->SetBorderMode(0); can5->SetBorderSize(2); can5->SetGridx(); can5->SetGridy(); can5->SetFrameFillColor(0); can5->SetFrameBorderMode(0); can5->SetFrameBorderMode(0); gPad->SetRightMargin(0.02); gPad->SetTopMargin(0.02);*/ r42 = (TH1D*)n4QS->Clone(); TPN_FourParticle[ch1_4]->Multiply(MRC_SC_4[4]); r42->Divide(TPN_FourParticle[ch1_4]); r42->GetXaxis()->SetRangeUser(0,0.08); r42->SetMinimum(0.5); r42->SetMaximum(14); r42->GetYaxis()->SetTitle("r_{4}{2}"); // //r42->Draw(); //ChaoticLimit_r42->Draw("same"); //TPN_FourParticle[ch1_4]->Draw(); //fTPNRejects4pion1->SetLineColor(2); //fTPNRejects4pion1->Draw("same"); } //double EA = exp(-pow(0.005 * 10/FmToGeV,2)/2.); //TF1 *C2mag = new TF1("C2mag","pow(1-x,2)*pow([0],2) + 2*x*(1-x)*[0]",0,1); //C2mag->FixParameter(0, EA); //TF1 *C4norm = new TF1("C4norm","( 6*pow(1-x,4)*pow([0],4) + 24*x*pow(1-x,3)*pow([0],3) + 4*(2*pow(1-x,3)*pow([0],3) + 6*x*pow(1-x,2)*pow([0],2)) + 3*pow(C2mag,2) + 6*C2mag + 1) / (6*pow(C2mag,2) + 8*pow(C2mag,1.5) + 3*pow(C2mag,2) + 6*C2mag + 1)",0,1); //C4norm->FixParameter(0, EA); //C4norm->Draw(); /*TF1 *C2mag = new TF1("C2mag","pow(1-[0],2)*exp(-pow(x*[1],2)/6.) + 2*[0]*(1-[0])*exp(-pow(x*[1],2)/12.)",0,0.2); C2mag->FixParameter(0, 0.25); C2mag->FixParameter(1, 10./FmToGeV); TF1 *C4norm = new TF1("C4norm","( 6*pow(1-[0],4)*exp(-pow(x*[1],2)/3.) + 24*[0]*pow(1-[0],3)*exp(-pow(x*[1],2)/4.) + 4*(2*pow(1-[0],3)*exp(-pow(x*[1],2)/4.) + 6*[0]*pow(1-[0],2)*exp(-pow(x*[1],2)/6.)) + 3*pow(C2mag,2) + 6*C2mag + 1) / (6*pow(C2mag,2) + 8*pow(C2mag,1.5) + 3*pow(C2mag,2) + 6*C2mag + 1)",0,0.2); C4norm->FixParameter(0, 0.25); C4norm->FixParameter(1, 10./FmToGeV);*/ //C4norm->Draw(); if(SaveToFile){ TString *savefilename = new TString("OutFiles/OutFile"); if(MCcase) savefilename->Append("MonteCarlo"); // if(SameCharge) savefilename->Append("SC"); else savefilename->Append("MC"); // if(!SameCharge) *savefilename += MixedCharge4pionType_def; // if(CHARGE==1) savefilename->Append("_Pos_"); else savefilename->Append("_Neg_"); // savefilename->Append("KT_"); *savefilename += EDbin+1; savefilename->Append("_M"); *savefilename += Mbin; savefilename->Append(".root"); TFile *savefile = new TFile(savefilename->Data(),"RECREATE"); // C2->Write("C2"); C2QS->Write("C2QS"); C3QS->Write("C3QS"); c3QS->Write("c3QS"); C4QS->Write("C4QS"); c4QS->Write("c4QS"); c4QS_RemovalStage1->Write("c4QS_RemovalStage1"); if(SameCharge) { r3->Write("r3"); r42->Write("r42"); c4QS_RemovalStage2->Write("c4QS_RemovalStage2"); TPFullWeight_ThreeParticle[ch1_3]->Write("C3QS_built"); TPFullWeight_FourParticle[ch1_4]->Write("C4QS_built"); // TPFullWeight_ThreeParticle_2D[ch1_3]->Write("C3QS_built2D"); TPFullWeight_FourParticle_2D[ch1_4]->Write("C4QS_built2D"); TPNegFullWeight_ThreeParticle_2D[ch1_3]->Write("C3QS_Negbuilt2D"); TPNegFullWeight_FourParticle_2D[ch1_4]->Write("C4QS_Negbuilt2D"); } // savefile->Close(); } } //________________________________________________________________________ void SetFSICorrelations(){ // read in 2-particle and 3-particle FSI correlations = K2 & K3 // 2-particle input histo from file is binned in qinv. 3-particle in qinv of each pair TFile *fsifile = new TFile("KFile.root","READ"); if(!fsifile->IsOpen()) { cout<<"No FSI file found!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"<<endl; } TH1D *temphistoSS[12]; TH1D *temphistoOS[12]; for(Int_t MB=0; MB<12; MB++) { TString *nameK2SS = new TString("K2ss_"); *nameK2SS += MB; temphistoSS[MB] = (TH1D*)fsifile->Get(nameK2SS->Data()); // TString *nameK2OS = new TString("K2os_"); *nameK2OS += MB; temphistoOS[MB] = (TH1D*)fsifile->Get(nameK2OS->Data()); // fFSIss[MB] = (TH1D*)temphistoSS[MB]->Clone(); fFSIos[MB] = (TH1D*)temphistoOS[MB]->Clone(); fFSIss[MB]->SetDirectory(0); fFSIos[MB]->SetDirectory(0); } // fsifile->Close(); } double Gamov(int chargeProduct, double qinv){ double arg = chargeProduct*2.*PI/(BohrR*qinv/2.); return arg/(exp(arg)-1); } void SetMomResCorrections(){ TString *momresfilename = new TString("MomResFile"); if(FileSetting==7) momresfilename->Append("_10percentIncrease"); if(CollisionType!=0) momresfilename->Append("_ppAndpPb"); momresfilename->Append(".root"); TFile *MomResFile = new TFile(momresfilename->Data(),"READ"); TString *proName[28]; for(int ii=0; ii<28; ii++){ proName[ii] = new TString("MRC_pro_"); *proName[ii] += ii; } MRC_SC_2[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_2_SC_term1"))->ProjectionY(proName[0]->Data(), RbinMRC, RbinMRC)); MRC_SC_2[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_2_SC_term2"))->ProjectionY(proName[1]->Data(), RbinMRC, RbinMRC)); MRC_SC_2[0]->SetDirectory(0); MRC_SC_2[1]->SetDirectory(0); // MRC_MC_2[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_2_MC_term1"))->ProjectionY(proName[2]->Data(), RbinMRC, RbinMRC)); MRC_MC_2[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_2_MC_term2"))->ProjectionY(proName[3]->Data(), RbinMRC, RbinMRC)); MRC_MC_2[0]->SetDirectory(0); MRC_MC_2[1]->SetDirectory(0); // MRC_SC_3[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_SC_term1"))->ProjectionY(proName[4]->Data(), RbinMRC, RbinMRC)); MRC_SC_3[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_SC_term2"))->ProjectionY(proName[5]->Data(), RbinMRC, RbinMRC)); MRC_SC_3[2] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_SC_term3"))->ProjectionY(proName[6]->Data(), RbinMRC, RbinMRC)); MRC_SC_3[0]->SetDirectory(0); MRC_SC_3[1]->SetDirectory(0); MRC_SC_3[2]->SetDirectory(0); // MRC_MC_3[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_MC_term1"))->ProjectionY(proName[7]->Data(), RbinMRC, RbinMRC)); MRC_MC_3[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_MC_term2"))->ProjectionY(proName[8]->Data(), RbinMRC, RbinMRC)); MRC_MC_3[2] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_MC_term3"))->ProjectionY(proName[9]->Data(), RbinMRC, RbinMRC)); MRC_MC_3[3] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_3_MC_term4"))->ProjectionY(proName[10]->Data(), RbinMRC, RbinMRC)); MRC_MC_3[0]->SetDirectory(0); MRC_MC_3[1]->SetDirectory(0); MRC_MC_3[2]->SetDirectory(0); MRC_MC_3[3]->SetDirectory(0); // MRC_SC_4[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_SC_term1"))->ProjectionY(proName[11]->Data(), RbinMRC, RbinMRC)); MRC_SC_4[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_SC_term2"))->ProjectionY(proName[12]->Data(), RbinMRC, RbinMRC)); MRC_SC_4[2] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_SC_term3"))->ProjectionY(proName[13]->Data(), RbinMRC, RbinMRC)); MRC_SC_4[3] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_SC_term4"))->ProjectionY(proName[14]->Data(), RbinMRC, RbinMRC)); MRC_SC_4[4] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_SC_term5"))->ProjectionY(proName[15]->Data(), RbinMRC, RbinMRC)); MRC_SC_4[0]->SetDirectory(0); MRC_SC_4[1]->SetDirectory(0); MRC_SC_4[2]->SetDirectory(0); MRC_SC_4[3]->SetDirectory(0); MRC_SC_4[4]->SetDirectory(0); // MRC_MC1_4[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term1"))->ProjectionY(proName[16]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term2"))->ProjectionY(proName[17]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[2] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term3"))->ProjectionY(proName[18]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[3] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term4"))->ProjectionY(proName[19]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[4] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term5"))->ProjectionY(proName[20]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[5] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC1_term6"))->ProjectionY(proName[21]->Data(), RbinMRC, RbinMRC)); MRC_MC1_4[0]->SetDirectory(0); MRC_MC1_4[1]->SetDirectory(0); MRC_MC1_4[2]->SetDirectory(0); MRC_MC1_4[3]->SetDirectory(0); MRC_MC1_4[4]->SetDirectory(0); MRC_MC1_4[5]->SetDirectory(0); // MRC_MC2_4[0] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term1"))->ProjectionY(proName[22]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[1] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term2"))->ProjectionY(proName[23]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[2] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term3"))->ProjectionY(proName[24]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[3] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term4"))->ProjectionY(proName[25]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[4] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term5"))->ProjectionY(proName[26]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[5] = (TH1D*)(((TH2D*)MomResFile->Get("MRC_4_MC2_term6"))->ProjectionY(proName[27]->Data(), RbinMRC, RbinMRC)); MRC_MC2_4[0]->SetDirectory(0); MRC_MC2_4[1]->SetDirectory(0); MRC_MC2_4[2]->SetDirectory(0); MRC_MC2_4[3]->SetDirectory(0); MRC_MC2_4[4]->SetDirectory(0); MRC_MC2_4[5]->SetDirectory(0); if(!MRC || MCcase_def){ for(int bin=1; bin<=MRC_SC_2[0]->GetNbinsX(); bin++){ MRC_SC_2[0]->SetBinContent(bin, 1.0); MRC_SC_2[1]->SetBinContent(bin, 1.0); MRC_MC_2[0]->SetBinContent(bin, 1.0); MRC_MC_2[1]->SetBinContent(bin, 1.0); } for(int bin=1; bin<=MRC_SC_3[0]->GetNbinsX(); bin++){ MRC_SC_3[0]->SetBinContent(bin, 1.0); MRC_SC_3[1]->SetBinContent(bin, 1.0); MRC_SC_3[2]->SetBinContent(bin, 1.0); MRC_MC_3[0]->SetBinContent(bin, 1.0); MRC_MC_3[1]->SetBinContent(bin, 1.0); MRC_MC_3[2]->SetBinContent(bin, 1.0); MRC_MC_3[3]->SetBinContent(bin, 1.0); } for(int bin=1; bin<=MRC_SC_4[0]->GetNbinsX(); bin++){ MRC_SC_4[0]->SetBinContent(bin, 1.0); MRC_SC_4[1]->SetBinContent(bin, 1.0); MRC_SC_4[2]->SetBinContent(bin, 1.0); MRC_SC_4[3]->SetBinContent(bin, 1.0); MRC_SC_4[4]->SetBinContent(bin, 1.0); MRC_MC1_4[0]->SetBinContent(bin, 1.0); MRC_MC1_4[1]->SetBinContent(bin, 1.0); MRC_MC1_4[2]->SetBinContent(bin, 1.0); MRC_MC1_4[3]->SetBinContent(bin, 1.0); MRC_MC1_4[4]->SetBinContent(bin, 1.0); MRC_MC1_4[5]->SetBinContent(bin, 1.0); // MRC_MC2_4[0]->SetBinContent(bin, 1.0); MRC_MC2_4[1]->SetBinContent(bin, 1.0); MRC_MC2_4[2]->SetBinContent(bin, 1.0); MRC_MC2_4[3]->SetBinContent(bin, 1.0); MRC_MC2_4[4]->SetBinContent(bin, 1.0); MRC_MC2_4[5]->SetBinContent(bin, 1.0); } } MomResFile->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); } float fact(float n){ return (n < 1.00001) ? 1 : fact(n - 1) * n; } //________________________________________________________________________________________ void SetMuonCorrections(){ TString *name = new TString("MuonCorrection"); if(FileSetting==8) name->Append("_92percent"); if(CollisionType!=0) name->Append("_ppAndpPb"); name->Append(".root"); TFile *MuonFile=new TFile(name->Data(),"READ"); TString *proName[22]; for(int ii=0; ii<22; ii++){ proName[ii] = new TString("MuonCorr_pro_"); *proName[ii] += ii; } C2muonCorrectionSC = (TH1D*)(((TH2D*)MuonFile->Get("C2muonCorrection_SC"))->ProjectionY(proName[0]->Data(), RbinMRC, RbinMRC)); C2muonCorrectionMC = (TH1D*)(((TH2D*)MuonFile->Get("C2muonCorrection_MC"))->ProjectionY(proName[1]->Data(), RbinMRC, RbinMRC)); WeightmuonCorrection = (TH1D*)(((TH2D*)MuonFile->Get("WeightmuonCorrection"))->ProjectionY(proName[2]->Data(), RbinMRC, RbinMRC)); C2muonCorrectionSC->SetDirectory(0); C2muonCorrectionMC->SetDirectory(0); WeightmuonCorrection->SetDirectory(0); // C3muonCorrectionSC[0] = (TH1D*)(((TH2D*)MuonFile->Get("C3muonCorrection_SC_term1"))->ProjectionY(proName[3]->Data(), RbinMRC, RbinMRC)); C3muonCorrectionSC[1] = (TH1D*)(((TH2D*)MuonFile->Get("C3muonCorrection_SC_term2"))->ProjectionY(proName[4]->Data(), RbinMRC, RbinMRC)); C3muonCorrectionMC[0] = (TH1D*)(((TH2D*)MuonFile->Get("C3muonCorrection_MC_term1"))->ProjectionY(proName[5]->Data(), RbinMRC, RbinMRC)); C3muonCorrectionMC[1] = (TH1D*)(((TH2D*)MuonFile->Get("C3muonCorrection_MC_term2"))->ProjectionY(proName[6]->Data(), RbinMRC, RbinMRC)); C3muonCorrectionMC[2] = (TH1D*)(((TH2D*)MuonFile->Get("C3muonCorrection_MC_term3"))->ProjectionY(proName[7]->Data(), RbinMRC, RbinMRC)); C3muonCorrectionSC[0]->SetDirectory(0); C3muonCorrectionSC[1]->SetDirectory(0); C3muonCorrectionMC[0]->SetDirectory(0); C3muonCorrectionMC[1]->SetDirectory(0); C3muonCorrectionMC[2]->SetDirectory(0); // C4muonCorrectionSC[0] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_SC_term1"))->ProjectionY(proName[8]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionSC[1] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_SC_term2"))->ProjectionY(proName[9]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionSC[2] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_SC_term3"))->ProjectionY(proName[10]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionSC[3] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_SC_term4"))->ProjectionY(proName[11]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionSC[0]->SetDirectory(0); C4muonCorrectionSC[1]->SetDirectory(0); C4muonCorrectionSC[2]->SetDirectory(0); C4muonCorrectionSC[3]->SetDirectory(0); // C4muonCorrectionMC1[0] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC1_term1"))->ProjectionY(proName[12]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC1[1] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC1_term2"))->ProjectionY(proName[13]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC1[2] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC1_term3"))->ProjectionY(proName[14]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC1[3] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC1_term4"))->ProjectionY(proName[15]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC1[4] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC1_term5"))->ProjectionY(proName[16]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC1[0]->SetDirectory(0); C4muonCorrectionMC1[1]->SetDirectory(0); C4muonCorrectionMC1[2]->SetDirectory(0); C4muonCorrectionMC1[3]->SetDirectory(0); C4muonCorrectionMC1[4]->SetDirectory(0); // C4muonCorrectionMC2[0] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC2_term1"))->ProjectionY(proName[17]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC2[1] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC2_term2"))->ProjectionY(proName[18]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC2[2] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC2_term3"))->ProjectionY(proName[19]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC2[3] = (TH1D*)(((TH2D*)MuonFile->Get("C4muonCorrection_MC2_term4"))->ProjectionY(proName[20]->Data(), RbinMRC, RbinMRC)); C4muonCorrectionMC2[4] = (TH1D*)C4muonCorrectionSC[3]->Clone(); C4muonCorrectionMC2[0]->SetDirectory(0); C4muonCorrectionMC2[1]->SetDirectory(0); C4muonCorrectionMC2[2]->SetDirectory(0); C4muonCorrectionMC2[3]->SetDirectory(0); C4muonCorrectionMC2[4]->SetDirectory(0); // // if(!MuonCorrection || MCcase_def){ for(int bin=1; bin<=C2muonCorrectionSC->GetNbinsX(); bin++){ C2muonCorrectionSC->SetBinContent(bin, 1.0); C2muonCorrectionMC->SetBinContent(bin, 1.0); WeightmuonCorrection->SetBinContent(bin, 1.0); } for(int bin=1; bin<=C3muonCorrectionSC[0]->GetNbinsX(); bin++){ C3muonCorrectionSC[0]->SetBinContent(bin, 1.0); C3muonCorrectionSC[1]->SetBinContent(bin, 1.0); C3muonCorrectionMC[0]->SetBinContent(bin, 1.0); C3muonCorrectionMC[1]->SetBinContent(bin, 1.0); C3muonCorrectionMC[2]->SetBinContent(bin, 1.0); } for(int bin=1; bin<=C4muonCorrectionSC[0]->GetNbinsX(); bin++){ C4muonCorrectionSC[0]->SetBinContent(bin, 1.0); C4muonCorrectionSC[1]->SetBinContent(bin, 1.0); C4muonCorrectionSC[2]->SetBinContent(bin, 1.0); C4muonCorrectionSC[3]->SetBinContent(bin, 1.0); C4muonCorrectionMC1[0]->SetBinContent(bin, 1.0); C4muonCorrectionMC1[1]->SetBinContent(bin, 1.0); C4muonCorrectionMC1[2]->SetBinContent(bin, 1.0); C4muonCorrectionMC1[3]->SetBinContent(bin, 1.0); C4muonCorrectionMC1[4]->SetBinContent(bin, 1.0); // C4muonCorrectionMC2[0]->SetBinContent(bin, 1.0); C4muonCorrectionMC2[1]->SetBinContent(bin, 1.0); C4muonCorrectionMC2[2]->SetBinContent(bin, 1.0); C4muonCorrectionMC2[3]->SetBinContent(bin, 1.0); C4muonCorrectionMC2[4]->SetBinContent(bin, 1.0); } } MuonFile->Close(); } //________________________________________________________________________ void SetFSIindex(Float_t R){ if(!MCcase_def){ if(CollisionType==0){ if(Mbin_def==0) fFSIindex = 0;//0-5% else if(Mbin_def==1) fFSIindex = 1;//5-10% else if(Mbin_def<=3) fFSIindex = 2;//10-20% else if(Mbin_def<=5) fFSIindex = 3;//20-30% else if(Mbin_def<=7) fFSIindex = 4;//30-40% else if(Mbin_def<=9) fFSIindex = 5;//40-50% else if(Mbin_def<=12) fFSIindex = 6;//40-50% else if(Mbin_def<=15) fFSIindex = 7;//40-50% else if(Mbin_def<=18) fFSIindex = 8;//40-50% else fFSIindex = 8;//90-100% }else fFSIindex = 9;// pp and pPb }else{// FSI binning for MC if(R>=10.) fFSIindex = 0; else if(R>=9.) fFSIindex = 1; else if(R>=8.) fFSIindex = 2; else if(R>=7.) fFSIindex = 3; else if(R>=6.) fFSIindex = 4; else if(R>=5.) fFSIindex = 5; else if(R>=4.) fFSIindex = 6; else if(R>=3.) fFSIindex = 7; else if(R>=2.) fFSIindex = 8; else fFSIindex = 9; } } //________________________________________________________________________ Float_t FSICorrelation(Int_t charge1, Int_t charge2, Float_t qinv){ // returns 2-particle Coulomb correlations = K2 Int_t qbinL = fFSIss[fFSIindex]->GetXaxis()->FindBin(qinv-fFSIss[fFSIindex]->GetXaxis()->GetBinWidth(1)/2.); Int_t qbinH = qbinL+1; if(qbinL <= 0) return 1.0; if(qbinH > fFSIss[fFSIindex]->GetNbinsX()) {// Scaled Gamov approximation int chargeproduct = 1; if(charge1!=charge2) { chargeproduct = -1; Float_t ScaleFac = (fFSIos[fFSIindex]->GetBinContent(fFSIos[fFSIindex]->GetNbinsX()-1) - 1); ScaleFac /= (Gamov(chargeproduct, fFSIos[fFSIindex]->GetXaxis()->GetBinCenter(fFSIos[fFSIindex]->GetNbinsX()-1)) - 1); return ( (Gamov(chargeproduct, qinv)-1)*ScaleFac + 1); }else{ Float_t ScaleFac = (fFSIss[fFSIindex]->GetBinContent(fFSIss[fFSIindex]->GetNbinsX()-1) - 1); ScaleFac /= (Gamov(chargeproduct, fFSIss[fFSIindex]->GetXaxis()->GetBinCenter(fFSIss[fFSIindex]->GetNbinsX()-1)) - 1); return ( (Gamov(chargeproduct, qinv)-1)*ScaleFac + 1); } } Float_t slope=0; if(charge1==charge2){ slope = fFSIss[fFSIindex]->GetBinContent(qbinL) - fFSIss[fFSIindex]->GetBinContent(qbinH); slope /= fFSIss[fFSIindex]->GetXaxis()->GetBinCenter(qbinL) - fFSIss[fFSIindex]->GetXaxis()->GetBinCenter(qbinH); return (slope*(qinv - fFSIss[fFSIindex]->GetXaxis()->GetBinCenter(qbinL)) + fFSIss[fFSIindex]->GetBinContent(qbinL)); }else { slope = fFSIos[fFSIindex]->GetBinContent(qbinL) - fFSIos[fFSIindex]->GetBinContent(qbinH); slope /= fFSIos[fFSIindex]->GetXaxis()->GetBinCenter(qbinL) - fFSIos[fFSIindex]->GetXaxis()->GetBinCenter(qbinH); return (slope*(qinv - fFSIos[fFSIindex]->GetXaxis()->GetBinCenter(qbinL)) + fFSIos[fFSIindex]->GetBinContent(qbinL)); } }
Plot_FourPion.C:1
Plot_FourPion.C:2
Plot_FourPion.C:3
Plot_FourPion.C:4
Plot_FourPion.C:5
Plot_FourPion.C:6
Plot_FourPion.C:7
Plot_FourPion.C:8
Plot_FourPion.C:9
Plot_FourPion.C:10
Plot_FourPion.C:11
Plot_FourPion.C:12
Plot_FourPion.C:13
Plot_FourPion.C:14
Plot_FourPion.C:15
Plot_FourPion.C:16
Plot_FourPion.C:17
Plot_FourPion.C:18
Plot_FourPion.C:19
Plot_FourPion.C:20
Plot_FourPion.C:21
Plot_FourPion.C:22
Plot_FourPion.C:23
Plot_FourPion.C:24
Plot_FourPion.C:25
Plot_FourPion.C:26
Plot_FourPion.C:27
Plot_FourPion.C:28
Plot_FourPion.C:29
Plot_FourPion.C:30
Plot_FourPion.C:31
Plot_FourPion.C:32
Plot_FourPion.C:33
Plot_FourPion.C:34
Plot_FourPion.C:35
Plot_FourPion.C:36
Plot_FourPion.C:37
Plot_FourPion.C:38
Plot_FourPion.C:39
Plot_FourPion.C:40
Plot_FourPion.C:41
Plot_FourPion.C:42
Plot_FourPion.C:43
Plot_FourPion.C:44
Plot_FourPion.C:45
Plot_FourPion.C:46
Plot_FourPion.C:47
Plot_FourPion.C:48
Plot_FourPion.C:49
Plot_FourPion.C:50
Plot_FourPion.C:51
Plot_FourPion.C:52
Plot_FourPion.C:53
Plot_FourPion.C:54
Plot_FourPion.C:55
Plot_FourPion.C:56
Plot_FourPion.C:57
Plot_FourPion.C:58
Plot_FourPion.C:59
Plot_FourPion.C:60
Plot_FourPion.C:61
Plot_FourPion.C:62
Plot_FourPion.C:63
Plot_FourPion.C:64
Plot_FourPion.C:65
Plot_FourPion.C:66
Plot_FourPion.C:67
Plot_FourPion.C:68
Plot_FourPion.C:69
Plot_FourPion.C:70
Plot_FourPion.C:71
Plot_FourPion.C:72
Plot_FourPion.C:73
Plot_FourPion.C:74
Plot_FourPion.C:75
Plot_FourPion.C:76
Plot_FourPion.C:77
Plot_FourPion.C:78
Plot_FourPion.C:79
Plot_FourPion.C:80
Plot_FourPion.C:81
Plot_FourPion.C:82
Plot_FourPion.C:83
Plot_FourPion.C:84
Plot_FourPion.C:85
Plot_FourPion.C:86
Plot_FourPion.C:87
Plot_FourPion.C:88
Plot_FourPion.C:89
Plot_FourPion.C:90
Plot_FourPion.C:91
Plot_FourPion.C:92
Plot_FourPion.C:93
Plot_FourPion.C:94
Plot_FourPion.C:95
Plot_FourPion.C:96
Plot_FourPion.C:97
Plot_FourPion.C:98
Plot_FourPion.C:99
Plot_FourPion.C:100
Plot_FourPion.C:101
Plot_FourPion.C:102
Plot_FourPion.C:103
Plot_FourPion.C:104
Plot_FourPion.C:105
Plot_FourPion.C:106
Plot_FourPion.C:107
Plot_FourPion.C:108
Plot_FourPion.C:109
Plot_FourPion.C:110
Plot_FourPion.C:111
Plot_FourPion.C:112
Plot_FourPion.C:113
Plot_FourPion.C:114
Plot_FourPion.C:115
Plot_FourPion.C:116
Plot_FourPion.C:117
Plot_FourPion.C:118
Plot_FourPion.C:119
Plot_FourPion.C:120
Plot_FourPion.C:121
Plot_FourPion.C:122
Plot_FourPion.C:123
Plot_FourPion.C:124
Plot_FourPion.C:125
Plot_FourPion.C:126
Plot_FourPion.C:127
Plot_FourPion.C:128
Plot_FourPion.C:129
Plot_FourPion.C:130
Plot_FourPion.C:131
Plot_FourPion.C:132
Plot_FourPion.C:133
Plot_FourPion.C:134
Plot_FourPion.C:135
Plot_FourPion.C:136
Plot_FourPion.C:137
Plot_FourPion.C:138
Plot_FourPion.C:139
Plot_FourPion.C:140
Plot_FourPion.C:141
Plot_FourPion.C:142
Plot_FourPion.C:143
Plot_FourPion.C:144
Plot_FourPion.C:145
Plot_FourPion.C:146
Plot_FourPion.C:147
Plot_FourPion.C:148
Plot_FourPion.C:149
Plot_FourPion.C:150
Plot_FourPion.C:151
Plot_FourPion.C:152
Plot_FourPion.C:153
Plot_FourPion.C:154
Plot_FourPion.C:155
Plot_FourPion.C:156
Plot_FourPion.C:157
Plot_FourPion.C:158
Plot_FourPion.C:159
Plot_FourPion.C:160
Plot_FourPion.C:161
Plot_FourPion.C:162
Plot_FourPion.C:163
Plot_FourPion.C:164
Plot_FourPion.C:165
Plot_FourPion.C:166
Plot_FourPion.C:167
Plot_FourPion.C:168
Plot_FourPion.C:169
Plot_FourPion.C:170
Plot_FourPion.C:171
Plot_FourPion.C:172
Plot_FourPion.C:173
Plot_FourPion.C:174
Plot_FourPion.C:175
Plot_FourPion.C:176
Plot_FourPion.C:177
Plot_FourPion.C:178
Plot_FourPion.C:179
Plot_FourPion.C:180
Plot_FourPion.C:181
Plot_FourPion.C:182
Plot_FourPion.C:183
Plot_FourPion.C:184
Plot_FourPion.C:185
Plot_FourPion.C:186
Plot_FourPion.C:187
Plot_FourPion.C:188
Plot_FourPion.C:189
Plot_FourPion.C:190
Plot_FourPion.C:191
Plot_FourPion.C:192
Plot_FourPion.C:193
Plot_FourPion.C:194
Plot_FourPion.C:195
Plot_FourPion.C:196
Plot_FourPion.C:197
Plot_FourPion.C:198
Plot_FourPion.C:199
Plot_FourPion.C:200
Plot_FourPion.C:201
Plot_FourPion.C:202
Plot_FourPion.C:203
Plot_FourPion.C:204
Plot_FourPion.C:205
Plot_FourPion.C:206
Plot_FourPion.C:207
Plot_FourPion.C:208
Plot_FourPion.C:209
Plot_FourPion.C:210
Plot_FourPion.C:211
Plot_FourPion.C:212
Plot_FourPion.C:213
Plot_FourPion.C:214
Plot_FourPion.C:215
Plot_FourPion.C:216
Plot_FourPion.C:217
Plot_FourPion.C:218
Plot_FourPion.C:219
Plot_FourPion.C:220
Plot_FourPion.C:221
Plot_FourPion.C:222
Plot_FourPion.C:223
Plot_FourPion.C:224
Plot_FourPion.C:225
Plot_FourPion.C:226
Plot_FourPion.C:227
Plot_FourPion.C:228
Plot_FourPion.C:229
Plot_FourPion.C:230
Plot_FourPion.C:231
Plot_FourPion.C:232
Plot_FourPion.C:233
Plot_FourPion.C:234
Plot_FourPion.C:235
Plot_FourPion.C:236
Plot_FourPion.C:237
Plot_FourPion.C:238
Plot_FourPion.C:239
Plot_FourPion.C:240
Plot_FourPion.C:241
Plot_FourPion.C:242
Plot_FourPion.C:243
Plot_FourPion.C:244
Plot_FourPion.C:245
Plot_FourPion.C:246
Plot_FourPion.C:247
Plot_FourPion.C:248
Plot_FourPion.C:249
Plot_FourPion.C:250
Plot_FourPion.C:251
Plot_FourPion.C:252
Plot_FourPion.C:253
Plot_FourPion.C:254
Plot_FourPion.C:255
Plot_FourPion.C:256
Plot_FourPion.C:257
Plot_FourPion.C:258
Plot_FourPion.C:259
Plot_FourPion.C:260
Plot_FourPion.C:261
Plot_FourPion.C:262
Plot_FourPion.C:263
Plot_FourPion.C:264
Plot_FourPion.C:265
Plot_FourPion.C:266
Plot_FourPion.C:267
Plot_FourPion.C:268
Plot_FourPion.C:269
Plot_FourPion.C:270
Plot_FourPion.C:271
Plot_FourPion.C:272
Plot_FourPion.C:273
Plot_FourPion.C:274
Plot_FourPion.C:275
Plot_FourPion.C:276
Plot_FourPion.C:277
Plot_FourPion.C:278
Plot_FourPion.C:279
Plot_FourPion.C:280
Plot_FourPion.C:281
Plot_FourPion.C:282
Plot_FourPion.C:283
Plot_FourPion.C:284
Plot_FourPion.C:285
Plot_FourPion.C:286
Plot_FourPion.C:287
Plot_FourPion.C:288
Plot_FourPion.C:289
Plot_FourPion.C:290
Plot_FourPion.C:291
Plot_FourPion.C:292
Plot_FourPion.C:293
Plot_FourPion.C:294
Plot_FourPion.C:295
Plot_FourPion.C:296
Plot_FourPion.C:297
Plot_FourPion.C:298
Plot_FourPion.C:299
Plot_FourPion.C:300
Plot_FourPion.C:301
Plot_FourPion.C:302
Plot_FourPion.C:303
Plot_FourPion.C:304
Plot_FourPion.C:305
Plot_FourPion.C:306
Plot_FourPion.C:307
Plot_FourPion.C:308
Plot_FourPion.C:309
Plot_FourPion.C:310
Plot_FourPion.C:311
Plot_FourPion.C:312
Plot_FourPion.C:313
Plot_FourPion.C:314
Plot_FourPion.C:315
Plot_FourPion.C:316
Plot_FourPion.C:317
Plot_FourPion.C:318
Plot_FourPion.C:319
Plot_FourPion.C:320
Plot_FourPion.C:321
Plot_FourPion.C:322
Plot_FourPion.C:323
Plot_FourPion.C:324
Plot_FourPion.C:325
Plot_FourPion.C:326
Plot_FourPion.C:327
Plot_FourPion.C:328
Plot_FourPion.C:329
Plot_FourPion.C:330
Plot_FourPion.C:331
Plot_FourPion.C:332
Plot_FourPion.C:333
Plot_FourPion.C:334
Plot_FourPion.C:335
Plot_FourPion.C:336
Plot_FourPion.C:337
Plot_FourPion.C:338
Plot_FourPion.C:339
Plot_FourPion.C:340
Plot_FourPion.C:341
Plot_FourPion.C:342
Plot_FourPion.C:343
Plot_FourPion.C:344
Plot_FourPion.C:345
Plot_FourPion.C:346
Plot_FourPion.C:347
Plot_FourPion.C:348
Plot_FourPion.C:349
Plot_FourPion.C:350
Plot_FourPion.C:351
Plot_FourPion.C:352
Plot_FourPion.C:353
Plot_FourPion.C:354
Plot_FourPion.C:355
Plot_FourPion.C:356
Plot_FourPion.C:357
Plot_FourPion.C:358
Plot_FourPion.C:359
Plot_FourPion.C:360
Plot_FourPion.C:361
Plot_FourPion.C:362
Plot_FourPion.C:363
Plot_FourPion.C:364
Plot_FourPion.C:365
Plot_FourPion.C:366
Plot_FourPion.C:367
Plot_FourPion.C:368
Plot_FourPion.C:369
Plot_FourPion.C:370
Plot_FourPion.C:371
Plot_FourPion.C:372
Plot_FourPion.C:373
Plot_FourPion.C:374
Plot_FourPion.C:375
Plot_FourPion.C:376
Plot_FourPion.C:377
Plot_FourPion.C:378
Plot_FourPion.C:379
Plot_FourPion.C:380
Plot_FourPion.C:381
Plot_FourPion.C:382
Plot_FourPion.C:383
Plot_FourPion.C:384
Plot_FourPion.C:385
Plot_FourPion.C:386
Plot_FourPion.C:387
Plot_FourPion.C:388
Plot_FourPion.C:389
Plot_FourPion.C:390
Plot_FourPion.C:391
Plot_FourPion.C:392
Plot_FourPion.C:393
Plot_FourPion.C:394
Plot_FourPion.C:395
Plot_FourPion.C:396
Plot_FourPion.C:397
Plot_FourPion.C:398
Plot_FourPion.C:399
Plot_FourPion.C:400
Plot_FourPion.C:401
Plot_FourPion.C:402
Plot_FourPion.C:403
Plot_FourPion.C:404
Plot_FourPion.C:405
Plot_FourPion.C:406
Plot_FourPion.C:407
Plot_FourPion.C:408
Plot_FourPion.C:409
Plot_FourPion.C:410
Plot_FourPion.C:411
Plot_FourPion.C:412
Plot_FourPion.C:413
Plot_FourPion.C:414
Plot_FourPion.C:415
Plot_FourPion.C:416
Plot_FourPion.C:417
Plot_FourPion.C:418
Plot_FourPion.C:419
Plot_FourPion.C:420
Plot_FourPion.C:421
Plot_FourPion.C:422
Plot_FourPion.C:423
Plot_FourPion.C:424
Plot_FourPion.C:425
Plot_FourPion.C:426
Plot_FourPion.C:427
Plot_FourPion.C:428
Plot_FourPion.C:429
Plot_FourPion.C:430
Plot_FourPion.C:431
Plot_FourPion.C:432
Plot_FourPion.C:433
Plot_FourPion.C:434
Plot_FourPion.C:435
Plot_FourPion.C:436
Plot_FourPion.C:437
Plot_FourPion.C:438
Plot_FourPion.C:439
Plot_FourPion.C:440
Plot_FourPion.C:441
Plot_FourPion.C:442
Plot_FourPion.C:443
Plot_FourPion.C:444
Plot_FourPion.C:445
Plot_FourPion.C:446
Plot_FourPion.C:447
Plot_FourPion.C:448
Plot_FourPion.C:449
Plot_FourPion.C:450
Plot_FourPion.C:451
Plot_FourPion.C:452
Plot_FourPion.C:453
Plot_FourPion.C:454
Plot_FourPion.C:455
Plot_FourPion.C:456
Plot_FourPion.C:457
Plot_FourPion.C:458
Plot_FourPion.C:459
Plot_FourPion.C:460
Plot_FourPion.C:461
Plot_FourPion.C:462
Plot_FourPion.C:463
Plot_FourPion.C:464
Plot_FourPion.C:465
Plot_FourPion.C:466
Plot_FourPion.C:467
Plot_FourPion.C:468
Plot_FourPion.C:469
Plot_FourPion.C:470
Plot_FourPion.C:471
Plot_FourPion.C:472
Plot_FourPion.C:473
Plot_FourPion.C:474
Plot_FourPion.C:475
Plot_FourPion.C:476
Plot_FourPion.C:477
Plot_FourPion.C:478
Plot_FourPion.C:479
Plot_FourPion.C:480
Plot_FourPion.C:481
Plot_FourPion.C:482
Plot_FourPion.C:483
Plot_FourPion.C:484
Plot_FourPion.C:485
Plot_FourPion.C:486
Plot_FourPion.C:487
Plot_FourPion.C:488
Plot_FourPion.C:489
Plot_FourPion.C:490
Plot_FourPion.C:491
Plot_FourPion.C:492
Plot_FourPion.C:493
Plot_FourPion.C:494
Plot_FourPion.C:495
Plot_FourPion.C:496
Plot_FourPion.C:497
Plot_FourPion.C:498
Plot_FourPion.C:499
Plot_FourPion.C:500
Plot_FourPion.C:501
Plot_FourPion.C:502
Plot_FourPion.C:503
Plot_FourPion.C:504
Plot_FourPion.C:505
Plot_FourPion.C:506
Plot_FourPion.C:507
Plot_FourPion.C:508
Plot_FourPion.C:509
Plot_FourPion.C:510
Plot_FourPion.C:511
Plot_FourPion.C:512
Plot_FourPion.C:513
Plot_FourPion.C:514
Plot_FourPion.C:515
Plot_FourPion.C:516
Plot_FourPion.C:517
Plot_FourPion.C:518
Plot_FourPion.C:519
Plot_FourPion.C:520
Plot_FourPion.C:521
Plot_FourPion.C:522
Plot_FourPion.C:523
Plot_FourPion.C:524
Plot_FourPion.C:525
Plot_FourPion.C:526
Plot_FourPion.C:527
Plot_FourPion.C:528
Plot_FourPion.C:529
Plot_FourPion.C:530
Plot_FourPion.C:531
Plot_FourPion.C:532
Plot_FourPion.C:533
Plot_FourPion.C:534
Plot_FourPion.C:535
Plot_FourPion.C:536
Plot_FourPion.C:537
Plot_FourPion.C:538
Plot_FourPion.C:539
Plot_FourPion.C:540
Plot_FourPion.C:541
Plot_FourPion.C:542
Plot_FourPion.C:543
Plot_FourPion.C:544
Plot_FourPion.C:545
Plot_FourPion.C:546
Plot_FourPion.C:547
Plot_FourPion.C:548
Plot_FourPion.C:549
Plot_FourPion.C:550
Plot_FourPion.C:551
Plot_FourPion.C:552
Plot_FourPion.C:553
Plot_FourPion.C:554
Plot_FourPion.C:555
Plot_FourPion.C:556
Plot_FourPion.C:557
Plot_FourPion.C:558
Plot_FourPion.C:559
Plot_FourPion.C:560
Plot_FourPion.C:561
Plot_FourPion.C:562
Plot_FourPion.C:563
Plot_FourPion.C:564
Plot_FourPion.C:565
Plot_FourPion.C:566
Plot_FourPion.C:567
Plot_FourPion.C:568
Plot_FourPion.C:569
Plot_FourPion.C:570
Plot_FourPion.C:571
Plot_FourPion.C:572
Plot_FourPion.C:573
Plot_FourPion.C:574
Plot_FourPion.C:575
Plot_FourPion.C:576
Plot_FourPion.C:577
Plot_FourPion.C:578
Plot_FourPion.C:579
Plot_FourPion.C:580
Plot_FourPion.C:581
Plot_FourPion.C:582
Plot_FourPion.C:583
Plot_FourPion.C:584
Plot_FourPion.C:585
Plot_FourPion.C:586
Plot_FourPion.C:587
Plot_FourPion.C:588
Plot_FourPion.C:589
Plot_FourPion.C:590
Plot_FourPion.C:591
Plot_FourPion.C:592
Plot_FourPion.C:593
Plot_FourPion.C:594
Plot_FourPion.C:595
Plot_FourPion.C:596
Plot_FourPion.C:597
Plot_FourPion.C:598
Plot_FourPion.C:599
Plot_FourPion.C:600
Plot_FourPion.C:601
Plot_FourPion.C:602
Plot_FourPion.C:603
Plot_FourPion.C:604
Plot_FourPion.C:605
Plot_FourPion.C:606
Plot_FourPion.C:607
Plot_FourPion.C:608
Plot_FourPion.C:609
Plot_FourPion.C:610
Plot_FourPion.C:611
Plot_FourPion.C:612
Plot_FourPion.C:613
Plot_FourPion.C:614
Plot_FourPion.C:615
Plot_FourPion.C:616
Plot_FourPion.C:617
Plot_FourPion.C:618
Plot_FourPion.C:619
Plot_FourPion.C:620
Plot_FourPion.C:621
Plot_FourPion.C:622
Plot_FourPion.C:623
Plot_FourPion.C:624
Plot_FourPion.C:625
Plot_FourPion.C:626
Plot_FourPion.C:627
Plot_FourPion.C:628
Plot_FourPion.C:629
Plot_FourPion.C:630
Plot_FourPion.C:631
Plot_FourPion.C:632
Plot_FourPion.C:633
Plot_FourPion.C:634
Plot_FourPion.C:635
Plot_FourPion.C:636
Plot_FourPion.C:637
Plot_FourPion.C:638
Plot_FourPion.C:639
Plot_FourPion.C:640
Plot_FourPion.C:641
Plot_FourPion.C:642
Plot_FourPion.C:643
Plot_FourPion.C:644
Plot_FourPion.C:645
Plot_FourPion.C:646
Plot_FourPion.C:647
Plot_FourPion.C:648
Plot_FourPion.C:649
Plot_FourPion.C:650
Plot_FourPion.C:651
Plot_FourPion.C:652
Plot_FourPion.C:653
Plot_FourPion.C:654
Plot_FourPion.C:655
Plot_FourPion.C:656
Plot_FourPion.C:657
Plot_FourPion.C:658
Plot_FourPion.C:659
Plot_FourPion.C:660
Plot_FourPion.C:661
Plot_FourPion.C:662
Plot_FourPion.C:663
Plot_FourPion.C:664
Plot_FourPion.C:665
Plot_FourPion.C:666
Plot_FourPion.C:667
Plot_FourPion.C:668
Plot_FourPion.C:669
Plot_FourPion.C:670
Plot_FourPion.C:671
Plot_FourPion.C:672
Plot_FourPion.C:673
Plot_FourPion.C:674
Plot_FourPion.C:675
Plot_FourPion.C:676
Plot_FourPion.C:677
Plot_FourPion.C:678
Plot_FourPion.C:679
Plot_FourPion.C:680
Plot_FourPion.C:681
Plot_FourPion.C:682
Plot_FourPion.C:683
Plot_FourPion.C:684
Plot_FourPion.C:685
Plot_FourPion.C:686
Plot_FourPion.C:687
Plot_FourPion.C:688
Plot_FourPion.C:689
Plot_FourPion.C:690
Plot_FourPion.C:691
Plot_FourPion.C:692
Plot_FourPion.C:693
Plot_FourPion.C:694
Plot_FourPion.C:695
Plot_FourPion.C:696
Plot_FourPion.C:697
Plot_FourPion.C:698
Plot_FourPion.C:699
Plot_FourPion.C:700
Plot_FourPion.C:701
Plot_FourPion.C:702
Plot_FourPion.C:703
Plot_FourPion.C:704
Plot_FourPion.C:705
Plot_FourPion.C:706
Plot_FourPion.C:707
Plot_FourPion.C:708
Plot_FourPion.C:709
Plot_FourPion.C:710
Plot_FourPion.C:711
Plot_FourPion.C:712
Plot_FourPion.C:713
Plot_FourPion.C:714
Plot_FourPion.C:715
Plot_FourPion.C:716
Plot_FourPion.C:717
Plot_FourPion.C:718
Plot_FourPion.C:719
Plot_FourPion.C:720
Plot_FourPion.C:721
Plot_FourPion.C:722
Plot_FourPion.C:723
Plot_FourPion.C:724
Plot_FourPion.C:725
Plot_FourPion.C:726
Plot_FourPion.C:727
Plot_FourPion.C:728
Plot_FourPion.C:729
Plot_FourPion.C:730
Plot_FourPion.C:731
Plot_FourPion.C:732
Plot_FourPion.C:733
Plot_FourPion.C:734
Plot_FourPion.C:735
Plot_FourPion.C:736
Plot_FourPion.C:737
Plot_FourPion.C:738
Plot_FourPion.C:739
Plot_FourPion.C:740
Plot_FourPion.C:741
Plot_FourPion.C:742
Plot_FourPion.C:743
Plot_FourPion.C:744
Plot_FourPion.C:745
Plot_FourPion.C:746
Plot_FourPion.C:747
Plot_FourPion.C:748
Plot_FourPion.C:749
Plot_FourPion.C:750
Plot_FourPion.C:751
Plot_FourPion.C:752
Plot_FourPion.C:753
Plot_FourPion.C:754
Plot_FourPion.C:755
Plot_FourPion.C:756
Plot_FourPion.C:757
Plot_FourPion.C:758
Plot_FourPion.C:759
Plot_FourPion.C:760
Plot_FourPion.C:761
Plot_FourPion.C:762
Plot_FourPion.C:763
Plot_FourPion.C:764
Plot_FourPion.C:765
Plot_FourPion.C:766
Plot_FourPion.C:767
Plot_FourPion.C:768
Plot_FourPion.C:769
Plot_FourPion.C:770
Plot_FourPion.C:771
Plot_FourPion.C:772
Plot_FourPion.C:773
Plot_FourPion.C:774
Plot_FourPion.C:775
Plot_FourPion.C:776
Plot_FourPion.C:777
Plot_FourPion.C:778
Plot_FourPion.C:779
Plot_FourPion.C:780
Plot_FourPion.C:781
Plot_FourPion.C:782
Plot_FourPion.C:783
Plot_FourPion.C:784
Plot_FourPion.C:785
Plot_FourPion.C:786
Plot_FourPion.C:787
Plot_FourPion.C:788
Plot_FourPion.C:789
Plot_FourPion.C:790
Plot_FourPion.C:791
Plot_FourPion.C:792
Plot_FourPion.C:793
Plot_FourPion.C:794
Plot_FourPion.C:795
Plot_FourPion.C:796
Plot_FourPion.C:797
Plot_FourPion.C:798
Plot_FourPion.C:799
Plot_FourPion.C:800
Plot_FourPion.C:801
Plot_FourPion.C:802
Plot_FourPion.C:803
Plot_FourPion.C:804
Plot_FourPion.C:805
Plot_FourPion.C:806
Plot_FourPion.C:807
Plot_FourPion.C:808
Plot_FourPion.C:809
Plot_FourPion.C:810
Plot_FourPion.C:811
Plot_FourPion.C:812
Plot_FourPion.C:813
Plot_FourPion.C:814
Plot_FourPion.C:815
Plot_FourPion.C:816
Plot_FourPion.C:817
Plot_FourPion.C:818
Plot_FourPion.C:819
Plot_FourPion.C:820
Plot_FourPion.C:821
Plot_FourPion.C:822
Plot_FourPion.C:823
Plot_FourPion.C:824
Plot_FourPion.C:825
Plot_FourPion.C:826
Plot_FourPion.C:827
Plot_FourPion.C:828
Plot_FourPion.C:829
Plot_FourPion.C:830
Plot_FourPion.C:831
Plot_FourPion.C:832
Plot_FourPion.C:833
Plot_FourPion.C:834
Plot_FourPion.C:835
Plot_FourPion.C:836
Plot_FourPion.C:837
Plot_FourPion.C:838
Plot_FourPion.C:839
Plot_FourPion.C:840
Plot_FourPion.C:841
Plot_FourPion.C:842
Plot_FourPion.C:843
Plot_FourPion.C:844
Plot_FourPion.C:845
Plot_FourPion.C:846
Plot_FourPion.C:847
Plot_FourPion.C:848
Plot_FourPion.C:849
Plot_FourPion.C:850
Plot_FourPion.C:851
Plot_FourPion.C:852
Plot_FourPion.C:853
Plot_FourPion.C:854
Plot_FourPion.C:855
Plot_FourPion.C:856
Plot_FourPion.C:857
Plot_FourPion.C:858
Plot_FourPion.C:859
Plot_FourPion.C:860
Plot_FourPion.C:861
Plot_FourPion.C:862
Plot_FourPion.C:863
Plot_FourPion.C:864
Plot_FourPion.C:865
Plot_FourPion.C:866
Plot_FourPion.C:867
Plot_FourPion.C:868
Plot_FourPion.C:869
Plot_FourPion.C:870
Plot_FourPion.C:871
Plot_FourPion.C:872
Plot_FourPion.C:873
Plot_FourPion.C:874
Plot_FourPion.C:875
Plot_FourPion.C:876
Plot_FourPion.C:877
Plot_FourPion.C:878
Plot_FourPion.C:879
Plot_FourPion.C:880
Plot_FourPion.C:881
Plot_FourPion.C:882
Plot_FourPion.C:883
Plot_FourPion.C:884
Plot_FourPion.C:885
Plot_FourPion.C:886
Plot_FourPion.C:887
Plot_FourPion.C:888
Plot_FourPion.C:889
Plot_FourPion.C:890
Plot_FourPion.C:891
Plot_FourPion.C:892
Plot_FourPion.C:893
Plot_FourPion.C:894
Plot_FourPion.C:895
Plot_FourPion.C:896
Plot_FourPion.C:897
Plot_FourPion.C:898
Plot_FourPion.C:899
Plot_FourPion.C:900
Plot_FourPion.C:901
Plot_FourPion.C:902
Plot_FourPion.C:903
Plot_FourPion.C:904
Plot_FourPion.C:905
Plot_FourPion.C:906
Plot_FourPion.C:907
Plot_FourPion.C:908
Plot_FourPion.C:909
Plot_FourPion.C:910
Plot_FourPion.C:911
Plot_FourPion.C:912
Plot_FourPion.C:913
Plot_FourPion.C:914
Plot_FourPion.C:915
Plot_FourPion.C:916
Plot_FourPion.C:917
Plot_FourPion.C:918
Plot_FourPion.C:919
Plot_FourPion.C:920
Plot_FourPion.C:921
Plot_FourPion.C:922
Plot_FourPion.C:923
Plot_FourPion.C:924
Plot_FourPion.C:925
Plot_FourPion.C:926
Plot_FourPion.C:927
Plot_FourPion.C:928
Plot_FourPion.C:929
Plot_FourPion.C:930
Plot_FourPion.C:931
Plot_FourPion.C:932
Plot_FourPion.C:933
Plot_FourPion.C:934
Plot_FourPion.C:935
Plot_FourPion.C:936
Plot_FourPion.C:937
Plot_FourPion.C:938
Plot_FourPion.C:939
Plot_FourPion.C:940
Plot_FourPion.C:941
Plot_FourPion.C:942
Plot_FourPion.C:943
Plot_FourPion.C:944
Plot_FourPion.C:945
Plot_FourPion.C:946
Plot_FourPion.C:947
Plot_FourPion.C:948
Plot_FourPion.C:949
Plot_FourPion.C:950
Plot_FourPion.C:951
Plot_FourPion.C:952
Plot_FourPion.C:953
Plot_FourPion.C:954
Plot_FourPion.C:955
Plot_FourPion.C:956
Plot_FourPion.C:957
Plot_FourPion.C:958
Plot_FourPion.C:959
Plot_FourPion.C:960
Plot_FourPion.C:961
Plot_FourPion.C:962
Plot_FourPion.C:963
Plot_FourPion.C:964
Plot_FourPion.C:965
Plot_FourPion.C:966
Plot_FourPion.C:967
Plot_FourPion.C:968
Plot_FourPion.C:969
Plot_FourPion.C:970
Plot_FourPion.C:971
Plot_FourPion.C:972
Plot_FourPion.C:973
Plot_FourPion.C:974
Plot_FourPion.C:975
Plot_FourPion.C:976
Plot_FourPion.C:977
Plot_FourPion.C:978
Plot_FourPion.C:979
Plot_FourPion.C:980
Plot_FourPion.C:981
Plot_FourPion.C:982
Plot_FourPion.C:983
Plot_FourPion.C:984
Plot_FourPion.C:985
Plot_FourPion.C:986
Plot_FourPion.C:987
Plot_FourPion.C:988
Plot_FourPion.C:989
Plot_FourPion.C:990
Plot_FourPion.C:991
Plot_FourPion.C:992
Plot_FourPion.C:993
Plot_FourPion.C:994
Plot_FourPion.C:995
Plot_FourPion.C:996
Plot_FourPion.C:997
Plot_FourPion.C:998
Plot_FourPion.C:999
Plot_FourPion.C:1000
Plot_FourPion.C:1001
Plot_FourPion.C:1002
Plot_FourPion.C:1003
Plot_FourPion.C:1004
Plot_FourPion.C:1005
Plot_FourPion.C:1006
Plot_FourPion.C:1007
Plot_FourPion.C:1008
Plot_FourPion.C:1009
Plot_FourPion.C:1010
Plot_FourPion.C:1011
Plot_FourPion.C:1012
Plot_FourPion.C:1013
Plot_FourPion.C:1014
Plot_FourPion.C:1015
Plot_FourPion.C:1016
Plot_FourPion.C:1017
Plot_FourPion.C:1018
Plot_FourPion.C:1019
Plot_FourPion.C:1020
Plot_FourPion.C:1021
Plot_FourPion.C:1022
Plot_FourPion.C:1023
Plot_FourPion.C:1024
Plot_FourPion.C:1025
Plot_FourPion.C:1026
Plot_FourPion.C:1027
Plot_FourPion.C:1028
Plot_FourPion.C:1029
Plot_FourPion.C:1030
Plot_FourPion.C:1031
Plot_FourPion.C:1032
Plot_FourPion.C:1033
Plot_FourPion.C:1034
Plot_FourPion.C:1035
Plot_FourPion.C:1036
Plot_FourPion.C:1037
Plot_FourPion.C:1038
Plot_FourPion.C:1039
Plot_FourPion.C:1040
Plot_FourPion.C:1041
Plot_FourPion.C:1042
Plot_FourPion.C:1043
Plot_FourPion.C:1044
Plot_FourPion.C:1045
Plot_FourPion.C:1046
Plot_FourPion.C:1047
Plot_FourPion.C:1048
Plot_FourPion.C:1049
Plot_FourPion.C:1050
Plot_FourPion.C:1051
Plot_FourPion.C:1052
Plot_FourPion.C:1053
Plot_FourPion.C:1054
Plot_FourPion.C:1055
Plot_FourPion.C:1056
Plot_FourPion.C:1057
Plot_FourPion.C:1058
Plot_FourPion.C:1059
Plot_FourPion.C:1060
Plot_FourPion.C:1061
Plot_FourPion.C:1062
Plot_FourPion.C:1063
Plot_FourPion.C:1064
Plot_FourPion.C:1065
Plot_FourPion.C:1066
Plot_FourPion.C:1067
Plot_FourPion.C:1068
Plot_FourPion.C:1069
Plot_FourPion.C:1070
Plot_FourPion.C:1071
Plot_FourPion.C:1072
Plot_FourPion.C:1073
Plot_FourPion.C:1074
Plot_FourPion.C:1075
Plot_FourPion.C:1076
Plot_FourPion.C:1077
Plot_FourPion.C:1078
Plot_FourPion.C:1079
Plot_FourPion.C:1080
Plot_FourPion.C:1081
Plot_FourPion.C:1082
Plot_FourPion.C:1083
Plot_FourPion.C:1084
Plot_FourPion.C:1085
Plot_FourPion.C:1086
Plot_FourPion.C:1087
Plot_FourPion.C:1088
Plot_FourPion.C:1089
Plot_FourPion.C:1090
Plot_FourPion.C:1091
Plot_FourPion.C:1092
Plot_FourPion.C:1093
Plot_FourPion.C:1094
Plot_FourPion.C:1095
Plot_FourPion.C:1096
Plot_FourPion.C:1097
Plot_FourPion.C:1098
Plot_FourPion.C:1099
Plot_FourPion.C:1100
Plot_FourPion.C:1101
Plot_FourPion.C:1102
Plot_FourPion.C:1103
Plot_FourPion.C:1104
Plot_FourPion.C:1105
Plot_FourPion.C:1106
Plot_FourPion.C:1107
Plot_FourPion.C:1108
Plot_FourPion.C:1109
Plot_FourPion.C:1110
Plot_FourPion.C:1111
Plot_FourPion.C:1112
Plot_FourPion.C:1113
Plot_FourPion.C:1114
Plot_FourPion.C:1115
Plot_FourPion.C:1116
Plot_FourPion.C:1117
Plot_FourPion.C:1118
Plot_FourPion.C:1119
Plot_FourPion.C:1120
Plot_FourPion.C:1121
Plot_FourPion.C:1122
Plot_FourPion.C:1123
Plot_FourPion.C:1124
Plot_FourPion.C:1125
Plot_FourPion.C:1126
Plot_FourPion.C:1127
Plot_FourPion.C:1128
Plot_FourPion.C:1129
Plot_FourPion.C:1130
Plot_FourPion.C:1131
Plot_FourPion.C:1132
Plot_FourPion.C:1133
Plot_FourPion.C:1134
Plot_FourPion.C:1135
Plot_FourPion.C:1136
Plot_FourPion.C:1137
Plot_FourPion.C:1138
Plot_FourPion.C:1139
Plot_FourPion.C:1140
Plot_FourPion.C:1141
Plot_FourPion.C:1142
Plot_FourPion.C:1143
Plot_FourPion.C:1144
Plot_FourPion.C:1145
Plot_FourPion.C:1146
Plot_FourPion.C:1147
Plot_FourPion.C:1148
Plot_FourPion.C:1149
Plot_FourPion.C:1150
Plot_FourPion.C:1151
Plot_FourPion.C:1152
Plot_FourPion.C:1153
Plot_FourPion.C:1154
Plot_FourPion.C:1155
Plot_FourPion.C:1156
Plot_FourPion.C:1157
Plot_FourPion.C:1158
Plot_FourPion.C:1159
Plot_FourPion.C:1160
Plot_FourPion.C:1161
Plot_FourPion.C:1162
Plot_FourPion.C:1163
Plot_FourPion.C:1164
Plot_FourPion.C:1165
Plot_FourPion.C:1166
Plot_FourPion.C:1167
Plot_FourPion.C:1168
Plot_FourPion.C:1169
Plot_FourPion.C:1170
Plot_FourPion.C:1171
Plot_FourPion.C:1172
Plot_FourPion.C:1173
Plot_FourPion.C:1174
Plot_FourPion.C:1175
Plot_FourPion.C:1176
Plot_FourPion.C:1177
Plot_FourPion.C:1178
Plot_FourPion.C:1179
Plot_FourPion.C:1180
Plot_FourPion.C:1181
Plot_FourPion.C:1182
Plot_FourPion.C:1183
Plot_FourPion.C:1184
Plot_FourPion.C:1185
Plot_FourPion.C:1186
Plot_FourPion.C:1187
Plot_FourPion.C:1188
Plot_FourPion.C:1189
Plot_FourPion.C:1190
Plot_FourPion.C:1191
Plot_FourPion.C:1192
Plot_FourPion.C:1193
Plot_FourPion.C:1194
Plot_FourPion.C:1195
Plot_FourPion.C:1196
Plot_FourPion.C:1197
Plot_FourPion.C:1198
Plot_FourPion.C:1199
Plot_FourPion.C:1200
Plot_FourPion.C:1201
Plot_FourPion.C:1202
Plot_FourPion.C:1203
Plot_FourPion.C:1204
Plot_FourPion.C:1205
Plot_FourPion.C:1206
Plot_FourPion.C:1207
Plot_FourPion.C:1208
Plot_FourPion.C:1209
Plot_FourPion.C:1210
Plot_FourPion.C:1211
Plot_FourPion.C:1212
Plot_FourPion.C:1213
Plot_FourPion.C:1214
Plot_FourPion.C:1215
Plot_FourPion.C:1216
Plot_FourPion.C:1217
Plot_FourPion.C:1218
Plot_FourPion.C:1219
Plot_FourPion.C:1220
Plot_FourPion.C:1221
Plot_FourPion.C:1222
Plot_FourPion.C:1223
Plot_FourPion.C:1224
Plot_FourPion.C:1225
Plot_FourPion.C:1226
Plot_FourPion.C:1227
Plot_FourPion.C:1228
Plot_FourPion.C:1229
Plot_FourPion.C:1230
Plot_FourPion.C:1231
Plot_FourPion.C:1232
Plot_FourPion.C:1233
Plot_FourPion.C:1234
Plot_FourPion.C:1235
Plot_FourPion.C:1236
Plot_FourPion.C:1237
Plot_FourPion.C:1238
Plot_FourPion.C:1239
Plot_FourPion.C:1240
Plot_FourPion.C:1241
Plot_FourPion.C:1242
Plot_FourPion.C:1243
Plot_FourPion.C:1244
Plot_FourPion.C:1245
Plot_FourPion.C:1246
Plot_FourPion.C:1247
Plot_FourPion.C:1248
Plot_FourPion.C:1249
Plot_FourPion.C:1250
Plot_FourPion.C:1251
Plot_FourPion.C:1252
Plot_FourPion.C:1253
Plot_FourPion.C:1254
Plot_FourPion.C:1255
Plot_FourPion.C:1256
Plot_FourPion.C:1257
Plot_FourPion.C:1258
Plot_FourPion.C:1259
Plot_FourPion.C:1260
Plot_FourPion.C:1261
Plot_FourPion.C:1262
Plot_FourPion.C:1263
Plot_FourPion.C:1264
Plot_FourPion.C:1265
Plot_FourPion.C:1266
Plot_FourPion.C:1267
Plot_FourPion.C:1268
Plot_FourPion.C:1269
Plot_FourPion.C:1270
Plot_FourPion.C:1271
Plot_FourPion.C:1272
Plot_FourPion.C:1273
Plot_FourPion.C:1274
Plot_FourPion.C:1275
Plot_FourPion.C:1276
Plot_FourPion.C:1277
Plot_FourPion.C:1278
Plot_FourPion.C:1279
Plot_FourPion.C:1280
Plot_FourPion.C:1281
Plot_FourPion.C:1282
Plot_FourPion.C:1283
Plot_FourPion.C:1284
Plot_FourPion.C:1285
Plot_FourPion.C:1286
Plot_FourPion.C:1287
Plot_FourPion.C:1288
Plot_FourPion.C:1289
Plot_FourPion.C:1290
Plot_FourPion.C:1291
Plot_FourPion.C:1292
Plot_FourPion.C:1293
Plot_FourPion.C:1294
Plot_FourPion.C:1295
Plot_FourPion.C:1296
Plot_FourPion.C:1297
Plot_FourPion.C:1298
Plot_FourPion.C:1299
Plot_FourPion.C:1300
Plot_FourPion.C:1301
Plot_FourPion.C:1302
Plot_FourPion.C:1303
Plot_FourPion.C:1304
Plot_FourPion.C:1305
Plot_FourPion.C:1306
Plot_FourPion.C:1307
Plot_FourPion.C:1308
Plot_FourPion.C:1309
Plot_FourPion.C:1310
Plot_FourPion.C:1311
Plot_FourPion.C:1312
Plot_FourPion.C:1313
Plot_FourPion.C:1314
Plot_FourPion.C:1315
Plot_FourPion.C:1316
Plot_FourPion.C:1317
Plot_FourPion.C:1318
Plot_FourPion.C:1319
Plot_FourPion.C:1320
Plot_FourPion.C:1321
Plot_FourPion.C:1322
Plot_FourPion.C:1323
Plot_FourPion.C:1324
Plot_FourPion.C:1325
Plot_FourPion.C:1326
Plot_FourPion.C:1327
Plot_FourPion.C:1328
Plot_FourPion.C:1329
Plot_FourPion.C:1330
Plot_FourPion.C:1331
Plot_FourPion.C:1332
Plot_FourPion.C:1333
Plot_FourPion.C:1334
Plot_FourPion.C:1335
Plot_FourPion.C:1336
Plot_FourPion.C:1337
Plot_FourPion.C:1338
Plot_FourPion.C:1339
Plot_FourPion.C:1340
Plot_FourPion.C:1341
Plot_FourPion.C:1342
Plot_FourPion.C:1343
Plot_FourPion.C:1344
Plot_FourPion.C:1345
Plot_FourPion.C:1346
Plot_FourPion.C:1347
Plot_FourPion.C:1348
Plot_FourPion.C:1349
Plot_FourPion.C:1350
Plot_FourPion.C:1351
Plot_FourPion.C:1352
Plot_FourPion.C:1353
Plot_FourPion.C:1354
Plot_FourPion.C:1355
Plot_FourPion.C:1356
Plot_FourPion.C:1357
Plot_FourPion.C:1358
Plot_FourPion.C:1359
Plot_FourPion.C:1360
Plot_FourPion.C:1361
Plot_FourPion.C:1362
Plot_FourPion.C:1363
Plot_FourPion.C:1364
Plot_FourPion.C:1365
Plot_FourPion.C:1366
Plot_FourPion.C:1367
Plot_FourPion.C:1368
Plot_FourPion.C:1369
Plot_FourPion.C:1370
Plot_FourPion.C:1371
Plot_FourPion.C:1372
Plot_FourPion.C:1373
Plot_FourPion.C:1374
Plot_FourPion.C:1375
Plot_FourPion.C:1376
Plot_FourPion.C:1377
Plot_FourPion.C:1378
Plot_FourPion.C:1379
Plot_FourPion.C:1380
Plot_FourPion.C:1381
Plot_FourPion.C:1382
Plot_FourPion.C:1383
Plot_FourPion.C:1384
Plot_FourPion.C:1385
Plot_FourPion.C:1386
Plot_FourPion.C:1387
Plot_FourPion.C:1388
Plot_FourPion.C:1389
Plot_FourPion.C:1390
Plot_FourPion.C:1391
Plot_FourPion.C:1392
Plot_FourPion.C:1393
Plot_FourPion.C:1394
Plot_FourPion.C:1395
Plot_FourPion.C:1396
Plot_FourPion.C:1397
Plot_FourPion.C:1398
Plot_FourPion.C:1399
Plot_FourPion.C:1400
Plot_FourPion.C:1401
Plot_FourPion.C:1402
Plot_FourPion.C:1403
Plot_FourPion.C:1404
Plot_FourPion.C:1405
Plot_FourPion.C:1406
Plot_FourPion.C:1407
Plot_FourPion.C:1408
Plot_FourPion.C:1409
Plot_FourPion.C:1410
Plot_FourPion.C:1411
Plot_FourPion.C:1412
Plot_FourPion.C:1413
Plot_FourPion.C:1414
Plot_FourPion.C:1415
Plot_FourPion.C:1416
Plot_FourPion.C:1417
Plot_FourPion.C:1418
Plot_FourPion.C:1419
Plot_FourPion.C:1420
Plot_FourPion.C:1421
Plot_FourPion.C:1422
Plot_FourPion.C:1423
Plot_FourPion.C:1424
Plot_FourPion.C:1425
Plot_FourPion.C:1426
Plot_FourPion.C:1427
Plot_FourPion.C:1428
Plot_FourPion.C:1429
Plot_FourPion.C:1430
Plot_FourPion.C:1431
Plot_FourPion.C:1432
Plot_FourPion.C:1433
Plot_FourPion.C:1434
Plot_FourPion.C:1435
Plot_FourPion.C:1436
Plot_FourPion.C:1437
Plot_FourPion.C:1438
Plot_FourPion.C:1439
Plot_FourPion.C:1440
Plot_FourPion.C:1441
Plot_FourPion.C:1442
Plot_FourPion.C:1443
Plot_FourPion.C:1444
Plot_FourPion.C:1445
Plot_FourPion.C:1446
Plot_FourPion.C:1447
Plot_FourPion.C:1448
Plot_FourPion.C:1449
Plot_FourPion.C:1450
Plot_FourPion.C:1451
Plot_FourPion.C:1452
Plot_FourPion.C:1453
Plot_FourPion.C:1454
Plot_FourPion.C:1455
Plot_FourPion.C:1456
Plot_FourPion.C:1457
Plot_FourPion.C:1458
Plot_FourPion.C:1459
Plot_FourPion.C:1460
Plot_FourPion.C:1461
Plot_FourPion.C:1462
Plot_FourPion.C:1463
Plot_FourPion.C:1464
Plot_FourPion.C:1465
Plot_FourPion.C:1466
Plot_FourPion.C:1467
Plot_FourPion.C:1468
Plot_FourPion.C:1469
Plot_FourPion.C:1470
Plot_FourPion.C:1471
Plot_FourPion.C:1472
Plot_FourPion.C:1473
Plot_FourPion.C:1474
Plot_FourPion.C:1475
Plot_FourPion.C:1476
Plot_FourPion.C:1477
Plot_FourPion.C:1478
Plot_FourPion.C:1479
Plot_FourPion.C:1480
Plot_FourPion.C:1481
Plot_FourPion.C:1482
Plot_FourPion.C:1483
Plot_FourPion.C:1484
Plot_FourPion.C:1485
Plot_FourPion.C:1486
Plot_FourPion.C:1487
Plot_FourPion.C:1488
Plot_FourPion.C:1489
Plot_FourPion.C:1490
Plot_FourPion.C:1491
Plot_FourPion.C:1492
Plot_FourPion.C:1493
Plot_FourPion.C:1494
Plot_FourPion.C:1495
Plot_FourPion.C:1496
Plot_FourPion.C:1497
Plot_FourPion.C:1498
Plot_FourPion.C:1499
Plot_FourPion.C:1500
Plot_FourPion.C:1501
Plot_FourPion.C:1502
Plot_FourPion.C:1503
Plot_FourPion.C:1504
Plot_FourPion.C:1505
Plot_FourPion.C:1506
Plot_FourPion.C:1507
Plot_FourPion.C:1508
Plot_FourPion.C:1509
Plot_FourPion.C:1510
Plot_FourPion.C:1511
Plot_FourPion.C:1512
Plot_FourPion.C:1513
Plot_FourPion.C:1514
Plot_FourPion.C:1515
Plot_FourPion.C:1516
Plot_FourPion.C:1517
Plot_FourPion.C:1518
Plot_FourPion.C:1519
Plot_FourPion.C:1520
Plot_FourPion.C:1521
Plot_FourPion.C:1522
Plot_FourPion.C:1523
Plot_FourPion.C:1524
Plot_FourPion.C:1525
Plot_FourPion.C:1526
Plot_FourPion.C:1527
Plot_FourPion.C:1528
Plot_FourPion.C:1529
Plot_FourPion.C:1530
Plot_FourPion.C:1531
Plot_FourPion.C:1532
Plot_FourPion.C:1533
Plot_FourPion.C:1534
Plot_FourPion.C:1535
Plot_FourPion.C:1536
Plot_FourPion.C:1537
Plot_FourPion.C:1538
Plot_FourPion.C:1539
Plot_FourPion.C:1540
Plot_FourPion.C:1541
Plot_FourPion.C:1542
Plot_FourPion.C:1543
Plot_FourPion.C:1544
Plot_FourPion.C:1545
Plot_FourPion.C:1546
Plot_FourPion.C:1547
Plot_FourPion.C:1548
Plot_FourPion.C:1549
Plot_FourPion.C:1550
Plot_FourPion.C:1551
Plot_FourPion.C:1552
Plot_FourPion.C:1553
Plot_FourPion.C:1554
Plot_FourPion.C:1555
Plot_FourPion.C:1556
Plot_FourPion.C:1557
Plot_FourPion.C:1558
Plot_FourPion.C:1559
Plot_FourPion.C:1560
Plot_FourPion.C:1561
Plot_FourPion.C:1562
Plot_FourPion.C:1563
Plot_FourPion.C:1564
Plot_FourPion.C:1565
Plot_FourPion.C:1566
Plot_FourPion.C:1567
Plot_FourPion.C:1568
Plot_FourPion.C:1569
Plot_FourPion.C:1570
Plot_FourPion.C:1571
Plot_FourPion.C:1572
Plot_FourPion.C:1573
Plot_FourPion.C:1574
Plot_FourPion.C:1575
Plot_FourPion.C:1576
Plot_FourPion.C:1577
Plot_FourPion.C:1578
Plot_FourPion.C:1579
Plot_FourPion.C:1580
Plot_FourPion.C:1581
Plot_FourPion.C:1582
Plot_FourPion.C:1583
Plot_FourPion.C:1584
Plot_FourPion.C:1585
Plot_FourPion.C:1586
Plot_FourPion.C:1587
Plot_FourPion.C:1588
Plot_FourPion.C:1589
Plot_FourPion.C:1590
Plot_FourPion.C:1591
Plot_FourPion.C:1592
Plot_FourPion.C:1593
Plot_FourPion.C:1594
Plot_FourPion.C:1595
Plot_FourPion.C:1596
Plot_FourPion.C:1597
Plot_FourPion.C:1598
Plot_FourPion.C:1599
Plot_FourPion.C:1600
Plot_FourPion.C:1601
Plot_FourPion.C:1602
Plot_FourPion.C:1603
Plot_FourPion.C:1604
Plot_FourPion.C:1605
Plot_FourPion.C:1606
Plot_FourPion.C:1607
Plot_FourPion.C:1608
Plot_FourPion.C:1609
Plot_FourPion.C:1610
Plot_FourPion.C:1611
Plot_FourPion.C:1612
Plot_FourPion.C:1613
Plot_FourPion.C:1614
Plot_FourPion.C:1615
Plot_FourPion.C:1616
Plot_FourPion.C:1617
Plot_FourPion.C:1618
Plot_FourPion.C:1619
Plot_FourPion.C:1620
Plot_FourPion.C:1621
Plot_FourPion.C:1622
Plot_FourPion.C:1623
Plot_FourPion.C:1624
Plot_FourPion.C:1625
Plot_FourPion.C:1626
Plot_FourPion.C:1627
Plot_FourPion.C:1628
Plot_FourPion.C:1629
Plot_FourPion.C:1630
Plot_FourPion.C:1631
Plot_FourPion.C:1632
Plot_FourPion.C:1633
Plot_FourPion.C:1634
Plot_FourPion.C:1635
Plot_FourPion.C:1636
Plot_FourPion.C:1637
Plot_FourPion.C:1638
Plot_FourPion.C:1639
Plot_FourPion.C:1640
Plot_FourPion.C:1641
Plot_FourPion.C:1642
Plot_FourPion.C:1643
Plot_FourPion.C:1644
Plot_FourPion.C:1645
Plot_FourPion.C:1646
Plot_FourPion.C:1647
Plot_FourPion.C:1648
Plot_FourPion.C:1649
Plot_FourPion.C:1650
Plot_FourPion.C:1651
Plot_FourPion.C:1652
Plot_FourPion.C:1653
Plot_FourPion.C:1654
Plot_FourPion.C:1655
Plot_FourPion.C:1656
Plot_FourPion.C:1657
Plot_FourPion.C:1658
Plot_FourPion.C:1659
Plot_FourPion.C:1660
Plot_FourPion.C:1661
Plot_FourPion.C:1662
Plot_FourPion.C:1663
Plot_FourPion.C:1664
Plot_FourPion.C:1665
Plot_FourPion.C:1666
Plot_FourPion.C:1667
Plot_FourPion.C:1668
Plot_FourPion.C:1669
Plot_FourPion.C:1670
Plot_FourPion.C:1671
Plot_FourPion.C:1672
Plot_FourPion.C:1673
Plot_FourPion.C:1674
Plot_FourPion.C:1675
Plot_FourPion.C:1676
Plot_FourPion.C:1677
Plot_FourPion.C:1678
Plot_FourPion.C:1679
Plot_FourPion.C:1680
Plot_FourPion.C:1681
Plot_FourPion.C:1682
Plot_FourPion.C:1683
Plot_FourPion.C:1684
Plot_FourPion.C:1685
Plot_FourPion.C:1686
Plot_FourPion.C:1687
Plot_FourPion.C:1688
Plot_FourPion.C:1689
Plot_FourPion.C:1690
Plot_FourPion.C:1691
Plot_FourPion.C:1692
Plot_FourPion.C:1693
Plot_FourPion.C:1694
Plot_FourPion.C:1695
Plot_FourPion.C:1696
Plot_FourPion.C:1697
Plot_FourPion.C:1698
Plot_FourPion.C:1699
Plot_FourPion.C:1700
Plot_FourPion.C:1701
Plot_FourPion.C:1702
Plot_FourPion.C:1703
Plot_FourPion.C:1704
Plot_FourPion.C:1705
Plot_FourPion.C:1706
Plot_FourPion.C:1707
Plot_FourPion.C:1708
Plot_FourPion.C:1709
Plot_FourPion.C:1710
Plot_FourPion.C:1711
Plot_FourPion.C:1712
Plot_FourPion.C:1713
Plot_FourPion.C:1714
Plot_FourPion.C:1715
Plot_FourPion.C:1716
Plot_FourPion.C:1717
Plot_FourPion.C:1718
Plot_FourPion.C:1719
Plot_FourPion.C:1720
Plot_FourPion.C:1721
Plot_FourPion.C:1722
Plot_FourPion.C:1723
Plot_FourPion.C:1724
Plot_FourPion.C:1725
Plot_FourPion.C:1726
Plot_FourPion.C:1727
Plot_FourPion.C:1728
Plot_FourPion.C:1729
Plot_FourPion.C:1730
Plot_FourPion.C:1731
Plot_FourPion.C:1732
Plot_FourPion.C:1733
Plot_FourPion.C:1734
Plot_FourPion.C:1735
Plot_FourPion.C:1736
Plot_FourPion.C:1737
Plot_FourPion.C:1738
Plot_FourPion.C:1739
Plot_FourPion.C:1740
Plot_FourPion.C:1741
Plot_FourPion.C:1742
Plot_FourPion.C:1743
Plot_FourPion.C:1744
Plot_FourPion.C:1745
Plot_FourPion.C:1746
Plot_FourPion.C:1747
Plot_FourPion.C:1748
Plot_FourPion.C:1749
Plot_FourPion.C:1750
Plot_FourPion.C:1751
Plot_FourPion.C:1752
Plot_FourPion.C:1753
Plot_FourPion.C:1754
Plot_FourPion.C:1755
Plot_FourPion.C:1756
Plot_FourPion.C:1757
Plot_FourPion.C:1758
Plot_FourPion.C:1759
Plot_FourPion.C:1760
Plot_FourPion.C:1761
Plot_FourPion.C:1762
Plot_FourPion.C:1763
Plot_FourPion.C:1764
Plot_FourPion.C:1765
Plot_FourPion.C:1766
Plot_FourPion.C:1767
Plot_FourPion.C:1768
Plot_FourPion.C:1769
Plot_FourPion.C:1770
Plot_FourPion.C:1771
Plot_FourPion.C:1772
Plot_FourPion.C:1773
Plot_FourPion.C:1774
Plot_FourPion.C:1775
Plot_FourPion.C:1776
Plot_FourPion.C:1777
Plot_FourPion.C:1778
Plot_FourPion.C:1779
Plot_FourPion.C:1780
Plot_FourPion.C:1781
Plot_FourPion.C:1782
Plot_FourPion.C:1783
Plot_FourPion.C:1784
Plot_FourPion.C:1785
Plot_FourPion.C:1786
Plot_FourPion.C:1787
Plot_FourPion.C:1788
Plot_FourPion.C:1789
Plot_FourPion.C:1790
Plot_FourPion.C:1791
Plot_FourPion.C:1792
Plot_FourPion.C:1793
Plot_FourPion.C:1794
Plot_FourPion.C:1795
Plot_FourPion.C:1796
Plot_FourPion.C:1797
Plot_FourPion.C:1798
Plot_FourPion.C:1799
Plot_FourPion.C:1800
Plot_FourPion.C:1801
Plot_FourPion.C:1802
Plot_FourPion.C:1803
Plot_FourPion.C:1804
Plot_FourPion.C:1805
Plot_FourPion.C:1806
Plot_FourPion.C:1807
Plot_FourPion.C:1808
Plot_FourPion.C:1809
Plot_FourPion.C:1810
Plot_FourPion.C:1811
Plot_FourPion.C:1812
Plot_FourPion.C:1813
Plot_FourPion.C:1814
Plot_FourPion.C:1815
Plot_FourPion.C:1816
Plot_FourPion.C:1817
Plot_FourPion.C:1818
Plot_FourPion.C:1819
Plot_FourPion.C:1820
Plot_FourPion.C:1821
Plot_FourPion.C:1822
Plot_FourPion.C:1823
Plot_FourPion.C:1824
Plot_FourPion.C:1825
Plot_FourPion.C:1826
Plot_FourPion.C:1827
Plot_FourPion.C:1828
Plot_FourPion.C:1829
Plot_FourPion.C:1830
Plot_FourPion.C:1831
Plot_FourPion.C:1832
Plot_FourPion.C:1833
Plot_FourPion.C:1834
Plot_FourPion.C:1835
Plot_FourPion.C:1836
Plot_FourPion.C:1837
Plot_FourPion.C:1838
Plot_FourPion.C:1839
Plot_FourPion.C:1840
Plot_FourPion.C:1841
Plot_FourPion.C:1842
Plot_FourPion.C:1843
Plot_FourPion.C:1844
Plot_FourPion.C:1845
Plot_FourPion.C:1846
Plot_FourPion.C:1847
Plot_FourPion.C:1848
Plot_FourPion.C:1849
Plot_FourPion.C:1850
Plot_FourPion.C:1851
Plot_FourPion.C:1852
Plot_FourPion.C:1853
Plot_FourPion.C:1854
Plot_FourPion.C:1855
Plot_FourPion.C:1856
Plot_FourPion.C:1857
Plot_FourPion.C:1858
Plot_FourPion.C:1859
Plot_FourPion.C:1860
Plot_FourPion.C:1861
Plot_FourPion.C:1862
Plot_FourPion.C:1863
Plot_FourPion.C:1864
Plot_FourPion.C:1865
Plot_FourPion.C:1866
Plot_FourPion.C:1867
Plot_FourPion.C:1868
Plot_FourPion.C:1869
Plot_FourPion.C:1870
Plot_FourPion.C:1871
Plot_FourPion.C:1872
Plot_FourPion.C:1873
Plot_FourPion.C:1874
Plot_FourPion.C:1875
Plot_FourPion.C:1876
Plot_FourPion.C:1877
Plot_FourPion.C:1878
Plot_FourPion.C:1879
Plot_FourPion.C:1880
Plot_FourPion.C:1881
Plot_FourPion.C:1882
Plot_FourPion.C:1883
Plot_FourPion.C:1884
Plot_FourPion.C:1885
Plot_FourPion.C:1886
Plot_FourPion.C:1887
Plot_FourPion.C:1888
Plot_FourPion.C:1889
Plot_FourPion.C:1890
Plot_FourPion.C:1891
Plot_FourPion.C:1892
Plot_FourPion.C:1893
Plot_FourPion.C:1894
Plot_FourPion.C:1895
Plot_FourPion.C:1896
Plot_FourPion.C:1897
Plot_FourPion.C:1898
Plot_FourPion.C:1899
Plot_FourPion.C:1900
Plot_FourPion.C:1901
Plot_FourPion.C:1902
Plot_FourPion.C:1903
Plot_FourPion.C:1904
Plot_FourPion.C:1905
Plot_FourPion.C:1906
Plot_FourPion.C:1907
Plot_FourPion.C:1908
Plot_FourPion.C:1909
Plot_FourPion.C:1910
Plot_FourPion.C:1911
Plot_FourPion.C:1912
Plot_FourPion.C:1913
Plot_FourPion.C:1914
Plot_FourPion.C:1915
Plot_FourPion.C:1916
Plot_FourPion.C:1917
Plot_FourPion.C:1918
Plot_FourPion.C:1919
Plot_FourPion.C:1920
Plot_FourPion.C:1921
Plot_FourPion.C:1922
Plot_FourPion.C:1923
Plot_FourPion.C:1924
Plot_FourPion.C:1925
Plot_FourPion.C:1926
Plot_FourPion.C:1927
Plot_FourPion.C:1928
Plot_FourPion.C:1929
Plot_FourPion.C:1930
Plot_FourPion.C:1931
Plot_FourPion.C:1932
Plot_FourPion.C:1933
Plot_FourPion.C:1934
Plot_FourPion.C:1935
Plot_FourPion.C:1936
Plot_FourPion.C:1937
Plot_FourPion.C:1938
Plot_FourPion.C:1939
Plot_FourPion.C:1940