AlgFarDetShieldPlankList.cxx

Go to the documentation of this file.

#include "AlgFarDetShieldPlankList.h"
#include "FarDetShieldPlankListHandle.h"
#include "FarDetShieldPlankHandle.h"

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

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

#include "CandData/CandRecord.h"
#include "Candidate/CandContext.h"
#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitListHandle.h"

ClassImp(AlgFarDetShieldPlankList)
 
CVSID("$Id: AlgFarDetShieldPlankList.cxx,v 1.5 2007/12/19 21:53:23 bspeak Exp $");

AlgFarDetShieldPlankList::AlgFarDetShieldPlankList() :
  fPairList(0)
{
  fPairList = new TObjArray();
}

AlgFarDetShieldPlankList::~AlgFarDetShieldPlankList()
{
  if(fPairList) delete fPairList;
}

void AlgFarDetShieldPlankList::RunAlg(AlgConfig& /* ac */, 
                                      CandHandle& ch, CandContext& cx)
{
  MSG("FarDetShieldPlankList",Msg::kDebug) << " *** AlgFarDetShieldPlankList::RunAlg(...) *** "<<endl;

  FarDetShieldPlankListHandle& myplanklist = dynamic_cast<FarDetShieldPlankListHandle&>(ch);

  Int_t DigitsBefore=0,DigitsAfter=0;

  Int_t i,j,kn,kp;
  Int_t veto,pln,match;
  TObjArray northlist,southlist;
  TObjArray tmpn,tmpp,tmp1;

  // Unpack CandContext
  CandRecord* candrec = (CandRecord*)(cx.GetCandRecord());
  const TObjArray* inputarray = dynamic_cast<const TObjArray*>(cx.GetDataIn());
  const CandDigitListHandle *cdlh = dynamic_cast<const CandDigitListHandle*>(inputarray->At(0));

  TIter digitr(cdlh->GetDaughterIterator());
  while(CandDigitHandle* cdh = dynamic_cast<CandDigitHandle*>(digitr())){
    if(cdh){
      veto = cdh->GetPlexSEIdAltL().IsVetoShield();

      //Find the lowest Plane value in the AltList to set as the key
      pln = cdh->GetPlexSEIdAltL().GetPlane();
      cdh->GetPlexSEIdAltL().SetFirst();
      while ( cdh->GetPlexSEIdAltL().IsValid() ) {
        if( cdh->GetPlexSEIdAltL().GetCurrentSEId().GetPlane() < pln )
          pln = cdh->GetPlexSEIdAltL().GetCurrentSEId().GetPlane();
        cdh->GetPlexSEIdAltL().Next();
      }
      cdh->GetPlexSEIdAltL().SetFirst();
      pln -= 500;

      if( veto && pln>0 && pln<500 ){
        fDigiList[pln].Add(cdh);
        MSG("FarDetShieldPlank",Msg::kVerbose) << " digit : " 
                                               << " " << cdh->GetPlexSEIdAltL().GetPlane() 
                                               << " " << cdh->GetPlexSEIdAltL().GetEnd()
                                               << " " << cdh->GetCharge(CalDigitType::kNone) 
                                               << " " << 1.0e9*cdh->GetSubtractedTime(CalTimeType::kNone) << endl;
        DigitsBefore++;
      }
    }
  }
    
  // Pair up CandDigits
  for(i=0;i<500;i++){
    if(1+fDigiList[i].GetLast()>0){
    
      for(j=0;j<1+fDigiList[i].GetLast();j++){
        CandDigitHandle* digit = (CandDigitHandle*)(fDigiList[i].At(j));
        const PlexSEIdAltL& pSEIdAltL = digit->GetPlexSEIdAltL();
        if(pSEIdAltL.GetEnd()==StripEnd::kNegative) northlist.Add(digit);
        if(pSEIdAltL.GetEnd()==StripEnd::kPositive) southlist.Add(digit);
      }

      for(kn=0;kn<1+northlist.GetLast();kn++){
        CandDigitHandle* dign = (CandDigitHandle*)(northlist.At(kn));
        const PlexSEIdAltL& nSEIdAltL = dign->GetPlexSEIdAltL();
        match=0;
        for(kp=0;kp<1+southlist.GetLast();kp++){
          CandDigitHandle* digp = (CandDigitHandle*)(southlist.At(kp));
          const PlexSEIdAltL& pSEIdAltL = digp->GetPlexSEIdAltL();
          if( 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))>-500.0  
           && 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))<500.0 ){
            for(PlexSEIdAltL::const_iterator itern = nSEIdAltL.begin();itern!=nSEIdAltL.end(); ++itern){
              for(PlexSEIdAltL::const_iterator iterp = pSEIdAltL.begin();iterp!=pSEIdAltL.end(); ++iterp){
                if((*itern).GetSEId().IsSameStrip((*iterp).GetSEId())) match=1;
              }
            }
          }
        }
        if(match) tmpn.Add(dign); else tmp1.Add(dign);
      }

      for(kp=0;kp<1+southlist.GetLast();kp++){
        CandDigitHandle* digp = (CandDigitHandle*)(southlist.At(kp));
        const PlexSEIdAltL& pSEIdAltL = digp->GetPlexSEIdAltL();
        match=0;
        for(kn=0;kn<1+northlist.GetLast();kn++){
          CandDigitHandle* dign = (CandDigitHandle*)(northlist.At(kn));
          const PlexSEIdAltL& nSEIdAltL = dign->GetPlexSEIdAltL();
          if( 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))>-500.0  
           && 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))<500.0 ){
            for(PlexSEIdAltL::const_iterator iterp = pSEIdAltL.begin();iterp!=pSEIdAltL.end(); ++iterp){
              for(PlexSEIdAltL::const_iterator itern = nSEIdAltL.begin();itern!=nSEIdAltL.end(); ++itern){
                if((*iterp).GetSEId().IsSameStrip((*itern).GetSEId())) match=1;
              }
            }
          }
        }
        if(match) tmpp.Add(digp); else tmp1.Add(digp);
      }

      if( 1+tmpn.GetLast()>0 && 1+tmpp.GetLast()>0 ){
        for(kn=0;kn<1+tmpn.GetLast();kn++){
          CandDigitHandle* dign = (CandDigitHandle*)(tmpn.At(kn));
          const PlexSEIdAltL& nSEIdAltL = dign->GetPlexSEIdAltL();   
          for(kp=0;kp<1+tmpp.GetLast();kp++){
            CandDigitHandle* digp = (CandDigitHandle*)(tmpp.At(kp));
            const PlexSEIdAltL& pSEIdAltL = digp->GetPlexSEIdAltL();
            if( 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))>-500.0  
             && 1.0e9*(dign->GetTime(CalTimeType::kNone)-digp->GetTime(CalTimeType::kNone))<500.0 ){
              match=0;
              for(PlexSEIdAltL::const_iterator itern = nSEIdAltL.begin();itern!=nSEIdAltL.end(); ++itern){
                for(PlexSEIdAltL::const_iterator iterp = pSEIdAltL.begin();iterp!=pSEIdAltL.end(); ++iterp){
                  if((*itern).GetSEId().IsSameStrip((*iterp).GetSEId())) match=1;
                }
              }
              if(match){
                TObjArray* myarray = new TObjArray();
                myarray->Add(dign); DigitsAfter++; 
                myarray->Add(digp); DigitsAfter++;
                fPairList->Add(myarray);
              }
            }
          }
        }
      }

      for(j=0;j<1+tmp1.GetLast();j++){
        CandDigitHandle* dig = (CandDigitHandle*)(tmp1.At(j));
        TObjArray* myarray = new TObjArray();
        myarray->Add(dig); DigitsAfter++; 
        fPairList->Add(myarray);
      }
      
      tmpn.Clear(); tmpp.Clear(); tmp1.Clear(); 
      northlist.Clear(); southlist.Clear();
    }
  }

  // Print out list digit pairs
  for(Int_t i=0; i<1+fPairList->GetLast();i++){
    TObjArray* temparray = (TObjArray*)(fPairList->At(i));
    MSG("FarDetShieldPlank",Msg::kVerbose) << " " << i << ":" << endl;

    CandDigitHandle* dig1 = (CandDigitHandle*)(temparray->At(0));
    if( dig1 ){
      MSG("FarDetShieldPlank",Msg::kVerbose) << "  (1) " <<  dig1->GetPlexSEIdAltL().GetPlane() << " " << dig1->GetPlexSEIdAltL().GetEnd() << " " << dig1->GetCharge(CalDigitType::kNone) << " " << 1.0e9*dig1->GetSubtractedTime(CalTimeType::kNone) << " " << "(" << dig1->GetChannelId().GetCrate() << "-" << dig1->GetChannelId().GetVarcId() << "-" << dig1->GetChannelId().GetVmm() << "-" << dig1->GetChannelId().GetVaAdcSel() << "-" << dig1->GetChannelId().GetVaChip() << "-" << dig1->GetChannelId().GetVaChannel() << ")" << endl;
    }

    CandDigitHandle* dig2 = (CandDigitHandle*)(temparray->At(1));
    if( dig2 ){
      MSG("FarDetShieldPlank",Msg::kVerbose) << "  (2) " <<  dig2->GetPlexSEIdAltL().GetPlane() << " " << dig2->GetPlexSEIdAltL().GetEnd() << " " << dig2->GetCharge(CalDigitType::kNone) << " " << 1.0e9*dig2->GetSubtractedTime(CalTimeType::kNone) << " " << "(" << dig2->GetChannelId().GetCrate() << "-" << dig2->GetChannelId().GetVarcId() << "-" << dig2->GetChannelId().GetVmm() << "-" << dig2->GetChannelId().GetVaAdcSel() << "-" << dig2->GetChannelId().GetVaChip() << "-" << dig2->GetChannelId().GetVaChannel() << ")" << endl;
    }
  }

  // Check on number digits
  MSG("FarDetShieldPlank",Msg::kDebug) << "  DIGITS BEFORE = " << DigitsBefore << "  DIGITS AFTER = " << DigitsAfter << endl;

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

  for(i=0;i<1+fPairList->GetLast();i++){
    TObjArray* mypair = (TObjArray*)(fPairList->At(i));
    cxx.SetDataIn(mypair);
    FarDetShieldPlankHandle myplank = FarDetShieldPlank::MakeCandidate(ahh,cxx);
    myplank.SetName("FarDetShieldPlank");
    myplank.SetTitle(TString("Created by FarDetShieldPlankList"));
    myplanklist.AddDaughterLink(myplank);
  }

  for(i=0;i<500;i++){
    fDigiList[i].Clear();
  }

  fPairList->Delete();

  return;
}


void AlgFarDetShieldPlankList::Trace(const char* /*c*/) const
{

}

---