EVE/alice-macros/tpc_hits_charge_split.C

   1 // $Id$
   2
   3 void tpc_hits_charge_split(const char *varexp    =
   4                         "TPC2.fArray.fR:TPC2.fArray.fFi:TPC2.fArray.fZ"
   5                         ":log(TPC2.fArray.fCharge)",
   6                         const char *selection = "TPC2.fArray.fR>80")
   7 {
   8   // Extracts 'major' TPC hits (not the compressed ones).
   9   // This gives ~2.5% of all hits.
  10
  11   AliRunLoader* rl =  Alieve::Event::AssertRunLoader();
  12   rl->LoadHits("TPC");
  13
  14   TTree* ht = rl->GetTreeH("TPC", false);
  15
  16   Reve::PointSetArray* l = new Reve::PointSetArray("TPC hits - Log-Charge Slices", "");
  17   l->SetSourceCS(TPointSelectorConsumer::TVT_RPhiZ);
  18   l->SetMarkerColor((Color_t)3);
  19   l->SetMarkerStyle(20); // full circle
  20   l->SetMarkerSize(.5);
  21
  22   gReve->AddRenderElement(l);
  23   l->InitBins("Log Charge", 20, 0, 5);
  24
  25   TPointSelector ps(ht, l, varexp, selection);
  26   ps.Select();
  27
  28   l->CloseBins();
  29
  30   gReve->Redraw3D();
  31 }