+ //Pad3: invariant mass fit
+ c4->cd(3);
+ TH1D *hvar24;
+ if (icasType == 0) hvar24 = (TH1D*) hlist->FindObject(Form("fHistMassXiMinus"));
+ else if (icasType == 1) hvar24 = (TH1D*) hlist->FindObject(Form("fHistMassXiPlus"));
+ else if (icasType == 2) hvar24 = (TH1D*) hlist->FindObject(Form("fHistMassOmegaMinus"));
+ else if (icasType == 3) hvar24 = (TH1D*) hlist->FindObject(Form("fHistMassOmegaPlus"));
+ hvar24->Draw("histo");
+ TF1* fitinvmassnoPID = new TF1("fitinvmassnoPID","gaus(0)+pol2(3)",lowlimmass,uplimmass);
+ fitinvmassnoPID->SetParName(0, "cnstntG");
+ fitinvmassnoPID->SetParName(1, "meanG");
+ fitinvmassnoPID->SetParName(2, "sigmaG");
+ fitinvmassnoPID->SetParLimits(0,0.,500000.);
+ if (icasType==0||icasType==1) {
+ fitinvmassnoPID->SetParameter(1, 1.32171);
+ fitinvmassnoPID->SetParLimits(1, 1.31,1.33);
+ fitinvmassnoPID->SetParLimits(2,0.001,0.005);
+ } else {
+ fitinvmassnoPID->SetParameter(1, 1.67245);
+ fitinvmassnoPID->SetParLimits(1, 1.664,1.68);
+ fitinvmassnoPID->SetParLimits(2,0.0008,0.006);
+ }
+ hvar24->Fit("fitinvmassnoPID","rimeN");
+ fitinvmassnoPID->SetLineColor(kRed);
+ fitinvmassnoPID->Draw("same");
+ Float_t meanGaussnoPID = fitinvmassnoPID->GetParameter(1);
+ Float_t sigmaGaussnoPID = fitinvmassnoPID->GetParameter(2);
+ cout<<"Mean: "<<meanGaussnoPID<<endl;
+ cout<<"Sigma: "<<sigmaGaussnoPID<<endl;
+ //Pad4: Text
+ c4->cd(4);
+ Float_t refwidth = 0.002;
+ if (icasType > 1) refwidth = 0.0025;
+ TPaveText *pave2 = new TPaveText(0.05,0.3,0.95,0.7);
+ pave2->SetFillColor(0);
+ pave2->SetTextSize(0.04);
+ pave2->SetTextAlign(12);
+ pave2->AddText("Invariant mass distribution #color[4]{WITHOUT} PID on dauther tracks");
+ if (icasType < 2) pave2->AddText("PDG mass: 1.32171 GeV/c^{2}");
+ else pave2->AddText("PDG mass: 1.67245 GeV/c^{2}");
+ pave2->AddText(Form("#color[1]{Mass from Fit: %.5f #pm %.5f GeV/c^{2}}",meanGaussnoPID,sigmaGaussnoPID));
+ if (sigmaGaussnoPID > refwidth - 0.0003 && sigmaGaussnoPID < refwidth + 0.0003) pave2->AddText("#color[3]{OK!! The width is compatible with standard.}");
+ else pave2->AddText("#color[2]{NOT OK!! Problem.}");
+ pave2->Draw();
+ cout<<" "<<refwidth - 0.0003<<"<"<<sigmaGaussnoPID<<"<"<<refwidth + 0.0003<<endl;
+ c4->SaveAs("fig_lf_Multistrange.pdf");
+