DSTUtil::ToF Namespace Reference

Functions
int	BarToArray (int bar)
int	ArrayToBar (int array)
void	GetDiscHitMap (const MIPPEventSummary *evt, map< int, int > &hitMap)
int	AdjacentHit (const MIPPEventSummary *evt, int pmtid)
bool	GetCleanTracks (const MIPPEventSummary *evt, std::vector< int > &trklist)
bool	GetCleanTracks (const MIPPEventSummary *evt, std::map< int, bool > &trktofmap)
bool	GetCleanTracks (const MIPPEventSummary *evt, const MIPPVertexSummary *vtx, std::vector< int > &trklist)
bool	GetCleanTracks (const MIPPEventSummary *evt, const MIPPVertexSummary *vtx, std::map< int, bool > &trktofmap)
void	GetXTalkOffset (const MIPPToFSummary *tof, const MIPPEventSummary *evt, double &OffsetTop, double &ResTop, double &OffsetBot, double &ResBot)
void	CorrectedLL (const MIPPToFSummary *itof, const MIPPEventSummary *evt, double topoff, double, double botoff, double, double &tmeas, double *CorrLL)
void	GetDiscHitMap (const MIPPEventSummary *evt, std::map< int, int > &hitmap)

---

Function Documentation

int DSTUtil::ToF::BarToArray

(

int

bar

)

Definition at line 1119 of file DSTUtil.cxx.

References TOFUtil::BarToArray().

Referenced by AdjacentHit(), TOFCalib2::Ana(), Ana(), TOFTempCorrector::Apply(), fcn(), GetData(), TOFTempCorrector::GraphAna(), InitParams(), main(), TOFTempCorrector::PreAna(), PreAna(), and ProcEvent().

{
  return TOFUtil::BarToArray(bar);
}

int DSTUtil::ToF::ArrayToBar

(

int

array

)

Definition at line 1126 of file DSTUtil.cxx.

References TOFUtil::ArrayToBar().

Referenced by Book(), DoGlobalFits(), FillAdjHitTable(), FillDistHists(), FillHistos(), GetData(), GetSeedOffset(), GetSeedSlew(), GetSeedSOL(), GetSeedTempCoeffs(), LoadTOFAdjHitTable(), main(), MakeEvtHistos(), MakeGlobalPlots(), MakeResolPlots(), SaveConsts(), and SetBranches().

{
  return TOFUtil::ArrayToBar(array);
}

void DSTUtil::ToF::GetDiscHitMap	(	const MIPPEventSummary *	evt,
		map< int, int > &	hitMap
	)

Definition at line 1134 of file DSTUtil.cxx.

References CMapTOF::BartoDisc(), cmap, MIPPEventSummary::GetToFPMT(), MIPPToFPMTSummary::id, CMapTOF::Init(), CMapTOF::Instance(), MIPPEventSummary::NToFPMT(), and MIPPEventSummary::run.

Referenced by GetXTalkOffset(), and main().

{
  hitMap.clear();

  int id, id2, crate, crate2, slot, slot2, chan, chan2;
  int ihist;

  CMapTOF& cmap = CMapTOF::Instance();
  cmap.Init((unsigned int)evt->run);

  for (int ipmt=0; ipmt < evt->NToFPMT(); ++ipmt) {
    id = evt->GetToFPMT(ipmt)->id;
    cmap.BartoDisc(id,crate,slot,chan);
    
    // determine surround hit pattern for this channel
    ihist = 0;
    for (int ipmt2=0; ipmt2 < evt->NToFPMT(); ++ipmt2) {
      id2 = evt->GetToFPMT(ipmt2)->id;
      if (id == id2) continue;
      cmap.BartoDisc(id2,crate2,slot2,chan2);
      if (crate == crate2 && slot == slot2) {
    if (chan2 == chan+2)  ihist |= 1;
    if (chan2 == chan+1)  ihist |= 2;
    if (chan2 == chan-1)  ihist |= 4;
    if (chan2 == chan-2)  ihist |= 8;
      }
    }
    hitMap[id] = ihist;
  }
}

int DSTUtil::ToF::AdjacentHit	(	const MIPPEventSummary *	evt,
		int	pmtid
	)

Definition at line 1171 of file DSTUtil.cxx.

References BarToArray(), MIPPEventSummary::GetToFPMT(), ibar, npmt, and MIPPEventSummary::NToFPMT().

Referenced by TOFCalib2::Ana(), Ana(), GetCleanTracks(), TOFTempCorrector::GraphAna(), main(), and TOFTempCorrector::PreAna().

                                                                  {

  int npmt = evt->NToFPMT();
  int ibar = DSTUtil::ToF::BarToArray(pmtid % 1000);

  int nhit = 0;

  for (int ipmt=0; ipmt< npmt; ++ipmt)  // loop over pmts
  {
    MIPPToFPMTSummary* pmt;
    pmt = evt->GetToFPMT(ipmt);
    int jbar = DSTUtil::ToF::BarToArray(pmt->id % 1000);

    if (ibar > 0 && ibar < 53) {
      if (jbar == ibar-1) ++nhit;
      if (jbar == ibar+1) ++nhit;
    }
    if (ibar == 0)
      if (jbar == 1) ++nhit;
    if (ibar == 53)
      if (jbar == 52) ++nhit;
  }
  
  return nhit;

}

bool DSTUtil::ToF::GetCleanTracks	(	const MIPPEventSummary *	evt,
		std::vector< int > &	trklist
	)

Definition at line 1200 of file DSTUtil.cxx.

References AdjacentHit(), MIPPEventSummary::GetToF(), MIPPEventSummary::GetToFPMT(), MIPPToFPMTSummary::id, MIPPEventSummary::NToF(), and MIPPEventSummary::NToFPMT().

Referenced by NuMIudst::Fill(), NuMIPID::Fill(), main(), ProcEvent(), and ProcPrevEvent().

{

  trklist.clear();

  if (! evt->NToFPMT() || ! evt->NToF()) return false;

  std::map<int,MIPPToFPMTSummary*> pmtlist;
  int id;

  //  cout << "there are " << evt->NToFPMT() << " ToF PMT summaries" << endl;
  for (int ipmt=0; ipmt < evt->NToFPMT(); ++ipmt) {
    id = evt->GetToFPMT(ipmt)->id;
    pmtlist[id] = evt->GetToFPMT(ipmt);
  }

  //  cout << "there are " << evt->NToF() << " ToF summaries" << endl;
  for (int itof=0; itof<evt->NToF(); ++itof) {
    MIPPToFSummary* tof = evt->GetToF(itof);
    if (! pmtlist[tof->tpmt]) continue;
    if (! pmtlist[tof->bpmt]) continue;

    if (pmtlist[tof->tpmt]->ntrk > 1 || 
    pmtlist[tof->bpmt]->ntrk > 1) continue;
    
    if (DSTUtil::ToF::AdjacentHit(evt,tof->tpmt) ||
    DSTUtil::ToF::AdjacentHit(evt,tof->bpmt)) continue;

    trklist.push_back(tof->itrk);
  }

  if (trklist.empty()) return false;

  return true;
}

bool DSTUtil::ToF::GetCleanTracks	(	const MIPPEventSummary *	evt,
		std::map< int, bool > &	trktofmap
	)

Definition at line 1239 of file DSTUtil.cxx.

References AdjacentHit(), MIPPEventSummary::GetToF(), MIPPEventSummary::GetToFPMT(), MIPPToFPMTSummary::id, MIPPEventSummary::NToF(), and MIPPEventSummary::NToFPMT().

{

  trktofmap.clear();

  if (! evt->NToFPMT() || ! evt->NToF()) return false;

  std::map<int,MIPPToFPMTSummary*> pmtlist;
  
  for (int ipmt=0; ipmt < evt->NToFPMT(); ++ipmt) {
    pmtlist[evt->GetToFPMT(ipmt)->id] = evt->GetToFPMT(ipmt);
  }

  for (int itof=0; itof<evt->NToF(); ++itof) {
    MIPPToFSummary* tof = evt->GetToF(itof);
    if (! pmtlist[tof->tpmt]) continue;
    if (! pmtlist[tof->bpmt]) continue;

    if (pmtlist[tof->tpmt]->ntrk > 1 || 
    pmtlist[tof->bpmt]->ntrk > 1) continue;
    
    if (DSTUtil::ToF::AdjacentHit(evt,tof->tpmt) ||
    DSTUtil::ToF::AdjacentHit(evt,tof->bpmt)) continue;

    trktofmap[tof->itrk]= true;
  }

  return true;
}

bool DSTUtil::ToF::GetCleanTracks	(	const MIPPEventSummary *	evt,
		const MIPPVertexSummary *	vtx,
		std::vector< int > &	trklist
	)

Definition at line 1272 of file DSTUtil.cxx.

References AdjacentHit(), MIPPEventSummary::GetToF(), MIPPEventSummary::GetToFPMT(), MIPPToFPMTSummary::id, MIPPEventSummary::NToF(), MIPPEventSummary::NToFPMT(), MIPPVertexSummary::ntrk, and MIPPVertexSummary::trkindex.

{

  trklist.clear();

  if (! evt->NToFPMT() || ! evt->NToF()) return false;

  std::map<int,MIPPToFPMTSummary*> pmtlist;
  
  for (int ipmt=0; ipmt < evt->NToFPMT(); ++ipmt)
    pmtlist[evt->GetToFPMT(ipmt)->id] = evt->GetToFPMT(ipmt);
  
  for (int itof=0; itof<evt->NToF(); ++itof) {
    MIPPToFSummary* tof = evt->GetToF(itof);
    if (! pmtlist[tof->tpmt]) continue;
    if (! pmtlist[tof->bpmt]) continue;
    
    if (pmtlist[tof->tpmt]->ntrk > 1 || 
    pmtlist[tof->bpmt]->ntrk > 1) continue;

    if (DSTUtil::ToF::AdjacentHit(evt,tof->tpmt) ||
    DSTUtil::ToF::AdjacentHit(evt,tof->bpmt)) continue;
    
    int ioff = vtx->trkindex;

    for (int itrk=0; itrk < vtx->ntrk; ++itrk)
      if (ioff+itrk == tof->itrk) {
    trklist.push_back(tof->itrk);
    break;
      }
    
  }

  if (trklist.empty()) return false;

  return true;
}

bool DSTUtil::ToF::GetCleanTracks	(	const MIPPEventSummary *	evt,
		const MIPPVertexSummary *	vtx,
		std::map< int, bool > &	trktofmap
	)

Definition at line 1314 of file DSTUtil.cxx.

References AdjacentHit(), MIPPEventSummary::GetToF(), MIPPEventSummary::GetToFPMT(), MIPPToFPMTSummary::id, MIPPEventSummary::NToF(), MIPPEventSummary::NToFPMT(), MIPPVertexSummary::ntrk, and MIPPVertexSummary::trkindex.

{

  trktofmap.clear();

  if (! evt->NToFPMT() || ! evt->NToF()) return false;

  std::map<int,MIPPToFPMTSummary*> pmtlist;
  
  for (int ipmt=0; ipmt < evt->NToFPMT(); ++ipmt)
    pmtlist[evt->GetToFPMT(ipmt)->id] = evt->GetToFPMT(ipmt);
  
  for (int itof=0; itof<evt->NToF(); ++itof) {
    MIPPToFSummary* tof = evt->GetToF(itof);
    if (! pmtlist[tof->tpmt]) continue;
    if (! pmtlist[tof->bpmt]) continue;
    
    if (pmtlist[tof->tpmt]->ntrk > 1 || 
    pmtlist[tof->bpmt]->ntrk > 1) continue;

    if (DSTUtil::ToF::AdjacentHit(evt,tof->tpmt) ||
    DSTUtil::ToF::AdjacentHit(evt,tof->bpmt)) continue;
    
    int ioff = vtx->trkindex;

    for (int itrk=0; itrk < vtx->ntrk; ++itrk)
      if (ioff+itrk == tof->itrk) {
    trktofmap[tof->itrk]=true;
    break;
      }
    
  }

  return true;
}

void DSTUtil::ToF::GetXTalkOffset	(	const MIPPToFSummary *	tof,
		const MIPPEventSummary *	evt,
		double &	OffsetTop,
		double &	ResTop,
		double &	OffsetBot,
		double &	ResBot
	)

Definition at line 1354 of file DSTUtil.cxx.

References adjhit, MIPPToFSummary::bpmt, GetDiscHitMap(), TOFAdjHit::Init(), TOFAdjHit::Instance(), TOFAdjHit::Offset(), TOFAdjHit::Res(), and MIPPToFSummary::tpmt.

{
  int chanidTop = tof->tpmt;
  int chanidBot = tof->bpmt;
  map<int,int> hitMap;
  DSTUtil::ToF::GetDiscHitMap(evt,hitMap);
  int hitsTop = hitMap[chanidTop];
  int hitsBot = hitMap[chanidBot];
  TOFAdjHit& adjhit = TOFAdjHit::Instance();
  adjhit.Init((unsigned int)14500);
  OffsetTop = adjhit.Offset(chanidTop,hitsTop);
  ResTop = adjhit.Res(chanidTop,hitsTop);
  OffsetBot = adjhit.Offset(chanidBot,hitsBot);
  ResBot = adjhit.Res(chanidBot,hitsBot);
}

void DSTUtil::ToF::CorrectedLL	(	const MIPPToFSummary *	itof,
		const MIPPEventSummary *	evt,
		double	topoff,
		double	topres,
		double	botoff,
		double	botres,
		double &	tmeas,
		double *	CorrLL
	)

Definition at line 1375 of file DSTUtil.cxx.

References TOFBarResol::BotRes(), MIPPToFSummary::bpmt, MIPPToFSummary::exptof, TOFBarResol::Init(), TOFBarResol::Instance(), PIDDef::kNPID, offset, MIPPToFSummary::res, MIPPEventSummary::run, MIPPToFSummary::tbot, MIPPToFSummary::tmeas, TOFBarResol::TopRes(), MIPPToFSummary::tpmt, and MIPPToFSummary::ttop.

{
  if(itof->ttop < 0 && itof->tbot < 0) {
    tmeas = -1e+9;
    for(int i=0; i<PIDDef::kNPID; ++i){
      CorrLL[i] = -1e+9;
    }
  }
  else{
    TOFBarResol& barRes = TOFBarResol::Instance();

    barRes.Init(evt->run);

    double offset = 0;

    double sumw = 0;

    double topresx = barRes.TopRes(itof->tpmt);
    double botresx = barRes.BotRes(itof->bpmt);

    if(itof->ttop > 0 ) {
      offset += topoff/topresx;
      sumw += 1./topresx;
    }
    if(itof->tbot > 0 ) {
      offset += botoff/botresx;
      sumw += 1./botresx;
    }

    offset = offset/sumw;

    double tmeasured = itof->tmeas + offset;

    for(int j=0; j<PIDDef::kNPID; ++j){
      double LL = (tmeasured - itof->exptof[j])/itof->res;
      LL *= -0.5*LL;
      LL -=log(2*TMath::Pi()*itof->res);
      CorrLL[j] = LL;

    }

    tmeas = tmeasured;
  }

}

void DSTUtil::ToF::GetDiscHitMap	(	const MIPPEventSummary *	evt,
		std::map< int, int > &	hitmap
	)

---