489 _amps1(amps1),
_amps2(amps2),
_ampsSum(ampsSum),
_r(r),
_delta(delta),
_gamma(gamma),
_tau(tau),
_dGamma(dGamma),
_dm(dm),
_eff_tag(eff_tag),
_w(w),
494 _r_dt =
new RooRealVar(
"dt",
"per-candidate time resolution estimate",0., 0.25);
495 _r_mistag =
new RooRealVar(
"mistag",
"mistag",0.);
496 _r_f =
new RooCategory(
"qf",
"qf");
497 _r_f->defineType(
"h+", +1);
498 _r_f->defineType(
"h-", -1);
499 _r_f->setRange(
"Range_p",
"h+");
500 _r_f->setRange(
"Range_m",
"h-");
501 _r_q =
new RooCategory(
"qt",
"qt");
502 _r_q->defineType(
"B+", +1);
503 _r_q->defineType(
"B-", -1) ;
504 _r_q->defineType(
"untagged", 0);
510 RooRealConstant::value(1.),
511 RooRealConstant::value(1.),
514 RooRealConstant::value(0.),
515 RooRealConstant::value(1.),
516 RooRealConstant::value(0.),
517 RooRealConstant::value(0.),
518 RooRealConstant::value(0.),
520 RooRealConstant::value(0.),
521 RooRealConstant::value(0.),
522 RooRealConstant::value(0.));
528 RooRealConstant::value(1.),
529 RooRealConstant::value(-1.),
532 RooRealConstant::value(0.),
533 RooRealConstant::value(1.),
534 RooRealConstant::value(0.),
535 RooRealConstant::value(0.),
536 RooRealConstant::value(0.),
538 RooRealConstant::value(0.),
539 RooRealConstant::value(0.),
540 RooRealConstant::value(0.));
546 RooRealConstant::value(-2.),
547 RooRealConstant::value(-2.),
550 RooRealConstant::value(0.),
551 RooRealConstant::value(1.),
552 RooRealConstant::value(0.),
553 RooRealConstant::value(0.),
554 RooRealConstant::value(0.),
556 RooRealConstant::value(0.),
557 RooRealConstant::value(0.),
558 RooRealConstant::value(0.));
564 RooRealConstant::value(2.),
565 RooRealConstant::value(-2.),
568 RooRealConstant::value(0.),
569 RooRealConstant::value(1.),
570 RooRealConstant::value(0.),
571 RooRealConstant::value(0.),
572 RooRealConstant::value(0.),
574 RooRealConstant::value(0.),
575 RooRealConstant::value(0.),
576 RooRealConstant::value(0.));
583 vector<double> myBinning = knot_positions.getVector();
586 vector<double> values = knot_values.getVector() ;
589 RooArgList tacc_list;
590 for(
int i= 0; i< values.size(); i++){
595 RooFormulaVar* coeff_last =
new RooFormulaVar((
"coeff_"+
anythingToString((
int)values.size()+1)).c_str(),(
"coeff_"+
anythingToString((
int)values.size()+1)).c_str(),
"@0 + ((@0-@1)/(@2-@3)) * (@4 - @2)", RooArgList(RooRealConstant::value(1.0), *tacc_list.find((
"coeff_"+
anythingToString((
int)values.size()-1)).c_str()) , RooRealConstant::value(myBinning[myBinning.size()-1]), RooRealConstant::value(myBinning[myBinning.size()-2]), RooRealConstant::value(
_r_t->getMax()) ));
597 tacc_list.add(*coeff_last);
605 for (
int i = 1; i<=h_spline->GetNbinsX(); i++) {
606 _r_t->setVal(h_spline->GetXaxis()->GetBinCenter(i));
607 h_spline->SetBinContent(i,
_spline->getVal());
610 TCanvas* c =
new TCanvas();
611 h_spline->SetLineColor(kRed);
612 h_spline->Draw(
"histc");
613 c->Print(
"spline.eps");
614 c->Print(
"spline.pdf");
617 _pdf_sigma_t =
new RooGenericPdf(
"pdf_sigma_t",
"pow(7. / @1, 7) / 720. * pow(@0, 6) * exp(-7. * @0 / @1)",RooArgList(*
_r_dt, RooRealConstant::value(0.04)));
DecRateCoeff_Bd * _cosh_coeff
NamedParameter< double > _min_TAU
AmpsPdfFlexiFast * _ampsSum
RooCubicSplineFun * _spline
DecRateCoeff_Bd * _sin_coeff
complex< double > _intAAbar
RooGaussEfficiencyModel * _efficiency
NamedParameter< double > _max_TAU
std::string anythingToString(const T &anything)
DecRateCoeff_Bd * _sinh_coeff
RooRealVar * _r_scale_mean_dt
cosh/sinh/cos/sin coefficients in decay rate equations
RooRealVar * _r_scale_sigma_dt
AmpsPdfFlexiFast * _amps1
AmpsPdfFlexiFast * _amps2
DecRateCoeff_Bd * _cos_coeff