BeamlineSimInput_module.cc
Go to the documentation of this file.
1  ////////////////////////////////////////////////////////////////////////////
2 /// \file BeamlineSimInput_module.cc
3 /// \module producer
4 /// \brief Module to add NOvASoft-formatted simulated beamline
5 /// information to the corresponding detector simulation.
6 /// \author Mike Wallbank (University of Cincinnati) <wallbank@fnal.gov>
7 /// Abhilash Y D (Syracuse University) <ayallapp@syr.edu>
8 /// \date January 2019
9 ////////////////////////////////////////////////////////////////////////////
10 
11 ////////////////////////////////////////////////////////////////////////////////////
12 /// What this module's all about.
13 ///
14 /// Beamline simulation:
15 /// Junting's beamline generator are used as input to the NOvA GEANT4 simulation,
16 /// which propagate the particles from the end of the beamline through the
17 /// detector.
18 /// The information stored in the beamline simulation files contain properties of
19 /// the particles as they traverse the beamline. This information isn't included
20 /// in the NOvASoft simulation files.
21 ///
22 /// BeamlineSimInput:
23 /// Designed to reconcile the information from the beamline simulation with the
24 /// detector simulation.
25 /// This module runs over the detector simulation files, which should contain
26 /// art::Events corresponding to the input particles from the end of the beamline.
27 /// It also takes the beamline output as input, reformats the data into the
28 /// NOvASoft data products, and adds them to the event.
29 ////////////////////////////////////////////////////////////////////////////////////
30 
31 // framework
32 #include "art/Framework/Core/EDAnalyzer.h"
33 #include "art/Framework/Core/EDProducer.h"
34 #include "art/Framework/Core/ModuleMacros.h"
35 #include "art/Framework/Principal/Event.h"
36 #include "fhiclcpp/ParameterSet.h"
37 #include "art/Framework/Principal/Handle.h"
38 #include "canvas/Persistency/Common/Ptr.h"
39 #include "canvas/Persistency/Common/PtrVector.h"
40 #include "art/Framework/Services/Registry/ServiceHandle.h"
41 #include "art/Framework/Services/Optional/TFileService.h"
42 #include "art/Framework/Services/Optional/TFileDirectory.h"
43 #include "canvas/Persistency/Common/FindManyP.h"
44 
45 // nova
46 #include "RawData/RawBeamline.h"
47 #include "BeamlineRecoBase/BeamlineDigit.h"
48 #include "BeamlineRecoBase/BeamlinePID.h"
49 #include "BeamlineRecoBase/MuonStackTrack.h"
50 #include "BeamlineRecoBase/Cherenkov.h"
51 #include "BeamlineRecoBase/ToF.h"
52 #include "BeamlineRecoBase/WCTrack.h"
53 #include "Utilities/AssociationUtil.h"
54 #include "nusimdata/SimulationBase/MCTruth.h"
55 
56 // stl
57 #include <iostream>
58 
59 // root
60 #include "TTree.h"
61 #include "TFile.h"
62 #include "TMath.h"
63 #include <TTreeReader.h>
64 #include <TTreeReaderValue.h>
65 #include <TTreeReaderArray.h>
66 #include <fstream>
67 #include <vector>
68 #include <utility>
69 
70 // -----------------------------------------------------------------------
71 namespace beamlinesim {
72 
73  class BeamlineSimInput : public art::EDProducer {
74 
75  public:
76 
77  BeamlineSimInput(const fhicl::ParameterSet& pset);
78 
79  void reconfigure(const fhicl::ParameterSet& pset);
80  void produce(art::Event& evt);
81  //void beginRun(art::Run &r);
82  void beginJob();
83  //void beginSubRun(const art::SubRun& subrun);
84  float_t bltof(float_t tof_us, float_t tof_ds, int presenttofus, int presenttofds) ;
85  TTreeReader reader;
86  TTree *tree;
87  int particle_no;
88  int event_no;
89  int final_index;
90  int total_particles;
91 
92  private:
93  // configuration
94  art::InputTag fRawToFDataLabel;
95 
96  // beamline file
97  TFile* fBeamlineSimFile;
98  TTree* fBeamlineSimTree;
99 
100  // not all beamline particles make an art::Event
101  // keep a track
102  int fBeamlineTreeIndex;
103 
104  // particle properties
105 
106  //std::vector<int> veventid;
107  //std::vector<float> vx;
108  //int feventid;
109 
110  //std::vector<int> visnoise;
111  std::vector<unsigned int> vPresentWire_Chamber1;
112  std::vector<unsigned int> vPresentWire_Chamber2;
113  std::vector<unsigned int> vPresentWire_Chamber3;
114  std::vector<unsigned int> vPresentWire_Chamber4;
115  std::vector<unsigned int> vPresenttofus;
116  std::vector<unsigned int> vPresenttofds;
117  std::vector<unsigned int> vPresentcherenkov;
118  std::vector<float_t> vXWire_Chamber1;
119  std::vector<float_t> vYWire_Chamber1;
120  std::vector<float_t> vZWire_Chamber1;
121  std::vector<float_t> vTWire_Chamber1;
122  std::vector<float_t> vXWire_Chamber2;
123  std::vector<float_t> vYWire_Chamber2;
124  std::vector<float_t> vZWire_Chamber2;
125  std::vector<float_t> vTWire_Chamber2;
126  std::vector<float_t> vXWire_Chamber3;
127  std::vector<float_t> vYWire_Chamber3;
128  std::vector<float_t> vZWire_Chamber3;
129  std::vector<float_t> vTWire_Chamber3;
130  std::vector<float_t> vXWire_Chamber4;
131  std::vector<float_t> vYWire_Chamber4;
132  std::vector<float_t> vZWire_Chamber4;
133  std::vector<float_t> vTWire_Chamber4;
134 // std::vector<float_t> vPxWire_Chamber1; //momentums
135 // std::vector<float_t> vPyWire_Chamber1;
136 // std::vector<float_t> vPzWire_Chamber1;
137 // std::vector<float_t> vPxWire_Chamber2;
138 // std::vector<float_t> vPyWire_Chamber2;
139 // std::vector<float_t> vPzWire_Chamber2;
140 // std::vector<float_t> vPxWire_Chamber3;
141 // std::vector<float_t> vPyWire_Chamber3;
142 // std::vector<float_t> vPzWire_Chamber3;
143 // std::vector<float_t> vPxWire_Chamber4;
144 // std::vector<float_t> vPyWire_Chamber4;
145 // std::vector<float_t> vPzWire_Chamber4;
146  std::vector<float_t> vttof_us;//tofs
147  std::vector<float_t> vttof_ds;
148  std::vector<float_t> vtcherenkov;
149  std::vector<float_t> vxnova;
150  std::vector<float_t> vynova;
151  std::vector<float_t> vznova;
152  std::vector<float_t> vtnova;
153  //std::vector<int> vspill_number;
154 
155 
156  //int fisnoise;
157  int fPresentWire_Chamber1;
158  int fPresentWire_Chamber2;
159  int fPresentWire_Chamber3;
160  int fPresentWire_Chamber4;
161  int fPresenttofus;
162  int fPresenttofds;
163  int fPresentcherenkov;
164  float_t fXWire_Chamber1;
165  float_t fYWire_Chamber1;
166  float_t fZWire_Chamber1;
167  float_t fTWire_Chamber1;
168  float_t fXWire_Chamber2;
169  float_t fYWire_Chamber2;
170  float_t fZWire_Chamber2;
171  float_t fTWire_Chamber2;
172  float_t fXWire_Chamber3;
173  float_t fYWire_Chamber3;
174  float_t fZWire_Chamber3;
175  float_t fTWire_Chamber3;
176  float_t fXWire_Chamber4;
177  float_t fYWire_Chamber4;
178  float_t fZWire_Chamber4;
179  float_t fTWire_Chamber4;
180 // float_t fPxWire_Chamber1; //momentums
181 // float_t fPyWire_Chamber1;
182 // float_t fPzWire_Chamber1;
183 // float_t fPxWire_Chamber2;
184 // float_t fPyWire_Chamber2;
185 // float_t fPzWire_Chamber2;
186 // float_t fPxWire_Chamber3;
187 // float_t fPyWire_Chamber3;
188 // float_t fPzWire_Chamber3;
189 // float_t fPxWire_Chamber4;
190 // float_t fPyWire_Chamber4;
191 // float_t fPzWire_Chamber4;
192  float_t fttof_us;//tofs
193  float_t fttof_ds;
194  float_t ftcherenkov;
195 // float_t fxnova;
196 // float_t fynova;
197 // float_t fznova;
198 // float_t ftnova;
199  //int fspill_number;
200  };
201 }
202 
203 // -----------------------------------------------------------------------
204 void beamlinesim::BeamlineSimInput::beginJob()
205 {
206  //beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/TestBeam/beamlineinput_new.root","READ");
207  //beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/junting/beamlinesimulation_wohoo.root","READ");
208  beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/junting/beamlinesimulation_wohoo.root","READ");
209  TTreeReader reader("tree", fBeamlineSimFile);
210  //TTreeReaderValue<int> feventid(reader, "event_id");
211  //TTreeReaderValue<float_t> fx(reader, "x");
212  //TTreeReaderValue<int> fisnoise(reader, "is_noise");
213  TTreeReaderValue<unsigned int> fPresentWire_Chamber1(reader, "present_wire_chamber_1");
214  TTreeReaderValue<unsigned int> fPresentWire_Chamber2(reader, "present_wire_chamber_2");
215  TTreeReaderValue<unsigned int> fPresentWire_Chamber3(reader, "present_wire_chamber_3");
216  TTreeReaderValue<unsigned int> fPresentWire_Chamber4(reader, "present_wire_chamber_4");
217  TTreeReaderValue<unsigned int> fPresenttofus(reader, "present_tof_us");
218  TTreeReaderValue<unsigned int> fPresenttofds(reader, "present_tof_ds");
219  TTreeReaderValue<unsigned int> fPresentcherenkov(reader, "present_cerenkov");
220  TTreeReaderValue<float> fXWire_Chamber1(reader, "x_wire_chamber_1");
221  TTreeReaderValue<float> fYWire_Chamber1(reader, "y_wire_chamber_1");
222  TTreeReaderValue<float> fZWire_Chamber1(reader, "z_wire_chamber_1");
223  TTreeReaderValue<float> fTWire_Chamber1(reader, "t_wire_chamber_1");
224  TTreeReaderValue<float> fXWire_Chamber2(reader, "x_wire_chamber_2");
225  TTreeReaderValue<float> fYWire_Chamber2(reader, "y_wire_chamber_2");
226  TTreeReaderValue<float> fZWire_Chamber2(reader, "z_wire_chamber_2");
227  TTreeReaderValue<float> fTWire_Chamber2(reader, "t_wire_chamber_2");
228  TTreeReaderValue<float> fXWire_Chamber3(reader, "x_wire_chamber_3");
229  TTreeReaderValue<float> fYWire_Chamber3(reader, "y_wire_chamber_3");
230  TTreeReaderValue<float> fZWire_Chamber3(reader, "z_wire_chamber_3");
231  TTreeReaderValue<float> fTWire_Chamber3(reader, "t_wire_chamber_3");
232  TTreeReaderValue<float> fXWire_Chamber4(reader, "x_wire_chamber_4");
233  TTreeReaderValue<float> fYWire_Chamber4(reader, "y_wire_chamber_4");
234  TTreeReaderValue<float> fZWire_Chamber4(reader, "z_wire_chamber_4");
235  TTreeReaderValue<float> fTWire_Chamber4(reader, "t_wire_chamber_4");
236 // TTreeReaderValue<float> fPxWire_Chamber1(reader, "px_wire_chamber_1");
237 // TTreeReaderValue<float> fPyWire_Chamber1(reader, "py_wire_chamber_1");
238 // TTreeReaderValue<float> fPzWire_Chamber1(reader, "pz_wire_chamber_1");
239 // TTreeReaderValue<float> fPxWire_Chamber2(reader, "px_wire_chamber_2");
240 // TTreeReaderValue<float> fPyWire_Chamber2(reader, "py_wire_chamber_2");
241 // TTreeReaderValue<float> fPzWire_Chamber2(reader, "pz_wire_chamber_2");
242 // TTreeReaderValue<float> fPxWire_Chamber3(reader, "px_wire_chamber_3");
243 // TTreeReaderValue<float> fPyWire_Chamber3(reader, "py_wire_chamber_3");
244 // TTreeReaderValue<float> fPzWire_Chamber3(reader, "pz_wire_chamber_3");
245 // TTreeReaderValue<float> fPxWire_Chamber4(reader, "px_wire_chamber_4");
246 // TTreeReaderValue<float> fPyWire_Chamber4(reader, "py_wire_chamber_4");
247 // TTreeReaderValue<float> fPzWire_Chamber4(reader, "pz_wire_chamber_4");
248  TTreeReaderValue<float> fttof_us(reader, "t_tof_us");
249  TTreeReaderValue<float> fttof_ds(reader, "t_tof_ds");
250  TTreeReaderValue<float> ftcherenkov(reader, "t_cerenkov");
251 // TTreeReaderValue<float> fxnova(reader, "x");
252 // TTreeReaderValue<float> fynova(reader, "y");
253 // TTreeReaderValue<float> fznova(reader, "z");
254 // TTreeReaderValue<float> ftnova(reader, "t");
255  //TTreeReaderValue<int> fspill_number(reader, "spill_number");
256 
257  while (reader.Next())
258  {
259  //veventid.push_back(*feventid);
260  //vx.push_back(*fx);
261 
262  //visnoise.push_back(*fisnoise);
263  vPresentWire_Chamber1.push_back(*fPresentWire_Chamber1);
264  vPresentWire_Chamber2.push_back(*fPresentWire_Chamber2);
265  vPresentWire_Chamber3.push_back(*fPresentWire_Chamber3);
266  vPresentWire_Chamber4.push_back(*fPresentWire_Chamber4);
267  vPresenttofus.push_back(*fPresenttofus);
268  vPresenttofds.push_back(*fPresenttofds);
269  vPresentcherenkov.push_back(*fPresentcherenkov);
270  vXWire_Chamber1.push_back(*fXWire_Chamber1);
271  vYWire_Chamber1.push_back(*fYWire_Chamber1);
272  vZWire_Chamber1.push_back(*fZWire_Chamber1);
273  vTWire_Chamber1.push_back(*fTWire_Chamber1);
274  vXWire_Chamber2.push_back(*fXWire_Chamber2);
275  vYWire_Chamber2.push_back(*fYWire_Chamber2);
276  vZWire_Chamber2.push_back(*fZWire_Chamber2);
277  vTWire_Chamber2.push_back(*fTWire_Chamber2);
278  vXWire_Chamber3.push_back(*fXWire_Chamber3);
279  vYWire_Chamber3.push_back(*fYWire_Chamber3);
280  vZWire_Chamber3.push_back(*fZWire_Chamber3);
281  vTWire_Chamber3.push_back(*fTWire_Chamber3);
282  vXWire_Chamber4.push_back(*fXWire_Chamber4);
283  vYWire_Chamber4.push_back(*fYWire_Chamber4);
284  vZWire_Chamber4.push_back(*fZWire_Chamber4);
285  vTWire_Chamber4.push_back(*fTWire_Chamber4);
286 // vPxWire_Chamber1.push_back(*fPxWire_Chamber1);
287 // vPyWire_Chamber1.push_back(*fPyWire_Chamber1);
288 // vPzWire_Chamber1.push_back(*fPzWire_Chamber1);
289 // vPxWire_Chamber2.push_back(*fPxWire_Chamber2);
290 // vPyWire_Chamber2.push_back(*fPyWire_Chamber2);
291 // vPzWire_Chamber2.push_back(*fPzWire_Chamber2);
292 // vPxWire_Chamber3.push_back(*fPxWire_Chamber3);
293 // vPyWire_Chamber3.push_back(*fPyWire_Chamber3);
294 // vPzWire_Chamber3.push_back(*fPzWire_Chamber3);
295 // vPxWire_Chamber4.push_back(*fPxWire_Chamber4);
296 // vPyWire_Chamber4.push_back(*fPyWire_Chamber4);
297 // vPzWire_Chamber4.push_back(*fPzWire_Chamber4);
298  vttof_us.push_back(*fttof_us);
299  vttof_ds.push_back(*fttof_ds);
300  vtcherenkov.push_back(*ftcherenkov);
301 // vxnova.push_back(*fxnova);
302 // vynova.push_back(*fynova);
303 // vznova.push_back(*fznova);
304 // vtnova.push_back(*ftnova);
305  //vspill_number.push_back(*fspill_number);
306  }
307 
308  //std::cout<<"eventid (1) " << veventid[1]<<std::endl;
309  //std::cout<<"eventid (2) " << veventid[2]<<std::endl;
310  //std::cout<<"x (1) " << vx[1]<<std::endl;
311  //std::cout<<"x (2) " << vx[2]<<std::endl;
312  return;
313 }
314 
315 beamlinesim::BeamlineSimInput::BeamlineSimInput(const fhicl::ParameterSet& pset) {
316  produces<std::vector<brb::ToF> >();
317  produces<std::vector<brb::Cherenkov> >();
318  produces<std::vector<rawdata::RawBeamlineWC> >();
319  total_particles = 0;
320  final_index = 0;
321  this->reconfigure(pset);
322 }
323 
324 // -----------------------------------------------------------------------
325 void beamlinesim::BeamlineSimInput::reconfigure(const fhicl::ParameterSet& pset) {
326 }
327 
328 // -----------------------------------------------------------------------
329 void beamlinesim::BeamlineSimInput::produce(art::Event& evt) {
330 
331  art::Handle<std::vector<simb::MCTruth> > mcTruthHandle;
332  std::vector<art::Ptr<simb::MCTruth> > mcTruths;
333  if (evt.getByLabel("generator", mcTruthHandle))
334  art::fill_ptr_vector(mcTruths, mcTruthHandle);
335 
336  std::cout << "There are " << mcTruths.size() << " mc truths in event " << evt.event() << std::endl;
337  particle_no=0;
338  event_no=evt.event();
339  std::cout<<"event number "<<event_no<<std::endl;
340  for (std::vector<art::Ptr<simb::MCTruth> >::const_iterator mcTruthIt = mcTruths.begin(); mcTruthIt != mcTruths.end(); ++mcTruthIt)
341  {
342  std::cout << " There are " << (*mcTruthIt)->NParticles() << " particles in this MCTruth." << std::endl;
343  particle_no=(*mcTruthIt)->NParticles();
344  }
345 
346  std::unique_ptr<std::vector<brb::ToF> > TOFptr(new std::vector<brb::ToF>);
347  std::unique_ptr<std::vector<brb::Cherenkov> > CHKVptr(new std::vector<brb::Cherenkov>);
348  std::unique_ptr<std::vector<rawdata::RawBeamlineWC> > WCptr(new std::vector<rawdata::RawBeamlineWC>);
349 
350  if(evt.event()!=6964)
351  {
352 
353  for (int i=0;i<particle_no;i++)
354  {
355  final_index = total_particles ;//+ i;
356  //total_particles=final_index;
357 // std::cout<<"putting particle number "<<i+1 <<" in the event"<<std::endl;
358 // std::cout<<"totalparticles is "<<total_particles<<std::endl;
359 // std::cout<<"final index is "<<final_index<<std::endl;
360  // reformat into novasoft data types
361  //Time of flight
362  brb::ToF tofobject;
363  std::pair <float,float> tofpair;
364  if (vPresenttofus[final_index]==1)
365  {
366  tofpair.first=vttof_us[final_index];
367  }
368  if ( vPresenttofds[final_index]==1)
369  {
370  tofpair.second=vttof_ds[final_index];
371  }
372  if ( vPresenttofus[final_index]==1 && vPresenttofds[final_index]==1)
373  {
374  tofobject.SetTime(vttof_us[final_index]-vttof_ds[final_index]);
375  tofobject.SetTimestamps(tofpair);
376  TOFptr->push_back(tofobject);
377  std::cout<<"the tof_us added to the event is "<< vttof_us[final_index] <<std::endl;
378  std::cout<<"the tof_ds added to the event is "<< vttof_ds[final_index] <<std::endl;
379  std::cout<<"the tof added to the event is "<< vttof_us[final_index]-vttof_ds[final_index] <<std::endl;
380  }
381  else
382  {
383  //std::cout<<"------------------------ didn't pass through both tof ------------------------------------------------------------------"<<std::endl;
384  }
385  //bool tofp=bool(vPresenttofus[final_index]) && bool(vPresenttofds[final_index]);
386  //std::cout<<"tof present up = "<< vPresenttofus[final_index] <<std::endl;
387  //std::cout<<"tof present up = "<< vPresenttofds[final_index] <<std::endl;
388 
389  //cherenkov
390 
391  brb::Cherenkov cherenkov_object;
392  if(vPresentcherenkov[final_index]==1){
393  cherenkov_object.SetTime(vtcherenkov[final_index]);
394  CHKVptr->push_back(cherenkov_object);
395  //std::cout<<"the cherenkov timestamp added to the event is "<< vtcherenkov[final_index] <<std::endl;
396  }
397  else
398  {
399  //std::cout<<"-------------------------- didn't pass through cherenkov-----------------------------------------------------------------"<<std::endl;
400  }
401 
402  //std::cout<<"cherenkov present = "<< vPresentcherenkov[final_index] <<std::endl;
403 
404  //wirechambers
405 
406  std::vector<rawdata::RawBeamlineWC> WCobject;
407 
408  std::vector<double> fchannelx;
409  std::vector<double> fchannely;
410  std::vector<uint64_t> fWCtime;
411  std::vector<unsigned int> fWCpresent;
412  std::vector<double> xWC;
413  std::vector<double> yWC;
414 
415  fchannelx.push_back(round((-(-491.51-vXWire_Chamber1[final_index])+64)/(TMath::Cos(16*TMath::Pi()/180))));
416  fchannelx.push_back(round((-(-910.57-vXWire_Chamber2[final_index])+64)/(TMath::Cos(16*TMath::Pi()/180))));
417  fchannelx.push_back(round(-(-1354.35-vXWire_Chamber3[final_index])+64));
418  fchannelx.push_back(round(-(-1354.35-vXWire_Chamber4[final_index])+64));
419 
420  fchannely.push_back(round(vYWire_Chamber1[final_index]+64));
421  fchannely.push_back(round(vYWire_Chamber2[final_index]+64));
422  fchannely.push_back(round(vYWire_Chamber3[final_index]+64));
423  fchannely.push_back(round(vYWire_Chamber4[final_index]+64));
424 
425  fWCtime.push_back(vTWire_Chamber1[final_index]);
426  fWCtime.push_back(vTWire_Chamber2[final_index]);
427  fWCtime.push_back(vTWire_Chamber3[final_index]);
428  fWCtime.push_back(vTWire_Chamber4[final_index]);
429 
430  fWCpresent.push_back(vPresentWire_Chamber1[final_index]);
431  fWCpresent.push_back(vPresentWire_Chamber2[final_index]);
432  fWCpresent.push_back(vPresentWire_Chamber3[final_index]);
433  fWCpresent.push_back(vPresentWire_Chamber4[final_index]);
434 
435  //just for debugging--------------------------------------------------
436  xWC.push_back(vXWire_Chamber1[final_index]);
437  xWC.push_back(vXWire_Chamber2[final_index]);
438  xWC.push_back(vXWire_Chamber3[final_index]);
439  xWC.push_back(vXWire_Chamber4[final_index]);
440 
441  yWC.push_back(vYWire_Chamber1[final_index]);
442  yWC.push_back(vYWire_Chamber2[final_index]);
443  yWC.push_back(vYWire_Chamber3[final_index]);
444  yWC.push_back(vYWire_Chamber4[final_index]);
445  //--mike helping me debug----------
446 // for (unsigned int wcNum = 0; wcNum < 4; ++wcNum) {
447 //
448 // rawdata::RawBeamlineWC::WCPulse xpulse(fchannelx[wcNum], fWCtime[wcNum]);
449 // rawdata::RawBeamlineWC::WCPulse ypulse(fchannely[wcNum], fWCtime[wcNum]);
450 //
451 // rawdata::RawBeamlineWC wc;//(wcNum, Form("WC%d", wcNum));
452 // wc.SetWCNumber(wcNum+1);
453 // //if (fWCpresent[wcNum]) {
454 // wc.AddXPulse(xpulse);
455 // wc.AddYPulse(ypulse);
456 // //}
457 //
458 // WCptr->push_back(wc);
459 //
460 // }
461  //------------------------------------
462 
463  for(unsigned int k=0;k<4;k++)
464  {
465  // Detector ID for this wire chamber
466  beamlinegeo::DetectorID detectorID(beamlinegeo::SystemID(beamlinegeo::BeamlineSystem::WC), k);
467 
468  rawdata::RawBeamlineWC::WCPulse xpulse(fchannelx[k],fWCtime[k]);
469  rawdata::RawBeamlineWC::WCPulse ypulse(fchannely[k],fWCtime[k]);
470  rawdata::RawBeamlineWC WCobject(detectorID);
471  WCobject.AddXPulse(xpulse);
472  WCobject.AddYPulse(ypulse);
473 
474  if (fWCpresent[k]==1) {
475 
476  std::cout<<"the x wire number for WC "<< k+1 <<" added to the event is "<< fchannelx[k]<<" and the x wire coordinate is "<<xWC[k] <<std::endl;
477  std::cout<<"the y wire number for WC "<< k+1 <<" added to the event is "<< fchannely[k] <<" and the y wire coordinate is "<<yWC[k] <<std::endl;
478  WCptr->push_back(WCobject);
479  }
480  else{
481  std::cout<<"-------------------- No WC "<<k+1<<" objects added ---------------------------"<<std::endl;
482  }
483 // std::cout<<"the x wire number for WC "<< k+1 <<" added to the event is "<< fchannelx[k] <<std::endl;
484 // std::cout<<"the x wire coordinate for WC "<<k+1<<" added to the event is "<<xWC[k] << std::endl;
485 // std::cout<<"the y wire number for WC "<< k+1 <<" added to the event is "<< fchannely[k] <<std::endl;
486 // std::cout<<"the y wire coordinate for WC "<<k+1<<" added to the event is "<<yWC[k] << std::endl;
487 // std::cout<<"wire chamber "<<k+1<<" present "<<fWCpresent[k] << std::endl;
488 
489  }
490 
491 
492  fchannelx.clear();
493  fchannely.clear();
494  fWCtime.clear();
495  fWCpresent.clear();
496  xWC.clear();
497  yWC.clear();
498  total_particles++;
499 
500  }
501 
502  evt.put(std::move(TOFptr));
503  evt.put(std::move(CHKVptr));
504  evt.put(std::move(WCptr));
505 
506  }
507  return;
508 
509 }
510 DEFINE_ART_MODULE(beamlinesim::BeamlineSimInput)
brb::ToF::SetTimestamps
void SetTimestamps(std::pair< float, float > timestamps)
Definition: ToF.cxx:47
brb::ToF::SetTime
void SetTime(float time)
Definition: ToF.cxx:41
beamlinesim::BeamlineSimInput::vPresentcherenkov
std::vector< unsigned int > vPresentcherenkov
Definition: BeamlineSimInput_module.cc:117
beamlinesim::BeamlineSimInput::vPresentWire_Chamber2
std::vector< unsigned int > vPresentWire_Chamber2
Definition: BeamlineSimInput_module.cc:112
beamlinesim::BeamlineSimInput::vPresentWire_Chamber3
std::vector< unsigned int > vPresentWire_Chamber3
Definition: BeamlineSimInput_module.cc:113
brb::Cherenkov::SetTime
void SetTime(float time)
Definition: Cherenkov.cxx:32
beamlinesim::BeamlineSimInput::vPresentWire_Chamber1
std::vector< unsigned int > vPresentWire_Chamber1
Definition: BeamlineSimInput_module.cc:111
DEFINE_ART_MODULE
DEFINE_ART_MODULE(TestTMapFile)
stan::math::round
fvar< T > round(const fvar< T > &x)
Definition: round.hpp:23
art::Event::put
ProductID put(std::unique_ptr< PROD > &&product)
Definition: Event.h:102
beamlinesim::BeamlineSimInput::BeamlineSimInput
BeamlineSimInput(const fhicl::ParameterSet &pset)
Definition: BeamlineSimInput_module.cc:315
BeamlineDigit.h
Encapsulation of reconstructed digitizer &#39;hits&#39;. Used for ToF PMTs and SiPMs, and Cherenkov and Muon ...
brb::ToF
Definition: ToF.h:25
beamlinesim::BeamlineSimInput::vPresentWire_Chamber4
std::vector< unsigned int > vPresentWire_Chamber4
Definition: BeamlineSimInput_module.cc:114
ToF.h
Encapsulation of reconstructed Time-of-Flight (ToF) information. Part of beamline reconstruction for ...
evt
int evt
Definition: FillTrainingTrees.h:22
WCTrack.h
Encapsulation of reconstructed Wire Chamber track. Part of beamline reconstruction for NOvA test beam...
art::Event::event
EventNumber_t event() const
Definition: Event.h:67
rawdata::RawBeamlineWC::AddXPulse
void AddXPulse(WCPulse pulse)
Definition: RawBeamline.cxx:199
om::cout
OStream cout
Definition: OStream.cxx:6
rawdata::RawBeamlineWC::AddYPulse
void AddYPulse(WCPulse pulse)
Definition: RawBeamline.cxx:205
beamlinesim::BeamlineSimInput::vPresenttofds
std::vector< unsigned int > vPresenttofds
Definition: BeamlineSimInput_module.cc:116
Event.h
beamlinesim::BeamlineSimInput::vPresenttofus
std::vector< unsigned int > vPresenttofus
Definition: BeamlineSimInput_module.cc:115
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
beamlinesim::BeamlineSimInput::reconfigure
void reconfigure(const fhicl::ParameterSet &pset)
Definition: BeamlineSimInput_module.cc:325
RawBeamline.h
Raw data definitions for beamline data used in NOvA test beam experiment.
Cherenkov.h
Encapsulation of &#39;reconstructed&#39; Cherenkov information. Part of beamline reconstruction for NOvA test...
art::fill_ptr_vector
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:464
BeamlinePID.h
Encapsulation of reconstructed PID information from detectors in the beamline (ToF, WCs, Cherenkov). Part of beamline reconstruction for NOvA test beam.
beamlinesim::BeamlineSimInput::bltof
float_t bltof(float_t tof_us, float_t tof_ds, int presenttofus, int presenttofds)
MuonStackTrack.h
Encapsulation of reconstructed track in the muon stack downstream of test beam detector. Part of beamline reconstruction for NOvA test beam.
Generated on Tue Dec 1 2020 05:14:22 for NOvA by   doxygen 1.8.11