kine_tracks.C
authormtadel <mtadel@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 Jul 2008 15:42:46 +0000 (15:42 +0000)
committermtadel <mtadel@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 Jul 2008 15:42:46 +0000 (15:42 +0000)
  Set default cut on particle pT, p to zero.

esd_tracks.C
  Set default cut on particle pT, p to zero.
  Increase max-radius of track-propagator to 520cm (to reach TOF,
  HMPID and PHOS).

visscan_init.C
  Add TOF clusters.

its_*.C
  Temporary workaround for a cint problem in root-5.20.

EVE/alice-macros/esd_tracks.C
EVE/alice-macros/its_digits.C
EVE/alice-macros/its_module_stepper.C
EVE/alice-macros/its_raw.C
EVE/alice-macros/kine_tracks.C
EVE/alice-macros/tof_clusters.C
EVE/macros/visscan_init.C

index 828c2c6..46d6641 100644 (file)
@@ -7,7 +7,7 @@
  * full copyright notice.                                                 *
  **************************************************************************/
 
-TEveTrack* esd_make_track(TEveTrackPropagator*   rnrStyle,
+TEveTrack* esd_make_track(TEveTrackPropagator*   trkProp,
                          Int_t                  index,
                          AliESDtrack*           at,
                          AliExternalTrackParam* tp=0)
@@ -29,7 +29,7 @@ TEveTrack* esd_make_track(TEveTrackPropagator*   rnrStyle,
   Double_t ep = at->GetP(), mc = at->GetMass();
   rt.fBeta = ep/TMath::Sqrt(ep*ep + mc*mc);
 
-  TEveTrack* track = new TEveTrack(&rt, rnrStyle);
+  TEveTrack* track = new TEveTrack(&rt, trkProp);
   //PH The line below is replaced waiting for a fix in Root
   //PH which permits to use variable siza arguments in CINT
   //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
@@ -65,7 +65,7 @@ Bool_t gkFixFailedITSExtr    = kFALSE;
 // Also show lines as generated by AliESDtrack.
 Bool_t gkMakeTrackParamLines = kFALSE;
 
-TEveTrackList* esd_tracks(Double_t min_pt=0.1, Double_t max_pt=100)
+TEveTrackList* esd_tracks(Double_t min_pt=0, Double_t max_pt=10000)
 {
   AliESDEvent* esd = AliEveEventManager::AssertESD();
 
@@ -75,8 +75,9 @@ TEveTrackList* esd_tracks(Double_t min_pt=0.1, Double_t max_pt=100)
 
   TEveTrackList* cont = new TEveTrackList("ESD Tracks");
   cont->SetMainColor(6);
-  TEveTrackPropagator* rnrStyle = cont->GetPropagator();
-  rnrStyle->SetMagField( 0.1*esd->GetMagneticField() );
+  TEveTrackPropagator* trkProp = cont->GetPropagator();
+  trkProp->SetMagField( 0.1*esd->GetMagneticField() );
+  trkProp->SetMaxR    ( 520 );
 
   gEve->AddElement(cont);
 
@@ -110,7 +111,7 @@ TEveTrackList* esd_tracks(Double_t min_pt=0.1, Double_t max_pt=100)
       innerTaken = kTRUE;
     }
 
-    TEveTrack* track = esd_make_track(rnrStyle, n, at, tp);
+    TEveTrack* track = esd_make_track(trkProp, n, at, tp);
     track->SetAttLineAttMarker(cont);
 
     if (!innerTaken) {
@@ -157,8 +158,9 @@ TEveTrackList* esd_tracks_from_array(TCollection* col, AliESDEvent* esd=0)
 
   TEveTrackList* cont = new TEveTrackList("ESD Tracks");
   cont->SetMainColor(6);
-  TEveTrackPropagator* rnrStyle = cont->GetPropagator();
-  rnrStyle->SetMagField( 0.1*esd->GetMagneticField() );
+  TEveTrackPropagator* trkProp = cont->GetPropagator();
+  trkProp->SetMagField( 0.1*esd->GetMagneticField() );
+  trkProp->SetMaxR    ( 520 );
 
   gEve->AddElement(cont);
 
@@ -176,7 +178,7 @@ TEveTrackList* esd_tracks_from_array(TCollection* col, AliESDEvent* esd=0)
     ++count;
     AliESDtrack* at = (AliESDtrack*) obj;
 
-    TEveTrack* track = esd_make_track(rnrStyle, count, at);
+    TEveTrack* track = esd_make_track(trkProp, count, at);
     track->SetAttLineAttMarker(cont);
     gEve->AddElement(track, cont);
   }
@@ -269,30 +271,35 @@ TEveElementList* esd_tracks_vertex_cut()
   tl[0] = new TEveTrackList("Sigma < 3");
   tc[0] = 0;
   tl[0]->GetPropagator()->SetMagField( 0.1*esd->GetMagneticField() );
+  tl[0]->GetPropagator()->SetMaxR    ( 520 );
   tl[0]->SetMainColor(3);
   gEve->AddElement(tl[0], cont);
 
   tl[1] = new TEveTrackList("3 < Sigma < 5");
   tc[1] = 0;
   tl[1]->GetPropagator()->SetMagField( 0.1*esd->GetMagneticField() );
+  tl[1]->GetPropagator()->SetMaxR    ( 520 );
   tl[1]->SetMainColor(7);
   gEve->AddElement(tl[1], cont);
 
   tl[2] = new TEveTrackList("5 < Sigma");
   tc[2] = 0;
   tl[2]->GetPropagator()->SetMagField( 0.1*esd->GetMagneticField() );
+  tl[2]->GetPropagator()->SetMaxR    ( 520 );
   tl[2]->SetMainColor(46);
   gEve->AddElement(tl[2], cont);
 
   tl[3] = new TEveTrackList("no ITS refit; Sigma < 5");
   tc[3] = 0;
   tl[3]->GetPropagator()->SetMagField( 0.1*esd->GetMagneticField() );
+  tl[3]->GetPropagator()->SetMaxR    ( 520 );
   tl[3]->SetMainColor(41);
   gEve->AddElement(tl[3], cont);
 
   tl[4] = new TEveTrackList("no ITS refit; Sigma > 5");
   tc[4] = 0;
   tl[4]->GetPropagator()->SetMagField( 0.1*esd->GetMagneticField() );
+  tl[4]->GetPropagator()->SetMaxR    ( 520 );
   tl[4]->SetMainColor(48);
   gEve->AddElement(tl[4], cont);
 
index 1aa1b99..da08664 100644 (file)
@@ -7,6 +7,8 @@
  * full copyright notice.                                                 *
  **************************************************************************/
 
+#include <vector> // bypass a cint problem in root-5.20
+
 // Load ITS digits.
 // Argument mode is a bitwise or determining which layers to import:
 //    1,  2 : SPD
index 34e7004..56d5f9d 100644 (file)
@@ -7,6 +7,8 @@
  * full copyright notice.                                                 *
  **************************************************************************/
 
+#include <vector> // bypass a cint problem in root-5.20
+
 class AliEveITSModuleStepper;
 
 void its_module_stepper(Int_t det = 0)
index 75c8f54..34c6b9e 100644 (file)
@@ -7,6 +7,8 @@
  * full copyright notice.                                                 *
  **************************************************************************/
 
+#include <vector> // bypass a cint problem in root-5.20
+
 // Load ITS raw-data.
 // Argument mode is a bitwise or determining which layers to import:
 //    1,  2 : SPD
index 173a193..d8ddf0d 100644 (file)
@@ -13,7 +13,7 @@
 #include "TParticlePDG.h"
 
 TEveTrackList*
-kine_tracks(Double_t min_pt  = 0.1,   Double_t min_p   = 0.2,
+kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
            Bool_t   pdg_col = kTRUE, Bool_t   recurse = kTRUE,
            Bool_t   use_track_refs = kTRUE)
 {
@@ -38,7 +38,7 @@ kine_tracks(Double_t min_pt  = 0.1,   Double_t min_p   = 0.2,
   Int_t N = stack->GetNtrack();
   for (Int_t i=0; i<N; ++i)
   {
-    if(stack->IsPhysicalPrimary(i))
+    if (stack->IsPhysicalPrimary(i))
     {
       TParticle* p = stack->Particle(i);
       if (p->Pt() < min_pt && p->P() < min_p) continue;
index 9389ff5..adedb1b 100644 (file)
@@ -45,7 +45,6 @@ TEvePointSet* tof_clusters(TEveElement* cont=0, Float_t maxR=390)
     if (!cTree->GetEvent(i)) continue;
 
     Int_t ncl=cl->GetEntriesFast();
-    cout<<" ncl = "<<ncl<<endl;
     while (ncl--) {
       AliCluster *c=(AliCluster*)cl->UncheckedAt(ncl);
       Float_t g[3]; //global coordinates
@@ -83,8 +82,8 @@ TEvePointSet* tof_clusters(TEveElement* cont=0, Float_t maxR=390)
   return clusters;
 }
 
-TEvePointSet* tof_clusters(Int_t selectedSector=-1,
-                          TEveElement* cont=0, Float_t maxR=390)
+TEvePointSet* tof_clusters_sec(Int_t selectedSector,
+                               TEveElement* cont=0, Float_t maxR=390)
 {
   AliEveEventManager::AssertGeometry();
 
@@ -109,7 +108,6 @@ TEvePointSet* tof_clusters(Int_t selectedSector=-1,
     if (!cTree->GetEvent(i)) continue;
 
     Int_t ncl=cl->GetEntriesFast();
-    cout<<" ncl = "<<ncl<<endl;
     while (ncl--) {
       AliCluster *c=(AliCluster*)cl->UncheckedAt(ncl);
       Float_t g[3]; //global coordinates
index 97dc164..4faf49e 100644 (file)
@@ -34,6 +34,7 @@ void visscan_init()
   TEveUtil::LoadMacro("its_clusters.C+");
   TEveUtil::LoadMacro("tpc_clusters.C+");
   TEveUtil::LoadMacro("trd_clusters.C+");
+  TEveUtil::LoadMacro("tof_clusters.C+");
 
   TEveLine::SetDefaultSmooth(1);
 
@@ -60,7 +61,7 @@ void visscan_init()
   if (gShowTRD) gGeomGentleTRD = geom_gentle_trd();
 
 
-  gROOT->ProcessLine(".L SplitGLView.C+g"); // !!!! debug-mode
+  gROOT->ProcessLine(".L SplitGLView.C+");
   TEveBrowser* browser = gEve->GetBrowser();
   browser->ExecPlugin("SplitGLView", 0, "new SplitGLView(gClient->GetRoot(), 600, 450, kTRUE)");
 
@@ -109,6 +110,13 @@ void on_new_event()
       trdc->SetMarkerStyle(4);
       trdc->SetMarkerSize(0.5);
     }
+
+    TEvePointSet* tofc = tof_clusters();
+    if (tofc) {
+      tofc->SetMarkerColor(kOrange);
+      tofc->SetMarkerStyle(4);
+      tofc->SetMarkerSize(0.5);
+    }
   }
   catch(TEveException& exc) {
     printf("Exception loading ITS/TPC clusters: %s\n", exc.Data());