HarpoSimStrips.cxx

Go to the documentation of this file.

#include "HarpoSimStrips.h"

#include "TMath.h"


ClassImp(HarpoSimStrips)


HarpoSimStripHit::HarpoSimStripHit() :
TObject(), 
  fStripId(-1), 
  fTime(-1), 
  fQ(-1)
{
}

HarpoSimStripHit::HarpoSimStripHit(Int_t stripId, Double_t t, Double_t q) :
TObject(), 
  fStripId(stripId), 
  fTime(t), 
  fQ(q)
{
}
  
HarpoSimStripHit::HarpoSimStripHit(const HarpoSimStripHit &hit) :
  TObject(hit), 
  fStripId(hit.fStripId), 
  fTime(hit.fTime), 
  fQ(hit.fQ)
{
}
 
  
HarpoSimStrips::HarpoSimStrips()
{

  Info("HarpoSimStrip","Creating");

  fStripHitsX = new TClonesArray("HarpoSimStripHit",fkMaxNhits);
  fStripHitsY = new TClonesArray("HarpoSimStripHit",fkMaxNhits);

  fNhitsX = 0;
  fNhitsY = 0;
  fTrackId = -1;
  fStripWidth = 1.;
}

HarpoSimStrips::HarpoSimStrips(const HarpoSimStrips &strips):
  TObject(strips), 
  fStripHitsX(strips.fStripHitsX), 
  fStripHitsY(strips.fStripHitsY),
  fNhitsX(strips.fNhitsX),
  fNhitsY(strips.fNhitsY),
  fTrackId(strips.fTrackId),
  fStripWidth(strips.fStripWidth)
{

  // fStripHitsX = new TClonesArray("HarpoSimStripHit",fkMaxNhits);
  // fStripHitsY = new TClonesArray("HarpoSimStripHit",fkMaxNhits);

}

HarpoSimStrips::~HarpoSimStrips()
{

  fStripHitsX->Delete();
  fStripHitsY->Delete();

}

Int_t HarpoSimStrips::Reset()
{
  // fStripHitsX->Delete();
  // fStripHitsY->Delete();
  fStripHitsX = new TClonesArray("HarpoSimStripHit",fkMaxNhits);
  fStripHitsY = new TClonesArray("HarpoSimStripHit",fkMaxNhits);
  fNhitsX = 0;
  fNhitsY = 0;
  return 1;
}

Int_t HarpoSimStrips::AddHit(Int_t plane, Int_t ch, Double_t t, Double_t q)
{

  if(plane==0){ 
    if(fNhitsX>=fkMaxNhits) return -1;
    new((*fStripHitsX)[fNhitsX]) HarpoSimStripHit(ch,t,q);
    fNhitsX++;
  }
  if(plane==1){
    if(fNhitsY>=fkMaxNhits) return -1;
    new((*fStripHitsY)[fNhitsY]) HarpoSimStripHit(ch,t,q);
    fNhitsY++;
  }

  return 1;
}


// Int_t HarpoSimStrips::AddCharge(Double_t charge, Double_t sigma, Double_t x, Double_t y, Double_t t)
// {

  // if(fNhitsX>=fkMaxNhits) return -1;

  // //Info("AddCharge","%d %g %g %g %g %g",fNhitsX,charge,sigma,x,y,t);
  

  // // Find Strip number (X coordinate)
  // Int_t iX = Int_t((x + 150.)/fStripWidth);
  // Int_t iY = Int_t((y + 150.)/fStripWidth);
  
  // Int_t nbins = Int_t(3*sigma/fStripWidth);
  // for(Int_t i = iX - nbins; i<iX+nbins+1; i++){
  //   Double_t Xi1 = fStripWidth*i-150;
  //   Double_t Xi2 = Xi1 + 0.379;
  //   Double_t Xi3 = Xi2 + 0.1;
  //   Double_t Xi4 = Xi2 + 0.421;
  //   Double_t chX = TMath::Abs(TMath::Erf((x-Xi1)/sigma) - TMath::Erf((x-Xi2)/sigma));
  //   Double_t chY = TMath::Abs(TMath::Erf((x-Xi3)/sigma) - TMath::Erf((x-Xi4)/sigma));
  //   chX *= charge;
  //   chY *= charge;
  //   if(chX>0.1 && fNhitsX<fkMaxNhits && i>=0 && i<NSTRIPSX){
  //     new((*fStripHitsX)[fNhitsX]) HarpoSimStripHit(i,t,chX);
  //     fNhitsX++;
  //   }
  //   for(Int_t j = iY - nbins; j<iY+nbins+1; j++){
  //     if(j<0 || j>=NSTRIPSY) continue;
  //     Double_t Yj1 = fStripWidth*j-150;
  //     Double_t Yj2 = Yj1 + 0.9;
  //     Double_t ch = TMath::Abs(TMath::Erf((y-Yj1)/sigma) - TMath::Erf((y-Yj2)/sigma));
  //     ch *= chY;
  //     if(ch>0.1 && fNhitsY<fkMaxNhits){
  //    new((*fStripHitsY)[fNhitsY]) HarpoSimStripHit(j,t,ch);
  //    fNhitsY++;
  //     }
  //   }
      
  // }
  
  // return 1;

// }