SimPixelTimeBucket.cxx

Go to the documentation of this file.

//
//
// SimPixelTimeBucket
//
// Description: This class holds all the data relevant to a
// single PMT pixel (all spots) at a given instant in time (a bucket).
// 
// Data members:
// fNSpots       number of spots in the pixel (1 or 8 typically)
// fDigiPE       a list of the digiPE that should have hit this pixel
// fDigiPEXtalk  a list of the digiPE that DID hit this pixel, after crosstalk.
// fTruth        a set of bits representing what happened.
// fPE           the number of digiPE that should have hit this pixel
// fPEXtalk      the number of digiPE that did hit this pixel, after crosstalk
// fCharge       the charge on the anode in this bucket, in Munits
// fTime         the earilest time of all the DigiPE in this bucket.

#include "DetSim/SimPixelTimeBucket.h"
#include "MessageService/MsgService.h"
#include "Conventions/Munits.h"

CVSID("$Id: SimPixelTimeBucket.cxx,v 1.11 2008/11/18 17:31:02 rhatcher Exp $");

ClassImp(SimPixelTimeBucket)

int SimPixelTimeBucket::fsTotalNumber = 0;

SimPixelTimeBucket::SimPixelTimeBucket()
{
  // Do not use!  Only for I/O.
  fsTotalNumber++;
}


SimPixelTimeBucket::SimPixelTimeBucket( Int_t nSpots ) :
  fNSpots(nSpots),
  fDigiPE(fNSpots+1),
  fDigiPEXtalk(fNSpots+1),
  fPE(fNSpots+1),
  fPEXtalk(fNSpots+1),
  fTruth(DigiSignal::kUnknown),
  fCharge(0),
  fTime(1e9)
{
  fsTotalNumber++;

  // Fill with zeros.
  for(int i=0;i<=fNSpots;i++) {
    fPE[i]=0;
    fPEXtalk[i]=0;
  }
}

SimPixelTimeBucket::~SimPixelTimeBucket()
{
  fsTotalNumber--;
}

void  SimPixelTimeBucket::SetCharge( Float_t charge )
{
  fCharge = charge;
  //if(charge<-100*Munits::fC) {
  //  cout << "Charge got set to " << fCharge/Munits::fC << endl;
  // }
}

void  SimPixelTimeBucket::AddCharge( Float_t charge )
{
  fCharge += charge;
  //if(fCharge<-100*Munits::fC) {
  //  cout << "Charge got added to " << fCharge/Munits::fC << " with increment " << charge/Munits::fC << endl;
  //}
}



Float_t SimPixelTimeBucket::GetTotalPE()
{
  float tot = 0;
  for(int i=1;i<=fNSpots;i++) {
    tot += fPE[i];
  }
  return tot;
}

Float_t SimPixelTimeBucket::GetTotalPEXtalk()
{
  float tot = 0;
  for(int i=1;i<=fNSpots;i++) {
    tot += fPEXtalk[i];
  }
  return tot;
}



void SimPixelTimeBucket::AddDigiPE(const DigiPE* digipe, Int_t spot)
{
  //
  // Add DigiPE to the bucket.  Use the provided spot number instead of the 
  // PlexPixelSpotId in the digipe, because the PlexPixelSpotID is simply wrong,
  // if accurately wrong.
  //
  
  //Int_t spot = digipe->GetPixelSpotId().GetSpot();
  if( (spot > fNSpots) || (spot < 1) ) {
    // Problem!
    MSG("DetSim",Msg::kWarning) << "AddDigiPE() Invalid Pixel Spot number!" << endl;
    return;
  }
  if(digipe->GetTime() < fTime) fTime = digipe->GetTime();
  
  // This adds the hit info and the pe pointer.
  fPE[spot] += 1;                    // True number of PE.. don't touch this again.
  
  fDigiPE[spot].insert(std::pair<Double_t,const DigiPE*>(digipe->GetTime(), digipe));
  
  // This is a real hit.
  
  fTruth |= digipe->GetSource();
}


DigiSignal* SimPixelTimeBucket::CreateSignal() 
{
  // Creates a new digiSignal with all the attendant data.
  // Caller becomes responsible for ownership.
  //
  DigiSignal* s = new DigiSignal();
  // Set up the charge
  s->SetCharge(fCharge);
  s->SetTruth(fTruth);

  // Add the PEs.
  for(UInt_t i = 0; i<fDigiPEXtalk.size(); i++) {
    PeList_t::iterator it = fDigiPEXtalk[i].begin();
    for( ; it!= fDigiPEXtalk[i].end() ; it++) {
      s->AddDigiPE(it->second);
    }
  }

  return s;
}

---