HLT/exa/runtracker.C

   1 /**
   2    Run this macro for cluster finder and track follower (see steering class
   3    AliLevel3.
   4 */
   5
   6 void runtracker(int minslice,int maxslice,char* path=0,char *rootfile=0,int nevent=1)
   7 {
   8   AliL3Logger l;
   9   l.Set(AliL3Logger::kAll);
  10   l.UseStdout();
  11   //l.UseStream();
  12
  13   char path_to_use[1024];
  14   if(!path)
  15     strcpy(path_to_use,"/usr/local/anders/data/hg_42105_s1-3/");
  16   //strcpy(path_to_use,"/prog/alice/data/Rawdata/6_patch/hg_1000_s1-3/");
  17   //strcpy(path_to_use,"/prog/alice/data/Rawdata/6_patch/1track_s1/");
  18   //strcpy(path_to_use,"/prog/alice/data/Rawdata/6_patch/pp/event_0/");
  19   //strcpy(path_to_use,"/prog/alice/data/Rawdata/1_patch/pp/pileups/event_1/");
  20   else strcpy(path_to_use,path);
  21
  22   //AliL3Transform::Init(path_to_use);
  23
  24   Int_t phi_segments,eta_segments,trackletlength,tracklength;
  25   Int_t rowscopetracklet,rowscopetrack;
  26   Double_t min_pt_fit,maxangle,goodDist,hitChi2Cut;
  27   Double_t goodHitChi2,trackChi2Cut,maxphi,maxeta;
  28
  29   for(Int_t ev=0; ev<nevent; ev++)
  30     {
  31       if(!rootfile)
  32         a = new AliLevel3();
  33       else
  34         a = new AliLevel3(rootfile);
  35
  36       a->Init(path_to_use);
  37
  38       phi_segments = 50;
  39       eta_segments = 100;
  40       trackletlength = 3;
  41       tracklength = 5;
  42       rowscopetracklet = 2;
  43       rowscopetrack = 2;
  44       min_pt_fit = 0;
  45       maxangle = 1.31;
  46       goodDist = 5;
  47       maxphi=100;
  48       maxeta=100;
  49       hitChi2Cut = 100;
  50       goodHitChi2 = 20;
  51       trackChi2Cut = 50;
  52
  53       //main vertex tracking parameters:
  54       a->SetTrackerParam(phi_segments,eta_segments,trackletlength,tracklength,
  55                          rowscopetracklet,rowscopetrack,
  56                          min_pt_fit,maxangle,goodDist,hitChi2Cut,
  57                          goodHitChi2,trackChi2Cut,50,maxphi,maxeta,kTRUE);
  58
  59       //a->DoRoi();    /*do region of interest*/
  60       //a->DoMc();     /*do monte carlo identification*/
  61       a->WriteFiles(); /*enable output*/
  62
  63       a->ProcessEvent(minslice,maxslice);
  64       //a->DoBench("benchmark_0");
  65
  66       delete a;
  67     }
  68 }
  69
  70