//_____________________________________________________________________________
AliGenBox::AliGenBox()
:AliGenerator(),
- fIpart(0)
+ fIpart(0),
+ fEtaMin(0),
+ fEtaMax(0)
{
//
// Default constructor
//_____________________________________________________________________________
AliGenBox::AliGenBox(Int_t npart)
:AliGenerator(npart),
- fIpart(kProton)
+ fIpart(kProton),
+ fEtaMin(0),
+ fEtaMax(0)
{
//
// Standard constructor
for(i=0;i<fNpart;i++) {
Rndm(random,3);
- theta=fThetaMin+random[0]*(fThetaMax-fThetaMin);
+
+ if (TestBit(kThetaRange)) {
+ theta = fThetaMin+random[0]*(fThetaMax-fThetaMin);
+ } else {
+ Float_t eta = fEtaMin+random[0]*(fEtaMax-fEtaMin);
+ theta = 2. * TMath::ATan(TMath::Exp(-eta));
+ }
+
+
if(TestBit(kMomentumRange)) {
pmom=fPMin+random[1]*(fPMax-fPMin);
pt=pmom*TMath::Sin(theta);
pt=fPtMin+random[1]*(fPtMax-fPtMin);
pmom=pt/TMath::Sin(theta);
}
+
phi=fPhiMin+random[2]*(fPhiMax-fPhiMin);
p[0] = pt*TMath::Cos(phi);
p[1] = pt*TMath::Sin(phi);
AliGenEventHeader* header = new AliGenEventHeader("BOX");
header->SetPrimaryVertex(fVertex);
- gAlice->SetGenEventHeader(header);
+ header->SetNProduced(fNpart);
+
+ // Passes header either to the container or to gAlice
+ if (fContainer) {
+ fContainer->AddHeader(header);
+ } else {
+ gAlice->SetGenEventHeader(header);
+ }
}
//_____________________________________________________________________________