AlgPassThru.cxx

Go to the documentation of this file.

#include <cassert>

#include "Calibrator/Calibrator.h"
#include "CandData/CandRecord.h"
#include "MuonRemoval/AlgPassThru.h"
#include "CandDigit/CandDigitHandle.h"
#include "Candidate/CandContext.h"
#include "Conventions/Munits.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawDigit.h"
#include "RawData/RawHeader.h"
#include "RawData/RawRecord.h"
#include "Validity/VldContext.h"

#include "Algorithm/AlgConfig.h"


ClassImp(AlgPassThru)

CVSID("$Id: AlgPassThru.cxx,v 1.2 2008/05/07 19:30:57 tjyang Exp $");


AlgPassThru::AlgPassThru()
{
}

AlgPassThru::~AlgPassThru()
{
}

void AlgPassThru::RunAlg(AlgConfig &ac , CandHandle &ch, CandContext &cx)
{
 
  assert(ch.InheritsFrom("CandDigitHandle"));

  //
  //Alg config
  //
  int cZeroWeights = 1;
  int cDoScale = 0;
  double cScaleFactor = 1.;
  double tmpf = 0;
  int tmpi = 0;
  if(ac.Get("ZeroWeights", tmpi)) cZeroWeights = tmpi;  
  if(ac.Get("doscale", tmpi)) cDoScale = tmpi;  
  if(ac.Get("scalefactor", tmpf)) cScaleFactor = tmpf;  

  //
  //Setup output digit
  //
  CandDigitHandle &digitout = (CandDigitHandle &) ch;
  
  //
  //Get the input data
  //
  const TObjArray* input = dynamic_cast<const TObjArray*>(cx.GetDataIn());
  if(input==NULL){
    MSG("RmMu", Msg::kError) << " Unable to get input data "<<endl;
    return;
  }
  
  const CandDigitHandle* digitin = dynamic_cast<const CandDigitHandle*> (input->At(0));
  assert(digitin);
  
  //copy over everything blindly
  RawChannelId channel(digitin->GetChannelId());
  digitout.SetChannelId(channel);
  digitout.SetTime(digitin->GetTime(CalTimeType::kNone));
  digitout.SetVaErrorBits(digitin->GetVaErrorBits());
  digitout.SetQieErrorBits(digitin->GetQieErrorBits());
  digitout.SetRawDigitIndex(digitin->GetRawDigitIndex());  

  if(cDoScale==0){
    digitout.SetCharge(digitin->GetCharge(CalDigitType::kNone));
  }else{    
    if(channel.GetElecType() == ElecType::kVA)
      digitout.SetVaErrorBits(digitin->GetVaErrorBits() | 0x80);
    else
      digitout.SetQieErrorBits(digitin->GetQieErrorBits() | 0x80);

    digitout.SetCharge(digitin->GetCharge(CalDigitType::kNone)*cScaleFactor);
  }


  if(cZeroWeights==0){
    //retain the demux solutions
    PlexSEIdAltL psalt(digitin->GetPlexSEIdAltL());
    digitout.SetPlexSEIdAltL(psalt);
  }else{
    //re-multiplex the solution
    const VldContext &vldc = *(cx.GetCandRecord()->GetVldContext());
    Calibrator& calibrator = Calibrator::Instance();
    calibrator.Reset(vldc);        
    PlexHandle ph(vldc);
    PlexSEIdAltL altlist = ph.GetSEIdAltL(channel, 
                                          &(calibrator), 
                                          (int)(digitin->GetCharge(CalDigitType::kNone)), 
                                          digitin->GetTime());
    digitout.SetPlexSEIdAltL(altlist);
  }
}

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

---