AlgFarDetDataQuality.cxx

Go to the documentation of this file.

#include "AlgFarDetDataQuality.h"

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandData/CandRecord.h"
#include "Candidate/CandContext.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h" 

#include "FarDetDataQualityHandle.h"
#include "FarDetDeadChipHandle.h"
#include "SnarlHeader.h"
#include "LiHeader.h"
#include "VarcErrHeader.h"
#include "SpillServHeader.h"
#include "VaChip.h"

#include "TObjArray.h"

//
// $Log $
//

ClassImp(AlgFarDetDataQuality)

CVSID("$Id: AlgFarDetDataQuality.cxx,v 1.2 2006/04/21 17:40:11 blake Exp $");

AlgFarDetDataQuality::AlgFarDetDataQuality()
{

}
  
AlgFarDetDataQuality::~AlgFarDetDataQuality()
{

}

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

  FarDetDataQualityHandle& fardet = dynamic_cast<FarDetDataQualityHandle&>(ch);
  const TObjArray* arr = dynamic_cast<const TObjArray*>(cx.GetDataIn()); 

  TObjArray* fChipList = new TObjArray(); 

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

  Int_t crate,varc,vmm;
  Int_t hotplnchips,coldplnchips,busyplnchips;
  Int_t hotshldchips,coldshldchips,busyshldchips;
  Int_t tpmtbox,tpmtnanosec;  
  Int_t id,time;
  Int_t tot,tot0,tot1,tot2;
  Int_t i,j;

  SnarlHeader* snarlheader = (SnarlHeader*)(arr->At(0)); 
  TObjArray* badlist = (TObjArray*)(arr->At(1));
  TObjArray* busylist = (TObjArray*)(arr->At(2));

  if( snarlheader->GetHitPlnChannels()>0 ){
    fardet.SetRun(snarlheader->GetRun());
    fardet.SetSubRun(snarlheader->GetSubRun());
    fardet.SetSnarl(snarlheader->GetSnarl());
    fardet.SetDate(snarlheader->GetDate());
    fardet.SetTime(snarlheader->GetTime());
    fardet.SetTimeFrame(snarlheader->GetTimeFrame());
    fardet.SetMicrosec(snarlheader->GetMicrosec());
    fardet.SetNanosec(snarlheader->GetNanosec());
    fardet.SetLiChannels(snarlheader->GetLiChannels());
    fardet.SetHitPlnChannels(snarlheader->GetHitPlnChannels());
    fardet.SetHitShldChannels(snarlheader->GetHitShldChannels());
    fardet.SetCrateMask(snarlheader->GetCrateMask());

    coldplnchips=0; hotplnchips=0; coldshldchips=0; hotshldchips=0;
    for(i=0;i<1+badlist->GetLast();i++){
      VaChip* mychip = (VaChip*)(badlist->At(i));
      if( mychip->GetVaChipStatus()==VaChip::kCold ){
        fChipList->Add(mychip);
        if(mychip->GetPlane()>=0) coldplnchips++; if(mychip->GetShld()>=0) coldshldchips++; 
      }
      if( mychip->GetVaChipStatus()==VaChip::kHot ){
        fChipList->Add(mychip);
        if(mychip->GetPlane()>=0) hotplnchips++; if(mychip->GetShld()>=0) hotshldchips++;
      }
    }

    busyplnchips=0; busyshldchips=0; 
    for(i=0;i<1+busylist->GetLast();i++){
      VaChip* mychip = (VaChip*)(busylist->At(i));
      if( mychip->GetVaChipStatus()==VaChip::kBusy ){ 
        tot=0; tot0=0; tot1=0; tot2=0;
        for(j=0;j<1+busylist->GetLast();j++){
          VaChip* mychip2 = (VaChip*)(busylist->At(j));
          crate=mychip2->GetCrate(); varc=mychip2->GetVarc(); vmm=mychip2->GetVmm();
          if(mychip->IsSameVmm(crate,varc,vmm)){
            if( mychip->GetPlane()>=0 ){
              if( mychip2->GetNanosec()-mychip->GetNanosec()>=-500
               && mychip2->GetNanosec()-mychip->GetNanosec()<500 
               && 3*mychip2->GetVaadc()+mychip2->GetVaChip()<=3*mychip->GetVaadc()+mychip->GetVaChip() ) tot0++;
              if( mychip2->GetNanosec()-mychip->GetNanosec()<-500
               && mychip2->GetNanosec()-mychip->GetNanosec()>=-5000 ) tot1++;
              if( mychip2->GetNanosec()-mychip->GetNanosec()<-5000
               && mychip2->GetNanosec()-mychip->GetNanosec()>=-10000 ) tot2++;
            }
            if( mychip->GetShld()>=0 ){
              if( mychip2->GetNanosec()-mychip->GetNanosec()>=-5
               && mychip2->GetNanosec()-mychip->GetNanosec()<5 
               && 3*mychip2->GetVaadc()+mychip2->GetVaChip()<=3*mychip->GetVaadc()+mychip->GetVaChip() ) tot0++;
              if( mychip2->GetNanosec()-mychip->GetNanosec()<-5
               && mychip2->GetNanosec()-mychip->GetNanosec()>=-5000 ) tot1++;
              if( mychip2->GetNanosec()-mychip->GetNanosec()<-5000
               && mychip2->GetNanosec()-mychip->GetNanosec()>=-10000 ) tot2++;
            }
          }
        }
        tot=tot0; if(tot1>=1) tot+=tot1; if(tot1+tot2>=2) tot+=tot2;
        if( mychip->GetNanosec()-snarlheader->GetNanosec()>-tot*5000 ){
          fChipList->Add(mychip);      
          if(mychip->GetPlane()>=0) busyplnchips++; if(mychip->GetShld()>=0) busyshldchips++; 
       }
      } 
    }

    if(1+fChipList->GetLast()>0){
      for(i=0;i<1+fChipList->GetLast();i++){
        VaChip* mychip = (VaChip*)(fChipList->At(i));
        cxx.SetDataIn(mychip);
        FarDetDeadChipHandle mydeadchip = FarDetDeadChip::MakeCandidate(ahh,cxx);
        mydeadchip.SetName(TString("FarDetDeadChipHandle"));
        mydeadchip.SetTitle(TString("Created by FarDetDataQualityModule"));
        fardet.AddDaughterLink(mydeadchip);
      }
    }

    fardet.SetDeadPlnChips(coldplnchips);
    fardet.SetHotPlnChips(hotplnchips);
    fardet.SetColdPlnChips(coldplnchips);
    fardet.SetBusyPlnChips(busyplnchips);

    fardet.SetDeadShldChips(coldshldchips);
    fardet.SetHotShldChips(hotshldchips);
    fardet.SetColdShldChips(coldshldchips);
    fardet.SetBusyShldChips(busyshldchips);
    
  }

  LiHeader* liheader = (LiHeader*)(arr->At(3));
  TObjArray* tpmtlist = (TObjArray*)(arr->At(4));
  if( liheader->GetTimeFrame()==snarlheader->GetTimeFrame() ){

    fardet.SetLiBox(liheader->GetBox());
    fardet.SetLiLed(liheader->GetLed());
    fardet.SetLiPulseHeight(liheader->GetPulseHeight());
    fardet.SetLiPulseWidth(liheader->GetPulseWidth());

    id=-1; time=-1000000000;
    tpmtbox=-1; tpmtnanosec=0;  
    for(i=0;i<1+tpmtlist->GetLast();i++){
      VaChip* mychip = (VaChip*)(tpmtlist->At(i));   
      if( mychip->GetNanosec()-snarlheader->GetNanosec()<1000
       && mychip->GetNanosec()-snarlheader->GetNanosec()>time ){
        id=i; time=mychip->GetNanosec()-snarlheader->GetNanosec();
      }
    }
    if(id>-1){
      VaChip* mychip = (VaChip*)(tpmtlist->At(id));   
      tpmtbox=mychip->GetVaChannel();
      tpmtnanosec=mychip->GetNanosec()-snarlheader->GetNanosec();
    }
  
    fardet.SetTpmtBox(tpmtbox);
    fardet.SetTpmtNanosec(tpmtnanosec);
  }

  VarcErrHeader* varcerrheader = (VarcErrHeader*)(arr->At(5));
  if( varcerrheader->GetTimeFrame()==snarlheader->GetTimeFrame() ){
    fardet.SetSparsErrors(varcerrheader->GetSparsErrors());
    fardet.SetEtcErrors(varcerrheader->GetEtcErrors());
  }

  SpillServHeader* spillservheader = (SpillServHeader*)(arr->At(6));
  if( spillservheader->GetTimeFrame()==snarlheader->GetTimeFrame() ){
    fardet.SetSpillStatus(spillservheader->GetSpillStatus());
    fardet.SetSpillType(spillservheader->GetSpillType());
    fardet.SetSpillGpsError(spillservheader->GetSpillGpsError());
  }

  delete fChipList;

  return;
}

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

}

---