ForwardEvolver.cc
No OneTemporary
Actions

Size

12 KB

Referenced Files

None

Subscribers

None

ForwardEvolver.cc
View Options

	// -- C++ --
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the ForwardEvolver class.
	//

	#include "ForwardEvolver.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"
	// #include "ThePEG/Interface/Parameter.h"
	#include "ThePEG/Interface/Reference.h"
	#include "Herwig++/Utilities/HwDebug.h"
	#include "ThePEG/Repository/EventGenerator.h"
	#include "ShowerParticle.h"
	#include "ShowerKinematics.h"
	#include "QtildaShowerKinematics1to2.h"
	#include "SplittingFunction.h"
	//#include "ShowerColourLine.h"
	#include "ThePEG/EventRecord/ColourLine.h"
	#include "ShowerConfig.h"
	#include "QQbarG.h"

	using namespace Herwig;
	using namespace ThePEG;

	ForwardEvolver::~ForwardEvolver() {}


	void ForwardEvolver::persistentOutput(PersistentOStream & os) const {
	os << _splittingGenerator;
	}


	void ForwardEvolver::persistentInput(PersistentIStream & is, int) {
	is >> _splittingGenerator;
	}


	ClassDescription<ForwardEvolver> ForwardEvolver::initForwardEvolver;
	// Definition of the static class description member.


	void ForwardEvolver::Init() {

	static ClassDocumentation<ForwardEvolver> documentation
	("This class is responsible for the forward showering of time-like particles",
	"It does also the special forward showering for decaying particles",
	"in which the angular ordering is reversed");

	static Reference<ForwardEvolver,SplittingGenerator>
	interfaceSplitGen("SplittingGenerator",
	"A reference to the SplittingGenerator object",
	&Herwig::ForwardEvolver::_splittingGenerator,
	false, false, true, false);

	}

	//-----------------------------------------------------------------------------

	bool ForwardEvolver::timeLikeShower(tPartCollHdlPtr ch,
	const tShowerVarsPtr showerVariables,
	//const tMECorrectionPtr meCorrectionPtr,
	tShowerParticlePtr particle,
	ShowerParticleVector & collecShoPar,
	const bool specialDecay )
	throw (Veto, Stop, Exception) {

	if ( HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << endl
	<< "ForwardEvolver::timeLikeShower() Evt #"
	<< generator()->currentEventNumber()
	<< "\t full __________________________"
	<< endl;
	}

	if ( HERWIG_DEBUG_LEVEL == HwDebug::minimal_Shower
	&& generator()->currentEventNumber() < 1000) {
	generator()->log() << "# event no " << generator()->currentEventNumber() << endl;
	}

	bool hasEmitted = false;
	bool vetoed = true;
	Energy ptMax;

	while (vetoed) { // ACHTUNG, just wrote this loop...

	if(HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower) {
	generator()->log() << "ForwardEvolver: attempt... _______________________________________________" << std::endl;
	}
	vetoed = false;
	ptMax = 0.0*GeV;
	if(abs(particle->id()) < 7) {
	if(particle->id() > 0) ptMax = showerVariables->largestPtQ();
	else ptMax = showerVariables->largestPtQbar();
	}
	particle->undecay();
	tShowerParticleVector particlesYetToShower;
	particlesYetToShower.push_back(particle);
	Branching fb;

	do {
	tShowerParticlePtr part = particlesYetToShower.back();
	particlesYetToShower.pop_back();

	if(HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower) {
	generator()->log() << "-- Yet " << particlesYetToShower.size() + 1
	<< ". next " << part->data().PDGName()
	// << " [" << part->number() << "]"
	<< endl;
	}

	if(HERWIG_DEBUG_LEVEL == HwDebug::minimal_Shower
	&& generator()->currentEventNumber() < 1000) {
	generator()->log() << "-<" << particlesYetToShower.size() + 1
	<< " " << part->data().PDGName();
	}
	//*LOOKHERE* update rhoD matrix of part ;

	fb=_splittingGenerator->chooseForwardBranching(ch, *part, specialDecay);

	if(HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower) {
	if(fb.first && fb.second ) {
	generator()->log() << " branching (int, Q_i -> Q_f) = ("
	<< fb.second->splittingFn()->interactionType()
	<< " "
	<< part->
	evolutionScales()[fb.second->splittingFn()->interactionType()]
	<< " -> "
	<< fb.first->qtilde()
	<< "). " << endl;
	} else {
	generator()->log() << " no branching." << endl;
	}
	}

	// accept it according to the showerVariables and soft correction;
	// soft correction:

	if(showerVariables->softMEC() && fb.first && fb.second &&
	fb.second->splittingFn()->interactionType() == ShowerIndex::QCD &&
	abs(part->id()) < 7) {

	double d_z = fb.first->z();
	Energy d_qt = fb.first->qtilde();
	Energy2 d_m2 = part->momentum().m2();
	Energy pPerp = (1.-d_z)sqrt( sqr(d_zd_qt) - d_m2);
	if ( part->id() == particle->id()
	&& HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << " Soft correction of "
	<< part->id()
	<< ": pT = "
	<< pPerp/GeV
	<< (pPerp > ptMax ? " > " : " < ")
	<< " pTmax = "
	<< ptMax/GeV
	<< " GeV"
	<< std::endl;
	}
	if(part->id() == particle->id() && pPerp > ptMax) {
	ptMax = pPerp;
	if ( HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << " Soft correction: is hardest so far."
	<< std::endl;
	}
	vetoed = MEVeto(particle, d_qt, d_z);
	}
	}
	if(vetoed) break;

	if(fb.first == ShoKinPtr() \|\| fb.second == tSudakovPtr()) {

	if ( HERWIG_DEBUG_LEVEL >= HwDebug::extreme_Shower ) {
	generator()->log() << "-- no further splitting."
	<< endl;
	}
	if ( HERWIG_DEBUG_LEVEL == HwDebug::minimal_Shower
	&& generator()->currentEventNumber() < 1000) {
	generator()->log() << endl;
	}
	//*LOOKHERE* rhoD propagation;
	} else {
	hasEmitted = true;
	_splittingGenerator->generateBranchingKinematics(ch, *part, fb);

	// Assign the splitting function and the shower kinematics
	// to the emitting particle.
	part->setShowerKinematics(fb.first);
	part->setSplittingFn(fb.second->splittingFn());

	// For the time being we are considering only 1->2 branching
	tSplittingFnPtr splitF = fb.second->splittingFn();
	if(splitF) {
	// Create the ShowerParticle objects for the two children of
	// the emitting particle; set the parent/child relationship;
	// add them to the collecShoPar and particlesYetToShower
	// collections. Remember that we our splitting functions
	// are associated only with particles (id>0) and not with
	// antiparticles (id<0)
	// (because we are assuming CP-conserving vertices):
	// therefore the signs of the decays products must be set by
	// hand explicitly; to simplify this, we assume that the
	// first product must always have the sign as the parent:
	// q->q+g, qbar->qbar+g, g->q+qbar, etc.
	// Notice that the momenta of the shower products is not
	// set: only at the end of the showering, during the
	// kinematics reconstruction such momenta are calculated and
	// set.

	ShowerParticlePtr showerProduct1;
	ShowerParticlePtr showerProduct2;
	if(part->data().id() > 0)
	showerProduct1 = new_ptr(ShowerParticle(getParticleData(fb.third[1])));
	else
	showerProduct1 = new_ptr(ShowerParticle(getParticleData(-fb.third[1])));
	showerProduct2 = new_ptr(ShowerParticle(getParticleData(fb.third[2])));

	const ShowerIndex::InteractionType interaction = splitF->interactionType();
	const Energy scale = part->showerKinematics()->qtilde();
	double zz = dynamic_ptr_cast<Ptr<QtildaShowerKinematics1to2>::
	transient_pointer>(part->showerKinematics())->z();
	// note that 1st child gets z, 2nd gets (1-z) by our convention.
	showerProduct1->setEvolutionScale(interaction, zz*scale);
	showerProduct2->setEvolutionScale(interaction, (1.-zz)*scale);
	showerProduct1->setInitiatesTLS(false);
	showerProduct2->setInitiatesTLS(false);

	ParticleVector theChildren;
	theChildren.push_back(showerProduct1);
	theChildren.push_back(showerProduct2);
	part->showerKinematics()->updateChildren(part, theChildren);

	// In the case of splittings which involves coloured particles,
	// set properly the colour flow of the branching.
	// Notice that the methods: ShowerColourLine::addColoured and
	// ShowerColourLine::addAntiColoured automatically set also,
	// respectively, the colourLine and antiColourLine of the
	// ShowerParticle object they received as argument.
	ShoColinePair parentShoColinePair =
	ShoColinePair(part->colourLine(), part->antiColourLine());
	ShoColinePair showerProduct1ShoColinePair = ShoColinePair();
	ShoColinePair showerProduct2ShoColinePair = ShoColinePair();
	splitF->colourConnection(parentShoColinePair,
	showerProduct1ShoColinePair,
	showerProduct2ShoColinePair);
	if ( showerProduct1ShoColinePair.first ) {
	showerProduct1ShoColinePair.first->addColoured( showerProduct1 );
	}
	if ( showerProduct1ShoColinePair.second ) {
	showerProduct1ShoColinePair.second->
	addAntiColoured( showerProduct1 );
	}
	if ( showerProduct2ShoColinePair.first ) {
	showerProduct2ShoColinePair.first->addColoured( showerProduct2 );
	}
	if ( showerProduct2ShoColinePair.second ) {
	showerProduct2ShoColinePair.second->
	addAntiColoured( showerProduct2 );
	}

	part->addChild(showerProduct1);
	part->addChild(showerProduct2);
	collecShoPar.insert(collecShoPar.end(), showerProduct1);
	collecShoPar.insert(collecShoPar.end(), showerProduct2);
	particlesYetToShower.push_back(showerProduct1);
	particlesYetToShower.push_back(showerProduct2);

	if ( HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << " Splitting: " << part->data().PDGName()
	<< " -> " << showerProduct1->data().PDGName()
	<< " + " << showerProduct2->data().PDGName()
	<< endl
	<< " \t \t \t \t colourLine \t antiColourLine " << endl;
	// Create a map of pointers to ShowerColourLine object,
	// in order to have a nice printing of colour lines,
	// numbered from 601 onwards, instead of printing
	// physical memory addresses.
	map<tShoColinePtr,int> mapShoColine;
	int countShoColine = 601;
	for ( int i=0; i<3; i++ ) {
	tShowerParticlePtr theParticle;
	if ( i == 0 ) {
	theParticle = part;
	generator()->log() << "\t parent: ";
	} else if ( i == 1 ) {
	theParticle = showerProduct1;
	generator()->log() << "\t product1: ";
	} else if ( i == 2 ) {
	theParticle = showerProduct2;
	generator()->log() << "\t product1: ";
	}
	generator()->log() << theParticle->data().PDGName() << "\t \t";
	if(theParticle->colourLine()) {
	if(mapShoColine.find(theParticle->colourLine()) == mapShoColine.end()) {
	mapShoColine[theParticle->colourLine()] = countShoColine++;
	}
	generator()->log() << mapShoColine[theParticle->colourLine()];
	} else {
	generator()->log() << "0";
	}
	generator()->log() << "\t \t";
	if ( theParticle->antiColourLine() ) {
	if ( mapShoColine.find( theParticle->antiColourLine() )
	== mapShoColine.end() ) {
	mapShoColine.insert(pair<tShoColinePtr,int>(theParticle->antiColourLine(),countShoColine++));
	}
	generator()->log() << mapShoColine.find(theParticle->
	antiColourLine())->second;
	} else {
	generator()->log() << "0";
	}
	generator()->log() << endl;
	} // for
	} // debug full

	if ( HERWIG_DEBUG_LEVEL == HwDebug::minimal_Shower
	&& generator()->currentEventNumber() < 1000) {
	generator()->log() << "->" << showerProduct1->data().PDGName()
	<< "+" << showerProduct2->data().PDGName()
	<< " ("
	<< fb.second->splittingFn()->interactionType()
	<< ", "
	<< part->evolutionScales()
	[fb.second->splittingFn()->interactionType()]
	<< ">"
	<< fb.first->qtilde()
	<< "; "
	<< fb.first->z()
	<< ")"
	<< endl;
	} // minimal debug
	} // if ( splitFun )
	} // if (shokin && sudakov) {...} else {
	} while (!particlesYetToShower.empty());
	} // while (vetoed)
	return hasEmitted;
	}



	bool ForwardEvolver::MEVeto(tcPPtr p, const Energy &q, const double &z) {
	bool veto = true;
	double weight = 0.;
	if ( HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << " MEVeto: "
	<< " from (id, E, m) = ("
	<< p->id() << ", "
	<< p->momentum().e() << ", "
	<< p->momentum().m() << ") ";
	}
	QQbarG qqg(2*p->momentum().e(), p->momentum().m());
	if (abs(p->id()) < 7) {
	if (p->id()>0) weight = qqg.qWeightX(q, z);
	else weight = qqg.qbarWeightX(q, z);
	veto = !rndbool(weight);
	}
	if ( HERWIG_DEBUG_LEVEL >= HwDebug::full_Shower ) {
	generator()->log() << (veto ? "VETO! ":"passed ")
	<< "w = " << weight << std::endl;
	}
	return veto;
	}

File Metadata

Mime Type: text/x-c
Expires: Tue, Sep 30, 5:55 AM (1 d, 15 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 6551760
Default Alt Text: ForwardEvolver.cc (12 KB)

ForwardEvolver.ccNo OneTemporaryActions

ForwardEvolver.ccView Options

File Metadata

Event Timeline

ForwardEvolver.cc
No OneTemporary
Actions

ForwardEvolver.cc
View Options