macros/graphics/Ed.C

Go to the documentation of this file.

// $Id: Ed.C,v 1.15 2010/04/29 19:16:36 fine Exp $
// *-- Author :    Valery Fine(fine@bnl.gov)   25/02/2009


class StEvent;
StEvent* event = 0;


//__________________________________________
void rd(int  hits=0, bool clear=false) 
{  
   // redraw the event
   if (event) {
      if (clear) gEventDisplay->Clear();
      if (hits) {
        if (hits & 0x1 ) gEventDisplay->Hits(event,1);
        if (hits & 0x2 ) gEventDisplay->EmcHits(event);
        if (hits & 0x4 ) gEventDisplay->FtpcHits(event,1);
 //       if (hits & 4) gEventDisplay->EmcHits(event,"eemc");
      }
      else gEventDisplay->Tracks(event);
      // gEventDisplay->Update();
   }
}
//__________________________________________
void skipe(int nEvents=1) {
   if(chain) chain->Skip(nEvents);
}

//__________________________________________

//__________________________________________
void ae(int tracks=-1, int  hits=-1) 
{
 // Advance till next "good" event
 // One may want to replace the "plain"  "if" clause below
 // with the full-fledged filter
   static int defaultTracks = 0;
   static int defaultHits   = 3;

   if (tracks != -1 ) defaultTracks = tracks;
   else tracks = defaultTracks;

   if (hits != -1 )   defaultHits = hits;
   else  hits = defaultHits;
   gEventDisplay->Clear();
 newevent:
     chain->MakeEvent();
     event = (StEvent*)chain->GetDataSet("StEvent");
     if (event && !event->trackNodes().empty())) {
         if (tracks) rd();     // Draw the tracks
         if (hits)  rd(hits); // Add the hits to the image
    } else {
        printf(" event is empty %p\n", event);
        goto newevent;
     }
 }

//__________________________________________
 void Ed(const char* file =
 "/star/institutions/bnl/fine/testfiles/st_physics_10169042_raw_4030001.event.root"
 , unsigned int nEvent=1, const char * detectorNames="TPC")
 {
   // Start application open the file provided.
   if ( gSystem->AccessPathName(file)) {
      cout << endl << endl 
           << "** Error ** : The input file: <"<< file << "> does not exist !!!" 
           << endl << endl
           <<  " Please select the existing one and re-call the function: " 
           << endl
           << endl << " root [0] Ed(\"new file ROOT file name\")" 
           << endl 
           << endl << "To draw the StEvent components with no detector geometry, use: "
           << endl
           << endl << " root [0] Ed(\"new file ROOT file name\",0)" 
           << endl
           << endl;
      return;
   }
   gROOT->Macro("Load.C"); 
   gSystem->Load("StDetectorDbMaker");
   gROOT->Macro(Form("bfc.C(%d,1,\"doevents\",\"%s\")",nEvent,file));
   delete gEventDisplay; // destroy the built-in display
   new StuDraw3DEvent(detectorNames); // create our own one (with no detector geometry)
//   new StuDraw3DEvent("TPC"); // create our own one (with TPC detector)
   gEventDisplay->SetBkColor(kBlack);
   printf("\n The display is ready!\n");
   printf(" call:\n");   
   printf("\t---\n");
   printf("\tae()\n");
   printf("\t---\n");
   printf("method to see the next event\n");
   ae();
   printf(" call:\n");   
   printf("\t---\n");
   printf("\tae()    - to draw   the default setting\n");
   printf("\tae(1,1) - to change the default and draw the tracks and its tpc hits\n");
   printf("\tae(0,2) - to change the default and draw the Bemc towers\n");
   printf("\t---\n");
   printf("method to see the next event\n");
 }

---