protoana::protonbeamana Class Reference

Inheritance diagram for protoana::protonbeamana:

Public Member Functions
	protonbeamana (fhicl::ParameterSet const &p)
	protonbeamana (protonbeamana const &)=delete
	protonbeamana (protonbeamana &&)=delete
protonbeamana &	operator= (protonbeamana const &)=delete
protonbeamana &	operator= (protonbeamana &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
void	endJob () override
void	reset ()
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
const art::InputTag	fTrackModuleLabel
const art::InputTag	fBeamModuleLabel
TTree *	fTree
int	run
int	subrun
int	event
double	evttime
int	fNactivefembs [6]
bool	beam_inst_valid
int	beam_inst_trigger
int	beam_inst_nTracks
int	beam_inst_nMomenta
std::vector< double >	beam_inst_TOF
std::vector< int >	beam_inst_PDG_candidates
std::vector< int >	beam_inst_TOF_Chan
std::vector< double >	beamPosx
std::vector< double >	beamPosy
std::vector< double >	beamPosz
std::vector< double >	beamDirx
std::vector< double >	beamDiry
std::vector< double >	beamDirz
std::vector< double >	beamMomentum
double	tof
std::vector< double >	tofs
std::vector< int >	ch_tofs
short	low_pressure_status
short	high_pressure_status
double	low_pressure
double	high_pressure
protoana::ProtoDUNEBeamlineUtils	fBeamlineUtils
protoana::ProtoDUNEEmptyEventFinder	fEmptyEventFinder
protoana::ProtoDUNEDataUtils	dataUtil
bool	reco_reconstructable_beam_event
std::string	fCalorimetryTag
std::string	fTrackerTag
std::string	fHitTag
std::string	fShowerTag
std::string	fPFParticleTag
fhicl::ParameterSet	fBeamPars
bool	fVerbose
geo::GeometryCore const *	fGeometry
std::vector< double >	primtrk_startx
std::vector< double >	primtrk_starty
std::vector< double >	primtrk_startz
std::vector< double >	primtrk_endx
std::vector< double >	primtrk_endy
std::vector< double >	primtrk_endz
std::vector< double >	primtrk_Dirx
std::vector< double >	primtrk_Diry
std::vector< double >	primtrk_Dirz
std::vector< double >	primtrklen
std::vector< double >	primtrkID
std::vector< int >	primtrk_trktag
int	fisprimarytrack
int	fisprimaryshower
std::vector< double >	primtrk_dqdx
std::vector< double >	primtrk_resrange
std::vector< double >	primtrk_dedx
std::vector< double >	primtrk_range
std::vector< double >	primtrk_ke
std::vector< double >	primtrk_hitx
std::vector< double >	primtrk_hity
std::vector< double >	primtrk_hitz
std::vector< double >	primtrk_pitch
std::vector< int >	primtrk_wid
std::vector< double >	primtrk_pt
std::vector< size_t >	primtrk_calo_hit_index
std::vector< int >	primtrk_ch
double	cosine_beam_primtrk
std::vector< int >	pdg_code
std::vector< int >	n_daughter
std::vector< int >	n_beamparticle
std::vector< int >	isPrimary
std::vector< int >	pfp_self
std::vector< int >	pfp_daughter
std::vector< double >	Zintersection
std::vector< double >	Yintersection
std::vector< double >	Xintersection
std::vector< double >	timeintersection

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 96 of file protonbeamana_module.cc.

Constructor & Destructor Documentation

protoana::protonbeamana::protonbeamana ( fhicl::ParameterSet const &  p )

explicit

Definition at line 240 of file protonbeamana_module.cc.

        :
                EDAnalyzer(p),
                //fSpacePointModuleLabel(p.get< art::InputTag >("SpacePointModuleLabel")),
                fTrackModuleLabel(p.get< art::InputTag >("TrackModuleLabel")),
                fBeamModuleLabel(p.get< art::InputTag >("BeamModuleLabel")),
                fBeamlineUtils(p.get<fhicl::ParameterSet>("BeamlineUtils")),
                fEmptyEventFinder(p.get< fhicl::ParameterSet >("EmptyEventFinder")),

                //fBeamModuleLabel(p.get<std::string>("BeamModuleLabel")),

                //fTimeDecoderModuleLabel(p.get< art::InputTag >("TimeDecoderModuleLabel")),

                dataUtil(p.get<fhicl::ParameterSet>("DataUtils")),
                fCalorimetryTag(p.get<std::string>("CalorimetryTag")),
                fTrackerTag(p.get<std::string>("TrackerTag")),
                fHitTag(p.get<std::string>("HitTag")),
                fShowerTag(p.get<std::string>("ShowerTag")),
                fPFParticleTag(p.get<std::string>("PFParticleTag")),
                //fGeneratorTag(p.get<std::string>("GeneratorTag")),
                fBeamPars(p.get<fhicl::ParameterSet>("BeamPars")),
                //fUseCERNCalibSelection(p.get<bool>("UseCERNCalibSelection")),
                fVerbose(p.get<bool>("Verbose"))
{
        //if (fSaveTrackInfo == false) fSaveCaloInfo = false;

}

protoana::protonbeamana::protonbeamana ( protonbeamana const &  )

delete

protoana::protonbeamana::protonbeamana ( protonbeamana &&  )

delete

Member Function Documentation

void protoana::protonbeamana::analyze ( art::Event const &  e )

overridevirtual

Use the pandora metadata to tell us if this is a beam particle or not

Implements art::EDAnalyzer.

Definition at line 268 of file protonbeamana_module.cc.

{
        //reset containers
        protoana::protonbeamana::reset();  


        //Access the Beam Event ======================================================================//
        auto beamHandle = evt.getValidHandle<std::vector<beam::ProtoDUNEBeamEvent>>("beamevent");
        std::vector<art::Ptr<beam::ProtoDUNEBeamEvent>> beamVec;
        if( beamHandle.isValid()){
                art::fill_ptr_vector(beamVec, beamHandle);
        }
        const beam::ProtoDUNEBeamEvent & beamEvent = *(beamVec.at(0)); //Should just have one
        beam_inst_trigger = beamEvent.GetTimingTrigger();
        /////////////////////////////////////////////////////////////

        //Check the quality of the event
        std::cout << "Timing Trigger: " << beamEvent.GetTimingTrigger() << std::endl; 
        std::cout << "Is Matched: "     << beamEvent.CheckIsMatched() << std::endl;

        if( !fBeamlineUtils.IsGoodBeamlineTrigger( evt ) ){
                std::cout << "Failed beam quality check!\n" << std::endl;
                beam_inst_valid = false;
                return; //returns, does nothing when !=IsGoodBeamlineTrigger
        }

        std::cout << "Passed quality check!" << std::endl << std::endl;

        int nTracks = beamEvent.GetBeamTracks().size();
        /////////////////////////////////////////////////////////////

        if (evt.isRealData()) {
                std::vector< int > pdg_cands = fBeamlineUtils.GetPID( beamEvent, 1. );
                beam_inst_PDG_candidates.insert( beam_inst_PDG_candidates.end(), pdg_cands.begin(), pdg_cands.end() ); //save info of incoming particle candidates
        }


        //Access momentum
        const std::vector< double > & momenta = beamEvent.GetRecoBeamMomenta();
        int nMomenta = momenta.size();
        std::cout << "Number of reconstructed momenta: " << nMomenta << std::endl;

        beam_inst_nTracks = nTracks;
        beam_inst_nMomenta = nMomenta;

        //if( momenta.size() > 0 ) 
                //std::cout << "Measured Momentum: " << momenta.at(0) << std::endl;
        ///////////////////////////////////////////////////////////// 

        //Access time of flight
        const std::vector< double > & the_tofs  = beamEvent.GetTOFs();
        const std::vector< int    > & the_chans = beamEvent.GetTOFChans();

        std::cout<<"run/subrun/event:"<<evt.run()<<"/"<<evt.subRun()<<"/"<<evt.id().event()<<std::endl; 
        std::cout << "Number of measured TOF: " << the_tofs.size() << std::endl;
        std::cout << "First TOF: "              << beamEvent.GetTOF()         << std::endl;
        std::cout << "First TOF Channel: "      << beamEvent.GetTOFChan()     << std::endl << std::endl;

        std::cout << "All (TOF, Channels): " << std::endl;
        for( size_t i = 0; i < the_tofs.size(); ++i ){
                std::cout << "\t(" << the_tofs[i] << ", " << the_chans[i] << ")" << std::endl;
                beam_inst_TOF.push_back(the_tofs[i]);
                beam_inst_TOF_Chan.push_back(the_chans[i]);
        }
        std::cout << std::endl;
        /////////////////////////////////////////////////////////////

        //Access Cerenkov info
        std::cout << "Cerenkov status, pressure:" << std::endl;
        std::cout << "C0: " << beamEvent.GetCKov0Status() << ", " << beamEvent.GetCKov0Pressure() << std::endl;
        std::cout << "C1: " << beamEvent.GetCKov1Status() << ", " << beamEvent.GetCKov1Pressure() << std::endl << std::endl;

        //if (beamEvent.GetBITrigger() == 1) {
        //ckov0status = beamEvent.GetCKov0Status(); //low_pressure_status
        //ckov1status = beamEvent.GetCKov1Status(); //high_pressure_status

        //ckov0pressure=beamEvent.GetCKov0Pressure();
        //ckov1pressure=beamEvent.GetCKov1Pressure();
        //}

        if(beamEvent.GetCKov0Pressure() < beamEvent.GetCKov1Pressure() ){
                high_pressure_status = beamEvent.GetCKov1Status();
                low_pressure_status = beamEvent.GetCKov0Status();

                high_pressure=beamEvent.GetCKov1Pressure();
                low_pressure=beamEvent.GetCKov0Pressure(); 
        }
        else{
                high_pressure_status = beamEvent.GetCKov0Status();
                low_pressure_status = beamEvent.GetCKov1Status();

                high_pressure=beamEvent.GetCKov0Pressure();
                low_pressure=beamEvent.GetCKov1Pressure(); 
        }

        ///////////////////////////////////////////////////////////// 

        double nom_beam_mon=fBeamPars.get<double>("Momentum");
        std::cout<<"Selected nominal beam momentum is: "<<nom_beam_mon<<" GeV/c"<<std::endl;

        //Access PID
        //std::vector< int > pids = fBeamlineUtils.GetPID( beamEvent, nom_beam_mon);
        //for( size_t i = 0; i < pids.size(); ++i ){
        //std::cout << "pid["<<i<<"]: "<< pids[i] << std::endl;
        //}

        PossibleParticleCands candidates = fBeamlineUtils.GetPIDCandidates( beamEvent, nom_beam_mon);
        std::cout << "electron " << candidates.electron << std::endl;
        std::cout << "muon "     << candidates.muon     << std::endl;
        std::cout << "pion "     << candidates.pion     << std::endl;
        std::cout << "kaon "     << candidates.kaon     << std::endl;
        std::cout << "proton "   << candidates.proton   << std::endl;

        std::cout << std::endl;
        //if (candidates.proton==0) {
        //std::cout<<"no proton!!!!!!"<<std::endl;
        //}
        /////////////////////////////////////////////////////////////

        //Manual Event Selection -----------------------------------------------------------------------------------------------------------------------------------------------//
        /*bool IsProton=false;
        //old Tof: (!fUseCERNCalibSelection&&tof>170.)
        if (nom_beam_mon==1.) { 
        if ((beamEvent.GetTOF()>110.&&beamEvent.GetTOF()<160.)&&low_pressure_status==0) {
        IsProton=true;
        std::cout<<"-->> This is a proton candidate..."<<std::endl;
        }
        }
        */

        //HY::For Carlo info
        std::vector<art::Ptr<recob::Track> > tracklist;
        auto trackListHandle = evt.getHandle< std::vector<recob::Track> >(fTrackModuleLabel);
        if (trackListHandle) art::fill_ptr_vector(tracklist, trackListHandle);
        else return;
        art::FindManyP<anab::Calorimetry> fmcal(trackListHandle, evt, fCalorimetryTag);
        art::FindManyP<recob::PFParticle> pfp_trk_assn(trackListHandle, evt, "pandoraTrack");
        auto allHits = evt.getValidHandle<std::vector<recob::Hit> >(fHitTag);

        //HY::Cluster info
        //art::Handle< std::vector<recob::PFParticle> > PFPListHandle;
        //std::vector<art::Ptr<recob::PFParticle> > pfplist;
        //if(evt.getByLabel("pandoraTrack",trackListHandle)) art::fill_ptr_vector(tracklist, trackListHandle);
        //if(evt.getByLabel("pandora",PFPListHandle)) art::fill_ptr_vector(pfplist, PFPListHandle);
        
        auto trackListHandle2 = evt.getHandle< std::vector<recob::Track> >("pandoraTrack");
        if (trackListHandle2) art::fill_ptr_vector(tracklist, trackListHandle2);

        std::vector<art::Ptr<recob::PFParticle> > pfplist;
        auto PFPListHandle = evt.getHandle< std::vector<recob::PFParticle> > ("pandora");
        if (PFPListHandle) art::fill_ptr_vector(pfplist, PFPListHandle);

        // to get information about the hits
        //art::Handle< std::vector<recob::Cluster> > clusterListHandle; // to get information about the hits
        //std::vector<art::Ptr<recob::Cluster>> clusterlist;
        //if(evt.getByLabel("pandora", clusterListHandle))
                //art::fill_ptr_vector(clusterlist, clusterListHandle);

        std::vector<art::Ptr<recob::Cluster>> clusterlist;
        auto clusterListHandle = evt.getHandle< std::vector<recob::Cluster> > ("pandora");
        if (clusterListHandle) art::fill_ptr_vector(clusterlist, clusterListHandle);

        art::FindManyP<recob::Cluster> fmcp(PFPListHandle,evt,"pandora");
        art::FindManyP<recob::Track> pftrack(PFPListHandle,evt,"pandoraTrack");
        std::cout<<"number of pfp_particles "<<pfplist.size()<<std::endl;
        //std::cout<<" size of pfParticles "<<pfParticles.size()<<std::endl;

        //HY::call geometry service
        fGeometry = &*(art::ServiceHandle<geo::Geometry>());



        // Implementation of required member function here.
        run = evt.run();
        subrun = evt.subRun();
        event = evt.id().event();
        art::Timestamp ts = evt.time();

         // Is this a reconstructable beam event?
        reco_reconstructable_beam_event = !fEmptyEventFinder.IsEmptyEvent(evt);
        std::cout<<"reco_reconstructable_beam_event:"<<reco_reconstructable_beam_event<<std::endl;

        //std::cout<<ts.timeHigh()<<" "<<ts.timeLow()<<std::endl;
        if (ts.timeHigh() == 0){
                //TTimeStamp tts(ts.timeLow());
                //evttime = tts.AsDouble();
                TTimeStamp ts2(ts.timeLow());
                evttime = ts2.AsDouble();
        }
        else{
                //TTimeStamp tts(ts.timeHigh(), ts.timeLow());
                //evttime = tts.AsDouble();
                TTimeStamp ts2(ts.timeHigh(), ts.timeLow());
                evttime = ts2.AsDouble();
        }

        // Get number of active fembs
        if(!evt.isRealData()){
                for(int ik=0; ik<6; ik++)
                        fNactivefembs[ik]=20;
        }
        else{
                for(int ik=0; ik<6; ik++)
                        fNactivefembs[ik]=dataUtil.GetNActiveFembsForAPA(evt,ik);
        }


        if (beamVec.size()&&candidates.proton==1) { //if beam proton
                if (beamEvent.GetTimingTrigger()==12) { //get beam timing trigger
                        if (beamEvent.CheckIsMatched()){ //if CheckIsMatched
                                //Get TOF info
                                if (beamEvent.GetTOFChan() != -1) {//if TOFChan == -1, then there was not a successful match, if it's 0, 1, 2, or 3, then there was a good match corresponding to the different pair-wise combinations of the upstream and downstream channels
                                        tof = beamEvent.GetTOF();
                                        for( size_t ii = 0; ii < the_tofs.size(); ++ii ){
                                                tofs.push_back(the_tofs[ii]);
                                                ch_tofs.push_back(the_chans[ii]);
                                        } 
                                }
                                //Get beam particle trajectory info
                                //auto & tracks = beaminfo[0]->GetBeamTracks();
                                auto & tracks = beamEvent.GetBeamTracks();
                                std::cout<<"ToF:"<<tof<<" [ns]"<<std::endl;
                                std::cout<<"beam trk size:"<<tracks.size()<<std::endl;
                                for (size_t i = 0; i<tracks.size(); ++i){
                                        beamPosx.push_back(tracks[i].End().X());
                                        beamPosy.push_back(tracks[i].End().Y());
                                        beamPosz.push_back(tracks[i].End().Z());
                                        beamDirx.push_back(tracks[i].StartDirection().X());
                                        beamDiry.push_back(tracks[i].StartDirection().Y());
                                        beamDirz.push_back(tracks[i].StartDirection().Z());

                                        std::cout<<"run/subrun/evt:"<<run<<"/"<<subrun<<"/"<<event<<std::endl;  
                                        std::cout<<"beamPosx/beamPosy/beamPosz:"<<tracks[i].End().X()<<"/"<<tracks[i].End().Y()<<"/"<<tracks[i].End().Z()<<std::endl;
                                        std::cout<<"beamDirx/beamDiry/beamDirz:"<<tracks[i].StartDirection().X()<<"/"<<tracks[i].StartDirection().Y()<<"/"<<tracks[i].StartDirection().Z()<<std::endl;
                                        //std::cout<<"beamDirx^2+beamDiry^2+beamDirz^2:"<<tracks[i].StartDirection().X()*tracks[i].StartDirection().X()+tracks[i].StartDirection().Y()*tracks[i].StartDirection().Y()+tracks[i].StartDirection().Z()*tracks[i].StartDirection().Z()<<std::endl;
                                }
                                //Get reconstructed beam momentum info
                                //auto & beammom = beaminfo[0]->GetRecoBeamMomenta();
                                //auto & beammom = beamEvent.GetRecoBeamMomenta();
                                //auto & beammom = beamEvent.GetRecoBeamMomenta();
                                //std::cout<<"==============================================================="<<std::endl;
                                std::cout<<"beam mom size:"<<momenta.size()<<std::endl;
                                for (size_t i = 0; i<momenta.size(); ++i){
                                        beamMomentum.push_back(momenta[i]);
                                        std::cout<<"beam mom["<<i<<"]:"<<momenta[i]<<" [GeV]"<<std::endl;
                                }
                                //std::cout<<"==============================================================="<<std::endl;

                                //put the track parameters in the cryo here ----------------------------------------------------------------------------//
                                /*
                                // Now we want to access the output from Pandora. This comes in the form of particle flow objects (recob::PFParticle).
                                // The primary PFParticles are those we want to consider and these PFParticles then have a hierarchy of daughters that
                                // describe the whole interaction of a given primary particle
                                //
                                //                     / daughter track
                                //                    /
                                //  primary track    /   
                                //  ---------------- ---- daughter track
                                //                   \
                                //                   /\-
                                //                   /\\-- daughter shower
                                //
                                // The above primary PFParticle will have links to three daughter particles, two track-like and one shower-like
                                */

                                std::cout<<"\n*******************************************************"<<std::endl;
                                std::cout<<"Moving on to the PFParticle section..."<<std::endl; 

                                // Get the PFParticle utility
                                protoana::ProtoDUNEPFParticleUtils pfpUtil;

                                // Get the track utility
                                protoana::ProtoDUNETrackUtils trackUtil;

                                // Get all of the PFParticles, by default from the "pandora" product
                                auto recoParticles = evt.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleTag);

                                // We'd like to find the beam particle. Pandora tries to do this for us, so let's use the PFParticle utility 
                                // to look for it. Pandora reconstructs slices containing one (or sometimes more) primary PFParticles. These    
                                // are tagged as either beam or cosmic for ProtoDUNE. This function automatically considers only those 
                                // PFParticles considered as primary

                                /// Use the pandora metadata to tell us if this is a beam particle or not
                                //bool isBeamParticle=dataUtil.IsBeamParticle(fPFParticleTag, evt, );
                                //bool IsBeamParticle(const recob::PFParticle &particle, art::Event const &evt, const std::string particleLabel) const;


                                //std::vector<const recob::PFParticle*> beamParticles = pfpUtil.GetPFParticlesFromBeamSlice(evt,fPFParticleTag);
                                auto beamParticles = pfpUtil.GetPFParticlesFromBeamSlice(evt,fPFParticleTag);
                                //HY:: From the new version, prepend recob::PFParticle* with "const"
                                if(beamParticles.size() == 0){
                                        std::cout << "We found no beam particles for this event... moving on" << std::endl;
                                        //return;
                                }

                                //unsigned int nPrimPFPartices=pfpUtil.GetNumberPrimaryPFParticle(evt,fPFParticleTag); //count all the particles
                                //std::cout<<"--> Number of all prim. PFPartices:"<<nPrimPFPartices<<std::endl;

                                //int tmp_counter=0;
                                n_beamparticle.push_back(beamParticles.size());
                                std::cout<<"we have "<<beamParticles.size()<<" beam particle(s)"<<std::endl;
                                for(const recob::PFParticle* particle : beamParticles){
                                        int nTrack=0;
                                        int nShower=0;

                                        // "particle" is the pointer to our beam particle. The recob::Track or recob::Shower object
                                        // of this particle might be more helpful. These return null pointers if not track-like / shower-like
                                        const recob::Track* thisTrack = pfpUtil.GetPFParticleTrack(*particle,evt,fPFParticleTag,fTrackerTag);
                                        const recob::Shower* thisShower = pfpUtil.GetPFParticleShower(*particle,evt,fPFParticleTag,fShowerTag);
                                        if(thisTrack != 0x0) { //thisTrack
                                                std::cout << "Beam particle is track-like" << std::endl;

                                                fisprimarytrack               = 1;
                                                fisprimaryshower              = 0;

                                                nTrack++;
                                                primtrk_trktag.push_back(1);

                                                std::vector<anab::Calorimetry> calovector = trackUtil.GetRecoTrackCalorimetry(*thisTrack, evt, fTrackerTag, fCalorimetryTag);
                                                for (auto & calo : calovector){
                                                        if (calo.PlaneID().Plane == 2){ //only collection plane
                                                                //std::cout<<"ck0"<<std::endl;
                                                                primtrk_range.push_back(calo.Range());
                                                                //std::cout<<"ck1"<<std::endl;
                                                                primtrk_ke.push_back(calo.KineticEnergy());
                                                                //std::cout<<"ck2"<<std::endl;
                                                                //std::cout<<"calo.Range():"<<calo.Range()<<std::endl;
                                                                //std::cout<<"calo.KineticEnergy():"<<calo.KineticEnergy()<<std::endl;
                                                                //std::cout<<"calo.dQdx().size():"<<calo.dQdx().size()<<std::endl;
                                                                for (size_t ihit = 0; ihit < calo.dQdx().size(); ++ihit){ //loop over hits
                                                                        //std::cout<<"ck3--0"<<std::endl;
                                                                        //std::cout<<"calo.dQdx["<<ihit<<"]:"<<calo.dQdx()[ihit]<<std::endl;
                                                                        //std::cout<<"calo.dEdx["<<ihit<<"]:"<<calo.dEdx()[ihit]<<std::endl;
                                                                        //std::cout<<"calo.DeadWireResRC["<<ihit<<"]:"<<calo.DeadWireResRC()[ihit]<<std::endl;
                                                                        primtrk_dqdx.push_back(calo.dQdx()[ihit]);
                                                                        primtrk_resrange.push_back(calo.ResidualRange()[ihit]);
                                                                        //primtrk_deadwireresrange.push_back(calo.DeadWireResRC()[ihit]);
                                                                        primtrk_dedx.push_back(calo.dEdx()[ihit]);
                                                                        primtrk_pitch.push_back(calo.TrkPitchVec()[ihit]);
                                                                        //std::cout<<"ck3--1"<<std::endl;

                                                                        const auto &primtrk_pos=(calo.XYZ())[ihit];
                                                                        primtrk_hitx.push_back(primtrk_pos.X());
                                                                        primtrk_hity.push_back(primtrk_pos.Y());
                                                                        primtrk_hitz.push_back(primtrk_pos.Z());
                                                                        //std::cout<<"ck3--2"<<std::endl;

                                                                        primtrk_calo_hit_index.push_back(calo.TpIndices()[ihit]);
                                                                        const recob::Hit & theHit = (*allHits)[ calo.TpIndices()[ihit] ];
                                                                        primtrk_wid.push_back(theHit.WireID().Wire);
                                                                        primtrk_pt.push_back(theHit.PeakTime());
                                                                        primtrk_ch.push_back(theHit.Channel()); 
                                                                        //std::cout<<"ck3--3"<<std::endl;
                                                                        //std::cout<<theHit.WireID().Wire<<std::endl;

                                                                        //double dedx_range=17.*pow(calo.ResidualRange()[ihit],-0.42);
                                                                        //if (dedx_range>=17.79&&dedx_range<=18.7&&thisTrack->ID()==38) {
                                                                        //if (dedx_range>=17.79&&dedx_range<=18.7) {
                                                                        //std::cout<<"dqdx="<<calo.dQdx()[ihit]<<"; resrange="<<calo.ResidualRange()[ihit]<<"; pitch="<<calo.TrkPitchVec()[ihit]<<std::endl;
                                                                        //std::cout<<"(X,Y,Z)="<<"("<<calo.XYZ()[ihit].X()<<","<<calo.XYZ()[ihit].Y()<<","<<calo.XYZ()[ihit].Z()<<")"<<std::endl;
                                                                        //std::cout<<"wid/pt/ch/triID:"<<theHit.WireID().Wire<<"/"<<theHit.PeakTime()<<"/"<<theHit.Channel()<<"/"<<thisTrack->ID()<<std::endl;
                                                                        //std::cout<<"wid/pt/ch/tpc:"<<theHit.WireID().Wire<<"/"<<theHit.PeakTime()<<"/"<<theHit.Channel()<<"/"<<theHit.WireID().TPC<<std::endl;
                                                                        //}
                                                                        //std::cout<<"(X,Y,Z)="<<"("<<primtrk_pos.X()<<","<<primtrk_pos.Y()<<","<<primtrk_pos.Z()<<")"<<std::endl;
                                                                        //std::cout<<"(X,Y,Z)="<<"("<<tmp_primtrk_hitx[ihit]<<","<<tmp_primtrk_hity[ihit]<<","<<tmp_primtrk_hitz[ihit]<<")"<<std::endl;
                                                                } //loop over hits
                                                                //std::cout<<"ck3"<<std::endl;
                                                                } //only collection plane
                                                        }

                                                        //HY::Here comes the start/end position of primary track
                                                        // Find the particle vertex. We need the tracker tag here because we need to do a bit of
                                                        // additional work if the PFParticle is track-like to find the vertex. 
                                                        const TVector3 vtx = pfpUtil.GetPFParticleVertex(*particle,evt,fPFParticleTag,fTrackerTag);
                                                        const TVector3 vtx_end(thisTrack->Trajectory().End().X(), thisTrack->Trajectory().End().Y(), thisTrack->Trajectory().End().Z());
                                                        primtrk_startx.push_back(vtx.X());
                                                        primtrk_starty.push_back(vtx.Y());
                                                        primtrk_startz.push_back(vtx.Z());
                                                        std::cout<<"vtx_X:"<<vtx.X()<<" ; vtx_Y:"<<vtx.Y()<<" ; vtx_Z:"<<vtx.Z()<<std::endl;

                                                        primtrk_endx.push_back(vtx_end.X());
                                                        primtrk_endy.push_back(vtx_end.Y());
                                                        primtrk_endz.push_back(vtx_end.Z());
                                                        std::cout<<"vtx_end_X:"<<vtx_end.X()<<" ; vtx_end_Y:"<<vtx_end.Y()<<" ; vtx_end_Z:"<<vtx_end.Z()<<std::endl;

                                                        if (vtx.Z()==vtx_end.Z()) { //warning message if Z_start=Z_end
                                                                std::cout<<"WARNING!! StartZ and EndZ are the same!!"<<std::endl;       
                                                        }

                                                        //Use Ajib's intersection calculation function
                                                        //[1]identify the beam track and tag other tracks
                                                        std::vector<float> Stw, Endw, Stt, Endt, Stwires, Endwires, Stticks, Endticks, TPCb, TPCcl;
                                                        Stw.clear(); Endw.clear(); Stt.clear();  Endt.clear();  Stwires.clear();  Endwires.clear(); Stticks.clear(); Endticks.clear(); TPCb.clear(); TPCcl.clear();
                                                        float den;
                                                        float numw, numt,wire_no,ticks_no;

                                                        for(size_t p1=0;p1<pfplist.size();p1++){
                                                                std::vector<art::Ptr<recob::Track>> trk=pftrack.at(p1);
                                                                std::vector<art::Ptr<recob::Cluster>> allClusters=fmcp.at(p1);
                                                                //std::cout<<"fprimaryID:"<<fprimaryID<<std::endl;
                                                                for(size_t c1=0;c1<allClusters.size();c1++){
                                                                        if(allClusters[c1]->Plane().Plane!=2) continue;
                                                                        if(trk.size() && int(trk[0].key())==(int)(thisTrack->ID())){
                                                                                Stw.push_back(allClusters[c1]->StartWire());
                                                                                Endw.push_back(allClusters[c1]->EndWire());
                                                                                Stt.push_back(allClusters[c1]->StartTick());
                                                                                Endt.push_back(allClusters[c1]->EndTick());
                                                                                TPCb.push_back(allClusters[c1]->Plane().TPC);
                                                                                //std::cout<<"\nst/end wire:"<<allClusters[c1]->StartWire()<<"/"<<allClusters[c1]->EndWire()<<std::endl;
                                                                        }
                                                                        else{
                                                                                Stwires.push_back(allClusters[c1]->StartWire());
                                                                                Endwires.push_back(allClusters[c1]->EndWire());
                                                                                Stticks.push_back(allClusters[c1]->StartTick());
                                                                                Endticks.push_back(allClusters[c1]->EndTick());
                                                                                TPCcl.push_back(allClusters[c1]->Plane().TPC);
                                                                        }
                                                                }
                                                        }
                                                        //[2]find interaction points if any (assuming all tracks are straight, find the interaction points)
                                                        for(size_t clt=0;clt<Stw.size();clt++){
                                                                for(size_t cl1=0;cl1<Stwires.size();cl1++){
                                                                        if(TPCcl[cl1]!=TPCb[clt]) continue;
                                                                        //std::cout<<"tpc are equal "<<std::endl;
                                                                        den=(Stw[clt]-Endw[clt])*(Stticks[cl1]-Endticks[cl1])-(Stt[clt]-Endt[clt])*(Stwires[cl1]-Endwires[cl1]);
                                                                        if(den==0) continue;
                                                                        numw=(Stw[clt]*Endt[clt]-Stt[clt]*Endw[clt])*(Stwires[cl1]-Endwires[cl1])-(Stw[clt]-Endw[clt])*(Stwires[cl1]*Endticks[cl1]-Stticks[cl1]*Endwires[cl1]);
                                                                        numt=(Stw[clt]*Endt[clt]-Stt[clt]*Endw[clt])*(Stticks[cl1]-Endticks[cl1])-(Stt[clt]-Endt[clt])*(Stwires[cl1]*Endticks[cl1]-Stticks[cl1]*Endwires[cl1]);
                                                                        wire_no=numw/den;
                                                                        //  std::cout<<"wireno and ticks not solution "<<wire_no<<"  "<<ticks_no<<std::endl;
                                                                        ticks_no=numt/den;
                                                                        if(((Stw[clt]<wire_no && Endw[clt]>wire_no)||(Stw[clt]>wire_no && Endw[clt]<wire_no))&&((Stt[clt]<ticks_no && Endt[clt]>ticks_no)||(Stt[clt]>ticks_no && Endt[clt]<ticks_no)) && ((Stwires[cl1]<wire_no && Endwires[cl1]>wire_no)||(Stwires[cl1]>wire_no && Endwires[cl1]<wire_no)) && ((Stticks[cl1]<ticks_no && Endticks[cl1]>ticks_no)||(Stticks[cl1]>ticks_no && Endticks[cl1]<ticks_no))) {
                                                                                std::cout<<"intersection wire and ticks are "<<std::round(wire_no)<<"  "<<ticks_no<<" Stw Endw StT EndT "<<Stwires[cl1]<<" "<<Endwires[cl1]<<" "<<Stticks[cl1]<<" "<<Endticks[cl1]<<std::endl;
                                                                                double xyzStart[3];
                                                                                double xyzEnd[3];
                                                                                unsigned int wireno=std::round(wire_no);
                                                                                geo::WireID wireid(0,TPCb[clt],2,wireno);
                                                                                fGeometry->WireEndPoints(0,TPCb[clt],2,wireno, xyzStart, xyzEnd);
                                                                                //fGeometry->WireEndPoints(wireid, xyzStart, xyzEnd);
                                                                                std::cout<<"Z position of intersection = "<<xyzStart[2]<<" "<<xyzEnd[2]<<"  "<<wireno<<std::endl;
                                                                                Zintersection.push_back(xyzStart[2]);
                                                                                Yintersection.push_back(xyzStart[1]);
                                                                                Xintersection.push_back(xyzStart[0]);
                                                                                timeintersection.push_back(ticks_no);
                                                                        }
                                                                }
                                                        }
                                                } //thisTrack
                                                if(thisShower != 0x0) { //thisShower 
                                                        std::cout << "Beam particle is shower-like" << std::endl;
                                                        fisprimarytrack               = 0;
                                                        fisprimaryshower              = 1;

                                                        nShower++;
                                                        primtrk_trktag.push_back(-1);
                                                } //thisShower

                                                //HY Add
                                                pdg_code.push_back(particle->PdgCode());
                                                n_daughter.push_back(particle->NumDaughters());
                                                isPrimary.push_back(particle->IsPrimary());
                                                pfp_self.push_back(particle->Self());
                                                //pfp_parent.push_back(particle->Parent());

                                                std::cout<<"pdg code:"<<particle->PdgCode()<<std::endl;
                                                std::cout<<"IsPrimary:"<<particle->IsPrimary()<<std::endl;
                                                std::cout<<"NumDaughters:"<<particle->NumDaughters()<<std::endl;
                                                std::cout<<"Self:"<<particle->Self()<<std::endl;        
                                                std::cout<<"Parent:"<<particle->Parent()<<std::endl;

                                                if ((particle->NumDaughters())>0) {
                                                        for (int ii=0; ii<(particle->NumDaughters());++ii) {
                                                                std::cout<<"Daughter["<<ii<<"]:"<<particle->Daughter(ii)<<std::endl;
                                                                pfp_daughter.push_back(particle->Daughter(ii));
                                                        }
                                                }
                                                else {
                                                        pfp_daughter.push_back(-99);
                                                }

                                                //int pdg_code=pfpUtil.PdgCode(*particle,evt,fPFParticleTag,fShowerTag);
                                                //std::cout<<"IsDaughters:"<<particle->Daughters()<<std::endl;  
                                                //HY::Add parameters for the primary tracks here ---------------//

                                                // Get track direction
                                                if (thisTrack) {
                                                        //pdg_code=thisTrack->PdgCode(); 
                                                        //std::cout<<"pdg code:"<<pdg_code<<std::endl;

                                                        auto trackdir = thisTrack->StartDirection();
                                                        std::cout<<"run/subrun/event:"<<run<<"/"<<subrun<<"/"<<event<<std::endl;        
                                                        std::cout<<"trkDirx/trkDiry/trkDirz:"<<trackdir.X()<<"/"<<trackdir.Y()<<"/"<<trackdir.Z()<<std::endl;
                                                        primtrk_Dirx.push_back(trackdir.X());
                                                        primtrk_Diry.push_back(trackdir.Y());
                                                        primtrk_Dirz.push_back(trackdir.Z());

                                                        primtrklen.push_back(thisTrack->Length()); //track length
                                                        std::cout<<"trk length: "<<thisTrack->Length()<<" [cm]"<<std::endl;
                                                        primtrkID.push_back(thisTrack->ID());
                                                        std::cout<<"trk ID: "<<thisTrack->ID()<<""<<std::endl; //HY::Fix me::trk ID seems wrong 
                                                        //std::cout<<"TrackId: "<<thisTrack->TrackId()<<std::endl;

                                                        //std::cout<<"trkDirx^2+trkDiry^2+trkDirz^2:"<<trackdir.X()*trackdir.X()+trackdir.Y()*trackdir.Y()+trackdir.Z()*trackdir.Z()<<std::endl;
                                                        //int nn=tracks.size()-1;
                                                        //cosine_beam_primtrk=tracks[nn].StartDirection().X()*trackdir.X()+tracks[nn].StartDirection().Y()*trackdir.Y()+tracks[nn].StartDirection().Z()*trackdir.Z();
                                                        if (tracks.size()){
                                                                cosine_beam_primtrk=tracks[0].StartDirection().X()*trackdir.X()+tracks[0].StartDirection().Y()*trackdir.Y()+tracks[0].StartDirection().Z()*trackdir.Z();
                                                        }
                                                        // fill a histogram of trackdir.X()*beamdir.X() + .....
                                                        // try to get calorimetry info of this track

                                                }
                                                // Now we can look for the interaction point of the particle if one exists, i.e where the particle
                                                // scatters off an argon nucleus. Shower-like objects won't have an interaction point, so we can
                                                // check this by making sure we get a sensible position
                                                const TVector3 interactionVtx = pfpUtil.GetPFParticleSecondaryVertex(*particle,evt,fPFParticleTag,fTrackerTag);


                                                // Let's get the daughter PFParticles... we can do this simply without the utility
                                                for(const int daughterID : particle->Daughters()){
                                                        // Daughter ID is the element of the original recoParticle vector
                                                        const recob::PFParticle *daughterParticle = &(recoParticles->at(daughterID));
                                                        std::cout << "Daughter " << daughterID << " has " << daughterParticle->NumDaughters() << " daughters" << std::endl;
                                                }

                                                // For actually studying the objects, it is easier to have the daughters in their track and shower forms.
                                                // We can use the utility to get a vector of track-like and a vector of shower-like daughters
                                                const std::vector<const recob::Track*> trackDaughters = pfpUtil.GetPFParticleDaughterTracks(*particle,evt,fPFParticleTag,fTrackerTag);  
                                                const std::vector<const recob::Shower*> showerDaughters = pfpUtil.GetPFParticleDaughterShowers(*particle,evt,fPFParticleTag,fShowerTag);  
                                                std::cout << "Beam particle has " << trackDaughters.size() << " track-like daughters and " << showerDaughters.size() << " shower-like daughters." << std::endl;

                                                std::cout<<"# total Tracks:"<<nTrack<<std::endl;
                                                std::cout<<"# total Showers:"<<nShower<<std::endl;
                                                //tmp_counter++;
                                        }
                                        //std::cout<<"tmp_counter:"<<tmp_counter<<"\n\n\n"<<std::endl;
                                        //'tmp_counter' should be the same as 'nBeamP'
                                        std::cout<<"*******************************************************"<<std::endl;

                                        //put the track parameters in the cryo here ----------------------------------------------------------------------------//


                                } //if CheckIsMatched
                        } //get beam timing trigger

                } //if beam proton

                fTree->Fill();
        }

void protoana::protonbeamana::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 818 of file protonbeamana_module.cc.

        {
                art::ServiceHandle<art::TFileService> tfs;
                fTree = tfs->make<TTree>("beamana","beam analysis tree");
                fTree->Branch("run",&run,"run/I");
                fTree->Branch("subrun",&subrun,"subrun/I");
                fTree->Branch("event",&event,"event/I");
                //fTree->Branch("trigger",&trigger,"trigger/I");
                fTree->Branch("evttime",&evttime,"evttime/D");
                fTree->Branch("reco_reconstructable_beam_event",&reco_reconstructable_beam_event);

                fTree->Branch("Nactivefembs",&fNactivefembs,"Nactivefembs[5]/I");

                fTree->Branch("beam_inst_valid", &beam_inst_valid);
                fTree->Branch("beam_inst_trigger", &beam_inst_trigger);
                fTree->Branch("beam_inst_nTracks", &beam_inst_nTracks);
                fTree->Branch("beam_inst_nMomenta", &beam_inst_nMomenta);

                fTree->Branch("beam_inst_PDG_candidates", &beam_inst_PDG_candidates);

                fTree->Branch("beam_inst_TOF", &beam_inst_TOF);
                fTree->Branch("beam_inst_TOF_Chan", &beam_inst_TOF_Chan);

                fTree->Branch("isprimarytrack", &fisprimarytrack, "isprimarytrack/I");
                fTree->Branch("isprimaryshower", &fisprimaryshower, "isprimaryshower/I");

                fTree->Branch("beamPosx",&beamPosx);
                fTree->Branch("beamPosy",&beamPosy);
                fTree->Branch("beamPosz",&beamPosz);
                fTree->Branch("beamDirx",&beamDirx);
                fTree->Branch("beamDiry",&beamDiry);
                fTree->Branch("beamDirz",&beamDirz);
                fTree->Branch("beamMomentum",&beamMomentum);
                fTree->Branch("tof", &tof, "tof/D");
                fTree->Branch("tofs",&tofs);
                fTree->Branch("ch_tofs",&ch_tofs);
                fTree->Branch("low_pressure_status", &low_pressure_status, "low_pressure_status/S");
                fTree->Branch("high_pressure_status", &high_pressure_status, "high_pressure_status/S");
                fTree->Branch("low_pressure", &low_pressure, "low_pressure/D");
                fTree->Branch("high_pressure", &high_pressure, "high_pressure/D");

                fTree->Branch("cosine_beam_primtrk", &cosine_beam_primtrk, "cosine_beam_primtrk/D");
                fTree->Branch("primtrk_startx",&primtrk_startx);
                fTree->Branch("primtrk_starty",&primtrk_starty);
                fTree->Branch("primtrk_startz",&primtrk_startz);

                fTree->Branch("primtrk_endx",&primtrk_endx);
                fTree->Branch("primtrk_endy",&primtrk_endy);
                fTree->Branch("primtrk_endz",&primtrk_endz);

                fTree->Branch("primtrk_Dirx",&primtrk_Dirx);
                fTree->Branch("primtrk_Diry",&primtrk_Diry);
                fTree->Branch("primtrk_Dirz",&primtrk_Dirz);
                fTree->Branch("primtrklen",&primtrklen);
                fTree->Branch("primtrkID",&primtrkID);

                fTree->Branch("primtrk_dqdx",&primtrk_dqdx);
                fTree->Branch("primtrk_dedx",&primtrk_dedx);
                fTree->Branch("primtrk_resrange",&primtrk_resrange);
                //fTree->Branch("primtrk_deadwireresrange",&primtrk_deadwireresrange);
                fTree->Branch("primtrk_range",&primtrk_range);
                fTree->Branch("primtrk_hitx",&primtrk_hitx);
                fTree->Branch("primtrk_hity",&primtrk_hity);
                fTree->Branch("primtrk_hitz",&primtrk_hitz);
                fTree->Branch("primtrk_ke",&primtrk_ke);
                fTree->Branch("primtrk_pitch",&primtrk_pitch);
                fTree->Branch("primtrk_wid",&primtrk_wid);
                fTree->Branch("primtrk_pt",&primtrk_pt);
                fTree->Branch("primtrk_calo_hit_index",&primtrk_calo_hit_index);
                fTree->Branch("primtrk_ch",&primtrk_ch);

                fTree->Branch("pdg_code", &pdg_code);
                fTree->Branch("n_beamparticle", &n_beamparticle);
                fTree->Branch("n_daughter", &n_daughter);
                fTree->Branch("isPrimary", &isPrimary);
                fTree->Branch("pfp_self", &pfp_self);
                //fTree->Branch("pfp_parent", &pfp_parent);
                fTree->Branch("pfp_daughter", &pfp_daughter);
                fTree->Branch("primtrk_trktag", &primtrk_trktag);

                fTree->Branch("Zintersection",&Zintersection);
                fTree->Branch("Yintersection",&Yintersection);
                fTree->Branch("Xintersection",&Xintersection);
                fTree->Branch("timeintersection",&timeintersection);

        }

void protoana::protonbeamana::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 905 of file protonbeamana_module.cc.

        {

        }

protonbeamana& protoana::protonbeamana::operator= ( protonbeamana const &  )

delete

protonbeamana& protoana::protonbeamana::operator= ( protonbeamana &&  )

delete

void protoana::protonbeamana::reset

(

)

Definition at line 910 of file protonbeamana_module.cc.

        {
                for(int k=0; k < 5; k++)
                        fNactivefembs[k] = -999;


                beam_inst_valid = true;
                beam_inst_trigger = -999;
                beam_inst_nTracks = -999;
                beam_inst_nMomenta = -999;
                reco_reconstructable_beam_event = false;

                beam_inst_PDG_candidates.clear();
                beam_inst_TOF.clear();
                beam_inst_TOF_Chan.clear();

                tof = -1;
                tofs.clear();
                ch_tofs.clear();

                low_pressure_status = -1;
                high_pressure_status = -1;
                low_pressure = -99;
                high_pressure = -99;

                beamPosx.clear();
                beamPosy.clear();
                beamPosz.clear();
                beamDirx.clear();
                beamDiry.clear();
                beamDirz.clear();
                beamMomentum.clear();

                fisprimarytrack = 0;
                fisprimaryshower = 0;

                primtrk_startx.clear();
                primtrk_starty.clear();
                primtrk_startz.clear();  

                primtrk_endx.clear();
                primtrk_endy.clear();
                primtrk_endz.clear();  

                primtrk_Dirx.clear();
                primtrk_Diry.clear();
                primtrk_Dirz.clear();
                primtrklen.clear();
                primtrkID.clear();
                primtrk_trktag.clear();

                //primtrk_dqdx->clear();
                //primtrk_resrange->clear();
                //primtrk_dedx->clear();
                primtrk_dqdx.clear();
                primtrk_resrange.clear();
                //primtrk_deadwireresrange.clear();
                primtrk_dedx.clear();
                primtrk_range.clear();
                primtrk_hitx.clear();
                primtrk_hity.clear();
                primtrk_hitz.clear();
                primtrk_ke.clear();
                primtrk_pitch.clear();
                primtrk_wid.clear();
                primtrk_pt.clear();
                primtrk_calo_hit_index.clear(); 
                primtrk_ch.clear();     

                pdg_code.clear();
                n_beamparticle.clear();
                n_daughter.clear();

                isPrimary.clear();
                pfp_self.clear();
                //pfp_parent.clear();
                pfp_daughter.clear();

                cosine_beam_primtrk=-99;

                Zintersection.clear();
                Yintersection.clear();
                Xintersection.clear();
                timeintersection.clear();


        }

Member Data Documentation

int protoana::protonbeamana::beam_inst_nMomenta

private

Definition at line 139 of file protonbeamana_module.cc.

int protoana::protonbeamana::beam_inst_nTracks

private

Definition at line 139 of file protonbeamana_module.cc.

std::vector< int > protoana::protonbeamana::beam_inst_PDG_candidates

private

Definition at line 141 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beam_inst_TOF

private

Definition at line 140 of file protonbeamana_module.cc.

std::vector< int > protoana::protonbeamana::beam_inst_TOF_Chan

private

Definition at line 141 of file protonbeamana_module.cc.

int protoana::protonbeamana::beam_inst_trigger

private

Definition at line 138 of file protonbeamana_module.cc.

bool protoana::protonbeamana::beam_inst_valid

private

Definition at line 137 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamDirx

private

Definition at line 147 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamDiry

private

Definition at line 148 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamDirz

private

Definition at line 149 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamMomentum

private

Definition at line 151 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamPosx

private

Definition at line 143 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamPosy

private

Definition at line 144 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::beamPosz

private

Definition at line 145 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::ch_tofs

private

Definition at line 155 of file protonbeamana_module.cc.

double protoana::protonbeamana::cosine_beam_primtrk

private

Definition at line 220 of file protonbeamana_module.cc.

protoana::ProtoDUNEDataUtils protoana::protonbeamana::dataUtil

private

Definition at line 164 of file protonbeamana_module.cc.

int protoana::protonbeamana::event

private

Definition at line 131 of file protonbeamana_module.cc.

double protoana::protonbeamana::evttime

private

Definition at line 133 of file protonbeamana_module.cc.

protoana::ProtoDUNEBeamlineUtils protoana::protonbeamana::fBeamlineUtils

private

Definition at line 162 of file protonbeamana_module.cc.

const art::InputTag protoana::protonbeamana::fBeamModuleLabel

private

Definition at line 123 of file protonbeamana_module.cc.

fhicl::ParameterSet protoana::protonbeamana::fBeamPars

private

Definition at line 176 of file protonbeamana_module.cc.

std::string protoana::protonbeamana::fCalorimetryTag

private

Definition at line 168 of file protonbeamana_module.cc.

protoana::ProtoDUNEEmptyEventFinder protoana::protonbeamana::fEmptyEventFinder

private

Definition at line 163 of file protonbeamana_module.cc.

geo::GeometryCore const* protoana::protonbeamana::fGeometry

private

Definition at line 181 of file protonbeamana_module.cc.

std::string protoana::protonbeamana::fHitTag

private

Definition at line 170 of file protonbeamana_module.cc.

int protoana::protonbeamana::fisprimaryshower

private

Definition at line 202 of file protonbeamana_module.cc.

int protoana::protonbeamana::fisprimarytrack

private

Definition at line 201 of file protonbeamana_module.cc.

int protoana::protonbeamana::fNactivefembs[6]

private

Definition at line 134 of file protonbeamana_module.cc.

std::string protoana::protonbeamana::fPFParticleTag

private

Definition at line 172 of file protonbeamana_module.cc.

std::string protoana::protonbeamana::fShowerTag

private

Definition at line 171 of file protonbeamana_module.cc.

std::string protoana::protonbeamana::fTrackerTag

private

Definition at line 169 of file protonbeamana_module.cc.

const art::InputTag protoana::protonbeamana::fTrackModuleLabel

private

Definition at line 120 of file protonbeamana_module.cc.

TTree* protoana::protonbeamana::fTree

private

Definition at line 127 of file protonbeamana_module.cc.

bool protoana::protonbeamana::fVerbose

private

Definition at line 178 of file protonbeamana_module.cc.

double protoana::protonbeamana::high_pressure

private

Definition at line 159 of file protonbeamana_module.cc.

short protoana::protonbeamana::high_pressure_status

private

Definition at line 157 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::isPrimary

private

Definition at line 225 of file protonbeamana_module.cc.

double protoana::protonbeamana::low_pressure

private

Definition at line 158 of file protonbeamana_module.cc.

short protoana::protonbeamana::low_pressure_status

private

Definition at line 156 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::n_beamparticle

private

Definition at line 224 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::n_daughter

private

Definition at line 223 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::pdg_code

private

Definition at line 222 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::pfp_daughter

private

Definition at line 228 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::pfp_self

private

Definition at line 226 of file protonbeamana_module.cc.

std::vector<size_t> protoana::protonbeamana::primtrk_calo_hit_index

private

Definition at line 217 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::primtrk_ch

private

Definition at line 218 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_dedx

private

Definition at line 208 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_Dirx

private

Definition at line 192 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_Diry

private

Definition at line 193 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_Dirz

private

Definition at line 194 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_dqdx

private

Definition at line 205 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_endx

private

Definition at line 188 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_endy

private

Definition at line 189 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_endz

private

Definition at line 190 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_hitx

private

Definition at line 211 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_hity

private

Definition at line 212 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_hitz

private

Definition at line 213 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_ke

private

Definition at line 210 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_pitch

private

Definition at line 214 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_pt

private

Definition at line 216 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_range

private

Definition at line 209 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_resrange

private

Definition at line 206 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_startx

private

Definition at line 184 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_starty

private

Definition at line 185 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrk_startz

private

Definition at line 186 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::primtrk_trktag

private

Definition at line 197 of file protonbeamana_module.cc.

std::vector<int> protoana::protonbeamana::primtrk_wid

private

Definition at line 215 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrkID

private

Definition at line 196 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::primtrklen

private

Definition at line 195 of file protonbeamana_module.cc.

bool protoana::protonbeamana::reco_reconstructable_beam_event

private

Definition at line 165 of file protonbeamana_module.cc.

int protoana::protonbeamana::run

private

Definition at line 129 of file protonbeamana_module.cc.

int protoana::protonbeamana::subrun

private

Definition at line 130 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::timeintersection

private

Definition at line 234 of file protonbeamana_module.cc.

double protoana::protonbeamana::tof

private

Definition at line 153 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::tofs

private

Definition at line 154 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::Xintersection

private

Definition at line 233 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::Yintersection

private

Definition at line 232 of file protonbeamana_module.cc.

std::vector<double> protoana::protonbeamana::Zintersection

private

Definition at line 231 of file protonbeamana_module.cc.

---

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

• protoduneana/protoduneana/singlephase/protonanalysis/protonbeamana_module.cc