- Int_t nparticles = gAlice->GetEvent(evnum);
- if (nparticles <= 0) return;
- printf("\n nparticles %d\n",nparticles);
-
- TClonesArray *particles = gAlice->Particles();
-
- sprintf(nameTH,"TreeH%d",evnum);
- printf("%s\n",nameTH);
- TTree *th = gAlice->TreeH();
- Int_t ntracks = (Int_t) th->GetEntries();
- if (ntracks<=0) return;
- // Start loop on tracks in the hits containers
- for (Int_t track=0; track<ntracks;track++) {
- gAlice->ResetHits();
- nbytes += th->GetEvent(track);
- particle=(TParticle*)particles->UncheckedAt(track);
- nhits =fHits->GetEntriesFast();
-
- for (hit=0;hit<nhits;hit++) {
- startHit = (AliSTARThit*)fHits->UncheckedAt(hit);
- pmt=startHit->fPmt;
- e=startHit->fEtot;
- x=startHit->X();
- y=startHit->Y();
- volume = startHit->fVolume;
- if(volume==1){
- timeright = startHit->fTime;
- if(timeright<besttimeright) {
- besttimeright=timeright;
- } //timeright
- }//time for right shoulder
- if(volume==2){
- timeleft = startHit->fTime;
- // printf("timeleft %f\n",timeleft);
- if(timeleft<besttimeleft) {
- besttimeleft=timeleft;
- } //timeleftbest
- }//time for left shoulder
- } //hit loop
- } //track loop