AlgFarDetStripList.cxx

Go to the documentation of this file.

#include "AlgFarDetStripList.h"
#include "FarDetStripHandle.h"
#include "FarDetStripListHandle.h"

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

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

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

#include "VaDigit.h"

ClassImp(AlgFarDetStripList)    

CVSID("$Id: AlgFarDetStripList.cxx,v 1.2 2006/07/13 12:10:34 blake Exp $");

AlgFarDetStripList::AlgFarDetStripList() :
  fPairList(0),
  fVaList(0)
{
  fPairList = new TObjArray();
  fVaList = new TObjArray();
}   

AlgFarDetStripList::~AlgFarDetStripList()
{
  if(fPairList) delete fPairList;
  if(fVaList) delete fVaList;
}
   
void AlgFarDetStripList::RunAlg(AlgConfig &ac, CandHandle &ch, CandContext &cx)
{
  MSG("FarDetStrip",Msg::kDebug) << " AlgFarDetStripList::RunAlg(...) " << endl;

  FarDetStripListHandle& mystriplist = dynamic_cast<FarDetStripListHandle&>(ch);

  Int_t fRemoveCrossTalk;
  fRemoveCrossTalk = ac.GetInt("RemoveCrossTalk");
  MSG("FarDetStrip",Msg::kVerbose) << " RemoveCrossTalk=" << fRemoveCrossTalk << endl;

  Int_t DigitsBefore=0,DigitsAfter=0;

  Double_t overlap;
  Int_t i,j,ie,iw,je,jw,ke,kw;
  Int_t pln,bin,match,xtalk;
  Int_t crate,varc,vmm,vaadc,vachip,vachannel;
  TObjArray tmpe,tmpw,tmps;
  TObjArray eastlist,westlist,tmplist;  
  TObjArray eastxtalklist,westxtalklist;

  // Unpack CandContext
  const CandDigitListHandle* cdlh = dynamic_cast<const CandDigitListHandle*>(cx.GetDataIn());
  TIter digitr(cdlh->GetDaughterIterator());
  while(CandDigitHandle* cdh = dynamic_cast<CandDigitHandle*>(digitr())){
    if(cdh){
      pln = cdh->GetPlexSEIdAltL().GetPlane();
      if( pln>0 && pln<500 ){
        RawChannelId rawch = cdh->GetChannelId();    
        crate = rawch.GetCrate();
        varc = rawch.GetVarcId();
        vmm = rawch.GetVmm();
        vaadc = rawch.GetVaAdcSel();
        vachip = rawch.GetVaChip();
        vachannel = rawch.GetVaChannel();
        if( crate<16 && varc<3 && vmm<6 && vaadc<2 && vachip<3 
         && vachannel>=2 && vachannel<=17 ){

          bin=-1;
          if(cdh->GetPlexSEIdAltL().GetEnd()==StripEnd::kPositive){ // west (odd)
            bin = 36*(crate/2)+12*(1+(1-varc))+2*(vmm)+(vaadc);
          }
          if(cdh->GetPlexSEIdAltL().GetEnd()==StripEnd::kNegative){ // east (even)
            bin = 36*(crate/2)+12*(1-(1-varc))+2*(vmm)+(vaadc);
          }
          if( bin>=0 && bin<300 ){
            fVfbList[bin].Add(cdh);

            CandDeMuxDigitHandle* cddh = dynamic_cast<CandDeMuxDigitHandle*>(cdh);
            xtalk=0;
            if( ( cddh->GetDeMuxDigitFlagWord()<8 )
             && ( (cddh->GetDeMuxDigitFlagWord() & CandDeMuxDigit::kXTalk)==(CandDeMuxDigit::kXTalk) ) ){
              xtalk=1; 
            }

            MSG("FarDetStrip",Msg::kVerbose) << " digit : "
                                             << " " << cdh->GetPlexSEIdAltL().GetPlane() 
                                             << " " << cdh->GetPlexSEIdAltL().GetEnd()
                                             << " " << cdh->GetCharge(CalDigitType::kNone) 
                                             << " " << 1.0e9*cdh->GetSubtractedTime(CalTimeType::kNone) 
                                             << " " << xtalk << endl;

            DigitsBefore++; 
          }
        }
      }
    } 
  }

  // Loop over Vfbs
  for(i=0;i<300;i++){
    if(1+fVfbList[i].GetLast()>0){

      // Deal with Crosstalk
      for(j=0;j<1+fVfbList[i].GetLast();j++){
        CandDeMuxDigitHandle* digit = (CandDeMuxDigitHandle*)(fVfbList[i].At(j));
        xtalk=0;
        if( fRemoveCrossTalk
         && ( digit->GetDeMuxDigitFlagWord()<8 )
         && ( (digit->GetDeMuxDigitFlagWord() & CandDeMuxDigit::kXTalk)==(CandDeMuxDigit::kXTalk) ) ){
          xtalk=1; 
        }
        if(digit->GetPlexSEIdAltL().GetEnd()==StripEnd::kPositive){ 
          if(!xtalk){
            VaDigit* vadigit = new VaDigit(digit);
            westlist.Add(vadigit); fVaList->Add(vadigit);
          } 
          else westxtalklist.Add(digit); 
        }
        if(digit->GetPlexSEIdAltL().GetEnd()==StripEnd::kNegative){ 
          if(!xtalk){
            VaDigit* vadigit = new VaDigit(digit);
            eastlist.Add(vadigit); fVaList->Add(vadigit);
          }
          else eastxtalklist.Add(digit); 
        }
      }

      for(iw=0;iw<1+westxtalklist.GetLast();iw++){
        CandDigitHandle* xtalkdigit = (CandDigitHandle*)(westxtalklist.At(iw));
        tmplist.Clear();
        for(jw=0;jw<1+westlist.GetLast();jw++){
          VaDigit* vadigit = (VaDigit*)(westlist.At(jw)); 
          if(vadigit->IsSamePmt(xtalkdigit)) tmplist.Add(vadigit);
        }
        if( 1+tmplist.GetLast()>0 ){
          overlap=0.0;
          for(jw=0;jw<1+tmplist.GetLast();jw++){
            VaDigit* vadigit = (VaDigit*)(tmplist.At(jw));
            overlap += vadigit->GetCrossTalk(xtalkdigit);
          }
          if( overlap>0.0){
            for(jw=0;jw<1+tmplist.GetLast();jw++){
              VaDigit* vadigit = (VaDigit*)(tmplist.At(jw));
              vadigit->AddCrossTalk(xtalkdigit,overlap);
            }
          }
        }
        DigitsAfter++; 
      }

      for(ie=0;ie<1+eastxtalklist.GetLast();ie++){
        CandDigitHandle* xtalkdigit = (CandDigitHandle*)(eastxtalklist.At(ie));
        tmplist.Clear();
        for(je=0;je<1+eastlist.GetLast();je++){
          VaDigit* vadigit = (VaDigit*)(eastlist.At(je));  
          if(vadigit->IsSamePmt(xtalkdigit)) tmplist.Add(vadigit);
        }
        if( 1+tmplist.GetLast()>0 ){
          overlap = 0.0;
          for(je=0;je<1+tmplist.GetLast();je++){
            VaDigit* vadigit = (VaDigit*)(tmplist.At(je));      
            overlap += vadigit->GetCrossTalk(xtalkdigit);
          }
          if( overlap>0.0 ){
            for(je=0;je<1+tmplist.GetLast();je++){
              VaDigit* vadigit = (VaDigit*)(tmplist.At(je));
              vadigit->AddCrossTalk(xtalkdigit,overlap);
            }
          }
        }
        DigitsAfter++;
      }

      // Match up Digits
      for(kw=0;kw<1+westlist.GetLast();kw++){
        VaDigit* digw = (VaDigit*)(westlist.At(kw));
        match=0;
        for(ke=0;ke<1+eastlist.GetLast();ke++){
          VaDigit* dige = (VaDigit*)(eastlist.At(ke));
          if( digw->IsSameStrip(dige) 
           && digw->IsSameTime(dige) ) match=1;
        }
        if(match) tmpw.Add(digw); else tmps.Add(digw);
      }

      for(ke=0;ke<1+eastlist.GetLast();ke++){
        VaDigit* dige = (VaDigit*)(eastlist.At(ke));
        match=0;
        for(kw=0;kw<1+westlist.GetLast();kw++){
          VaDigit* digw = (VaDigit*)(westlist.At(kw));
          if( dige->IsSameStrip(digw) 
           && dige->IsSameTime(digw) ) match=1;
        }
        if(match) tmpe.Add(dige); else tmps.Add(dige);
      }

      if( 1+tmpe.GetLast()>0 && 1+tmpw.GetLast()>0 ){
        for(ke=0;ke<1+tmpe.GetLast();ke++){
          VaDigit* dige = (VaDigit*)(tmpe.At(ke));
          for(kw=0;kw<1+tmpw.GetLast();kw++){
            VaDigit* digw = (VaDigit*)(tmpw.At(kw));
            if( dige->IsSameStrip(digw) 
             && dige->IsSameTime(digw) ){
              TObjArray* myarray = new TObjArray();
              myarray->Add(dige); DigitsAfter++;
              myarray->Add(digw); DigitsAfter++;

              fPairList->Add(myarray);
            }
          }
        }
      }

      for(j=0;j<1+tmps.GetLast();j++){
        VaDigit* dig = (VaDigit*)(tmps.At(j));
        TObjArray* myarray = new TObjArray();
        myarray->Add(dig); DigitsAfter++; 
        fPairList->Add(myarray);
      }
      
      tmpe.Clear(); tmpw.Clear(); tmps.Clear();
      eastxtalklist.Clear(); westxtalklist.Clear(); 
      eastlist.Clear(); westlist.Clear();
    }
  }

  // Digit Check
  MSG("FarDetStrip",Msg::kDebug) << "  DIGITS BEFORE = " << DigitsBefore << "  DIGITS AFTER = " << DigitsAfter << endl;

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

  for(i=0;i<1+fPairList->GetLast();i++){
    TObjArray* myarray = (TObjArray*)(fPairList->At(i));

    cxx.SetDataIn(myarray);
    FarDetStripHandle mystrip = FarDetStrip::MakeCandidate(ahh,cxx);
    mystrip.SetName("FarDetStrip");
    mystrip.SetTitle(TString("Created by AlgFarDetStripList"));

    mystriplist.AddDaughterLink(mystrip);
  }

  for(i=0;i<300;i++){
    fVfbList[i].Clear();
  }

  fPairList->Delete();
  fVaList->Delete();

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

}

---