BeamTOF Class Reference

#include <BeamTOF.h>

Inheritance diagram for BeamTOF:

List of all members.

Public Member Functions
	BeamTOF (const char *version)
	~BeamTOF ()
JobCResult	Reco ()
JobCResult	Ana (const EDMEventHandle &evt)
void	Update (const CfgConfig &c)
Private Attributes
int	fDebug
int	fNEvent
TH1F *	beamtof
CMapTOF *	cmaptof
std::vector< const TOFDigit * >	tofdigits
std::vector< const TOFDigit * >	t0digits

---

Detailed Description

Definition at line 15 of file BeamTOF.h.

---

Constructor & Destructor Documentation

BeamTOF::BeamTOF

(

const char *

version

)

Definition at line 29 of file BeamTOF.cxx.

References beamtof, cmaptof, fDebug, CMapTOF::Instance(), name, and JobCModule::SetCfgVersion().

                                    : 
   JobCModule("BeamTOF")
{
  // Do this first to load default configuration
   this->SetCfgVersion(version);

   if (fDebug) {
     std::cout <<"Debug value is: " << fDebug << std::endl;
   }

  // Initialize the histograms
  char name[80], title[256];
  sprintf(title, "Time of Flight for Beam Tracks");
  sprintf(name, "beamtof");
  // Histos here
  beamtof = new TH1F(name, title, 100, 0, 30); //20 50

  // Init the connection map
  cmaptof = &CMapTOF::Instance();
}

BeamTOF::~BeamTOF

(

)

Definition at line 52 of file BeamTOF.cxx.

References beamtof.

{
  // Write out the histos
  if(beamtof->GetEntries() != 0) beamtof->Write();
  delete beamtof;
}

---

Member Function Documentation

JobCResult BeamTOF::Ana

(

const EDMEventHandle &

evt

)

[virtual]

Reimplemented from JobCModule.

Definition at line 61 of file BeamTOF.cxx.

References TOFDigit::ADC(), beamtof, EDMTriggerId::BeamTrig(), CENTERCHAN, TOFDigit::Channel(), EDMEventHeader::Event(), fDebug, fNEvent, EDMDataBucket::Get(), EDMEventHandle::Header(), JobCModule::kFailed, JobCModule::kPassed, EDMEventHandle::Raw(), TOFBarHit::SetChannel(), TOFBarHit::SetLowerADC(), TOFBarHit::SetLowerTime(), TOFBarHit::SetUpperADC(), TOFBarHit::SetUpperTime(), t0digits, TOFDigit::TDC(), tofdigits, and EDMEventHeader::TriggerId().

{
  // Just to make sure...
  fNEvent=evt.Header().Event(); 
  if(!evt.Header().TriggerId().BeamTrig()) {
    if(fDebug) std::cout << "Not a beam trigger" << std::endl;
    return JobCModule::kFailed;
  }
  
  // Clear the list and fill it
  tofdigits.clear();
  t0digits.clear();

  // Some useful information 
  std::cout << "***********************************" << std::endl;
  std::cout << "Starting event " << fNEvent <<std::endl;
  // Read in tof digits 
  try { 
    evt.Raw().Get("./tof", tofdigits);
  } catch(EDMException e) {
    std::cerr << "Failed to get TOF Digits for event " 
      << fNEvent << std::endl;
    return JobCModule::kFailed;
  }
  if(fDebug) {
    std::cout << "read in tof digits " << tofdigits.size() << std::endl; 
  }
  
  // Read in t0 digits
  try {
    evt.Raw().Get("./t0", t0digits);
  } catch(EDMException e) {
    std::cerr << "Failed to get t0 digits for event "
      << fNEvent << std::endl;
    return JobCModule::kFailed;
  }    
  if(fDebug) {
    std::cout << "read in t0 digits " << t0digits.size() << std::endl;
  }

  // I'm assuming all four pmts fire.
  // T00 = 500*, T01 = 501*, TBD = 503*
  double t0sum = 0;
  for(unsigned int j=0; j<t0digits.size(); j++) {
    if(t0digits[j]->Channel() < 5015 && t0digits[j]->Channel() > 5010) {
      switch(t0digits[j]->Channel()) {
        case 5011: //5011: 300 150 5031: 700 300
          if(t0digits[j]->ADC() > 300 || t0digits[j]->ADC() < 150) 
            return JobCModule::kFailed;
          break;
        case 5012: //5012: 275 125   5032: 600 250
          if(t0digits[j]->ADC() > 275 || t0digits[j]->ADC() < 125) 
            return JobCModule::kFailed;
          break;
        case 5013: //5013: 350 180 5033: 800 300
          if(t0digits[j]->ADC() > 350 || t0digits[j]->ADC() < 180) 
            return JobCModule::kFailed;
          break;
        case 5014: //5014: 360 150  5034: 800 300
          if(t0digits[j]->ADC() > 360 || t0digits[j]->ADC() < 150) 
            return JobCModule::kFailed;
          break;
        default:
          break;
      }    
      t0sum+=t0digits[j]->TDC();
      std::cout << "t0 fired " << j-4 << " times" << std::endl;
    }
  }
  
  // Divide sum by 4
  double t0time = t0sum/4;

  // Convert to ns using that damned connection map
  t0time = t0time / (1000/30); //T01: 30   TBD: 60 
  //cmaptof->TDCBinSize();
  //if(fDebug) 
  std::cout << "T0 time: " << t0time << std::endl;
  // Put the tofdigits in a map so they can be indexed by channel
  std::map<int, const TOFDigit*> tofmap;
  for(unsigned int j=0; j<tofdigits.size(); j++) {
    tofmap[tofdigits[j]->Channel()] = tofdigits[j];
  }  
  
  for(unsigned int upper=0;upper<tofdigits.size();upper++) {
    // Check to see if the upper PMT has a hit
    // Skip over lower hits and add those in
    // When doing upper hits, also do calibration bars.
    if(tofdigits[upper]->Channel() != 2501 
        && tofdigits[upper]->Channel() != 2502)
      if(tofdigits[upper]->Channel()>1000) continue; 

    // Use the loop to make sure we account for all the hits in the wall
    // Setup upper and lower pmt pointers.
    int upperchan=tofdigits[upper]->Channel();
    int lowerchan=upperchan+1000;
    
    if(fDebug>1) std::cout << "Upper TDC: " 
      << tofdigits[upper]->TDC() << std::endl;
    
    const TOFDigit* upperpmt=tofmap[upperchan];
    const TOFDigit* lowerpmt=tofmap[lowerchan];
    
    if(fDebug>1) std::cout << "Upper channel: " << upperchan 
      << " Lower channel: " << lowerchan << std::endl;
    
    // See if Upper or Lower Tube was hit:
    // If both upper and lower have no hit--continue to next bar
    if (upperpmt->TDC() > 2047 && lowerpmt->TDC() > 2047) continue;
    if(upperpmt->TDC() == 0 && lowerpmt->TDC() == 0) {
      if(fDebug) std::cout << "Both TDCs == 0!" << std::endl;
      break;
    }
    
    // Create a bar hit and use that
    TOFBarHit bar;
    bar.SetChannel(upperchan);
    // Don't care about ADC's so this'll be commented 
    // until it happens that we do 
    bar.SetUpperADC(upperpmt->ADC());
    bar.SetLowerADC(lowerpmt->ADC());
    
    // Check for TDC overflow
    if(upperpmt->TDC() > 2047) { // overflow
      bar.SetUpperTime(-10000.);
    } else {
      // Assertively do nothing.
    }

    if(lowerpmt->TDC() > 2047) { // overflow
      bar.SetLowerTime(-10000.);
    } // else { also do nothing }
 
    // Could probably use bar to hold times...  but since I don't care about
    // cleanliness in this module, I won't! ;)
     
    if(lowerpmt->TDC() > 2047 or upperpmt->TDC() > 2047) {
      if(fDebug) std::cout << "TDC value too large!" << std::endl;
      return JobCModule::kFailed;
    } else {
      // A few useful numbers
      //double diff = upperpmt->TDC() - lowerpmt->TDC();
      double sum = lowerpmt->TDC() + upperpmt->TDC();
      double toftime = sum/2 / (1000/50); // tof time in ns
      std::cout << "tof time: " << toftime << std::endl;
      double tof = 0;
      
      // Do things for the center channel
      std::cout << upperpmt->Channel() << std::endl;
      if(upperpmt->Channel() == CENTERCHAN) {
        if(CENTERCHAN == 320) {
          if(upperpmt->ADC() > 600 || upperpmt->ADC() < 400) {
            std::cout << "failing on upper: " << upperpmt->ADC() << std::endl;
            return JobCModule::kFailed;
          }
          if(lowerpmt->ADC() > 120 || lowerpmt->ADC() < 85) {
            std::cout << "failing on lower: " << lowerpmt->ADC() << std::endl;
            return JobCModule::kFailed;
          }
        }
        else if(CENTERCHAN == 326) {
          // 310 230
          if(upperpmt->ADC() > 380 || upperpmt->ADC() < 180) {
            std::cout << "failing on upper: " << upperpmt->ADC() << std::endl;
            return JobCModule::kFailed;
          }
          if(lowerpmt->ADC() > 145 || lowerpmt->ADC() < 90) {
            std::cout << "failing on lower: " << lowerpmt->ADC() << std::endl;
            return JobCModule::kFailed;
          }
        }
        tof = toftime - t0time;
        //if(fDebug)
        std::cout << "Time of flight: " << tof << std::endl;
        beamtof->Fill(tof);
      }
    }
  }
  
  return JobCModule::kPassed;
} 

JobCResult BeamTOF::Reco

(

)

Definition at line 245 of file BeamTOF.cxx.

References JobCModule::kPassed.

                         {
  return JobCModule::kPassed;
}

void BeamTOF::Update

(

const CfgConfig &

c

)

[virtual]

Implements CfgObserver.

Definition at line 252 of file BeamTOF.cxx.

References fDebug.

{
  c("Debug").Get(fDebug);
}

---

Member Data Documentation

TH1F* BeamTOF::beamtof [private]

Definition at line 30 of file BeamTOF.h.

Referenced by Ana(), BeamTOF(), and ~BeamTOF().

CMapTOF* BeamTOF::cmaptof [private]

Definition at line 31 of file BeamTOF.h.

Referenced by BeamTOF().

int BeamTOF::fDebug [private]

Definition at line 27 of file BeamTOF.h.

Referenced by Ana(), BeamTOF(), and Update().

int BeamTOF::fNEvent [private]

Definition at line 28 of file BeamTOF.h.

Referenced by Ana().

std::vector<const TOFDigit*> BeamTOF::t0digits [private]

Definition at line 33 of file BeamTOF.h.

Referenced by Ana().

std::vector<const TOFDigit*> BeamTOF::tofdigits [private]

Definition at line 32 of file BeamTOF.h.

Referenced by Ana().

---

The documentation for this class was generated from the following files:

• BeamTOF.h
• BeamTOF.cxx

---