1 //example script on what to do with the star events
3 // root readStarEventSimple.C
5 void readStarEventsSimple()
7 gSystem->Load("libTree.so");
8 gSystem->Load("libVMC.so");
9 gSystem->Load("libPhysics.so");
10 gSystem->Load("libPWGflowBase");
12 Int_t maxNumberOfEvents = 1000;
14 //define reference particles
15 AliStarTrackCuts* rpCuts = AliStarTrackCuts::StandardCuts();
17 //define particles of interest
18 AliStarTrackCuts* poiCuts = AliStarTrackCuts::StandardCuts();
19 poiCuts->SetPtMin(1.0);
22 AliStarEventCuts* starEventCuts = AliStarEventCuts::StandardCuts();
25 AliStarEventReader starReader("/data/alice3/jthomas/testData/") ;
26 while ( starReader.GetNextEvent() ) // Get next event
28 AliStarEvent* starEvent = starReader.GetEvent();
29 if ( !starEventCuts->PassesCuts(starEvent) ) continue; // Test if the event is good
31 AliFlowEventSimple* flowEvent = new AliFlowEventStar(starEvent,rpCuts,poiCuts); // make a flow event from a star event (aka "the magic")
33 /////analysis here////////////////
37 //////////////////////////////////
39 //starEvent->Print("all");
45 if (i>maxNumberOfEvents) break;