+ nameRfl+=i;
+ nameMassNocutsS="hMassS_";
+ nameMassNocutsS+=i;
+ nameMassNocutsB="hMassB_";
+ nameMassNocutsB+=i;
+
+ //histograms of cut variable distributions
+ if(fFillVarHists){
+ if(fReadMC){
+ // dca
+ namedistr="hdcaS_";
+ namedistr+=i;
+ TH1F *hdcaS = new TH1F(namedistr.Data(), "DCA distribution;dca [cm]",200,0.,0.1);
+ // impact parameter
+ namedistr="hd0piS_";
+ namedistr+=i;
+ TH1F *hd0piS = new TH1F(namedistr.Data(), "Impact parameter distribution (pions);d0(#pi) [cm]",200,-0.1,0.1);
+
+ namedistr="hd0KS_";
+ namedistr+=i;
+ TH1F *hd0KS = new TH1F(namedistr.Data(), "Impact parameter distribution (kaons);d0(K) [cm]",200,-0.1,0.1);
+ namedistr="hd0d0S_";
+ namedistr+=i;
+ TH1F *hd0d0S = new TH1F(namedistr.Data(), "d_{0}#timesd_{0} distribution;d_{0}#timesd_{0} [cm^{2}]",200,-0.001,0.001);
+
+ //decay lenght
+ namedistr="hdeclS_";
+ namedistr+=i;
+ TH1F *hdeclengthS = new TH1F(namedistr.Data(), "Decay Length^{2} distribution;Decay Length^{2} [cm]",200,0,0.015);
+
+ namedistr="hnormdeclS_";
+ namedistr+=i;
+ TH1F *hnormdeclengthS = new TH1F(namedistr.Data(), "Normalized Decay Length^{2} distribution;(Decay Length/Err)^{2} ",200,0,12.);
+
+ namedistr="hdeclxyS_";
+ namedistr+=i;
+ TH1F* hdeclxyS=new TH1F(namedistr.Data(),"Decay Length XY distribution;Decay Length XY [cm]",200,0,0.15);
+ namedistr="hnormdeclxyS_";
+ namedistr+=i;
+ TH1F* hnormdeclxyS=new TH1F(namedistr.Data(),"Normalized decay Length XY distribution;Decay Length XY/Err",200,0,6.);
+
+ namedistr="hdeclxyd0d0S_";
+ namedistr+=i;
+ TH2F* hdeclxyd0d0S=new TH2F(namedistr.Data(),"Correlation decay Length XY - d_{0}#times d_{0};Decay Length XY [cm];d_{0}#times d_{0}[cm^{2}]",200,0,0.15,200,-0.001,0.001);
+
+ namedistr="hnormdeclxyd0d0S_";
+ namedistr+=i;
+ TH2F* hnormdeclxyd0d0S=new TH2F(namedistr.Data(),"Correlation normalized decay Length XY - d_{0}#times d_{0};Decay Length XY/Err;d_{0}#times d_{0}[cm^{2}]",200,0,6,200,-0.001,0.001);
+
+ // costhetapoint
+ namedistr="hcosthetapointS_";
+ namedistr+=i;
+ TH1F *hcosthetapointS = new TH1F(namedistr.Data(), "cos#theta_{Point} distribution;cos#theta_{Point}",200,0,1.);
+
+ namedistr="hcosthetapointxyS_";
+ namedistr+=i;
+ TH1F *hcosthetapointxyS = new TH1F(namedistr.Data(), "cos#theta_{Point} XYdistribution;cos#theta_{Point}",300,0.,1.);
+
+ TH1F* tmpMS = new TH1F(nameMassNocutsS.Data(),"D^{0} invariant mass; M [GeV]; Entries",300,1.5648,2.1648); //range (MD0-300MeV, mD0 + 300MeV)
+ tmpMS->Sumw2();
+
+
+ fDistr->Add(hdcaS);
+
+ fDistr->Add(hd0piS);
+ fDistr->Add(hd0KS);
+
+ fDistr->Add(hd0d0S);
+
+ fDistr->Add(hcosthetapointS);
+
+ fDistr->Add(hcosthetapointxyS);
+
+ fDistr->Add(hdeclengthS);
+
+ fDistr->Add(hnormdeclengthS);
+
+ fDistr->Add(hdeclxyS);
+
+ fDistr->Add(hnormdeclxyS);
+
+ fDistr->Add(hdeclxyd0d0S);
+ fDistr->Add(hnormdeclxyd0d0S);
+
+ fDistr->Add(tmpMS);
+ }
+
+ // dca
+ namedistr="hdcaB_";
+ namedistr+=i;
+ TH1F *hdcaB = new TH1F(namedistr.Data(), "DCA distribution;dca [cm]",200,0.,0.1);
+
+ // impact parameter
+ namedistr="hd0B_";
+ namedistr+=i;
+ TH1F *hd0B = new TH1F(namedistr.Data(), "Impact parameter distribution (both);d0 [cm]",200,-0.1,0.1);
+
+ namedistr="hd0d0B_";
+ namedistr+=i;
+ TH1F *hd0d0B = new TH1F(namedistr.Data(), "d_{0}#timesd_{0} distribution;d_{0}#timesd_{0} [cm^{2}]",200,-0.001,0.001);
+
+ //decay lenght
+ namedistr="hdeclB_";
+ namedistr+=i;
+ TH1F *hdeclengthB = new TH1F(namedistr.Data(), "Decay Length^{2} distribution;Decay Length^{2} [cm^{2}]",200,0,0.015);
+
+ namedistr="hnormdeclB_";
+ namedistr+=i;
+ TH1F *hnormdeclengthB = new TH1F(namedistr.Data(), "Normalized Decay Length distribution;(Decay Length/Err)^{2} ",200,0,12.);
+
+ namedistr="hdeclxyB_";
+ namedistr+=i;
+ TH1F* hdeclxyB=new TH1F(namedistr.Data(),"Decay Length XY distribution;Decay Length XY [cm]",200,0,0.15);
+ namedistr="hnormdeclxyB_";
+ namedistr+=i;
+ TH1F* hnormdeclxyB=new TH1F(namedistr.Data(),"Normalized decay Length XY distribution;Decay Length XY/Err",200,0,6.);
+
+ namedistr="hdeclxyd0d0B_";
+ namedistr+=i;
+ TH2F* hdeclxyd0d0B=new TH2F(namedistr.Data(),"Correlation decay Length XY - d_{0}#times d_{0};Decay Length XY [cm];d_{0}#times d_{0}[cm^{2}]",200,0,0.15,200,-0.001,0.001);
+
+ namedistr="hnormdeclxyd0d0B_";
+ namedistr+=i;
+ TH2F* hnormdeclxyd0d0B=new TH2F(namedistr.Data(),"Correlation normalized decay Length XY - d_{0}#times d_{0};Decay Length XY/Err;d_{0}#times d_{0}[cm^{2}]",200,0,6,200,-0.001,0.001);
+
+ // costhetapoint
+ namedistr="hcosthetapointB_";
+ namedistr+=i;
+ TH1F *hcosthetapointB = new TH1F(namedistr.Data(), "cos#theta_{Point} distribution;cos#theta_{Point}",200,0,1.);
+
+ namedistr="hcosthetapointxyB_";
+ namedistr+=i;
+ TH1F *hcosthetapointxyB = new TH1F(namedistr.Data(), "cos#theta_{Point} XY distribution;cos#theta_{Point} XY",300,0.,1.);
+
+ TH1F* tmpMB = new TH1F(nameMassNocutsB.Data(),"D^{0} invariant mass; M [GeV]; Entries",300,1.5648,2.1648); //range (MD0-300MeV, mD0 + 300MeV)
+ tmpMB->Sumw2();
+
+
+
+ fDistr->Add(hdcaB);
+
+ fDistr->Add(hd0B);
+
+ fDistr->Add(hd0d0B);
+
+ fDistr->Add(hcosthetapointB);
+
+ fDistr->Add(hcosthetapointxyB);
+
+ fDistr->Add(hdeclengthB);
+
+ fDistr->Add(hnormdeclengthB);
+
+ fDistr->Add(hdeclxyB);