StRoot: macros/examples/geant.C Source File

00001 // $Id: geant.C,v 1.14 2001/05/30 18:24:59 perev Exp $
00002 //=======================================================================
00003 // owner: Pavel Nevski
00004 // what it does: 
00005 //=======================================================================
00006 //
00007 //#define gtrack
00008 TBrowser *b = 0;
00009 class St_geant_Maker;
00010 St_geant_Maker *geant=0;
00011 void geant(const Int_t Nevents=1,
00012            const Char_t *fzfile ="/star/rcf/simu/cocktail/hadronic/default/lowdensity/year_1h/half_field/hadronic_on/Gstardata/hc_lowdensity.400_evts.fz")
00013 {
00014   gROOT->LoadMacro("bfc.C");
00015   bfc(0,"fzin sim_T globT gen_T",fzfile);
00016   Int_t i=0;
00017   for (Int_t i =1; i <= Nevents; i++){
00018     chain->Clear();
00019     if (chain->Make(i)>=kStEOF) break;
00020     St_g2t_track *track = (St_g2t_track *) chain->FindObject("g2t_track");
00021      g2t_track_st *trk = track->GetTable();
00022      for (Int_t j = 0; j < track->GetNRows(); j++,trk++) {
00023        if (TMath::Abs(trk->charge) < 0.5) continue;
00024        if (! trk->eg_label) continue;
00025        if (trk->n_tpc_hit < 26) continue;
00026        hyp = -1;
00027        if (trk->ge_pid == 2 || trk->ge_pid == 3) hyp = 3;
00028        if (trk->ge_pid == 8 || trk->ge_pid == 9) hyp = 2;
00029        if (trk->ge_pid ==11 || trk->ge_pid ==12) hyp = 1;
00030        if (trk->ge_pid ==14 || trk->ge_pid ==15) hyp = 0;
00031        Double_t pT =  trk->charge*trk->pt;
00032        Double_t Eta = trk->eta;
00033        printf("eg: %d ge: %d pT %f eta %f  Ntpc %d\n",
00034        trk->eg_label,trk->ge_pid,pT,Eta,trk->n_tpc_hit);
00035      }
00036     printf ("=========================================== Done with Event no. %d\n",i);
00037   }
00038 }