AlgShieldPlankList.cxx

Go to the documentation of this file.

//
// $Id: AlgShieldPlankList.cxx,v 1.3 2007/02/04 05:26:51 rhatcher Exp $ 
//
// AlgShieldPlankList
//
// Package: ShieldPlank
//
// Super class for shield strips
//
// Contact: Caius Howcroft<howcroft@hep.phy.cam.ac.uk>
//
//

#include <cmath>
#include <cassert>

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "Algorithm/AlgConfig.h"
#include "Candidate/CandContext.h"
#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitListHandle.h"

#include "ShieldPlank/AlgShieldPlankList.h"
#include "ShieldPlank/CandShieldPlankList.h"
#include "ShieldPlank/CandShieldPlankListHandle.h"

#include "ShieldPlank/CandShieldPlank.h"
#include "ShieldPlank/CandShieldPlankList.h"

#include "Conventions/Mphysical.h"
#include "Conventions/Munits.h"
#include "Conventions/PlaneView.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "Navigation/NavKey.h"
#include "Navigation/NavSet.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "Validity/VldContext.h"

static NavKey StripSRKeyFromPSEId(const CandDigitHandle *);
 
CVSID("$Id: AlgShieldPlankList.cxx,v 1.3 2007/02/04 05:26:51 rhatcher Exp $");
//______________________________________________________________________
AlgShieldPlankList::AlgShieldPlankList()
{
}

//______________________________________________________________________
AlgShieldPlankList::~AlgShieldPlankList()
{
}

//______________________________________________________________________
void AlgShieldPlankList::RunAlg(AlgConfig& /* ac */, 
                                CandHandle& ch, CandContext& cx)
{

  MSG("ShieldPlank",Msg::kDebug) << "Start ALGSHIELDSTRIPLIST"<<endl;
  assert(cx.GetDataIn());

  if (!(cx.GetDataIn()->InheritsFrom("CandDigitListHandle"))) {
    MSG("ShieldPlank",Msg::kError) << "AlgShieldPlankList wasnt passed a CandDigitList"<<endl;
    return;
  }

  //======================================================
  //Set up all pointer junk
  //
  CandShieldPlankListHandle& csslh = (CandShieldPlankListHandle&) ch;
  CandRecord *candrec = cx.GetCandRecord();
  assert(candrec);
  const VldContext *vldcptr = candrec->GetVldContext();
  assert(vldcptr);
  VldContext vldc = *vldcptr;
  

  CandContext cxx(this,cx.GetMom());
  cxx.SetCandRecord(candrec);

  AlgFactory &af = AlgFactory::GetInstance();
  AlgHandle ah = af.GetAlgHandle("AlgShieldPlank","default");
  
  
  const CandDigitListHandle *cdlh =
    dynamic_cast<const CandDigitListHandle*>(cx.GetDataIn());
     


  //=====================================================
  // Sort by Plane and Strip
  //
  MSG("ShieldPlank",Msg::kDebug) << "Sort by strip and plane:"<<cdlh->GetNDaughters()<<endl;
  CandDigitHandleItr cdhItr(cdlh->GetDaughterIterator());
  CandDigitHandleKeyFunc *cdhKf = cdhItr.CreateKeyFunc();
  cdhKf->SetFun(StripSRKeyFromPSEId);
  cdhItr.GetSet()->AdoptSortKeyFunc(cdhKf);
  cdhKf = 0;
  MSG("ShieldPlank",Msg::kDebug) << "Done sort"<<endl;
  
  
  //=====================================================
  //
  //
  CandDigitHandle *prev = 0;
  CandDigitHandle *first = 0;
  TObjArray cdhAr;
  cdhAr.Clear();
  while (CandDigitHandle *curr = cdhItr()) {
    PlexSEIdAltL pseid = curr->GetPlexSEIdAltL();

    if(!pseid.IsVetoShield())continue;
    MSG("ShieldPlank",Msg::kDebug) << "considering digit, plane " << pseid.GetPlane() << " plexsize " << pseid.GetSize() << endl;
    if (prev) {
      bool issameplank = 0;
      const PlexSEIdAltL& pSEIdAltL = prev->GetPlexSEIdAltL();
      const PlexSEIdAltL& cSEIdAltL = curr->GetPlexSEIdAltL();
      for(PlexSEIdAltL::const_iterator iterp = pSEIdAltL.begin();
          iterp!=pSEIdAltL.end(); ++iterp)
        {
           for(PlexSEIdAltL::const_iterator iterc = cSEIdAltL.begin();
               iterc!=cSEIdAltL.end(); ++iterc)
             {
               if((*iterp).GetSEId().IsSameStrip( (*iterc).GetSEId() ))
                 {
                   //found match
                   MSG("ShieldPlank",Msg::kDebug) << "MATCH " <<(*iterp).GetSEId().GetPlane() << "=="<<(*iterc).GetSEId().GetPlane()<< " Strip:"<<  (*iterp).GetSEId().GetStrip()<<" == "<<(*iterc).GetSEId().GetStrip()<<endl;            
                   issameplank = 1;
                   break;
                 }
             }
           if(issameplank) break;
        }

      
      if (issameplank) {
        cdhAr.Add(curr);
      } else {
        MSG("ShieldPlank",Msg::kDebug) << "Making shield strip from "<< cdhAr.GetLast()<<" digits"<<endl; 
        cxx.SetDataIn(&cdhAr);

        CandShieldPlankHandle csh = CandShieldPlank::MakeCandidate(ah,cxx);
        csslh.AddDaughterLink(csh);
        cdhAr.Clear();
        cdhAr.Add(curr);
        first = curr;
      }
    } else {
      cdhAr.Add(curr);
      first = curr;
    }
    prev = curr;
  }
  
  if (prev) {
    cxx.SetDataIn(&cdhAr);
    MSG("ShieldPlank",Msg::kDebug) << "Making shield strip from "<< cdhAr.GetLast()<<" digits"<<endl; 
    CandShieldPlankHandle csh = CandShieldPlank::MakeCandidate(ah,cxx);
    csslh.AddDaughterLink(csh);
  }

  //
  //Set the Abs time of shield plank list to the abs time of the digit list
  //
  double abstime(cdlh->GetAbsTime());
  if(!finite(abstime)) abstime = -99e99;
  csslh.SetAbsTime(abstime);
  
}

NavKey StripSRKeyFromPSEId(const CandDigitHandle *cdh)
{
  //junk function stolen from roy
  return const_cast<CandDigitHandle *>(cdh)->GetPlexSEIdAltL().GetBestSEId().GetEncoded();
}

//______________________________________________________________________
void AlgShieldPlankList::Trace(const char* /*c*/) const
{
  //I have no idea what this does! more pointless code?
}

ClassImp(AlgShieldPlankList)

---