FullTimePdf Class Reference

#include <FullTimePdf.h>

Inheritance diagram for FullTimePdf:

Public Member Functions
void	parametersChanged ()
void	beginFit ()
void	endFit ()
double	un_normalised (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent &evt)
double	getNorm (IDalitzEvent &evt)
double	getSampledPdfVal (IDalitzEvent &evt)
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 (double &eta_SS)
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)
double	getCalibratedResolution (double dt)
double	getCalibratedResolution (double &dt, double &scale_sigma_dt, double &scale_sigma_2_dt)
virtual double	getVal_withPs (IDalitzEvent &evt)
virtual double	getVal_noPs (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent *evt)
virtual double	getVal_withPs (IDalitzEvent *evt)
virtual double	getVal_noPs (IDalitzEvent *evt)
RooDataSet *	sampleEvents (int N=10000)
void	generateBkgToys (int N, DalitzEventList &eventListData)
DalitzEventList	generateToysRooFit (int N=10000, int run=-1, int trigger=-1)
DalitzEventList	generateToys (int N=10000, int run=-1, int trigger=-1)
void	saveEventListToFile (DalitzEventList &eventList, string name="toys.root")
DalitzEvent	generateWeightedEvent ()
TH1D *	plotSpline ()
	FullTimePdf (const MINT::FitParameter &C, const MINT::FitParameter &D, const MINT::FitParameter &D_bar, const MINT::FitParameter &S, const MINT::FitParameter &S_bar, const MINT::FitParameter &k, 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="")
Public Member Functions inherited from MINT::PdfBase< IDalitzEvent >
	PdfBase ()
	PdfBase (const PdfBase< IDalitzEvent > &)
virtual double	RealVal (IDalitzEvent &evt)
virtual double	getNewVal (IDalitzEvent &evt)
virtual	~PdfBase ()
Public Member Functions inherited from MINT::IPdf< IDalitzEvent >
virtual void	Gradient (IDalitzEvent &evt, std::vector< double > &grad, MINT::MinuitParameterSet *mps)
virtual bool	useAnalyticGradient ()
Public Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
virtual	~IReturnRealForEvent ()

Protected Attributes
const FitParameter &	_C
	fit parameters More...
const FitParameter &	_D
const FitParameter &	_D_bar
const FitParameter &	_S
const FitParameter &	_S_bar
const FitParameter &	_k
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
TimePdfMaster *	_timePdfMaster
DalitzEventPattern	_pat
NamedParameter< double >	_min_TAU
NamedParameter< double >	_max_TAU
NamedParameter< double >	_min_TAUERR
NamedParameter< double >	_max_TAUERR

Additional Inherited Members
Protected Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
	IReturnRealForEvent ()

Detailed Description

Definition at line 79 of file FullTimePdf.h.

Constructor & Destructor Documentation

FullTimePdf()

FullTimePdf::FullTimePdf ( const MINT::FitParameter &  C, const MINT::FitParameter &  D, const MINT::FitParameter &  D_bar, const MINT::FitParameter &  S, const MINT::FitParameter &  S_bar, const MINT::FitParameter &  k, 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

Definition at line 554 of file FullTimePdf.h.

                 :
    _C(C),
    _D(D),
    _D_bar(D_bar),
    _S(S),
    _S_bar(S_bar),
    _k(k),
    _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)
    {
        _timePdfMaster = new TimePdfMaster(_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,
                                          _p0_os, _p1_os, _delta_p0_os, _delta_p1_os, 
                                          _avg_eta_os, _tageff_os, _tageff_asym_os, 
                                          _p0_ss, p1_ss, _delta_p0_ss, _delta_p1_ss, 
                                          _avg_eta_ss, _tageff_ss, _tageff_asym_ss, 
                                          _production_asym, _detection_asym,_marginalPdfsPrefix);

    _timePdfMaster->setAllFitParameters();

        _timePdfMaster->setCP_coeff(1., 1.,
                                   _C,-_C,
                                   _k * _D, _k * _D_bar,
                                   _k * _S, _k * _S_bar  );  

        NamedParameter<int> EventPattern("Event Pattern", 521, 321, 211, -211, 443);
        _pat = DalitzEventPattern(EventPattern.getVector());
    }

Member Function Documentation

beginFit()

void FullTimePdf::beginFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 144 of file FullTimePdf.h.

                   {
        _timePdfMaster->listFitParDependencies();
    }

endFit()

void FullTimePdf::endFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 147 of file FullTimePdf.h.

                 {
    }

generateBkgToys()

void FullTimePdf::generateBkgToys ( int  N, DalitzEventList &  eventListData  )

inline

Definition at line 276 of file FullTimePdf.h.

                                                               {

    DalitzEventList eventList,eventListDataSideband;

    for(int i=0; i< eventListData.size(); i++)
    {   
        DalitzEvent evt = eventListData[i];
        if(abs(evt.getValueFromVector(9) -5370) > 60)eventListDataSideband.Add(evt);
    }
    int N_sample = eventListDataSideband.size();

    vector<int> b_indices;
    while( b_indices.size() < N )b_indices.push_back(TMath::Nint(gRandom->Uniform(0,N_sample)));
    sort(b_indices.begin(), b_indices.end());
    
    TRandom3 rndm;
    for(int i=0; i< N; i++)
    {   
        DalitzEvent evt = eventListDataSideband[b_indices[i]];
        eventList.Add(evt);
    }
    saveEventListToFile(eventList);
    }

generateToys()

DalitzEventList FullTimePdf::generateToys ( int  N = 10000, int  run = -1, int  trigger = -1  )

inline

Estimate max val

Safety check on the maxmimum generated height

Hit-and-miss

Definition at line 336 of file FullTimePdf.h.

                                                                                {

    time_t startTime = time(0);

    cout << "Generating " << N << " events" << endl;
    DalitzEventList eventList;

    vector<double> vals;    
    for(int i = 0; i < 100000; i++){
        DalitzEvent evt = generateWeightedEvent();
        double val = getVal(evt)/evt.getGeneratorPdfRelativeToPhaseSpace(); //_timePdfMaster->get_marginalPdfs_product(evt);
        vals.push_back(val);    
    }

    cout << "Now calculating maximum val " << vals.size() << endl;
    double amax,pmax;
    generalisedPareto_estimateMaximum(vals,0.999,amax,pmax);
    
    double pdf_max = 1.;
    if(!TMath::IsNaN(pmax) && pmax > 0 && pmax < 100 * amax)pdf_max = pmax;
    else if(!TMath::IsNaN(amax))pdf_max = amax;
    // for safety
    pdf_max *= 1.5; 

    cout << "pdf_max " << pdf_max << endl;

    int N_gen = 0;
    int N_tot = 0;
    while(true){
            DalitzEvent evt = generateWeightedEvent();
            double pdfVal = getVal(evt)/evt.getGeneratorPdfRelativeToPhaseSpace(); // /_timePdfMaster->get_marginalPdfs_product(evt);
            
            const double height = gRandom->Uniform(0,pdf_max);
            
            if( pdfVal > pdf_max ){
                std::cout << "ERROR: PDF above determined maximum." << std::endl;
                std::cout << pdfVal << " > " << pdf_max << std::endl;
                pdf_max = pdf_max * 2.;
            }
            
            if( height < pdfVal ) { 
                evt.setValueInVector(7, run);
                evt.setValueInVector(8, trigger);
                eventList.Add(evt);
                N_gen++;
                //if (0ul == (N_gen % 500ul)) cout << "Generated event " << N_gen << "/" << N << endl;
            }       
            N_tot ++;
            if(N_gen == N)break;
    }

    cout << " Done. " << " Total time since start " << (time(0) - startTime)/60.0 << " min." << endl;
    cout << "Generated " << N_gen << " events ! Efficiecy = " << (double)N_gen/(double)N_tot << endl;

//  saveEventListToFile(eventList);
    return eventList;
   }

generateToysRooFit()

DalitzEventList FullTimePdf::generateToysRooFit ( int  N = 10000, int  run = -1, int  trigger = -1  )

inline

Definition at line 300 of file FullTimePdf.h.

                                                                                      {

    cout << "Generating " << N << " events" << endl;
    DalitzEventList eventList;

    RooDataSet* sample = sampleEvents(N);
    
    for(int i = 0; i < sample->numEntries(); i++){
    
        RooArgSet* sample_set= (RooArgSet*)sample->get(i);
    
        double t_MC = ((RooRealVar*)sample_set->find("t"))->getVal() ;
        double dt_MC = ((RooRealVar*)sample_set->find("dt"))->getVal() ;
        double eta_OS_MC = ((RooRealVar*)sample_set->find("eta_OS"))->getVal() ;
        double eta_SS_MC = ((RooRealVar*)sample_set->find("eta_SS"))->getVal() ;

        int f_MC = ((RooCategory*)sample_set->find("qf"))->getIndex() ;
        int q_OS_MC = ((RooCategory*)sample_set->find("q_OS"))->getIndex() ;
        int q_SS_MC = ((RooCategory*)sample_set->find("q_SS"))->getIndex() ;    
    
        DalitzEvent evt(_pat,gRandom);
        evt.setValueInVector(0, t_MC);
        evt.setValueInVector(1, dt_MC);
        evt.setValueInVector(2, f_MC);
        evt.setValueInVector(3, q_OS_MC);
        evt.setValueInVector(4, eta_OS_MC);
        evt.setValueInVector(5, q_SS_MC);
        evt.setValueInVector(6, eta_SS_MC);
        evt.setValueInVector(7, run);
        evt.setValueInVector(8, trigger);

        eventList.Add(evt);
    }
    return eventList;
   }

generateWeightedEvent()

DalitzEvent FullTimePdf::generateWeightedEvent ( )

inline

Definition at line 506 of file FullTimePdf.h.

                                       {

    while(true){
        double t_MC = gRandom->Exp(1./_Gamma);
//      double t_MC = gRandom->Exp(1./_Gamma)+_min_TAU;
                if(t_MC > _max_TAU || t_MC < _min_TAU)continue;

        DalitzEvent evt(_pat,gRandom);

        vector<double> marginal_vals = _timePdfMaster->getRandom_marginalVals();
        double dt_MC = marginal_vals[0] ;
        double eta_OS_MC = marginal_vals[1] ;
        double eta_SS_MC = marginal_vals[2] ;

        int f_MC = (gRandom->Uniform() > 0.5) ? 1 : -1;     
    
        // true flavor
        int q_MC = (gRandom->Uniform() > 0.5) ? 1 : -1;

            int q_SS_MC = (gRandom->Uniform() > 2./3.) ? 0 : q_MC ;
            int q_OS_MC = (gRandom->Uniform() > 2./3.) ? 0 : q_MC ;
 
        q_OS_MC = (gRandom->Uniform() < 0.5) ? - q_OS_MC : q_OS_MC;
        q_SS_MC = (gRandom->Uniform() < 0.5) ? - q_SS_MC : q_SS_MC;
        
        eta_OS_MC = (q_OS_MC == 0) ? 0.5 : eta_OS_MC;
        eta_SS_MC = (q_SS_MC == 0) ? 0.5 : eta_SS_MC;

        if(f_MC<0)evt.CP_conjugateYourself();
        evt.setValueInVector(0, t_MC);
        evt.setValueInVector(1, dt_MC);
        evt.setValueInVector(2, f_MC);
        evt.setValueInVector(3, q_OS_MC);
        evt.setValueInVector(4, eta_OS_MC);
        evt.setValueInVector(5, q_SS_MC);
        evt.setValueInVector(6, eta_SS_MC);

        evt.setGeneratorPdfRelativeToPhaseSpace(_Gamma * exp(-t_MC*_Gamma) / ( ( exp(-_min_TAU*_Gamma) - exp(-_max_TAU * _Gamma) )));
//      evt.setGeneratorPdfRelativeToPhaseSpace(_Gamma * ( exp(-t_MC*_Gamma) + _min_TAU )  / ( ( 1. - exp(-_max_TAU * _Gamma) ) + _Gamma * _min_TAU * ( _max_TAU - 0.)));

        return evt;
    }
    }

getCalibratedMistag()

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

inline

Definition at line 244 of file FullTimePdf.h.

                                                                                                                               {
        return _timePdfMaster->getCalibratedMistag(eta, avg_eta, p0, p1, delta_p0, delta_p1);
    }

getCalibratedMistag_OS() [1/3]

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

inline

Definition at line 220 of file FullTimePdf.h.

                                                                     {
        return _timePdfMaster->getCalibratedMistag_OS(evt);
    }

getCalibratedMistag_OS() [2/3]

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

inline

Definition at line 224 of file FullTimePdf.h.

                                                                  {
        return _timePdfMaster->getCalibratedMistag_OS(eta_OS);
    }

getCalibratedMistag_OS() [3/3]

std::pair<double, double> FullTimePdf::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 228 of file FullTimePdf.h.

                                                                                                                                                             {
        return _timePdfMaster->getCalibratedMistag_OS(evt, avg_eta_os, p0_os, p1_os, delta_p0_os, delta_p1_os);
    }

getCalibratedMistag_SS() [1/3]

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

inline

Definition at line 232 of file FullTimePdf.h.

                                                                     {
        return _timePdfMaster->getCalibratedMistag_SS(evt);
    }

getCalibratedMistag_SS() [2/3]

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

inline

Definition at line 236 of file FullTimePdf.h.

                                                                  {
        return _timePdfMaster->getCalibratedMistag_SS(eta_SS);
    }

getCalibratedMistag_SS() [3/3]

std::pair<double, double> FullTimePdf::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 240 of file FullTimePdf.h.

                                                                                                                                                             {
        return _timePdfMaster->getCalibratedMistag_SS(evt, avg_eta_ss, p0_ss, p1_ss, delta_p0_ss, delta_p1_ss);
    }

getCalibratedResolution() [1/2]

double FullTimePdf::getCalibratedResolution ( double  dt )

inline

Definition at line 248 of file FullTimePdf.h.

                                             {
        return _timePdfMaster->getCalibratedResolution(dt);
    }

getCalibratedResolution() [2/2]

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

inline

Definition at line 252 of file FullTimePdf.h.

                                                                                              {
        return _timePdfMaster->getCalibratedResolution(dt,scale_sigma_dt,scale_sigma_2_dt);
    }

getNorm()

double FullTimePdf::getNorm ( IDalitzEvent &  evt )

inline

Definition at line 192 of file FullTimePdf.h.

                                     {

        _timePdfMaster->setAllObservablesAndFitParameters(evt);
        
        // C,Cbar,D,Dbar,S,Sbar
        _timePdfMaster->setCP_coeff(1.,
                   1.,
                                   _C,
                                   -_C,
                                   _k * _D,
                                   _k * _D_bar,
                                   _k * _S,
                                   _k * _S_bar
                                   );

        double norm =
        _timePdfMaster->get_cosh_term_Integral(evt)
        +  _timePdfMaster->get_cos_term_Integral(evt)
        +  _timePdfMaster->get_sinh_term_Integral(evt)
        +  _timePdfMaster->get_sin_term_Integral(evt);
        
        return norm;
    }

getSampledPdfVal()

double FullTimePdf::getSampledPdfVal ( IDalitzEvent &  evt )

inline

Definition at line 216 of file FullTimePdf.h.

                                                  {
    return _timePdfMaster->getSamplingPdfVal(evt);
    }

getVal() [1/2]

virtual double FullTimePdf::getVal ( IDalitzEvent &  evt )

inlinevirtual

Implements MINT::PdfBase< IDalitzEvent >.

Definition at line 179 of file FullTimePdf.h.

                                            {
        
        double val = un_normalised(evt);
        
        double norm =
        _timePdfMaster->get_cosh_term_Integral(evt)
        +  _timePdfMaster->get_cos_term_Integral(evt)
        +  _timePdfMaster->get_sinh_term_Integral(evt)
        +  _timePdfMaster->get_sin_term_Integral(evt);
        
        return val/norm;
    }

getVal() [2/2]

virtual double FullTimePdf::getVal ( IDalitzEvent *  evt )

inlinevirtual

Definition at line 259 of file FullTimePdf.h.

                                            {
        if(0 == evt) return 0;
        return getVal(*evt);
    }

getVal_noPs() [1/2]

virtual double FullTimePdf::getVal_noPs ( IDalitzEvent &  evt )

inlinevirtual

Definition at line 257 of file FullTimePdf.h.

{return getVal(evt);}

getVal_noPs() [2/2]

virtual double FullTimePdf::getVal_noPs ( IDalitzEvent *  evt )

inlinevirtual

Definition at line 267 of file FullTimePdf.h.

                                                 {
        if(0 == evt) return 0;
        return getVal_noPs(*evt);
    }

getVal_withPs() [1/2]

virtual double FullTimePdf::getVal_withPs ( IDalitzEvent &  evt )

inlinevirtual

Definition at line 256 of file FullTimePdf.h.

{return getVal(evt);}

getVal_withPs() [2/2]

virtual double FullTimePdf::getVal_withPs ( IDalitzEvent *  evt )

inlinevirtual

Definition at line 263 of file FullTimePdf.h.

                                                   {
        if(0 == evt) return 0;
        return getVal_withPs(*evt);
    }

parametersChanged()

void FullTimePdf::parametersChanged ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 142 of file FullTimePdf.h.

                            {
    }

plotSpline()

TH1D* FullTimePdf::plotSpline ( )

inline

Definition at line 550 of file FullTimePdf.h.

                       {
     return _timePdfMaster->plotSpline();
    }

sampleEvents()

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

inline

Definition at line 272 of file FullTimePdf.h.

                                           {
    return _timePdfMaster->sampleEvents(N);
    }

saveEventListToFile()

void FullTimePdf::saveEventListToFile ( DalitzEventList &  eventList, string  name = "toys.root"  )

inline

Definition at line 397 of file FullTimePdf.h.

                                                                                   {

        TFile* out = new TFile(name.c_str(),"RECREATE");
        TTree* tree = new TTree("DecayTree","DecayTree");
            
        double t,dt;
        int Ds_ID;
        int q_OS,q,q_SS;
        double eta_OS;
        double eta_SS;
        double sw,w;
        int run,trigger;
        int year;
        
        double K[4];
        double pip[4];
        double pim[4];
        double Ds[4];
        double mB;

            TBranch* br_mB = tree->Branch( "Bs_DTF_MM", &mB, "Bs_DTF_MM/D" );
            TBranch* br_sw = tree->Branch( "N_Bs_sw", &sw, "N_Bs_sw/D" );
            TBranch* br_w = tree->Branch( "weight", &w, "weight/D" );

            TBranch* br_t = tree->Branch( "Bs_BsDTF_TAU", &t, "Bs_BsDTF_TAU/D" );
            TBranch* br_dt = tree->Branch( "Bs_BsDTF_TAUERR", &dt, "Bs_BsDTF_TAUERR/D" );

        TBranch* br_Ds_ID = tree->Branch("Ds_ID",&Ds_ID,"Ds_ID/I");

        TBranch* br_q_OS =tree->Branch("OS_Combination_DEC",&q_OS,"OS_Combination_DEC/I");
        TBranch* br_eta_OS =  tree->Branch("OS_Combination_PROB",&eta_OS,"OS_Combination_PROB/D");
        TBranch* br_q_SS = tree->Branch("SS_Kaon_DEC",&q_SS,"SS_Kaon_DEC/I");
        TBranch* br_eta_SS = tree->Branch("SS_Kaon_PROB",&eta_SS,"SS_Kaon_PROB/D");

        TBranch* br_run = tree->Branch("run",&run,"run/I");
            TBranch* br_year = tree->Branch( "year", &year, "year/I" );
        TBranch* br_trigger = tree->Branch("TriggerCat",&trigger,"TriggerCat/I");
        
        TBranch* br_K0 = tree->Branch("BsDTF_Kplus_PX",&K[0],"BsDTF_Kplus_PX/D");
        TBranch* br_K1 =tree->Branch("BsDTF_Kplus_PY",&K[1],"BsDTF_Kplus_PY/D");
        TBranch* br_K2 =tree->Branch("BsDTF_Kplus_PZ",&K[2],"BsDTF_Kplus_PZ/D");
        TBranch* br_K3 =tree->Branch("BsDTF_Kplus_PE",&K[3],"BsDTF_Kplus_PE/D");
    
            TBranch* br_pip0 =tree->Branch("BsDTF_piplus_PX",&pip[0],"BsDTF_piplus_PX/D");
        TBranch* br_pip1 =tree->Branch("BsDTF_piplus_PY",&pip[1],"BsDTF_piplus_PY/D");
        TBranch* br_pip2 =tree->Branch("BsDTF_piplus_PZ",&pip[2],"BsDTF_piplus_PZ/D");
        TBranch* br_pip3 =tree->Branch("BsDTF_piplus_PE",&pip[3],"BsDTF_piplus_PE/D");
    
        TBranch* br_pim0 =tree->Branch("BsDTF_piminus_PX",&pim[0],"BsDTF_piminus_PX/D");
        TBranch* br_pim1 =tree->Branch("BsDTF_piminus_PY",&pim[1],"BsDTF_piminus_PY/D");
        TBranch* br_pim2 =tree->Branch("BsDTF_piminus_PZ",&pim[2],"BsDTF_piminus_PZ/D");
        TBranch* br_pim3 =tree->Branch("BsDTF_piminus_PE",&pim[3],"BsDTF_piminus_PE/D");
    
        TBranch* br_Ds0 =tree->Branch("BsDTF_Ds_PX",&Ds[0],"BsDTF_Ds_PX/D");
        TBranch* br_Ds1 =tree->Branch("BsDTF_Ds_PY",&Ds[1],"BsDTF_Ds_PY/D");
        TBranch* br_Ds2 =tree->Branch("BsDTF_Ds_PZ",&Ds[2],"BsDTF_Ds_PZ/D");
        TBranch* br_Ds3 =tree->Branch("BsDTF_Ds_PE",&Ds[3],"BsDTF_Ds_PE/D");

        for(int i= 0; i < eventList.size(); i++){

        t = eventList[i].getValueFromVector(0);
        dt = eventList[i].getValueFromVector(1);

        Ds_ID = - eventList[i].getValueFromVector(2);
        
        q_OS = eventList[i].getValueFromVector(3);
        eta_OS = eventList[i].getValueFromVector(4);

        q_SS = eventList[i].getValueFromVector(5);
        eta_SS = eventList[i].getValueFromVector(6);

        run = eventList[i].getValueFromVector(7);
        if(run == 2) year = 16;
        else if(run == 3) year = 17;
        else year = 12;
        trigger = eventList[i].getValueFromVector(8);

        mB = eventList[i].p(0).M(); //eventList[i].getValueFromVector(9);
        sw = eventList[i].getWeight();
        w = eventList[i].getWeight();

        Ds[0] = eventList[i].p(1).Px()/MeV;
        Ds[1] = eventList[i].p(1).Py()/MeV;
        Ds[2] = eventList[i].p(1).Pz()/MeV;
        Ds[3] = eventList[i].p(1).E()/MeV;

        K[0] = eventList[i].p(2).Px()/MeV;
        K[1] = eventList[i].p(2).Py()/MeV;
        K[2] = eventList[i].p(2).Pz()/MeV;
        K[3] = eventList[i].p(2).E()/MeV;

        pip[0] = eventList[i].p(3).Px()/MeV;
        pip[1] = eventList[i].p(3).Py()/MeV;
        pip[2] = eventList[i].p(3).Pz()/MeV;
        pip[3] = eventList[i].p(3).E()/MeV;

        pim[0] = eventList[i].p(4).Px()/MeV;
        pim[1] = eventList[i].p(4).Py()/MeV;
        pim[2] = eventList[i].p(4).Pz()/MeV;
        pim[3] = eventList[i].p(4).E()/MeV;

        tree->Fill();
         }

         tree->Write();
         out->Write();
         out->Close();
    }

un_normalised()

double FullTimePdf::un_normalised ( IDalitzEvent &  evt )

inline

Definition at line 150 of file FullTimePdf.h.

                                                  {
        
        //const double t = (double) evt.getValueFromVector(0);
        //if(t < _min_TAU || t > _max_TAU )return 0.;

        _timePdfMaster->setAllObservablesAndFitParameters(evt);
        
        // C,Cbar,D,Dbar,S,Sbar
        _timePdfMaster->setCP_coeff(1.,
                   1.,
                                   _C,
                                   -_C,
                                   _k * _D,
                                   _k * _D_bar,
                                   _k * _S,
                                   _k * _S_bar
                                   );
        
        double val =
        (
         _timePdfMaster->get_cosh_term_Val(evt)
         +  _timePdfMaster->get_cos_term_Val(evt)
         +  _timePdfMaster->get_sinh_term_Val(evt)
         +  _timePdfMaster->get_sin_term_Val(evt)
         );// * _timePdfMaster->get_marginalPdfs_product(evt);
        
        return val;
    }

Member Data Documentation

_avg_eta_os

const MINT::FitParameter& FullTimePdf::_avg_eta_os

protected

Definition at line 114 of file FullTimePdf.h.

_avg_eta_ss

const MINT::FitParameter& FullTimePdf::_avg_eta_ss

protected

Definition at line 122 of file FullTimePdf.h.

_C

const FitParameter& FullTimePdf::_C

protected

fit parameters

Definition at line 83 of file FullTimePdf.h.

_c0

const MINT::FitParameter& FullTimePdf::_c0

protected

Definition at line 99 of file FullTimePdf.h.

_c1

const MINT::FitParameter& FullTimePdf::_c1

protected

Definition at line 100 of file FullTimePdf.h.

_c2

const MINT::FitParameter& FullTimePdf::_c2

protected

Definition at line 101 of file FullTimePdf.h.

_c3

const MINT::FitParameter& FullTimePdf::_c3

protected

Definition at line 102 of file FullTimePdf.h.

_c4

const MINT::FitParameter& FullTimePdf::_c4

protected

Definition at line 103 of file FullTimePdf.h.

_c5

const MINT::FitParameter& FullTimePdf::_c5

protected

Definition at line 104 of file FullTimePdf.h.

_c6

const MINT::FitParameter& FullTimePdf::_c6

protected

Definition at line 105 of file FullTimePdf.h.

_c7

const MINT::FitParameter& FullTimePdf::_c7

protected

Definition at line 106 of file FullTimePdf.h.

_c8

const MINT::FitParameter& FullTimePdf::_c8

protected

Definition at line 107 of file FullTimePdf.h.

_c9

const MINT::FitParameter& FullTimePdf::_c9

protected

Definition at line 108 of file FullTimePdf.h.

_D

const FitParameter& FullTimePdf::_D

protected

Definition at line 84 of file FullTimePdf.h.

_D_bar

const FitParameter& FullTimePdf::_D_bar

protected

Definition at line 85 of file FullTimePdf.h.

_delta_p0_os

const MINT::FitParameter& FullTimePdf::_delta_p0_os

protected

Definition at line 112 of file FullTimePdf.h.

_delta_p0_ss

const MINT::FitParameter& FullTimePdf::_delta_p0_ss

protected

Definition at line 120 of file FullTimePdf.h.

_delta_p1_os

const MINT::FitParameter& FullTimePdf::_delta_p1_os

protected

Definition at line 113 of file FullTimePdf.h.

_delta_p1_ss

const MINT::FitParameter& FullTimePdf::_delta_p1_ss

protected

Definition at line 121 of file FullTimePdf.h.

_detection_asym

const MINT::FitParameter& FullTimePdf::_detection_asym

protected

Definition at line 127 of file FullTimePdf.h.

_dGamma

const MINT::FitParameter& FullTimePdf::_dGamma

protected

Definition at line 91 of file FullTimePdf.h.

_dm

const MINT::FitParameter& FullTimePdf::_dm

protected

Definition at line 92 of file FullTimePdf.h.

_Gamma

const MINT::FitParameter& FullTimePdf::_Gamma

protected

Definition at line 90 of file FullTimePdf.h.

_k

const FitParameter& FullTimePdf::_k

protected

Definition at line 88 of file FullTimePdf.h.

_marginalPdfsPrefix

string FullTimePdf::_marginalPdfsPrefix

protected

Definition at line 129 of file FullTimePdf.h.

_max_TAU

NamedParameter<double> FullTimePdf::_max_TAU

protected

Definition at line 137 of file FullTimePdf.h.

_max_TAUERR

NamedParameter<double> FullTimePdf::_max_TAUERR

protected

Definition at line 139 of file FullTimePdf.h.

_min_TAU

NamedParameter<double> FullTimePdf::_min_TAU

protected

Definition at line 136 of file FullTimePdf.h.

_min_TAUERR

NamedParameter<double> FullTimePdf::_min_TAUERR

protected

Definition at line 138 of file FullTimePdf.h.

_offset_sigma_dt

const MINT::FitParameter& FullTimePdf::_offset_sigma_dt

protected

Definition at line 94 of file FullTimePdf.h.

_p0_os

const MINT::FitParameter& FullTimePdf::_p0_os

protected

Definition at line 110 of file FullTimePdf.h.

_p0_ss

const MINT::FitParameter& FullTimePdf::_p0_ss

protected

Definition at line 118 of file FullTimePdf.h.

_p1_os

const MINT::FitParameter& FullTimePdf::_p1_os

protected

Definition at line 111 of file FullTimePdf.h.

_p1_ss

const MINT::FitParameter& FullTimePdf::_p1_ss

protected

Definition at line 119 of file FullTimePdf.h.

_pat

DalitzEventPattern FullTimePdf::_pat

protected

Definition at line 133 of file FullTimePdf.h.

_production_asym

const MINT::FitParameter& FullTimePdf::_production_asym

protected

Definition at line 126 of file FullTimePdf.h.

_S

const FitParameter& FullTimePdf::_S

protected

Definition at line 86 of file FullTimePdf.h.

_S_bar

const FitParameter& FullTimePdf::_S_bar

protected

Definition at line 87 of file FullTimePdf.h.

_scale_mean_dt

const MINT::FitParameter& FullTimePdf::_scale_mean_dt

protected

Definition at line 95 of file FullTimePdf.h.

_scale_sigma_2_dt

const MINT::FitParameter& FullTimePdf::_scale_sigma_2_dt

protected

Definition at line 97 of file FullTimePdf.h.

_scale_sigma_dt

const MINT::FitParameter& FullTimePdf::_scale_sigma_dt

protected

Definition at line 96 of file FullTimePdf.h.

_tageff_asym_os

const MINT::FitParameter& FullTimePdf::_tageff_asym_os

protected

Definition at line 116 of file FullTimePdf.h.

_tageff_asym_ss

const MINT::FitParameter& FullTimePdf::_tageff_asym_ss

protected

Definition at line 124 of file FullTimePdf.h.

_tageff_os

const MINT::FitParameter& FullTimePdf::_tageff_os

protected

Definition at line 115 of file FullTimePdf.h.

_tageff_ss

const MINT::FitParameter& FullTimePdf::_tageff_ss

protected

Definition at line 123 of file FullTimePdf.h.

_timePdfMaster

TimePdfMaster* FullTimePdf::_timePdfMaster

protected

Definition at line 132 of file FullTimePdf.h.

---

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

• Mint/FullTimePdf.h