TimePdfMaster Class Reference

#include <TimePdfMaster.h>

Public Member Functions
	TimePdfMaster (const MINT::FitParameter &Gamma, const MINT::FitParameter &dGamma, const MINT::FitParameter &dm, const MINT::FitParameter &offset_sigma_dt, const MINT::FitParameter &scale_mean_dt, const MINT::FitParameter &scale_sigma_dt, const MINT::FitParameter &scale_sigma_2_dt, const MINT::FitParameter &c0, const MINT::FitParameter &c1, const MINT::FitParameter &c2, const MINT::FitParameter &c3, const MINT::FitParameter &c4, const MINT::FitParameter &c5, const MINT::FitParameter &c6, const MINT::FitParameter &c7, const MINT::FitParameter &c8, const MINT::FitParameter &c9, const MINT::FitParameter &p0_os, const MINT::FitParameter &p1_os, const MINT::FitParameter &delta_p0_os, const MINT::FitParameter &delta_p1_os, const MINT::FitParameter &avg_eta_os, const MINT::FitParameter &tageff_os, const MINT::FitParameter &tageff_asym_os, const MINT::FitParameter &p0_ss, const MINT::FitParameter &p1_ss, const MINT::FitParameter &delta_p0_ss, const MINT::FitParameter &delta_p1_ss, const MINT::FitParameter &avg_eta_ss, const MINT::FitParameter &tageff_ss, const MINT::FitParameter &tageff_asym_ss, const MINT::FitParameter &production_asym, const MINT::FitParameter &detection_asym, string marginalPdfsPrefix)
RooDataSet *	sampleEvents (int N=10000)
vector< double >	getRandom_marginalVals ()
void	fillProtoData (double dt, int f, int q_OS, double eta_OS, int q_SS, double eta_SS)
double	getSamplingPdfVal (IDalitzEvent &evt)
double	get_cosh_term_Val (IDalitzEvent &evt)
double	get_cosh_term_Integral (IDalitzEvent &evt)
double	get_sinh_term_Val (IDalitzEvent &evt)
double	get_sinh_term_Integral (IDalitzEvent &evt)
double	get_cos_term_Val (IDalitzEvent &evt)
double	get_cos_term_Integral (IDalitzEvent &evt)
double	get_sin_term_Val (IDalitzEvent &evt)
double	get_sin_term_Integral (IDalitzEvent &evt)
double	get_cosh_coeff_Val (IDalitzEvent &evt)
double	get_cos_coeff_Val (IDalitzEvent &evt)
double	get_sinh_coeff_Val (IDalitzEvent &evt)
double	get_sin_coeff_Val (IDalitzEvent &evt)
double	get_marginalPdfs_Val (IDalitzEvent &evt)
double	get_marginalPdfs_product (IDalitzEvent &evt)
double	get_spline_Val (IDalitzEvent &evt)
double	get_tau_Val ()
double	get_dm_Val ()
double	get_dGamma_Val ()
std::pair< double, double >	getCalibratedMistag_OS (IDalitzEvent &evt)
std::pair< double, double >	getCalibratedMistag_OS (double &eta_OS)
std::pair< double, double >	getCalibratedMistag_OS (IDalitzEvent &evt, double &avg_eta_os, double &p0_os, double &p1_os, double &delta_p0_os, double &delta_p1_os)
std::pair< double, double >	getCalibratedMistag_SS (IDalitzEvent &evt)
std::pair< double, double >	getCalibratedMistag_SS (IDalitzEvent &evt, double &avg_eta_ss, double &p0_ss, double &p1_ss, double &delta_p0_ss, double &delta_p1_ss)
std::pair< double, double >	getCalibratedMistag (double eta, double avg_eta, double p0, double p1, double delta_p0, double delta_p1)
std::pair< double, double >	getCalibratedMistag_SS (double &eta_SS)
double	getCalibratedResolution (double &dt)
double	getCalibratedResolution (double &dt, double &scale_sigma_dt, double &scale_sigma_2_dt)
void	listFitParDependencies ()
void	setAllFitParameters ()
void	setAllObservables (IDalitzEvent &evt)
void	setAllObservablesToMean (IDalitzEvent &evt)
void	setCP_coeff (double norm, double norm_bar, double C, double C_bar, double D, double D_bar, double S, double S_bar)
vector< double >	getCP_coeff ()
void	setAllObservablesAndFitParameters (IDalitzEvent &evt)
TH1D *	plotSpline ()
virtual	~TimePdfMaster ()

Protected Attributes
const MINT::FitParameter &	_Gamma
const MINT::FitParameter &	_dGamma
const MINT::FitParameter &	_dm
const MINT::FitParameter &	_offset_sigma_dt
const MINT::FitParameter &	_scale_mean_dt
const MINT::FitParameter &	_scale_sigma_dt
const MINT::FitParameter &	_scale_sigma_2_dt
const MINT::FitParameter &	_c0
const MINT::FitParameter &	_c1
const MINT::FitParameter &	_c2
const MINT::FitParameter &	_c3
const MINT::FitParameter &	_c4
const MINT::FitParameter &	_c5
const MINT::FitParameter &	_c6
const MINT::FitParameter &	_c7
const MINT::FitParameter &	_c8
const MINT::FitParameter &	_c9
const MINT::FitParameter &	_p0_os
const MINT::FitParameter &	_p1_os
const MINT::FitParameter &	_delta_p0_os
const MINT::FitParameter &	_delta_p1_os
const MINT::FitParameter &	_avg_eta_os
const MINT::FitParameter &	_tageff_os
const MINT::FitParameter &	_tageff_asym_os
const MINT::FitParameter &	_p0_ss
const MINT::FitParameter &	_p1_ss
const MINT::FitParameter &	_delta_p0_ss
const MINT::FitParameter &	_delta_p1_ss
const MINT::FitParameter &	_avg_eta_ss
const MINT::FitParameter &	_tageff_ss
const MINT::FitParameter &	_tageff_asym_ss
const MINT::FitParameter &	_production_asym
const MINT::FitParameter &	_detection_asym
RooRealVar *	_r_t
	Cast of MINT parameters to B2DX parameters. More...
RooRealVar *	_r_dt
RooRealVar *	_r_dt_scaled
RooCategory *	_r_q_OS
RooRealVar *	_r_eta_OS
RooCategory *	_r_q_SS
RooRealVar *	_r_eta_SS
RooCategory *	_r_f
RooCategory *	_r_run
RooCategory *	_r_trigger
RooRealVar *	_r_tau
RooRealVar *	_r_dGamma
RooRealVar *	_r_dm
RooRealVar *	_r_scale_mean_dt
RooRealVar *	_r_offset_sigma_dt
RooRealVar *	_r_scale_sigma_dt
RooRealVar *	_r_scale_sigma_2_dt
RooRealVar *	_r_c0
RooRealVar *	_r_c1
RooRealVar *	_r_c2
RooRealVar *	_r_c3
RooRealVar *	_r_c4
RooRealVar *	_r_c5
RooRealVar *	_r_c6
RooRealVar *	_r_c7
RooRealVar *	_r_c8
RooRealVar *	_r_c9
RooFormulaVar *	_coeff_last
RooRealVar *	_r_p0_os
RooRealVar *	_r_p1_os
RooRealVar *	_r_delta_p0_os
RooRealVar *	_r_delta_p1_os
RooRealVar *	_r_avg_eta_os
RooRealVar *	_r_tageff_os
RooRealVar *	_r_tageff_asym_os
RooRealVar *	_r_p0_ss
RooRealVar *	_r_p1_ss
RooRealVar *	_r_delta_p0_ss
RooRealVar *	_r_delta_p1_ss
RooRealVar *	_r_avg_eta_ss
RooRealVar *	_r_tageff_ss
RooRealVar *	_r_tageff_asym_ss
RooRealVar *	_r_production_asym
RooRealVar *	_r_detection_asym
RooCubicSplineFun *	_spline
RooProduct *	_splinePdf
RooGaussEfficiencyModel *	_efficiency
RooRealVar *	_r_norm
RooRealVar *	_r_norm_bar
RooRealVar *	_r_C
RooRealVar *	_r_C_bar
RooRealVar *	_r_D
RooRealVar *	_r_D_bar
RooRealVar *	_r_S
RooRealVar *	_r_S_bar
DecRateCoeff_Bd *	_cos_coeff
DecRateCoeff_Bd *	_cosh_coeff
DecRateCoeff_Bd *	_sin_coeff
DecRateCoeff_Bd *	_sinh_coeff
TimePdfIntegrator *	_cosh_term
TimePdfIntegrator *	_sinh_term
TimePdfIntegrator *	_cos_term
TimePdfIntegrator *	_sin_term
string	_marginalPdfsPrefix
TFile *	_f_pdfs
TH1D *	_h_dt
RooDataHist *	_r_h_dt
RooAbsPdf *	_pdf_sigma_t
TH1D *	_h_q_f
TH1D *	_h_q_OS
TH1D *	_h_q_SS
RooAbsPdf *	_pdf_eta_OS
TH1D *	_h_eta_OS
RooDataHist *	_r_h_eta_OS
RooAbsPdf *	_pdf_eta_SS
TH1D *	_h_eta_SS
RooDataHist *	_r_h_eta_SS
NamedParameter< double >	_min_TAU
NamedParameter< double >	_max_TAU
NamedParameter< double >	_min_TAUERR
NamedParameter< double >	_max_TAUERR
RooBDecay *	_samplingPdf_t
RooBDecay *	_fitPdf_t
RooGaussModel *	_resmodel
RooEffProd *	_samplingPdf_t_eff
RooProdPdf *	_samplingPdf
RooProdPdf *	_fitPdf
RooDataSet *	_protoData
RooMCStudy *	_sampleGen

Detailed Description

Definition at line 52 of file TimePdfMaster.h.

Constructor & Destructor Documentation

TimePdfMaster()

TimePdfMaster::TimePdfMaster ( const MINT::FitParameter &  Gamma, const MINT::FitParameter &  dGamma, const MINT::FitParameter &  dm, const MINT::FitParameter &  offset_sigma_dt, const MINT::FitParameter &  scale_mean_dt, const MINT::FitParameter &  scale_sigma_dt, const MINT::FitParameter &  scale_sigma_2_dt, const MINT::FitParameter &  c0, const MINT::FitParameter &  c1, const MINT::FitParameter &  c2, const MINT::FitParameter &  c3, const MINT::FitParameter &  c4, const MINT::FitParameter &  c5, const MINT::FitParameter &  c6, const MINT::FitParameter &  c7, const MINT::FitParameter &  c8, const MINT::FitParameter &  c9, const MINT::FitParameter &  p0_os, const MINT::FitParameter &  p1_os, const MINT::FitParameter &  delta_p0_os, const MINT::FitParameter &  delta_p1_os, const MINT::FitParameter &  avg_eta_os, const MINT::FitParameter &  tageff_os, const MINT::FitParameter &  tageff_asym_os, const MINT::FitParameter &  p0_ss, const MINT::FitParameter &  p1_ss, const MINT::FitParameter &  delta_p0_ss, const MINT::FitParameter &  delta_p1_ss, const MINT::FitParameter &  avg_eta_ss, const MINT::FitParameter &  tageff_ss, const MINT::FitParameter &  tageff_asym_ss, const MINT::FitParameter &  production_asym, const MINT::FitParameter &  detection_asym, string  marginalPdfsPrefix  )

inline

Init B2DX parameters

Acceptance

Decay rate coefficients

Time pdf integrators

Marginal pdfs

Definition at line 211 of file TimePdfMaster.h.

                                                                                                                              : 
        _Gamma(Gamma),
        _dGamma(dGamma),
        _dm(dm),
        _offset_sigma_dt(offset_sigma_dt),    
        _scale_mean_dt(scale_mean_dt),
        _scale_sigma_dt(scale_sigma_dt),
        _scale_sigma_2_dt(scale_sigma_2_dt),
        _c0(c0),
        _c1(c1),
        _c2(c2),
        _c3(c3),
        _c4(c4),
        _c5(c5),
        _c6(c6),
        _c7(c7),
        _c8(c8),
        _c9(c9),
        _p0_os(p0_os),
        _p1_os(p1_os),
        _delta_p0_os(delta_p0_os),
        _delta_p1_os(delta_p1_os),
        _avg_eta_os(avg_eta_os),
        _tageff_os(tageff_os),
        _tageff_asym_os(tageff_asym_os),
        _p0_ss(p0_ss),
        _p1_ss(p1_ss),
        _delta_p0_ss(delta_p0_ss),
        _delta_p1_ss(delta_p1_ss),
        _avg_eta_ss(avg_eta_ss),
        _tageff_ss(tageff_ss),
        _tageff_asym_ss(tageff_asym_ss),
        _production_asym(production_asym),
        _detection_asym(detection_asym),
        _marginalPdfsPrefix(marginalPdfsPrefix),
        _min_TAU("min_TAU", 0.4),
        _max_TAU("max_TAU", 10.),
        _min_TAUERR("min_TAUERR", 0.),
        _max_TAUERR("max_TAUERR", 0.1)
    {
        
        // Observables
        _r_t = new RooRealVar("t", "t", _min_TAU, _max_TAU);
        _r_dt = new RooRealVar("dt", "dt",_min_TAUERR,_max_TAUERR);
        _r_f = new RooCategory("qf", "qf");
        _r_f->defineType("h+", +1);
        _r_f->defineType("h-", -1);
        _r_f->setRange("Range_p","h+");
        _r_f->setRange("Range_m","h-");
        _r_q_OS = new RooCategory("q_OS", "q_OS");
        _r_q_OS->defineType("B+", +1);
        _r_q_OS->defineType("B-", -1) ;
        _r_q_OS->defineType("untagged", 0);
        _r_eta_OS = new RooRealVar("eta_OS", "eta_OS",0.,0.5);
        _r_q_SS = new RooCategory("q_SS", "q_SS");
        _r_q_SS->defineType("B+", +1);
        _r_q_SS->defineType("B-", -1) ;
        _r_q_SS->defineType("untagged", 0);
        _r_eta_SS = new RooRealVar("eta_SS", "eta_SS",0.,0.5);
        _r_run = new RooCategory("run", "run");
        _r_run->defineType("Run1", 1);
        _r_run->defineType("Run2", 2);
        _r_trigger = new RooCategory("trigger", "trigger");
        _r_trigger->defineType("t0", 0);
        _r_trigger->defineType("t1", 1);

        // Fit parameters
        _r_tau = new RooRealVar("tau", "tau",1./_Gamma);
        _r_dGamma = new RooRealVar("dGamma", "dGamma",_dGamma);
        _r_dm = new RooRealVar("dm", "dm",_dm);
        
        _r_scale_mean_dt = new RooRealVar("scale_mean_dt", "scale_mean_dt", _scale_mean_dt);
        _r_offset_sigma_dt = new RooRealVar("offset_sigma_dt", "offset_sigma_dt", _offset_sigma_dt);
        _r_scale_sigma_dt = new RooRealVar("scale_sigma_dt", "scale_sigma_dt", _scale_sigma_dt);
    _r_scale_sigma_2_dt = new RooRealVar("scale_sigma_2_dt", "scale_sigma_2_dt", _scale_sigma_2_dt);      

        _r_c0 = new RooRealVar("coeff_0", "coeff_0",_c0);
        _r_c1 = new RooRealVar("coeff_1", "coeff_1",_c1);
        _r_c2 = new RooRealVar("coeff_2", "coeff_2",_c2);
        _r_c3 = new RooRealVar("coeff_3", "coeff_3",_c3);
        _r_c4 = new RooRealVar("coeff_4", "coeff_4",_c4);
        _r_c5 = new RooRealVar("coeff_5", "coeff_5",_c5);
        _r_c6 = new RooRealVar("coeff_6", "coeff_6",_c6);
        _r_c7 = new RooRealVar("coeff_7", "coeff_7",_c7);
        _r_c8 = new RooRealVar("coeff_8", "coeff_8",_c8);
        _r_c9 = new RooRealVar("coeff_9", "coeff_9",_c9);
        
        _r_p0_os = new RooRealVar("p0_os", "p0_os",_p0_os);
        _r_p1_os = new RooRealVar("p1_os", "p1_os",_p1_os);
        _r_delta_p0_os = new RooRealVar("delta_p0_os", "delta_p0_os",_delta_p0_os);
        _r_delta_p1_os = new RooRealVar("delta_p1_os", "delta_p1_os",_delta_p1_os);
        _r_avg_eta_os = new RooRealVar("avg_eta_os", "avg_eta_os",_avg_eta_os);
        _r_tageff_os = new RooRealVar("tageff_os", "tageff_os",_tageff_os);
        _r_tageff_asym_os = new RooRealVar("tageff_asym_os", "tageff_asym_os",_tageff_asym_os);
        _r_p0_ss = new RooRealVar("p0_ss", "p0_ss",_p0_ss);
        _r_p1_ss = new RooRealVar("p1_ss", "p1_ss",_p1_ss);
        _r_delta_p0_ss = new RooRealVar("delta_p0_ss", "delta_p0_ss",_delta_p0_ss);
        _r_delta_p1_ss = new RooRealVar("delta_p1_ss", "delta_p1_ss",_delta_p1_ss);
        _r_avg_eta_ss = new RooRealVar("avg_eta_ss", "avg_eta_ss",_avg_eta_ss);
        _r_tageff_ss = new RooRealVar("tageff_ss", "tageff_ss",_tageff_ss);
        _r_tageff_asym_ss= new RooRealVar("tageff_asym_ss", "tageff_asym_ss",_tageff_asym_ss);
        
        _r_production_asym = new RooRealVar("production_asym", "production_asym",_production_asym);
        _r_detection_asym = new RooRealVar("detection_asym", "detection_asym",_detection_asym);
        
        vector<RooRealVar*> v_coeff;
        v_coeff.push_back(_r_c0);
        v_coeff.push_back(_r_c1);
        v_coeff.push_back(_r_c2);
        v_coeff.push_back(_r_c3);
        v_coeff.push_back(_r_c4);
        v_coeff.push_back(_r_c5);
        v_coeff.push_back(_r_c6);
        v_coeff.push_back(_r_c7);
        v_coeff.push_back(_r_c8);
        v_coeff.push_back(_r_c9);
 
        //SPLINE KNOTS
        NamedParameter<double> knot_positions("knot_positions", 0.5, 1., 1.5, 2., 3., 6., 9.5);
        vector<double> myBinning = knot_positions.getVector();
        
        RooArgList tacc_list;
        for(int i= 0; i<= myBinning.size(); i++){
            tacc_list.add(*v_coeff[i]);
        }
        
        _coeff_last = new RooFormulaVar(("coeff_"+anythingToString((int)myBinning.size()+1)).c_str(),("coeff_"+anythingToString((int)myBinning.size()+1)).c_str(), "@0 + ((@0-@1)/(@2-@3)) * (@4 - @2)", RooArgList(RooRealConstant::value(1.0), *tacc_list.find(("coeff_"+anythingToString((int)myBinning.size()-1)).c_str())  , RooRealConstant::value(myBinning[myBinning.size()-1]), RooRealConstant::value(myBinning[myBinning.size()-2]), RooRealConstant::value(_r_t->getMax()) ));
        
        tacc_list.add(*_coeff_last);    
        
        // Cubic spline function
        _spline = new RooCubicSplineFun("spline", "spline", *_r_t, myBinning, tacc_list);        
        //_splinePdf = new RooCubicSplinePdf("splinePdf", "splinePdf", *_r_t, myBinning, tacc_list);
    _splinePdf = new RooProduct("splinePdf", "splinePdf", RooArgList(*_spline,RooRealConstant::value(0.1)));

        //_efficiency = new RooGaussEfficiencyModel("resmodel", "resmodel", *_r_t, *_spline, RooRealConstant::value(0.), *_r_dt, *_r_scale_mean_dt, *_r_scale_sigma_dt );
        _r_dt_scaled = (RooRealVar*) new RooFormulaVar( "r_dt_scaled","r_dt_scaled", "@0+@1*@2+@3*@2*@2",RooArgList(*_r_offset_sigma_dt,*_r_scale_sigma_dt,*_r_dt,*_r_scale_sigma_2_dt));
        _efficiency = new RooGaussEfficiencyModel("resmodel", "resmodel", *_r_t, *_spline, RooRealConstant::value(0.), *_r_dt_scaled, *_r_scale_mean_dt, RooRealConstant::value(1.) );
        
        // CP coefficients
        _r_norm = new RooRealVar("norm", "norm",1);
        _r_norm_bar = new RooRealVar("norm_bar", "norm_bar",1);
        _r_C = new RooRealVar("C", "C",1);
        _r_C_bar = new RooRealVar("Cbar", "Cbar",-1);
        //_r_C_bar= (RooRealVar*) new RooFormulaVar("Cbar","-1. * @0",RooArgList(*_r_C));
        _r_D= new RooRealVar("D", "D",0);
        _r_D_bar= new RooRealVar("Dbar", "Dbar",0);
        _r_S= new RooRealVar("S", "S",0);
        _r_S_bar= new RooRealVar("Sbar", "Sbar",0);
        
        _cosh_coeff = new DecRateCoeff_Bd("cosh_coeff",
                                          "cosh_coeff",
                                          DecRateCoeff_Bd::kCosh ,
                                          *_r_f,
                                          *_r_norm,
                                          *_r_norm_bar,
                                          *_r_q_OS,
                                          *_r_eta_OS,
                                          *_r_p0_os,
                                          *_r_p1_os,
                                          *_r_delta_p0_os,
                                          *_r_delta_p1_os,
                                          *_r_avg_eta_os,
                                          *_r_tageff_os,
                                          *_r_tageff_asym_os,
                                          *_r_q_SS,
                                          *_r_eta_SS,
                                          *_r_p0_ss,
                                          *_r_p1_ss,
                                          *_r_delta_p0_ss,
                                          *_r_delta_p1_ss,
                                          *_r_avg_eta_ss,
                                          *_r_tageff_ss,
                                          *_r_tageff_asym_ss,
                                          *_r_production_asym,
                                          *_r_detection_asym );
        
        _cos_coeff = new DecRateCoeff_Bd("cos_coeff",
                                         "cos_coeff",
                                         DecRateCoeff_Bd::kCos ,
                                         *_r_f,
                                         *_r_C,
                                         *_r_C_bar,
                                         *_r_q_OS,
                                         *_r_eta_OS,
                                         *_r_p0_os,
                                         *_r_p1_os,
                                         *_r_delta_p0_os,
                                         *_r_delta_p1_os,
                                         *_r_avg_eta_os,
                                         *_r_tageff_os,
                                         *_r_tageff_asym_os,
                                         *_r_q_SS,
                                         *_r_eta_SS,
                                         *_r_p0_ss,
                                         *_r_p1_ss,
                                         *_r_delta_p0_ss,
                                         *_r_delta_p1_ss,
                                         *_r_avg_eta_ss,
                                         *_r_tageff_ss,
                                         *_r_tageff_asym_ss,
                                         *_r_production_asym,
                                         *_r_detection_asym );
        
        _sinh_coeff = new DecRateCoeff_Bd("sinh_coeff",
                                          "sinh_coeff",
                                          DecRateCoeff_Bd::kSinh ,
                                          *_r_f,
                                          *_r_D,
                                          *_r_D_bar,
                                          *_r_q_OS,
                                          *_r_eta_OS,
                                          *_r_p0_os,
                                          *_r_p1_os,
                                          *_r_delta_p0_os,
                                          *_r_delta_p1_os,
                                          *_r_avg_eta_os,
                                          *_r_tageff_os,
                                          *_r_tageff_asym_os,
                                          *_r_q_SS,
                                          *_r_eta_SS,
                                          *_r_p0_ss,
                                          *_r_p1_ss,
                                          *_r_delta_p0_ss,
                                          *_r_delta_p1_ss,
                                          *_r_avg_eta_ss,
                                          *_r_tageff_ss,
                                          *_r_tageff_asym_ss,
                                          *_r_production_asym,
                                          *_r_detection_asym );
        
        _sin_coeff = new DecRateCoeff_Bd("sin_coeff",
                                         "sin_coeff",
                                         DecRateCoeff_Bd::kSin,
                                         *_r_f,
                                         *_r_S,
                                         *_r_S_bar,
                                         *_r_q_OS,
                                         *_r_eta_OS,
                                         *_r_p0_os,
                                         *_r_p1_os,
                                         *_r_delta_p0_os,
                                         *_r_delta_p1_os,
                                         *_r_avg_eta_os,
                                         *_r_tageff_os,
                                         *_r_tageff_asym_os,
                                         *_r_q_SS,
                                         *_r_eta_SS,
                                         *_r_p0_ss,
                                         *_r_p1_ss,
                                         *_r_delta_p0_ss,
                                         *_r_delta_p1_ss,
                                         *_r_avg_eta_ss,
                                         *_r_tageff_ss,
                                         *_r_tageff_asym_ss,
                                         *_r_production_asym,
                                         *_r_detection_asym );
        
        _cosh_term = new TimePdfIntegrator(coshBasis,_efficiency,
                                           _Gamma,_dGamma,_dm,
                                           _offset_sigma_dt, _scale_mean_dt, _scale_sigma_dt, _scale_sigma_2_dt
                                           ,_c0, _c1, _c2
                                           ,_c3, _c4, _c5
                                           ,_c6, _c7, _c8
                                           ,_c9);
        _cos_term = new TimePdfIntegrator(cosBasis,_efficiency,
                                          _Gamma,_dGamma,_dm,
                                          _offset_sigma_dt, _scale_mean_dt, _scale_sigma_dt, _scale_sigma_2_dt
                                          ,_c0, _c1, _c2
                                          ,_c3, _c4, _c5
                                          ,_c6, _c7, _c8
                                          ,_c9);
        _sinh_term = new TimePdfIntegrator(sinhBasis,_efficiency,
                                           _Gamma,_dGamma,_dm,
                                           _offset_sigma_dt, _scale_mean_dt, _scale_sigma_dt, _scale_sigma_2_dt
                                           ,_c0, _c1, _c2
                                           ,_c3, _c4, _c5
                                           ,_c6, _c7, _c8
                                           ,_c9);
        _sin_term = new TimePdfIntegrator(sinBasis,_efficiency,
                                          _Gamma,_dGamma,_dm,
                                          _offset_sigma_dt, _scale_mean_dt, _scale_sigma_dt, _scale_sigma_2_dt
                                          ,_c0, _c1, _c2
                                          ,_c3, _c4, _c5
                                          ,_c6, _c7, _c8
                                          ,_c9);
        
    cout << (string)_marginalPdfsPrefix << endl;
    if((string)_marginalPdfsPrefix == "Uniform"){
        _pdf_sigma_t = (RooAbsPdf*) (new RooUniform("pdf_sigma_t","pdf_sigma_t",*_r_dt));
        _pdf_eta_OS = (RooAbsPdf*) (new RooUniform("pdf_eta_OS","pdf_eta_OS",*_r_eta_OS));
        _pdf_eta_SS = (RooAbsPdf*) (new RooUniform("pdf_eta_SS","pdf_eta_SS",*_r_eta_SS));
    }
    else {
        _f_pdfs = new TFile("Mistag_pdfs.root","OPEN");
    
        _h_dt = new TH1D( *((TH1D*) _f_pdfs->Get(("h_dt_norm_"+(string)_marginalPdfsPrefix).c_str())));
        // RooHistPdf doesn't like negative or 0 bins (can happen due to sweights), set them to a small positive number
        _h_dt->Smooth();
        for(int i= 1 ; i<=_h_dt->GetNbinsX(); i++){
            if(_h_dt->GetBinContent(i) <= 0.)_h_dt->SetBinContent(i,0.000000001*_h_dt->GetMaximum());
        }
        _r_h_dt = new RooDataHist("r_h_dt","r_h_dt",*_r_dt,_h_dt);
        _pdf_sigma_t = (RooAbsPdf*) (new RooHistPdf("pdf_sigma_t","pdf_sigma_t",*_r_dt,*_r_h_dt));
        
        _h_eta_OS = new TH1D( *((TH1D*) _f_pdfs->Get(("h_w_OS_norm_"+(string)_marginalPdfsPrefix).c_str())));
        _h_eta_OS->Smooth();
        for(int i= 1 ; i<=_h_eta_OS->GetNbinsX(); i++){
            if(_h_eta_OS->GetBinContent(i) <= 0.)_h_eta_OS->SetBinContent(i,0.000000001*_h_eta_OS->GetMaximum());
        }
        _r_h_eta_OS = new RooDataHist("r_eta_OS","r_eta_OS",*_r_eta_OS,_h_eta_OS);
        _pdf_eta_OS = (RooAbsPdf*) (new RooHistPdf("pdf_eta_OS","pdf_eta_OS",*_r_eta_OS,*_r_h_eta_OS));
        
        _h_eta_SS = new TH1D( *((TH1D*) _f_pdfs->Get(("h_w_SS_norm_"+(string)_marginalPdfsPrefix).c_str())));
        _h_eta_SS->Smooth();
        for(int i= 1 ; i<=_h_eta_SS->GetNbinsX(); i++){
            if(_h_eta_SS->GetBinContent(i) <= 0.)_h_eta_SS->SetBinContent(i,0.000000001*_h_eta_SS->GetMaximum());
        }
        _r_h_eta_SS = new RooDataHist("r_eta_SS","r_eta_SS",*_r_eta_SS,_h_eta_SS);
        _pdf_eta_SS = (RooAbsPdf*) (new RooHistPdf("pdf_eta_SS","pdf_eta_SS",*_r_eta_SS,*_r_h_eta_SS));

//      _h_q_f = new TH1D( *((TH1D*) _f_pdfs->Get(("h_q_f_norm_"+(string)_marginalPdfsPrefix).c_str())));
//      _h_q_OS = new TH1D( *((TH1D*) _f_pdfs->Get(("h_q_OS_norm_"+(string)_marginalPdfsPrefix).c_str())));
//      _h_q_SS = new TH1D( *((TH1D*) _f_pdfs->Get(("h_q_SS_norm_"+(string)_marginalPdfsPrefix).c_str())));
        }

    _resmodel = new RooGaussModel("resmodel", "resmodel", *_r_t,  RooRealConstant::value(0.), *_r_dt_scaled);              

    _samplingPdf_t = new RooBDecay(("samplingPdf_t"+(string)_marginalPdfsPrefix).c_str(), ("samplingPdf_t"+(string)_marginalPdfsPrefix).c_str(),
                *_r_t,*_r_tau, *_r_dGamma, 
                                *_cosh_coeff,*_sinh_coeff,*_cos_coeff,*_sin_coeff,
                *_r_dm, *_resmodel, RooBDecay::SingleSided); 

    _samplingPdf_t_eff = new RooEffProd(("samplingPdf_t_eff"+(string)_marginalPdfsPrefix).c_str(), ("samplingPdf_t_eff"+(string)_marginalPdfsPrefix).c_str(),
                         *_samplingPdf_t, *_splinePdf);
    
        _samplingPdf = new RooProdPdf(("samplingPdf"+(string)_marginalPdfsPrefix).c_str(),("samplingPdf"+(string)_marginalPdfsPrefix).c_str(),
                RooArgSet(*_pdf_sigma_t,*_pdf_eta_OS,*_pdf_eta_SS),Conditional(RooArgSet(*_samplingPdf_t_eff),RooArgSet(*_r_t,*_r_f,*_r_q_OS,*_r_q_SS)));


    _fitPdf_t = new RooBDecay(("fitPdf_t"+(string)_marginalPdfsPrefix).c_str(), ("fitPdf_t"+(string)_marginalPdfsPrefix).c_str(),
                *_r_t,*_r_tau, *_r_dGamma, 
                                *_cosh_coeff,*_sinh_coeff,*_cos_coeff,*_sin_coeff,
                *_r_dm, *_efficiency, RooBDecay::SingleSided); 

        _fitPdf = new RooProdPdf(("fitPdf"+(string)_marginalPdfsPrefix).c_str(),("fitPdf"+(string)_marginalPdfsPrefix).c_str(),
                RooArgSet(*_pdf_sigma_t,*_pdf_eta_OS,*_pdf_eta_SS),Conditional(RooArgSet(*_fitPdf_t),RooArgSet(*_r_t,*_r_f,*_r_q_OS,*_r_q_SS)));

    _sampleGen = 0;
    }

~TimePdfMaster()

virtual TimePdfMaster::~TimePdfMaster ( )

inlinevirtual

Definition at line 879 of file TimePdfMaster.h.

                           {
       // Plot acceptance
       TH1D* h_spline = plotSpline();       
       TCanvas* c = new TCanvas();
       h_spline->SetLineColor(kRed);
       h_spline->Draw("histc");
       c->Print("spline.eps");
       //c->Print("spline.pdf");
   }

Member Function Documentation

fillProtoData()

void TimePdfMaster::fillProtoData ( double  dt, int  f, int  q_OS, double  eta_OS, int  q_SS, double  eta_SS  )

inline

Definition at line 602 of file TimePdfMaster.h.

                                                                                          {
        _r_dt->setVal(dt);        
        _r_q_SS->setIndex(q_SS);
        _r_q_OS->setIndex(q_OS);
        _r_f->setIndex(f);
        _r_eta_OS->setVal(eta_OS);
        _r_eta_SS->setVal(eta_SS);
        _protoData->add(RooArgSet(*_r_dt,*_r_f,*_r_q_OS,*_r_eta_OS,*_r_q_SS,*_r_eta_SS));
    }

get_cos_coeff_Val()

double TimePdfMaster::get_cos_coeff_Val ( IDalitzEvent &  evt )

inline

Definition at line 653 of file TimePdfMaster.h.

                                               {
        return _cos_coeff->evaluate();
    }

get_cos_term_Integral()

double TimePdfMaster::get_cos_term_Integral ( IDalitzEvent &  evt )

inline

Definition at line 637 of file TimePdfMaster.h.

                                                   {
        return _cos_coeff->analyticalIntegral(2) * _cos_term->getVal(evt).imag();
    }

get_cos_term_Val()

double TimePdfMaster::get_cos_term_Val ( IDalitzEvent &  evt )

inline

Definition at line 633 of file TimePdfMaster.h.

                                              {
        return _cos_coeff->evaluate() * _cos_term->getVal(evt).real();
    }

get_cosh_coeff_Val()

double TimePdfMaster::get_cosh_coeff_Val ( IDalitzEvent &  evt )

inline

Definition at line 649 of file TimePdfMaster.h.

                                                {
        return _cosh_coeff->evaluate();
    }

get_cosh_term_Integral()

double TimePdfMaster::get_cosh_term_Integral ( IDalitzEvent &  evt )

inline

Definition at line 621 of file TimePdfMaster.h.

                                                    {
        return _cosh_coeff->analyticalIntegral(2) * _cosh_term->getVal(evt).imag();
    }

get_cosh_term_Val()

double TimePdfMaster::get_cosh_term_Val ( IDalitzEvent &  evt )

inline

Definition at line 617 of file TimePdfMaster.h.

                                               {
        return _cosh_coeff->evaluate() * _cosh_term->getVal(evt).real();
    }

get_dGamma_Val()

double TimePdfMaster::get_dGamma_Val ( )

inline

Definition at line 711 of file TimePdfMaster.h.

                           {
        return (double)_dGamma;
    }

get_dm_Val()

double TimePdfMaster::get_dm_Val ( )

inline

Definition at line 707 of file TimePdfMaster.h.

                       {
        return (double)_dm;
    }

get_marginalPdfs_product()

double TimePdfMaster::get_marginalPdfs_product ( IDalitzEvent &  evt )

inline

Definition at line 695 of file TimePdfMaster.h.

                                                      {
        return  _pdf_sigma_t->getVal(RooArgSet(*_r_dt)) * _pdf_eta_OS->getVal(RooArgSet(*_r_eta_OS)) * _pdf_eta_SS->getVal(RooArgSet(*_r_eta_SS));
    }

get_marginalPdfs_Val()

double TimePdfMaster::get_marginalPdfs_Val ( IDalitzEvent &  evt )

inline

Definition at line 665 of file TimePdfMaster.h.

                                                  {
        const int q_OS = (int)evt.getValueFromVector(3);
        const int q_SS = (int)evt.getValueFromVector(5);
        
        /*
        // check norm 
        TRandom3 r(0);
        double sum1 = 0.;
        double sum2 = 0.;
        double sum3 = 0.;

        for (int i=0; i < 1000000; i++) {
            _r_eta_OS->setVal(r.Uniform(0,0.5));
            _r_eta_SS->setVal(r.Uniform(0,0.5));
            _r_dt->setVal(r.Uniform(0,0.1));
            sum1 += _pdf_eta_OS->getVal(RooArgSet(*_r_eta_OS));
            sum2 += _pdf_eta_SS->getVal(RooArgSet(*_r_eta_SS));
            sum3 += _pdf_sigma_t->getVal(RooArgSet(*_r_dt));
        }
        cout << "norm = " << 0.5 * sum1/1000000. << endl;
        cout << "norm = " << 0.5 * sum2/1000000. << endl;
        cout << "norm = " << 0.1 * sum3/1000000. << endl;
        throw "";
        */
        
        return  _pdf_sigma_t->getVal(RooArgSet(*_r_dt))
                * ( abs(q_OS) * _pdf_eta_OS->getVal(RooArgSet(*_r_eta_OS)) + ( 1. - abs(q_OS)) * 1./0.5 )
                * ( abs(q_SS) * _pdf_eta_SS->getVal(RooArgSet(*_r_eta_SS)) + ( 1. - abs(q_SS)) * 1./0.5 );
    }

get_sin_coeff_Val()

double TimePdfMaster::get_sin_coeff_Val ( IDalitzEvent &  evt )

inline

Definition at line 661 of file TimePdfMaster.h.

                                               {
        return _sin_coeff->evaluate();
    }

get_sin_term_Integral()

double TimePdfMaster::get_sin_term_Integral ( IDalitzEvent &  evt )

inline

Definition at line 645 of file TimePdfMaster.h.

                                                   {
        return _sin_coeff->analyticalIntegral(2) * _sin_term->getVal(evt).imag();
    }

get_sin_term_Val()

double TimePdfMaster::get_sin_term_Val ( IDalitzEvent &  evt )

inline

Definition at line 641 of file TimePdfMaster.h.

                                              {
        return _sin_coeff->evaluate() * _sin_term->getVal(evt).real();
    }

get_sinh_coeff_Val()

double TimePdfMaster::get_sinh_coeff_Val ( IDalitzEvent &  evt )

inline

Definition at line 657 of file TimePdfMaster.h.

                                                {
        return _sinh_coeff->evaluate();
    }

get_sinh_term_Integral()

double TimePdfMaster::get_sinh_term_Integral ( IDalitzEvent &  evt )

inline

Definition at line 629 of file TimePdfMaster.h.

                                                    {
        return _sinh_coeff->analyticalIntegral(2) * _sinh_term->getVal(evt).imag();
    }

get_sinh_term_Val()

double TimePdfMaster::get_sinh_term_Val ( IDalitzEvent &  evt )

inline

Definition at line 625 of file TimePdfMaster.h.

                                               {
        return _sinh_coeff->evaluate() * _sinh_term->getVal(evt).real();
    }

get_spline_Val()

double TimePdfMaster::get_spline_Val ( IDalitzEvent &  evt )

inline

Definition at line 699 of file TimePdfMaster.h.

                                            {
        return  _spline->getVal();
    }

get_tau_Val()

double TimePdfMaster::get_tau_Val ( )

inline

Definition at line 703 of file TimePdfMaster.h.

                        {
        return (double) 1./_Gamma;
    }

getCalibratedMistag()

std::pair<double, double> TimePdfMaster::getCalibratedMistag ( double  eta, double  avg_eta, double  p0, double  p1, double  delta_p0, double  delta_p1  )

inline

Definition at line 735 of file TimePdfMaster.h.

                                                                                                                               {
        return _cosh_coeff->calibrate(eta, avg_eta, p0, p1, delta_p0, delta_p1);
    }

getCalibratedMistag_OS() [1/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_OS ( IDalitzEvent &  evt )

inline

Definition at line 715 of file TimePdfMaster.h.

                                                                     {
        return _cosh_coeff->calibrate(evt.getValueFromVector(4), _avg_eta_os, _p0_os, _p1_os, _delta_p0_os, _delta_p1_os);
    }

getCalibratedMistag_OS() [2/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_OS ( double &  eta_OS )

inline

Definition at line 719 of file TimePdfMaster.h.

                                                                  {
        return _cosh_coeff->calibrate(eta_OS, _avg_eta_os, _p0_os, _p1_os, _delta_p0_os, _delta_p1_os);
    }

getCalibratedMistag_OS() [3/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_OS ( IDalitzEvent &  evt, double &  avg_eta_os, double &  p0_os, double &  p1_os, double &  delta_p0_os, double &  delta_p1_os  )

inline

Definition at line 723 of file TimePdfMaster.h.

                                                                                                                                                             {
        return _cosh_coeff->calibrate(evt.getValueFromVector(4), avg_eta_os, p0_os, p1_os, delta_p0_os, delta_p1_os);
    }

getCalibratedMistag_SS() [1/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_SS ( IDalitzEvent &  evt )

inline

Definition at line 727 of file TimePdfMaster.h.

                                                                     {
        return _cosh_coeff->calibrate(evt.getValueFromVector(6), _avg_eta_ss, _p0_ss, _p1_ss, _delta_p0_ss, _delta_p1_ss);
    }

getCalibratedMistag_SS() [2/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_SS ( IDalitzEvent &  evt, double &  avg_eta_ss, double &  p0_ss, double &  p1_ss, double &  delta_p0_ss, double &  delta_p1_ss  )

inline

Definition at line 731 of file TimePdfMaster.h.

                                                                                                                                                             {
        return _cosh_coeff->calibrate(evt.getValueFromVector(6), avg_eta_ss, p0_ss, p1_ss, delta_p0_ss, delta_p1_ss);
    }

getCalibratedMistag_SS() [3/3]

std::pair<double, double> TimePdfMaster::getCalibratedMistag_SS ( double &  eta_SS )

inline

Definition at line 739 of file TimePdfMaster.h.

                                                                  {
        return _cosh_coeff->calibrate(eta_SS, _avg_eta_ss, _p0_ss, _p1_ss, _delta_p0_ss, _delta_p1_ss);
    }

getCalibratedResolution() [1/2]

double TimePdfMaster::getCalibratedResolution ( double &  dt )

inline

Definition at line 743 of file TimePdfMaster.h.

                                              {
        return _offset_sigma_dt + _scale_sigma_dt * dt + _scale_sigma_2_dt * dt *dt;
    }

getCalibratedResolution() [2/2]

double TimePdfMaster::getCalibratedResolution ( double &  dt, double &  scale_sigma_dt, double &  scale_sigma_2_dt  )

inline

Definition at line 747 of file TimePdfMaster.h.

                                                                                              {
        return _offset_sigma_dt + scale_sigma_dt * dt + scale_sigma_2_dt * dt *dt;
    }

getCP_coeff()

vector<double> TimePdfMaster::getCP_coeff ( )

inline

Definition at line 840 of file TimePdfMaster.h.

                                {

    vector<double> vals;
        vals.push_back(_r_norm->getVal());
        vals.push_back(_r_norm_bar->getVal());
        vals.push_back(_r_C->getVal());
        vals.push_back(_r_C_bar->getVal());
        vals.push_back(_r_D->getVal());
        vals.push_back(_r_D_bar->getVal());
        vals.push_back(_r_S->getVal());
        vals.push_back(_r_S_bar->getVal());
    
    return vals;
    }

getRandom_marginalVals()

vector<double> TimePdfMaster::getRandom_marginalVals ( )

inline

Definition at line 594 of file TimePdfMaster.h.

                                           {
    vector<double> vals;
    vals.push_back(_h_dt->GetRandom());
    vals.push_back(_h_eta_OS->GetRandom());
    vals.push_back(_h_eta_SS->GetRandom());
    return vals;
    }

getSamplingPdfVal()

double TimePdfMaster::getSamplingPdfVal ( IDalitzEvent &  evt )

inline

Definition at line 612 of file TimePdfMaster.h.

                                               {
    setAllObservablesAndFitParameters(evt);
    return _fitPdf->getVal(RooArgSet(*_r_t,*_r_dt,*_r_eta_OS,*_r_eta_SS,*_r_f,*_r_q_OS,*_r_q_SS));
    }

listFitParDependencies()

void TimePdfMaster::listFitParDependencies ( )

inline

Definition at line 751 of file TimePdfMaster.h.

                                {
       std::cout << "TimePDfIntegrator depends on these fitParameters:" << std::endl;
       _cosh_term->listFitParDependencies(std::cout);
       //std::cout << "TimePDfMaster depends on these fitParameters (sinh):" << std::endl;
       //_sinh_term->listFitParDependencies(std::cout);
       //std::cout << "" << std::endl;
   }

plotSpline()

TH1D* TimePdfMaster::plotSpline ( )

inline

Definition at line 860 of file TimePdfMaster.h.

                     {
       setAllFitParameters();
       TH1D *h_spline = new TH1D("", ";t (ps);Efficiency (a.u.)", 100, _min_TAU, _max_TAU);
       for (int i = 1; i<=h_spline->GetNbinsX(); i++) {
           _r_t->setVal(h_spline->GetXaxis()->GetBinCenter(i));
           h_spline->SetBinContent(i,_spline->getVal());
       }
       TCanvas* c = new TCanvas();
       h_spline->SetLineColor(kRed);
       h_spline->SetLineWidth(5);

       h_spline->SetMinimum(0.);
       h_spline->SetMaximum(1.2);

       h_spline->Draw("histc");
       c->Print("spline.eps");
       return h_spline;
   }

sampleEvents()

RooDataSet* TimePdfMaster::sampleEvents ( int  N = 10000 )

inline

Definition at line 577 of file TimePdfMaster.h.

                                           {
        if(_sampleGen == 0){
            _protoData = new RooDataSet("protoData","protoData",RooArgSet(*_r_dt,*_r_eta_OS,*_r_eta_SS));
        int N_proto = N > 100000 ? N : 100000;
        for(int i = 0 ; i < N; i++){
                _r_dt->setVal(_h_dt->GetRandom());        
                _r_eta_OS->setVal(_h_eta_OS->GetRandom());
                _r_eta_SS->setVal(_h_eta_SS->GetRandom());
                _protoData->add(RooArgSet(*_r_dt,*_r_eta_OS,*_r_eta_SS));
        }
        _sampleGen = new RooMCStudy(*_samplingPdf,RooArgSet(*_r_t,*_r_q_OS,*_r_q_SS,*_r_f),ProtoData(*_protoData,kFALSE,kTRUE));
    }
    _sampleGen->generate(1,N,kTRUE);
    RooDataSet* data = (RooDataSet*)_sampleGen->genData(0);
    return data;
    }

setAllFitParameters()

void TimePdfMaster::setAllFitParameters ( )

inline

Definition at line 759 of file TimePdfMaster.h.

                             {       
       _r_tau->setVal(1./_Gamma);
       _r_dGamma->setVal(_dGamma);
       _r_dm->setVal(_dm);
       
       _r_scale_mean_dt->setVal(_scale_mean_dt);
       _r_offset_sigma_dt->setVal(_offset_sigma_dt);
       _r_scale_sigma_dt->setVal(_scale_sigma_dt);
       _r_scale_sigma_2_dt->setVal(_scale_sigma_2_dt);
       _r_c0->setVal(_c0);
       _r_c1->setVal(_c1);
       _r_c2->setVal(_c2);
       _r_c3->setVal(_c3);
       _r_c4->setVal(_c4);
       _r_c5->setVal(_c5);
       _r_c6->setVal(_c6);
       _r_c7->setVal(_c7);
       _r_c8->setVal(_c8);
       _r_c9->setVal(_c9);
       
       _r_p0_os->setVal(_p0_os);
       _r_p1_os->setVal(_p1_os);
       _r_delta_p0_os->setVal(_delta_p0_os);
       _r_delta_p1_os->setVal(_delta_p1_os);
       _r_avg_eta_os->setVal(_avg_eta_os);
       _r_tageff_os->setVal(_tageff_os);
       _r_tageff_asym_os->setVal(_tageff_asym_os);
       _r_p0_ss->setVal(_p0_ss);
       _r_p1_ss->setVal(_p1_ss);
       _r_delta_p0_ss->setVal(_delta_p0_ss);
       _r_delta_p1_ss->setVal(_delta_p1_ss);
       _r_avg_eta_ss->setVal(_avg_eta_ss);
       _r_tageff_ss->setVal(_tageff_ss);
       _r_tageff_asym_ss->setVal(_tageff_asym_ss);
       
       _r_production_asym->setVal(_production_asym);
       _r_detection_asym->setVal(_detection_asym);
   }

setAllObservables()

void TimePdfMaster::setAllObservables ( IDalitzEvent &  evt )

inline

Definition at line 798 of file TimePdfMaster.h.

                                            {
       const double t = (double) evt.getValueFromVector(0);
       const double dt = (double) evt.getValueFromVector(1);
       const int f = (int)evt.getValueFromVector(2);
       const int q_OS = (int)evt.getValueFromVector(3);
       const double eta_OS = (double) evt.getValueFromVector(4);
       const int q_SS = (int)evt.getValueFromVector(5);
       const double eta_SS = (double) evt.getValueFromVector(6);
       
       _r_t->setVal(t);
       _r_dt->setVal(dt);
       _r_f->setIndex(f);
       _r_q_OS->setIndex(q_OS);
       _r_eta_OS->setVal(eta_OS);
       _r_q_SS->setIndex(q_SS);
       _r_eta_SS->setVal(eta_SS);
    }

setAllObservablesAndFitParameters()

void TimePdfMaster::setAllObservablesAndFitParameters ( IDalitzEvent &  evt )

inline

Definition at line 855 of file TimePdfMaster.h.

                                                             {
        setAllObservables(evt);
        setAllFitParameters();
    }

setAllObservablesToMean()

void TimePdfMaster::setAllObservablesToMean ( IDalitzEvent &  evt )

inline

Definition at line 816 of file TimePdfMaster.h.

                                                  {
       _r_t->setVal(1./_Gamma);
       if((string)_marginalPdfsPrefix == "Uniform")_r_dt->setVal(1./_Gamma/100.);
       else _r_dt->setVal(_h_dt->GetMean());
       _r_f->setIndex((int)evt.getValueFromVector(2));
       _r_q_OS->setIndex(0);
       if((string)_marginalPdfsPrefix == "Uniform")_r_eta_OS->setVal(0.);
       else _r_eta_OS->setVal(_h_eta_OS->GetMean());
       _r_q_SS->setIndex(0);
       if((string)_marginalPdfsPrefix == "Uniform")_r_eta_SS->setVal(0.);
       else _r_eta_SS->setVal(_h_eta_SS->GetMean());
    }

setCP_coeff()

void TimePdfMaster::setCP_coeff ( double  norm, double  norm_bar, double  C, double  C_bar, double  D, double  D_bar, double  S, double  S_bar  )

inline

Definition at line 829 of file TimePdfMaster.h.

                                                                                                                     {
        _r_norm->setVal(norm);
        _r_norm_bar->setVal(norm_bar);
        _r_C->setVal(C);
        _r_C_bar->setVal(C_bar);
        _r_D->setVal(D);
        _r_D_bar->setVal(D_bar);
        _r_S->setVal(S);
        _r_S_bar->setVal(S_bar);
    }

Member Data Documentation

_avg_eta_os

const MINT::FitParameter& TimePdfMaster::_avg_eta_os

protected

Definition at line 81 of file TimePdfMaster.h.

_avg_eta_ss

const MINT::FitParameter& TimePdfMaster::_avg_eta_ss

protected

Definition at line 89 of file TimePdfMaster.h.

_c0

const MINT::FitParameter& TimePdfMaster::_c0

protected

Definition at line 66 of file TimePdfMaster.h.

_c1

const MINT::FitParameter& TimePdfMaster::_c1

protected

Definition at line 67 of file TimePdfMaster.h.

_c2

const MINT::FitParameter& TimePdfMaster::_c2

protected

Definition at line 68 of file TimePdfMaster.h.

_c3

const MINT::FitParameter& TimePdfMaster::_c3

protected

Definition at line 69 of file TimePdfMaster.h.

_c4

const MINT::FitParameter& TimePdfMaster::_c4

protected

Definition at line 70 of file TimePdfMaster.h.

_c5

const MINT::FitParameter& TimePdfMaster::_c5

protected

Definition at line 71 of file TimePdfMaster.h.

_c6

const MINT::FitParameter& TimePdfMaster::_c6

protected

Definition at line 72 of file TimePdfMaster.h.

_c7

const MINT::FitParameter& TimePdfMaster::_c7

protected

Definition at line 73 of file TimePdfMaster.h.

_c8

const MINT::FitParameter& TimePdfMaster::_c8

protected

Definition at line 74 of file TimePdfMaster.h.

_c9

const MINT::FitParameter& TimePdfMaster::_c9

protected

Definition at line 75 of file TimePdfMaster.h.

_coeff_last

RooFormulaVar* TimePdfMaster::_coeff_last

protected

Definition at line 127 of file TimePdfMaster.h.

_cos_coeff

DecRateCoeff_Bd* TimePdfMaster::_cos_coeff

protected

Definition at line 163 of file TimePdfMaster.h.

_cos_term

TimePdfIntegrator* TimePdfMaster::_cos_term

protected

Definition at line 171 of file TimePdfMaster.h.

_cosh_coeff

DecRateCoeff_Bd* TimePdfMaster::_cosh_coeff

protected

Definition at line 164 of file TimePdfMaster.h.

_cosh_term

TimePdfIntegrator* TimePdfMaster::_cosh_term

protected

Definition at line 169 of file TimePdfMaster.h.

_delta_p0_os

const MINT::FitParameter& TimePdfMaster::_delta_p0_os

protected

Definition at line 79 of file TimePdfMaster.h.

_delta_p0_ss

const MINT::FitParameter& TimePdfMaster::_delta_p0_ss

protected

Definition at line 87 of file TimePdfMaster.h.

_delta_p1_os

const MINT::FitParameter& TimePdfMaster::_delta_p1_os

protected

Definition at line 80 of file TimePdfMaster.h.

_delta_p1_ss

const MINT::FitParameter& TimePdfMaster::_delta_p1_ss

protected

Definition at line 88 of file TimePdfMaster.h.

_detection_asym

const MINT::FitParameter& TimePdfMaster::_detection_asym

protected

Definition at line 94 of file TimePdfMaster.h.

_dGamma

const MINT::FitParameter& TimePdfMaster::_dGamma

protected

Definition at line 58 of file TimePdfMaster.h.

_dm

const MINT::FitParameter& TimePdfMaster::_dm

protected

Definition at line 59 of file TimePdfMaster.h.

_efficiency

RooGaussEfficiencyModel* TimePdfMaster::_efficiency

protected

Definition at line 150 of file TimePdfMaster.h.

_f_pdfs

TFile* TimePdfMaster::_f_pdfs

protected

Definition at line 176 of file TimePdfMaster.h.

_fitPdf

RooProdPdf* TimePdfMaster::_fitPdf

protected

Definition at line 206 of file TimePdfMaster.h.

_fitPdf_t

RooBDecay* TimePdfMaster::_fitPdf_t

protected

Definition at line 202 of file TimePdfMaster.h.

_Gamma

const MINT::FitParameter& TimePdfMaster::_Gamma

protected

Definition at line 57 of file TimePdfMaster.h.

_h_dt

TH1D* TimePdfMaster::_h_dt

protected

Definition at line 178 of file TimePdfMaster.h.

_h_eta_OS

TH1D* TimePdfMaster::_h_eta_OS

protected

Definition at line 187 of file TimePdfMaster.h.

_h_eta_SS

TH1D* TimePdfMaster::_h_eta_SS

protected

Definition at line 191 of file TimePdfMaster.h.

_h_q_f

TH1D* TimePdfMaster::_h_q_f

protected

Definition at line 182 of file TimePdfMaster.h.

_h_q_OS

TH1D* TimePdfMaster::_h_q_OS

protected

Definition at line 183 of file TimePdfMaster.h.

_h_q_SS

TH1D* TimePdfMaster::_h_q_SS

protected

Definition at line 184 of file TimePdfMaster.h.

_marginalPdfsPrefix

string TimePdfMaster::_marginalPdfsPrefix

protected

Definition at line 175 of file TimePdfMaster.h.

_max_TAU

NamedParameter<double> TimePdfMaster::_max_TAU

protected

Definition at line 196 of file TimePdfMaster.h.

_max_TAUERR

NamedParameter<double> TimePdfMaster::_max_TAUERR

protected

Definition at line 198 of file TimePdfMaster.h.

_min_TAU

NamedParameter<double> TimePdfMaster::_min_TAU

protected

Definition at line 195 of file TimePdfMaster.h.

_min_TAUERR

NamedParameter<double> TimePdfMaster::_min_TAUERR

protected

Definition at line 197 of file TimePdfMaster.h.

_offset_sigma_dt

const MINT::FitParameter& TimePdfMaster::_offset_sigma_dt

protected

Definition at line 61 of file TimePdfMaster.h.

_p0_os

const MINT::FitParameter& TimePdfMaster::_p0_os

protected

Definition at line 77 of file TimePdfMaster.h.

_p0_ss

const MINT::FitParameter& TimePdfMaster::_p0_ss

protected

Definition at line 85 of file TimePdfMaster.h.

_p1_os

const MINT::FitParameter& TimePdfMaster::_p1_os

protected

Definition at line 78 of file TimePdfMaster.h.

_p1_ss

const MINT::FitParameter& TimePdfMaster::_p1_ss

protected

Definition at line 86 of file TimePdfMaster.h.

_pdf_eta_OS

RooAbsPdf* TimePdfMaster::_pdf_eta_OS

protected

Definition at line 186 of file TimePdfMaster.h.

_pdf_eta_SS

RooAbsPdf* TimePdfMaster::_pdf_eta_SS

protected

Definition at line 190 of file TimePdfMaster.h.

_pdf_sigma_t

RooAbsPdf* TimePdfMaster::_pdf_sigma_t

protected

Definition at line 180 of file TimePdfMaster.h.

_production_asym

const MINT::FitParameter& TimePdfMaster::_production_asym

protected

Definition at line 93 of file TimePdfMaster.h.

_protoData

RooDataSet* TimePdfMaster::_protoData

protected

Definition at line 207 of file TimePdfMaster.h.

_r_avg_eta_os

RooRealVar* TimePdfMaster::_r_avg_eta_os

protected

Definition at line 133 of file TimePdfMaster.h.

_r_avg_eta_ss

RooRealVar* TimePdfMaster::_r_avg_eta_ss

protected

Definition at line 140 of file TimePdfMaster.h.

_r_C

RooRealVar* TimePdfMaster::_r_C

protected

Definition at line 155 of file TimePdfMaster.h.

_r_c0

RooRealVar* TimePdfMaster::_r_c0

protected

Definition at line 117 of file TimePdfMaster.h.

_r_c1

RooRealVar* TimePdfMaster::_r_c1

protected

Definition at line 118 of file TimePdfMaster.h.

_r_c2

RooRealVar* TimePdfMaster::_r_c2

protected

Definition at line 119 of file TimePdfMaster.h.

_r_c3

RooRealVar* TimePdfMaster::_r_c3

protected

Definition at line 120 of file TimePdfMaster.h.

_r_c4

RooRealVar* TimePdfMaster::_r_c4

protected

Definition at line 121 of file TimePdfMaster.h.

_r_c5

RooRealVar* TimePdfMaster::_r_c5

protected

Definition at line 122 of file TimePdfMaster.h.

_r_c6

RooRealVar* TimePdfMaster::_r_c6

protected

Definition at line 123 of file TimePdfMaster.h.

_r_c7

RooRealVar* TimePdfMaster::_r_c7

protected

Definition at line 124 of file TimePdfMaster.h.

_r_c8

RooRealVar* TimePdfMaster::_r_c8

protected

Definition at line 125 of file TimePdfMaster.h.

_r_c9

RooRealVar* TimePdfMaster::_r_c9

protected

Definition at line 126 of file TimePdfMaster.h.

_r_C_bar

RooRealVar* TimePdfMaster::_r_C_bar

protected

Definition at line 156 of file TimePdfMaster.h.

_r_D

RooRealVar* TimePdfMaster::_r_D

protected

Definition at line 157 of file TimePdfMaster.h.

_r_D_bar

RooRealVar* TimePdfMaster::_r_D_bar

protected

Definition at line 158 of file TimePdfMaster.h.

_r_delta_p0_os

RooRealVar* TimePdfMaster::_r_delta_p0_os

protected

Definition at line 131 of file TimePdfMaster.h.

_r_delta_p0_ss

RooRealVar* TimePdfMaster::_r_delta_p0_ss

protected

Definition at line 138 of file TimePdfMaster.h.

_r_delta_p1_os

RooRealVar* TimePdfMaster::_r_delta_p1_os

protected

Definition at line 132 of file TimePdfMaster.h.

_r_delta_p1_ss

RooRealVar* TimePdfMaster::_r_delta_p1_ss

protected

Definition at line 139 of file TimePdfMaster.h.

_r_detection_asym

RooRealVar* TimePdfMaster::_r_detection_asym

protected

Definition at line 145 of file TimePdfMaster.h.

_r_dGamma

RooRealVar* TimePdfMaster::_r_dGamma

protected

Definition at line 109 of file TimePdfMaster.h.

_r_dm

RooRealVar* TimePdfMaster::_r_dm

protected

Definition at line 110 of file TimePdfMaster.h.

_r_dt

RooRealVar* TimePdfMaster::_r_dt

protected

Definition at line 98 of file TimePdfMaster.h.

_r_dt_scaled

RooRealVar* TimePdfMaster::_r_dt_scaled

protected

Definition at line 99 of file TimePdfMaster.h.

_r_eta_OS

RooRealVar* TimePdfMaster::_r_eta_OS

protected

Definition at line 101 of file TimePdfMaster.h.

_r_eta_SS

RooRealVar* TimePdfMaster::_r_eta_SS

protected

Definition at line 103 of file TimePdfMaster.h.

_r_f

RooCategory* TimePdfMaster::_r_f

protected

Definition at line 104 of file TimePdfMaster.h.

_r_h_dt

RooDataHist* TimePdfMaster::_r_h_dt

protected

Definition at line 179 of file TimePdfMaster.h.

_r_h_eta_OS

RooDataHist* TimePdfMaster::_r_h_eta_OS

protected

Definition at line 188 of file TimePdfMaster.h.

_r_h_eta_SS

RooDataHist* TimePdfMaster::_r_h_eta_SS

protected

Definition at line 192 of file TimePdfMaster.h.

_r_norm

RooRealVar* TimePdfMaster::_r_norm

protected

Definition at line 153 of file TimePdfMaster.h.

_r_norm_bar

RooRealVar* TimePdfMaster::_r_norm_bar

protected

Definition at line 154 of file TimePdfMaster.h.

_r_offset_sigma_dt

RooRealVar* TimePdfMaster::_r_offset_sigma_dt

protected

Definition at line 113 of file TimePdfMaster.h.

_r_p0_os

RooRealVar* TimePdfMaster::_r_p0_os

protected

Definition at line 129 of file TimePdfMaster.h.

_r_p0_ss

RooRealVar* TimePdfMaster::_r_p0_ss

protected

Definition at line 136 of file TimePdfMaster.h.

_r_p1_os

RooRealVar* TimePdfMaster::_r_p1_os

protected

Definition at line 130 of file TimePdfMaster.h.

_r_p1_ss

RooRealVar* TimePdfMaster::_r_p1_ss

protected

Definition at line 137 of file TimePdfMaster.h.

_r_production_asym

RooRealVar* TimePdfMaster::_r_production_asym

protected

Definition at line 144 of file TimePdfMaster.h.

_r_q_OS

RooCategory* TimePdfMaster::_r_q_OS

protected

Definition at line 100 of file TimePdfMaster.h.

_r_q_SS

RooCategory* TimePdfMaster::_r_q_SS

protected

Definition at line 102 of file TimePdfMaster.h.

_r_run

RooCategory* TimePdfMaster::_r_run

protected

Definition at line 105 of file TimePdfMaster.h.

_r_S

RooRealVar* TimePdfMaster::_r_S

protected

Definition at line 159 of file TimePdfMaster.h.

_r_S_bar

RooRealVar* TimePdfMaster::_r_S_bar

protected

Definition at line 160 of file TimePdfMaster.h.

_r_scale_mean_dt

RooRealVar* TimePdfMaster::_r_scale_mean_dt

protected

Definition at line 112 of file TimePdfMaster.h.

_r_scale_sigma_2_dt

RooRealVar* TimePdfMaster::_r_scale_sigma_2_dt

protected

Definition at line 115 of file TimePdfMaster.h.

_r_scale_sigma_dt

RooRealVar* TimePdfMaster::_r_scale_sigma_dt

protected

Definition at line 114 of file TimePdfMaster.h.

_r_t

RooRealVar* TimePdfMaster::_r_t

protected

Cast of MINT parameters to B2DX parameters.

Definition at line 97 of file TimePdfMaster.h.

_r_tageff_asym_os

RooRealVar* TimePdfMaster::_r_tageff_asym_os

protected

Definition at line 135 of file TimePdfMaster.h.

_r_tageff_asym_ss

RooRealVar* TimePdfMaster::_r_tageff_asym_ss

protected

Definition at line 142 of file TimePdfMaster.h.

_r_tageff_os

RooRealVar* TimePdfMaster::_r_tageff_os

protected

Definition at line 134 of file TimePdfMaster.h.

_r_tageff_ss

RooRealVar* TimePdfMaster::_r_tageff_ss

protected

Definition at line 141 of file TimePdfMaster.h.

_r_tau

RooRealVar* TimePdfMaster::_r_tau

protected

Definition at line 108 of file TimePdfMaster.h.

_r_trigger

RooCategory* TimePdfMaster::_r_trigger

protected

Definition at line 106 of file TimePdfMaster.h.

_resmodel

RooGaussModel* TimePdfMaster::_resmodel

protected

Definition at line 203 of file TimePdfMaster.h.

_sampleGen

RooMCStudy* TimePdfMaster::_sampleGen

protected

Definition at line 208 of file TimePdfMaster.h.

_samplingPdf

RooProdPdf* TimePdfMaster::_samplingPdf

protected

Definition at line 205 of file TimePdfMaster.h.

_samplingPdf_t

RooBDecay* TimePdfMaster::_samplingPdf_t

protected

Definition at line 201 of file TimePdfMaster.h.

_samplingPdf_t_eff

RooEffProd* TimePdfMaster::_samplingPdf_t_eff

protected

Definition at line 204 of file TimePdfMaster.h.

_scale_mean_dt

const MINT::FitParameter& TimePdfMaster::_scale_mean_dt

protected

Definition at line 62 of file TimePdfMaster.h.

_scale_sigma_2_dt

const MINT::FitParameter& TimePdfMaster::_scale_sigma_2_dt

protected

Definition at line 64 of file TimePdfMaster.h.

_scale_sigma_dt

const MINT::FitParameter& TimePdfMaster::_scale_sigma_dt

protected

Definition at line 63 of file TimePdfMaster.h.

_sin_coeff

DecRateCoeff_Bd* TimePdfMaster::_sin_coeff

protected

Definition at line 165 of file TimePdfMaster.h.

_sin_term

TimePdfIntegrator* TimePdfMaster::_sin_term

protected

Definition at line 172 of file TimePdfMaster.h.

_sinh_coeff

DecRateCoeff_Bd* TimePdfMaster::_sinh_coeff

protected

Definition at line 166 of file TimePdfMaster.h.

_sinh_term

TimePdfIntegrator* TimePdfMaster::_sinh_term

protected

Definition at line 170 of file TimePdfMaster.h.

_spline

RooCubicSplineFun* TimePdfMaster::_spline

protected

Definition at line 148 of file TimePdfMaster.h.

_splinePdf

RooProduct* TimePdfMaster::_splinePdf

protected

Definition at line 149 of file TimePdfMaster.h.

_tageff_asym_os

const MINT::FitParameter& TimePdfMaster::_tageff_asym_os

protected

Definition at line 83 of file TimePdfMaster.h.

_tageff_asym_ss

const MINT::FitParameter& TimePdfMaster::_tageff_asym_ss

protected

Definition at line 91 of file TimePdfMaster.h.

_tageff_os

const MINT::FitParameter& TimePdfMaster::_tageff_os

protected

Definition at line 82 of file TimePdfMaster.h.

_tageff_ss

const MINT::FitParameter& TimePdfMaster::_tageff_ss

protected

Definition at line 90 of file TimePdfMaster.h.

---

The documentation for this class was generated from the following file:

• Mint/TimePdfMaster.h