MappedDalitzBWArea Class Reference

#include <MappedDalitzBWArea.h>

Public Member Functions
	MappedDalitzBWArea ()
	MappedDalitzBWArea (const DalitzEventPattern &pat, TRandom *rnd=gRandom)
	MappedDalitzBWArea (const DalitzEventPattern &pat, const MINT::counted_ptr< IGenFct > &oneFunction, TRandom *rnd=gRandom)
	MappedDalitzBWArea (const DalitzEventPattern &pat, const MINT::counted_ptr< IGenFct > &oneFunction, const MINT::counted_ptr< IGenFct > &otherFunction, TRandom *rnd=gRandom)
	MappedDalitzBWArea (const DalitzEventPattern &pat, const std::vector< MINT::counted_ptr< IGenFct > > &functions, TRandom *rnd=gRandom)
	MappedDalitzBWArea (const MappedDalitzBWArea &other)
virtual	~MappedDalitzBWArea ()
bool	checkIntegration () const
MINT::counted_ptr< DalitzEvent >	tryEventForOwner () const
bool	isInside (const DalitzEvent &evt) const
bool	isInside (const DalitzCoordinate &dc) const
bool	isInside (const std::vector< DalitzCoordinate > &dcList) const
double	genValue (const DalitzEvent &evt) const
double	size () const
double	integral () const
void	setUnWeightPs (bool doSo=true)
bool	unWeightPs ()
void	print (std::ostream &os=std::cout) const
MINT::counted_ptr< IGenFct >	fct_for_s (const std::vector< int > &indices) const
MINT::counted_ptr< IGenFct >	fct_for_s (int i, int j, int k=-1) const
MappedDalitzBWArea &	operator= (const MappedDalitzBWArea &other)
bool	setRnd (TRandom *rnd=gRandom)

Private Member Functions
void	setup ()
void	findMapping ()
void	applyLimits ()
bool	allConsecutive () const
bool	allStandardised () const
std::vector< TLorentzVector > &	mapP4 (const DalitzEvent &evt, const Permutation &mapping, std::vector< TLorentzVector > &p4) const
std::vector< TLorentzVector >	mapP4 (const DalitzEvent &evt, const Permutation &mapping) const

Private Attributes
TRandom *	_rnd
DalitzEventPattern	_pat
DalitzEventPattern	_mappedPat
DalitzBWArea	_area
std::vector< MINT::counted_ptr< IGenFct > >	_limits
Permutation	_mapping
Permutation	_inverseMapping

Detailed Description

Definition at line 21 of file MappedDalitzBWArea.h.

Constructor & Destructor Documentation

MappedDalitzBWArea() [1/6]

MappedDalitzBWArea::MappedDalitzBWArea

(

)

Definition at line 24 of file MappedDalitzBWArea.cpp.

  : _rnd(0)
  , _pat()
  , _mappedPat()
  , _area()
  , _limits()
  , _mapping()
  , _inverseMapping()
{
  setup();
}

MappedDalitzBWArea() [2/6]

MappedDalitzBWArea::MappedDalitzBWArea	(	const DalitzEventPattern &	pat,
		TRandom *	rnd = gRandom
	)

Definition at line 35 of file MappedDalitzBWArea.cpp.

  : _rnd(rnd)
  , _pat(pat)
  , _mappedPat(pat)
  , _area(pat)
  , _limits()
  , _mapping(pat.size())
  , _inverseMapping(pat.size())
{
  setup();
}

MappedDalitzBWArea() [3/6]

MappedDalitzBWArea::MappedDalitzBWArea	(	const DalitzEventPattern &	pat,
		const MINT::counted_ptr< IGenFct > &	oneFunction,
		TRandom *	rnd = gRandom
	)

Definition at line 47 of file MappedDalitzBWArea.cpp.

  : _rnd(rnd)
  , _pat(pat)
  , _mappedPat(pat)
  , _limits(1)
  , _mapping(pat.size())
  , _inverseMapping(pat.size())
{
  _limits[0] = oneLimit;
  setup();
}

MappedDalitzBWArea() [4/6]

MappedDalitzBWArea::MappedDalitzBWArea	(	const DalitzEventPattern &	pat,
		const MINT::counted_ptr< IGenFct > &	oneFunction,
		const MINT::counted_ptr< IGenFct > &	otherFunction,
		TRandom *	rnd = gRandom
	)

Definition at line 61 of file MappedDalitzBWArea.cpp.

  : _rnd(rnd)
  , _pat(pat)
  , _mappedPat(pat)
  , _limits(2)
  , _mapping(pat.size())
  , _inverseMapping(pat.size())
{
  _limits[0] = oneLimit;
  _limits[1] = otherLimit;
  /*
  cout << " just made a MappedDalitzBWArea with limits :"
       << oneLimit << ", " << otherLimit << " or "
       << _limits[0] << ", " << _limits[1] << endl;
  */
  setup();
  /*
  cout << " ..after setup: "
       << _limits[0] << ", " << _limits[1] << endl;
  */
}

MappedDalitzBWArea() [5/6]

MappedDalitzBWArea::MappedDalitzBWArea	(	const DalitzEventPattern &	pat,
		const std::vector< MINT::counted_ptr< IGenFct > > &	functions,
		TRandom *	rnd = gRandom
	)

Definition at line 87 of file MappedDalitzBWArea.cpp.

  : _rnd(rnd)
  , _pat(pat)
  , _mappedPat(pat)
  , _limits(lim)
  , _mapping(pat.size())
  , _inverseMapping(pat.size())
{
  setup();
}

MappedDalitzBWArea() [6/6]

MappedDalitzBWArea::MappedDalitzBWArea

(

const MappedDalitzBWArea &

other

)

Definition at line 101 of file MappedDalitzBWArea.cpp.

  : _rnd(other._rnd)
  , _pat(other._pat)
  , _mappedPat(other._mappedPat)
  , _area(other._area)
  , _limits(other._limits)
  , _mapping(other._mapping)
  , _inverseMapping(other._inverseMapping)
{
}

~MappedDalitzBWArea()

MappedDalitzBWArea::~MappedDalitzBWArea ( )

virtual

Definition at line 123 of file MappedDalitzBWArea.cpp.

                                       {
}

Member Function Documentation

allConsecutive()

bool MappedDalitzBWArea::allConsecutive ( ) const

private

Definition at line 126 of file MappedDalitzBWArea.cpp.

                                             {
  for(unsigned int i=0; i < _limits.size(); i++){
    DalitzCoordinate newCoordinates = (_limits[i])->getCoordinate().mapMe(_mapping);
    /*
    cout << "mapped new coordinates: " 
     << i << ": " << newCoordinates << endl;
    */
    if(! (newCoordinates.I_am_Consecutive())) return false;
  }
  return true;
}

allStandardised()

bool MappedDalitzBWArea::allStandardised ( ) const

private

Definition at line 138 of file MappedDalitzBWArea.cpp.

                                              {
  // such that only s123, s12,
  // or             s12, s34
  // are filled (or s123, or s12 for single 
  // limits).
  // Helps with generating later.

  if(_limits.size() > 2 || _limits.empty() ){
    return true; // no standards for that, all OK
  }
  DalitzCoordinate s123(1,2,3);
  DalitzCoordinate s12(1,2);
  DalitzCoordinate s34(3,4);
  if(_limits.size() == 1){
    if(_limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s123)) return true;
    if(_limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s12)) return true;
  }
  if(_limits.size() == 2){
    if(_limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s123)
       && _limits[1]->getCoordinate().mapMe(_mapping).sameIndices(s12)) return true;
    if(_limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s12)
       && _limits[1]->getCoordinate().mapMe(_mapping).sameIndices(s34)) return true;
    if(_limits[1]->getCoordinate().mapMe(_mapping).sameIndices(s123)
       && _limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s12)) return true;
    if(_limits[1]->getCoordinate().mapMe(_mapping).sameIndices(s12)
       && _limits[0]->getCoordinate().mapMe(_mapping).sameIndices(s34)) return true;
  }
  return false;
}

applyLimits()

void MappedDalitzBWArea::applyLimits ( )

private

Definition at line 235 of file MappedDalitzBWArea.cpp.

                                    {
  bool dbThis=false;
  for(unsigned int i=0; i < _limits.size(); i++){
    if(dbThis){
      cout << "mapped co-ordinates " 
       << (_limits[i])->getCoordinate().mapMe(_mapping)
       << endl;
    }
    _area.setFcn((_limits[i])->getCoordinate().mapMe(_mapping), _limits[i]);
  }
}

checkIntegration()

bool MappedDalitzBWArea::checkIntegration ( ) const

inline

Definition at line 94 of file MappedDalitzBWArea.h.

                               {
    return _area.checkIntegration();}

fct_for_s() [1/2]

MINT::counted_ptr<IGenFct> MappedDalitzBWArea::fct_for_s

(

const std::vector< int > &

indices

)

const

fct_for_s() [2/2]

MINT::counted_ptr<IGenFct> MappedDalitzBWArea::fct_for_s	(	int	i,
		int	j,
		int	k = -1
	)		const

findMapping()

void MappedDalitzBWArea::findMapping ( )

private

Definition at line 168 of file MappedDalitzBWArea.cpp.

                                    {
  bool dbThis=false;
  _mapping.makeUnity();
  if(allConsecutive() && allStandardised()) return;

  DalitzEventPattern tmpPat(_pat);
  tmpPat[0] =1;
  for(unsigned int i=1; i< _pat.size(); i++){
    tmpPat[i]=0;
  }

  Permutator ptor(tmpPat);// this finds all mappings
  // leaving the 0th particle alone (the mother) and
  // swapping all the daughters - the pattern
  // is set up to achieve this ('particles' 1,2,3,4,.. are identical)

  bool foundOne=false;
  for(unsigned int i=0; i< ptor.size(); i++){
    _mapping = ptor[i];
    if(dbThis) cout << " trying out mapping\n " << _mapping << endl;
    if(allConsecutive() && allStandardised()){
      foundOne=true;
      break;
    }
  }
  if(! foundOne){
    cout << " MappedDalitzBWArea::findMapping: WARNING!!"
     << " didn't find a mapping!!!"
     << endl;
    throw "shit!";
  }
  _inverseMapping = _mapping.getInverse();
  for(unsigned int i=0; i < _pat.size(); i++){
    int mappedIndex = _mapping[i];
    if(mappedIndex < 0 || mappedIndex >= (int) _mappedPat.size()){
      cout << "ERROR in MappedDalitzBWArea::findMapping()"
       << "\n   Index out of range: " << mappedIndex
       << " not in [0, " << _mappedPat.size()-1
       << endl;
      throw "index out of range.";
    }
    if(dbThis) {
      cout << " mapped index of " << i 
       << " = " << mappedIndex << endl;
    }
    _mappedPat[mappedIndex] = _pat[i];
  }

  if(dbThis){
    cout << " using the following mapping:\n "
     << _mapping << endl;
    cout << "_mappedPat " << _mappedPat << endl;
  }
  return;
}

genValue()

double MappedDalitzBWArea::genValue

(

const DalitzEvent &

evt

)

const

Definition at line 307 of file MappedDalitzBWArea.cpp.

                                                               {
  bool dbThis=false;

  double returnVal = _area.genValue(&evt, _inverseMapping);
  if(dbThis){
    vector<TLorentzVector> p4 =  mapP4(evt, _mapping );
    DalitzEvent mappedEvt(_mappedPat, p4);
    isInside(evt);
    cout << " mapping: " << _mapping << endl;    
    cout << "old style: " << _area.genValue(&mappedEvt)
     << "\n new style: " << _area.genValue(&evt, _inverseMapping)
     << "\n MappedDalitzBWArea::genValue returning " << returnVal << endl;
  }
  return returnVal;
}

integral()

double MappedDalitzBWArea::integral ( ) const

inline

Definition at line 107 of file MappedDalitzBWArea.h.

                         {
    return _area.integral();
  }

isInside() [1/3]

bool MappedDalitzBWArea::isInside

(

const DalitzEvent &

evt

)

const

Definition at line 323 of file MappedDalitzBWArea.cpp.

                                                             {
  bool dbThis=false;

  if(dbThis){
    vector<TLorentzVector> p4 = 
      mapP4(evt, _mapping );
    DalitzEvent mappedEvt(_mappedPat, p4);
    cout << " old style: " <<  _area.isInside(mappedEvt);
    cout << " new style: " << _area.isInside(evt, _inverseMapping);
  }
  return _area.isInside(evt, _inverseMapping);
}

isInside() [2/3]

bool MappedDalitzBWArea::isInside

(

const DalitzCoordinate &

dc

)

const

Definition at line 336 of file MappedDalitzBWArea.cpp.

                                                                 {
  return _area.isInside(dc.mapMe(_mapping));
}

isInside() [3/3]

bool MappedDalitzBWArea::isInside

(

const std::vector< DalitzCoordinate > &

dcList

)

const

Definition at line 339 of file MappedDalitzBWArea.cpp.

                                                                                {
  //  cout << " inside this: " << *this << endl;
  for(unsigned int i=0; i < dcList.size(); i++){
    if(! isInside(dcList[i])){
      //     cout << "Not inside: " << dcList[i]
      //       << "\n\tmapped: " << dcList[i].mapMe(_mapping)
      //       << endl;
      return false;
    }
  }
  return true;

}

mapP4() [1/2]

std::vector<TLorentzVector>& MappedDalitzBWArea::mapP4 ( const DalitzEvent &  evt, const Permutation &  mapping, std::vector< TLorentzVector > &  p4  ) const

inlineprivate

Definition at line 49 of file MappedDalitzBWArea.h.

                            {
    unsigned int n = evt.eventPattern().size();
    p4.resize(n); 
    for(unsigned int i=0; i < n; i++){
      int mappedIndex = mapping[i];
      if(mappedIndex < 0 || mappedIndex + 1 > (int) n){
    std::cout << "ERROR in MappedDalitzBWArea::mapP4()"
          << "\n   Index out of range: " << mappedIndex
          << " not in [0, " <<   n-1
          << std::endl;
    throw "index out of range.";
      }
      p4[mappedIndex] = evt.p(i);
    }
    return p4;
  }

mapP4() [2/2]

std::vector<TLorentzVector> MappedDalitzBWArea::mapP4 ( const DalitzEvent &  evt, const Permutation &  mapping  ) const

inlineprivate

Definition at line 68 of file MappedDalitzBWArea.h.

                          {
    std::vector<TLorentzVector> p4(evt.eventPattern().size());
    return mapP4(evt, mapping, p4);
  }

operator=()

MappedDalitzBWArea & MappedDalitzBWArea::operator=

(

const MappedDalitzBWArea &

other

)

Definition at line 112 of file MappedDalitzBWArea.cpp.

                                                                                {
  _rnd = other._rnd;
  _pat = other._pat;
  _mappedPat = other._mappedPat;
  _area = other._area;
  _limits = other._limits;
  _mapping = other._mapping;
  _inverseMapping = other._inverseMapping;
  return *this;
}

print()

void MappedDalitzBWArea::print

(

std::ostream &

os = std::cout

)

const

Definition at line 353 of file MappedDalitzBWArea.cpp.

                                                  {
  os << " MappedDalitzBWArea "
     << " limits: ";
  for(unsigned int i=0; i < _limits.size(); i++){
    if(0 == _limits[i]){
      os << "\n  ZERO-PTR LIMIT!" << endl;
    }else{
      os << "\n  " << (_limits[i])->getCoordinate() << endl;
    }
  }

  os << "\n  Mapping:\n " << _mapping
     << "\n area: " << _area;
}

setRnd()

bool MappedDalitzBWArea::setRnd

(

TRandom *

rnd = gRandom

)

Definition at line 368 of file MappedDalitzBWArea.cpp.

                                           {
  _rnd = rnd;
  _area.setRnd(_rnd);
  return true;
}

setUnWeightPs()

void MappedDalitzBWArea::setUnWeightPs ( bool  doSo = true )

inline

Definition at line 110 of file MappedDalitzBWArea.h.

{_area.setUnWeightPs(doSo);}

setup()

void MappedDalitzBWArea::setup ( )

private

Definition at line 224 of file MappedDalitzBWArea.cpp.

                              {
  bool dbThis=false;
  sort(_limits.begin(), _limits.end());

  findMapping();
  if(dbThis) cout << " mappedPat " << _mappedPat << endl;
  _area.setPattern(_mappedPat);
  _area.setRnd(_rnd);
  applyLimits();
}

size()

double MappedDalitzBWArea::size ( ) const

inline

Definition at line 105 of file MappedDalitzBWArea.h.

                     {
    return _area.size();}

tryEventForOwner()

counted_ptr< DalitzEvent > MappedDalitzBWArea::tryEventForOwner

(

)

const

Definition at line 278 of file MappedDalitzBWArea.cpp.

                                                                   {
  bool dbThis=false;
  if(dbThis){ 
    cout << "Hello from MappedDalitzBWArea::tryEventForOwner()" 
     << endl;
  }
  counted_ptr<DalitzEvent> evt( _area.tryEventForOwner(_inverseMapping));

  /*
  if(dbThis && 0 != evt){
    counted_ptr<DalitzEvent> mappedEvt( _area.tryEventForOwner());
    if(0 == mappedEvt){
      return mappedEvt;
    }
    static vector<TLorentzVector> p4;
    mapP4( *mappedEvt, _inverseMapping, p4);
    
    counted_ptr<DalitzEvent> evtCheck(new DalitzEvent(_pat, p4));
    evtCheck->setWeight(mappedEvt->getWeight());

    cout << "compare new, old event and weight: "
     << "\n\t new: " << *evt
     << "\n\t old: " << *evtCheck
     << endl;
  }
  */
  return evt;
}

unWeightPs()

bool MappedDalitzBWArea::unWeightPs ( )

inline

Definition at line 111 of file MappedDalitzBWArea.h.

{return _area.unWeightPs();}

Member Data Documentation

_area

DalitzBWArea MappedDalitzBWArea::_area

private

Definition at line 26 of file MappedDalitzBWArea.h.

_inverseMapping

Permutation MappedDalitzBWArea::_inverseMapping

private

Definition at line 30 of file MappedDalitzBWArea.h.

_limits

std::vector<MINT::counted_ptr<IGenFct> > MappedDalitzBWArea::_limits

private

Definition at line 28 of file MappedDalitzBWArea.h.

_mappedPat

DalitzEventPattern MappedDalitzBWArea::_mappedPat

private

Definition at line 24 of file MappedDalitzBWArea.h.

_mapping

Permutation MappedDalitzBWArea::_mapping

private

Definition at line 30 of file MappedDalitzBWArea.h.

_pat

DalitzEventPattern MappedDalitzBWArea::_pat

private

Definition at line 24 of file MappedDalitzBWArea.h.

_rnd

TRandom* MappedDalitzBWArea::_rnd

private

Definition at line 22 of file MappedDalitzBWArea.h.

---

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

• Mint/MappedDalitzBWArea.h
• src/Mojito/BreitWignerMC/MappedDalitzBWArea.cpp