PWG2/FLOW/macros/readStarEventsSimple.C

   1 //example script on what to do with the star events
   2 //run e.g. like this:
   3 //                    root readStarEventSimple.C
   4
   5 void  readStarEventsSimple()
   6 {
   7   gSystem->Load("libTree.so");
   8   gSystem->Load("libVMC.so");
   9   gSystem->Load("libPhysics.so");
  10   gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  11   gSystem->AddIncludePath("-I$ALICE_ROOT/PWG2/FLOW/AliFlowCommon");
  12   gSystem->Load("libSTEERBase");
  13   gSystem->Load("libESD");
  14   gSystem->Load("libAOD");
  15   gSystem->Load("libANALYSIS");
  16   gSystem->Load("libANALYSISalice");
  17   gSystem->Load("libCORRFW");
  18   gSystem->Load("libPWG2flowCommon");
  19   gSystem->Load("libPWG2flowTasks");
  20   gROOT->LoadMacro("AliStarTrack.cxx++");
  21   gROOT->LoadMacro("AliStarEvent.cxx++");
  22   gROOT->LoadMacro("AliStarEventReader.cxx++");
  23   gROOT->LoadMacro("AliFlowEventStar.cxx++");
  24
  25   AliStarEventReader*  starReader = new AliStarEventReader( "/data/alice3/jthomas/testData/") ;
  26
  27   while ( starReader->GetNextEvent() )                                // Get next event
  28   {
  29     AliStarEvent* starEvent = starReader->GetEvent();
  30     if ( !starReader->AcceptEvent(starEvent) ) continue;              // Test if the event is good
  31
  32     AliFlowEventSimple* flowEvent = new AliFlowEventStar(starEvent);  // make a flow event from a star event (aka "the magic")
  33
  34     /////analysis here////////////////
  35
  36
  37
  38     //////////////////////////////////
  39
  40     starEvent->Print("all");
  41     flowEvent->Print("all");
  42
  43     delete flowEvent;
  44     break;
  45   }
  46
  47   delete starReader ;
  48 }