All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Private Attributes | List of all members
sbn::NuMIRetriever Class Reference
Inheritance diagram for sbn::NuMIRetriever:

Public Member Functions

 NuMIRetriever (fhicl::ParameterSet const &p)
 
 NuMIRetriever (NuMIRetriever const &)=delete
 
 NuMIRetriever (NuMIRetriever &&)=delete
 
NuMIRetriever & operator= (NuMIRetriever const &)=delete
 
NuMIRetriever & operator= (NuMIRetriever &&)=delete
 
void produce (art::Event &e) override
 
void beginSubRun (art::SubRun &sr) override
 
void endSubRun (art::SubRun &sr) override
 

Private Attributes

std::vector< sbn::NuMISpillInfo > fOutbeamInfos
 
double fTimePad
 
double fBFPEpsilion
 
std::string fURL
 
std::string raw_data_label_
 
std::string fDeviceUsedForTiming
 
int TotalBeamSpills
 
art::ServiceHandle
< ifbeam_ns::IFBeam > 		ifbeam_handle
 
std::unique_ptr
< ifbeam_ns::BeamFolder > 		bfp
 

Detailed Description

Definition at line 40 of file NuMIRetriever_module.cc.

Constructor & Destructor Documentation

sbn::NuMIRetriever::NuMIRetriever ( fhicl::ParameterSet const &  p)
explicit

Definition at line 67 of file NuMIRetriever_module.cc.

68  : EDProducer{p},
69  fTimePad(p.get<double>("TimePadding", 0.5)), //epsilon in seconds, buffer time to look for spills
70  fBFPEpsilion(p.get<double>("BFPEpsilon", 0.02)), // 20 ms, tuned for BNB, check for NuMI here might need to be larger
71  raw_data_label_(p.get<std::string>("raw_data_label")),
72  fDeviceUsedForTiming(p.get<std::string>("DeviceUsedForTiming")),
73  bfp(ifbeam_handle->getBeamFolder(p.get<std::string>("Bundle"), p.get<std::string>("URL"), p.get<double>("TimeWindow")))
74 {
75 
76  bfp->set_epsilon(fBFPEpsilion);
77  bfp->setValidWindow(500.);
78  produces<std::vector<sbn::NuMISpillInfo>, art::InSubRun>();
79  TotalBeamSpills = 0;
80 }
p
pdgs p
Definition: selectors.fcl:22
sbn::NuMIRetriever::bfp
std::unique_ptr< ifbeam_ns::BeamFolder > bfp
Definition: NuMIRetriever_module.cc:64
sbn::NuMIRetriever::fDeviceUsedForTiming
std::string fDeviceUsedForTiming
Definition: NuMIRetriever_module.cc:61
sbn::NuMIRetriever::ifbeam_handle
art::ServiceHandle< ifbeam_ns::IFBeam > ifbeam_handle
Definition: NuMIRetriever_module.cc:63
sbn::NuMIRetriever::NuMIRetriever ( NuMIRetriever const &  )
delete
sbn::NuMIRetriever::NuMIRetriever ( NuMIRetriever &&  )
delete

Member Function Documentation

void sbn::NuMIRetriever::beginSubRun ( art::SubRun &  sr)
override

Definition at line 236 of file NuMIRetriever_module.cc.

237 {
238  return;
239 }
void sbn::NuMIRetriever::endSubRun ( art::SubRun &  sr)
override

Definition at line 241 of file NuMIRetriever_module.cc.

242 {
243  // We will add all of the BNBSpillInfo data-products to the
244  // art::SubRun so it persists
245  // currently this is ~2.7 kB/event or ~0.07 kB/spill
246 
247  std::cout << "Total number of DAQ Spills : " << TotalBeamSpills << std::endl;
248  std::cout << "Total number of Selected Spills : " << fOutbeamInfos.size() << std::endl;
249 
250  auto p = std::make_unique< std::vector< sbn::NuMISpillInfo > >(fOutbeamInfos);
251 
252  sr.put(std::move(p), art::subRunFragment());
253 
254  // Clear the (now old) infos
255  fOutbeamInfos.clear();
256 
257  return;
258 }
p
pdgs p
Definition: selectors.fcl:22
sbn::NuMIRetriever::fOutbeamInfos
std::vector< sbn::NuMISpillInfo > fOutbeamInfos
Definition: NuMIRetriever_module.cc:55
cout
BEGIN_PROLOG could also be cout
Definition: messageservice.fcl:10
NuMIRetriever& sbn::NuMIRetriever::operator= ( NuMIRetriever const &  )
delete
NuMIRetriever& sbn::NuMIRetriever::operator= ( NuMIRetriever &&  )
delete
void sbn::NuMIRetriever::produce ( art::Event &  e)
override

Definition at line 82 of file NuMIRetriever_module.cc.

83 {
84 
85  // If this is the first event in the run, then ignore it
86  // We do not currently have the ability to figure out the first
87  // spill that the DAQ was sensitive to, so don't try to save any
88  // spill information
89  //
90  // TODO: long-term goal -- can we fix this?
91  if (e.event() == 1) return;
92 
93  int gate_type = 0;
94  art::Handle< std::vector<artdaq::Fragment> > raw_data_ptr;
95  e.getByLabel(raw_data_label_, "ICARUSTriggerV2", raw_data_ptr);
96  auto const & raw_data = (*raw_data_ptr);
97 
98  double t_current_event = 0;
99  double t_previous_event = 0;
100  double number_of_gates_since_previous_event = 0;
101 
102  for(auto raw_datum : raw_data){
103 
104  uint64_t artdaq_ts = raw_datum.timestamp();
105  icarus::ICARUSTriggerV2Fragment frag(raw_datum);
106  std::string data = frag.GetDataString();
107  char *buffer = const_cast<char*>(data.c_str());
108  icarus::ICARUSTriggerInfo datastream_info = icarus::parse_ICARUSTriggerV2String(buffer);
109  gate_type = datastream_info.gate_type;
110  number_of_gates_since_previous_event = frag.getDeltaGatesNuMI();
111 
112  t_current_event = static_cast<double>(artdaq_ts)/(1000000000.); //check this offset...
113  if(gate_type == 2)
114  t_previous_event = (static_cast<double>(frag.getLastTimestampNuMI()))/(1000000000.);
115  else
116  t_previous_event = (static_cast<double>(frag.getLastTimestampOther()))/(1000000000.);
117 
118  }
119 
120  std::cout << std::setprecision(19) << "Previous : " << t_previous_event << ", Current : " << t_current_event << std::endl;
121  //We only want to process NuMI gates, i.e. type 2
122  if(gate_type == 2)
123  {
124  // Keep track of the number of beam gates the DAQ thinks
125  // are in this job
126  TotalBeamSpills += number_of_gates_since_previous_event;
127 
128  // These lines get everything primed within the IFBeamDB
129  // They seem redundant but they are needed
130  try{auto cur_vec_temp = bfp->GetNamedVector((t_previous_event)-fTimePad,"E:HP121[]");} catch (WebAPIException &we) {}
131 
132  try{auto cur_vec_temp_2 = bfp->GetNamedVector((t_current_event)+fTimePad,"E:VP121[]");} catch (WebAPIException &we) {}
133  try{auto packed_MTGTDS_temp = bfp->GetNamedVector((t_current_event)+fTimePad, "E:MTGTDS[]");} catch(WebAPIException &we) {}
134  std::vector<double> times_temps = bfp->GetTimeList(fDeviceUsedForTiming);
135 
136  int spill_count = 0;
137  // Iterating through each of the beamline times
138  for (size_t i = 0; i < times_temps.size(); i++) {
139 
140  // Only continue if these times are matched to our DAQ time
141  // plus or minus some time padding, currently using 3.3 ms
142  // which is half the Booster Rep Rate
143  if(e.event() != 1){//We already addressed the "first event" above
144  if(times_temps[i] > t_current_event){continue;}
145  if(times_temps[i] <= t_previous_event){continue;}
146  }
147 
148  //count found spills
149  spill_count++;
150 
151  //initialize all devices found in NuMISpillInfo.h in sbnobj
152  double HRNDIR = -1.;
153  double NSLINA = -1.;
154  double NSLINB = -1.;
155  double NSLINC = -1.;
156  double NSLIND = -1.;
157  double TOR101 = -1.;
158  double TORTGT = -1.;
159  double TR101D = -1.;
160  double TRTGTD = -1.;
161  std::vector< double > HP121;
162  std::vector< double > VP121;
163  std::vector< double > HPTGT;
164  std::vector< double > VPTGT;
165  std::vector< double > HITGT;
166  std::vector< double > VITGT;
167  std::vector< double > MTGTDS;
168  double TRTGTD_time = -1.;
169  std::cout << "Grabbing IFBeam info!" << std::endl;
170  try{bfp->GetNamedData(times_temps[i], "E:TRTGTD@",&TRTGTD,&TRTGTD_time);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
171  try{bfp->GetNamedData(times_temps[i], "E:TR101D",&TR101D);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
172  try{bfp->GetNamedData(times_temps[i], "E:HRNDIR",&HRNDIR);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
173  try{bfp->GetNamedData(times_temps[i], "E:NSLINA",&NSLINA);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
174  try{bfp->GetNamedData(times_temps[i], "E:NSLINB",&NSLINB);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
175  try{bfp->GetNamedData(times_temps[i], "E:NSLINC",&NSLINC);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
176  try{bfp->GetNamedData(times_temps[i], "E:NSLIND",&NSLIND);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
177  try{bfp->GetNamedData(times_temps[i], "E:TOR101",&TOR101);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
178  try{bfp->GetNamedData(times_temps[i], "E:TORTGT",&TORTGT);}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
179  // BPM Positions and Intensities - they each have 7 elements
180  // First is an average value of a few batches (often 2,3,4)
181  // used for auto-tuning, so we should disregard it
182 
183  try{HP121 = bfp->GetNamedVector(times_temps[i], "E:HP121[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " << "got exception: " << we.what() << "\n";}
184  try{VP121 = bfp->GetNamedVector(times_temps[i], "E:VP121[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
185  try{HPTGT = bfp->GetNamedVector(times_temps[i], "E:HPTGT[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
186  try{VPTGT = bfp->GetNamedVector(times_temps[i], "E:VPTGT[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
187  try{HITGT = bfp->GetNamedVector(times_temps[i], "E:HITGT[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
188  try{VITGT = bfp->GetNamedVector(times_temps[i], "E:VITGT[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
189  try{MTGTDS = bfp->GetNamedVector(times_temps[i], "E:MTGTDS[]");}catch (WebAPIException &we) {std::cout << "At time : " << times_temps[i] << " " <<"got exception: " << we.what() << "\n";}
190 
191  std::cout << "Finished getting IFBeam info" << std::endl;
192  std::cout << "BFP Time: " << times_temps[i] << " TOROID Time: " << TRTGTD_time << " TOROID COUNT: " << TRTGTD << std::endl;
193  unsigned long int time_closest_int = (int) TRTGTD_time;
194  double time_closest_ns = (TRTGTD_time - time_closest_int)*1000000000;
195 
196  sbn::NuMISpillInfo NuMIbeamInfo;
197  NuMIbeamInfo.TORTGT = TORTGT*1e12; //include factor of 1e12 protons in POT calculation
198  NuMIbeamInfo.TOR101 = TOR101*1e12; //include factor of 1e12 protons in POT calculation
199  NuMIbeamInfo.TRTGTD = TRTGTD*1e12; //include factor of 1e12 protons in POT calculation
200  NuMIbeamInfo.TR101D = TR101D*1e12; //include factor of 1e12 protons in POT calculation
201  NuMIbeamInfo.HRNDIR = HRNDIR;
202  NuMIbeamInfo.NSLINA = NSLINA;
203  NuMIbeamInfo.NSLINB = NSLINB;
204  NuMIbeamInfo.NSLINC = NSLINC;
205  NuMIbeamInfo.NSLIND = NSLIND;
206 
207  NuMIbeamInfo.HP121 = HP121;
208  NuMIbeamInfo.VP121 = VP121;
209  NuMIbeamInfo.HPTGT = HPTGT;
210  NuMIbeamInfo.VPTGT = VPTGT;
211  NuMIbeamInfo.HITGT = HITGT;
212  NuMIbeamInfo.VITGT = VITGT;
213  NuMIbeamInfo.MTGTDS = MTGTDS;
214 
215  NuMIbeamInfo.time = times_temps[i];
216  NuMIbeamInfo.event = e.event();
217  NuMIbeamInfo.spill_time_s = time_closest_int;
218  NuMIbeamInfo.spill_time_ns = time_closest_ns;
219  // Save the Number of DAQ Gates in the first saved spill
220  if (spill_count == 1) {
221  NuMIbeamInfo.daq_gates = number_of_gates_since_previous_event;
222  }
223  else {
224  NuMIbeamInfo.daq_gates = 0;
225  }
226 
227  fOutbeamInfos.push_back(NuMIbeamInfo);
228  }
229  if(spill_count > number_of_gates_since_previous_event)
230  std::cout << "Event Spills : " << spill_count << ", DAQ Spills : " << number_of_gates_since_previous_event << " \t \t ::: WRONG!"<< std::endl;
231  else
232  std::cout << "Event Spills : " << spill_count << ", DAQ Spills : " << number_of_gates_since_previous_event << std::endl;
233  }
234 }
sbn::NuMISpillInfo::TOR101
float TOR101
Other monitored POT at Start.
Definition: NuMISpillInfo.h:30
sbn::NuMISpillInfo::HP121
std::vector< double > HP121
List of Horizontal position at start per-bunch.
Definition: NuMISpillInfo.h:15
sbn::NuMISpillInfo::TR101D
float TR101D
Best monitored POT at Start.
Definition: NuMISpillInfo.h:28
sbn::NuMISpillInfo::TRTGTD
float TRTGTD
Best monitored POT at Target.
Definition: NuMISpillInfo.h:27
sbn::NuMISpillInfo::time
float time
Time of device used to lookup spill.
Definition: NuMISpillInfo.h:31
sbn::NuMISpillInfo::daq_gates
unsigned int daq_gates
Definition: NuMISpillInfo.h:37
sbn::NuMISpillInfo::HPTGT
std::vector< double > HPTGT
List of Horizontal position at target per-bunch.
Definition: NuMISpillInfo.h:17
sbn::NuMIRetriever::bfp
std::unique_ptr< ifbeam_ns::BeamFolder > bfp
Definition: NuMIRetriever_module.cc:64
sbn::NuMISpillInfo::TORTGT
float TORTGT
Other monitored POT at Target.
Definition: NuMISpillInfo.h:29
sbn::NuMISpillInfo::spill_time_s
unsigned long int spill_time_s
The IFDB Beam Spill Time, unit sec.
Definition: NuMISpillInfo.h:33
sbn::NuMISpillInfo::VITGT
std::vector< double > VITGT
List of Vertical-Monitor intensity at tartget per bunch.
Definition: NuMISpillInfo.h:20
sbn::NuMIRetriever::fDeviceUsedForTiming
std::string fDeviceUsedForTiming
Definition: NuMIRetriever_module.cc:61
sbn::NuMISpillInfo::MTGTDS
std::vector< double > MTGTDS
Definition: NuMISpillInfo.h:21
sbn::NuMISpillInfo::HITGT
std::vector< double > HITGT
List of Horizontal-Monitor intensity at target per-bunch.
Definition: NuMISpillInfo.h:19
e
do i e
Definition: updateVersionList.sh:68
sbn::NuMISpillInfo::spill_time_ns
unsigned long int spill_time_ns
The IFDB Beam Spill Time, unit nsec.
Definition: NuMISpillInfo.h:34
sbn::NuMIRetriever::fOutbeamInfos
std::vector< sbn::NuMISpillInfo > fOutbeamInfos
Definition: NuMIRetriever_module.cc:55
sbn::NuMISpillInfo::event
unsigned int event
Definition: NuMISpillInfo.h:36
sbn::NuMISpillInfo::VP121
std::vector< double > VP121
List of Vertical position at start per-bunch.
Definition: NuMISpillInfo.h:16
sbn::NuMISpillInfo::VPTGT
std::vector< double > VPTGT
List of Vertical position at target per-bunch.
Definition: NuMISpillInfo.h:18
cout
BEGIN_PROLOG could also be cout
Definition: messageservice.fcl:10

Member Data Documentation

std::unique_ptr<ifbeam_ns::BeamFolder> sbn::NuMIRetriever::bfp
private

Definition at line 64 of file NuMIRetriever_module.cc.

double sbn::NuMIRetriever::fBFPEpsilion
private

Definition at line 57 of file NuMIRetriever_module.cc.

std::string sbn::NuMIRetriever::fDeviceUsedForTiming
private

Definition at line 61 of file NuMIRetriever_module.cc.

std::vector< sbn::NuMISpillInfo > sbn::NuMIRetriever::fOutbeamInfos
private

Definition at line 55 of file NuMIRetriever_module.cc.

double sbn::NuMIRetriever::fTimePad
private

Definition at line 56 of file NuMIRetriever_module.cc.

std::string sbn::NuMIRetriever::fURL
private

Definition at line 58 of file NuMIRetriever_module.cc.

art::ServiceHandle<ifbeam_ns::IFBeam> sbn::NuMIRetriever::ifbeam_handle
private

Definition at line 63 of file NuMIRetriever_module.cc.

std::string sbn::NuMIRetriever::raw_data_label_
private

Definition at line 60 of file NuMIRetriever_module.cc.

int sbn::NuMIRetriever::TotalBeamSpills
private

Definition at line 62 of file NuMIRetriever_module.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.5