StripKG.cxx

Go to the documentation of this file.

// Cluster3D package
//
// StripKG.cxx
//
// StripKG - a strip
//
// Author:  Katarzyna Grzelak K.Grzelak1@physics.ox.ac.uk
//

#include "Cluster3D/StripKG.h"
#include "Algorithm/AlgConfig.h"
#include "MessageService/MsgService.h" // needed for CVS
#include "RecoBase/CandStripHandle.h"     
#include "UgliGeometry/UgliGeomHandle.h"   // UgliGeom handling

ClassImp(StripKG)

CVSID("$Id: StripKG.cxx,v 1.2 2007/03/01 17:44:17 rhatcher Exp $");

StripKG::StripKG():handle(0),index(0),strip(0),plane(0),planeview(0),
digit(0),tpos(0.),zpos(0.),rawE(0.),rawW(0.),
siglinE(0.),siglinW(0.),sigcorE(0.),sigcorW(0.),
peE(0.),peW(0.),timerawbE(0.),timerawbW(0.),timeraweE(0.),
timeraweW(0.),timeE(0.),timeW(0.),timephE(0.),
timephW(0.),hwidth(0.),hlength(0.),hthick(0.),clusternb(999),xtalk(false),
ambig(0){
//
//  Purpose:  Default constructor
//
  MSG("Cluster3D", Msg::kVerbose)
    << "Default constructor " << endl;
}
//............................................................................

void StripKG::Print(){
//
//  Purpose:  Prints contents of the strip object.
//
//  Arguments: n/a
//
//  Return:  n/a
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//

  cout << "Content of the Strip object" << endl;
  cout << "index= " << index 
       << "strip= " << strip
       << "plane= " << plane
       << "planeview= " << planeview
       << "digit= " << digit  
       << endl;
}

//............................................................................

void StripKG::Set(CandStripHandle* stp, UgliGeomHandle *ugh, Int_t lp){
//
//  Purpose: Set StripKG data members 
//
//  Arguments: 
//           stp   in  - strip handle   
//           ugh   in  - UgliGeom handle
//            lp   in  - serial strip number
//              
//  Return:  n/a
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============

  // Light propagation velocity in a clear fiber in m/ns
  Double_t cClear = 0.56*Munits::c_light*Munits::ns; 
  // Light propagation velocity in a Wls  fiber in m/ns
  Double_t cWls = 0.73*Munits::c_light*Munits::ns; 
  // Light propagation velocity in a scintillator in m/ns
  //Double_t cScint = 0.56*Munits::c_light*Munits::ns; 

  // Two ends of strip
  StripEnd::StripEnd_t stpend[2] = {StripEnd::kNegative,StripEnd::kPositive};

  handle=stp; 
  index=lp;
  strip=stp->GetStrip();  
  plane=stp->GetPlane();
  planeview=(Int_t)(stp->GetPlaneView());
  digit=stp->GetNDigit();
  tpos=stp->GetTPos();
  zpos=stp->GetZPos();
  rawE=stp->GetCharge(CalDigitType::kNone,stpend[0]);
  rawW=stp->GetCharge(CalDigitType::kNone,stpend[1]);
  siglinE=stp->GetCharge(CalDigitType::kSigLin,stpend[0]);
  siglinW=stp->GetCharge(CalDigitType::kSigLin,stpend[1]);
  sigcorE=stp->GetCharge(CalDigitType::kSigCorr,stpend[0]);
  sigcorW=stp->GetCharge(CalDigitType::kSigCorr,stpend[1]);
  peE=stp->GetCharge(CalDigitType::kPE,stpend[0]);
  peW=stp->GetCharge(CalDigitType::kPE,stpend[1]);
  timerawbE=stp->GetBegTime(stpend[0])/Munits::ns;
  timerawbW=stp->GetBegTime(stpend[1])/Munits::ns; 
  timeraweE=stp->GetEndTime(stpend[0])/Munits::ns;
  timeraweW=stp->GetEndTime(stpend[1])/Munits::ns;

  PlexStripEndId stripid = stp->GetStripEndId();
  UgliStripHandle stph = ugh->GetStripHandle(stripid);

  Double_t tClearE = stph.ClearFiber(stpend[0])/cClear;
  Double_t tWlsE = stph.WlsPigtail(stpend[0])/cWls;
  timeE=timerawbE-tWlsE-tClearE;

  Double_t tClearW = stph.ClearFiber(stpend[1])/cClear;
  Double_t tWlsW = stph.WlsPigtail(stpend[1])/cWls;
  timeW=timerawbW-tWlsW-tClearW;

  timephE=stp->GetTime(stpend[0])/Munits::ns;
  timephW=stp->GetTime(stpend[1])/Munits::ns;

  hwidth=stph.GetHalfWidth();
  hlength=stph.GetHalfLength();
  hthick=stph.GetHalfThickness();
}

//............................................................................

void StripKG::TimePos(AlgConfig &ac, Bool_t oneStp, Float_t &dt, Float_t &timpos, Float_t &halflength) const{
//
//  Purpose:  Return hit position from timing information.
//
//  Arguments: 
//           ac           in  -  AlgConfig 
//           oneStp       in  -  true = one hit per plane => use timing info
//                               false = more than one hit => do not use timing
//           dt          out  -  time difference [ns] 
//           timpos      out  -  position from time [strips]
//           halflength  out  -  half of the length of the cell [strips]
//              
//  Return:  n/a
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============
//
// Fit results: mean strip & velocity in the scintillator in strips/ns
// avStp = mean strip/length of half of the strip in strips
// vScint/2 = velocity in the scintillator [strip/ns]
// dtRes = time resolution [ns] 
// widthStp = strip width [m]

  const Double_t vScint = ac.GetDouble("vLightInScint");
  const Double_t t1t2Res = ac.GetDouble("t1-t2ResolutionInNs");

  Double_t dtRes;

  if(oneStp==true){
    dtRes = 3.*t1t2Res;
  }
  else {
    dtRes = t1t2Res;   
  } 

  Float_t avStp = 96.0;

  dt=timeE-timeW; 

  if(Timing()==2 && planeview==2 && oneStp==true){
    timpos=avStp-vScint*dt;
    halflength=dtRes*vScint;
  }
  else if(Timing()==2 && planeview==3 && oneStp==true){
    timpos=avStp+vScint*dt;
    halflength=dtRes*vScint;
  }
  else{
    timpos=avStp;
    halflength=avStp;
  }
}

//............................................................................

Int_t StripKG::Timing() const{   
//
//  Purpose:  Return 
//
//  Arguments: 
//              
//  Return:  n/a
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============

  Int_t nbEnd;
  if(timeE>-1000. && timeW>-1000.){
    nbEnd=2; 
  }
  else if(timeE>-1000. || timeW>-1000.){
    nbEnd=1;
  }
  else{
    nbEnd=0;
  }
  return nbEnd;
}

//............................................................................

Bool_t StripKG::XTalk(AlgConfig &ac){
//
//  Purpose:  Is it a crosstalk hit ?
//
//  Arguments: 
//           ac         in  -  AlgConfig
//              
//  Return:  Bool_t crosstalk hit or not ?
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============

  const Double_t peCut = ac.GetDouble("xTalkPeCut");
  MSG("Cluster3D", Msg::kVerbose)
   << "pe cut = " << peCut << endl;

  if(digit==1 && peE<peCut && peW<peCut){
    xtalk=true;
    return true;    
  }
  else{
    return false;
    xtalk=false;
  }
}

//............................................................................

CandStripHandle* StripKG::GetHandle() const{
//
//  Purpose:  return handle to a strip
//
//  Arguments: n/a
//              
//  Return:  CandStripHandle* handle 
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============
  
  return handle;
}

//............................................................................

---