GfxMCVecList.cxx

Go to the documentation of this file.

#include "GfxMCVecList.h"
#include "GfxMCVecListCfg.h"
#include "GfxProxy.h"
#include "ViewState.h"
#include "GfxMCVec.h"
#include <Midad/Base/Mint.h>
#include "GfxMCVecListMenu.h"

#include <MinosObjectMap/MomNavigator.h>
#include <Record/SimSnarlRecord.h>
#include <Record/RecArrayAllocator.h>
#include <TParticle.h>
#include <TClonesArray.h>
#include <DataUtil/TruthHelper.h>
#include <RecoBase/CandSliceHandle.h>
#include <string>
using namespace std;

//#include "<Midad/Base/pdg_kludge.h"

#include <sigc++/sigc++.h>
#include <sigc++/class_slot.h>
using namespace SigC;


#include <Buttons.h> // part of ROOT

static GfxProxy<GfxMCVecList> 
gsGfxMCVecListProxy("MCVecList", new GfxMCVecListMenu());

GfxMCVecList::GfxMCVecList()
    :fCfg(0)
      ,fHideVecs(false)
{
    fCfg = new GfxMCVecListCfg();
    fCfg->modified_signal.connect
        (SigC::slot_class(*this,&GfxMCVecList::ReConfigure));
    //   pdg_kludge();
    RecArrayAllocator& a = RecArrayAllocator::Instance();
    fVecsCA = a.GetArray("GfxMCVec");
} 

GfxMCVecList::~GfxMCVecList()
{
 RecArrayAllocator& a = RecArrayAllocator::Instance();
 if (fVecsCA) { a.ReleaseArray(fVecsCA); fVecsCA = 0; }
}


void GfxMCVecList::Init(PageDisplay& pd, EVD& mp)
{

  fEVD= &mp;
  fPageDisplay= &pd;
  ViewState* vs = this->GetViewState();
  vs->spatial_metric.connect(slot_class(*this,&GfxMCVecList::ReConfigure));
}


void GfxMCVecList::Configure(Mint& mint)
{

    if (fVecs.size()) this->ClearVecList();
    if (fHideVecs) return;
    const MomNavigator* mom = mint.GetJint().GetMom();
    if (!mom) return;
    
    SimSnarlRecord* ssr = 
      dynamic_cast<SimSnarlRecord*>(mom->GetFragment("SimSnarlRecord"));
    if (!ssr) return;
      const TClonesArray* ctca = 
        dynamic_cast<const TClonesArray*>
        (ssr->FindComponent("TClonesArray","StdHep"));
    if (!ctca) return;

    float evntP=0;
    int ind, siz = ctca->GetEntriesFast();
    for (ind=0; ind < siz; ++ind) {
      TParticle* part = dynamic_cast<TParticle*>((*ctca)[ind]);
      if (!part) {
        cerr << "Non particle: " << (*ctca)[ind]->GetName() << endl;
        continue;
      }
     
      if(part->P()>evntP && part->P()<100.0 )evntP=part->P(); 
    }
    const CandSliceHandle * slc = fEVD->GetCandSlice();
    Int_t neuId= -1;
    Int_t nextneuId = -1;
    if(slc){
      CandSliceHandle slch = *slc;
      TruthHelper truthHelp(mom);
      neuId=truthHelp.GetBestSliceNeuMatch(slch);
      nextneuId=truthHelp.GetNextNeuId(neuId);
    }
    CandEventHandle * ev = fEVD->GetCandEvent();
    if(ev){
      TruthHelper truthHelp(mom);
      neuId=truthHelp.GetBestEventNeuMatch(*ev);
      nextneuId=truthHelp.GetNextNeuId(neuId);
    }
    for (ind=0; ind < siz; ++ind) {      
      TParticle* part = dynamic_cast<TParticle*>((*ctca)[ind]);
      if (!part) {
        cerr << "Non particle: " << (*ctca)[ind]->GetName() << endl;
        continue;
      }
      if(neuId<0 || (ind >=neuId && ind < nextneuId)){
        if(abs(part->GetPdgCode())<19 ||       // elementary particles
           abs(part->GetPdgCode())==211 ||    // charged pions
           abs(part->GetPdgCode())==2212      // proton
           ){
          int last = fVecsCA->GetLast()+1;
          new ((*fVecsCA)[last]) GfxMCVec(*part,evntP,*this);
          GfxMCVec* gs = dynamic_cast<GfxMCVec*>(fVecsCA->UncheckedAt(last));
          fVecs.push_back(gs);
        }
      }
    }
    this->ReConfigure();

}

void GfxMCVecList::ReConfigure()
{
  if (fHideVecs) return;
  int ind=0;
    for (MCVecList_t::iterator it=fVecs.begin(); it != fVecs.end(); ++it) {
        (*it)->Configure(ind);
        ind++;
    }
}

void GfxMCVecList::ClearVecList() 
{
    fVecs.clear();
    fVecsCA->Clear("C");
}

void GfxMCVecList::Draw(Option_t* option)
{
//    cerr << "GfxMCVecList::Draw("<<option<<")\n";
    for (MCVecList_t::iterator it=fVecs.begin(); it != fVecs.end(); ++it) {
        (*it)->Draw(option);
        
    }
}

void GfxMCVecList::ExecuteEvent(int event, int /* px */, int /* py */, GfxMCVec* MCVec)
{
    switch (event) {
    case kMouseEnter: {
        text_info.emit(Form("%s: %s  P=%8.6f",
                            "MCVec", MCVec->GetParticle().GetName(),
                            MCVec->GetParticle().P()));
        break;
    }
    default: // nothing
        break;
    }

}



void GfxMCVecList::HideVecs(bool hide)
{
    fHideVecs = hide;
    this->ReConfigure();
    if(GetGraphicsView()==GraphicsView::kXY)fEVD->UpdateCanvas();

}

---