]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EVGEN/AliGenHIJINGpara.cxx
added the QA checkable-chekers using TTask and TFolder
[u/mrichter/AliRoot.git] / EVGEN / AliGenHIJINGpara.cxx
index 3efb3a7627d621fb2b47ddeedd8d9af539d4dd02..c7044d974dd7812c5200c83714ac102e4c5bf5e0 100644 (file)
 
 /*
 $Log$
+Revision 1.5  2000/12/21 16:24:06  morsch
+Coding convention clean-up
+
+Revision 1.4  2000/11/30 07:12:50  alibrary
+Introducing new Rndm and QA classes
+
+Revision 1.3  2000/10/02 21:28:06  fca
+Removal of useless dependecies via forward declarations
+
+Revision 1.2  2000/07/11 18:24:55  fca
+Coding convention corrections + few minor bug fixes
+
 Revision 1.1  2000/06/09 20:20:30  morsch
 Same class as previously in AliSimpleGen.cxx
 All coding rule violations except RS3 corrected (AM)
 
 */
-///////////////////////////////////////////////////////////////////
-//                                                               //
-//    Generate the final state of the interaction as the input   //
-//    to the MonteCarlo                                          //
-//
+
+// Parameterisation of pi and K, eta and pt distributions
+// used for the ALICE TDRs.
+// eta: according to HIJING (shadowing + quenching)
+// pT : according to CDF measurement at 1.8 TeV
+// Author: andreas.morsch@cern.ch
+
+
 //Begin_Html
 /*
 <img src="picts/AliGeneratorClass.gif">
@@ -41,6 +56,7 @@ All coding rule violations except RS3 corrected (AM)
 ///////////////////////////////////////////////////////////////////
 
 #include "AliGenHIJINGpara.h"
+#include "TF1.h"
 #include "AliRun.h"
 #include "AliConst.h"
 #include "AliPDG.h"
@@ -221,7 +237,8 @@ void AliGenHIJINGpara::Init()
     Float_t phiFrac = (fPhiMax-fPhiMin)/2/TMath::Pi();
     fParentWeight = Float_t(fNpart)/intETASel*ptFrac*phiFrac;
     
-    printf("\n The number of particles in the selected kinematic region corresponds to %f percent of a full event\n ", 100.*fParentWeight);
+    printf("%s: The number of particles in the selected kinematic region corresponds to %f percent of a full event\n ", 
+          ClassName(),100.*fParentWeight);
     
 }
 
@@ -253,7 +270,7 @@ void AliGenHIJINGpara::Generate()
     //
     for (j=0;j<3;j++) origin[j]=fOrigin[j];
     if(fVertexSmear==kPerEvent) {
-       gMC->Rndm(random,6);
+       Rndm(random,6);
        for (j=0;j<3;j++) {
            origin[j]+=fOsigma[j]*TMath::Cos(2*random[2*j]*TMath::Pi())*
                TMath::Sqrt(-2*TMath::Log(random[2*j+1]));
@@ -261,7 +278,7 @@ void AliGenHIJINGpara::Generate()
     }
     for(i=0;i<fNpart;i++) {
        while(1) {
-           gMC->Rndm(random,3);
+           Rndm(random,3);
            if(random[0]<kBorne) {
                part=kPions[Int_t (random[1]*3)];
                ptf=fPtpi;
@@ -282,13 +299,13 @@ void AliGenHIJINGpara::Generate()
            p[1]=pt*TMath::Sin(phi);
            p[2]=pl;
            if(fVertexSmear==kPerTrack) {
-               gMC->Rndm(random,6);
+               Rndm(random,6);
                for (j=0;j<3;j++) {
                    origin[j]=fOrigin[j]+fOsigma[j]*TMath::Cos(2*random[2*j]*TMath::Pi())*
                        TMath::Sqrt(-2*TMath::Log(random[2*j+1]));
                }
            }
-           gAlice->SetTrack(fTrackIt,-1,part,p,origin,polar,0,"Primary",nt,fParentWeight);
+           gAlice->SetTrack(fTrackIt,-1,part,p,origin,polar,0,kPPrimary,nt,fParentWeight);
            break;
        }
     }