HcalShape Class Reference

shaper for Hcal (not for HF) More...

#include <HcalShape.h>

Inheritance diagram for HcalShape:

List of all members.

Public Member Functions
void	computeShape ()
void	display () const
	HcalShape (const HcalShape &d)
	HcalShape ()
virtual double	operator() (double time) const
virtual double	timeToRise () const
virtual	~HcalShape ()
Private Attributes
int	nbin_
std::vector< float >	nt_

---

Detailed Description

shaper for Hcal (not for HF)

Definition at line 15 of file HcalShape.h.

---

Constructor & Destructor Documentation

HcalShape::HcalShape

(

)

Definition at line 4 of file HcalShape.cc.

References computeShape().

: nbin_(256),
  nt_(nbin_, 0.)
{
   computeShape();
}

HcalShape::HcalShape

(

const HcalShape &

d

)

Definition at line 11 of file HcalShape.cc.

                                     :
  CaloVShape(d),
  nbin_(d.nbin_),
  nt_(d.nt_)
{
}

virtual HcalShape::~HcalShape

(

)

[inline, virtual]

Definition at line 23 of file HcalShape.h.

{}

---

Member Function Documentation

void HcalShape::computeShape

(

)

Definition at line 24 of file HcalShape.cc.

References funct::exp(), i, j, gen::k, nbin_, lumiNorm::norm, nt_, and tmax.

Referenced by HcalShape().

{

  // pulse shape time constant_s in ns
  const float ts1  = 8.;          // scintillation time constant_s : 1,2,3
  const float ts2  = 10.;           
  const float ts3  = 29.3;         
  const float thpd = 4.;          // HPD current_ collection drift time
  const float tpre = 9.;          // preamp time constant_ (refit on TB04 data)
  
  const float wd1 = 2.;           // relative weights of decay exponent_s 
  const float wd2 = 0.7;
  const float wd3 = 1.;
  
  // pulse shape componnt_s over a range of time 0 ns to 255 ns in 1 ns steps
  std::vector<float> nth(nbin_,0.0);   // zeroing HPD drift shape
  std::vector<float> ntp(nbin_,0.0);   // zeroing Binkley preamp shape
  std::vector<float> ntd(nbin_,0.0);   // zeroing Scintillator decay shape

  int i,j,k;
  float norm;

  // HPD starts at I and rises to 2I in thpd of time
  norm=0.0;
  for(j=0;j<thpd && j<nbin_;j++){
    nth[j] = 1.0 + j/thpd;
    norm += nth[j];
  }
  // normalize integrated current_ to 1.0
  for(j=0;j<thpd && j<nbin_;j++){
    nth[j] /= norm;
  }
  
  // Binkley shape over 6 time constant_s
  norm=0.0;
  for(j=0;j<6*tpre && j<nbin_;j++){
    ntp[j] = j*exp(-(j*j)/(tpre*tpre));
    norm += ntp[j];
  }
  // normalize pulse area to 1.0
  for(j=0;j<6*tpre && j<nbin_;j++){
    ntp[j] /= norm;
  }

// ignore stochastic variation of photoelectron emission
// <...>

// effective tile plus wave-length shifter decay time over 4 time constant_s
  int tmax = 6 * static_cast<int>(ts3);
 
  norm=0.0;
  for(j=0;j<tmax && j<nbin_;j++){
    ntd[j] = wd1 * exp(-j/ts1) + 
      wd2 * exp(-j/ts2) + 
      wd3 * exp(-j/ts3) ; 
    norm += ntd[j];
  }
  // normalize pulse area to 1.0
  for(j=0;j<tmax && j<nbin_;j++){
    ntd[j] /= norm;
  }
  
  int t1,t2,t3,t4;
  for(i=0;i<tmax && i<nbin_;i++){
    t1 = i;
    //    t2 = t1 + top*rand;
    // ignoring jitter from optical path length
    t2 = t1;
    for(j=0;j<thpd && j<nbin_;j++){
      t3 = t2 + j;
      for(k=0;k<4*tpre && k<nbin_;k++){       // here "4" is set deliberately,
 t4 = t3 + k;                         // as in test fortran toy MC ...
 if(t4<nbin_){                         
   int ntb=t4;                        
   nt_[ntb] += ntd[i]*nth[j]*ntp[k];
        }
      }
    }
  }
  
  // normalize for 1 GeV pulse height
  norm = 0.;
  for(i=0;i<nbin_;i++){
    norm += nt_[i];
  }

  //cout << " Convoluted SHAPE ==============  " << endl;
  for(i=0; i<nbin_; i++){
    nt_[i] /= norm;
    //  cout << " shape " << i << " = " << nt_[i] << endl;   
  }

}

void HcalShape::display

(

)

const [inline]

Definition at line 27 of file HcalShape.h.

{}

double HcalShape::operator()

(

double

time

)

const [virtual]

Implements CaloVShape.

Definition at line 118 of file HcalShape.cc.

References nbin_, and nt_.

{

  // return pulse amplitude for request time in ns
  int jtime = static_cast<int>(time_+0.5);
  if(jtime>=0 && jtime<nbin_){
    return nt_[jtime];
  } else {
    return 0.0;
  }

}

double HcalShape::timeToRise

(

)

const [virtual]

Implements CaloVShape.

Definition at line 19 of file HcalShape.cc.

{
   return 0. ;
}

---

Member Data Documentation

int HcalShape::nbin_ [private]

Definition at line 33 of file HcalShape.h.

Referenced by computeShape(), and operator()().

std::vector<float> HcalShape::nt_ [private]

Definition at line 34 of file HcalShape.h.

Referenced by computeShape(), and operator()().