public:
// Constructor.
- SigmaMultiparton() {}
+ SigmaMultiparton() : nChan(0),
+ sigmaTsum(0), sigmaUsum(0),
+ pickOther(0), pickedU(0),
+ rndmPtr(0x0) {}
// Destructor.
~SigmaMultiparton() {
public:
// Constructor.
- MultipartonInteractions() : bIsSet(false) {}
+ MultipartonInteractions() :
+ allowRescatter(false), allowDoubleRes(false), canVetoMPI(false),
+ pTmaxMatch(0), alphaSorder(0), alphaEMorder(0), bProfile(0), processLevel(0),
+ bSelScale(0), rescatterMode(0), nQuarkIn(0), nSample(0), enhanceScreening(0),
+ alphaSvalue(0), Kfactor(0), pT0Ref(0), ecmRef(0), ecmPow(0), pTmin(0), coreRadius(0),
+ coreFraction(0), expPow(0), ySepResc(0), deltaYResc(0), sigmaPomP(0), mPomP(0), pPomP(0),
+ mMaxPertDiff(0), mMinPertDiff(0),
+ a1(0), a0now(0), a02now(0), bstepNow(0), a2max(0), b2now(0), enhanceBmax(0), enhanceBnow(0),
+ id1Save(0), id2Save(0),
+ pT2Save(0), x1Save(0), x2Save(0), sHatSave(0), tHatSave(0), uHatSave(0),
+ alpSsave(0), alpEMsave(0), pT2FacSave(0), pT2RenSave(0), xPDF1nowSave(0),
+ xPDF2nowSave(0),
+ dSigmaDtSelSave(0x0),
+ vsc1(0), vsc2(0),
+ hasBaryonBeams(false), hasLowPow(false), globalRecoilFSR(false),
+ iDiffSys(0), nMaxGlobalRecoilFSR(0),
+ eCM(0), sCM(0), pT0(0), pT20(0), pT2min(0), pTmax(0), pT2max(0), pT20R(0), pT20minR(0),
+ pT20maxR(0), pT20min0maxR(0), pT2maxmin(0), sigmaND(0), pT4dSigmaMax(0),
+ pT4dProbMax(0), dSigmaApprox(0), sigmaInt(0),
+ zeroIntCorr(0), normOverlap(0), nAvg(0), kNow(0), normPi(0), bAvg(0), bDiv(0),
+ probLowB(0), radius2B(0), radius2C(0), fracA(0), fracB(0), fracC(0), fracAhigh(0),
+ fracBhigh(0), fracChigh(0), fracABChigh(0), expRev(0), cDiv(0), cMax(0),
+ bIsSet(false), bSetInFirst(false), isAtLowB(false), pickOtherSel(false),
+ id1(0), id2(0), i1Sel(0), i2Sel(0), id1Sel(0), id2Sel(0),
+ bNow(0), enhanceB(0), pT2(0), pT2shift(0), pT2Ren(0), pT2Fac(0), x1(0), x2(0), xT(0), xT2(0),
+ tau(0), y(0), sHat(0), tHat(0), uHat(0), alpS(0), alpEM(0), xPDF1now(0), xPDF2now(0),
+ dSigmaSum(0), x1Sel(0), x2Sel(0), sHatSel(0), tHatSel(0), uHatSel(0),
+ nStep(0), iStepFrom(0), iStepTo(0),
+ eCMsave(0), eStepSize(0), eStepSave(0), eStepFrom(0), eStepTo(0),
+ infoPtr(0x0),
+ rndmPtr(0x0),
+ beamAPtr(0x0),
+ beamBPtr(0x0),
+ couplingsPtr(0x0),
+ partonSystemsPtr(0x0),
+ sigmaTotPtr(0x0),
+ userHooksPtr(0x0),
+ sigma2gg(), sigma2qg(), sigma2qqbarSame(), sigma2qq(),
+ sigma2Sel(0x0),
+ dSigmaDtSel(0x0),
+ alphaS(),
+ alphaEM() {
+ for (int i=0; i<101; i++) sudExpPT[i] = 0;
+ for (int i=0; i<5; i++) {
+ pT0Save[i] = 0;
+ pT4dSigmaMaxSave[i] = 0; pT4dProbMaxSave[i] = 0; sigmaIntSave[i] = 0;
+
+ for (int j=0; j<101; j++) sudExpPTSave[i][j] = 0;
+ zeroIntCorrSave[i] = 0; normOverlapSave[i] = 0;
+ kNowSave[i] = 0; bAvgSave[i] = 0; bDivSave[i] = 0; probLowBSave[i] = 0;
+ fracAhighSave[i] = 0; fracBhighSave[i] = 0; fracChighSave[i] = 0;
+ fracABChighSave[i] = 0; cDivSave[i] = 0; cMaxSave[i] = 0;
+ }
+ }
// Initialize the generation process for given beams.
bool init( bool doMPIinit, int iDiffSysIn, Info* infoPtrIn,