Update master to aliroot
[u/mrichter/AliRoot.git] / EPOS / TEpos.cxx
index 8f7c0bb..b894541 100644 (file)
@@ -1,18 +1,18 @@
-/*
- *###################################################################
- *#        EPOS 1.67     K. WERNER, T. PIEROG, S. PORTEBOEUF.       #
- *#                      Contact: werner@subatech.in2p3.fr          #
- *###################################################################
- *
- * TEpos.cxx
- * 
- * Wraper class for interfacing EPOS model, derived from ROOT's TGenerator.
- * It generates temporary input file for the model, providing user with
- * ability to add his/hers own lines to the input.
- * Output is read directly from common blocks.
- *
- *      Author: Piotr Ostrowski, postrow@if.pw.edu.pl
- */
+//
+//###################################################################
+//#        EPOS 1.67     K. WERNER, T. PIEROG, S. PORTEBOEUF.       #
+//#                      Contact: werner@subatech.in2p3.fr          #
+//###################################################################
+//
+// TEpos.cxx
+// 
+// Wraper class for interfacing EPOS model, derived from ROOT's TGenerator.
+// It generates temporary input file for the model, providing user with
+// ability to add his/hers own lines to the input.
+// Output is read directly from common blocks.
+//
+//      Author: Piotr Ostrowski, postrow@if.pw.edu.pl
+//
 
 
 #include <TClonesArray.h>
@@ -75,14 +75,37 @@ TEpos::~TEpos() {
        delete fIdConverter;
 }
 
-TEpos& TEpos::operator=(const TEpos&) {
-       if (!fIdConverter) {
-               fIdConverter = new AliGenEposIsajetToPdgConverter();
-       }
-       return *this;
+TEpos& TEpos::operator=(const TEpos& epos) {
+  //operator=
+  if (this != &epos) {
+    if (fIdConverter) {
+       delete fIdConverter;
+       fIdConverter = new AliGenEposIsajetToPdgConverter(*epos.fIdConverter);
+    } else{
+       fIdConverter = 0;
+    }
+    fLaproj = epos.fLaproj;
+    fMaproj = epos.fMaproj;
+    fLatarg = epos.fLatarg;
+    fMatarg = epos.fMatarg;
+    fBminim = epos.fBminim;
+    fBmaxim = epos.fBmaxim;
+    fPhimin = epos.fPhimin;
+    fPhimax = epos.fPhimax;
+    fEcms   = epos.fEcms;
+    fSplitting = epos.fSplitting;
+    for(UInt_t i = 0; i < epos.fNoDecays.size(); ++i) {
+      fNoDecays.push_back(epos.fNoDecays[i]);
+    }
+    for(UInt_t i = 0; i < epos.fExtraInputLines.size(); ++i) {
+      fExtraInputLines.push_back(epos.fExtraInputLines[i]);
+    }
+  }
+  return *this;
 }
 
 void TEpos::Initialize() {
+  // Generates input file and prepares EPOS to read from it.
        Int_t nopeno = 0;
        GenerateInputFile();
        aaset(0);
@@ -105,16 +128,17 @@ void TEpos::Initialize() {
                IniModel(appli.model);
        ebin.nrebin = 1;
        ainit();
+       aseed(2);
 }
 
 void TEpos::GenerateEvent() {
-       cseed.seedj = gRandom->Rndm() * 1e10;
-       aseed(2);
+//     cseed.seedj = gRandom->Rndm() * 1e10;
        Int_t n = 1;
        evgen(n);
 }
 
 Int_t TEpos::ImportParticles(TClonesArray *particles, Option_t *) {
+  //Fills provided ClonesArray with generated particles
        particles->Clear();
        if (!cevt.nevt) return 0;
        Int_t numpart = cptl.nptl;
@@ -141,6 +165,7 @@ Int_t TEpos::ImportParticles(TClonesArray *particles, Option_t *) {
 }
 
 TObjArray*  TEpos::ImportParticles(Option_t *) {
+  //Creates new particle array
        fParticles->Clear();
        if (!cevt.nevt) return NULL;
        Int_t numpart = cptl.nptl;
@@ -176,15 +201,16 @@ void TEpos::AddExtraInputLine(const char *line) {
 }
 
 void TEpos::GenerateInputFile() {
+  // Generate input file in EPOS format
        ofstream file(GetInputFileName(), ios_base::out | ios_base::trunc);
        char epo[256];
        char *epoEnv = getenv("EPO");
        if (epoEnv) {
-               strcpy(epo, epoEnv);
+               strncpy(epo, epoEnv, 255);
        } else {
-               strcpy(epo, getenv("ALICE_ROOT"));
+               strncpy(epo, getenv("ALICE_ROOT"), 255);
        }
-       strcat(epo, "/EPOS/epos167");
+       strncat(epo, "/EPOS/epos167", 255);
 
        file << "fname pathnx " << epo << "/" << endl;
        file << "fname histo none" << endl;
@@ -251,25 +277,25 @@ void TEpos::GenerateInputFile() {
        file.close();
 }
 
-Float_t TEpos::GetBimevt() { return cevt.bimevt; }
-Float_t TEpos::GetPhievt() { return cevt.phievt; }
-Int_t TEpos::GetKolevt() { return cevt.kolevt; }
-Int_t TEpos::GetKoievt() { return cevt.koievt; }
-Float_t TEpos::GetPmxevt() { return cevt.pmxevt; }
-Float_t TEpos::GetEgyevt() { return cevt.egyevt; }
-Int_t TEpos::GetNpjevt() { return cevt.npjevt; }
-Int_t TEpos::GetNtgevt() { return cevt.ntgevt; }
-Int_t TEpos::GetNpnevt() { return cevt.npnevt; }
-Int_t TEpos::GetNppevt() { return cevt.nppevt; }
-Int_t TEpos::GetNtnevt() { return cevt.ntnevt; }
-Int_t TEpos::GetNtpevt() { return cevt.ntpevt; }
-Int_t TEpos::GetJpnevt() { return cevt.jpnevt; }
-Int_t TEpos::GetJppevt() { return cevt.jppevt; }
-Int_t TEpos::GetJtnevt() { return cevt.jtnevt; }
-Int_t TEpos::GetJtpevt() { return cevt.jtpevt; }
-Float_t TEpos::GetXbjevt() { return cevt.xbjevt; }
-Float_t TEpos::GetQsqevt() { return cevt.qsqevt; }
-Int_t TEpos::GetNglevt() { return cevt.nglevt; }
-Float_t TEpos::GetZppevt() { return cevt.zppevt; }
-Float_t TEpos::GetZptevt() { return cevt.zptevt; }
+Float_t TEpos::GetBimevt() const { return cevt.bimevt; }
+Float_t TEpos::GetPhievt() const { return cevt.phievt; }
+Int_t TEpos::GetKolevt() const { return cevt.kolevt; }
+Int_t TEpos::GetKoievt() const { return cevt.koievt; }
+Float_t TEpos::GetPmxevt() const { return cevt.pmxevt; }
+Float_t TEpos::GetEgyevt() const { return cevt.egyevt; }
+Int_t TEpos::GetNpjevt() const { return cevt.npjevt; }
+Int_t TEpos::GetNtgevt() const { return cevt.ntgevt; }
+Int_t TEpos::GetNpnevt() const { return cevt.npnevt; }
+Int_t TEpos::GetNppevt() const { return cevt.nppevt; }
+Int_t TEpos::GetNtnevt() const { return cevt.ntnevt; }
+Int_t TEpos::GetNtpevt() const { return cevt.ntpevt; }
+Int_t TEpos::GetJpnevt() const { return cevt.jpnevt; }
+Int_t TEpos::GetJppevt() const { return cevt.jppevt; }
+Int_t TEpos::GetJtnevt() const { return cevt.jtnevt; }
+Int_t TEpos::GetJtpevt() const { return cevt.jtpevt; }
+Float_t TEpos::GetXbjevt() const { return cevt.xbjevt; }
+Float_t TEpos::GetQsqevt() const { return cevt.qsqevt; }
+Int_t TEpos::GetNglevt() const { return cevt.nglevt; }
+Float_t TEpos::GetZppevt() const { return cevt.zppevt; }
+Float_t TEpos::GetZptevt() const { return cevt.zptevt; }