NtpSRNDAPPlaneHistory.cxx

Go to the documentation of this file.

// Program name: NtpSRNDAPPlaneHistory.cxx              
//                                                  
// Package: CandNtupleSR/Module
//                                                                    
// Purpose: store all strips from one plane
//          in the same snarl for AP calculation
//                                                                    
// Contact: Tobias Raufer

#include <iostream>

#include "TMath.h"

#include "CandNtupleSR/Module/NtpSRNDAPPlaneHistory.h"
#include "MessageService/MsgService.h"

using std::cout;
using std::endl;
using std::multimap;

  CVSID("$Id: NtpSRNDAPPlaneHistory.cxx,v 1.1 2005/11/05 23:53:51 schubert Exp $");

//......................................................................
NtpSRNDAPPlaneHistory::NtpSRNDAPPlaneHistory() {
  this->Reset();
  fPlane=-1;
}

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

NtpSRNDAPPlaneHistory::NtpSRNDAPPlaneHistory(Int_t plane)
{
  MSG("NtpSRNDAPPlaneHistory",Msg::kVerbose) 
    <<"Running NtpSRNDAPPlaneHistory constructor..."<<std::endl;
  this->Reset();
  fPlane = plane;
}

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

NtpSRNDAPPlaneHistory::NtpSRNDAPPlaneHistory(Int_t plane, 
                                   multimap<Double_t,Double_t> stripList) :
  fPlane(plane),
  fStripList(stripList)
{
  MSG("NtpSRNDAPPlaneHistory",Msg::kVerbose)    
    <<"Running NtpSRNDAPPlaneHistory constructor for plane "
    << plane <<std::endl;
}

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

NtpSRNDAPPlaneHistory::~NtpSRNDAPPlaneHistory() {
  MSG("NtpSRNDAPPlaneHistory",Msg::kVerbose)   
    <<"Running NtpSRNDAPPlaneHistory destructor..."<<std::endl;
  this->Reset();
}

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

void NtpSRNDAPPlaneHistory::Reset() {
  fStripList.clear();
}

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

void NtpSRNDAPPlaneHistory::AddStrip(Double_t time, Double_t adc) {
  
  fStripList.insert(std::pair<Double_t,Double_t>(time,adc));
}

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

Int_t NtpSRNDAPPlaneHistory::GetNStrips() {
  return fStripList.size();
}

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

Double_t NtpSRNDAPPlaneHistory::APPrediction(Double_t time) {

  Double_t result = 0;

  multimap<Double_t,Double_t>::iterator pos;
  for (pos = fStripList.begin(); pos!=fStripList.end(); ++pos) {
    if (pos->first  < (time-100e-9)) {
      Double_t deltaT = time-pos->first;
      result += pos->second*(1.67e-2*TMath::Exp(-deltaT/343e-9) 
                             + 8.5e-4*TMath::Exp(-deltaT/2.4e-6));
    }
  }

  MSG("NDAPHistory",Msg::kVerbose) << "AP prediction for time "
                                << time << ": " << result << endl;
  
  return result;
}

---