AlgFarDetStrip.cxx

Go to the documentation of this file.

#include "AlgFarDetStrip.h"
#include "FarDetStripHandle.h"

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

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

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

#include "Plex/PlexHandle.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "UgliGeometry/UgliStripHandle.h"
#include "Validity/VldContext.h"
#include "Conventions/SimFlag.h"

#include "TObjArray.h"
#include "VaDigit.h"

ClassImp(AlgFarDetStrip)    

CVSID("$Id: AlgFarDetStrip.cxx,v 1.3 2006/04/21 16:57:41 chapman Exp $");

AlgFarDetStrip::AlgFarDetStrip()
{
  MSG("FarDetStrip",Msg::kVerbose) << " AlgFarDetStrip::AlgFarDetStrip() " << endl;
  
}   

AlgFarDetStrip::~AlgFarDetStrip()
{
  MSG("FarDetStrip",Msg::kVerbose) << " AlgFarDetStrip::~AlgFarDetStrip() " << endl;
  
}
   
void AlgFarDetStrip::RunAlg(AlgConfig& /* ac */, 
                            CandHandle& ch, CandContext& cx)
{
  MSG("FarDetStrip",Msg::kVerbose) << " AlgFarDetStrip::RunAlg(...) " << endl;

  FarDetStripHandle& mystrip = dynamic_cast<FarDetStripHandle&>(ch);

  Double_t corr=1.0;
  Int_t i;

  const TObjArray* myarray = dynamic_cast<const TObjArray*>(cx.GetDataIn());
  CandRecord* candrec = (CandRecord*)(cx.GetCandRecord());
  VldContext* vldc = (VldContext*)(candrec->GetVldContext());
  UgliGeomHandle ugh(*vldc);
  PlexHandle ph(*vldc);

  Int_t myXtalk[3];
  Double_t myRawTime[3];
  Double_t myCalTime[3];
  Double_t myChargeNone[2];
  Double_t myChargeAdc[2];
  Double_t myChargePE[2];
  Double_t myChargeSigLin[2];
  Double_t myChargeSigCorr[2];
  Double_t myChargePEcorr[2];
  Double_t myGreenFibre[2];
  Double_t myWlsPigtail[2];
  Double_t myClearFibre[2];
  Double_t myHalfLength;

  for(i=0;i<3;i++){
    myRawTime[i] = 0.0;
    myCalTime[i] = 0.0;
    myXtalk[i] = 0;
  }

  for(i=0;i<2;i++){
    myChargeNone[i] = 0.0;
    myChargeAdc[i] = 0.0;
    myChargePE[i] = 0.0;
    myChargeSigLin[i] = 0.0;
    myChargeSigCorr[i] = 0.0;
    myChargePEcorr[i] = 0.0;
    myGreenFibre[i] = 0.0;
    myWlsPigtail[i] = 0.0;
    myClearFibre[i] = 0.0;
  }

  myHalfLength = 0.0;

 
  for(i=0;i<1+myarray->GetLast();i++){
    VaDigit* vadigit = (VaDigit*)(myarray->At(i));
    corr = (vadigit->GetChargeAdc()+vadigit->GetCrossTalkChargeAdc())/(vadigit->GetChargeAdc());

    CandDigitHandle* digit = (CandDigitHandle*)(vadigit->GetCandDigitHandle());
    PlexStripEndId pSEId = digit->GetPlexSEIdAltL().GetBestSEId();

    if(!i){
      mystrip.SetBestStripEndId(pSEId);  
      UgliStripHandle striphandle = ugh.GetStripHandle(pSEId);
      UgliScintPlnHandle scintplnhandle = ugh.GetScintPlnHandle(pSEId);

      myGreenFibre[0] = striphandle.GetHalfLength();
      myWlsPigtail[0] = striphandle.WlsPigtail(StripEnd::kNegative);
      myClearFibre[0] = striphandle.ClearFiber(StripEnd::kNegative);

      myGreenFibre[1] = striphandle.GetHalfLength();
      myWlsPigtail[1] = striphandle.WlsPigtail(StripEnd::kPositive);
      myClearFibre[1] = striphandle.ClearFiber(StripEnd::kPositive);

      myHalfLength = striphandle.GetHalfLength();

      mystrip.SetPlane(pSEId.GetPlane());
      mystrip.SetStrip(pSEId.GetStrip());
      mystrip.SetTPos(striphandle.GetTPos());
      mystrip.SetZPos(scintplnhandle.GetZ0());
    }

    if(pSEId.GetEnd()==StripEnd::kNegative){
      CandDeMuxDigitHandle* demuxdigit = (CandDeMuxDigitHandle*)(digit);
      if( (demuxdigit->GetDeMuxDigitFlagWord()<8)
       && ( (demuxdigit->GetDeMuxDigitFlagWord() & CandDeMuxDigit::kXTalk)==(CandDeMuxDigit::kXTalk) ) ){
        myXtalk[0]=1; myXtalk[2]=1;
      }
      myRawTime[0] = digit->GetSubtractedTime(CalTimeType::kNone);
      myCalTime[0] = digit->GetSubtractedTime(CalTimeType::kT0);
      myChargeNone[0] = digit->GetCharge(CalDigitType::kNone);
      myChargeAdc[0] = corr*(digit->GetCharge(CalDigitType::kNone));
      myChargePE[0] = corr*(digit->GetCharge(CalDigitType::kPE));
      myChargeSigLin[0] = corr*(digit->GetCharge(CalDigitType::kSigLin));
      myChargeSigCorr[0] = corr*(digit->GetCharge(CalDigitType::kSigCorr));
      myChargePEcorr[0] = myChargeSigCorr[0]/67.5;
      
      MSG("FarDetStrip",Msg::kVerbose) << " -ve : " 
                                       << " Qraw=" << myChargeNone[0] 
                                       << " Qadc=" << myChargeAdc[0]
                                       << " Qpe=" << myChargePE[0]
                                       << " QpeCorr=" << myChargePEcorr[0] << endl;
    }

    if(pSEId.GetEnd()==StripEnd::kPositive){
      CandDeMuxDigitHandle* demuxdigit = (CandDeMuxDigitHandle*)(digit);
      if( (demuxdigit->GetDeMuxDigitFlagWord()<8)
       && ( (demuxdigit->GetDeMuxDigitFlagWord() & CandDeMuxDigit::kXTalk)==(CandDeMuxDigit::kXTalk) ) ){
        myXtalk[1]=1; myXtalk[2]=1;
      }
      myRawTime[1] = digit->GetSubtractedTime(CalTimeType::kNone);
      myCalTime[1] = digit->GetSubtractedTime(CalTimeType::kT0);
      myChargeNone[1] = digit->GetCharge(CalDigitType::kNone);
      myChargeAdc[1] = corr*(digit->GetCharge(CalDigitType::kNone));
      myChargePE[1] = corr*(digit->GetCharge(CalDigitType::kPE));
      myChargeSigLin[1] = corr*(digit->GetCharge(CalDigitType::kSigLin));
      myChargeSigCorr[1] = corr*(digit->GetCharge(CalDigitType::kSigCorr));
      myChargePEcorr[1] = myChargeSigCorr[1]/67.5;
      
      MSG("FarDetStrip",Msg::kVerbose) << " +ve : " 
                                       << " Qraw=" << myChargeNone[1] 
                                       << " Qadc=" << myChargeAdc[1]
                                       << " Qpe=" << myChargePE[1]
                                       << " QpeCorr=" << myChargePEcorr[1] << endl;
                                          
    }
    
    mystrip.AddDaughterLink(*digit);
  }

  Double_t Sqrt=0.0,Sqct=0.0,Sq=0.0;
  for(i=0;i<2;i++){
    Sqrt+=myRawTime[i]*myChargePE[i];
    Sqct+=myCalTime[i]*myChargePE[i]; 
    Sq+=myChargePE[i];
  }
  if(Sq>0.0){
    myRawTime[2] = Sqrt/Sq;
    myCalTime[2] = Sqct/Sq;
  }

  mystrip.SetTime(CalTimeType::kNone,myRawTime);
  mystrip.SetTime(CalTimeType::kT0,myCalTime);
  mystrip.SetCalTimeType(CalTimeType::kT0);

  mystrip.SetCharge(CalDigitType::kNone,myChargeAdc);
  mystrip.SetCharge(CalDigitType::kPE,myChargePE);
  mystrip.SetCharge(CalDigitType::kSigLin,myChargeSigLin);
  mystrip.SetCharge(CalDigitType::kSigCorr,myChargeSigCorr);

  mystrip.SetFibreLengths(myGreenFibre,myWlsPigtail,myClearFibre,myHalfLength);
  mystrip.SetXtalk(myXtalk);

  MSG("FarDetStrip",Msg::kDebug) << "  CANDSTRIP : "
                                 << " (" << mystrip.GetNDaughters() << ")"
                                 << " " << mystrip.GetPlane()
                                 << " " << mystrip.GetStrip() 
                                 << " " << mystrip.GetCharge()
                                 << " " << 1.0e9*mystrip.GetTime() << endl;
}
   
void AlgFarDetStrip::Trace(const char* /*c*/ ) const
{

}

        

---