systematics/cutResolution/VertexResolution.C

0001 void VertexResolution()
0002 {
0003   TFile* KS0DataFile = new TFile("/sphenix/user/rosstom/TrackingAnalysis/Physics_Val_TF/data/analysisMacros/KShort_fullDataset_finalCuts_0p2pTCut_rapidity1p0Cut_BCOcut_.root");
0004   TTree* KS0DataTree = (TTree*)KS0DataFile->Get("DecayTree");
0005
0006   TFile* XiDataFile = new TFile("/sphenix/user/rosstom/TrackingAnalysis/Physics_Val_TF/data/analysisMacros/Ximinus_fullDataset_finalCuts_0p2pTCut_rapidity1p0Cut_BCOcut_.root");
0007   TTree* XiDataTree = (TTree*)XiDataFile->Get("DecayTree");
0008
0009   TH2F *h_PV_x_y_KS0 = new TH2F("h_PV_x_y_KS0",";Primary Vertex x [cm];Primary Vertex y [cm]",50,-0.2,0.2,50,-0.2,0.2);
0010   TH2F *h_PV_x_y_Xi = new TH2F("h_PV_x_y_Xi",";Primary Vertex x [cm];Primary Vertex y [cm]",50,-0.2,0.2,50,-0.2,0.2);
0011   TH2F *h_PV_z_r_KS0 = new TH2F("h_PV_z_r_KS0",";Primary Vertex z [cm];Primary Vertex r [cm]",50,-10,10,50,0,0.2);
0012   TH2F *h_PV_z_r_Xi = new TH2F("h_PV_z_r_Xi",";Primary Vertex z [cm];Primary Vertex r [cm]",50,-10,10,50,0,0.2);
0013
0014   TH2F *h_PVcov_x_y_KS0 = new TH2F("h_PVcov_x_y_KS0",";Primary Vertex #sigma_{x};Primary Vertex #sigma_{y}",50,0,0.01,50,0,0.01);
0015   TH2F *h_PVcov_x_y_Xi = new TH2F("h_PVcov_x_y_Xi",";Primary Vertex #sigma_{x};Primary Vertex #sigma_{y}",50,0,0.01,50,0,0.01);
0016   TH2F *h_PVcov_x_z_KS0 = new TH2F("h_PVcov_x_z_KS0",";Primary Vertex #sigma_{x};Primary Vertex #sigma_{z}",50,0,0.01,50,0,0.01);
0017   TH2F *h_PVcov_x_z_Xi = new TH2F("h_PVcov_x_z_Xi",";Primary Vertex #sigma_{x};Primary Vertex #sigma_{z}",50,0,0.01,50,0,0.01);
0018
0019   TH1F *h_PVvolume_r_KS0 = new TH1F("h_PVvolume_r_KS0",";Primary Vertex Volume Radius [cm];Entries",50,0,0.01);
0020   TH1F *h_PVvolume_r_Xi = new TH1F("h_PVvolume_r_Xi",";Primary Vertex Volume Radius [cm];Entries",50,0,0.01);
0021   TH1F *h_PVcov_r_KS0 = new TH1F("h_PVcov_r_KS0",";Primary Vertex #sigma_{R};Entries",50,0,0.01);
0022   TH1F *h_PVcov_r_Xi = new TH1F("h_PVcov_r_Xi",";Primary Vertex #sigma_{R};Entries",50,0,0.01);
0023   TH1F *h_PVcov_rxy_KS0 = new TH1F("h_PVcov_rxy_KS0",";Primary Vertex #sigma_{Rxy};Entries",50,0,0.01);
0024   TH1F *h_PVcov_rxy_Xi = new TH1F("h_PVcov_rxy_Xi",";Primary Vertex #sigma_{Rxy};Entries",50,0,0.01);
0025
0026   TH2F *h_PVcov_r_nTracks_KS0 = new TH2F("h_PVcov_r_nTracks_KS0",";Number of Vertex Tracks;Primary Vertex #sigma_{R}",10,0,10,50,0,0.01);
0027
0028   TH2F *h_SV_x_y_KS0 = new TH2F("h_SV_x_y_KS0",";Secondary Vertex x [cm];Secondary Vertex y [cm]",50,-5,5,50,-5,5);
0029   TH2F *h_SV_x_y_Xi = new TH2F("h_SV_x_y_Xi",";Secondary Vertex x [cm];Secondary Vertex y [cm]",50,-5,5,50,-5,5);
0030   TH2F *h_SV_z_r_KS0 = new TH2F("h_SV_z_r_KS0",";Secondary Vertex z [cm];Secondary Vertex r [cm]",50,-10,10,50,0,5);
0031   TH2F *h_SV_z_r_Xi = new TH2F("h_SV_z_r_Xi",";Secondary Vertex z [cm];Secondary Vertex r [cm]",50,-10,10,50,0,5);
0032
0033   TH2F *h_SVcov_x_y_KS0 = new TH2F("h_SVcov_x_y_KS0",";Secondary Vertex #sigma_{x};Secondary Vertex #sigma_{y}",50,0,0.004,50,0,0.004);
0034   TH2F *h_SVcov_x_y_Xi = new TH2F("h_SVcov_x_y_Xi",";Secondary Vertex #sigma_{x};Secondary Vertex #sigma_{y}",50,0,0.004,50,0,0.004);
0035   TH2F *h_SVcov_x_z_KS0 = new TH2F("h_SVcov_x_z_KS0",";Secondary Vertex #sigma_{x};Secondary Vertex #sigma_{z}",50,0,0.004,50,0,0.004);
0036   TH2F *h_SVcov_x_z_Xi = new TH2F("h_SVcov_x_z_Xi",";Secondary Vertex #sigma_{x};Secondary Vertex #sigma_{z}",50,0,0.004,50,0,0.004);
0037
0038   TH1F *h_SVvolume_r_KS0 = new TH1F("h_SVvolume_r_KS0",";Secondary Vertex Volume Radius [cm];Entries",50,0,0.01);
0039   TH1F *h_SVvolume_r_Xi = new TH1F("h_SVvolume_r_Xi",";Secondary Vertex Volume Radius [cm];Entries",50,0,0.01);
0040   TH1F *h_SVcov_r_KS0 = new TH1F("h_SVcov_r_KS0",";Secondary Vertex #sigma_{R};Entries",50,0,0.01);
0041   TH1F *h_SVcov_r_Xi = new TH1F("h_SVcov_r_Xi",";Secondary Vertex #sigma_{R};Entries",50,0,0.01);
0042   TH1F *h_SVcov_rxy_KS0 = new TH1F("h_SVcov_rxy_KS0",";Secondary Vertex #sigma_{Rxy};Entries",50,0,0.01);
0043   TH1F *h_SVcov_rxy_Xi = new TH1F("h_SVcov_rxy_Xi",";Secondary Vertex #sigma_{Rxy};Entries",50,0,0.01);
0044
0045   TH2F *h_Lambda_x_y_Xi = new TH2F("h_Lambda_x_y_Xi",";#Lambda^{0} Vertex x [cm];#Lambda^{0} Vertex y [cm]",50,-5,5,50,-5,5);
0046   TH2F *h_Lambda_z_r_Xi = new TH2F("h_Lambda_z_r_Xi",";#Lambda^{0} Vertex z [cm];#Lambda^{0} Vertex r [cm]",50,-10,10,50,0,5);
0047
0048   TH2F *h_Lambdacov_x_y_Xi = new TH2F("h_Lambdacov_x_y_Xi",";#Lambda^{0} Vertex #sigma_{x};#Lambda^{0} Vertex #sigma_{y}",50,0,0.004,50,0,0.004);
0049   TH2F *h_Lambdacov_x_z_Xi = new TH2F("h_Lambdacov_x_z_Xi",";#Lambda^{0} Vertex #sigma_{x};#Lambda^{0} Vertex #sigma_{z}",50,0,0.004,50,0,0.004);
0050
0051   TH1F *h_Lambdavolume_r_Xi = new TH1F("h_Lambdavolume_r_Xi",";#Lambda^{0} Vertex Volume Radius [cm];Entries",50,0,0.01);
0052   TH1F *h_Lambdacov_r_Xi = new TH1F("h_Lambdacov_r_Xi",";#Lambda^{0} Vertex #sigma_{R};Entries",50,0,0.01);
0053   TH1F *h_Lambdacov_rxy_Xi = new TH1F("h_Lambdacov_rxy_Xi",";#Lambda^{0} Vertex #sigma_{Rxy};Entries",50,0,0.01);
0054
0055   TH1F *h_pion_track1_r_KS0 = new TH1F("h_pion_track1_r_KS0",";#pi^{#pm} Track 1 #sigma_{R};Entries",50,0,0.01);
0056   TH1F *h_pion_track2_r_KS0 = new TH1F("h_pion_track2_r_KS0",";#pi^{#pm} Track 1 #sigma_{R};Entries",50,0,0.01);
0057   TH1F *h_pion_track1_rxy_KS0 = new TH1F("h_pion_track1_rxy_KS0",";#pi^{#pm} Track 1 #sigma_{Rxy};Entries",50,0,0.01);
0058   TH1F *h_pion_track2_rxy_KS0 = new TH1F("h_pion_track2_rxy_KS0",";#pi^{#pm} Track 1 #sigma_{Rxy};Entries",50,0,0.01);
0059
0060   TH1F *h_pion_track_1_r_Xi = new TH1F("h_pion_track_1_r_Xi",";#pi^{#pm} Track 1 #sigma_{R};Entries",50,0,0.01);
0061   TH1F *h_proton_track_2_r_Xi = new TH1F("h_proton_track_2_r_Xi",";p^{#pm} Track 2 #sigma_{R};Entries",50,0,0.01);
0062   TH1F *h_pion_track_3_r_Xi = new TH1F("h_pion_track_3_r_Xi",";#pi^{#pm} Track 3 #sigma_{R};Entries",50,0,0.01);
0063   TH1F *h_pion_track_1_rxy_Xi = new TH1F("h_pion_track_1_rxy_Xi",";#pi^{#pm} Track 1 #sigma_{Rxy};Entries",50,0,0.01);
0064   TH1F *h_proton_track_2_rxy_Xi = new TH1F("h_proton_track_2_rxy_Xi",";p^{#pm} Track 2 #sigma_{Rxy};Entries",50,0,0.01);
0065   TH1F *h_pion_track_3_rxy_Xi = new TH1F("h_pion_track_3_rxy_Xi",";#pi^{#pm} Track 3 #sigma_{Rxy};Entries",50,0,0.01);
0066
0067   std::string Xi_charge_cut = "Lambda0_track_1_charge == track_3_charge";
0068
0069   /*
0070   KS0DataTree->Draw("primary_vertex_y:primary_vertex_x>>h_PV_x_y_KS0");
0071   XiDataTree->Draw("primary_vertex_y:primary_vertex_x>>h_PV_x_y_Xi",Xi_charge_cut.c_str());
0072   KS0DataTree->Draw("sqrt(pow(primary_vertex_x,2) + pow(primary_vertex_y,2)):primary_vertex_z>>h_PV_z_r_KS0");
0073   XiDataTree->Draw("sqrt(pow(primary_vertex_x,2) + pow(primary_vertex_y,2)):primary_vertex_z>>h_PV_z_r_Xi",Xi_charge_cut.c_str());
0074
0075   KS0DataTree->Draw("sqrt(primary_vertex_Covariance[2]):sqrt(primary_vertex_Covariance[0])>>h_PVcov_x_y_KS0");
0076   XiDataTree->Draw("sqrt(primary_vertex_Covariance[2]):sqrt(primary_vertex_Covariance[0])>>h_PVcov_x_y_Xi",Xi_charge_cut.c_str());
0077   KS0DataTree->Draw("sqrt(primary_vertex_Covariance[5]):sqrt(primary_vertex_Covariance[0])>>h_PVcov_x_z_KS0");
0078   XiDataTree->Draw("sqrt(primary_vertex_Covariance[5]):sqrt(primary_vertex_Covariance[0])>>h_PVcov_x_z_Xi",Xi_charge_cut.c_str());
0079
0080   KS0DataTree->Draw("pow((3*primary_vertex_volume)/(4*pi),(1./3.))>>h_PVvolume_r_KS0");
0081   XiDataTree->Draw("pow((3*primary_vertex_volume)/(4*pi),(1./3.))>>h_PVvolume_r_Xi",Xi_charge_cut.c_str());
0082   */
0083   KS0DataTree->Draw("pow(primary_vertex_Covariance[0] + primary_vertex_Covariance[2] + primary_vertex_Covariance[5],1/2)>>h_PVcov_r_KS0");
0084   /*
0085   XiDataTree->Draw("pow(primary_vertex_Covariance[0] + primary_vertex_Covariance[2] + primary_vertex_Covariance[5],1/2)>>h_PVcov_r_Xi",Xi_charge_cut.c_str());
0086   KS0DataTree->Draw("pow(primary_vertex_Covariance[0] + primary_vertex_Covariance[2],1/2)>>h_PVcov_rxy_KS0");
0087   */
0088   XiDataTree->Draw("pow(primary_vertex_Covariance[0] + primary_vertex_Covariance[2],1/2)>>h_PVcov_rxy_Xi",Xi_charge_cut.c_str());
0089   /*
0090   KS0DataTree->Draw("pow(primary_vertex_Covariance[0] + primary_vertex_Covariance[2] + primary_vertex_Covariance[5],1/2):nTracksOfVertex>>h_PVcov_r_nTracks_KS0");
0091
0092   KS0DataTree->Draw("K_S0_y:K_S0_x>>h_SV_x_y_KS0");
0093   XiDataTree->Draw("Ximinus_y:Ximinus_x>>h_SV_x_y_Xi",Xi_charge_cut.c_str());
0094   KS0DataTree->Draw("sqrt(pow(K_S0_x,2) + pow(K_S0_y,2)):K_S0_z>>h_SV_z_r_KS0");
0095   XiDataTree->Draw("sqrt(pow(Ximinus_x,2) + pow(Ximinus_y,2)):Ximinus_z>>h_SV_z_r_Xi",Xi_charge_cut.c_str());
0096
0097   KS0DataTree->Draw("sqrt(K_S0_Covariance[2]):sqrt(K_S0_Covariance[0])>>h_SVcov_x_y_KS0");
0098   XiDataTree->Draw("sqrt(Ximinus_Covariance[2]):sqrt(Ximinus_Covariance[0])>>h_SVcov_x_y_Xi",Xi_charge_cut.c_str());
0099   KS0DataTree->Draw("sqrt(K_S0_Covariance[5]):sqrt(K_S0_Covariance[0])>>h_SVcov_x_z_KS0");
0100   XiDataTree->Draw("sqrt(Ximinus_Covariance[5]):sqrt(Ximinus_Covariance[0])>>h_SVcov_x_z_Xi",Xi_charge_cut.c_str());
0101
0102   KS0DataTree->Draw("pow((3*K_S0_vertex_volume)/(4*pi),(1./3.))>>h_SVvolume_r_KS0");
0103   XiDataTree->Draw("pow((3*Ximinus_vertex_volume)/(4*pi),(1./3.))>>h_SVvolume_r_Xi",Xi_charge_cut.c_str());
0104   KS0DataTree->Draw("pow(K_S0_Covariance[0] + K_S0_Covariance[2] + K_S0_Covariance[5],1/2)>>h_SVcov_r_KS0");
0105   XiDataTree->Draw("pow(Ximinus_Covariance[0] + Ximinus_Covariance[2] + Ximinus_Covariance[5],1/2)>>h_SVcov_r_Xi",Xi_charge_cut.c_str());
0106   KS0DataTree->Draw("pow(K_S0_Covariance[0] + K_S0_Covariance[2],1/2)>>h_SVcov_rxy_KS0");
0107   XiDataTree->Draw("pow(Ximinus_Covariance[0] + Ximinus_Covariance[2],1/2)>>h_SVcov_rxy_Xi",Xi_charge_cut.c_str());
0108
0109   XiDataTree->Draw("Lambda0_y:Lambda0_x>>h_Lambda_x_y_Xi",Xi_charge_cut.c_str());
0110   XiDataTree->Draw("sqrt(pow(Lambda0_x,2) + pow(Lambda0_y,2)):Lambda0_z>>h_Lambda_z_r_Xi",Xi_charge_cut.c_str());
0111
0112   XiDataTree->Draw("sqrt(Lambda0_Covariance[2]):sqrt(Lambda0_Covariance[0])>>h_Lambdacov_x_y_Xi",Xi_charge_cut.c_str());
0113   XiDataTree->Draw("sqrt(Lambda0_Covariance[5]):sqrt(Lambda0_Covariance[0])>>h_Lambdacov_x_z_Xi",Xi_charge_cut.c_str());
0114
0115   XiDataTree->Draw("pow((3*Lambda0_vertex_volume)/(4*pi),(1./3.))>>h_Lambdavolume_r_Xi",Xi_charge_cut.c_str());
0116   XiDataTree->Draw("pow(Lambda0_Covariance[0]+Lambda0_Covariance[2]+Lambda0_Covariance[5],1/2)>>h_Lambdacov_r_Xi",Xi_charge_cut.c_str());
0117   XiDataTree->Draw("pow(Lambda0_Covariance[0]+Lambda0_Covariance[2],1/2)>>h_Lambdacov_rxy_Xi",Xi_charge_cut.c_str());
0118
0119   KS0DataTree->Draw("sqrt(track_1_Covariance[0]+track_1_Covariance[2]+track_1_Covariance[5])>>h_pion_track1_r_KS0");
0120   KS0DataTree->Draw("sqrt(track_2_Covariance[0]+track_2_Covariance[2]+track_2_Covariance[5])>>h_pion_track2_r_KS0");
0121   KS0DataTree->Draw("sqrt(track_1_Covariance[0]+track_1_Covariance[2])>>h_pion_track1_rxy_KS0");
0122   KS0DataTree->Draw("sqrt(track_2_Covariance[0]+track_2_Covariance[2])>>h_pion_track2_rxy_KS0");
0123
0124   XiDataTree->Draw("sqrt(Lambda0_track_1_Covariance[0]+Lambda0_track_1_Covariance[2]+Lambda0_track_1_Covariance[5])>>h_pion_track_1_r_Xi",Xi_charge_cut.c_str());
0125   XiDataTree->Draw("sqrt(Lambda0_track_2_Covariance[0]+Lambda0_track_2_Covariance[2]+Lambda0_track_2_Covariance[5])>>h_proton_track_2_r_Xi",Xi_charge_cut.c_str());
0126   XiDataTree->Draw("sqrt(track_3_Covariance[0]+track_3_Covariance[2]+track_3_Covariance[5])>>h_pion_track_3_r_Xi",Xi_charge_cut.c_str());
0127   XiDataTree->Draw("sqrt(Lambda0_track_1_Covariance[0]+Lambda0_track_1_Covariance[2])>>h_pion_track_1_rxy_Xi",Xi_charge_cut.c_str());
0128   XiDataTree->Draw("sqrt(Lambda0_track_2_Covariance[0]+Lambda0_track_2_Covariance[2])>>h_proton_track_2_rxy_Xi",Xi_charge_cut.c_str());
0129   XiDataTree->Draw("sqrt(track_3_Covariance[0]+track_3_Covariance[2])>>h_pion_track_3_rxy_Xi",Xi_charge_cut.c_str());
0130   */
0131
0132   TF1* fitFunc1 = new TF1("gaus1", "gaus", 0.002,0.0075);
0133   fitFunc1->SetLineColor(kRed);
0134   fitFunc1->SetParName(0, "Const"); fitFunc1->SetParameter(0, 5e5); fitFunc1->SetParLimits(0, 0, 1e6);
0135   fitFunc1->SetParName(1, "Mean"); fitFunc1->SetParameter(1, 0.0049); fitFunc1->SetParLimits(1, 0.001, 0.01);
0136   fitFunc1->SetParName(2, "Width"); fitFunc1->SetParameter(2, 0.001); fitFunc1->SetParLimits(2, 1e-4, 0.01);
0137   h_PVcov_r_KS0->Fit(fitFunc1, "ER");
0138   TF1* fitFunc2 = new TF1("gaus2", "gaus", 0.001,0.007);
0139   fitFunc2->SetLineColor(kRed);
0140   fitFunc2->SetParName(0, "Const"); fitFunc2->SetParameter(0, 5.5e3); fitFunc2->SetParLimits(0, 0, 1e4);
0141   fitFunc2->SetParName(1, "Mean"); fitFunc2->SetParameter(1, 0.0039); fitFunc2->SetParLimits(1, 0.001, 0.01);
0142   fitFunc2->SetParName(2, "Width"); fitFunc2->SetParameter(2, 0.001); fitFunc2->SetParLimits(2, 1e-4, 0.01);
0143   h_PVcov_rxy_Xi->Fit(fitFunc2, "ER");
0144   /*
0145   TF1* fitFunc3 = new TF1("gaus3", "gaus", 0.0,0.003);
0146   fitFunc3->SetLineColor(kRed);
0147   fitFunc3->SetParName(0, "Const"); fitFunc3->SetParameter(0, 9e5); fitFunc3->SetParLimits(0, 0, 5e6);
0148   fitFunc3->SetParName(1, "Mean"); fitFunc3->SetParameter(1, 0.001); fitFunc3->SetParLimits(1, 1e-4, 0.003);
0149   fitFunc3->SetParName(2, "Width"); fitFunc3->SetParameter(2, 0.001); fitFunc3->SetParLimits(2, 1e-4, 0.01);
0150   h_SVvolume_r_KS0->Fit(fitFunc3, "ER");
0151   TF1* fitFunc4 = new TF1("gaus4", "gaus", 0.0,0.005);
0152   fitFunc4->SetLineColor(kRed);
0153   fitFunc4->SetParName(0, "Const"); fitFunc4->SetParameter(0, 3.5e3); fitFunc4->SetParLimits(0, 0, 1e4);
0154   fitFunc4->SetParName(1, "Mean"); fitFunc4->SetParameter(1, 0.002); fitFunc4->SetParLimits(1, 0.001, 0.004);
0155   fitFunc4->SetParName(2, "Width"); fitFunc4->SetParameter(2, 0.001); fitFunc4->SetParLimits(2, 1e-4, 0.01);
0156   h_SVvolume_r_Xi->Fit(fitFunc4, "ER");
0157   TF1* fitFunc5 = new TF1("gaus5", "gaus", 0.0,0.005);
0158   fitFunc5->SetLineColor(kRed);
0159   fitFunc5->SetParName(0, "Const"); fitFunc5->SetParameter(0, 3.5e3); fitFunc5->SetParLimits(0, 0, 1e4);
0160   fitFunc5->SetParName(1, "Mean"); fitFunc5->SetParameter(1, 0.002); fitFunc5->SetParLimits(1, 0.001, 0.004);
0161   fitFunc5->SetParName(2, "Width"); fitFunc5->SetParameter(2, 0.001); fitFunc5->SetParLimits(2, 1e-4, 0.01);
0162   h_Lambdavolume_r_Xi->Fit(fitFunc5, "ER");
0163
0164   double fitMean1 = fitFunc1->GetParameter(1);
0165   double fitWidth1 = fitFunc1->GetParameter(2);
0166   double fitMean2 = fitFunc2->GetParameter(1);
0167   double fitWidth2 = fitFunc2->GetParameter(2);
0168   double fitMean3 = fitFunc3->GetParameter(1);
0169   double fitWidth3 = fitFunc3->GetParameter(2);
0170   double fitMean4 = fitFunc4->GetParameter(1);
0171   double fitWidth4 = fitFunc4->GetParameter(2);
0172   double fitMean5 = fitFunc5->GetParameter(1);
0173   double fitWidth5 = fitFunc5->GetParameter(2);
0174
0175   std::cout << "K_S0 PV mean and sigma : " << fitMean1 << "+/-" << fitWidth1 << std::endl;
0176   std::cout << "Xi PV mean and sigma : " << fitMean2 << "+/-" << fitWidth2 << std::endl;
0177   std::cout << "K_S0 SV mean and sigma : " << fitMean3 << "+/-" << fitWidth3 << std::endl;
0178   std::cout << "Xi SV mean and sigma : " << fitMean4 << "+/-" << fitWidth4 << std::endl;
0179   std::cout << "Xi Lambda mean and sigma : " << fitMean5 << "+/-" << fitWidth5 << std::endl;
0180   */
0181
0182   TFile* fout = new TFile("VertexDistributionHistograms_temp.root","RECREATE");
0183   fout->cd();
0184
0185   /*
0186   h_PV_x_y_KS0->Write();
0187   h_PV_x_y_Xi->Write();
0188   h_PV_z_r_KS0->Write();
0189   h_PV_z_r_Xi->Write();
0190   h_PVcov_x_y_KS0->Write();
0191   h_PVcov_x_y_Xi->Write();
0192   h_PVcov_x_z_KS0->Write();
0193   h_PVcov_x_z_Xi->Write();
0194   */
0195   h_PVcov_r_KS0->Write();
0196   /*
0197   h_PVcov_r_Xi->Write();
0198   h_PVcov_rxy_KS0->Write();
0199   */
0200   h_PVcov_rxy_Xi->Write();
0201   /*
0202   h_PVvolume_r_KS0->Write();
0203   h_PVvolume_r_Xi->Write();
0204   h_SV_x_y_KS0->Write();
0205   h_SV_x_y_Xi->Write();
0206   h_SV_z_r_KS0->Write();
0207   h_SV_z_r_Xi->Write();
0208   h_SVcov_x_y_KS0->Write();
0209   h_SVcov_x_y_Xi->Write();
0210   h_SVcov_x_z_KS0->Write();
0211   h_SVcov_x_z_Xi->Write();
0212   h_SVcov_r_KS0->Write();
0213   h_SVcov_r_Xi->Write();
0214   h_SVcov_rxy_KS0->Write();
0215   h_SVcov_rxy_Xi->Write();
0216   h_SVvolume_r_KS0->Write();
0217   h_SVvolume_r_Xi->Write();
0218   h_Lambda_x_y_Xi->Write();
0219   h_Lambda_z_r_Xi->Write();
0220   h_Lambdacov_x_y_Xi->Write();
0221   h_Lambdacov_x_z_Xi->Write();
0222   h_Lambdavolume_r_Xi->Write();
0223   h_Lambdacov_r_Xi->Write();
0224   h_Lambdacov_rxy_Xi->Write();
0225   h_PVcov_r_nTracks_KS0->Write();
0226   h_pion_track1_r_KS0->Write();
0227   h_pion_track2_r_KS0->Write();
0228   h_pion_track1_rxy_KS0->Write();
0229   h_pion_track2_rxy_KS0->Write();
0230   h_pion_track_1_r_Xi->Write();
0231   h_proton_track_2_r_Xi->Write();
0232   h_pion_track_3_r_Xi->Write();
0233   h_pion_track_1_rxy_Xi->Write();
0234   h_proton_track_2_rxy_Xi->Write();
0235   h_pion_track_3_rxy_Xi->Write();
0236   */
0237   fout->Close();
0238
0239   std::cout << "Completed" << std::endl;
0240 }