hooks for PMD flow analysis
[u/mrichter/AliRoot.git] / 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 }