/*
$Log$
+Revision 1.11 1999/10/05 17:15:45 fca
+Minor syntax for the Alpha OSF
+
Revision 1.10 1999/10/01 09:24:40 fca
Protect against no current file in FinishEvent
printf("Dummy version of muon step -- it should never happen!!\n");
/*
const Float_t kRaddeg = 180/TMath::Pi();
- AliMC* pMC = AliMC::GetMC();
Int_t nsec, ipart;
Float_t x[4], p[4];
Float_t pt, th0, th2;
char proc[5];
if(fAccCut) {
- if((nsec=pMC->NSecondaries())>0) {
- pMC->ProdProcess(proc);
- if((pMC->TrackPid()==443 || pMC->TrackPid()==553) && !strcmp(proc,"DCAY")) {
+ if((nsec=gMC->NSecondaries())>0) {
+ gMC->ProdProcess(proc);
+ if((gMC->TrackPid()==443 || gMC->TrackPid()==553) && !strcmp(proc,"DCAY")) {
//
// Check angular acceptance
// --- and have muons from resonance decays in the wanted window ---
if(nsec != 2) {
printf(" AliMUON::StepManager: Strange resonance Decay into %d particles\n",nsec);
- pMC->StopEvent();
+ gMC->StopEvent();
} else {
- pMC->GetSecondary(0,ipart,x,p);
+ gMC->GetSecondary(0,ipart,x,p);
pt = TMath::Sqrt(p[0]*p[0]+p[1]*p[1]);
th0 = TMath::ATan2(pt,p[2])*kRaddeg;
- pMC->GetSecondary(1,ipart,x,p);
+ gMC->GetSecondary(1,ipart,x,p);
pt = TMath::Sqrt(p[0]*p[0]+p[1]*p[1]);
th2 = TMath::ATan2(pt,p[2])*kRaddeg;
if(!(fAccMin < th0 && th0 < fAccMax) ||
!(fAccMin < th2 && th2 < fAccMax))
- pMC->StopEvent();
+ gMC->StopEvent();
}
}
}
File=new TFile(fFileName);
cout<<"I have opened "<<fFileName<<" file "<<endl;
fHits2 = new TClonesArray("AliMUONhit",1000 );
- fParticles2 = new TClonesArray("GParticle",1000);
+ fParticles2 = new TClonesArray("TParticle",1000);
first=kFALSE;
}
File->cd();
}
// TClonesArray &Partarray = *(gAlice->Particles());
TClonesArray *PartArray = gAlice->Particles();
- GParticle *Part;
+ TParticle *Part;
// Import the Kine and Hits Trees for the event evNumber in the file
// Di-Muon Loop
Float_t pt1,pt2;
- GParticle* Muon1, *Muon2;
+ TParticle* Muon1, *Muon2;
/*
Combinator->ResetRange();
for (Combinator->FirstMuonPairSelected(Muon1,Muon2);
* $Id$
*
* $Log$
+* Revision 1.3 1999/10/05 17:15:45 fca
+* Minor syntax for the Alpha OSF
+*
* Revision 1.2 1999/08/06 14:12:30 fca
* Remove several warnings
*
$,XMBEG,YMBEG,ZMBEG,XMEND,YMEND,ZMEND
$,BV(MAXFLD)
+ CHARACTER*255 CHDIR
+
ISXFMAP = 3
IF(ISXFMAP.EQ.1) THEN
10000 FORMAT(' *SXFMAP* Magnetic field map flag: ',I5
$ ,'; Reading magnetic field map data ')
*
- OPEN(77,FILE='data/field01.dat',
+ CHDIR=' '
+ CALL GETENVF('ALICE_ROOT',CHDIR)
+ LEND=LNBLNK(CHDIR)
+ OPEN(77,FILE=CHDIR(1:LEND)//'/data/field01.dat',
$ FORM='FORMATTED',STATUS='OLD')
READ(77,*) NX,NY,NZ,DX,DY,DZ,XMBEG,YMBEG,ZMBEG
PRINT*,'NX,NY,NZ,DX,DY,DZ,XMBEG,YMBEG,ZMBEG',