TimeDependentGeneratorOld Class Reference

#include <TimeDependentGeneratorOld.h>

Classes
class	GenTimeEvent
class	GenTimePoint

Public Types
typedef std::list< GenTimePoint >	GenList
typedef std::map< int, GenList >	GenMap
typedef std::pair< std::complex< double >, std::complex< double > >	AmpPair

Public Member Functions
	TimeDependentGeneratorOld (const std::string &name, const bool overwrite, TRandom3 *rndm, double precision, const DalitzEventPattern &pattern, double width, double deltam, double deltagamma, double qoverp, double phi, double tmax, int ntimepoints, const bool saveIntegEvents=true, double tmin=0., TH1F *h_efficiency=NULL, float resWidth=0.05, bool addExpEffects=false)
AmpPair	amplitude_coefficients (const int tag, const double decaytime)
int	generate_tag () const
double	generate_decay_time (const int tag) const
MINT::counted_ptr< IDalitzEvent >	generate_dalitz_event (const int tag, const double decaytime) const
MINT::counted_ptr< IDalitzEvent >	generate_event () const
const std::map< int, SplineGenerator >	time_generators () const

Static Public Member Functions
static DalitzEventPattern	anti (DalitzEventPattern pat)

Private Attributes
const std::string	m_name
TRandom3 *	m_rndm
const DalitzEventPattern	m_pattern
const DalitzEventPattern	m_cppattern
const double	m_width
const double	m_deltam
const double	m_deltagamma
const double	m_qoverp
const double	m_phi
const double	m_tmax
const double	m_tmin
const int	m_ntimepoints
GenMap	m_genmap
std::map< int, SplineGenerator >	m_timegenerators
double	m_tagintegralfrac
double	m_precision
TH1F *	m_h_efficiency
TSpline3	m_efficiencyFit
float	m_resWidth
bool	m_addExpEffects

Detailed Description

Definition at line 17 of file TimeDependentGeneratorOld.h.

Member Typedef Documentation

AmpPair

typedef std::pair<std::complex<double>, std::complex<double> > TimeDependentGeneratorOld::AmpPair

Definition at line 34 of file TimeDependentGeneratorOld.h.

GenList

typedef std::list<GenTimePoint> TimeDependentGeneratorOld::GenList

Definition at line 32 of file TimeDependentGeneratorOld.h.

GenMap

typedef std::map<int, GenList> TimeDependentGeneratorOld::GenMap

Definition at line 33 of file TimeDependentGeneratorOld.h.

Constructor & Destructor Documentation

TimeDependentGeneratorOld()

TimeDependentGeneratorOld::TimeDependentGeneratorOld	(	const std::string &	name,
		const bool	overwrite,
		TRandom3 *	rndm,
		double	precision,
		const DalitzEventPattern &	pattern,
		double	width,
		double	deltam,
		double	deltagamma,
		double	qoverp,
		double	phi,
		double	tmax,
		int	ntimepoints,
		const bool	saveIntegEvents = true,
		double	tmin = 0.,
		TH1F *	h_efficiency = NULL,
		float	resWidth = 0.05,
		bool	addExpEffects = false
	)

Definition at line 98 of file TimeDependentGeneratorOld.cpp.

                                                                                   :
  m_name(name),
  m_rndm(rndm),
  m_pattern(pattern),
  m_cppattern(anti(pattern)),
  m_width(width),
  m_deltam(deltam),
  m_deltagamma(deltagamma),
  m_qoverp(qoverp),
  m_phi(phi),
  m_tmax(tmax),
  m_tmin(tmin),
  m_ntimepoints(ntimepoints),
  m_genmap(),
  m_timegenerators(),
  m_tagintegralfrac(0.),
  m_precision(precision),
  m_h_efficiency(h_efficiency),
  m_resWidth(resWidth),
  m_addExpEffects(addExpEffects),
  m_efficiencyFit()
{
  // If overwrite is true and the integrators directory exists, delete it.
  if(overwrite && exists(name)){
    cout << "Deleting previous integrators in directory " << name << endl ;
    string cmd("rm -rf " + name) ;
    system(cmd.c_str()) ;
  }

  if(!exists(name)){
    string cmd("mkdir -p " + name) ;
    system(cmd.c_str()) ;
  }
    
  const DalitzEventPattern* patterns[] = {&m_cppattern, &m_pattern} ;
  double sampleinterval = m_ntimepoints > 1 ? (m_tmax - m_tmin)/(m_ntimepoints-1) : 0. ;
  // Loop over flavours.
  for(int tag = -1 ; tag <= 1 ; tag += 2) {
    m_genmap[tag] = GenList() ;
    const DalitzEventPattern* evtpat = patterns[(tag+1)/2] ;
    const DalitzEventPattern* antipat = patterns[((tag+1)/2 + 1) % 2] ;
    // cout << "evtpat " ;
    // evtpat->print() ;
    // cout << " antipat " ;
    // antipat->print() ;
    // cout << endl ;
    vector<double> times ;
    vector<double> integrals ;
    // Loop over decay time sample points.
    for(int i = 0 ; i < m_ntimepoints ; ++i){
      double decaytime = m_tmin + i * sampleinterval ;
      AmpPair amps = amplitude_coefficients(tag, decaytime) ;
      FitAmpSum* model(new FitAmpSum(*evtpat)) ;
      *model *= amps.first ;
      if(amps.second != complex<double>(0., 0.)){
    FitAmpSum antimodel(*antipat) ;
    antimodel *= amps.second ;
    model->add(antimodel) ;
      }
      SignalGenerator* generator = new SignalGenerator(*evtpat, model) ;
      //cout << "Time point " << i << endl ;
      //model->print() ;
      //auto evt = generator->newEvent() ;
      //model->printAllAmps(*evt) ;
      ostringstream fname ;
      fname << m_name << "/tag_" << tag << "_decaytime_" << decaytime ;
      double integral(0.) ;
      // Calculate the integral if necessary.
      if(!exists(fname.str())){
    const string eventsFile(fname.str() + "_events.root") ;
    DalitzPdfSaveInteg dalitz(*evtpat, model, m_precision, fname.str(),
                  eventsFile, "topUp", fname.str()) ;
    integral = dalitz.getIntegralValue() ;
    dalitz.saveIntegrator(fname.str()) ;
    if(!saveIntegEvents){
      string cmd("rm " + eventsFile) ;
      system(cmd.c_str()) ;
    }
      }
      // Else retrive the integral from a file.
      else {
    auto intcalc = model->makeIntegrationCalculator() ;
    intcalc->retrieve(fname.str()) ;
    integral = intcalc->integral() ;
      }
      cout << "Make generator with tag " << tag << ", decay time " << decaytime
       << ", coeffprod " << amps.first.real()
       << " + " << amps.first.imag() << " j, "
       << " coeffmix " << amps.second.real() << " + " << amps.second.imag()
       << " j, integral " << integral << endl ;
      if(integral == 0.){
    complex<double> coeff = amps.first + amps.second ;
    integral = sqrt(coeff.real() * coeff.real() + coeff.imag() * coeff.imag()) ;
      }
      m_genmap[tag].push_back(GenTimePoint(decaytime, model, integral, generator)) ;
      times.push_back(decaytime) ;
      integrals.push_back(integral) ;
    }
    // Make a spline interpolator of the decay time distributions, to be used to generate
    // the decay times.
    ostringstream splinename ;
    splinename << "timespline_tag_" << tag ;
    string splinenamestr = splinename.str() ;
    TSpline3 timespline(splinenamestr.c_str(), &times[0], &integrals[0], times.size()) ;
    timespline.SetName(splinenamestr.c_str()) ;
    m_timegenerators.insert(make_pair(tag, SplineGenerator(rndm, timespline))) ;
  }
  // Calculate the integrated CP asymmetry.
  double integminus = m_timegenerators.find(-1)->second.integral() ;
  double integplus = m_timegenerators.find(1)->second.integral() ;
  cout << "Integrated CP asymmetry is " << (integplus - integminus)/(integplus + integminus) << endl ;
  m_tagintegralfrac = integminus/(integminus + integplus) ;
  double tauminus = m_timegenerators.find(-1)->second.mean() ;
  double tauplus = m_timegenerators.find(1)->second.mean() ;
  cout << "Tau minus: " << tauminus << endl ;
  cout << "Tau plus: " << tauplus << endl ;
  cout << "AGamma: " << (tauminus - tauplus)/(tauminus + tauplus) << endl ;
  
  if( m_h_efficiency != NULL){
    m_efficiencyFit = TSpline3(m_h_efficiency) ;
  }
}

Member Function Documentation

amplitude_coefficients()

TimeDependentGeneratorOld::AmpPair TimeDependentGeneratorOld::amplitude_coefficients	(	const int	tag,
		const double	decaytime
	)

Definition at line 230 of file TimeDependentGeneratorOld.cpp.

                                                                                       {
  double coeff = exp(-decaytime * 0.5 * (m_width + 0.5 * m_deltagamma)) ;
  complex<double> expterm = exp(complex<double>(0.5 * m_deltagamma * decaytime, m_deltam * decaytime)) ;
  complex<double> plusterm = 1. + expterm ;
  complex<double> minusterm = 1. - expterm ;
  complex<double> coeffprod = coeff * plusterm ;
  complex<double> coeffmix = pow(polar(m_qoverp, m_phi), tag) * coeff * minusterm ;
  return AmpPair(coeffprod, coeffmix) ;
}

anti()

DalitzEventPattern TimeDependentGeneratorOld::anti ( DalitzEventPattern  pat )

static

Definition at line 78 of file TimeDependentGeneratorOld.cpp.

                                                                         {
  pat[0].antiThis() ;
  return pat ;
}

generate_dalitz_event()

MINT::counted_ptr< IDalitzEvent > TimeDependentGeneratorOld::generate_dalitz_event	(	const int	tag,
		const double	decaytime
	)		const

Definition at line 257 of file TimeDependentGeneratorOld.cpp.

                                                                                                                          {
  const GenList& genlist = m_genmap.find(tag)->second ;
  GenList::const_iterator igen = genlist.begin() ;
  while(igen->decaytime < decaytime && igen != genlist.end())
    ++igen ;
  if(igen == genlist.end()){
    cerr << "TimeDependentGeneratorOld::generate_dalitz_event: ERROR: Got impossible decay time: "
     << decaytime << " (tmax = " << m_tmax << ")" << endl ;
    return MINT::counted_ptr<IDalitzEvent>(0) ;
  }
  // Unlikely, but best to check.
  if(decaytime == igen->decaytime)
    return igen->generator->newEvent() ;
  GenList::const_iterator igenprev(igen) ;
  --igen ;
  // Pick between the generators either side of the decay time according to how close they are to it.
  double gensel = m_rndm->Rndm() ;
  if(gensel < 1. - (decaytime - igenprev->decaytime)/(igen->decaytime - igenprev->decaytime))
    return igenprev->generator->newEvent() ;
  return igen->generator->newEvent() ;
}

generate_decay_time()

double TimeDependentGeneratorOld::generate_decay_time

(

const int

tag

)

const

Definition at line 249 of file TimeDependentGeneratorOld.cpp.

                                                                         {
  double decaytime = m_tmax + 1. ;
  while(decaytime > m_tmax)
    decaytime = m_timegenerators.find(tag)->second.gen_random() ;
  return decaytime ; 
}

generate_event()

MINT::counted_ptr< IDalitzEvent > TimeDependentGeneratorOld::generate_event

(

)

const

Definition at line 280 of file TimeDependentGeneratorOld.cpp.

                                                                              {
  int tag = generate_tag() ;
  double decaytime = generate_decay_time(tag) ;
  double smeareddecaytime = -999. ;

  smeareddecaytime = decaytime + m_rndm->Gaus(0, m_resWidth) ;

  if ( (m_h_efficiency != NULL) && (m_addExpEffects) ){
    int i = 0 ;
    float efficiency = 0 ;
    int maxiter = 100000 ;
    while(true){
      if( smeareddecaytime > m_efficiencyFit.GetXmax() ){
        efficiency = m_efficiencyFit.Eval( m_efficiencyFit.GetXmax() ) ;
      } 
      else if( smeareddecaytime > m_efficiencyFit.GetXmin() ){
    efficiency = m_efficiencyFit.Eval( smeareddecaytime ) ;
      }
      else{
        efficiency = m_efficiencyFit.Eval( m_efficiencyFit.GetXmin() ) ;
      }      

      if(m_rndm->Rndm() < efficiency){
        break;
      }
 
      tag = generate_tag() ;
      decaytime = generate_decay_time(tag) ;
      smeareddecaytime =  decaytime + m_rndm->Gaus(0, m_resWidth);

      i += 1 ;
      if( i > maxiter ){ 
        cout << "WARNING: Decay time generation limit exceeded." << endl ; 
        break ;
      }
    }
  }
  MINT::counted_ptr<IDalitzEvent> evt = generate_dalitz_event(tag, decaytime) ;
  return MINT::counted_ptr<IDalitzEvent>(new GenTimeEvent(*evt, tag, decaytime, smeareddecaytime)) ;
}

generate_tag()

int TimeDependentGeneratorOld::generate_tag

(

)

const

Definition at line 241 of file TimeDependentGeneratorOld.cpp.

                                                  {
  double rndm = m_rndm->Rndm() ;
  if(rndm < m_tagintegralfrac)
    return -1 ;
  return 1 ;
}

time_generators()

const map< int, SplineGenerator > TimeDependentGeneratorOld::time_generators

(

)

const

Definition at line 322 of file TimeDependentGeneratorOld.cpp.

                                                                                 {
  return m_timegenerators;
}

Member Data Documentation

m_addExpEffects

bool TimeDependentGeneratorOld::m_addExpEffects

private

Definition at line 126 of file TimeDependentGeneratorOld.h.

m_cppattern

const DalitzEventPattern TimeDependentGeneratorOld::m_cppattern

private

Definition at line 103 of file TimeDependentGeneratorOld.h.

m_deltagamma

const double TimeDependentGeneratorOld::m_deltagamma

private

Definition at line 107 of file TimeDependentGeneratorOld.h.

m_deltam

const double TimeDependentGeneratorOld::m_deltam

private

Definition at line 106 of file TimeDependentGeneratorOld.h.

m_efficiencyFit

TSpline3 TimeDependentGeneratorOld::m_efficiencyFit

private

Definition at line 123 of file TimeDependentGeneratorOld.h.

m_genmap

GenMap TimeDependentGeneratorOld::m_genmap

private

Definition at line 115 of file TimeDependentGeneratorOld.h.

m_h_efficiency

TH1F* TimeDependentGeneratorOld::m_h_efficiency

private

Definition at line 122 of file TimeDependentGeneratorOld.h.

m_name

const std::string TimeDependentGeneratorOld::m_name

private

Definition at line 100 of file TimeDependentGeneratorOld.h.

m_ntimepoints

const int TimeDependentGeneratorOld::m_ntimepoints

private

Definition at line 113 of file TimeDependentGeneratorOld.h.

m_pattern

const DalitzEventPattern TimeDependentGeneratorOld::m_pattern

private

Definition at line 102 of file TimeDependentGeneratorOld.h.

m_phi

const double TimeDependentGeneratorOld::m_phi

private

Definition at line 109 of file TimeDependentGeneratorOld.h.

m_precision

double TimeDependentGeneratorOld::m_precision

private

Definition at line 120 of file TimeDependentGeneratorOld.h.

m_qoverp

const double TimeDependentGeneratorOld::m_qoverp

private

Definition at line 108 of file TimeDependentGeneratorOld.h.

m_resWidth

float TimeDependentGeneratorOld::m_resWidth

private

Definition at line 125 of file TimeDependentGeneratorOld.h.

m_rndm

TRandom3* TimeDependentGeneratorOld::m_rndm

private

Definition at line 101 of file TimeDependentGeneratorOld.h.

m_tagintegralfrac

double TimeDependentGeneratorOld::m_tagintegralfrac

private

Definition at line 119 of file TimeDependentGeneratorOld.h.

m_timegenerators

std::map<int, SplineGenerator> TimeDependentGeneratorOld::m_timegenerators

private

Definition at line 117 of file TimeDependentGeneratorOld.h.

m_tmax

const double TimeDependentGeneratorOld::m_tmax

private

Definition at line 111 of file TimeDependentGeneratorOld.h.

m_tmin

const double TimeDependentGeneratorOld::m_tmin

private

Definition at line 112 of file TimeDependentGeneratorOld.h.

m_width

const double TimeDependentGeneratorOld::m_width

private

Definition at line 105 of file TimeDependentGeneratorOld.h.

---

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

• Mint/TimeDependentGeneratorOld.h
• src/Utils/TimeDependentGeneratorOld.cpp