/*
$Log$
+Revision 1.5 1999/09/29 09:24:12 fca
+Introduction of the Copyright and cvs Log
+
*/
+#include <iostream.h>
+
#include "AliGenExtFile.h"
#include "AliGenMUONlib.h"
#include "AliMC.h"
/*
$Log$
+Revision 1.13 1999/11/04 11:30:31 fca
+Correct the logics for SetForceDecay
+
Revision 1.12 1999/11/03 17:43:20 fca
New version from G.Martinez & A.Morsch
// particle type
Int_t Ipart = fIpParaFunc();
fChildWeight=(fPythia->GetBraPart(Ipart))*fParentWeight;
- Float_t am=fPythia->GetPMAS(fPythia->LuComp(Ipart),1);
+ Float_t am=fPythia->GetPMAS(fPythia->Lucomp(Ipart),1);
gMC->Rndm(random,2);
//
// phi
/*
$Log$
+Revision 1.12 1999/11/03 17:43:20 fca
+New version from G.Martinez & A.Morsch
+
Revision 1.11 1999/09/29 09:24:14 fca
Introduction of the Copyright and cvs Log
*/
fPythia->SetCKIN(4,fPtHardMax);
fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc);
fPythia->ForceDecay(fForceDecay);
- fPythia->LuList(0);
- fPythia->PyStat(2);
+ fPythia->Lulist(0);
+ fPythia->Pystat(2);
// Parent and Children Selection
switch (fProcess)
{
}
while(1)
{
- fPythia->PyEvnt();
+ fPythia->Pyevnt();
fTrials++;
fPythia->ImportParticles(particles,"All");
Int_t np = particles->GetEntriesFast();
/*
$Log$
+Revision 1.5 1999/11/03 17:43:20 fca
+New version from G.Martinez & A.Morsch
+
Revision 1.4 1999/09/29 09:24:14 fca
Introduction of the Copyright and cvs Log
fEcms = energy;
fStrucFunc = strucfunc;
// don't decay p0
- SetMDCY(LuComp(111),1,0);
+ SetMDCY(Lucomp(111),1,0);
// select structure function
SetMSTP(52,2);
SetMSTP(51,strucfunc);
//
// force decay of particle into products with multiplicity mult
- Int_t kc=LuComp(particle);
+ Int_t kc=Lucomp(particle);
SetMDCY(kc,1,1);
Int_t ifirst=GetMDCY(kc,2);
Int_t ilast=ifirst+GetMDCY(kc,3)-1;
//
// phi-> mu+mu- and phi -> e+e-
// clone the original phi
- kc = LuComp(333);
+ kc = Lucomp(333);
nkc = 41;
for (i=1;i<=3;i++) {
// Vector mesons
//
// phi clone for dilepton decay-channel
- kc=LuComp(41);
+ kc=Lucomp(41);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(113,"Phi",3,mass,0,tlife);
fGPCode[kc][0]=113;
fGPCode[kc][1]=113;
// J/Psi
- kc=LuComp(443);
+ kc=Lucomp(443);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(114,"J/Psi",3,mass,0,tlife);
fGPCode[kc][0]=114;
fGPCode[kc][1]=114;
// psi prime
- kc=LuComp(30443);
+ kc=Lucomp(30443);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(115,"Psi'",3,mass,0,tlife);
fGPCode[kc][0]=115;
fGPCode[kc][1]=115;
// upsilon(1s)
- kc=LuComp(553);
+ kc=Lucomp(553);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(116,"Upsilon",3,mass,0,tlife);
fGPCode[kc][0]=116;
fGPCode[kc][1]=116;
// upsilon(2s)
- kc=LuComp(30553);
+ kc=Lucomp(30553);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(117,"Upsilon'",3,mass,0,tlife);
fGPCode[kc][0]=117;
fGPCode[kc][1]=117;
// upsilon(3s)
- kc=LuComp(30553);
+ kc=Lucomp(30553);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(118,"Upsilon''",3,mass,0,tlife);
// charmed mesons
//
// D^+/-
- kc=LuComp(411);
+ kc=Lucomp(411);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(119,"D^+",3,mass, 1,tlife);
fGPCode[kc][0]=119;
fGPCode[kc][1]=120;
// D^0
- kc=LuComp(421);
+ kc=Lucomp(421);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(121,"D^0",3,mass,0,tlife);
fGPCode[kc][0]=121;
fGPCode[kc][1]=122;
// D_s
- kc=LuComp(431);
+ kc=Lucomp(431);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(123,"D_s^+",3,mass, 1,tlife);
fGPCode[kc][0]=123;
fGPCode[kc][1]=124;
// Lambda_c
- kc=LuComp(4122);
+ kc=Lucomp(4122);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(125,"Lambda_c+",3,mass, 1,tlife);
//
// beauty mesons
// B_0
- kc=LuComp(511);
+ kc=Lucomp(511);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(127,"B^0",3,mass, 0,tlife);
fGPCode[kc][0]=127;
fGPCode[kc][1]=128;
// B^+-
- kc=LuComp(521);
+ kc=Lucomp(521);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(129,"B^+",3,mass, 1,tlife);
fGPCode[kc][0]=129;
fGPCode[kc][1]=130;
// B_s
- kc=LuComp(531);
+ kc=Lucomp(531);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(131,"B_s",3,mass, 0,tlife);
fGPCode[kc][0]=131;
fGPCode[kc][1]=132;
// Lambda_b
- kc=LuComp(5122);
+ kc=Lucomp(5122);
mass =GetPMAS(kc,1);
tlife=GetPMAS(kc,4);
gMC->Gspart(133,"Lambda_b",3,mass, 0,tlife);
// set up correspondance between standard GEANT particle codes
// and PYTHIA kf
- kc=LuComp(22); // gamma
+ kc=Lucomp(22); // gamma
fGPCode[kc][0]=1;
fGPCode[kc][1]=1;
- kc=LuComp(11); // positron
+ kc=Lucomp(11); // positron
fGPCode[kc][0]=2;
fGPCode[kc][1]=3;
- kc=LuComp(12); // neutrino
+ kc=Lucomp(12); // neutrino
fGPCode[kc][0]=4;
fGPCode[kc][1]=4;
- kc=LuComp(13); // muon
+ kc=Lucomp(13); // muon
fGPCode[kc][0]=5;
fGPCode[kc][1]=6;
- kc=LuComp(111); // pi0
+ kc=Lucomp(111); // pi0
fGPCode[kc][0]=7;
fGPCode[kc][1]=7;
- kc=LuComp(211); // pi+
+ kc=Lucomp(211); // pi+
fGPCode[kc][0]=8;
fGPCode[kc][1]=9;
- kc=LuComp(130); // K0 short
+ kc=Lucomp(130); // K0 short
fGPCode[kc][0]=10;
fGPCode[kc][1]=10;
- kc=LuComp(321); // K+/-
+ kc=Lucomp(321); // K+/-
fGPCode[kc][0]=11;
fGPCode[kc][1]=12;
- kc=LuComp(2112); // neutron/anti-neutron
+ kc=Lucomp(2112); // neutron/anti-neutron
fGPCode[kc][0]=13;
fGPCode[kc][1]=25;
- kc=LuComp(2212); // proton/anti-proton
+ kc=Lucomp(2212); // proton/anti-proton
fGPCode[kc][0]=14;
fGPCode[kc][1]=15;
- kc=LuComp(310); // K0 short
+ kc=Lucomp(310); // K0 short
fGPCode[kc][0]=16;
fGPCode[kc][1]=16;
- kc=LuComp(221); // eta
+ kc=Lucomp(221); // eta
fGPCode[kc][0]=17;
fGPCode[kc][1]=17;
- kc=LuComp(3122); // lambda
+ kc=Lucomp(3122); // lambda
fGPCode[kc][0]=18;
fGPCode[kc][1]=18;
- kc=LuComp(3222); // sigma+/antisigma+
+ kc=Lucomp(3222); // sigma+/antisigma+
fGPCode[kc][0]=19;
fGPCode[kc][1]=29;
- kc=LuComp(3212); // sigma0/antisigma0
+ kc=Lucomp(3212); // sigma0/antisigma0
fGPCode[kc][0]=20;
fGPCode[kc][1]=28;
- kc=LuComp(3112); // sigma-/antisigma-
+ kc=Lucomp(3112); // sigma-/antisigma-
fGPCode[kc][0]=21;
fGPCode[kc][1]=27;
- kc=LuComp(3322); // xsi0-/antixsi0
+ kc=Lucomp(3322); // xsi0-/antixsi0
fGPCode[kc][0]=22;
fGPCode[kc][1]=30;
- kc=LuComp(3312); // xsi-/antixsi+
+ kc=Lucomp(3312); // xsi-/antixsi+
fGPCode[kc][0]=23;
fGPCode[kc][1]=31;
- kc=LuComp(3334); // omega/antiomega
+ kc=Lucomp(3334); // omega/antiomega
fGPCode[kc][0]=24;
fGPCode[kc][1]=32;
}
Int_t AliPythia::GetGeantCode(Int_t kf)
{
- Int_t kc=LuComp(TMath::Abs(kf));
+ Int_t kc=Lucomp(TMath::Abs(kf));
return (kf > 0) ? fGPCode[kc][0] : fGPCode[kc][1];
}
Float_t AliPythia::GetBraPart(Int_t kf)
{
- Int_t kc=LuComp(TMath::Abs(kf));
+ Int_t kc=Lucomp(TMath::Abs(kf));
return fBraPart[kc];
}
// convert to compressed code and print result (for debugging only)
virtual Int_t CheckedLuComp(Int_t kf)
{
- Int_t kc=LuComp(kf);
+ Int_t kc=Lucomp(kf);
printf("\n Lucomp kf,kc %d %d",kf,kc);
return kc;
}
/*
$Log$
+Revision 1.8 1999/09/29 09:24:23 fca
+Introduction of the Copyright and cvs Log
+
*/
#include <TSlider.h>
#include <TSliderBox.h>
#include <TGaxis.h>
-#include <TGXW.h>
+#include <TVirtualX.h>
#include <TMath.h>
#include <TMatrix.h>
#include <X3DBuffer.h>
if (ysize < 100) ysize = 750;
Int_t xsize = Int_t(size*830./ysize);
fCanvas = new TCanvas("Canvas", "MUON Clusters Display",14,47,xsize,ysize);
- fCanvas->SetEditable(kIsNotEditable);
+ fCanvas->SetEditable(kFALSE);
fCanvas->ToggleEventStatus();
// Create main display pad
switch (event) {
case kButton1Down:
- gGXW->SetLineColor(-1);
+ gVirtualX->SetLineColor(-1);
gPad->TAttLine::Modify(); //Change line attributes only if necessary
x0 = gPad->AbsPixeltoX(px);
y0 = gPad->AbsPixeltoY(py);
return;
case kButton1Motion:
- if (linedrawn) gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ if (linedrawn) gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
pxold = px;
pyold = py;
linedrawn = 1;
- gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
return;
case kButton1Up:
/*
$Log$
+Revision 1.5 1999/09/29 09:24:29 fca
+Introduction of the Copyright and cvs Log
+
*/
#include <TSlider.h>
#include <TSliderBox.h>
#include <TGaxis.h>
-#include <TGXW.h>
+#include <TVirtualX.h>
#include <TMath.h>
#include <X3DBuffer.h>
if (ysize < 100) ysize = 750;
Int_t xsize = Int_t(size*830./ysize);
fCanvas = new TCanvas("Canvas", "RICH Clusters Display",14,47,xsize,ysize);
- fCanvas->SetEditable(kIsNotEditable);
+ fCanvas->SetEditable(kFALSE);
fCanvas->ToggleEventStatus();
// Create main display pad
switch (event) {
case kButton1Down:
- gGXW->SetLineColor(-1);
+ gVirtualX->SetLineColor(-1);
gPad->TAttLine::Modify(); //Change line attributes only if necessary
x0 = gPad->AbsPixeltoX(px);
y0 = gPad->AbsPixeltoY(py);
return;
case kButton1Motion:
- if (linedrawn) gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ if (linedrawn) gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
pxold = px;
pyold = py;
linedrawn = 1;
- gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
return;
case kButton1Up:
/*
$Log$
+Revision 1.5 1999/09/29 09:24:29 fca
+Introduction of the Copyright and cvs Log
+
*/
#include <TSlider.h>
#include <TSliderBox.h>
#include <TGaxis.h>
-#include <TGXW.h>
+#include <TVirtualX.h>
#include <TMath.h>
#include <X3DBuffer.h>
if (ysize < 100) ysize = 750;
Int_t xsize = Int_t(size*830./ysize);
fCanvas = new TCanvas("Canvas", "ALICE Event Display",14,47,xsize,ysize);
- fCanvas->SetEditable(kIsNotEditable);
+ fCanvas->SetEditable(kFALSE);
fCanvas->ToggleEventStatus();
// Create main display pad
switch (event) {
case kButton1Down:
- gGXW->SetLineColor(-1);
+ gVirtualX->SetLineColor(-1);
gPad->TAttLine::Modify(); //Change line attributes only if necessary
x0 = gPad->AbsPixeltoX(px);
y0 = gPad->AbsPixeltoY(py);
return;
case kButton1Motion:
- if (linedrawn) gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ if (linedrawn) gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
pxold = px;
pyold = py;
linedrawn = 1;
- gGXW->DrawBox(px0, py0, pxold, pyold, TGXW::kHollow);
+ gVirtualX->DrawBox(px0, py0, pxold, pyold, TVirtualX::kHollow);
return;
case kButton1Up:
# Current Version
-VERSION = 3.02pre
+VERSION = 3.02
# Which makedepend
MAKEDEPEND = $(ROOTSYS)/bin/rmkdepend
# Root libraries
-ROOTLIBS = -L$(ROOTSYS)/lib -lNew -lBase -lCint -lClib -lCont -lFunc \
- -lGraf -lGraf3d -lHist -lHtml -lMeta -lMinuit -lNet \
- -lPostscript -lProof -lRint -lTree -lUnix -lZip -lMatrix \
- -lMeta -lPhysics
+ROOTLIBS = $(shell $(ROOTSYS)/bin/root-config --glibs) -lMinuit -lHtml
# Alice libraries
-lITS -lPHOSdummy -lPHOS -lCASTOR -lRICH -lZDC -lSTRUCT \
-lTGeant3 -lSTART -lCPV -lEVGEN -lSTEER
-# Graphic libraries
-
-ROOTGLIBS = -lGpad -lGui -lGX11 -lX3d -lXpm
-
# CERN Libraries
CLIBS = -lminicern
G3LIBS = -lgeant321
+
SHLD = $(LD)
SOFLAGS = -b /usr/lib/libm.sl
-LIBS = $(ROOTLIBS) -lNew -lm
-GLIBS = $(ROOTLIBS) $(ROOTGLIBS) -L/usr/lib/X11R6 -L/usr/local/lib \
- -lX11 -lm -ldld
+LIBS = $(ROOTLIBS)
+GLIBS = $(ROOTLIBS) -L/usr/lib/X11R6 -lX11
SHLD = $(LD)
SOFLAGS = -b /usr/lib/libm.sl
-LIBS = $(ROOTLIBS) -lNew -lm
-GLIBS = $(ROOTLIBS) $(ROOTGLIBS) -L/usr/lib/X11R6 -L/usr/local/lib \
- -lX11 -lm -ldld
+LIBS = $(ROOTLIBS)
+GLIBS = $(ROOTLIBS) -L/usr/lib/X11R6 -lX11
SOFLAGS = -Wl,-soname,$(notdir $@) -shared
SHLIB = -lg2c
-LIBS = $(ROOTLIBS) -lm -ldl -rdynamic
-GLIBS = $(ROOTLIBS) $(ROOTGLIBS) -L/usr/X11R6/lib -lX11 -lg2c -ldl -lcrypt -rdynamic
+LIBS = $(ROOTLIBS)
+GLIBS = $(ROOTLIBS) -L/usr/X11R6/lib -lX11 -lg2c -lcrypt