Inheritance diagram for ProtoDUNEelectronWireAna:

Public Member Functions
	ProtoDUNEelectronWireAna (fhicl::ParameterSet const &p)

	ProtoDUNEelectronWireAna (ProtoDUNEelectronWireAna const &)=delete

	ProtoDUNEelectronWireAna (ProtoDUNEelectronWireAna &&)=delete

ProtoDUNEelectronWireAna &	operator= (ProtoDUNEelectronWireAna const &)=delete

ProtoDUNEelectronWireAna &	operator= (ProtoDUNEelectronWireAna &&)=delete

void	analyze (art::Event const &evt) override

void	beginJob () override

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 Member Functions
void	Initialize ()

Private Attributes
protoana::ProtoDUNEPFParticleUtils	pfpUtil

protoana::ProtoDUNEShowerUtils	showerUtil

protoana::ProtoDUNETruthUtils	truthUtil

geo::GeometryCore const *	fGeometry = &*(art::ServiceHandle<geo::Geometry>())

std::string	fShowerTag

std::string	fPFParticleTag

std::string	fWireTag

TTree *	fTree

unsigned int	frun

unsigned int	fsubrun

unsigned int	fevent

std::vector< int >	fprimary_Shower_wire_w

std::vector< float >	fprimary_Shower_wire_ch

std::vector< float >	fprimary_Shower_wire_X

std::vector< float >	fprimary_Shower_wire_Z

std::vector< float >	fprimary_Shower_wire_Y

std::vector< float >	fprimary_Shower_MCwire_E

std::vector< int >	fprimary_Shower_MCwire_w

float	fprimaryStartPosition [3]

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 54 of file ProtoDUNEelectronWireAna_module.cc.

Constructor & Destructor Documentation

ProtoDUNEelectronWireAna::ProtoDUNEelectronWireAna ( fhicl::ParameterSet const & p )

explicit

Definition at line 119 of file ProtoDUNEelectronWireAna_module.cc.

   : EDAnalyzer{p},
   fShowerTag(p.get<std::string>("ShowerTag")),
   fPFParticleTag(p.get<std::string>("PFParticleTag")),
   fWireTag(p.get<std::string>("WireTag"))
 {
   // Call appropriate consumes<>() for any products to be retrieved by this module.
 }

ProtoDUNEelectronWireAna::ProtoDUNEelectronWireAna ( ProtoDUNEelectronWireAna const & )

delete

ProtoDUNEelectronWireAna::ProtoDUNEelectronWireAna ( ProtoDUNEelectronWireAna && )

delete

Member Function Documentation

void ProtoDUNEelectronWireAna::analyze ( art::Event const & evt )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 128 of file ProtoDUNEelectronWireAna_module.cc.

 {
 
   Initialize(); 
   frun = evt.run();
   fsubrun = evt.subRun();
   fevent  = evt.id().event(); 
  
   //check for reco pandora stuff
   auto recoParticleHandle = evt.getHandle<std::vector<recob::PFParticle>>(fPFParticleTag);
   if( !recoParticleHandle ) return;
 
   //std::vector<int> hit_w;
   std::map<int,std::vector<double>> hit_w_and_t1;
   std::map<int,std::vector<double>> hit_w_and_t2;
   std::map<int,std::vector<double>> hit_w_and_y;
   std::map<int,std::vector<double>> hit_w_and_x;
   std::vector<const recob::PFParticle*> pfParticles = pfpUtil.GetPFParticlesFromBeamSlice(evt,fPFParticleTag);
   const simb::MCParticle* mcparticle = NULL;
   bool doAna = false;
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(evt);
   auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataFor(evt, clockData);
   for(const recob::PFParticle* particle : pfParticles){
      const recob::Shower* thisShower = pfpUtil.GetPFParticleShower(*particle,evt,fPFParticleTag,fShowerTag);
      if( thisShower == 0x0 ) continue;
      if( !evt.isRealData() ){
        mcparticle = truthUtil.GetMCParticleFromRecoShower(clockData, *thisShower, evt, "pandoraShower");
        if( abs(mcparticle->PdgCode()) != 11 ) return;
      }
      fprimaryStartPosition[0] = thisShower->ShowerStart().X();
      fprimaryStartPosition[1] = thisShower->ShowerStart().Y();
      fprimaryStartPosition[2] = thisShower->ShowerStart().Z();
      doAna = true; 
      const std::vector<const recob::Hit*> sh_hits = showerUtil.GetRecoShowerHits(*thisShower, evt, fShowerTag);
      art::FindManyP<recob::Wire> wFromHits(sh_hits,evt,"hitpdune");
      art::FindManyP<recob::SpacePoint> spFromShowerHits(sh_hits,evt,fPFParticleTag);
      for( size_t j=0; j<sh_hits.size() && j<MAXHits; ++j){
        if( sh_hits[j]->WireID().Plane != 2 ) continue;
        std::vector<art::Ptr<recob::Wire>> wires = wFromHits.at(j);
        //hit_w.push_back(wires[0]->Channel());
        hit_w_and_t1[wires[0]->Channel()].push_back(sh_hits[j]->PeakTimeMinusRMS(5.0));
        hit_w_and_t2[wires[0]->Channel()].push_back(sh_hits[j]->PeakTimePlusRMS(5.0));
        //std::cout<<wires[0]->Channel()<<" "<<sh_hits[j]->PeakTime()<<" "<<sh_hits[j]->PeakTimePlusRMS(5.0)<<" "<<sh_hits[j]->PeakTimeMinusRMS(5.0)<<std::endl;
        hit_w_and_x[wires[0]->Channel()].push_back(detProp.ConvertTicksToX(sh_hits[j]->PeakTime(),sh_hits[j]->WireID().Plane,sh_hits[j]->WireID().TPC,0));
        std::vector<art::Ptr<recob::SpacePoint>> sp = spFromShowerHits.at(j); 
        if(!sp.empty()){
           hit_w_and_y[wires[0]->Channel()].push_back(sp[0]->XYZ()[1]);
        }
        else hit_w_and_y[wires[0]->Channel()].push_back(fprimaryStartPosition[1]); //use vtx if no sp 
      }   
      break;
   }
  
   //one wire can have various hits so lets remove duplicate wires
   //std::sort(hit_w.begin(),hit_w.end());
   //hit_w.erase(std::unique(hit_w.begin(),hit_w.end()), hit_w.end()); 
   
   if( doAna ){
     auto const& wires = evt.getValidHandle<std::vector<recob::Wire> >(fWireTag);
     //auto const& wires = evt.getValidHandle<std::vector<recob::Wire> >("wclsdatasp:gauss"); //new reco 
     auto w1 = hit_w_and_t1.begin();
     auto w2 = hit_w_and_t2.begin();
     auto x  = hit_w_and_x.begin();
     auto y  = hit_w_and_y.begin(); 
     while( w1 != hit_w_and_t1.end()){
       int it_w = w1->first;
       int n_hits = w1->second.size();
       // std::cout<<it_w<<" "<<n_hits<<" "<<w1->second[n_hits-1]<<" "<<w2->second[n_hits-1]<<std::endl;
       double t1 =  w1->second[0]; //first hit
       double t2 =  w2->second[n_hits-1]; //last hit
       double x_w, y_w;
       if( x->second.size() < 1 ){
         x_w = (x->second[0]-x->second[n_hits-1])/2.0;
         y_w = (y->second[0]-y->second[n_hits-1])/2.0;
       }
       else {
         x_w = x->second[0];
         y_w = y->second[0];
       }
       for(auto & wire : * wires){
          int channel_no = wire.Channel();
          int plane = fGeometry->View(wire.Channel()); 
          if( plane != 2 ) continue;
          std::vector< geo::WireID > wireID= fGeometry->ChannelToWire(channel_no);
          const geo::WireGeo* pwire = fGeometry->WirePtr(wireID.at(0)); //for collection plane there is one wire per channel
          TVector3 xyzWire = pwire->GetCenter<TVector3>();
          if( it_w == channel_no ){  
            double charge =0.0;
            for(size_t i = 0; i < wire.Signal().size(); ++i){
               if( i > t1 && i < t2 ) charge += wire.Signal()[i];
            }   
            fprimary_Shower_wire_Z.push_back(xyzWire.Z()); 
            fprimary_Shower_wire_Y.push_back(y_w); 
            fprimary_Shower_wire_X.push_back(x_w); 
            fprimary_Shower_wire_ch.push_back(charge);
            fprimary_Shower_wire_w.push_back(channel_no);
            
            break;
          }   
       }// all recob::wire
       w1 ++; w2 ++;
       x ++; y ++;
     }//wire from shower 
 
     //look at MC info if available
     if(!evt.isRealData()){ 
       auto mcTruths = evt.getValidHandle<std::vector<simb::MCTruth>>("generator");
       const simb::MCParticle* geantGoodParticle = truthUtil.GetGeantGoodParticle((*mcTruths)[0],evt);
       int trackid = mcparticle->TrackId();
       if( mcparticle->TrackId() == geantGoodParticle->TrackId() ){
         if( doAna ){   //we have a reco shower
           auto simchannelHandle = evt.getHandle< std::vector<sim::SimChannel> >("largeant");
           if (simchannelHandle){
             for(auto const& simchannel : (*simchannelHandle)){
                if(fGeometry->View(simchannel.Channel()) != 2) continue;
                auto const& alltimeslices = simchannel.TDCIDEMap();
                double EperCh=0;
                for(auto const& tslice : alltimeslices){
                   auto const& simide = tslice.second;
                   // Loop over energy deposits
                   for(auto const& eDep : simide){
                      if(eDep.trackID == trackid || eDep.trackID == -trackid){
                        EperCh += eDep.energy;
                      }
                   }
                } 
                if( EperCh== 0 ) continue;//save only channeles with ID
                fprimary_Shower_MCwire_w.push_back(simchannel.Channel());
                fprimary_Shower_MCwire_E.push_back(EperCh);
             }
           }
         }//do analysis of MC
       }//is MC same as beam?
     }//is MC?
   }//shower
 
 
   fTree->Fill();
 }

void ProtoDUNEelectronWireAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 98 of file ProtoDUNEelectronWireAna_module.cc.

 {
 
   art::ServiceHandle<art::TFileService> tfs;
   fTree = tfs->make<TTree>("WireAna","Recob::Wire");
   fTree->Branch("run",&frun,"run/I");
   fTree->Branch("subrun",&fsubrun,"subrun/I");
   fTree->Branch("event",&fevent,"event/I");
   
   fTree->Branch("primaryStartPosition",          &fprimaryStartPosition,         "primaryStartPosition[3]/f");
   fTree->Branch("primary_Shower_wire_w",&fprimary_Shower_wire_w);
   fTree->Branch("primary_Shower_wire_ch",&fprimary_Shower_wire_ch);
   fTree->Branch("primary_Shower_wire_X",&fprimary_Shower_wire_X);
   fTree->Branch("primary_Shower_wire_Z",&fprimary_Shower_wire_Z);
   fTree->Branch("primary_Shower_wire_Y",&fprimary_Shower_wire_Y);
   fTree->Branch("primary_Shower_MCwire_w",&fprimary_Shower_MCwire_w);
   fTree->Branch("primary_Shower_MCwire_E",&fprimary_Shower_MCwire_E);
  
 }

void ProtoDUNEelectronWireAna::Initialize ( void )

private

Definition at line 267 of file ProtoDUNEelectronWireAna_module.cc.

                                          {
   frun =-999;
   fsubrun =-999;
   fevent =-999; 
   fprimary_Shower_wire_w.clear();
   fprimary_Shower_wire_Y.clear();
   fprimary_Shower_wire_X.clear();
   fprimary_Shower_wire_Z.clear();
   fprimary_Shower_wire_ch.clear();
   fprimary_Shower_MCwire_w.clear();
   fprimary_Shower_MCwire_E.clear();
   for(int k=0; k < 3; k++){
     fprimaryStartPosition[k] = -999.0;
   }
 }