////////////////////////////////////////////////////////////////////////////////
-#include <Riostream.h>
+#include <iostream.h>
#include "TROOT.h"
#include "TCanvas.h"
#include "AliGenGeVSim.h"
#include "AliGeVSimParticle.h"
+#include "AliGenGeVSimEventHeader.h"
ClassImp(AliGenGeVSim);
if (psi < 0 || psi > 360 )
Error ("AliGenGeVSim", "Reaction plane angle ( %d )out of range [0..360]", psi);
- fPsi = psi * 2 * TMath::Pi() / 360. ;
+ fPsi = psi * TMath::Pi() / 180. ;
fIsMultTotal = isMultTotal;
// initialization
form = 0;
form = (TF1 *)gROOT->GetFunction("gevsimPsi");
- if (form) fPsi = form->Eval(gAlice->GetEvNumber());
+ if (form) fPsi = form->Eval(gAlice->GetEvNumber()) * TMath::Pi() / 180;
form = 0;
form = (TF1 *)gROOT->GetFunction("gevsimPsiRndm");
- if (form) fPsi = form->GetRandom();
+ if (form) fPsi = form->GetRandom() * TMath::Pi() / 180;
+
+
+ cout << "Psi = " << fPsi << "\t" << (Int_t)(fPsi*180./TMath::Pi()) << endl;
fPhiFormula->SetParameter(0, fPsi);
}
fPtYHist = (TH2D*)gROOT->FindObject(buff);
}
- if (!fPtYHist) Error(where, msg[3], pdg);
+ if (!fPtYHist) Error(where, msg[4], pdg);
}
}
}
}
+ // prepare and store header
+
+ AliGenGeVSimEventHeader *header = new AliGenGeVSimEventHeader("GeVSim header");
+ TArrayF eventVertex(3,orgin);
+
+ header->SetPrimaryVertex(eventVertex);
+ header->SetEventPlane(fPsi);
+ header->SetEllipticFlow(fPhiFormula->GetParameter(2));
+
+ gAlice->SetGenEventHeader(header);
+
delete v;
}