DSTUtil Namespace Reference

Namespaces
namespace	Bflux
namespace	Calo
namespace	Ckov
namespace	Global
namespace	GoodRun
namespace	NuMI
namespace	Prsc
namespace	RICH
namespace	RootUtil
namespace	Target
namespace	ToF
namespace	Trig
Functions
int	HavePileUpCommon (const MIPPEventSummary *evt)
	Functions which can be used to select (reject) pile-up events independent of beam particle momentum.
int	HavePileUp60 (const MIPPEventSummary *evt)
	Functions which can be used to select (reject) pile-up events in 60 GeV/c momentum.
int	HavePileUp120 (const MIPPEventSummary *evt)
	Functions which can be used to select (reject) pile-up events in 120 GeV/c momentum.
bool	GoodBeam (const MIPPRunSummary *rs, const MIPPEventSummary *evt, const MIPPBeamTrkSummary *fBeam)
bool	GoodVertex (const MIPPRunSummary *runsum, const MIPPVertexSummary *vtx, float rWin, float zWin)
bool	GoodTrack (const MIPPTrackSummary *trk)
void	SigHandle (int sig)
	Default signal handler.
TTree *	GetTreeByName (TFile *const f, std::string const fTreeDir, std::string const fTreeName)
	Given a MIPP histogram file, return a tree with name fTreeName in the first directory with name starting with fTreeDir.
TTree *	GetEvtTree (TFile *f)
	Given a MIPP histogram file, return DST tree.
TTree *	GetSpillTree (TFile *f)
	Given a MIPP histogram file, return DST tree.
void	GetFileList (int argc, char **argv, std::vector< std::string > &list)
	Parses command-line arguments to extract names of valid ROOT files.
void	AddToFileList (char *arg, std::vector< std::string > &list)
	Parses arg and add to list if it is a ROOT file or add all ROOT fileslisted in file arg to list.
MIPPRunSummary *	GetRunSummary (TFile *f)
	Gets the first entry from run summary tree.
MIPPRunSummary *	GetRunSummaryVerbose (TFile *f)
	like DSTUtil::GetRunSummary(TFile* f), but print some info about this run to stdout.
MIPPRunSummary *	GetRunSummary (TFile *f, const int &subrun)
	Gets entry for subrun from run summary tree.
TTree *	GetMCMatchTree (TFile *f)
	Given a MIPP histogram file, return MCMatch tree.
TTree *	GetVeeTree (TFile *f, string &ntuple)
TTree *	GetVeeTree (TFile *f, std::string &ntuple)

---

Function Documentation

int DSTUtil::HavePileUpCommon

(

const MIPPEventSummary *

evt

)

Functions which can be used to select (reject) pile-up events independent of beam particle momentum.

Definition at line 39 of file DSTUtil.cxx.

References MIPPBeamTrkSummary::avgt, MIPPEventSummary::GetBTrk(), MIPPEventSummary::GetT0(), RIMCConfig::kMCRunOffset, MIPPEventSummary::NBTrk(), MIPPEventSummary::NT0(), MIPPT0Summary::Qtot(), and MIPPEventSummary::run.

Referenced by HavePileUp120(), and HavePileUp60().

                                                        {
  // No pile up on MC events
  if (evt->run >= RIMCConfig::kMCRunOffset) return 0;
  
  // Avoid multiple beam tracks
  if (evt->NBTrk() > 1) return 1;
  
  // Make sure beam track time is ok, if we have a beam track at all
  MIPPBeamTrkSummary* trk = evt->GetBTrk(0);
  if (trk && (trk->avgt < -5 || trk->avgt > 20)) return 2;
  
  // Check T0 summed charge
  float maxT0ADC = 1.8 * 4 * 3;// 3 counters with 4 tubes each
  float minT0ADC = 0.5 * 4 * 3;// 3 counters with 4 tubes each
  if (evt->NT0() == 5) {
    float totADC = 0;
    totADC += evt->GetT0(0)->Qtot(); // T00
    totADC += evt->GetT0(1)->Qtot(); // T01
    totADC += evt->GetT0(3)->Qtot(); // TBD
    
    if (totADC > maxT0ADC || totADC < minT0ADC) return 3; 
  }
  return 0;
}

int DSTUtil::HavePileUp60

(

const MIPPEventSummary *

evt

)

Functions which can be used to select (reject) pile-up events in 60 GeV/c momentum.

Definition at line 65 of file DSTUtil.cxx.

References HavePileUpCommon().

Referenced by main().

                                                    {
  return HavePileUpCommon(evt);
}

int DSTUtil::HavePileUp120

(

const MIPPEventSummary *

evt

)

Functions which can be used to select (reject) pile-up events in 120 GeV/c momentum.

Definition at line 70 of file DSTUtil.cxx.

References MIPPEventSummary::bckov1adc, MIPPEventSummary::bckov1tdc, MIPPEventSummary::bckov2adc, MIPPEventSummary::bckov2tdc, MIPPEventSummary::GetT0(), MIPPEventSummary::GetTrk(), HavePileUpCommon(), MIPPEventSummary::NT0(), MIPPEventSummary::NTrk(), MIPPEventSummary::t0, and t0.

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

                                                     {
  // Check common PileUp criteria
  int fCommon = HavePileUpCommon(evt);
  if(fCommon != 0) return fCommon;

  // Check beam Cherenkov summed charge
  if (evt->bckov1adc[0] + evt->bckov2adc[0] > 330 * 2) return 4;

  // Check T0 time differences
  float minT0diff[5] = {-1.5, -0.5, -0.2, -6, -1};
  float maxT0diff[5] = {1.5, 0.5, 0.25, 7, 0.8};
  for (int i = 0; i < evt->NT0(); ++i) {
    MIPPT0Summary* t0 = evt->GetT0(i);
    if (t0->timeavg - evt->t0 < minT0diff[i] || 
    t0->timeavg - evt->t0 > maxT0diff[i]) return 5;
  }

  // Check BCkov TDC's
  // Require upper inner PMT have a valid signal
  if (evt->bckov1tdc[0] < 800 || evt->bckov1tdc[0] > 1040) return 6;
  
  // Require upper outer PMT have no signal
  if (evt->bckov1tdc[1] > 1 && evt->bckov1tdc[1] < 2000) return 6;

  // Require downstream inner PMT have a valid signal
  if (evt->bckov1tdc[0] < 700 || evt->bckov2tdc[0] > 950) return 6;
  
  // Require downstream outer PMT have no signal
  if (evt->bckov1tdc[1] > 1 && evt->bckov1tdc[1] < 2000) return 6;

  // Check for out-of-time tracks
  float minT = -15;
  float maxT = 55;
  for (int i = 0; i < evt->NTrk(); ++i) {
    MIPPTrackSummary* trk = evt->GetTrk(i);
    if (trk->avgt < minT || trk->avgt > maxT) return 7;
  }

  return 0;
}

bool DSTUtil::GoodBeam	(	const MIPPRunSummary *	rs,
		const MIPPEventSummary *	evt,
		const MIPPBeamTrkSummary *	fBeam
	)

Definition at line 121 of file DSTUtil.cxx.

References MIPPBeamTrkSummary::dxdzbc3, MIPPBeamTrkSummary::dydzbc3, MIPPEventSummary::GetBTrk(), MIPPEventSummary::NBTrk(), MIPPRunSummary::tgta, MIPPRunSummary::tgtPos, and MIPPBeamTrkSummary::xbc3.

{
  if (rs->tgta > 1000) return false;     // do not apply it for NuMI target

  int beamPosFlag = 0;
  for(int ibeam=0; ibeam<evt->NBTrk(); ++ibeam) {

    fBeam = evt->GetBTrk(ibeam);

    // Recalculate beam track position at center of target
    double dzz = rs->tgtPos[2] - fBeam->xbc3[2];   
    double dxpos = (fBeam->xbc3[0] + fBeam->dxdzbc3 * dzz) - rs->tgtPos[0];
    double dypos = (fBeam->xbc3[1] + fBeam->dydzbc3 * dzz) - rs->tgtPos[1];
    double xyp = sqrt(dxpos*dxpos + (dypos-0.1)*(dypos-0.1));
    if(xyp < 1.5) beamPosFlag++;
  } // end ibeam loop

  if(beamPosFlag > 0) return true;
  return false;
}

bool DSTUtil::GoodVertex	(	const MIPPRunSummary *	runsum,
		const MIPPVertexSummary *	vtx,
		float	rWin = 1.,
		float	zWin = 1.
	)

Definition at line 152 of file DSTUtil.cxx.

References MIPPVertexSummary::intrkindex, MIPPVertexSummary::ntrk, runsum, MIPPRunSummary::tgtPos, MIPPRunSummary::tgtSize, and MIPPVertexSummary::x.

{
  if (!vtx) return false;
  if (vtx->intrkindex < 0) return false; // a good vertex requires a beam track
  // if no beam track, then this could be a V0

  if (vtx->ntrk < 2) return false;

  const float* pos = vtx->x;
  float dx = pos[0] - runsum->tgtPos[0];
  float dy = pos[1] - runsum->tgtPos[1];
  float dr = sqrt(dx*dx+dy*dy);
  float dz = pos[2] - runsum->tgtPos[2];

  if (fabs(dr) > rWin*runsum->tgtSize[0]) return false;
  if (fabs(dz) > zWin*runsum->tgtSize[2]) return false;

  return true;
}

bool DSTUtil::GoodTrack

(

const MIPPTrackSummary *

trk

)

Definition at line 178 of file DSTUtil.cxx.

References MIPPTrackSummary::gof, and MIPPTrackSummary::ntpchit.

Referenced by ProcEvent().

{
  if (trk->gof < 1.e-3) return false;
  if (trk->ntpchit < 20 || trk->ntpchit >= 100) return false;
  return true;
}

void DSTUtil::SigHandle

(

int

sig

)

Default signal handler.

Definition at line 188 of file DSTUtil.cxx.

References gStop.

Referenced by AddToFileList(), GetFileList(), and main().

{
  ++gStop;
  if (gStop > 2) abort();
  cout << "Stopping on signal " << sig << endl;
}

TTree * DSTUtil::GetTreeByName	(	TFile *const	f,
		std::string const	fTreeDir,
		std::string const	fTreeName
	)

Given a MIPP histogram file, return a tree with name fTreeName in the first directory with name starting with fTreeDir.

(Don't have to know if we deal with FillEventSummary_pass4 or FillEventSummary_mc, ...)

Definition at line 860 of file DSTUtil.cxx.

References tree.

Referenced by GetEvtTree(), GetMCMatchTree(), GetRunSummary(), GetSpillTree(), GetVeeTree(), and main().

{
  TList* l = f->GetListOfKeys();
  for (int i = 0; i < l->GetEntries(); ++i) {
    string fNm(l->At(i)->GetName());
    if (fNm.substr(0,fTreeDir.size()) == fTreeDir) {
      fNm += "/"+fTreeName;
      TTree* tree = dynamic_cast<TTree*>(f->Get(fNm.c_str()));
      if (tree) return tree;
    }
  }
  cout << "Couldn't find " << fTreeName << " in " << f->GetName()
       << endl;
  return 0;
}

TTree * DSTUtil::GetEvtTree

(

TFile *

f

)

Given a MIPP histogram file, return DST tree.

This function will return the tree from the first DST folder that it finds

Definition at line 882 of file DSTUtil.cxx.

References GetTreeByName().

Referenced by TOFTempCorrector::init(), main(), ValiMC::SetupDSTTree(), and GlobalPid::SetupDSTTree().

{
  return GetTreeByName(f, "FillEventSummary", "fNtTree");
}

TTree * DSTUtil::GetSpillTree

(

TFile *

f

)

Given a MIPP histogram file, return DST tree.

This function will return the tree from the first DST folder that it finds

Definition at line 891 of file DSTUtil.cxx.

References GetTreeByName().

Referenced by DSTUtil::Bflux::GetBeamFluxForRun(), and main().

{
  return GetTreeByName(f, "FillEventSummary", "fSpillTree");
}

void DSTUtil::GetFileList	(	int	argc,
		char **	argv,
		std::vector< std::string > &	list
	)

Parses command-line arguments to extract names of valid ROOT files.

The valid arguments are either ROOT file name or name of file which contains ROOT file names

Definition at line 901 of file DSTUtil.cxx.

References AddToFileList(), and SigHandle().

Referenced by main(), and ParseArgs().

{
  signal(SIGINT, SigHandle);
  //  struct stat buf;

  for (int i = 1; i < argc; ++i) {
    AddToFileList(argv[i],list);
  }
}

void DSTUtil::AddToFileList	(	char *	arg,
		std::vector< std::string > &	list
	)

Parses arg and add to list if it is a ROOT file or add all ROOT fileslisted in file arg to list.

Definition at line 916 of file DSTUtil.cxx.

References SigHandle().

Referenced by GetFileList(), and main().

{
  signal(SIGINT, SigHandle);
  struct stat buf;

  int len = strlen(arg);
  if (len > 5 && strcmp(".root", &arg[len-5]) == 0) {
    // Should be a valid ROOT file name
    if (strncmp("root:/", arg, 6) == 0 ||
        strncmp("dcap:/", arg, 6) == 0 ||
        stat(arg, &buf) == 0) {
      cout << "Adding file " << arg << endl;
      list.push_back(arg);
    }
    else cout << "File " << arg << " doesn't exist?" << endl;
  }
  else {
    // This must be a text file containing ROOT file names
    ifstream ifi(arg);
    std::string str;
    while (ifi >> str) {
      if (str.size() < 5) continue;
      if (str.substr(str.size() - 5) != ".root") continue;
      cout << "Adding file " << str << endl;
      list.push_back(str);
    }
  }
}

MIPPRunSummary * DSTUtil::GetRunSummary

(

TFile *

f

)

Gets the first entry from run summary tree.

User is responsible for deleting the new object that gets created

Definition at line 949 of file DSTUtil.cxx.

References GetTreeByName(), and tree.

Referenced by DSTUtil::Bflux::GetBeamFluxForRun(), GetRunSummaryVerbose(), TOFTempCorrector::init(), main(), ValiMC::SetupRunSummary(), and GlobalPid::SetupRunSummary().

{
  TTree* tree = GetTreeByName(f, "FillEventSummary", "fRunTree");
  if (tree && tree->GetEntries() > 0) {
    MIPPRunSummary* rs = new MIPPRunSummary;
    tree->SetBranchAddress("fRunInfo", &rs);
    tree->GetEntry(0);
    if (tree->GetEntries() > 1) {
      cout << "DSTUtil::GetRunSummary: Warning: The file " 
           << f->GetName()
           << " contains more than 1 RunSummary. The first is returned"
           << "\n (Multiple subruns processed into one DST file.)"
           << endl;
    }
    return rs;
  }
  return 0;  
}

MIPPRunSummary * DSTUtil::GetRunSummaryVerbose

(

TFile *

f

)

like DSTUtil::GetRunSummary(TFile* f), but print some info about this run to stdout.

Definition at line 972 of file DSTUtil.cxx.

References GetRunSummary(), and MIPPRunSummary::Print().

{
  MIPPRunSummary* rs = GetRunSummary(f);
  if (rs) rs->Print();
  return rs;
}

MIPPRunSummary * DSTUtil::GetRunSummary	(	TFile *	f,
		const int &	subrun
	)

Gets entry for subrun from run summary tree.

User is responsible for deleting the new object that gets created

Definition at line 983 of file DSTUtil.cxx.

References GetTreeByName(), MIPPRunSummary::Print(), MIPPRunSummary::subrun, and tree.

{
  TTree* tree = GetTreeByName(f, "FillEventSummary", "fRunTree");
  if (!tree) return 0;
  MIPPRunSummary* rs = new MIPPRunSummary;
  tree->SetBranchAddress("fRunInfo", &rs);
  int fEntry = -1;
  for (int i=0; i!=tree->GetEntries(); ++i) {
    tree->GetEntry(i);
    if (subrun == rs->subrun) {
      if (fEntry == -1) fEntry = i;
      else {
        cout << "DSTUtil::GetRunSummary: Warning: The file " 
             << f->GetName()
             << " contains more than 1 RunSummary for subrun " 
             << subrun
             << ". The first is returned.\n (This really should "
             << "not happen. You should check this DST.)" << endl;
      }
    }
  }
  if (fEntry>=0) {
    tree->GetEntry(fEntry);
    rs->Print();
    return rs;
  }
  //don't have this subrun
  delete rs;
  return 0;  
}

TTree * DSTUtil::GetMCMatchTree

(

TFile *

f

)

Given a MIPP histogram file, return MCMatch tree.

This function will return the tree from the first DST folder that it finds

Definition at line 1018 of file DSTUtil.cxx.

References GetTreeByName().

Referenced by main().

{
  return GetTreeByName(f, "MCMatch", "fMatchesNt");
}

TTree* DSTUtil::GetVeeTree	(	TFile *	f,
		string &	ntuple
	)

Definition at line 1023 of file DSTUtil.cxx.

References GetTreeByName().

Referenced by plot_MCTree(), and plot_RCTree().

{
  return GetTreeByName(f, "VFinder", ntuple);
}

TTree* DSTUtil::GetVeeTree	(	TFile *	f,
		std::string &	ntuple
	)

---