# define pyclus pyclus_
# define pycell pycell_
# define pyshow pyshow_
+# define pyshowq pyshowq_
# define pyrobo pyrobo_
# define pyquen pyquen_
# define pyevnw pyevnw_
# define pyjoin pyjoin_
+# define qpygin0 qpygin0_
# define type_of_call
#else
# define pyclus PYCLUS
# define pycell PYCELL
+# define pyshow PYSHOW
+# define pyshowq PYSHOWQ
# define pyrobo PYROBO
# define pyquen PYQUEN
# define pyevnw PYEVNW
# define pyjoin PYJOIN
+# define qpygin0 QPYGIN0
# define type_of_call _stdcall
#endif
extern "C" void type_of_call pyclus(Int_t & );
extern "C" void type_of_call pycell(Int_t & );
extern "C" void type_of_call pyshow(Int_t &, Int_t &, Double_t &);
+extern "C" void type_of_call pyshowq(Int_t &, Int_t &, Double_t &);
+extern "C" void type_of_call qpygin0();
extern "C" void type_of_call pyrobo(Int_t &, Int_t &, Double_t &, Double_t &, Double_t &, Double_t &, Double_t &);
extern "C" void type_of_call pyquen(Double_t &, Int_t &, Double_t &);
extern "C" void type_of_call pyevnw();
pythia.Copy(*this);
}
-void AliPythia6::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc)
+void AliPythia6::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t /*tune*/)
{
// Initialise the process to generate
if (!AliPythiaRndm::GetPythiaRandom())
SetMSUB(94,1); // double diffraction
SetMSUB(95,1); // low pt production
- AtlasTuning_MC09();
+ AtlasTuningMC09();
break;
case kPyMbWithDirectPhoton:
// With parton showers on we are generating "Z inclusive process"
SetMSTP(61,1); //Initial QCD & QED showers on
SetMSTP(71,1); //Final QCD & QED showers on
-
+ break;
+ case kPyMBRSingleDiffraction:
+ case kPyMBRDoubleDiffraction:
+ case kPyMBRCentralDiffraction:
break;
}
pyjoin(npart, ipart);
}
+void AliPythia6::Pyshowq(Int_t ip1, Int_t ip2, Double_t qmax)
+{
+// Call qPythia showering
+//
+ pyshowq(ip1, ip2, qmax);
+}
+
+void AliPythia6::Qpygin0()
+{
+ //position of the hard scattering in the nuclear overlapping area.
+ //just for qpythia.
+ qpygin0();
+}
+
void AliPythia6::Pyclus(Int_t& njet)
{
// Call Pythia clustering algorithm
Int_t itype = (qPdg[j] == 21) ? 2 : 1;
- Double_t eloss = fQuenchingWeights->GetELossRandomKFast(itype, int0[j], int1[j], eq[j]);
+ // Double_t eloss = fQuenchingWeights->GetELossRandomKFast(itype, int0[j], int1[j], eq[j]);
+ Double_t eloss = fQuenchingWeights->GetELossRandomK(itype, int0[j], int1[j], eq[j]);
if (TMath::Abs(yq[j]) > 2.5 || eq[j] < 10.) {
fZQuench[j] = 0.;
SetPARP(67,1); // Regulates Initial State Radiation
}
-void AliPythia6::AtlasTuning_MC09()
+void AliPythia6::AtlasTuningMC09()
{
//
// Configuration for the ATLAS tuning