/*
$Log$
+Revision 1.3 2003/03/31 15:47:16 morsch
+Unhandled cases added.
+
+Revision 1.2 2003/03/27 09:28:38 morsch
+Charm production process configuration for pPb and pp added
+(A. Dainese, N. Carrer)
+
Revision 1.1 2003/03/15 15:00:48 morsch
Classed imported from EVGEN.
fFragmentation = kTRUE;
// Default counting mode
fCountMode = kCountAll;
+ // Pycel
+ SetPycellParameters();
}
AliGenPythia::AliGenPythia(const AliGenPythia & Pythia)
// Destructor
}
+void AliGenPythia::SetPycellParameters(Float_t etamax, Int_t neta, Int_t nphi,
+ Float_t thresh, Float_t etseed, Float_t minet, Float_t r)
+{
+// Set pycell parameters
+ fPycellEtaMax = etamax;
+ fPycellNEta = neta;
+ fPycellNPhi = nphi;
+ fPycellThreshold = thresh;
+ fPycellEtSeed = etseed;
+ fPycellMinEtJet = minet;
+ fPycellMaxRadius = r;
+}
+
+
+
void AliGenPythia::SetEventListRange(Int_t eventFirst, Int_t eventLast)
{
// Set a range of event numbers, for which a table
case kPyCharm:
case kPyCharmUnforced:
case kPyCharmPbPbMNR:
+ case kPyCharmppMNR:
+ case kPyCharmpPbMNR:
fParentSelect[0] = 411;
fParentSelect[1] = 421;
fParentSelect[2] = 431;
break;
}
//
+//
+// JetFinder for Trigger
+//
+// Configure detector (EMCAL like)
+//
+ fPythia->SetPARU(51, fPycellEtaMax);
+ fPythia->SetMSTU(51, fPycellNEta);
+ fPythia->SetMSTU(52, fPycellNPhi);
+//
+// Configure Jet Finder
+//
+ fPythia->SetPARU(58, fPycellThreshold);
+ fPythia->SetPARU(52, fPycellEtSeed);
+ fPythia->SetPARU(53, fPycellMinEtJet);
+ fPythia->SetPARU(54, fPycellMaxRadius);
+ fPythia->SetMSTU(54, 2);
+//
// This counts the total number of calls to Pyevnt() per run.
fTrialsRun = 0;
fQ = 0.;
// Calls the Pythia jet finding algorithm to find jets in the current event
//
//
-// Configure detector (EMCAL like)
-//
- fPythia->SetPARU(51,2.);
- fPythia->SetMSTU(51,Int_t(96 * 2./0.7));
- fPythia->SetMSTU(52,3 * 144);
-//
-// Configure Jet Finder
-//
- fPythia->SetPARU(58, eCellMin);
- fPythia->SetPARU(52, eCellSeed);
- fPythia->SetPARU(53, eMin);
- fPythia->SetPARU(54, rMax);
- fPythia->SetMSTU(54, 2);
//
// Save jets
Int_t n = fPythia->GetN();
// Loading jets from common block
//
} else {
-//
-// Configure detector (EMCAL like)
-//
- fPythia->SetPARU(51,2.);
- fPythia->SetMSTU(51,Int_t(96 * 2./0.7));
- fPythia->SetMSTU(52,3 * 144);
-//
-// Configure Jet Finder
-//
- fPythia->SetPARU(58, 0.0);
- fPythia->SetPARU(52, 4.0);
- fPythia->SetPARU(53, 10.0);
- fPythia->SetPARU(54, 1.0);
- fPythia->SetMSTU(54, 2);
+
//
// Run Jet Finder
fPythia->Pycell(nJets);