Page MenuHomeHEPForge
Files F7878235

No OneTemporary
Actions

View Options

diff --git a/UnderlyingEvent/MPIXSecReweighter.cc b/UnderlyingEvent/MPIXSecReweighter.cc
--- a/UnderlyingEvent/MPIXSecReweighter.cc
+++ b/UnderlyingEvent/MPIXSecReweighter.cc
@@ -1,98 +1,98 @@
// -*- C++ -*-
//
// MPIXSecReweighter.cc is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the MPIXSecReweighter class.
//
#include "MPIXSecReweighter.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/Repository/UseRandom.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/Utilities/DescribeClass.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "MPIHandler.h"
using namespace Herwig;
MPIXSecReweighter::MPIXSecReweighter()
: sumWeights(0.), xSec(ZERO) {}
MPIXSecReweighter::~MPIXSecReweighter() {}
void MPIXSecReweighter::
handle(EventHandler & eh, const tPVector &,
const Hint &) {
if ( MPIHandler::currentHandler() &&
dynamic_cast<StandardEventHandler*>(&eh) ) {
double corr = 1.;
StandardEventHandler& seh =
dynamic_cast<StandardEventHandler&>(eh);
CrossSection mpiXSec = MPIHandler::currentHandler()->inelasticXSec();
double weight = seh.currentEvent()->weight();
if ( weight == 0. )
return;
CrossSection next = seh.integratedXSecNoReweight();
assert(next != ZERO);
if ( xSec != ZERO ) {
corr = (mpiXSec/next)*( 1. +(sumWeights/weight)*(1.-next/xSec) );
sumWeights += weight;
xSec = next;
} else {
xSec = next;
sumWeights = weight;
corr = mpiXSec/xSec;
}
seh.reweight(corr);
}
}
IBPtr MPIXSecReweighter::clone() const {
return new_ptr(*this);
}
IBPtr MPIXSecReweighter::fullclone() const {
return new_ptr(*this);
}
// If needed, insert default implementations of virtual function defined
// in the InterfacedBase class here (using ThePEG-interfaced-impl in Emacs).
void MPIXSecReweighter::persistentOutput(PersistentOStream &) const {}
void MPIXSecReweighter::persistentInput(PersistentIStream &, int) {}
// *** Attention *** The following static variable is needed for the type
// description system in ThePEG. Please check that the template arguments
// are correct (the class and its base class), and that the constructor
// arguments are correct (the class name and the name of the dynamically
// loadable library where the class implementation can be found).
DescribeClass<MPIXSecReweighter,StepHandler>
- describeHerwigMPIXSecReweighter("Herwig::MPIXSecReweighter", "SimpleKTCut.so HwMPI.so");
+ describeHerwigMPIXSecReweighter("Herwig::MPIXSecReweighter", "JetCuts.so SimpleKTCut.so HwMPI.so");
void MPIXSecReweighter::Init() {
static ClassDocumentation<MPIXSecReweighter> documentation
("MPIXSecReweighter sets up the proper minimum bias cross section.");
}

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 5:29 PM (1 d, 14 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805390
Default Alt Text
(3 KB)

Event Timeline