beamlinesim::BeamlineSimInput Class Reference

Inheritance diagram for beamlinesim::BeamlineSimInput:

Public Types
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	BeamlineSimInput (const fhicl::ParameterSet &pset)
void	reconfigure (const fhicl::ParameterSet &pset)
void	produce (art::Event &evt)
void	beginJob ()
float_t	bltof (float_t tof_us, float_t tof_ds, int presenttofus, int presenttofds)
template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const
template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const
bool	modifiesEvent () const
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)
base_engine_t &	createEngine (seed_t seed)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Public Attributes
TTreeReader	reader
TTree *	tree
int	particle_no
int	event_no
int	final_index
int	total_particles

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Private Attributes
art::InputTag	fRawToFDataLabel
TFile *	fBeamlineSimFile
TTree *	fBeamlineSimTree
int	fBeamlineTreeIndex
std::vector< unsigned int >	vPresentWire_Chamber1
std::vector< unsigned int >	vPresentWire_Chamber2
std::vector< unsigned int >	vPresentWire_Chamber3
std::vector< unsigned int >	vPresentWire_Chamber4
std::vector< unsigned int >	vPresenttofus
std::vector< unsigned int >	vPresenttofds
std::vector< unsigned int >	vPresentcherenkov
std::vector< float_t >	vXWire_Chamber1
std::vector< float_t >	vYWire_Chamber1
std::vector< float_t >	vZWire_Chamber1
std::vector< float_t >	vTWire_Chamber1
std::vector< float_t >	vXWire_Chamber2
std::vector< float_t >	vYWire_Chamber2
std::vector< float_t >	vZWire_Chamber2
std::vector< float_t >	vTWire_Chamber2
std::vector< float_t >	vXWire_Chamber3
std::vector< float_t >	vYWire_Chamber3
std::vector< float_t >	vZWire_Chamber3
std::vector< float_t >	vTWire_Chamber3
std::vector< float_t >	vXWire_Chamber4
std::vector< float_t >	vYWire_Chamber4
std::vector< float_t >	vZWire_Chamber4
std::vector< float_t >	vTWire_Chamber4
std::vector< float_t >	vttof_us
std::vector< float_t >	vttof_ds
std::vector< float_t >	vtcherenkov
std::vector< float_t >	vxnova
std::vector< float_t >	vynova
std::vector< float_t >	vznova
std::vector< float_t >	vtnova
int	fPresentWire_Chamber1
int	fPresentWire_Chamber2
int	fPresentWire_Chamber3
int	fPresentWire_Chamber4
int	fPresenttofus
int	fPresenttofds
int	fPresentcherenkov
float_t	fXWire_Chamber1
float_t	fYWire_Chamber1
float_t	fZWire_Chamber1
float_t	fTWire_Chamber1
float_t	fXWire_Chamber2
float_t	fYWire_Chamber2
float_t	fZWire_Chamber2
float_t	fTWire_Chamber2
float_t	fXWire_Chamber3
float_t	fYWire_Chamber3
float_t	fZWire_Chamber3
float_t	fTWire_Chamber3
float_t	fXWire_Chamber4
float_t	fYWire_Chamber4
float_t	fZWire_Chamber4
float_t	fTWire_Chamber4
float_t	fttof_us
float_t	fttof_ds
float_t	ftcherenkov

Detailed Description

Definition at line 73 of file BeamlineSimInput_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

beamlinesim::BeamlineSimInput::BeamlineSimInput

(

const fhicl::ParameterSet &

pset

)

Definition at line 315 of file BeamlineSimInput_module.cc.

References final_index, reconfigure(), and total_particles.

                                                                           {
    produces<std::vector<brb::ToF> >();
    produces<std::vector<brb::Cherenkov> >();
    produces<std::vector<rawdata::RawBeamlineWC> >();
    total_particles = 0;
    final_index = 0;
    this->reconfigure(pset);
}

Member Function Documentation

void beamlinesim::BeamlineSimInput::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 204 of file BeamlineSimInput_module.cc.

References fBeamlineSimFile, fPresentcherenkov, fPresenttofds, fPresenttofus, fPresentWire_Chamber1, fPresentWire_Chamber2, fPresentWire_Chamber3, fPresentWire_Chamber4, ftcherenkov, fttof_ds, fttof_us, fTWire_Chamber1, fTWire_Chamber2, fTWire_Chamber3, fTWire_Chamber4, fXWire_Chamber1, fXWire_Chamber2, fXWire_Chamber3, fXWire_Chamber4, fYWire_Chamber1, fYWire_Chamber2, fYWire_Chamber3, fYWire_Chamber4, fZWire_Chamber1, fZWire_Chamber2, fZWire_Chamber3, fZWire_Chamber4, reader, vPresentcherenkov, vPresenttofds, vPresenttofus, vPresentWire_Chamber1, vPresentWire_Chamber2, vPresentWire_Chamber3, vPresentWire_Chamber4, vtcherenkov, vttof_ds, vttof_us, vTWire_Chamber1, vTWire_Chamber2, vTWire_Chamber3, vTWire_Chamber4, vXWire_Chamber1, vXWire_Chamber2, vXWire_Chamber3, vXWire_Chamber4, vYWire_Chamber1, vYWire_Chamber2, vYWire_Chamber3, vYWire_Chamber4, vZWire_Chamber1, vZWire_Chamber2, vZWire_Chamber3, and vZWire_Chamber4.

{
    //beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/TestBeam/beamlineinput_new.root","READ");
    //beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/junting/beamlinesimulation_wohoo.root","READ");
    beamlinesim::BeamlineSimInput::fBeamlineSimFile = new TFile("/nova/app/users/ayallapp/testbeam/test_testbeam/junting/beamlinesimulation_wohoo.root","READ");
    TTreeReader reader("tree", fBeamlineSimFile);
    //TTreeReaderValue<int> feventid(reader, "event_id");
    //TTreeReaderValue<float_t> fx(reader, "x");
    //TTreeReaderValue<int> fisnoise(reader, "is_noise");
    TTreeReaderValue<unsigned int> fPresentWire_Chamber1(reader, "present_wire_chamber_1");
    TTreeReaderValue<unsigned int> fPresentWire_Chamber2(reader, "present_wire_chamber_2");
    TTreeReaderValue<unsigned int> fPresentWire_Chamber3(reader, "present_wire_chamber_3");
    TTreeReaderValue<unsigned int> fPresentWire_Chamber4(reader, "present_wire_chamber_4");
    TTreeReaderValue<unsigned int> fPresenttofus(reader, "present_tof_us");
    TTreeReaderValue<unsigned int> fPresenttofds(reader, "present_tof_ds");
    TTreeReaderValue<unsigned int> fPresentcherenkov(reader, "present_cerenkov");
    TTreeReaderValue<float> fXWire_Chamber1(reader, "x_wire_chamber_1");
    TTreeReaderValue<float> fYWire_Chamber1(reader, "y_wire_chamber_1");
    TTreeReaderValue<float> fZWire_Chamber1(reader, "z_wire_chamber_1");
    TTreeReaderValue<float> fTWire_Chamber1(reader, "t_wire_chamber_1");
    TTreeReaderValue<float> fXWire_Chamber2(reader, "x_wire_chamber_2");
    TTreeReaderValue<float> fYWire_Chamber2(reader, "y_wire_chamber_2");
    TTreeReaderValue<float> fZWire_Chamber2(reader, "z_wire_chamber_2");
    TTreeReaderValue<float> fTWire_Chamber2(reader, "t_wire_chamber_2");
    TTreeReaderValue<float> fXWire_Chamber3(reader, "x_wire_chamber_3");
    TTreeReaderValue<float> fYWire_Chamber3(reader, "y_wire_chamber_3");
    TTreeReaderValue<float> fZWire_Chamber3(reader, "z_wire_chamber_3");
    TTreeReaderValue<float> fTWire_Chamber3(reader, "t_wire_chamber_3");
    TTreeReaderValue<float> fXWire_Chamber4(reader, "x_wire_chamber_4");
    TTreeReaderValue<float> fYWire_Chamber4(reader, "y_wire_chamber_4");
    TTreeReaderValue<float> fZWire_Chamber4(reader, "z_wire_chamber_4");
    TTreeReaderValue<float> fTWire_Chamber4(reader, "t_wire_chamber_4");
//    TTreeReaderValue<float> fPxWire_Chamber1(reader, "px_wire_chamber_1");
//    TTreeReaderValue<float> fPyWire_Chamber1(reader, "py_wire_chamber_1");
//    TTreeReaderValue<float> fPzWire_Chamber1(reader, "pz_wire_chamber_1");
//    TTreeReaderValue<float> fPxWire_Chamber2(reader, "px_wire_chamber_2");
//    TTreeReaderValue<float> fPyWire_Chamber2(reader, "py_wire_chamber_2");
//    TTreeReaderValue<float> fPzWire_Chamber2(reader, "pz_wire_chamber_2");
//    TTreeReaderValue<float> fPxWire_Chamber3(reader, "px_wire_chamber_3");
//    TTreeReaderValue<float> fPyWire_Chamber3(reader, "py_wire_chamber_3");
//    TTreeReaderValue<float> fPzWire_Chamber3(reader, "pz_wire_chamber_3");
//    TTreeReaderValue<float> fPxWire_Chamber4(reader, "px_wire_chamber_4");
//    TTreeReaderValue<float> fPyWire_Chamber4(reader, "py_wire_chamber_4");
//    TTreeReaderValue<float> fPzWire_Chamber4(reader, "pz_wire_chamber_4");
    TTreeReaderValue<float> fttof_us(reader, "t_tof_us");
    TTreeReaderValue<float> fttof_ds(reader, "t_tof_ds");
    TTreeReaderValue<float> ftcherenkov(reader, "t_cerenkov");
//    TTreeReaderValue<float> fxnova(reader, "x");
//    TTreeReaderValue<float> fynova(reader, "y");
//    TTreeReaderValue<float> fznova(reader, "z");
//    TTreeReaderValue<float> ftnova(reader, "t");
    //TTreeReaderValue<int> fspill_number(reader, "spill_number");
    
    while (reader.Next())
    {
        //veventid.push_back(*feventid);
        //vx.push_back(*fx);
        
        //visnoise.push_back(*fisnoise);
        vPresentWire_Chamber1.push_back(*fPresentWire_Chamber1);
        vPresentWire_Chamber2.push_back(*fPresentWire_Chamber2);
        vPresentWire_Chamber3.push_back(*fPresentWire_Chamber3);
        vPresentWire_Chamber4.push_back(*fPresentWire_Chamber4);
        vPresenttofus.push_back(*fPresenttofus);
        vPresenttofds.push_back(*fPresenttofds);
        vPresentcherenkov.push_back(*fPresentcherenkov);
        vXWire_Chamber1.push_back(*fXWire_Chamber1);
        vYWire_Chamber1.push_back(*fYWire_Chamber1);
        vZWire_Chamber1.push_back(*fZWire_Chamber1);
        vTWire_Chamber1.push_back(*fTWire_Chamber1);
        vXWire_Chamber2.push_back(*fXWire_Chamber2);
        vYWire_Chamber2.push_back(*fYWire_Chamber2);
        vZWire_Chamber2.push_back(*fZWire_Chamber2);
        vTWire_Chamber2.push_back(*fTWire_Chamber2);
        vXWire_Chamber3.push_back(*fXWire_Chamber3);
        vYWire_Chamber3.push_back(*fYWire_Chamber3);
        vZWire_Chamber3.push_back(*fZWire_Chamber3);
        vTWire_Chamber3.push_back(*fTWire_Chamber3);
        vXWire_Chamber4.push_back(*fXWire_Chamber4);
        vYWire_Chamber4.push_back(*fYWire_Chamber4);
        vZWire_Chamber4.push_back(*fZWire_Chamber4);
        vTWire_Chamber4.push_back(*fTWire_Chamber4);
//        vPxWire_Chamber1.push_back(*fPxWire_Chamber1);
//        vPyWire_Chamber1.push_back(*fPyWire_Chamber1);
//        vPzWire_Chamber1.push_back(*fPzWire_Chamber1);
//        vPxWire_Chamber2.push_back(*fPxWire_Chamber2);
//        vPyWire_Chamber2.push_back(*fPyWire_Chamber2);
//        vPzWire_Chamber2.push_back(*fPzWire_Chamber2);
//        vPxWire_Chamber3.push_back(*fPxWire_Chamber3);
//        vPyWire_Chamber3.push_back(*fPyWire_Chamber3);
//        vPzWire_Chamber3.push_back(*fPzWire_Chamber3);
//        vPxWire_Chamber4.push_back(*fPxWire_Chamber4);
//        vPyWire_Chamber4.push_back(*fPyWire_Chamber4);
//        vPzWire_Chamber4.push_back(*fPzWire_Chamber4);
        vttof_us.push_back(*fttof_us);
        vttof_ds.push_back(*fttof_ds);
        vtcherenkov.push_back(*ftcherenkov);
//        vxnova.push_back(*fxnova);
//        vynova.push_back(*fynova);
//        vznova.push_back(*fznova);
//        vtnova.push_back(*ftnova);
        //vspill_number.push_back(*fspill_number);
    }
    
    //std::cout<<"eventid (1) " << veventid[1]<<std::endl;
    //std::cout<<"eventid (2) " << veventid[2]<<std::endl;
    //std::cout<<"x (1) " << vx[1]<<std::endl;
    //std::cout<<"x (2) " << vx[2]<<std::endl;
    return;
}

float_t beamlinesim::BeamlineSimInput::bltof	(	float_t	tof_us,
		float_t	tof_ds,
		int	presenttofus,
		int	presenttofds
	)

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 146 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 161 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 171 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed )

inherited

Referenced by evgen::CosmicPionGen::CosmicPionGen(), evgen::CosmicsGen::CosmicsGen(), fuzz::FuzzyKVertex::FuzzyKVertex(), g4n::G4Gen::G4Gen(), photrans::ImprovedTransport::ImprovedTransport(), art::MixFilter< T >::initEngine_(), evgen::IsoGen::IsoGen(), rsim::MakeNoiseSpectrumFile::MakeNoiseSpectrumFile(), photrans::PhotonSystRewriter::PhotonSystRewriter(), photrans::PhotonTransport::PhotonTransport(), rsim::ReadoutSim::ReadoutSim(), evgen::SingleGen::SingleGen(), and evgen::SupernovaGen::SupernovaGen().

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

CurrentProcessingContext const* art::EDProducer::currentContext ( ) const

protectedinherited

seed_t art::EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset, char const  key[] = "seed", seed_t const  implicit_seed = -1  )

inherited

Referenced by art::MixFilter< T >::initEngine_().

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const &  instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

Referenced by skim::NueSkimmer::CopyMichelSlice(), and skim::NueSkimmer::CopyMichelTrack().

  {
    return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                               instanceName);
  }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID ( ModuleDescription const &  moduleDescription, std::string const &  instanceName  ) const

inherited

Definition at line 56 of file ProducerBase.h.

References art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

  {
    auto const& pd =
      get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
    return pd.productID();
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it )

inherited

Definition at line 189 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 204 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it )

inherited

Definition at line 214 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID(), and string.

    {
      return true;
    }

static cet::exempt_ptr<Consumer> art::Consumer::non_module_context ( )

staticinherited

Referenced by art::RootOutput::beginSubRun(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EmptyEvent::readRun_(), and art::EmptyEvent::readSubRun_().

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

void beamlinesim::BeamlineSimInput::produce ( art::Event &  evt )

virtual

Implements art::EDProducer.

Definition at line 329 of file BeamlineSimInput_module.cc.

References rawdata::RawBeamlineWC::AddXPulse(), rawdata::RawBeamlineWC::AddYPulse(), om::cout, DEFINE_ART_MODULE(), allTimeWatchdog::endl, art::Event::event(), event_no, art::fill_ptr_vector(), final_index, art::DataViewImpl::getByLabel(), MECModelEnuComparisons::i, particle_no, art::Event::put(), stan::math::round(), brb::Cherenkov::SetTime(), brb::ToF::SetTime(), brb::ToF::SetTimestamps(), total_particles, vPresentcherenkov, vPresenttofds, vPresenttofus, vPresentWire_Chamber1, vPresentWire_Chamber2, vPresentWire_Chamber3, vPresentWire_Chamber4, vtcherenkov, vttof_ds, vttof_us, vTWire_Chamber1, vTWire_Chamber2, vTWire_Chamber3, vTWire_Chamber4, vXWire_Chamber1, vXWire_Chamber2, vXWire_Chamber3, vXWire_Chamber4, vYWire_Chamber1, vYWire_Chamber2, vYWire_Chamber3, vYWire_Chamber4, and beamlinegeo::WC.

                                                       {
    
    art::Handle<std::vector<simb::MCTruth> > mcTruthHandle;
    std::vector<art::Ptr<simb::MCTruth> > mcTruths;
    if (evt.getByLabel("generator", mcTruthHandle))
        art::fill_ptr_vector(mcTruths, mcTruthHandle);
    
    std::cout << "There are " << mcTruths.size() << " mc truths in event " << evt.event() << std::endl;
    particle_no=0;
    event_no=evt.event();
    std::cout<<"event number "<<event_no<<std::endl;
    for (std::vector<art::Ptr<simb::MCTruth> >::const_iterator mcTruthIt = mcTruths.begin(); mcTruthIt != mcTruths.end(); ++mcTruthIt)
    {
        std::cout << "  There are " << (*mcTruthIt)->NParticles() << " particles in this MCTruth." << std::endl;
        particle_no=(*mcTruthIt)->NParticles();
    }
    
    std::unique_ptr<std::vector<brb::ToF> > TOFptr(new std::vector<brb::ToF>);
    std::unique_ptr<std::vector<brb::Cherenkov> > CHKVptr(new std::vector<brb::Cherenkov>);
    std::unique_ptr<std::vector<rawdata::RawBeamlineWC> > WCptr(new std::vector<rawdata::RawBeamlineWC>);
    
    if(evt.event()!=6964)
    {
    
    for (int i=0;i<particle_no;i++)
    {
        final_index = total_particles ;//+ i;
        //total_particles=final_index;
//        std::cout<<"putting particle number "<<i+1 <<" in the event"<<std::endl;
//        std::cout<<"totalparticles is "<<total_particles<<std::endl;
//        std::cout<<"final index is "<<final_index<<std::endl;
        // reformat into novasoft data types
        //Time of flight
        brb::ToF tofobject;
        std::pair <float,float> tofpair;
        if (vPresenttofus[final_index]==1)
        {
          tofpair.first=vttof_us[final_index];
        }
        if ( vPresenttofds[final_index]==1)
        {
          tofpair.second=vttof_ds[final_index];
        }
        if ( vPresenttofus[final_index]==1 && vPresenttofds[final_index]==1)
        {
            tofobject.SetTime(vttof_us[final_index]-vttof_ds[final_index]);
            tofobject.SetTimestamps(tofpair);
            TOFptr->push_back(tofobject);
            std::cout<<"the tof_us added to the event is "<< vttof_us[final_index] <<std::endl;
            std::cout<<"the tof_ds added to the event is "<< vttof_ds[final_index] <<std::endl;
            std::cout<<"the tof added to the event is "<< vttof_us[final_index]-vttof_ds[final_index] <<std::endl;
        }
        else
        {
            //std::cout<<"------------------------ didn't pass through both tof ------------------------------------------------------------------"<<std::endl;
        }
        //bool tofp=bool(vPresenttofus[final_index]) && bool(vPresenttofds[final_index]);
        //std::cout<<"tof present up =  "<< vPresenttofus[final_index] <<std::endl;
        //std::cout<<"tof present up =  "<< vPresenttofds[final_index] <<std::endl;
        
        //cherenkov
        
        brb::Cherenkov cherenkov_object;
        if(vPresentcherenkov[final_index]==1){
        cherenkov_object.SetTime(vtcherenkov[final_index]);
        CHKVptr->push_back(cherenkov_object);
        //std::cout<<"the cherenkov timestamp added to the event is "<< vtcherenkov[final_index] <<std::endl;
        }
        else
        {
            //std::cout<<"-------------------------- didn't pass through cherenkov-----------------------------------------------------------------"<<std::endl;
        }
        
        //std::cout<<"cherenkov present  =  "<< vPresentcherenkov[final_index] <<std::endl;
        
        //wirechambers
        
        std::vector<rawdata::RawBeamlineWC> WCobject;

        std::vector<double> fchannelx;
        std::vector<double> fchannely;
        std::vector<uint64_t> fWCtime;
        std::vector<unsigned int> fWCpresent;
        std::vector<double> xWC;
        std::vector<double> yWC;

        fchannelx.push_back(round((-(-491.51-vXWire_Chamber1[final_index])+64)/(TMath::Cos(16*TMath::Pi()/180))));
        fchannelx.push_back(round((-(-910.57-vXWire_Chamber2[final_index])+64)/(TMath::Cos(16*TMath::Pi()/180))));
        fchannelx.push_back(round(-(-1354.35-vXWire_Chamber3[final_index])+64));
        fchannelx.push_back(round(-(-1354.35-vXWire_Chamber4[final_index])+64));
        
        fchannely.push_back(round(vYWire_Chamber1[final_index]+64));
        fchannely.push_back(round(vYWire_Chamber2[final_index]+64));
        fchannely.push_back(round(vYWire_Chamber3[final_index]+64));
        fchannely.push_back(round(vYWire_Chamber4[final_index]+64));
        
        fWCtime.push_back(vTWire_Chamber1[final_index]);
        fWCtime.push_back(vTWire_Chamber2[final_index]);
        fWCtime.push_back(vTWire_Chamber3[final_index]);
        fWCtime.push_back(vTWire_Chamber4[final_index]);
        
        fWCpresent.push_back(vPresentWire_Chamber1[final_index]);
        fWCpresent.push_back(vPresentWire_Chamber2[final_index]);
        fWCpresent.push_back(vPresentWire_Chamber3[final_index]);
        fWCpresent.push_back(vPresentWire_Chamber4[final_index]);
        
        //just for debugging--------------------------------------------------
        xWC.push_back(vXWire_Chamber1[final_index]);
        xWC.push_back(vXWire_Chamber2[final_index]);
        xWC.push_back(vXWire_Chamber3[final_index]);
        xWC.push_back(vXWire_Chamber4[final_index]);
       
        yWC.push_back(vYWire_Chamber1[final_index]);
        yWC.push_back(vYWire_Chamber2[final_index]);
        yWC.push_back(vYWire_Chamber3[final_index]);
        yWC.push_back(vYWire_Chamber4[final_index]);
  //--mike helping me debug----------
//        for (unsigned int wcNum = 0; wcNum < 4; ++wcNum) {
//
//            rawdata::RawBeamlineWC::WCPulse xpulse(fchannelx[wcNum], fWCtime[wcNum]);
//            rawdata::RawBeamlineWC::WCPulse ypulse(fchannely[wcNum], fWCtime[wcNum]);
//
//            rawdata::RawBeamlineWC wc;//(wcNum, Form("WC%d", wcNum));
//            wc.SetWCNumber(wcNum+1);
//            //if (fWCpresent[wcNum]) {
//                wc.AddXPulse(xpulse);
//                wc.AddYPulse(ypulse);
//            //}
//
//            WCptr->push_back(wc);
//
//        }
  //------------------------------------
        
        for(unsigned int k=0;k<4;k++)
        {
          // Detector ID for this wire chamber
          beamlinegeo::DetectorID detectorID(beamlinegeo::SystemID(beamlinegeo::BeamlineSystem::WC), k);

            rawdata::RawBeamlineWC::WCPulse xpulse(fchannelx[k],fWCtime[k]);
            rawdata::RawBeamlineWC::WCPulse ypulse(fchannely[k],fWCtime[k]);
            rawdata::RawBeamlineWC WCobject(detectorID);
            WCobject.AddXPulse(xpulse);
            WCobject.AddYPulse(ypulse);
            
            if (fWCpresent[k]==1) {

                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;
                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;
                WCptr->push_back(WCobject);
            }
            else{
            std::cout<<"-------------------- No WC "<<k+1<<" objects added ---------------------------"<<std::endl;
            }
//            std::cout<<"the x wire number for WC "<< k+1 <<" added to the event is "<< fchannelx[k] <<std::endl;
//            std::cout<<"the x wire coordinate for WC "<<k+1<<" added to the event is "<<xWC[k] << std::endl;
//            std::cout<<"the y wire number for WC "<< k+1 <<" added to the event is "<< fchannely[k] <<std::endl;
//            std::cout<<"the y wire coordinate for WC "<<k+1<<" added to the event is "<<yWC[k] << std::endl;
//            std::cout<<"wire chamber  "<<k+1<<" present "<<fWCpresent[k] << std::endl;
            
        }
        
       
        fchannelx.clear();
        fchannely.clear();
        fWCtime.clear();
        fWCpresent.clear();
        xWC.clear();
        yWC.clear();
        total_particles++;
        
    }
    
    evt.put(std::move(TOFptr));
    evt.put(std::move(CHKVptr));
    evt.put(std::move(WCptr));
        
    }
    return;
    
}

void beamlinesim::BeamlineSimInput::reconfigure

(

const fhicl::ParameterSet &

pset

)

Definition at line 325 of file BeamlineSimInput_module.cc.

Referenced by BeamlineSimInput().

                                                                           {
}

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Referenced by art::RootOutput::endJob().

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

Member Data Documentation

int beamlinesim::BeamlineSimInput::event_no

Definition at line 88 of file BeamlineSimInput_module.cc.

Referenced by produce().

TFile* beamlinesim::BeamlineSimInput::fBeamlineSimFile

private

Definition at line 97 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

TTree* beamlinesim::BeamlineSimInput::fBeamlineSimTree

private

Definition at line 98 of file BeamlineSimInput_module.cc.

int beamlinesim::BeamlineSimInput::fBeamlineTreeIndex

private

Definition at line 102 of file BeamlineSimInput_module.cc.

int beamlinesim::BeamlineSimInput::final_index

Definition at line 89 of file BeamlineSimInput_module.cc.

Referenced by BeamlineSimInput(), and produce().

int beamlinesim::BeamlineSimInput::fPresentcherenkov

private

Definition at line 163 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresenttofds

private

Definition at line 162 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresenttofus

private

Definition at line 161 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresentWire_Chamber1

private

Definition at line 157 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresentWire_Chamber2

private

Definition at line 158 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresentWire_Chamber3

private

Definition at line 159 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::fPresentWire_Chamber4

private

Definition at line 160 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

art::InputTag beamlinesim::BeamlineSimInput::fRawToFDataLabel

private

Definition at line 94 of file BeamlineSimInput_module.cc.

float_t beamlinesim::BeamlineSimInput::ftcherenkov

private

Definition at line 194 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fttof_ds

private

Definition at line 193 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fttof_us

private

Definition at line 192 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fTWire_Chamber1

private

Definition at line 167 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fTWire_Chamber2

private

Definition at line 171 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fTWire_Chamber3

private

Definition at line 175 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fTWire_Chamber4

private

Definition at line 179 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fXWire_Chamber1

private

Definition at line 164 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fXWire_Chamber2

private

Definition at line 168 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fXWire_Chamber3

private

Definition at line 172 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fXWire_Chamber4

private

Definition at line 176 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fYWire_Chamber1

private

Definition at line 165 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fYWire_Chamber2

private

Definition at line 169 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fYWire_Chamber3

private

Definition at line 173 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fYWire_Chamber4

private

Definition at line 177 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fZWire_Chamber1

private

Definition at line 166 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fZWire_Chamber2

private

Definition at line 170 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fZWire_Chamber3

private

Definition at line 174 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

float_t beamlinesim::BeamlineSimInput::fZWire_Chamber4

private

Definition at line 178 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::particle_no

Definition at line 87 of file BeamlineSimInput_module.cc.

Referenced by produce().

TTreeReader beamlinesim::BeamlineSimInput::reader

Definition at line 85 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

int beamlinesim::BeamlineSimInput::total_particles

Definition at line 90 of file BeamlineSimInput_module.cc.

Referenced by BeamlineSimInput(), and produce().

TTree* beamlinesim::BeamlineSimInput::tree

Definition at line 86 of file BeamlineSimInput_module.cc.

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresentcherenkov

private

Definition at line 117 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresenttofds

private

Definition at line 116 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresenttofus

private

Definition at line 115 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresentWire_Chamber1

private

Definition at line 111 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresentWire_Chamber2

private

Definition at line 112 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresentWire_Chamber3

private

Definition at line 113 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<unsigned int> beamlinesim::BeamlineSimInput::vPresentWire_Chamber4

private

Definition at line 114 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vtcherenkov

private

Definition at line 148 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vtnova

private

Definition at line 152 of file BeamlineSimInput_module.cc.

std::vector<float_t> beamlinesim::BeamlineSimInput::vttof_ds

private

Definition at line 147 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vttof_us

private

Definition at line 146 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vTWire_Chamber1

private

Definition at line 121 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vTWire_Chamber2

private

Definition at line 125 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vTWire_Chamber3

private

Definition at line 129 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vTWire_Chamber4

private

Definition at line 133 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vxnova

private

Definition at line 149 of file BeamlineSimInput_module.cc.

std::vector<float_t> beamlinesim::BeamlineSimInput::vXWire_Chamber1

private

Definition at line 118 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vXWire_Chamber2

private

Definition at line 122 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vXWire_Chamber3

private

Definition at line 126 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vXWire_Chamber4

private

Definition at line 130 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vynova

private

Definition at line 150 of file BeamlineSimInput_module.cc.

std::vector<float_t> beamlinesim::BeamlineSimInput::vYWire_Chamber1

private

Definition at line 119 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vYWire_Chamber2

private

Definition at line 123 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vYWire_Chamber3

private

Definition at line 127 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vYWire_Chamber4

private

Definition at line 131 of file BeamlineSimInput_module.cc.

Referenced by beginJob(), and produce().

std::vector<float_t> beamlinesim::BeamlineSimInput::vznova

private

Definition at line 151 of file BeamlineSimInput_module.cc.

std::vector<float_t> beamlinesim::BeamlineSimInput::vZWire_Chamber1

private

Definition at line 120 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

std::vector<float_t> beamlinesim::BeamlineSimInput::vZWire_Chamber2

private

Definition at line 124 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

std::vector<float_t> beamlinesim::BeamlineSimInput::vZWire_Chamber3

private

Definition at line 128 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

std::vector<float_t> beamlinesim::BeamlineSimInput::vZWire_Chamber4

private

Definition at line 132 of file BeamlineSimInput_module.cc.

Referenced by beginJob().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/TestBeam/BeamlineSimInput_module.cc