NuBeam Class Reference

#include <NuBeam.h>

List of all members.

Public Member Functions
	NuBeam ()
	~NuBeam ()
Bool_t	IsGoodSpillAndFillPot (const NtpBDLiteRecord *recBD, const NuConfig &config, NuEvent &nu) const

---

Constructor & Destructor Documentation

NuBeam::NuBeam (   )

Definition at line 35 of file NuBeam.cxx. References MSG. { MSG("NuBeam",Msg::kDebug) <<"Running NuBeam Constructor..."<<endl; MSG("NuBeam",Msg::kDebug) <<"Finished NuBeam Constructor"<<endl; }

NuBeam::~NuBeam (   )

Definition at line 47 of file NuBeam.cxx. References MSG. { MSG("NuBeam",Msg::kDebug) <<"Running NuBeam Destructor..."<<endl; MSG("NuBeam",Msg::kDebug) <<"Finished NuBeam Destructor"<<endl; }

---

Member Function Documentation

Bool_t NuBeam::IsGoodSpillAndFillPot (  const NtpBDLiteRecord *  recBD, const NuConfig &  config, NuEvent &  nu )  const

This function makes plots/counts the POTs NOTE: in order to get POT counting right when the coil is bad then the coil information has to be obtained before running this function Definition at line 59 of file NuBeam.cxx. References BeamType::AsString(), BeamMonSpill::BeamType(), BeamMonSpill::BpmAtTarget(), NtpBDLiteRecord::bpmint, NtpBDLiteRecord::bposx, NtpBDLiteRecord::bposy, NtpBDLiteRecord::bwidx, NtpBDLiteRecord::bwidy, NuEvent::coilIsOk, BMSpillAna::Config(), NuConfig::detector, RecRecordImp< T >::GetHeader(), BMSpillAna::GetSpill(), BeamMonSpill::GetStatusBits(), BeamDataLiteHeader::GetTimeDiffStreamSpill(), NuEvent::goodBeamSntp, NtpBDLiteRecord::horncur, NuConfig::hornCurrent, NuEvent::isGoodDataQuality, Registry::LockValues(), MAXMSG, NuConfig::mcVersion, MSG, NuEvent::pot, Registry::PrettyPrint(), BMSpillAna::Print(), NuUtilities::ReturnConventionsBeamType(), NuConfig::runPeriod, BMSpillAna::SelectSpill(), Registry::Set(), BMSpillAna::SetSpill(), NtpBDLiteRecord::tor101, NtpBDLiteRecord::tortgt, NtpBDLiteRecord::tr101d, NtpBDLiteRecord::trtgtd, Registry::UnLockValues(), BMSpillAna::UseDatabaseCuts(), and NuEvent::useDBForDataQuality. Referenced by NuAnalysis::ChargeSeparationOneSnarl(), NuAnalysis::ChargeSignCut(), NuAnalysis::Efficiencies(), NuAnalysis::EnergySpect(), NuAnalysis::LIRejectionTest(), NuAnalysis::MakeFullDST(), NuAnalysis::N_1(), and NuAnalysis::NuMuBarAppearance(). { static TH1F* hPottor101=0; static TH1F* hPottr101d=0; static TH1F* hPottortgt=0; static TH1F* hPottrtgtd=0; //good pot static TH1F* hPotGoodtor101=0; static TH1F* hPotGoodtr101d=0; static TH1F* hPotGoodtortgt=0; static TH1F* hPotGoodtrtgtd=0; //bad pot static TH1F* hPotBadtor101=0; static TH1F* hPotBadtr101d=0; static TH1F* hPotBadtortgt=0; static TH1F* hPotBadtrtgtd=0; //horncur, etc static TH1F* hHorncur=0; static TH1F* hHorncurBad=0; static TH1F* hBwidx=0; static TH1F* hBwidxBad=0; static TH1F* hBwidy=0; static TH1F* hBwidyBad=0; static TH1F* hBposx=0; static TH1F* hBposxBad=0; static TH1F* hBposy=0; static TH1F* hBposyBad=0; static TH1F* hBpmint=0; static TH1F* hBpmintBad=0; static TH1F* hAvBposx=0; static TH1F* hAvBposxBad=0; static TH1F* hAvBposy=0; static TH1F* hAvBposyBad=0; static TH1F* hDeltaBposx=0; static TH1F* hDeltaBposxBad=0; static TH1F* hDeltaBposy=0; static TH1F* hDeltaBposyBad=0; static TH1F* hBDTimeDiff=0; static TH1F* hBDTimeDiffBad=0; static TH1F* hBDTimeDiffBadBad=0; static TH1F* hBeamType=0; static TH1F* hBeamTypeBad=0; static TH1F* hTargetIn=0; static TH1F* hTargetInBad=0; static TH1F* hSpillsPerFile=0; if (!hPottor101){ MAXMSG("NuBeam",Msg::kDebug,1) <<"Creating POT plots..."<<endl; //these 1st four POT histograms also take into account the coil //the PotGood histos below are for all spills where the beam //was good (regardless of the coil) hPottor101=new TH1F("hPottor101","hPottor101",100000,-10,90); hPottor101->GetXaxis()->SetTitle("POT (x10^{12})"); hPottor101->GetXaxis()->CenterTitle(); hPottor101->GetYaxis()->SetTitle(""); hPottor101->GetYaxis()->CenterTitle(); hPottor101->SetFillColor(0); hPottor101->SetLineColor(1); //hPottor101->SetBit(TH1::kCanRebin); hPottr101d=new TH1F("hPottr101d","hPottr101d",100000,-10,90); hPottr101d->GetXaxis()->SetTitle("POT (x10^{12})"); hPottr101d->GetXaxis()->CenterTitle(); hPottr101d->GetYaxis()->SetTitle(""); hPottr101d->GetYaxis()->CenterTitle(); hPottr101d->SetFillColor(0); hPottr101d->SetLineColor(1); //hPottr101d->SetBit(TH1::kCanRebin); hPottortgt=new TH1F("hPottortgt","hPottortgt",100000,-10,90); hPottortgt->GetXaxis()->SetTitle("POT (x10^{12})"); hPottortgt->GetXaxis()->CenterTitle(); hPottortgt->GetYaxis()->SetTitle(""); hPottortgt->GetYaxis()->CenterTitle(); hPottortgt->SetFillColor(0); hPottortgt->SetLineColor(1); //hPottortgt->SetBit(TH1::kCanRebin); hPottrtgtd=new TH1F("hPottrtgtd","hPottrtgtd",100000,-10,90); hPottrtgtd->GetXaxis()->SetTitle("POT (x10^{12})"); hPottrtgtd->GetXaxis()->CenterTitle(); hPottrtgtd->GetYaxis()->SetTitle(""); hPottrtgtd->GetYaxis()->CenterTitle(); hPottrtgtd->SetFillColor(0); hPottrtgtd->SetLineColor(1); //hPottrtgtd->SetBit(TH1::kCanRebin); //Good POT hPotGoodtor101=new TH1F("hPotGoodtor101","hPotGoodtor101", 100000,-10,90); hPotGoodtor101->GetXaxis()->SetTitle("POT (x10^{12})"); hPotGoodtor101->GetXaxis()->CenterTitle(); hPotGoodtor101->GetYaxis()->SetTitle(""); hPotGoodtor101->GetYaxis()->CenterTitle(); hPotGoodtor101->SetFillColor(0); hPotGoodtor101->SetLineColor(1); //hPotGoodtor101->SetBit(TH1::kCanRebin); hPotGoodtr101d=new TH1F("hPotGoodtr101d","hPotGoodtr101d", 100000,-10,90); hPotGoodtr101d->GetXaxis()->SetTitle("POT (x10^{12})"); hPotGoodtr101d->GetXaxis()->CenterTitle(); hPotGoodtr101d->GetYaxis()->SetTitle(""); hPotGoodtr101d->GetYaxis()->CenterTitle(); hPotGoodtr101d->SetFillColor(0); hPotGoodtr101d->SetLineColor(1); //hPotGoodtr101d->SetBit(TH1::kCanRebin); hPotGoodtortgt=new TH1F("hPotGoodtortgt","hPotGoodtortgt", 100000,-10,90); hPotGoodtortgt->GetXaxis()->SetTitle("POT (x10^{12})"); hPotGoodtortgt->GetXaxis()->CenterTitle(); hPotGoodtortgt->GetYaxis()->SetTitle(""); hPotGoodtortgt->GetYaxis()->CenterTitle(); hPotGoodtortgt->SetFillColor(0); hPotGoodtortgt->SetLineColor(1); //hPotGoodtortgt->SetBit(TH1::kCanRebin); hPotGoodtrtgtd=new TH1F("hPotGoodtrtgtd","hPotGoodtrtgtd", 100000,-10,90); hPotGoodtrtgtd->GetXaxis()->SetTitle("POT (x10^{12})"); hPotGoodtrtgtd->GetXaxis()->CenterTitle(); hPotGoodtrtgtd->GetYaxis()->SetTitle(""); hPotGoodtrtgtd->GetYaxis()->CenterTitle(); hPotGoodtrtgtd->SetFillColor(0); hPotGoodtrtgtd->SetLineColor(1); //hPotGoodtrtgtd->SetBit(TH1::kCanRebin); //bad POT hPotBadtor101=new TH1F("hPotBadtor101","hPotBadtor101", 100000,-10,90); hPotBadtor101->GetXaxis()->SetTitle("POT (x10^{12})"); hPotBadtor101->GetXaxis()->CenterTitle(); hPotBadtor101->GetYaxis()->SetTitle(""); hPotBadtor101->GetYaxis()->CenterTitle(); hPotBadtor101->SetFillColor(0); hPotBadtor101->SetLineColor(1); //hPotBadtor101->SetBit(TH1::kCanRebin); hPotBadtr101d=new TH1F("hPotBadtr101d","hPotBadtr101d", 100000,-10,90); hPotBadtr101d->GetXaxis()->SetTitle("POT (x10^{12})"); hPotBadtr101d->GetXaxis()->CenterTitle(); hPotBadtr101d->GetYaxis()->SetTitle(""); hPotBadtr101d->GetYaxis()->CenterTitle(); hPotBadtr101d->SetFillColor(0); hPotBadtr101d->SetLineColor(1); //hPotBadtr101d->SetBit(TH1::kCanRebin); hPotBadtortgt=new TH1F("hPotBadtortgt","hPotBadtortgt", 100000,-10,90); hPotBadtortgt->GetXaxis()->SetTitle("POT (x10^{12})"); hPotBadtortgt->GetXaxis()->CenterTitle(); hPotBadtortgt->GetYaxis()->SetTitle(""); hPotBadtortgt->GetYaxis()->CenterTitle(); hPotBadtortgt->SetFillColor(0); hPotBadtortgt->SetLineColor(1); //hPotBadtortgt->SetBit(TH1::kCanRebin); hPotBadtrtgtd=new TH1F("hPotBadtrtgtd","hPotBadtrtgtd", 100000,-10,90); hPotBadtrtgtd->GetXaxis()->SetTitle("POT (x10^{12})"); hPotBadtrtgtd->GetXaxis()->CenterTitle(); hPotBadtrtgtd->GetYaxis()->SetTitle(""); hPotBadtrtgtd->GetYaxis()->CenterTitle(); hPotBadtrtgtd->SetFillColor(0); hPotBadtrtgtd->SetLineColor(1); //hPotBadtrtgtd->SetBit(TH1::kCanRebin); //horn current hHorncur=new TH1F("hHorncur","hHorncur", 1000,-250,250); hHorncur->GetXaxis()->SetTitle("Horn Current"); hHorncur->GetXaxis()->CenterTitle(); hHorncur->GetYaxis()->SetTitle(""); hHorncur->GetYaxis()->CenterTitle(); hHorncur->SetFillColor(0); hHorncur->SetLineColor(1); //hHorncur->SetBit(TH1::kCanRebin); hHorncurBad=new TH1F("hHorncurBad","hHorncurBad", 1000,-250,250); hHorncurBad->GetXaxis()->SetTitle("Horn Current"); hHorncurBad->GetXaxis()->CenterTitle(); hHorncurBad->GetYaxis()->SetTitle(""); hHorncurBad->GetYaxis()->CenterTitle(); hHorncurBad->SetFillColor(0); hHorncurBad->SetLineColor(1); //hHorncurBad->SetBit(TH1::kCanRebin); hBwidx=new TH1F("hBwidx","hBwidx", 1000,-0.002,0.01); hBwidx->GetXaxis()->SetTitle("Beam Width x"); hBwidx->GetXaxis()->CenterTitle(); hBwidx->GetYaxis()->SetTitle(""); hBwidx->GetYaxis()->CenterTitle(); hBwidx->SetFillColor(0); hBwidx->SetLineColor(1); //hBwidx->SetBit(TH1::kCanRebin); hBwidxBad=new TH1F("hBwidxBad","hBwidxBad", 1000,-0.002,0.01); hBwidxBad->GetXaxis()->SetTitle("Beam Width x"); hBwidxBad->GetXaxis()->CenterTitle(); hBwidxBad->GetYaxis()->SetTitle(""); hBwidxBad->GetYaxis()->CenterTitle(); hBwidxBad->SetFillColor(0); hBwidxBad->SetLineColor(1); //hBwidxBad->SetBit(TH1::kCanRebin); hBwidy=new TH1F("hBwidy","hBwidy", 1000,-0.002,0.01); hBwidy->GetXaxis()->SetTitle("Beam Width y"); hBwidy->GetXaxis()->CenterTitle(); hBwidy->GetYaxis()->SetTitle(""); hBwidy->GetYaxis()->CenterTitle(); hBwidy->SetFillColor(0); hBwidy->SetLineColor(1); //hBwidy->SetBit(TH1::kCanRebin); hBwidyBad=new TH1F("hBwidyBad","hBwidyBad", 1000,-0.002,0.01); hBwidyBad->GetXaxis()->SetTitle("Beam Width y"); hBwidyBad->GetXaxis()->CenterTitle(); hBwidyBad->GetYaxis()->SetTitle(""); hBwidyBad->GetYaxis()->CenterTitle(); hBwidyBad->SetFillColor(0); hBwidyBad->SetLineColor(1); //hBwidyBad->SetBit(TH1::kCanRebin); //beam intensity per batch hBpmint=new TH1F("hBpmint","hBpmint", 2010,-10,2000); hBpmint->GetXaxis()->SetTitle("Beam intensity"); hBpmint->GetXaxis()->CenterTitle(); hBpmint->GetYaxis()->SetTitle(""); hBpmint->GetYaxis()->CenterTitle(); hBpmint->SetFillColor(0); hBpmint->SetLineColor(1); //hBpmint->SetBit(TH1::kCanRebin); hBpmintBad=new TH1F("hBpmintBad","hBpmintBad", 2010,-10,2000); hBpmintBad->GetXaxis()->SetTitle("Beam intensity"); hBpmintBad->GetXaxis()->CenterTitle(); hBpmintBad->GetYaxis()->SetTitle(""); hBpmintBad->GetYaxis()->CenterTitle(); hBpmintBad->SetFillColor(0); hBpmintBad->SetLineColor(1); //hBpmintBad->SetBit(TH1::kCanRebin); hBposx=new TH1F("hBposx","hBposx", 5000,-0.01,0.01); hBposx->GetXaxis()->SetTitle("Beam Position x"); hBposx->GetXaxis()->CenterTitle(); hBposx->GetYaxis()->SetTitle(""); hBposx->GetYaxis()->CenterTitle(); hBposx->SetFillColor(0); hBposx->SetLineColor(1); //hBposx->SetBit(TH1::kCanRebin); hBposxBad=new TH1F("hBposxBad","hBposxBad", 5000,-0.01,0.01); hBposxBad->GetXaxis()->SetTitle("Beam Position x"); hBposxBad->GetXaxis()->CenterTitle(); hBposxBad->GetYaxis()->SetTitle(""); hBposxBad->GetYaxis()->CenterTitle(); hBposxBad->SetFillColor(0); hBposxBad->SetLineColor(1); //hBposxBad->SetBit(TH1::kCanRebin); hBposy=new TH1F("hBposy","hBposy", 5000,-0.01,0.01);// hBposy->GetXaxis()->SetTitle("Beam Position y"); hBposy->GetXaxis()->CenterTitle(); hBposy->GetYaxis()->SetTitle(""); hBposy->GetYaxis()->CenterTitle(); hBposy->SetFillColor(0); hBposy->SetLineColor(1); //hBposy->SetBit(TH1::kCanRebin); hBposyBad=new TH1F("hBposyBad","hBposyBad", 5000,-0.01,0.01); hBposyBad->GetXaxis()->SetTitle("Beam Position y"); hBposyBad->GetXaxis()->CenterTitle(); hBposyBad->GetYaxis()->SetTitle(""); hBposyBad->GetYaxis()->CenterTitle(); hBposyBad->SetFillColor(0); hBposyBad->SetLineColor(1); //hBposyBad->SetBit(TH1::kCanRebin); hAvBposx=new TH1F("hAvBposx","hAvBposx", 5000,-0.01,0.01); hAvBposx->GetXaxis()->SetTitle("Av. Beam Position x"); hAvBposx->GetXaxis()->CenterTitle(); hAvBposx->GetYaxis()->SetTitle(""); hAvBposx->GetYaxis()->CenterTitle(); hAvBposx->SetFillColor(0); hAvBposx->SetLineColor(1); //hAvBposx->SetBit(TH1::kCanRebin); hAvBposxBad=new TH1F("hAvBposxBad","hAvBposxBad", 5000,-0.01,0.01); hAvBposxBad->GetXaxis()->SetTitle("Av. Beam Position x"); hAvBposxBad->GetXaxis()->CenterTitle(); hAvBposxBad->GetYaxis()->SetTitle(""); hAvBposxBad->GetYaxis()->CenterTitle(); hAvBposxBad->SetFillColor(0); hAvBposxBad->SetLineColor(1); //hAvBposxBad->SetBit(TH1::kCanRebin); hAvBposy=new TH1F("hAvBposy","hAvBposy", 5000,-0.01,0.01);// hAvBposy->GetXaxis()->SetTitle("Av. Beam Position y"); hAvBposy->GetXaxis()->CenterTitle(); hAvBposy->GetYaxis()->SetTitle(""); hAvBposy->GetYaxis()->CenterTitle(); hAvBposy->SetFillColor(0); hAvBposy->SetLineColor(1); //hAvBposy->SetBit(TH1::kCanRebin); hAvBposyBad=new TH1F("hAvBposyBad","hAvBposyBad", 5000,-0.01,0.01); hAvBposyBad->GetXaxis()->SetTitle("Av. Beam Position y"); hAvBposyBad->GetXaxis()->CenterTitle(); hAvBposyBad->GetYaxis()->SetTitle(""); hAvBposyBad->GetYaxis()->CenterTitle(); hAvBposyBad->SetFillColor(0); hAvBposyBad->SetLineColor(1); //hAvBposyBad->SetBit(TH1::kCanRebin); hDeltaBposx=new TH1F("hDeltaBposx","hDeltaBposx", 5000,-0.01,0.01);// hDeltaBposx->GetXaxis()->SetTitle("Delta Beam Position x"); hDeltaBposx->GetXaxis()->CenterTitle(); hDeltaBposx->GetYaxis()->SetTitle(""); hDeltaBposx->GetYaxis()->CenterTitle(); hDeltaBposx->SetFillColor(0); hDeltaBposx->SetLineColor(1); //hDeltaBposx->SetBit(TH1::kCanRebin); hDeltaBposxBad=new TH1F("hDeltaBposxBad","hDeltaBposxBad", 5000,-0.01,0.01); hDeltaBposxBad->GetXaxis()->SetTitle("Delta Beam Position x"); hDeltaBposxBad->GetXaxis()->CenterTitle(); hDeltaBposxBad->GetYaxis()->SetTitle(""); hDeltaBposxBad->GetYaxis()->CenterTitle(); hDeltaBposxBad->SetFillColor(0); hDeltaBposxBad->SetLineColor(1); //hDeltaBposxBad->SetBit(TH1::kCanRebin); hDeltaBposy=new TH1F("hDeltaBposy","hDeltaBposy", 5000,-0.01,0.01);// hDeltaBposy->GetXaxis()->SetTitle("Delta Beam Position y"); hDeltaBposy->GetXaxis()->CenterTitle(); hDeltaBposy->GetYaxis()->SetTitle(""); hDeltaBposy->GetYaxis()->CenterTitle(); hDeltaBposy->SetFillColor(0); hDeltaBposy->SetLineColor(1); //hDeltaBposy->SetBit(TH1::kCanRebin); hDeltaBposyBad=new TH1F("hDeltaBposyBad","hDeltaBposyBad", 5000,-0.01,0.01); hDeltaBposyBad->GetXaxis()->SetTitle("Delta Beam Position y"); hDeltaBposyBad->GetXaxis()->CenterTitle(); hDeltaBposyBad->GetYaxis()->SetTitle(""); hDeltaBposyBad->GetYaxis()->CenterTitle(); hDeltaBposyBad->SetFillColor(0); hDeltaBposyBad->SetLineColor(1); //hDeltaBposyBad->SetBit(TH1::kCanRebin); //timediff hBDTimeDiff=new TH1F("hBDTimeDiff","hBDTimeDiff", 5000,-2.5,2.5); hBDTimeDiff->GetXaxis()->SetTitle("Beam/ND Time Difference (s)"); hBDTimeDiff->GetXaxis()->CenterTitle(); hBDTimeDiff->GetYaxis()->SetTitle(""); hBDTimeDiff->GetYaxis()->CenterTitle(); hBDTimeDiff->SetFillColor(0); hBDTimeDiff->SetLineColor(1); //hBDTimeDiff->SetBit(TH1::kCanRebin); hBDTimeDiffBad=new TH1F("hBDTimeDiffBad","hBDTimeDiffBad", 5000,-2.5,2.5); hBDTimeDiffBad->GetXaxis()->SetTitle("Beam/ND Time Difference (s)"); hBDTimeDiffBad->GetXaxis()->CenterTitle(); hBDTimeDiffBad->GetYaxis()->SetTitle(""); hBDTimeDiffBad->GetYaxis()->CenterTitle(); hBDTimeDiffBad->SetFillColor(0); hBDTimeDiffBad->SetLineColor(1); //hBDTimeDiffBad->SetBit(TH1::kCanRebin); hBDTimeDiffBadBad=new TH1F("hBDTimeDiffBadBad","hBDTimeDiffBadBad", 5000,-2.5,1000); hBDTimeDiffBadBad->GetXaxis()-> SetTitle("Beam/ND Time Difference (s)"); hBDTimeDiffBadBad->GetXaxis()->CenterTitle(); hBDTimeDiffBadBad->GetYaxis()->SetTitle(""); hBDTimeDiffBadBad->GetYaxis()->CenterTitle(); hBDTimeDiffBadBad->SetFillColor(0); hBDTimeDiffBadBad->SetLineColor(1); //hBDTimeDiffBadBad->SetBit(TH1::kCanRebin); hBeamType=new TH1F("hBeamType","hBeamType", 100,-10,90);// hBeamType->GetXaxis()->SetTitle("Beam Type"); hBeamType->GetXaxis()->CenterTitle(); hBeamType->GetYaxis()->SetTitle(""); hBeamType->GetYaxis()->CenterTitle(); hBeamType->SetFillColor(0); hBeamType->SetLineColor(1); //hBeamType->SetBit(TH1::kCanRebin); hBeamTypeBad=new TH1F("hBeamTypeBad","hBeamTypeBad", 100,-10,90); hBeamTypeBad->GetXaxis()->SetTitle("Beam Type"); hBeamTypeBad->GetXaxis()->CenterTitle(); hBeamTypeBad->GetYaxis()->SetTitle(""); hBeamTypeBad->GetYaxis()->CenterTitle(); hBeamTypeBad->SetFillColor(0); hBeamTypeBad->SetLineColor(1); //hBeamTypeBad->SetBit(TH1::kCanRebin); hTargetIn=new TH1F("hTargetIn","hTargetIn", 100,-10,90);// hTargetIn->GetXaxis()->SetTitle("Target In"); hTargetIn->GetXaxis()->CenterTitle(); hTargetIn->GetYaxis()->SetTitle(""); hTargetIn->GetYaxis()->CenterTitle(); hTargetIn->SetFillColor(0); hTargetIn->SetLineColor(1); //hTargetIn->SetBit(TH1::kCanRebin); hTargetInBad=new TH1F("hTargetInBad","hTargetInBad", 100,-10,90); hTargetInBad->GetXaxis()->SetTitle("Target In"); hTargetInBad->GetXaxis()->CenterTitle(); hTargetInBad->GetYaxis()->SetTitle(""); hTargetInBad->GetYaxis()->CenterTitle(); hTargetInBad->SetFillColor(0); hTargetInBad->SetLineColor(1); //hTargetInBad->SetBit(TH1::kCanRebin); hSpillsPerFile=new TH1F("hSpillsPerFile","hSpillsPerFile", 10000,0,10000); hSpillsPerFile->GetXaxis()->SetTitle("hSpillsPerFile"); hSpillsPerFile->GetXaxis()->CenterTitle(); hSpillsPerFile->GetYaxis()->SetTitle(""); hSpillsPerFile->GetYaxis()->CenterTitle(); hSpillsPerFile->SetFillColor(0); hSpillsPerFile->SetLineColor(1); //hSpillsPerFile->SetBit(TH1::kCanRebin); } if (!recBD) { //this is the MC case //this is the pot for carrot MC at LE-10 Float_t potPerSpill=24.2;//x1e12 but ntp assumes e12 Float_t spillsPerFile=400; //This function, for e.g., returns the value 24.2, since PoT/snarl //for Daikon L010185 Near MC is 2.42*1013. //Convert to canonical BeamType in the case of intensity //weighted MC NuUtilities utils; BeamType::BeamType_t convBeamType = utils.ReturnConventionsBeamType(config); MAXMSG("NuBeam",Msg::kInfo,1) << "Beam type to be used for PoT counting: " << BeamType::AsString(convBeamType) << endl; potPerSpill=MCInfo::GetMCPoT (static_cast<Detector::Detector_t>(config.detector), static_cast<BeamType::BeamType_t>(convBeamType), static_cast<ReleaseType::Release_t>(config.mcVersion)); spillsPerFile=MCInfo::GetNoSnarlPerFile (static_cast<Detector::Detector_t>(config.detector), static_cast<BeamType::BeamType_t>(convBeamType), static_cast<ReleaseType::Release_t>(config.mcVersion), config.runPeriod); //remove the FD factor of 1e8 so that it can be stored in the histo //also just put a 1 for spills per file, since it's pot/file really if (config.detector==Detector::kFar) { if (potPerSpill < 0) { MAXMSG("NuBeam",Msg::kWarning,50) << "*** Warning: potPerSpill undefined for beam type: " << BeamType::AsString(convBeamType) << " (" << convBeamType << ") in the far detector." << endl << "ASSUMING 6.5e20 per file, which is standard for Daikon. " << endl << "This is dangerous - make sure the above beam type is added to DataUtil/MCInfo::GetMCPoT() ASAP." << endl; potPerSpill = 6.5e8; } potPerSpill/=1e8; if (spillsPerFile<=0) { MAXMSG("NuBeam",Msg::kInfo,1) <<"Overriding spills per file (="<<spillsPerFile <<"), now set to 1"<<endl; spillsPerFile=1; } } //sanity check: histogram overflow Int_t n=hSpillsPerFile->GetNbinsX(); if (hSpillsPerFile->GetBinContent(n+1)>0){ MSG("NuBeam",Msg::kError) <<"n bins="<<n <<", overflow="<<hSpillsPerFile->GetBinContent(n+1)<<endl; } MAXMSG("NuBeam",Msg::kInfo,1) <<"No NtpBDLiteRecord, so using potPerSpill="<<potPerSpill <<"x10^12 ND (10^20 FD)"<<" and spillsPerFile="<<spillsPerFile <<" for beamType="<<convBeamType <<", hornCurrent="<<config.hornCurrent <<", mcVersion="<<config.mcVersion<<endl; //fill histos hPottor101->Fill(potPerSpill); hPottr101d->Fill(potPerSpill); hPottortgt->Fill(potPerSpill); hPottrtgtd->Fill(potPerSpill); hSpillsPerFile->Fill(spillsPerFile); //write out the potPerSpill nu.pot=potPerSpill; nu.goodBeamSntp=true;//always true for MC //don't fill the bad histos ever in this case } //this is the data case else { MAXMSG("NuBeam",Msg::kInfo,1) <<"Using NtpBDLiteRecord for POT"<<endl; const NtpBDLiteRecord& ntpBD=*recBD; static BeamMonSpill spill; static BMSpillAna bmsa; //tell it to use the database defined values of the cuts bmsa.UseDatabaseCuts(); Bool_t changeBeamCuts=false; //change default values (only once though) static Registry* pr=0; if (!pr && changeBeamCuts){ MAXMSG("NuBeam",Msg::kInfo,1) <<"Changing default values of BMSpillAna"<<endl; pr=new Registry(); Registry& r=*pr;//get a reference r.UnLockValues(); //r.Set("TargetIn",1); //r.Set("BeamType",-1);//don't care if ME,LE,etc //r.Set("PosTgtXMin",-2.0*Munits::mm); //have to make sure that zero is not included r.Set("PosTgtXMax",-0.0000001*Munits::mm); r.Set("PosTgtYMin", 0.0000001*Munits::mm); //r.Set("PosTgtYMax", 2.0*Munits::mm); r.LockValues(); bmsa.Config(r); //print the registry r.PrettyPrint(cout); } //give it the spill bmsa.SetSpill(ntpBD,spill); Double_t xmean=0; Double_t ymean=0; Double_t xrms=0; Double_t yrms=0; bmsa.GetSpill().BpmAtTarget(xmean,ymean,xrms,yrms); Int_t beamType=spill.GetStatusBits().beam_type; Int_t targetIn=spill.GetStatusBits().target_in; //write to the object whether the beam is good nu.goodBeamSntp=bmsa.SelectSpill(); //set the pot value nu.pot=ntpBD.tortgt; if (nu.pot<0) nu.pot=0;//bad spills can be <0 but don't want to sum Float_t tor101=ntpBD.tor101; if (tor101<0) tor101=0; Float_t tr101d=ntpBD.tr101d; if (tr101d<0) tr101d=0; Float_t tortgt=ntpBD.tortgt; if (tortgt<0) tortgt=0; Float_t trtgtd=ntpBD.trtgtd; if (trtgtd<0) trtgtd=0; //check if a good spill if (!nu.goodBeamSntp) { MAXMSG("NuBeam",Msg::kWarning, 20) <<"NuBeam: BAD SPILL:" << endl << " > Maximum time diffrence (s): " << ntpBD.GetHeader().GetTimeDiffStreamSpill() << endl << " > Spill intensity (1e12 pot): " << trtgtd << endl << " > Horn Current (kA): " << ntpBD.horncur << endl << " > fBeamType: " << BeamType::AsString(static_cast<BeamType::BeamType_t>(spill.BeamType())) <<" ("<<spill.BeamType()<<")"<<endl << " > Horizontal beam position (mm): " << xmean << endl << " > Vertical beam position (mm): " << ymean << endl << " > Horizontal beam width a (mm): " << ntpBD.bwidx << endl << " > Vertical beam width a (mm): " << ntpBD.bwidy << endl; //for the Bad plots to have any meaning you have to make sure //that the spill was in time, otherwise you just plot //data from spills ages ago. if (ntpBD.GetHeader().GetTimeDiffStreamSpill()>1){ hBDTimeDiffBadBad->Fill(ntpBD.GetHeader(). GetTimeDiffStreamSpill()); } else { hPotBadtor101->Fill(tor101); hPotBadtr101d->Fill(tr101d); hPotBadtortgt->Fill(tortgt); hPotBadtrtgtd->Fill(trtgtd); hHorncurBad->Fill(ntpBD.horncur); hBwidxBad->Fill(ntpBD.bwidx); hBwidyBad->Fill(ntpBD.bwidy); hBDTimeDiffBad->Fill(ntpBD.GetHeader(). GetTimeDiffStreamSpill()); hAvBposxBad->Fill(xmean); hAvBposyBad->Fill(ymean); for (Int_t batch=0;batch<6;batch++){ hBpmintBad->Fill(ntpBD.bpmint[batch]); //check if batch existed in spill if (ntpBD.bpmint[batch]>1) { hBposxBad->Fill(ntpBD.bposx[batch]); hBposyBad->Fill(ntpBD.bposy[batch]); //calc delta from mean position //should probably use some sort of intensity weighting hDeltaBposxBad->Fill(ntpBD.bposx[batch]-xmean); hDeltaBposyBad->Fill(ntpBD.bposy[batch]-ymean); } } hTargetInBad->Fill(targetIn); hBeamTypeBad->Fill(beamType); } return nu.goodBeamSntp; } MAXMSG("NuBeam",Msg::kInfo,1) <<"Found first good spill"<<", trtgtd="<<trtgtd <<", timeDiff="<<ntpBD.GetHeader().GetTimeDiffStreamSpill() <<", beamType="<<beamType<<", tgtIn="<<targetIn<<endl; static Bool_t firstTime=true; if (firstTime) { bmsa.Print(); cout<<endl; firstTime=false; } //perform sanity check that hornCurrent is //within +/- 4% of the config one if (!(fabs(ntpBD.horncur)>fabs(config.hornCurrent)*0.96 && fabs(ntpBD.horncur)<fabs(config.hornCurrent)*1.04)){ MSG("NuBeam",Msg::kError) <<"HornCurrent problem for BeamType (sntp)="<<beamType <<", event hornCur="<<ntpBD.horncur <<" config.hornCurrent="<<config.hornCurrent <<endl; TH1F* hError=dynamic_cast<TH1F*>(gROOT->FindObject("hError")); if (hError) hError->Fill(1);//set an error flag } //these histograms are the ones used for POT counting //have to check the coil and data quality as well if (nu.goodBeamSntp && nu.coilIsOk && (nu.isGoodDataQuality || !nu.useDBForDataQuality)) { hPottor101->Fill(tor101); hPottr101d->Fill(tr101d); hPottortgt->Fill(tortgt); hPottrtgtd->Fill(trtgtd); } //these are just for when the beam is good (regardless of coil) hPotGoodtor101->Fill(tor101); hPotGoodtr101d->Fill(tr101d); hPotGoodtortgt->Fill(tortgt); hPotGoodtrtgtd->Fill(trtgtd); hHorncur->Fill(ntpBD.horncur); hBwidx->Fill(ntpBD.bwidx); hBwidy->Fill(ntpBD.bwidy); hBDTimeDiff->Fill(ntpBD.GetHeader().GetTimeDiffStreamSpill()); hAvBposx->Fill(xmean); hAvBposy->Fill(ymean); for (Int_t batch=0;batch<6;batch++){ hBpmint->Fill(ntpBD.bpmint[batch]); //check if batch existed in spill if (ntpBD.bpmint[batch]>1) { hBposx->Fill(ntpBD.bposx[batch]); hBposy->Fill(ntpBD.bposy[batch]); //calc delta from mean position //should probably use some sort of intensity weighting hDeltaBposx->Fill(ntpBD.bposx[batch]-xmean); hDeltaBposy->Fill(ntpBD.bposy[batch]-ymean); } } hTargetIn->Fill(targetIn); hBeamType->Fill(beamType); } return nu.goodBeamSntp; }

---

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

• NuBeam.h
• NuBeam.cxx

---