Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtPropSLPole.cpp
| Show First 20 Lines • Show All 194 Lines • ▼ Show 20 Lines | double EvtPropSLPole::calBreitWigner( EvtParticle* pmeson, EvtPoint1D point ) | ||||
| } | } | ||||
| _massMax = _mass + maxdelta; | _massMax = _mass + maxdelta; | ||||
| if ( _massMin < 0. ) | if ( _massMin < 0. ) | ||||
| _massMin = 0.; | _massMin = 0.; | ||||
| EvtParticle* par = pmeson->getParent(); | EvtParticle* par = pmeson->getParent(); | ||||
| double maxMass = -1.; | double maxMass = -1.; | ||||
| if ( par != 0 ) { | if ( par != nullptr ) { | ||||
| if ( par->hasValidP4() ) | if ( par->hasValidP4() ) | ||||
| maxMass = par->mass(); | maxMass = par->mass(); | ||||
| for ( size_t i = 0; i < par->getNDaug(); i++ ) { | for ( size_t i = 0; i < par->getNDaug(); i++ ) { | ||||
| EvtParticle* tDaug = par->getDaug( i ); | EvtParticle* tDaug = par->getDaug( i ); | ||||
| if ( pmeson != tDaug ) | if ( pmeson != tDaug ) | ||||
| maxMass -= EvtPDL::getMinMass( tDaug->getId() ); | maxMass -= EvtPDL::getMinMass( tDaug->getId() ); | ||||
| } | } | ||||
| } | } | ||||
| EvtId* dauId = 0; | EvtId* dauId = nullptr; | ||||
| double* dauMasses = 0; | double* dauMasses = nullptr; | ||||
| size_t nDaug = pmeson->getNDaug(); | size_t nDaug = pmeson->getNDaug(); | ||||
| if ( nDaug > 0 ) { | if ( nDaug > 0 ) { | ||||
| dauId = new EvtId[nDaug]; | dauId = new EvtId[nDaug]; | ||||
| dauMasses = new double[nDaug]; | dauMasses = new double[nDaug]; | ||||
| for ( size_t j = 0; j < nDaug; j++ ) { | for ( size_t j = 0; j < nDaug; j++ ) { | ||||
| dauId[j] = pmeson->getDaug( j )->getId(); | dauId[j] = pmeson->getDaug( j )->getId(); | ||||
| dauMasses[j] = pmeson->getDaug( j )->mass(); | dauMasses[j] = pmeson->getDaug( j )->mass(); | ||||
| } | } | ||||
| } | } | ||||
| EvtId* parId = 0; | EvtId* parId = nullptr; | ||||
| EvtId* othDaugId = 0; | EvtId* othDaugId = nullptr; | ||||
| EvtParticle* tempPar = pmeson->getParent(); | EvtParticle* tempPar = pmeson->getParent(); | ||||
| if ( tempPar ) { | if ( tempPar ) { | ||||
| parId = new EvtId( tempPar->getId() ); | parId = new EvtId( tempPar->getId() ); | ||||
| if ( tempPar->getNDaug() == 2 ) { | if ( tempPar->getNDaug() == 2 ) { | ||||
| if ( tempPar->getDaug( 0 ) == pmeson ) | if ( tempPar->getDaug( 0 ) == pmeson ) | ||||
| othDaugId = new EvtId( tempPar->getDaug( 1 )->getId() ); | othDaugId = new EvtId( tempPar->getDaug( 1 )->getId() ); | ||||
| else | else | ||||
| othDaugId = new EvtId( tempPar->getDaug( 0 )->getId() ); | othDaugId = new EvtId( tempPar->getDaug( 0 )->getId() ); | ||||
| ▲ Show 20 Lines • Show All 326 Lines • Show Last 20 Lines | |||||