Helpers.cxx

Go to the documentation of this file.

#include <iostream>
using std::cout;
using std::endl;

#include "TClonesArray.h"
#include "CalDetSI/Helpers.h"
#include "Navigation/NavKey.h"
#include "Conventions/Detector.h"
#include "Conventions/ElecType.h"
#include "Conventions/StripEnd.h"
#include "RerootExodus/RerootExodus.h"

NavKey CalHelpers::KeyFromPlane(const CandDigitHandle* cdh)
{
     return GetPlane(*cdh);
}

NavKey CalHelpers::KeyFromStrip(const CandDigitHandle* cdh)
{
     return GetBestStrip(*cdh);
}

NavKey CalHelpers::KeyFromEnd(const CandDigitHandle* cdh)
{
     return GetBestEnd(*cdh);
}

NavKey CalHelpers::KeyFromTime(const CandDigitHandle* cdh)
{
     return cdh->GetPlexSEIdAltL().GetBestItem().GetTime();
}

NavKey CalHelpers::KeyFromDigitTime(const CandDigitHandle* cdh)
{
     return cdh->GetTime();
}


NavKey CalHelpers::KeyFromCrate(const CandDigitHandle* cdh)
{
     return cdh->GetChannelId().GetCrate();
}

NavKey CalHelpers::KeyFromVarc(const CandDigitHandle* cdh)
{
     return cdh->GetChannelId().GetVarcId();
}


CalHelpers::MCDigit::MCDigit():
  
  fMCPlane(0),
  fMCExtr(0),
  fMCCell(0),
  fMCTPos(0.),
  fMCRAWPE(0.),
  fMCCorrPE(0.),
  fMCTDC(0.),
  fMCTubePixel(0),
  fMCSignalPE(0.),
  fMCInitialTDC(0.),
  fMCAveDistTrue(0.),
  fMCHitBits(0),
  fMCCellPos(0),
  fMCPlanePos(0),
  fMCSumETrue(0)
{}

CalHelpers::MCDigit::MCDigit(const PlexStripEndId pse, REROOT_FLSDigit flsdig)
{  
  fMCPlane = flsdig.IPln();
  fMCExtr = flsdig.IExtr();
  fMCCell = flsdig.ICell(); 
  fMCTPos = flsdig.TPos();
  fMCHitBits = flsdig.HitBits(); 
  fMCCellPos = flsdig.CellPos();
  fMCPlanePos = flsdig.PlanePos(); 

  if(pse.GetEnd()==StripEnd::kNegative){
    fMCRAWPE = flsdig.RawA();
    fMCCorrPE = flsdig.CorrA();
    fMCTDC = flsdig.TDCA();
    fMCTubePixel = flsdig.TubePixelA(); 
    fMCSignalPE = flsdig.SignalPEA();
    fMCInitialTDC = flsdig.InitialTDCA();
    fMCAveDistTrue = flsdig.AveDistTrueA();
  }
  else{
    fMCRAWPE = flsdig.RawB();
    fMCCorrPE = flsdig.CorrB();
    fMCTDC = flsdig.TDCB();
    fMCTubePixel = flsdig.TubePixelB(); 
    fMCSignalPE = flsdig.SignalPEB();
    fMCInitialTDC = flsdig.InitialTDCB();
    fMCAveDistTrue = flsdig.AveDistTrueB();
  }
  // added by MAK
  fMCSumETrue=flsdig.SumETrue();
  if(fMCRAWPE<=0){
      cout<<"MCDigit Constructor: fMCRAWPE <= 0"<<endl
          <<"flsdig.RawA(): "<<flsdig.RawA()<<endl
          <<"flsdig.RawB(): "<<flsdig.RawB()<<endl
          <<"Plane: "<<pse.GetPlane()<<endl
          <<"Strip: "<<pse.GetStrip()<<endl
          <<"Strip End: "<<pse.GetEnd()<<endl
          <<endl;
  }
}

CalHelpers::MCDigit::MCDigit(const MCDigit &mcd):
  fMCPlane(mcd.fMCPlane),
  fMCExtr(mcd.fMCExtr),
  fMCCell(mcd.fMCCell),
  fMCTPos(mcd.fMCTPos),
  fMCRAWPE(mcd.fMCRAWPE),
  fMCCorrPE(mcd.fMCCorrPE),
  fMCTDC(mcd.fMCTDC),
  fMCTubePixel(mcd.fMCTubePixel),
  fMCSignalPE(mcd.fMCSignalPE),
  fMCInitialTDC(mcd.fMCInitialTDC),
  fMCAveDistTrue(mcd.fMCAveDistTrue),
  fMCHitBits(mcd.fMCHitBits),
  fMCCellPos(mcd.fMCCellPos),
  fMCPlanePos(mcd.fMCPlanePos),
  fMCSumETrue(mcd.fMCSumETrue)
{}

  

CalHelpers::MCDigit::~MCDigit() {}

CalHelpers::MCDigitMap::MCDigitMap():
   fMCDigitMap()
{}

CalHelpers::MCDigitMap::~MCDigitMap() {}

void CalHelpers::MCDigitMap::MakeMCDigitMap()
{

  static const TClonesArray *mdcl = RerootExodus::GetFLSDigitList();

  for(int i=0;i<mdcl->GetEntries();i++){
    const REROOT_FLSDigit *dig = dynamic_cast<const REROOT_FLSDigit *>(mdcl->At(i));
    if(dig == 0){
      continue;
    }
    PlexStripEndId psea = RerootExodus::PECAB2SEId(dig->IPln(), 
                                                  dig->IExtr(), dig->ICell(), 0);
    // for debugging
    cout<<"Gave iab = 0"<<endl;
    cout<<"Got back stripend: "<<psea.GetEnd()<<endl;

    MCDigit mcdig1(psea, *dig);
    fMCDigitMap[psea]=mcdig1;
    PlexStripEndId pseb = RerootExodus::PECAB2SEId(dig->IPln(), 
                                                  dig->IExtr(), dig->ICell(), 1);
    // for debugging
    cout<<"Gave iab = 1"<<endl;
    cout<<"Got back stripend: "<<pseb.GetEnd()<<endl;

    MCDigit mcdig2(pseb, *dig);
    fMCDigitMap[pseb]=mcdig2;

  }
}

---