DAQQuickClustering Class Reference

#include <DAQQuickClustering_module.h>

Inheritance diagram for DAQQuickClustering:

Public Member Functions
	DAQQuickClustering (fhicl::ParameterSet const &p)
	DAQQuickClustering (DAQQuickClustering const &)=delete
	DAQQuickClustering (DAQQuickClustering &&)=delete
DAQQuickClustering &	operator= (DAQQuickClustering const &)=delete
DAQQuickClustering &	operator= (DAQQuickClustering &&)=delete
void	analyze (art::Event const &evt) override
void	reconfigure (fhicl::ParameterSet const &p)
void	beginJob () override
void	endJob () 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	ResetVariables ()
void	trigger (std::vector< cluster > &vec_Clusters, unsigned int const &config)
void	clusterCut (std::vector< cluster > &vec_Clusters, unsigned int const &config)
void	clusterChannels (std::vector< recoHit > &vec_Hits, std::vector< ClusterHitsInTime > &vec_ChannelCluster, unsigned int const &config)
void	makeConfigGraph ()
void	FillMyMaps (std::map< int, simb::MCParticle > &MyMap, art::FindManyP< simb::MCParticle > Assn, art::ValidHandle< std::vector< simb::MCTruth > > Hand)
PType	WhichParType (int TrID)
bool	InMyMap (int TrID, std::map< int, simb::MCParticle > ParMap)

Private Attributes
std::vector< int >	cut_AdjChanTolerance
std::vector< int >	cut_HitsInWindow
std::vector< int >	cut_MinChannels
std::vector< int >	cut_MinChanWidth
std::vector< float >	cut_TimeWindowSize
std::vector< float >	cut_TotalADC
int	firstCatch
int	secondCatch
int	thirdCatch
std::string	fRawDigitLabel
std::string	fHitLabel
std::string	fGEANTLabel
std::string	fMARLLabel
std::map< int, simb::MCParticle >	MarlParts
std::string	fAPALabel
std::map< int, simb::MCParticle >	APAParts
std::string	fCPALabel
std::map< int, simb::MCParticle >	CPAParts
std::string	fAr39Label
std::map< int, simb::MCParticle >	Ar39Parts
std::string	fNeutLabel
std::map< int, simb::MCParticle >	NeutParts
std::string	fKrypLabel
std::map< int, simb::MCParticle >	KrypParts
std::string	fPlonLabel
std::map< int, simb::MCParticle >	PlonParts
std::string	fRdonLabel
std::map< int, simb::MCParticle >	RdonParts
std::string	fAr42Label
std::map< int, simb::MCParticle >	Ar42Parts
std::map< int, PType >	trkIDToPType
unsigned int	NConfigs
unsigned int	NCuts
double	detectorScaling
int	Run
int	SubRun
int	Event
std::map< int, std::vector< double > >	map_EventToMC
TTree *	t_Output_unusedhits
TTree *	t_Output_clusteredhits
TH1D *	h_ENu_MC
TH1D *	h_MarlTime_MC
TH1D *	h_TimeElapsed
int	out_Cluster
int	out_Event
int	out_Config
int	out_StartChan
int	out_EndChan
int	out_ChanWidth
int	out_NChan
int	out_Type
int	out_NHits
float	out_SumADC
float	out_FirstTimeHit
float	out_LastTimeHit
float	out_TimeWidth
double	out_ENu
double	out_ENu_Lep
double	out_MarlTime
std::vector< int >	out_HitView
std::vector< int >	out_GenType
std::vector< int >	out_HitChan
std::vector< double >	out_HitTime
std::vector< double >	out_HitSADC
std::vector< double >	out_HitRMS
int	NTotHits
int	NColHits
int	NIndHits
int	nHitsNotBackTracked
int	HitView [nMaxHits]
int	HitSize [nMaxHits]
int	HitTPC [nMaxHits]
int	HitChan [nMaxHits]
float	HitTime [nMaxHits]
float	HitRMS [nMaxHits]
float	HitSADC [nMaxHits]
float	HitInt [nMaxHits]
float	HitPeak [nMaxHits]
int	GenType [nMaxHits]
float	Hit_X [nMaxHits]
float	Hit_Y [nMaxHits]
float	Hit_Z [nMaxHits]
float	Hit_Energy [nMaxHits]
float	Hit_NumElectrons [nMaxHits]
int	NCorrespondingIDEs [nMaxHits]
int	TotGen_Marl
int	TotGen_APA
int	TotGen_CPA
int	TotGen_Ar39
int	TotGen_Neut
int	TotGen_Kryp
int	TotGen_Plon
int	TotGen_Rdon
int	TotGen_Ar42
int	VertexChan
int	Nu_Type
int	Nu_Lep_Type
int	Mode
int	CCNC
int	HitNucleon
int	Target
std::vector< int >	MarlSample
std::vector< double >	MarlTime
std::vector< double >	MarlWeight
double	ENu
double	ENu_Lep
double	VertX
double	VertY
double	VertZ
double	VertexT
double	Px
double	Py
double	Pz
art::ServiceHandle< geo::Geometry >	geo
art::ServiceHandle< cheat::BackTrackerService >	bt_serv
art::ServiceHandle< cheat::ParticleInventoryService >	pi_serv

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 193 of file DAQQuickClustering_module.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 97 of file DAQQuickClustering_module.cc.

                                                                 :EDAnalyzer(p){

   this->reconfigure(p);

}

DAQQuickClustering::DAQQuickClustering ( DAQQuickClustering const &  )

delete

DAQQuickClustering::DAQQuickClustering ( DAQQuickClustering &&  )

delete

Member Function Documentation

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

overridevirtual

Implements art::EDAnalyzer.

Definition at line 407 of file DAQQuickClustering_module.cc.

{
  ResetVariables();
  
  Run    = evt.run();
  SubRun = evt.subRun();
  Event  = evt.event();

  //GET INFORMATION ABOUT THE DETECTOR'S GEOMETRY.
  auto const* geo = lar::providerFrom<geo::Geometry>();

// GET THE RECO HITS.
  auto reco_hits = evt.getValidHandle<std::vector<recob::Hit> >(fHitLabel);
  
  //LIFT OUT THE MARLEY PARTICLES.
  auto MarlTrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fMARLLabel);
  art::FindManyP<simb::MCParticle> MarlAssn(MarlTrue,evt,fGEANTLabel);
  FillMyMaps(MarlParts, MarlAssn, MarlTrue);
  TotGen_Marl = MarlParts.size();

  //SUPERNOVA TRUTH.
  art::FindManyP<sim::SupernovaTruth> SNTruth(MarlTrue, evt, fMARLLabel);

  double Px_(0), Py_(0), Pz_(0), Pnorm(1);
  for(unsigned int i = 0; i < MarlTrue->size(); i++)
  {
    Nu_Type     = MarlTrue->at(i).GetNeutrino().Nu().PdgCode();
    ENu         = MarlTrue->at(i).GetNeutrino().Nu().E();
    Mode        = MarlTrue->at(i).GetNeutrino().Mode();
    CCNC        = MarlTrue->at(i).GetNeutrino().CCNC();
    Target      = MarlTrue->at(i).GetNeutrino().Target();
    HitNucleon  = MarlTrue->at(i).GetNeutrino().HitNuc();
    Nu_Lep_Type = MarlTrue->at(i).GetNeutrino().Lepton().PdgCode(); 
    VertX       = MarlTrue->at(i).GetNeutrino().Lepton().Vx(); 
    VertY       = MarlTrue->at(i).GetNeutrino().Lepton().Vy(); 
    VertZ       = MarlTrue->at(i).GetNeutrino().Lepton().Vz(); 
    Px_         = MarlTrue->at(i).GetNeutrino().Lepton().Px();
    Py_         = MarlTrue->at(i).GetNeutrino().Lepton().Py();
    Pz_         = MarlTrue->at(i).GetNeutrino().Lepton().Pz();
    ENu_Lep     = MarlTrue->at(i).GetNeutrino().Lepton().E();
    for (unsigned int j = 0; j < SNTruth.at(i).size(); j++) 
    {
      const sim::SupernovaTruth ThisTr = (*SNTruth.at(i).at(j));
      MarlTime  .push_back(ThisTr.SupernovaTime);
      MarlWeight.push_back(ThisTr.Weight);
      MarlSample.push_back(ThisTr.SamplingMode);
    }
  }
  Pnorm = std::sqrt(Px_*Px_+Py_*Py_+Pz_*Pz_);
  Px = Px_/Pnorm;
  Py = Py_/Pnorm;
  Pz = Pz_/Pnorm;
  bool caught = false;
  double Vertex[3] = {VertX, VertY, VertZ};
  geo::WireID WireID;
  geo::PlaneID Plane(geo->FindTPCAtPosition(Vertex),geo::kZ);
  try
  {
    WireID = geo->NearestWireID(Vertex, Plane);
  }
  catch(...)
  {
    caught = true;
  }

  if(caught)
    VertexChan = -1; 
  else
    VertexChan = geo->PlaneWireToChannel(WireID);
  
  //CM/MICROSECOND.
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
  double drift_velocity = detProp.DriftVelocity(detProp.Efield(),detProp.Temperature());
  //CM/TICK
  drift_velocity = drift_velocity*0.5;
  VertexT = VertX/drift_velocity;

  auto APATrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fAPALabel);
  art::FindManyP<simb::MCParticle> APAAssn(APATrue,evt,fGEANTLabel);
  FillMyMaps( APAParts, APAAssn, APATrue );
  TotGen_APA = APAParts.size();

  auto CPATrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fCPALabel);
  art::FindManyP<simb::MCParticle> CPAAssn(CPATrue,evt,fGEANTLabel);
  FillMyMaps( CPAParts, CPAAssn, CPATrue );
  TotGen_CPA = CPAParts.size();

  auto Ar39True = evt.getValidHandle<std::vector<simb::MCTruth> >(fAr39Label);
  art::FindManyP<simb::MCParticle> Ar39Assn(Ar39True,evt,fGEANTLabel);
  FillMyMaps( Ar39Parts, Ar39Assn, Ar39True );
  TotGen_Ar39 = Ar39Parts.size();

  auto NeutTrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fNeutLabel);
  art::FindManyP<simb::MCParticle> NeutAssn(NeutTrue,evt,fGEANTLabel);
  FillMyMaps( NeutParts, NeutAssn, NeutTrue );
  TotGen_Neut = NeutParts.size();

  auto KrypTrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fKrypLabel);
  art::FindManyP<simb::MCParticle> KrypAssn(KrypTrue,evt,fGEANTLabel);
  FillMyMaps( KrypParts, KrypAssn, KrypTrue );
  TotGen_Kryp = KrypParts.size();

  auto PlonTrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fPlonLabel);
  art::FindManyP<simb::MCParticle> PlonAssn(PlonTrue,evt,fGEANTLabel);
  FillMyMaps( PlonParts, PlonAssn, PlonTrue );
  TotGen_Plon = PlonParts.size();

  auto RdonTrue = evt.getValidHandle<std::vector<simb::MCTruth> >(fRdonLabel);
  art::FindManyP<simb::MCParticle> RdonAssn(RdonTrue,evt,fGEANTLabel);
  FillMyMaps( RdonParts, RdonAssn, RdonTrue );
  TotGen_Rdon = RdonParts.size();

  auto Ar42True = evt.getValidHandle<std::vector<simb::MCTruth> >(fAr42Label);
  art::FindManyP<simb::MCParticle> Ar42Assn(Ar42True,evt,fGEANTLabel);
  FillMyMaps( Ar42Parts, Ar42Assn, Ar42True );
  TotGen_Ar42 = Ar42Parts.size();
  
  std::map<PType, std::map< int, simb::MCParticle >&> PTypeToMap{
      { kMarl, MarlParts},
      {kAPA, APAParts},
      {kCPA, CPAParts},
      {kAr39, Ar39Parts},
      {kNeut, NeutParts},
      {kKryp, KrypParts},
      {kPlon, PlonParts},
      {kRdon, RdonParts},
      {kAr42, Ar42Parts}
  };

  for(auto const& it : PTypeToMap){
      const PType p=it.first;
      auto const& m=it.second;
      for(auto const& it2 : m){
          trkIDToPType.insert(std::make_pair(it2.first, p));
      }
  }
  
  //std::cout << "THE EVENTS NUMBER IS: " << Event << std::endl;

  std::vector< recob::Hit > ColHits_Marl;
  std::vector< recob::Hit > ColHits_CPA;
  std::vector< recob::Hit > ColHits_APA;
  std::vector< recob::Hit > ColHits_Ar39;
  std::vector< recob::Hit > ColHits_Neut;
  std::vector< recob::Hit > ColHits_Kryp;
  std::vector< recob::Hit > ColHits_Plon;
  std::vector< recob::Hit > ColHits_Rdon;
  std::vector< recob::Hit > ColHits_Oth;
  std::vector< recob::Hit > ColHits_Ar42;

  NTotHits = reco_hits->size();
  int colHitCount(0);
  int LoopHits = std::min( NTotHits, nMaxHits );
  //if(LoopHits != NTotHits)
  // std::cerr << "---- There are " << NTotHits
  //           << " hits in the event, but array is of size " 
  //           << nMaxHits << ", so looping over first " << LoopHits
  //           << " hits." << std::endl;

  for(int hit = 0; hit < LoopHits; ++hit) 
  {
    recob::Hit const& ThisHit = reco_hits->at(hit);  
    if(ThisHit.View() == geo::kU || ThisHit.View() == geo::kV) ++NIndHits;
    else                                                       ++NColHits;
  }

  for(int hit = 0; hit < LoopHits; ++hit) 
  {
    recob::Hit const& ThisHit = reco_hits->at(hit);  

    if (ThisHit.View() == 2) {
    std::vector< sim::TrackIDE > ThisHitIDE; 
    //GETTING HOLD OF THE SIM::IDEs.
    std::vector<const sim::IDE*> ThisSimIDE;
    try
    {
      ThisHitIDE = bt_serv->HitToTrackIDEs(clockData,  ThisHit );
    }
    catch(...)
    {
      //std::cout << "FIRST CATCH" << std::endl;
      firstCatch++;
      try
      {
        ThisSimIDE = bt_serv->HitToSimIDEs_Ps(clockData, ThisHit);
      }
      catch(...)
      {
        //std::cout << "SECOND CATCH" << std::endl; 
        secondCatch++;
        continue;
      }
      continue;
    }
    try
    {
      ThisSimIDE = bt_serv->HitToSimIDEs_Ps(clockData, ThisHit);
    }
    catch(...)
    {
      //std::cout << "THIRD CATCH" << std::endl; 
      thirdCatch++;
      continue;
    }
    
    HitView[colHitCount] = ThisHit.View();
    HitSize[colHitCount] = ThisHit.EndTick() - ThisHit.StartTick();
    HitTPC [colHitCount] = ThisHit.WireID().TPC;
    HitChan[colHitCount] = ThisHit.Channel();
    HitTime[colHitCount] = ThisHit.PeakTime();
    HitRMS [colHitCount] = ThisHit.RMS();
    HitSADC[colHitCount] = ThisHit.SummedADC();
    HitInt [colHitCount] = ThisHit.Integral();
    HitPeak[colHitCount] = ThisHit.PeakAmplitude();
    NCorrespondingIDEs[colHitCount] = ThisHitIDE.size();
    
    if(ThisHitIDE.size()==0)
      nHitsNotBackTracked++;

    //WHICH PARTICLE CONTRIBUTED MOST TO THE HIT.
    int    MainTrID = -1;
    double TopEFrac = -DBL_MAX;
    for (size_t ideL=0; ideL < ThisHitIDE.size(); ++ideL) 
    {
      if ( ThisHitIDE[ideL].energyFrac > TopEFrac ) 
      {
        TopEFrac = ThisHitIDE[ideL].energyFrac;
        MainTrID = ThisHitIDE[ideL].trackID;
      }
    }

    PType ThisPType      = WhichParType( MainTrID );
    GenType[colHitCount] = ThisPType;

    if(MainTrID == -1)
    {
      Hit_X[colHitCount]            = -1;
      Hit_Y[colHitCount]            = -1;
      Hit_Z[colHitCount]            = -1;
      Hit_Energy[colHitCount]       = -1;
      Hit_NumElectrons[colHitCount] = -1;
    }
    else
    {
      for(unsigned int i = 0; i < ThisSimIDE.size(); i++)
      {
        if(ThisSimIDE.at(i)->trackID==MainTrID)
        {
          Hit_X[colHitCount]            = ThisSimIDE.at(i)->x;
          Hit_Y[colHitCount]            = ThisSimIDE.at(i)->y;
          Hit_Z[colHitCount]            = ThisSimIDE.at(i)->z;
          Hit_Energy[colHitCount]       = ThisSimIDE.at(i)->energy;
          Hit_NumElectrons[colHitCount] = ThisSimIDE.at(i)->numElectrons;
          break;
        }
      }
    }
    switch (ThisPType){
    case kUnknown: { ColHits_Oth .push_back( ThisHit ); break; }
    case kMarl:    { ColHits_Marl.push_back( ThisHit ); break; }
    case kAPA:     { ColHits_APA .push_back( ThisHit ); break; }
    case kCPA:     { ColHits_CPA .push_back( ThisHit ); break; }
    case kAr39:    { ColHits_Ar39.push_back( ThisHit ); break; }
    case kNeut:    { ColHits_Neut.push_back( ThisHit ); break; }
    case kKryp:    { ColHits_Kryp.push_back( ThisHit ); break; }
    case kPlon:    { ColHits_Plon.push_back( ThisHit ); break; }
    case kRdon:    { ColHits_Rdon.push_back( ThisHit ); break; }
    case kAr42:    { ColHits_Ar42.push_back( ThisHit ); break; }
    default: break;
    }

    colHitCount++;
    }
  } 

  std::vector<int> vec_ClusterCount(NConfigs);
  
  h_ENu_MC     ->Fill(1000*ENu);
  h_MarlTime_MC->Fill(MarlTime.back());
  
  map_EventToMC[Event] = {ENu, ENu_Lep, MarlTime.back()};
  
  //MAKE RECOHIT OBJECTS EVENTWISE FROM THE TREE.
  std::vector<recoHit> vec_Hits;
  for(int j = 0; j < NColHits; j++)
  {
    recoHit hit(Event,      HitView[j], GenType[j],
                HitChan[j], HitTime[j], HitSADC[j], HitRMS[j]);
    vec_Hits.push_back(hit);
  }

  for(unsigned int j = 0; j < 1; j++)
  {
    std::vector<ClusterHitsInTime> vec_ChannelCluster;
    std::vector<cluster>           vec_Clusters;
    TStopwatch *timeElapsed = new TStopwatch();
    clusterChannels(vec_Hits, vec_ChannelCluster, j);

    for(unsigned int k = 0; k < vec_ChannelCluster.size(); k++)
    {
      std::vector<cluster> vec_Temp = vec_ChannelCluster.at(k).getClusterVector();
      //std::cout << " vec_temp.size() " << vec_Temp.size() << std::endl;
      for(unsigned int l = 0; l < vec_Temp.size(); l++)
        vec_Clusters.push_back(vec_Temp.at(l));
    }
    clusterCut(vec_Clusters, j);
    trigger   (vec_Clusters, j);
    h_TimeElapsed->Fill(timeElapsed->RealTime()*1000);
    delete timeElapsed;
    timeElapsed = NULL;
    //std::cout << " vec_Cluster.size() " << vec_Clusters.size() << std::endl;

    //FILL THE OUTPUT TREE.
    for(unsigned int k  = 0; k  < vec_Clusters.size(); k++)
    {
      if(vec_Clusters.at(k).getTriggerFlag()==1)
      {
        out_Config       = j;
        out_Cluster      = vec_ClusterCount.at(j);
        out_Event        = vec_Clusters.at(k).getEvent();
        out_StartChan    = vec_Clusters.at(k).getStartChan();
        out_EndChan      = vec_Clusters.at(k).getEndChan();
        out_ChanWidth    = vec_Clusters.at(k).getChanWidth();
        out_NChan        = vec_Clusters.at(k).getNChan();
        out_Type         = vec_Clusters.at(k).getType();
        out_NHits        = vec_Clusters.at(k).getNHits();
        out_SumADC       = vec_Clusters.at(k).getHitSADC();
        out_FirstTimeHit = vec_Clusters.at(k).getFirstTimeHit();
        out_LastTimeHit  = vec_Clusters.at(k).getLastTimeHit();
        out_TimeWidth    = vec_Clusters.at(k).getTimeWidth();
        out_ENu          = map_EventToMC[vec_Clusters.at(k).getEvent()].at(0);
        out_ENu_Lep      = map_EventToMC[vec_Clusters.at(k).getEvent()].at(1);
        out_MarlTime     = map_EventToMC[vec_Clusters.at(k).getEvent()].at(2);
        for(unsigned int l = 0; l < vec_Clusters.at(k).getHits().size(); l++)
        {
          out_HitView.push_back(vec_Clusters.at(k).getHits().at(l).getHitView());  
          out_GenType.push_back(vec_Clusters.at(k).getHits().at(l).getGenType());  
          out_HitChan.push_back(vec_Clusters.at(k).getHits().at(l).getHitChan());  
          out_HitTime.push_back(vec_Clusters.at(k).getHits().at(l).getHitTime());  
          out_HitSADC.push_back(vec_Clusters.at(k).getHits().at(l).getHitSADC());  
          out_HitRMS .push_back(vec_Clusters.at(k).getHits().at(l).getHitRMS ());  
        }
        t_Output_clusteredhits->Fill();
        out_HitView.clear(); out_GenType.clear(); out_HitChan.clear(); 
        out_HitTime.clear(); out_HitSADC.clear(); out_HitRMS.clear();

        vec_ClusterCount.at(j)++;
      }
    }
    vec_ChannelCluster.clear();
    vec_Clusters.clear();
  }
  vec_Hits.clear();
}

void DAQQuickClustering::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 293 of file DAQQuickClustering_module.cc.

{
     
  art::ServiceHandle<art::TFileService> tfs;

  makeConfigGraph();

  t_Output_unusedhits = tfs->make<TTree>("DAQQuickClustering_unusedhits",
                              "DAQQuickClustering_unusedhits");

  t_Output_unusedhits->Branch("Cluster",      &out_Cluster,     "Cluster/I");
  t_Output_unusedhits->Branch("Event",        &out_Event,       "Event/I");
  t_Output_unusedhits->Branch("Config",       &out_Config,      "Config/I");
  t_Output_unusedhits->Branch("StartChan",    &out_StartChan,   "StartChan/I");
  t_Output_unusedhits->Branch("EndChan",      &out_EndChan,     "EndChan/I");
  t_Output_unusedhits->Branch("ChanWidth",    &out_ChanWidth,   "ChanWidth/I");
  t_Output_unusedhits->Branch("NChan",        &out_NChan,       "NChan/I");
  t_Output_unusedhits->Branch("Type",         &out_Type,        "Type/I");
  t_Output_unusedhits->Branch("NHits",        &out_NHits,       "NHits/I");
  t_Output_unusedhits->Branch("SumADC",       &out_SumADC,      "SumADC/F");
  t_Output_unusedhits->Branch("FirstTimeHit", &out_FirstTimeHit,"FirstTimeHit/F");
  t_Output_unusedhits->Branch("LastTimeHit",  &out_LastTimeHit, "LastTimeHit/F");
  t_Output_unusedhits->Branch("TimeWidth",    &out_TimeWidth,   "TimeWidth/F");
  t_Output_unusedhits->Branch("ENu",          &out_ENu,         "ENu/D");
  t_Output_unusedhits->Branch("ENu_Lep",      &out_ENu_Lep,     "ENu_Lep/D");
  t_Output_unusedhits->Branch("MarlTime",     &out_MarlTime,    "MarlTime/D");
  t_Output_unusedhits->Branch("HitView",      &out_HitView);
  t_Output_unusedhits->Branch("GenType",      &out_GenType);
  t_Output_unusedhits->Branch("HitChan",      &out_HitChan);
  t_Output_unusedhits->Branch("HitTime",      &out_HitTime);
  t_Output_unusedhits->Branch("HitSADC",      &out_HitSADC);
  t_Output_unusedhits->Branch("HitRMS",       &out_HitRMS);

  t_Output_clusteredhits = tfs->make<TTree>("DAQQuickClustering_clusteredhits",
                                            "DAQQuickClustering_clusteredhits");

  t_Output_clusteredhits->Branch("Cluster",      &out_Cluster,     "Cluster/I");
  t_Output_clusteredhits->Branch("Event",        &out_Event,       "Event/I");
  t_Output_clusteredhits->Branch("Config",       &out_Config,      "Config/I");
  t_Output_clusteredhits->Branch("StartChan",    &out_StartChan,   "StartChan/I");
  t_Output_clusteredhits->Branch("EndChan",      &out_EndChan,     "EndChan/I");
  t_Output_clusteredhits->Branch("ChanWidth",    &out_ChanWidth,   "ChanWidth/I");
  t_Output_clusteredhits->Branch("NChan",        &out_NChan,       "NChan/I");
  t_Output_clusteredhits->Branch("Type",         &out_Type,        "Type/I");
  t_Output_clusteredhits->Branch("NHits",        &out_NHits,       "NHits/I");
  t_Output_clusteredhits->Branch("SumADC",       &out_SumADC,      "SumADC/F");
  t_Output_clusteredhits->Branch("FirstTimeHit", &out_FirstTimeHit,"FirstTimeHit/F");
  t_Output_clusteredhits->Branch("LastTimeHit",  &out_LastTimeHit, "LastTimeHit/F");
  t_Output_clusteredhits->Branch("TimeWidth",    &out_TimeWidth,   "TimeWidth/F");
  t_Output_clusteredhits->Branch("ENu",          &out_ENu,         "ENu/D");
  t_Output_clusteredhits->Branch("ENu_Lep",      &out_ENu_Lep,     "ENu_Lep/D");
  t_Output_clusteredhits->Branch("MarlTime",     &out_MarlTime,    "MarlTime/D");
  t_Output_clusteredhits->Branch("HitView",      &out_HitView);
  t_Output_clusteredhits->Branch("GenType",      &out_GenType);
  t_Output_clusteredhits->Branch("HitChan",      &out_HitChan);
  t_Output_clusteredhits->Branch("HitTime",      &out_HitTime);
  t_Output_clusteredhits->Branch("HitSADC",      &out_HitSADC);
  t_Output_clusteredhits->Branch("HitRMS",       &out_HitRMS);
  h_ENu_MC      = tfs->make<TH1D>("h_ENu_MC","h_ENu_MC",35,0,50);
  h_MarlTime_MC = tfs->make<TH1D>("h_MarlTime_MC","h_MarlTime_MC",100,-0.1,10.5);
  h_TimeElapsed = tfs->make<TH1D>("h_TimeElapsed", "h_TimeElapsed", 50,0,0.5);
  h_ENu_MC->Sumw2();
  std::cout << "Finished beginJob" << std::endl;
}

void DAQQuickClustering::clusterChannels ( std::vector< recoHit > &  vec_Hits, std::vector< ClusterHitsInTime > &  vec_ChannelCluster, unsigned int const &  config  )

private

Definition at line 193 of file DAQQuickClustering_module.cc.

{
  //HERE IT IS ASSUMED THAT THE HITS APPEAR SEQUENTIALLY BY CHANNEL.

  for(unsigned int i = 0; i < vec_Hits.size()-1; i++)
  {
    std::vector<recoHit> vec_TempHits;
    double width = std::abs(vec_Hits.at(i  ).getHitChan()-
                            vec_Hits.at(i+1).getHitChan());
    //std::cout << cut_AdjChanTolerance.at(config) << std::endl;
    if(width <= cut_AdjChanTolerance.at(config))
    {
      int  channelCount = 1;
      
      vec_TempHits.push_back(vec_Hits.at(i  ));
      vec_TempHits.push_back(vec_Hits.at(i+1));

      while((i+channelCount+1)<vec_Hits.size() &&
            std::abs(vec_Hits.at(i+channelCount).getHitChan()-vec_Hits.at(i+channelCount+1).getHitChan())<=cut_AdjChanTolerance.at(config))
      {
        vec_TempHits.push_back(vec_Hits.at(i + channelCount + 1));
        channelCount++;
      }

      i = i + channelCount;
      ClusterHitsInTime temp(cut_TimeWindowSize[config]);
      temp.DoIt(vec_TempHits);
      
      vec_ChannelCluster.push_back(temp);
      //std::cout << " vec_ChannelCluster.size() " << vec_ChannelCluster.size() << std::endl;
            
    }
  }

  return;
}

void DAQQuickClustering::clusterCut ( std::vector< cluster > &  vec_Clusters, unsigned int const &  config  )

private

Definition at line 234 of file DAQQuickClustering_module.cc.

{
  //REMEMBER WE NEED BOTH A MAXIMUM AND MINIMUM CHANNEL WIDTH DUE TO COSMICS.

  for(std::vector<cluster>::iterator it_Clusters=vec_Clusters.begin();
      it_Clusters!=vec_Clusters.end();)
  {
    bool fail_MinChannels(false), fail_ChanWidth(false), fail_TotalADC(false);
    if(it_Clusters->getNChan()<cut_MinChannels.at(config))
      fail_MinChannels = true;

    if(it_Clusters->getChanWidth()<cut_MinChanWidth.at(config))
      fail_ChanWidth = true;

    if(it_Clusters->getHitSADC()<cut_TotalADC.at(config))
      fail_TotalADC = true;

    if(fail_MinChannels == true || fail_ChanWidth == true || fail_TotalADC == true)
      it_Clusters = vec_Clusters.erase(it_Clusters);
    else
      it_Clusters++;
  }
  return;
}

void DAQQuickClustering::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 384 of file DAQQuickClustering_module.cc.

{
  std::cout << "Job ended." << std::endl; 
  std::cerr << "firstCatch  " << firstCatch  << std::endl; 
  std::cerr << "secondCatch " << secondCatch << std::endl; 
  std::cerr << "thirdCatch  " << thirdCatch  << std::endl; 
}

void DAQQuickClustering::FillMyMaps ( std::map< int, simb::MCParticle > &  MyMap, art::FindManyP< simb::MCParticle >  Assn, art::ValidHandle< std::vector< simb::MCTruth > >  Hand  )

private

Definition at line 393 of file DAQQuickClustering_module.cc.

{
  for ( size_t L1=0; L1 < Hand->size(); ++L1 ) {
    for ( size_t L2=0; L2 < Assn.at(L1).size(); ++L2 ) {
      const simb::MCParticle ThisPar = (*Assn.at(L1).at(L2));
      MyMap[ThisPar.TrackId()] = ThisPar;
    }
  }
  return;
}

bool DAQQuickClustering::InMyMap ( int  TrID, std::map< int, simb::MCParticle >  ParMap  )

private

Definition at line 371 of file DAQQuickClustering_module.cc.

{
     
  std::map<int, simb::MCParticle>::iterator ParIt;
  ParIt = ParMap.find(TrID);

  if(ParIt != ParMap.end()) return true;
  else                       return false;

}

void DAQQuickClustering::makeConfigGraph ( )

private

Definition at line 273 of file DAQQuickClustering_module.cc.

{
  art::ServiceHandle<art::TFileService> tfs;

  for(unsigned int i = 0; i < NConfigs; i++)
  {
    TGraph* g_config = tfs->make<TGraph>(NCuts);
    g_config->SetName(Form("g_ConfigDefinitions%i",i));
    g_config->SetPoint(0, 1, cut_AdjChanTolerance[i]);
    g_config->SetPoint(1, 2, cut_HitsInWindow    [i]);
    g_config->SetPoint(2, 3, cut_MinChannels     [i]);
    g_config->SetPoint(3, 4, cut_MinChanWidth    [i]);
    g_config->SetPoint(4, 5, cut_TimeWindowSize  [i]);
    g_config->SetPoint(5, 6, cut_TotalADC        [i]);
  }
 
}

DAQQuickClustering& DAQQuickClustering::operator= ( DAQQuickClustering const &  )

delete

DAQQuickClustering& DAQQuickClustering::operator= ( DAQQuickClustering &&  )

delete

void DAQQuickClustering::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 105 of file DAQQuickClustering_module.cc.

{

  fRawDigitLabel = p.get<std::string>("RawDigitLabel");
  fHitLabel      = p.get<std::string>("HitLabel"     );
  
  fGEANTLabel    = p.get<std::string>("GEANT4Label"  );
  fMARLLabel     = p.get<std::string>("MARLEYLabel"  );
  fAPALabel      = p.get<std::string>("APALabel"     );
  fCPALabel      = p.get<std::string>("CPALabel"     );
  fAr39Label     = p.get<std::string>("Argon39Label" );
  fNeutLabel     = p.get<std::string>("NeutronLabel" );
  fKrypLabel     = p.get<std::string>("KryptonLabel" );
  fPlonLabel     = p.get<std::string>("PoloniumLabel");
  fRdonLabel     = p.get<std::string>("RadonLabel"   );
  fAr42Label     = p.get<std::string>("Argon42Label" );
                                      

  cut_AdjChanTolerance = p.get<std::vector<int>>  ("AdjChanTolerance");
  cut_HitsInWindow     = p.get<std::vector<int>>  ("HitsInWindow");
  cut_MinChannels      = p.get<std::vector<int>>  ("MinChannels");
  cut_MinChanWidth     = p.get<std::vector<int>>  ("MinChanWidth");
  cut_TimeWindowSize   = p.get<std::vector<float>>("TimeWindowSize");
  cut_TotalADC         = p.get<std::vector<float>>("TotalADC");

  NConfigs = cut_AdjChanTolerance.size();
  NCuts    = 6;

  detectorScaling = p.get<double>("detectorScaling");

} // Reconfigure

void DAQQuickClustering::ResetVariables ( )

private

Definition at line 139 of file DAQQuickClustering_module.cc.

{
  MarlParts.clear(); APAParts .clear(); CPAParts .clear(); Ar39Parts.clear();
  NeutParts.clear(); KrypParts.clear(); PlonParts.clear(); RdonParts.clear();
  Ar42Parts.clear();
  trkIDToPType.clear();
  Run = SubRun = Event = -1;

  TotGen_Marl = TotGen_APA  = TotGen_CPA  = TotGen_Ar39 = 0;
  TotGen_Neut = TotGen_Kryp = TotGen_Plon = TotGen_Rdon = 0;
  TotGen_Ar42 = 0;

  out_HitView.clear();
  out_GenType.clear();
  out_HitChan.clear();
  out_HitTime.clear();
  out_HitSADC.clear();
  out_HitRMS .clear(); 
   
  NTotHits = NColHits = NIndHits = 0; 
  nHitsNotBackTracked = 0;
  for (int hh=0; hh<nMaxHits; ++hh)
  {
    HitView[hh] = HitSize[hh] = HitChan[hh] = GenType[hh] = 0;
    HitTime[hh] = HitRMS [hh] = HitSADC[hh] = 0;
    HitInt [hh] = HitPeak[hh] = HitTPC [hh] = 0;
    NCorrespondingIDEs[hh] = 0;
    Hit_X[hh] = 0;                  
    Hit_Y[hh] = 0;                  
    Hit_Z[hh] = 0;                  
    Hit_Energy[hh] = 0;             
    Hit_NumElectrons[hh] = 0;
  } 

  MarlSample.clear();
  MarlTime  .clear();
  MarlWeight.clear();
  ENu = 0;
  ENu_Lep = 0;
  VertX = 0;
  VertY = 0;
  VertZ = 0;
  VertexT = 0;
  Px    = 0;
  Py    = 0;
  Pz    = 0;
  VertexChan = 0;
  map_EventToMC.clear();
   

}

void DAQQuickClustering::trigger ( std::vector< cluster > &  vec_Clusters, unsigned int const &  config  )

private

Definition at line 261 of file DAQQuickClustering_module.cc.

{
  
  for(unsigned int i = 0; i < vec_Clusters.size(); i++)
    if(vec_Clusters.at(i).getNHits()>=cut_HitsInWindow.at(config))
      vec_Clusters.at(i).setTriggerFlag(1);

  return;
}

PType DAQQuickClustering::WhichParType ( int  TrID )

private

Definition at line 359 of file DAQQuickClustering_module.cc.

{
  PType ThisPType = kUnknown;
  auto const& it=trkIDToPType.find(TrID);
  if(it!=trkIDToPType.end()){
    ThisPType=it->second;
  }
  
  return ThisPType;
}

Member Data Documentation

std::map< int, simb::MCParticle > DAQQuickClustering::APAParts

private

Definition at line 240 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::Ar39Parts

private

Definition at line 242 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::Ar42Parts

private

Definition at line 247 of file DAQQuickClustering_module.h.

art::ServiceHandle<cheat::BackTrackerService> DAQQuickClustering::bt_serv

private

Definition at line 346 of file DAQQuickClustering_module.h.

int DAQQuickClustering::CCNC

private

Definition at line 329 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::CPAParts

private

Definition at line 241 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::cut_AdjChanTolerance

private

Definition at line 225 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::cut_HitsInWindow

private

Definition at line 226 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::cut_MinChannels

private

Definition at line 227 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::cut_MinChanWidth

private

Definition at line 228 of file DAQQuickClustering_module.h.

std::vector<float> DAQQuickClustering::cut_TimeWindowSize

private

Definition at line 229 of file DAQQuickClustering_module.h.

std::vector<float> DAQQuickClustering::cut_TotalADC

private

Definition at line 230 of file DAQQuickClustering_module.h.

double DAQQuickClustering::detectorScaling

private

Definition at line 255 of file DAQQuickClustering_module.h.

double DAQQuickClustering::ENu

private

Definition at line 335 of file DAQQuickClustering_module.h.

double DAQQuickClustering::ENu_Lep

private

Definition at line 336 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Event

private

Definition at line 259 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fAPALabel

private

Definition at line 240 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fAr39Label

private

Definition at line 242 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fAr42Label

private

Definition at line 247 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fCPALabel

private

Definition at line 241 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fGEANTLabel

private

Definition at line 238 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fHitLabel

private

Definition at line 237 of file DAQQuickClustering_module.h.

int DAQQuickClustering::firstCatch

private

Definition at line 232 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fKrypLabel

private

Definition at line 244 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fMARLLabel

private

Definition at line 239 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fNeutLabel

private

Definition at line 243 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fPlonLabel

private

Definition at line 245 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fRawDigitLabel

private

Definition at line 236 of file DAQQuickClustering_module.h.

std::string DAQQuickClustering::fRdonLabel

private

Definition at line 246 of file DAQQuickClustering_module.h.

int DAQQuickClustering::GenType[nMaxHits]

private

Definition at line 308 of file DAQQuickClustering_module.h.

art::ServiceHandle<geo::Geometry> DAQQuickClustering::geo

private

Definition at line 345 of file DAQQuickClustering_module.h.

TH1D* DAQQuickClustering::h_ENu_MC

private

Definition at line 268 of file DAQQuickClustering_module.h.

TH1D* DAQQuickClustering::h_MarlTime_MC

private

Definition at line 269 of file DAQQuickClustering_module.h.

TH1D* DAQQuickClustering::h_TimeElapsed

private

Definition at line 270 of file DAQQuickClustering_module.h.

float DAQQuickClustering::Hit_Energy[nMaxHits]

private

Definition at line 312 of file DAQQuickClustering_module.h.

float DAQQuickClustering::Hit_NumElectrons[nMaxHits]

private

Definition at line 313 of file DAQQuickClustering_module.h.

float DAQQuickClustering::Hit_X[nMaxHits]

private

Definition at line 309 of file DAQQuickClustering_module.h.

float DAQQuickClustering::Hit_Y[nMaxHits]

private

Definition at line 310 of file DAQQuickClustering_module.h.

float DAQQuickClustering::Hit_Z[nMaxHits]

private

Definition at line 311 of file DAQQuickClustering_module.h.

int DAQQuickClustering::HitChan[nMaxHits]

private

Definition at line 302 of file DAQQuickClustering_module.h.

float DAQQuickClustering::HitInt[nMaxHits]

private

Definition at line 306 of file DAQQuickClustering_module.h.

int DAQQuickClustering::HitNucleon

private

Definition at line 330 of file DAQQuickClustering_module.h.

float DAQQuickClustering::HitPeak[nMaxHits]

private

Definition at line 307 of file DAQQuickClustering_module.h.

float DAQQuickClustering::HitRMS[nMaxHits]

private

Definition at line 304 of file DAQQuickClustering_module.h.

float DAQQuickClustering::HitSADC[nMaxHits]

private

Definition at line 305 of file DAQQuickClustering_module.h.

int DAQQuickClustering::HitSize[nMaxHits]

private

Definition at line 300 of file DAQQuickClustering_module.h.

float DAQQuickClustering::HitTime[nMaxHits]

private

Definition at line 303 of file DAQQuickClustering_module.h.

int DAQQuickClustering::HitTPC[nMaxHits]

private

Definition at line 301 of file DAQQuickClustering_module.h.

int DAQQuickClustering::HitView[nMaxHits]

private

Definition at line 299 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::KrypParts

private

Definition at line 244 of file DAQQuickClustering_module.h.

std::map<int,std::vector<double> > DAQQuickClustering::map_EventToMC

private

Definition at line 261 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::MarlParts

private

Definition at line 239 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::MarlSample

private

Definition at line 332 of file DAQQuickClustering_module.h.

std::vector<double> DAQQuickClustering::MarlTime

private

Definition at line 333 of file DAQQuickClustering_module.h.

std::vector<double> DAQQuickClustering::MarlWeight

private

Definition at line 334 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Mode

private

Definition at line 328 of file DAQQuickClustering_module.h.

int DAQQuickClustering::NColHits

private

Definition at line 296 of file DAQQuickClustering_module.h.

unsigned int DAQQuickClustering::NConfigs

private

Definition at line 252 of file DAQQuickClustering_module.h.

int DAQQuickClustering::NCorrespondingIDEs[nMaxHits]

private

Definition at line 314 of file DAQQuickClustering_module.h.

unsigned int DAQQuickClustering::NCuts

private

Definition at line 253 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::NeutParts

private

Definition at line 243 of file DAQQuickClustering_module.h.

int DAQQuickClustering::nHitsNotBackTracked

private

Definition at line 298 of file DAQQuickClustering_module.h.

int DAQQuickClustering::NIndHits

private

Definition at line 297 of file DAQQuickClustering_module.h.

int DAQQuickClustering::NTotHits

private

Definition at line 295 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Nu_Lep_Type

private

Definition at line 327 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Nu_Type

private

Definition at line 326 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_ChanWidth

private

Definition at line 277 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_Cluster

private

Definition at line 272 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_Config

private

Definition at line 274 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_EndChan

private

Definition at line 276 of file DAQQuickClustering_module.h.

double DAQQuickClustering::out_ENu

private

Definition at line 285 of file DAQQuickClustering_module.h.

double DAQQuickClustering::out_ENu_Lep

private

Definition at line 286 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_Event

private

Definition at line 273 of file DAQQuickClustering_module.h.

float DAQQuickClustering::out_FirstTimeHit

private

Definition at line 282 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::out_GenType

private

Definition at line 289 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::out_HitChan

private

Definition at line 290 of file DAQQuickClustering_module.h.

std::vector<double> DAQQuickClustering::out_HitRMS

private

Definition at line 293 of file DAQQuickClustering_module.h.

std::vector<double> DAQQuickClustering::out_HitSADC

private

Definition at line 292 of file DAQQuickClustering_module.h.

std::vector<double> DAQQuickClustering::out_HitTime

private

Definition at line 291 of file DAQQuickClustering_module.h.

std::vector<int> DAQQuickClustering::out_HitView

private

Definition at line 288 of file DAQQuickClustering_module.h.

float DAQQuickClustering::out_LastTimeHit

private

Definition at line 283 of file DAQQuickClustering_module.h.

double DAQQuickClustering::out_MarlTime

private

Definition at line 287 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_NChan

private

Definition at line 278 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_NHits

private

Definition at line 280 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_StartChan

private

Definition at line 275 of file DAQQuickClustering_module.h.

float DAQQuickClustering::out_SumADC

private

Definition at line 281 of file DAQQuickClustering_module.h.

float DAQQuickClustering::out_TimeWidth

private

Definition at line 284 of file DAQQuickClustering_module.h.

int DAQQuickClustering::out_Type

private

Definition at line 279 of file DAQQuickClustering_module.h.

art::ServiceHandle<cheat::ParticleInventoryService> DAQQuickClustering::pi_serv

private

Definition at line 347 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::PlonParts

private

Definition at line 245 of file DAQQuickClustering_module.h.

double DAQQuickClustering::Px

private

Definition at line 341 of file DAQQuickClustering_module.h.

double DAQQuickClustering::Py

private

Definition at line 342 of file DAQQuickClustering_module.h.

double DAQQuickClustering::Pz

private

Definition at line 343 of file DAQQuickClustering_module.h.

std::map< int, simb::MCParticle > DAQQuickClustering::RdonParts

private

Definition at line 246 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Run

private

Definition at line 257 of file DAQQuickClustering_module.h.

int DAQQuickClustering::secondCatch

private

Definition at line 233 of file DAQQuickClustering_module.h.

int DAQQuickClustering::SubRun

private

Definition at line 258 of file DAQQuickClustering_module.h.

TTree* DAQQuickClustering::t_Output_clusteredhits

private

Definition at line 266 of file DAQQuickClustering_module.h.

TTree* DAQQuickClustering::t_Output_unusedhits

private

Definition at line 265 of file DAQQuickClustering_module.h.

int DAQQuickClustering::Target

private

Definition at line 331 of file DAQQuickClustering_module.h.

int DAQQuickClustering::thirdCatch

private

Definition at line 234 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_APA

private

Definition at line 316 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Ar39

private

Definition at line 318 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Ar42

private

Definition at line 323 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_CPA

private

Definition at line 317 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Kryp

private

Definition at line 320 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Marl

private

Definition at line 315 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Neut

private

Definition at line 319 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Plon

private

Definition at line 321 of file DAQQuickClustering_module.h.

int DAQQuickClustering::TotGen_Rdon

private

Definition at line 322 of file DAQQuickClustering_module.h.

std::map<int, PType> DAQQuickClustering::trkIDToPType

private

Definition at line 250 of file DAQQuickClustering_module.h.

int DAQQuickClustering::VertexChan

private

Definition at line 325 of file DAQQuickClustering_module.h.

double DAQQuickClustering::VertexT

private

Definition at line 340 of file DAQQuickClustering_module.h.

double DAQQuickClustering::VertX

private

Definition at line 337 of file DAQQuickClustering_module.h.

double DAQQuickClustering::VertY

private

Definition at line 338 of file DAQQuickClustering_module.h.

double DAQQuickClustering::VertZ

private

Definition at line 339 of file DAQQuickClustering_module.h.

---

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

• duneana/duneana/DAQSimAna/DAQQuickClustering_module.h
• duneana/duneana/DAQSimAna/DAQQuickClustering_module.cc