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->Load("libPWG2flowCommon");
  11
  12   Int_t maxNumberOfEvents = 1000;
  13
  14   //define reference particles
  15   AliStarTrackCuts* rpCuts = AliStarTrackCuts::StandardCuts();
  16
  17   //define particles of interest
  18   AliStarTrackCuts* poiCuts = AliStarTrackCuts::StandardCuts();
  19   poiCuts->SetPtMin(1.0);
  20
  21   //define event cuts
  22   AliStarEventCuts* starEventCuts = AliStarEventCuts::StandardCuts();
  23
  24   Int_t i=0;
  25   AliStarEventReader starReader("/data/alice3/jthomas/testData/") ;
  26   while ( starReader.GetNextEvent() )                                // Get next event
  27   {
  28     AliStarEvent* starEvent = starReader.GetEvent();
  29     if ( !starEventCuts->PassesCuts(starEvent) ) continue;              // Test if the event is good
  30
  31     AliFlowEventSimple* flowEvent = new AliFlowEventStar(starEvent,rpCuts,poiCuts);  // make a flow event from a star event (aka "the magic")
  32
  33     /////analysis here////////////////
  34
  35
  36
  37     //////////////////////////////////
  38
  39     //starEvent->Print("all");
  40     flowEvent->Print();
  41
  42     delete flowEvent;
  43
  44     i++;
  45     if (i>maxNumberOfEvents) break;
  46   }
  47   delete rpCuts;
  48   delete poiCuts;
  49   delete starEventCuts;
  50 }