AlgFarDetSliceList.cxx

Go to the documentation of this file.

#include "AlgFarDetSliceList.h"
#include "FarDetSliceHandle.h"
#include "FarDetSliceListHandle.h"

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "Algorithm/AlgConfig.h"

#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"

#include "Candidate/CandContext.h"
#include "RecoBase/CandStripHandle.h"
#include "RecoBase/CandStripListHandle.h"

#include "Fabrication/FabPlnInstallLookup.h"

ClassImp(AlgFarDetSliceList)    

CVSID("$Id: AlgFarDetSliceList.cxx,v 1.2 2008/08/25 13:03:10 blake Exp $");

AlgFarDetSliceList::AlgFarDetSliceList() :
  fSliceList(0)
{
  fSliceList = new TObjArray();
}   

AlgFarDetSliceList::~AlgFarDetSliceList()
{
  if(fSliceList) delete fSliceList;
}
   
void AlgFarDetSliceList::RunAlg(AlgConfig &ac, CandHandle &ch, CandContext &cx)
{
  MSG("FarDetSlice",Msg::kDebug) << " AlgFarDetSliceList::RunAlg(...) " << endl;

  FarDetSliceListHandle& myslicelist = dynamic_cast<FarDetSliceListHandle&>(ch);

  // Unpack AlgConfig
  Int_t pln;
  Int_t fMinPlane = 0;
  Int_t fMaxPlane = 485;
  Int_t fCheckInstall = 0;
  fMinPlane = ac.GetInt("MinPlane");
  fMaxPlane = ac.GetInt("MaxPlane");
  fCheckInstall = ac.GetInt("CheckInstall");

  MSG("FarDetSlice",Msg::kDebug) << " AlgFarDetSliceList Config: " << endl
                                 << "  MinPlane=" << fMinPlane 
                                 << "  MaxPlane=" << fMaxPlane 
                                 << "  CheckInstall=" << fCheckInstall << endl;

  // Check installation status
  if( fCheckInstall ){

    // Find maximum installed plane
    MSG("FarDetSlice",Msg::kDebug) << " Checking Installation Status: " << endl;
    Int_t fMaxInstallPlane = 485, fTmpInstallPlane = 485;
    CandRecord* candrec = (CandRecord*)(cx.GetCandRecord());
    VldContext *vldc = (VldContext*)(candrec->GetVldContext());
    FabPlnInstallLookup fMaxPlaneLookup(*vldc);
    PlexPlaneId fMaxPlaneId = fMaxPlaneLookup.GetLastNormalPlane();
    if(fMaxPlaneId.IsValid()) fTmpInstallPlane = fMaxPlaneId.GetPlane();
    if(fTmpInstallPlane<fMaxInstallPlane) fMaxInstallPlane = fTmpInstallPlane;
    MSG("FarDetSlice",Msg::kDebug) << "  Maximum Installed Plane=" << fMaxInstallPlane << endl;
  
    // Reset maximum plane
    if( fMaxInstallPlane<fMaxPlane ){
      MSG("FarDetSlice",Msg::kDebug) << "   ... resetting MaxPlane to: " << fMaxInstallPlane << endl;
      fMaxPlane=fMaxInstallPlane;
    }

    // Reset minimum plane
    if( fMaxInstallPlane<fMinPlane ){
      MSG("FarDetSlice",Msg::kDebug) << "   ... resetting MinPlane to: " << fMaxInstallPlane << endl;
      fMinPlane=fMaxInstallPlane;
    }
  }

  // Unpack CandContext
  const CandStripListHandle* cslh = dynamic_cast<const CandStripListHandle*>(cx.GetDataIn());
  TIter strpitr(cslh->GetDaughterIterator());
  while(CandStripHandle* csh = dynamic_cast<CandStripHandle*>(strpitr())){
    if(csh){
      pln = csh->GetPlane();
      if( pln>0 && pln<500 
       && pln>fMinPlane && pln<fMaxPlane ){
        fSliceList->Add(csh);
      }
    } 
  }

  // Make CandSlices
  AlgFactory &af = AlgFactory::GetInstance();
  AlgHandle ahh = af.GetAlgHandle("AlgFarDetSlice","default");
  CandContext cxx(this,cx.GetMom());
  cxx.SetCandRecord(cx.GetCandRecord());

  cxx.SetDataIn(fSliceList);
  FarDetSliceHandle myslice = FarDetSlice::MakeCandidate(ahh,cxx);
  myslice.SetName("FarDetSlice");
  myslice.SetTitle(TString("Created by AlgFarDetSliceList"));
  myslicelist.AddDaughterLink(myslice);
  
  fSliceList->Clear();

}
   
void AlgFarDetSliceList::Trace(const char* /*c*/ ) const
{

}

---