pdunedp::Purity Class Reference

Inheritance diagram for pdunedp::Purity:

Classes
struct	Config

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	Purity (Parameters const &config)
	Purity (Purity const &)=delete
	Purity (Purity &&)=delete
Purity &	operator= (Purity const &)=delete
Purity &	operator= (Purity &&)=delete
void	analyze (art::Event const &e) override
void	MakeDataProduct ()
void	beginJob () override
void	endJob () override
void	Clear ()
void	FillEventHitsTree (std::vector< recob::Hit > hits)
double	GetCharge (std::vector< recob::Hit > hits)
double	GetCharge (const std::vector< art::Ptr< recob::Hit > > hits)
bool	IsCrossing (TVector3 Start, TVector3 End)
bool	StitchTracks (recob::Track Track1, recob::Track Track2, TVector3 &Edge1, TVector3 &Edge2)
void	Make_dEdx (std::vector< double > &dEdx, std::vector< double > &range, const std::vector< pdunedp::bHitInfo > &hits, recob::Track mip, int Flip, unsigned int plane)
void	FillTajectoryGraph (std::map< size_t, recob::Track > MipCandidate, art::FindManyP< recob::Hit > HitsTrk)
void	FillPurityHist (recob::Track track, std::vector< art::Ptr< recob::Hit >> hits)
void	FindMipInfo (recob::Track mip, std::vector< art::Ptr< recob::Hit >> vhits, const std::vector< const recob::TrackHitMeta *, std::allocator< const recob::TrackHitMeta * > > vmeta)
double	GetCorrectedCharge (recob::Track trk, double charge, unsigned int plane)
double	GetCorrection (recob::Track trk, unsigned int plane)
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
calo::CalorimetryAlg	fCalorimetryAlg
art::ServiceHandle< geo::Geometry >	geom
art::InputTag	fCalWireModuleLabel
art::InputTag	fHitModuleLabel
art::InputTag	fClusterModuleLabel
art::InputTag	fTrackModuleLabel
double	fTotalGain
double	fECut
double	fLength
double	fStitchAngle
double	fStitchDistance
std::vector< double >	fVolCut
int	fNumOfBins
double	fADCtoCharge
double	fMinDx
int	Events =0
int	SkipEvents =0
int	BadEvents =0
int	fRun
int	fSubRun
int	fEvent
int	fNtotHits
int	fNtotTracks
int	fNHitsTrack
int	fTrkNum
int	fNumOfMips
int	fMipNum
int	fMipIndex
int	fPlane
int	fIndex
int	fStartTick
int	fEndTick
double	fGodnessOfFit
double	fChargeIntegral
double	fPeakAmplitude
double	fTrackLength
double	fChi2Ndof
double	fHitsCharge
double	fTrackCharge
double	fHitsTrackCharge
double	fMipCharge
double	fMipLength
double	fMipPhi
double	fMipTheta
int	fNHitsMip
double	fCosX
double	fCosY
double	fCosZ
double	fSpacePointX
double	fSpacePointY
double	fSpacePointZ
double	fSummedADC
double	fSummedADC_corrected
double	fADCIntegral_corrected
double	fADCIntegral
double	fADCPeak
double	fADCPeak_corrected
double	fCorrectedCharge
double	fCharge
double	fDrift
int	fWire
double	fCorrection
double	fdQds
double	fdEdx
double	fRange
int	fBin
double	fSummedCharge
int	fEntries
double	fADCToElectrons
double	fElectronCharge = 1.60217662e-4
double	fMipChargeCm = 10
std::map< size_t, std::vector< pdunedp::bHitInfo >[3] >	fTrk2InfoMap
std::map< size_t, recob::Track >	TrackList
std::map< int, int >	goodevents
TTree *	fTree
TTree *	fTreeTrk
TTree *	fTreeMip
TTree *	fTreeHitsMip
TTree *	fTreeCalib
TTree *	fTreeHits
TTree *	fTreePurity
TTree *	fTreePurityMean
TH1D *	htbin [3][100]
TH1D *	htbin_singlehits [3][100]
TH1D *	htbin_num [3][100]
TH2D *	hTrkTrajectory_0
TH2D *	hTrkTrajectory_1

Additional Inherited Members
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 76 of file Purity_module.cc.

Member Typedef Documentation

using pdunedp::Purity::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 139 of file Purity_module.cc.

Constructor & Destructor Documentation

pdunedp::Purity::Purity ( Parameters const &  config )

explicit

Definition at line 232 of file Purity_module.cc.

                                               : EDAnalyzer(config),
  fCalorimetryAlg(config().CalorimetryAlg()),
  fCalWireModuleLabel(config().CalWireModuleLabel()),
  fHitModuleLabel(config().HitModuleLabel()),
  fClusterModuleLabel(config().ClusterModuleLabel()),
  fTrackModuleLabel(config().TrackModuleLabel()),
  fTotalGain(config().TotalGain()),
  fECut(config().EnergyCut()),
  fLength(config().Length()),
  fStitchAngle(config().StitchAngle()),
  fStitchDistance(config().StitchDistance()),
  fVolCut(config().VolCut()),
  fNumOfBins(config().NumOfBins()),
  fADCtoCharge(config().ADCtoCharge()),
  fMinDx(config().MinDx())
{
  this->MakeDataProduct();
}

pdunedp::Purity::Purity ( Purity const &  )

delete

pdunedp::Purity::Purity ( Purity &&  )

delete

Member Function Documentation

void pdunedp::Purity::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 389 of file Purity_module.cc.

                                             {
  fRun = e.run();
  fSubRun = e.subRun();
  fEvent = e.id().event();   Events++;

  Clear();

  //Require data product handles
  auto CalwireHandle = e.getValidHandle<std::vector<recob::Wire> >(fCalWireModuleLabel);
  auto HitHandle = e.getValidHandle<std::vector<recob::Hit> >(fHitModuleLabel);
  auto ClusterHandle = e.getValidHandle<std::vector<recob::Cluster> >(fClusterModuleLabel);
  auto TrackHandle = e.getValidHandle<std::vector<recob::Track> >(fTrackModuleLabel);

  //check data products exists
  if(!CalwireHandle || !HitHandle || !ClusterHandle || !TrackHandle){
    BadEvents++;
    return;
  }

  //first selection of data based on ADC counts (separate energetic showers)
  float SumWireADC=0;
  for(auto const& Calwire : *CalwireHandle){
    for(float ADC : Calwire.Signal()){
      SumWireADC += ADC;
    }
  }
  float Edep = SumWireADC*fADCToElectrons*fElectronCharge;
  float Efrac = Edep/(fMipChargeCm*350*fTotalGain);
  mf::LogVerbatim("pdunedp::Purity")<< "Energy: " << Edep << " " << Efrac;
  if(Efrac > fECut){
    mf::LogVerbatim("pdunedp::Purity") << "Too much energy energy deposited! Not a mip";
    SkipEvents++;
    return;
  }

  //Get Hits charge
  fNtotHits = HitHandle->size();
  fHitsCharge = GetCharge(*HitHandle);

  FillEventHitsTree(*HitHandle);

  fNtotTracks = (int)TrackHandle->size();
  if(fNtotTracks == 0){
    mf::LogVerbatim("pdunedp::Purity") << "Event has no 3D tracks";
    SkipEvents++;
    return;
  }

  art::FindManyP< recob::Hit > HitsFromTrack(TrackHandle, e, fTrackModuleLabel);
  art::FindManyP< recob::Hit, recob::TrackHitMeta > TrackHitMeta(TrackHandle, e, fTrackModuleLabel);

  double ChargeTrk=0;
  //int skippedTrk=0;
  for(size_t t=0; t<(size_t)fNtotTracks; t++){
    //retrive information for every track (not mip selected yet)
    auto track = TrackHandle->at(t);
    TrackList[t] = track;
    fTrackLength = track.Length();
    fChi2Ndof = track.Chi2PerNdof();
    fNHitsTrack = HitsFromTrack.at(t).size();
    fTrkNum = (int)t;
    fTrackCharge = GetCharge(HitsFromTrack.at(t));
    ChargeTrk+= GetCharge(HitsFromTrack.at(t));
    fTreeTrk->Fill();
  }
  fHitsTrackCharge = ChargeTrk; //<-- Sum charge from all event

  std::map<size_t, recob::Track > fMipCandidate;
  /*Try to select mips. Particle passing trougth the detector (Selecting a fiducial
  volume to keep into account the inefficiencies of the lems ).
  Tracks that doesn't meet the criteria of crossing the detector after a first pass
  are stitched togheter and the crossing criteria are evaluated again.
  Stitch is done only on two consecuitve tracks. TODO: Stitch more than two tracks*/

  std::map<size_t, recob::Track >::iterator It;
  for(It = TrackList.begin(); It !=TrackList.end(); It++){
    //first check if the track is already inside the candidate vector
    if (fMipCandidate.find(It->first) != fMipCandidate.end()){
      continue;
    }

    TVector3 Start = (It->second).Vertex<TVector3>(); TVector3 End = (It->second).End<TVector3>();
    if(!IsCrossing(Start, End) && (TrackList.size() > 1)){
      //It is not crossing, let's see if it can be stitched to something else
      for(std::map<size_t, recob::Track >::iterator It2=TrackList.begin(); It2!=TrackList.end(); It2++){
        if(It->first == It2->first){ continue; }
        TVector3 newStart; TVector3 newEnd;
        if(!StitchTracks(It->second, It2->second, newStart, newEnd)){
          continue; //Not stitching: continue searcing
        }else{
          if(IsCrossing(newStart, newEnd)){
            fMipCandidate[It->first] = It->second;
            fMipCandidate[It2->first] = It2->second;
          }else{
            break; //Stop the loop: we stitched and still not crossing...
          }
        }
      }//end while
    }else if(!IsCrossing(Start, End) && (TrackList.size() == 1)){
      continue;
    }else{
      fMipCandidate[It->first] = It->second;
    }
  }

  if(!fMipCandidate.size()){
    mf::LogVerbatim("pdunedp::Purity") << "No mips selected! ";
    SkipEvents++;
    return;
  }
  //print the list of mips:
  mf::LogVerbatim("pdunedp::Purity") << "List of mips candidates in event " << fEvent;
  size_t num=0;
  for(std::map<size_t, recob::Track >::iterator mip = fMipCandidate.begin(); mip !=fMipCandidate.end(); mip++){
    mf::LogVerbatim("pdunedp::Purity") << "Track: " << mip->first;
    fMipNum = (int)num++;
    fMipIndex = (int)mip->first;
    fMipCharge = GetCharge(HitsFromTrack.at(mip->first));
    fMipLength = (mip->second).Length();
    fMipPhi = (mip->second).Phi();      //polar angle
    fMipTheta = (mip->second).Theta();  //Azimutal angle
    fNHitsMip = (int)HitsFromTrack.at(mip->first).size();
    auto vhit = TrackHitMeta.at(mip->first);
    auto vmeta = TrackHitMeta.data(mip->first);
    FindMipInfo(mip->second, vhit, vmeta);
    FillPurityHist(mip->second, HitsFromTrack.at(mip->first));
    fTreeMip->Fill();
  }

  for (auto const & trkEntry : fTrk2InfoMap)
  {
    //check flip
    int Flip = 0; // tag track if it is reconstructed in the opposite direction
    if (fMipCandidate[trkEntry.first].End().Z() < fMipCandidate[trkEntry.first].Vertex().Z()){ Flip = 1; }

    for(unsigned int plane=0; plane < geom->Nplanes(0); plane++){ //TODO<<--Extend to different geometries
      std::vector< double > dEdx, range;
      auto const & info = trkEntry.second;
      Make_dEdx(dEdx, range, info[plane], fMipCandidate[trkEntry.first], Flip, plane);
      for (size_t i = 0; i < dEdx.size(); ++i){
        fdEdx = dEdx[i];
        fPlane = plane;
        fIndex = trkEntry.first;
        fRange = range[i];
        fTreeCalib->Fill();
      }
    }
  }

  FillTajectoryGraph(fMipCandidate, HitsFromTrack);
  goodevents[fSubRun]=fEvent;
  fTree->Fill();
}//end analyzer

void pdunedp::Purity::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 251 of file Purity_module.cc.

                            {
  fADCToElectrons = 1./dp->ElectronsToADC();
  //auto simChannelExtract = &*art::ServiceHandle<detsim::DPhaseSimChannelExtractService>();
  //fTotalGain = simChannelExtract->GainPerView()*2;
}

void pdunedp::Purity::Clear

(

)

Definition at line 1005 of file Purity_module.cc.

                         {
  TrackList.clear();
  fTrk2InfoMap.clear();
}

void pdunedp::Purity::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 1010 of file Purity_module.cc.

                          {
  mf::LogVerbatim("pdunedp::Purity") << "====== Run " << fRun << " summary ======";
  mf::LogVerbatim("pdunedp::Purity") << "Total Events: " << Events;
  mf::LogVerbatim("pdunedp::Purity") << "Bad Events: " << BadEvents;
  mf::LogVerbatim("pdunedp::Purity") << "Skipped Events: " << SkipEvents;
  mf::LogVerbatim("pdunedp::Purity") << "selected Events: " << Events-BadEvents-SkipEvents <<" \nList: ";
  for(std::map<int, int>::iterator goodevent = goodevents.begin(); goodevent != goodevents.end(); goodevent++){
      mf::LogVerbatim("pdunedp::Purity") << "SubRun: " << goodevent->first << " Event: " << goodevent->second ; //<----write this list on file (?)
  }
  mf::LogVerbatim("pdunedp::Purity") << "===========================";

}

void pdunedp::Purity::FillEventHitsTree

(

std::vector< recob::Hit >

hits

)

Definition at line 573 of file Purity_module.cc.

                                                             {
  //Fill a tree with additionals informations about the Event
  if(!hits.size()){ return;}

  for(auto hit : hits){
    unsigned short plane = hit.WireID().Plane;
    int wire = hit.WireID().Wire;
    double dqadc = hit.Integral();
    double tdrift = hit.PeakTime();
    if (!std::isnormal(dqadc) || (dqadc < 0)) continue;
    double dq = dqadc*fADCtoCharge;

    fPlane = plane;
    fStartTick = hit.StartTick();
    fEndTick = hit.EndTick();
    fADCIntegral = dqadc;
    fPeakAmplitude = hit.PeakAmplitude();
    fSummedADC = hit.SummedADC();
    fGodnessOfFit= hit.GoodnessOfFit();
    fChargeIntegral = dq;
    fDrift = tdrift;
    fWire = wire;
    fTreeHits->Fill();
  }
    return;
}

void pdunedp::Purity::FillPurityHist	(	recob::Track	track,
		std::vector< art::Ptr< recob::Hit >>	hits
	)

Definition at line 872 of file Purity_module.cc.

                                                                                    {
  /*Function intended to read the hits from a track and fill the histograms that can be
  used for further purity analysis*/
  mf::LogVerbatim("pdunedp::Purity") << "Start purity for track: " << fMipIndex;

  //determine the size (in ticks of each bin)
 for(int pl=0; pl<(int)geom->Nplanes(0); pl++){
  int TicksPerBin = detProp.NumberTimeSamples()/fNumOfBins;

  double charge[100]; double num[100];
  for (int nb=0; nb<fNumOfBins; nb++){ charge[nb]=0; num[nb] =0; } //init the charge bins to 0

  if(!hits.size()){
    mf::LogError("pdunedp::Purity") << "The track has no hit associated!";
    return;
  }
  int nfound =0; int nstart; int nstop;

  for (size_t nh=0; nh<hits.size(); nh++){

      double dqadc = hits[nh]->Integral(); //ADCxticks
      double htime = hits[nh]->PeakTime(); //ticks
      unsigned short plane = hits[nh]->WireID().Plane;
      if(plane != pl) {continue;}
      //double conv = fCalorimetryAlg.ElectronsFromADCArea(dqadc, plane)*fElectronCharge;
      double conv = fADCtoCharge;
      double dq = dqadc*conv;
      dq = GetCorrectedCharge(track, dq, plane);

      //dq = GetCorrectedCharge(track, dq, plane); //corrected charge wrt to the track angle

      if (nh==0){
        nstart=0;
        nstop=fNumOfBins;
      }
      if (nh>0){
        nstart=std::max(0,nfound-5);
        nstop =std::min(nfound+5,fNumOfBins);
      }

      for(int nb=nstart; nb<nstop; nb++) {
        double tmin=TicksPerBin*nb;
        double tmax=TicksPerBin*(nb+1);
        //mf::LogVerbatim("pdunedp::Purity") << "nb " << nb << " " << tmin << " " <<tmax;

        if ( (htime>=tmin) && (htime<tmax) ){
          htbin_singlehits[pl][nb]->Fill (dq);
          fBin = nb;
          fCharge = dq;
          fPlane = pl;
          fTreePurity->Fill();
          charge[nb]+=dq;
          num[nb]++;
          //mf::LogVerbatim("pdunedp::Purity") << "charge" << charge[nb];
          nfound=nb;
          //mf::LogVerbatim("pdunedp::Purity") << "nfound" << nb;

       }
     } //close loop on time bins
  } //close loop on hit on tracks

  // fill histograms for fit
  //select first and last time interval with charge deposition >0
  int ilast=-1000;
  int ifirst=1000;

  for (int nb=0; nb<fNumOfBins; nb++) {
    if (charge[nb]>0 && num[nb] >0) {
      //mf::LogVerbatim("pdunedp::Purity") << "nb " << nb << " Norm: " << charge[nb];
          if (nb<ifirst) ifirst=nb;
                if (nb>ilast)  ilast=nb;
        }
  }

  //mf::LogVerbatim("pdunedp::Purity") << "first: " << ifirst << " last " << ilast;

  for (int nb=0; nb<ilast-ifirst-1; nb++){
      mf::LogVerbatim("pdunedp::Purity") << "plane " << pl << " nb " << nb+ifirst+1<< " charge bin " << charge[nb+ifirst+1]/num[nb+ifirst+1];
      //htbin_norm[nb]->Fill(charge[nb+ifirst+1]/charge[ifirst+1]); //fill the bin with charge normalized to first bin
      htbin[pl][nb]->Fill( charge[nb+ifirst+1]/num[nb+ifirst+1] );
      htbin_num[pl][nb]->Fill( num[nb+ifirst+1] );
  }

  /*
  for (int nb=0; nb<ilast-ifirst; nb++){
      //mf::LogVerbatim("pdunedp::Purity") << "nb " << nb+ifirst+1<< " Norm: " << charge[ifirst+1] << " charge bin " << charge[nb+ifirst+1];
      htbin_corr[nb]->Fill( charge[nb+ifirst]*angle_corr ); //fill the bin with charge normalized to first bin
  }
  */
 }
 return;
}

void pdunedp::Purity::FillTajectoryGraph	(	std::map< size_t, recob::Track >	MipCandidate,
		art::FindManyP< recob::Hit >	HitsTrk
	)

Definition at line 746 of file Purity_module.cc.

                                                                          {
  //Merge the 2D hit position view in a single graph to evaluate defects in mip selection
  //<<--TODO Generic geometry
  std::map<size_t, recob::Track >::iterator It;
  for(It = MipCandidate.begin(); It !=MipCandidate.end(); It++){
    auto hits = HitsTrk.at(It->first);
    if(!hits.size()){ continue; }
    for(auto const  hit : hits){
      unsigned short plane = hit->WireID().Plane;
      switch (plane) {
        case 0:
          hTrkTrajectory_0->Fill(hit->WireID().Wire, hit->PeakTime());
          break;
        case 1:
          hTrkTrajectory_1->Fill(hit->WireID().Wire, hit->PeakTime());
          break;
        default:
          mf::LogError("pdunedp::Purity") << "Invalid Plane;";
          break;
      }
    }
  }
  return;
}

void pdunedp::Purity::FindMipInfo	(	recob::Track	mip,
		std::vector< art::Ptr< recob::Hit >>	vhits,
		const std::vector< const recob::TrackHitMeta *, std::allocator< const recob::TrackHitMeta * > >	vmeta
	)

Definition at line 772 of file Purity_module.cc.

                                                                                           {
  /*This function is intended to caluclate the most important quantities from a mip
  and fill a tree for further analysis. There will be an entry for every hit in the mip.
  Stitched mips are considered independently. T0 is assumed to be 0, as consequence of
  the selection done*/

  if(!vhits.size()){return;}
  for(size_t h =0; h< vhits.size(); h++){
    unsigned int plane= vhits[h]->WireID().Plane;
    int wire= vhits[h]->WireID().Wire;
    //double pitch = geom->WirePitch(vhits[h]->View());
    double dQadc = vhits[h]->Integral();
    if (!std::isnormal(dQadc) || (dQadc < 0)) continue;
    size_t idx = vmeta[h]->Index();
    double dx =  vmeta[h]->Dx();

    auto TrajPoint3D = mip.TrajectoryPoint(idx);

    double dQ = dQadc*fADCtoCharge;
    fTrk2InfoMap[fMipIndex][plane].emplace_back(idx, dx, dQ, wire);

    TVector3 Dir = mip.Vertex<TVector3>() - mip.End<TVector3>();
    //Direction cosines of the track
    double CosX = 1./Dir.Unit().X();
    double CosY = 1./Dir.Unit().Y();
    double CosZ = 1./Dir.Unit().Z();

    double dQ_corr = GetCorrectedCharge(mip, dQ, plane);
    double PeakTime = vhits[h]->PeakTime();

    //Fill Purity tree
    fWire   = wire;
    fPlane  = plane;
    fCosX   = CosX;
    fCosY   = CosY;
    fCosZ   = CosZ;
    fADCIntegral_corrected = GetCorrectedCharge(mip, dQadc, plane); //corrected value
    fADCPeak = vhits[h]->PeakAmplitude();
    fADCPeak_corrected = GetCorrectedCharge(mip, vhits[h]->PeakAmplitude(), plane);//corrected value
    fSummedADC = vhits[h]->SummedADC();
    fSummedADC_corrected = GetCorrectedCharge(mip, fSummedADC, plane);//corrected value
    fADCIntegral    = dQadc;
    fCharge = dQ;
    fCorrectedCharge = dQ_corr;
    fCorrection = dQ_corr/dQ;
    fSpacePointX = TrajPoint3D.position.X();
    fSpacePointY = TrajPoint3D.position.Y();
    fSpacePointZ = TrajPoint3D.position.Z();
    if(dx>0){
        fdQds = dQ/dx;
    }
    fDrift = PeakTime;
    fStartTick = vhits[h]->StartTick();
    fEndTick = vhits[h]->EndTick();
    fTreeHitsMip->Fill();
  }
  return;
}

double pdunedp::Purity::GetCharge

(

std::vector< recob::Hit >

hits

)

Definition at line 558 of file Purity_module.cc.

                                                       {
  //It returns the uncalibrated charge in the detector given a list of hits (summed on both views)
  if(!hits.size()){ return 0.0;}

  double charge=0;
  for(auto hit : hits){
    //unsigned short plane = hit.WireID().Plane;
    double dqadc = hit.Integral();
    if (!std::isnormal(dqadc) || (dqadc < 0)) continue;
    //charge += dqadc*fCalorimetryAlg.ElectronsFromADCArea(dqadc, plane)*fElectronCharge;
    charge+= dqadc*fADCtoCharge;
  }
    return charge;
}

double pdunedp::Purity::GetCharge

(

const std::vector< art::Ptr< recob::Hit > >

hits

)

Definition at line 543 of file Purity_module.cc.

                                                                       {
  //It returns the uncalibrated charge in the detector given a list of hits (summed on both views)
  if(!hits.size()){ return 0.0;}

  double charge=0;
  for(auto const  hit : hits){
    //unsigned short plane = hit->WireID().Plane;
    double dqadc = hit->Integral();
    if (!std::isnormal(dqadc) || (dqadc < 0)) continue;
    //charge += dqadc*fCalorimetryAlg.ElectronsFromADCArea(dqadc, plane)*fElectronCharge;
    charge+= dqadc*fADCtoCharge;
  }
    return charge;
}

double pdunedp::Purity::GetCorrectedCharge	(	recob::Track	trk,
		double	charge,
		unsigned int	plane
	)

Definition at line 965 of file Purity_module.cc.

                                                                                         {
  /*Returns the corrected charge value corrected for the angle between the track direction and the wire pitch*/
  double angle =0.;
  double charge_corr =0.;

  switch (plane) {
    case 0:
      angle = atan( fabs( (trk.End().X()-trk.Vertex().X())/(trk.End().Y()-trk.Vertex().Y()) ) );
      break;
    case 1:
      angle = atan( fabs( (trk.End().X()-trk.Vertex().X())/(trk.End().Z()-trk.Vertex().Z()) ) );
      break;
    default:
      mf::LogError("pdunedp::Purity") << "Invalid plane number!";
      break;
  }
  if(fabs(cos(angle))>0.){ charge_corr = charge*fabs(cos(angle)); }
  return charge_corr;
}

double pdunedp::Purity::GetCorrection	(	recob::Track	trk,
		unsigned int	plane
	)

Definition at line 985 of file Purity_module.cc.

                                                                     {
  /*Returns the corrected charge value corrected for the angle between the track direction and the wire pitch*/
  double angle =0.;
  double correction =0.;

  switch (plane) {
    case 0:
      angle = atan( fabs( (trk.End().X()-trk.Vertex().X())/(trk.End().Y()-trk.Vertex().Y()) ) );
      break;
    case 1:
      angle = atan( fabs( (trk.End().X()-trk.Vertex().X())/(trk.End().Z()-trk.Vertex().Z()) ) );
      break;
    default:
      mf::LogError("pdunedp::Purity") << "Invalid plane number!";
      break;
  }
  if(fabs(cos(angle))>0.){ correction = fabs(cos(angle)); }
  return correction;
}

bool pdunedp::Purity::IsCrossing	(	TVector3	Start,
		TVector3	End
	)

Definition at line 654 of file Purity_module.cc.

                                                          {
  bool isCrossing = false;
  //define the boundaries (for the moment consider only 1 TPC) //<--TODO Generalisaztion to Multiple TPCs
  double vtx[3] = {0, 0, 150}; //<<--TODO Generalisaztion
  art::ServiceHandle<geo::Geometry> geom;
  geo::TPCID idtpc = geom->FindTPCAtPosition(vtx);
  if (geom->HasTPC(idtpc)) // <----Assuming there is only one TPC
  {
    /*
    const geo::TPCGeo& tpcgeo = geom->GetElement(idtpc);
    double minx = tpcgeo.MinX() + fVolCut[0]; double maxx = tpcgeo.MaxX() - fVolCut[1];
    double miny = tpcgeo.MinY() + fVolCut[2]; double maxy = tpcgeo.MaxY() - fVolCut[3];
    double minz = tpcgeo.MinZ() + fVolCut[4]; double maxz = tpcgeo.MaxZ() - fVolCut[5];

  //check if both start and end are both outside the fiducial volume (in at least one direction )
    if(   (   ((minx > Start.X()) || (maxx < Start.X()))
           || ((miny > Start.Y()) || (maxy < Start.Y()))
           || ((minz > Start.Z()) || (maxz < Start.Z())) ) //condiotion on start
        &&(   ((minx > End.X()  ) || (maxx < End.X()  ))
           || ((miny > End.Y()  ) || (maxy < End.Y()  ))
           || ((minz > End.Z()  ) || (maxz < End.Z()  )) ) //condiotion on end
      ){
        isCrossing=true;
      }else{
        isCrossing=false;
      }
  //from the point of view of X only particles that goes trougth the whole detector can be
  //considered particles on trigger (check all the possible cases)
    TVector3 Diff = Start-End;
    double ProjX = fabs(Diff.Unit().X());
    if( acos(1./ProjX) < 30*(TMath::Pi()/180) ){ //<<-- only on very vertical tracks
      if((maxx < Start.X())||(maxx < End.X())){
        if((minx < Start.X())||(minx < End.X()))
          isCrossing=false;
      }
      if((minx > Start.X())||(minx > End.X())){
        if((maxx > Start.X())||(maxx > End.X()))
          isCrossing=false;
      }
    }
  //Mips must be above certain length (for not be confused with random hadronic processes)
    if( (Diff).Mag() < fLength){
      isCrossing=false;
    }
    */

    //Cuts used for the purity analysis: vertical track, from anode to cathode contained within some fiducial volume in z
    //boundaries are minx -50 maxx +50
    //               miny -50 maxy +50
    //               minz 0   maxz +300
    const geo::TPCGeo& tpcgeo = geom->GetElement(idtpc);
    double minx = tpcgeo.MinX() + fVolCut[0]; double maxx = tpcgeo.MaxX() - fVolCut[1];
    double minz = tpcgeo.MinZ() + fVolCut[4]; double maxz = tpcgeo.MaxZ() - fVolCut[5];

    //first cut is applied on the length
    if( (End-Start).Mag() > fLength ){

      //check if the track is flipped on Z
      if( (End.Z() - Start.Z()) > 0 ){
        if( Start.Z() > minz && End.Z() < maxz){
          //check if the track is flipped on X
          if(End.X() > Start.X()){
            if((Start.X() <= minx && End.X() >= maxx)){
              isCrossing = true;
            }
          }else if(End.X() < Start.X()){
            if((End.X() <= minx && Start.X() >= maxx)){
              isCrossing = true;
            }
          }
        }
      }else if((End.Z() - Start.Z()) < 0){
        if( End.Z() > minz && Start.Z() < maxz){
          //check if the track is flipped on X
          if(End.X() > Start.X()){
            if((Start.X() <= minx && End.X() >= maxx)){
              isCrossing = true;
            }
          }else if(End.X() < Start.X()){
            if((End.X() <= minx && Start.X() >= maxx)){
              isCrossing = true;
            }
          }
        }
      }else{
        mf::LogError("pdunedp::Purity") << "Track starts and ends in the same Z coordinate";
      }
    }//end fLength
  }//end has tpc
  return isCrossing;
}

void pdunedp::Purity::Make_dEdx	(	std::vector< double > &	dEdx,
		std::vector< double > &	range,
		const std::vector< pdunedp::bHitInfo > &	hits,
		recob::Track	mip,
		int	Flip,
		unsigned int	plane
	)

Definition at line 832 of file Purity_module.cc.

                                                                                                                    {
  if (!hits.size()) return;

        dEdx.clear(); range.clear();

        double rmax = mip.Length();

        int i0 = hits.size() - 1; int i1 = -1; int di = -1;
        if (Flip) {i0 = 0; i1 = hits.size(); di = 1;}

        double de = 0.0;
        double dx = 0.0;
        double r0 = 0.0; double r1 = 0.0; double r = 0.0;

        double minDx = 0.1; // can be a parameter

        while ((i0 != i1) && (r < rmax))
        {
                dx = 0.0; de = 0.0;
                while ((i0 != i1) && (dx <= minDx))
                {
                        de += hits[i0].dE;
                        dx += hits[i0].dx;
                        i0 += di;
                }

                r0 = r1;
                r1 += dx;
                r = 0.5 * (r0 + r1);

                if ((de > 0.0) && (dx > 0.0) && (r < rmax))
                {
                        dEdx.push_back(de/dx);
                        range.push_back(r*GetCorrection(mip, plane));
                }
        }
  return;
}

void pdunedp::Purity::MakeDataProduct

(

)

Definition at line 257 of file Purity_module.cc.

                                   {

  //Tfile Services
  art::ServiceHandle<art::TFileService> tfs;

  //one entry for every event
  fTree = tfs->make<TTree>("Event", "LAr purity analysis information");
  fTree->Branch("fRun", &fRun,"fRun/I");
  fTree->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTree->Branch("fEvent", &fEvent, "fEvent/I");
  fTree->Branch("fNtotHits", &fNtotHits, "fNtotHits/I");
  fTree->Branch("fNtotTracks", &fNtotTracks, "fNtotTracks/I");
  fTree->Branch("fNumOfMips", &fNumOfMips, "fNumOfMips/I");
  fTree->Branch("fHitsCharge", &fHitsCharge, "fHitsCharge/D");
  fTree->Branch("fHitsTrackCharge", &fHitsTrackCharge, "fHitsTrackCharge/D");

  //branch for purity measuremtns usign mips: one entry per mip
  fTreeHits =tfs->make<TTree>("HitsEvent","Info on hits in event");
  fTreeHits->Branch("fRun", &fRun,"fRun/I");
  fTreeHits->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreeHits->Branch("fEvent", &fEvent, "fEvent/I");
  fTreeHits->Branch("fPlane", &fPlane, "fPlane/I");
  fTreeHits->Branch("fADCIntegral", &fADCIntegral, "fADCIntegral/D");
  fTreeHits->Branch("fChargeIntegral", &fChargeIntegral, "fChargeIntegral/D");
  fTreeHits->Branch("fPeakAmplitude", &fPeakAmplitude, "fPeakAmplitude/D");
  fTreeHits->Branch("fSummedADC", &fSummedADC, "fSummedADC/D");
  fTreeHits->Branch("fGodnessOfFit", &fGodnessOfFit, "fGodnessOfFit/D");
  fTreeHits->Branch("fDrift", &fDrift, "fDrift/D");
  fTreeHits->Branch("fWire", &fWire, "fWire/I");
  fTreeHits->Branch("fStartTick", &fStartTick, "fStartTick/I");
  fTreeHits->Branch("fEndTick", &fEndTick, "fEndTick/I");

  //one entry for every track
  fTreeTrk = tfs->make<TTree>("TrkInfo", "Information on tracks");
  fTreeTrk->Branch("fRun", &fRun,"fRun/I");
  fTreeTrk->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreeTrk->Branch("fEvent", &fEvent, "fEvent/I");
  fTreeTrk->Branch("fTrkNum", &fTrkNum, "fTrkNum/I");
  fTreeTrk->Branch("fChi2Ndof", &fChi2Ndof, "fChi2Ndof/D");
  fTreeTrk->Branch("fTrackLength", &fTrackLength, "fTrackLength/D");
  fTreeTrk->Branch("fTrackCharge", &fTrackCharge, "fTrackCharge/D");
  fTreeTrk->Branch("fNHitsTrack", &fNHitsTrack, "fNHitsTrack/I");

  //one entry for every selected mip (stitched mips are considered separately)
  fTreeMip = tfs->make<TTree>("MipInfo", "Information on mips");
  fTreeMip->Branch("fRun", &fRun,"fRun/I");
  fTreeMip->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreeMip->Branch("fEvent", &fEvent, "fEvent/I");
  fTreeMip->Branch("fMipIndex", &fMipIndex, "fMipIndex/I");
  fTreeMip->Branch("fMipNum", &fMipNum, "fMipNum/I");
  fTreeMip->Branch("fMipCharge", &fMipCharge, "fMipcharge/D");
  fTreeMip->Branch("fMipLength", &fMipLength, "fMipLength/D");
  fTreeMip->Branch("fMipLength", &fMipLength, "fMipLength/D");
  fTreeMip->Branch("fMipPhi", &fMipPhi, "fMipPhi/D");
  fTreeMip->Branch("fMipTheta", &fMipTheta, "fMipTheta/D");
  fTreeMip->Branch("fNHitsMip", &fNHitsMip, "fNHitsMip/I");

  //branch for purity measuremtns usign mips: one entry per mip
  fTreeHitsMip =tfs->make<TTree>("HitsMip","Info hits in mips");
  fTreeHitsMip->Branch("fRun", &fRun,"fRun/I");
  fTreeHitsMip->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreeHitsMip->Branch("fEvent", &fEvent, "fEvent/I");
  fTreeHitsMip->Branch("fMipIndex", &fMipIndex, "fMipIndex/I");
  fTreeHitsMip->Branch("fMipNum", &fMipNum, "fMipNum/I");
  fTreeHitsMip->Branch("fPlane", &fPlane, "fPlane/I");
  fTreeHitsMip->Branch("fCosX", &fCosX, "fCosX/D");
  fTreeHitsMip->Branch("fCosY", &fCosY, "fCosY/D");
  fTreeHitsMip->Branch("fCosZ", &fCosZ, "fCosZ/D");
  fTreeHitsMip->Branch("fADCIntegral", &fADCIntegral, "fADCIntegral/D");
  fTreeHitsMip->Branch("fADCIntegral_corrected", &fADCIntegral_corrected, "fADCIntegral_corrected/D");
  fTreeHitsMip->Branch("fADCPeak", &fADCPeak, "fADCPeak/D");
  fTreeHitsMip->Branch("fADCPeak_corrected", &fADCPeak_corrected, "fADCPeak_corrected/D");
  fTreeHitsMip->Branch("fCorrection", &fCorrection, "fCorrection/D");
  fTreeHitsMip->Branch("fSummedADC", &fSummedADC, "fSummedADC/D");
  fTreeHitsMip->Branch("fSummedADC_corrected", &fSummedADC_corrected, "fSummedADC_corrected/D");
  fTreeHitsMip->Branch("fCharge", &fCharge, "fCharge/D");
  fTreeHitsMip->Branch("fCorrectedCharge", &fCorrectedCharge, "fCorrectedCharge/D");
  fTreeHitsMip->Branch("fDrift", &fDrift, "fDrift/D");
  fTreeHitsMip->Branch("fStartTick", &fStartTick, "fStartTick/I");
  fTreeHitsMip->Branch("fEndTick", &fEndTick, "fEndTick/I");
  fTreeHitsMip->Branch("fWire", &fWire, "fWire/I");
  fTreeHitsMip->Branch("fSpacePointX", &fSpacePointX, "fSpacePointX/D");
  fTreeHitsMip->Branch("fSpacePointY", &fSpacePointY, "fSpacePointY/D");
  fTreeHitsMip->Branch("fSpacePointZ", &fSpacePointZ, "fSpacePointZ/D");
  fTreeHitsMip->Branch("fdQds", &fdQds, "fdQds/D");

  fTreeCalib = tfs->make<TTree>("Calibration", "dE/dx info");
  fTreeCalib->Branch("fRun", &fRun, "fRun/I");
  fTreeCalib->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreeCalib->Branch("fEvent", &fEvent, "fEvent/I");
  fTreeCalib->Branch("fIndex", &fIndex, "fIndex/I");//
  fTreeCalib->Branch("fdEdx", &fdEdx, "fdEdx/D");
  fTreeCalib->Branch("fRange", &fRange, "fRange/D");
  fTreeCalib->Branch("fPlane", &fPlane, "fPlane/I");

  //fill branch with purity measurements
  fTreePurity =tfs->make<TTree>("PurityHit","Corrected charge for purity analysis");
  fTreePurity->Branch("fRun", &fRun,"fRun/I");
  fTreePurity->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreePurity->Branch("fEvent", &fEvent, "fEvent/I");
  fTreePurity->Branch("fPlane", &fPlane, "fPlane/I");
  fTreePurity->Branch("fCharge", &fCharge, "fCharge/D");
  fTreePurity->Branch("fDrift", &fDrift, "fDrift/D");
  fTreePurity->Branch("fBin", &fBin, "fBin/I");
  fTreePurity->Branch("fMipIndex", &fMipIndex, "fMipIndex/I");

  fTreePurityMean =tfs->make<TTree>("PurityMean","Summed values for every drift bin");
  fTreePurityMean->Branch("fRun", &fRun,"fRun/I");
  fTreePurityMean->Branch("fSubRun", &fSubRun,"fSubRun/I");
  fTreePurityMean->Branch("fEvent", &fEvent, "fEvent/I");
  fTreePurityMean->Branch("fPlane", &fPlane, "fPlane/I");
  fTreePurityMean->Branch("fBin", &fBin, "fBin/I");
  fTreePurityMean->Branch("fSummedCharge", &fSummedCharge, "fSummedCharge/D");
  fTreePurityMean->Branch("fEntries", &fEntries, "fEntries/I");

  //TODO<<--Generalize these histograms to different geometries
  hTrkTrajectory_0 = tfs->make<TH2D>("hTrkTrajectory_0", "Selected mips hit position view 0;Channel;Ticks", 320, 0, 319, 1667, 0, 1666);
  hTrkTrajectory_1 = tfs->make<TH2D>("hTrkTrajectory_1", "Selected mips hit position view 1;Channel;Ticks", 960, 0, 959, 1667, 0, 1666);

  for(int plane = 0; plane < (int)geom->Nplanes(0); plane++){
    for (int nb=0; nb <fNumOfBins; nb++){
      std::string histname_singlehits = "histname_singlehits"+std::to_string(nb)+"_view"+std::to_string(plane);
      std::string histname_average = "histname_average"+std::to_string(nb)+"_view"+std::to_string(plane);
      std::string histname_num = "histname_numhits"+std::to_string(nb)+"_view"+std::to_string(plane);

      htbin_singlehits[plane][nb] = tfs->make<TH1D>(histname_singlehits.c_str(), ";Single hits charge distribution (fC)", 100, 0, 200);
      htbin[plane][nb] = tfs->make<TH1D>(histname_average.c_str(), ";Average charge (fC)", 100, 0, 200);
      htbin_num[plane][nb] = tfs->make<TH1D>(histname_num.c_str(), ";Number of bins per track", 100, 0, 200);
    }
  }
}

Purity& pdunedp::Purity::operator= ( Purity const &  )

delete

Purity& pdunedp::Purity::operator= ( Purity &&  )

delete

bool pdunedp::Purity::StitchTracks	(	recob::Track	Track1,
		recob::Track	Track2,
		TVector3 &	Edge1,
		TVector3 &	Edge2
	)

Definition at line 600 of file Purity_module.cc.

                                                                                                        {
  //Attempt to sticth two tracks togheter. return the non stitched verteces of the two tracks.
  //  mf::LogVerbatim("pdunedp::Purity") << "Doing Stitch";
  bool Stitch = false;

  double Dx = sqrt(pow(Track1.End().X()-Track2.Vertex().X(),2) + pow(Track1.End().Y()-Track2.Vertex().Y(),2)+pow(Track1.End().Z()-Track2.Vertex().Z(),2));
  double Angle= (180.0/3.14159)*Track1.EndDirection<TVector3>().Angle(Track2.VertexDirection<TVector3>());
  int Criteria = 1;

  if(Dx > sqrt(pow(Track1.End().X()-Track2.End().X(),2) + pow(Track1.End().Y()-Track2.End().Y(),2) + pow(Track1.End().Z()-Track2.End().Z(),2)))
        {
            Dx= sqrt(pow(Track1.End().X()-Track2.End().X(),2) + pow(Track1.End().Y()-Track2.End().Y(),2) + pow(Track1.End().Z()-Track2.End().Z(),2));
            Angle = 180. - (180.0/3.14159)*Track1.EndDirection<TVector3>().Angle(Track2.EndDirection<TVector3>());
            Criteria= 2;
         }

  if(Dx > sqrt(pow(Track1.Vertex().X()-Track2.End().X(),2) + pow(Track1.Vertex().Y()-Track2.End().Y(),2) + pow(Track1.Vertex().Z()-Track2.End().Z(),2)))
  {
            Dx = sqrt(pow(Track1.Vertex().X()-Track2.End().X(),2) + pow(Track1.Vertex().Y()-Track2.End().Y(),2) + pow(Track1.Vertex().Z()-Track2.End().Z(),2));
            Angle = (180.0/3.14159)*Track1.VertexDirection<TVector3>().Angle(Track2.EndDirection<TVector3>());
            Criteria = 3;
         }

  if(Dx > sqrt(pow(Track1.Vertex().X()-Track2.Vertex().X(),2) + pow(Track1.Vertex().Y()-Track2.Vertex().Y(),2) + pow(Track1.Vertex().Z()-Track2.Vertex().Z(),2)))
        {
            Dx = sqrt(pow(Track1.Vertex().X()-Track2.Vertex().X(),2) + pow(Track1.Vertex().Y()-Track2.Vertex().Y(),2) + pow(Track1.Vertex().Z()-Track2.Vertex().Z(),2));
            Angle = 180. - (180.0/3.14159)*Track1.VertexDirection<TVector3>().Angle(Track2.VertexDirection<TVector3>());
            Criteria = 4;
        }

  //check the stitching criteria
  if( (fStitchAngle > Angle) && (fStitchDistance > Dx)){
    Stitch = true;
    switch (Criteria) {
      case 1:
        Edge1 = Track1.Vertex<TVector3>(); Edge2 = Track2.End<TVector3>();
        break;
      case 2:
        Edge1 = Track1.Vertex<TVector3>(); Edge2 = Track2.Vertex<TVector3>();
        break;
      case 3:
        Edge1 = Track1.End<TVector3>(); Edge2 = Track2.Vertex<TVector3>();
        break;
      case 4:
        Edge1 = Track1.End<TVector3>(); Edge2 = Track2.End<TVector3>();
        break;
      default:
        mf::LogError("pdunedp::Purity") << "Unknown error!";
        break;
    }
  }
  return Stitch;
}

Member Data Documentation

int pdunedp::Purity::BadEvents =0

private

Definition at line 187 of file Purity_module.cc.

int pdunedp::Purity::Events =0

private

Definition at line 187 of file Purity_module.cc.

double pdunedp::Purity::fADCIntegral

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fADCIntegral_corrected

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fADCPeak

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fADCPeak_corrected

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fADCtoCharge

private

Definition at line 184 of file Purity_module.cc.

double pdunedp::Purity::fADCToElectrons

private

Definition at line 216 of file Purity_module.cc.

int pdunedp::Purity::fBin

private

Definition at line 211 of file Purity_module.cc.

calo::CalorimetryAlg pdunedp::Purity::fCalorimetryAlg

private

Definition at line 171 of file Purity_module.cc.

art::InputTag pdunedp::Purity::fCalWireModuleLabel

private

Definition at line 174 of file Purity_module.cc.

double pdunedp::Purity::fCharge

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fChargeIntegral

private

Definition at line 197 of file Purity_module.cc.

double pdunedp::Purity::fChi2Ndof

private

Definition at line 200 of file Purity_module.cc.

art::InputTag pdunedp::Purity::fClusterModuleLabel

private

Definition at line 176 of file Purity_module.cc.

double pdunedp::Purity::fCorrectedCharge

private

Definition at line 208 of file Purity_module.cc.

double pdunedp::Purity::fCorrection

private

Definition at line 209 of file Purity_module.cc.

double pdunedp::Purity::fCosX

private

Definition at line 205 of file Purity_module.cc.

double pdunedp::Purity::fCosY

private

Definition at line 205 of file Purity_module.cc.

double pdunedp::Purity::fCosZ

private

Definition at line 205 of file Purity_module.cc.

double pdunedp::Purity::fdEdx

private

Definition at line 210 of file Purity_module.cc.

double pdunedp::Purity::fdQds

private

Definition at line 209 of file Purity_module.cc.

double pdunedp::Purity::fDrift

private

Definition at line 209 of file Purity_module.cc.

double pdunedp::Purity::fECut

private

Definition at line 180 of file Purity_module.cc.

double pdunedp::Purity::fElectronCharge = 1.60217662e-4

private

Definition at line 217 of file Purity_module.cc.

int pdunedp::Purity::fEndTick

private

Definition at line 195 of file Purity_module.cc.

int pdunedp::Purity::fEntries

private

Definition at line 213 of file Purity_module.cc.

int pdunedp::Purity::fEvent

private

Definition at line 188 of file Purity_module.cc.

double pdunedp::Purity::fGodnessOfFit

private

Definition at line 196 of file Purity_module.cc.

art::InputTag pdunedp::Purity::fHitModuleLabel

private

Definition at line 175 of file Purity_module.cc.

double pdunedp::Purity::fHitsCharge

private

Definition at line 201 of file Purity_module.cc.

double pdunedp::Purity::fHitsTrackCharge

private

Definition at line 203 of file Purity_module.cc.

int pdunedp::Purity::fIndex

private

Definition at line 192 of file Purity_module.cc.

double pdunedp::Purity::fLength

private

Definition at line 181 of file Purity_module.cc.

double pdunedp::Purity::fMinDx

private

Definition at line 185 of file Purity_module.cc.

double pdunedp::Purity::fMipCharge

private

Definition at line 204 of file Purity_module.cc.

double pdunedp::Purity::fMipChargeCm = 10

private

Definition at line 218 of file Purity_module.cc.

int pdunedp::Purity::fMipIndex

private

Definition at line 191 of file Purity_module.cc.

double pdunedp::Purity::fMipLength

private

Definition at line 204 of file Purity_module.cc.

int pdunedp::Purity::fMipNum

private

Definition at line 191 of file Purity_module.cc.

double pdunedp::Purity::fMipPhi

private

Definition at line 204 of file Purity_module.cc.

double pdunedp::Purity::fMipTheta

private

Definition at line 204 of file Purity_module.cc.

int pdunedp::Purity::fNHitsMip

private

Definition at line 204 of file Purity_module.cc.

int pdunedp::Purity::fNHitsTrack

private

Definition at line 189 of file Purity_module.cc.

int pdunedp::Purity::fNtotHits

private

Definition at line 189 of file Purity_module.cc.

int pdunedp::Purity::fNtotTracks

private

Definition at line 189 of file Purity_module.cc.

int pdunedp::Purity::fNumOfBins

private

Definition at line 183 of file Purity_module.cc.

int pdunedp::Purity::fNumOfMips

private

Definition at line 190 of file Purity_module.cc.

double pdunedp::Purity::fPeakAmplitude

private

Definition at line 198 of file Purity_module.cc.

int pdunedp::Purity::fPlane

private

Definition at line 192 of file Purity_module.cc.

double pdunedp::Purity::fRange

private

Definition at line 210 of file Purity_module.cc.

int pdunedp::Purity::fRun

private

Definition at line 188 of file Purity_module.cc.

double pdunedp::Purity::fSpacePointX

private

Definition at line 206 of file Purity_module.cc.

double pdunedp::Purity::fSpacePointY

private

Definition at line 206 of file Purity_module.cc.

double pdunedp::Purity::fSpacePointZ

private

Definition at line 206 of file Purity_module.cc.

int pdunedp::Purity::fStartTick

private

Definition at line 194 of file Purity_module.cc.

double pdunedp::Purity::fStitchAngle

private

Definition at line 182 of file Purity_module.cc.

double pdunedp::Purity::fStitchDistance

private

Definition at line 182 of file Purity_module.cc.

int pdunedp::Purity::fSubRun

private

Definition at line 188 of file Purity_module.cc.

double pdunedp::Purity::fSummedADC

private

Definition at line 207 of file Purity_module.cc.

double pdunedp::Purity::fSummedADC_corrected

private

Definition at line 207 of file Purity_module.cc.

double pdunedp::Purity::fSummedCharge

private

Definition at line 213 of file Purity_module.cc.

double pdunedp::Purity::fTotalGain

private

Definition at line 179 of file Purity_module.cc.

double pdunedp::Purity::fTrackCharge

private

Definition at line 202 of file Purity_module.cc.

double pdunedp::Purity::fTrackLength

private

Definition at line 200 of file Purity_module.cc.

art::InputTag pdunedp::Purity::fTrackModuleLabel

private

Definition at line 177 of file Purity_module.cc.

TTree* pdunedp::Purity::fTree

private

Definition at line 225 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreeCalib

private

Definition at line 225 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreeHits

private

Definition at line 226 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreeHitsMip

private

Definition at line 225 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreeMip

private

Definition at line 225 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreePurity

private

Definition at line 226 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreePurityMean

private

Definition at line 226 of file Purity_module.cc.

TTree* pdunedp::Purity::fTreeTrk

private

Definition at line 225 of file Purity_module.cc.

std::map< size_t, std::vector< pdunedp::bHitInfo >[3] > pdunedp::Purity::fTrk2InfoMap

private

Definition at line 220 of file Purity_module.cc.

int pdunedp::Purity::fTrkNum

private

Definition at line 190 of file Purity_module.cc.

std::vector<double> pdunedp::Purity::fVolCut

private

Definition at line 183 of file Purity_module.cc.

int pdunedp::Purity::fWire

private

Definition at line 209 of file Purity_module.cc.

art::ServiceHandle<geo::Geometry> pdunedp::Purity::geom

private

Definition at line 172 of file Purity_module.cc.

std::map<int, int> pdunedp::Purity::goodevents

private

Definition at line 223 of file Purity_module.cc.

TH1D* pdunedp::Purity::htbin[3][100]

private

Definition at line 228 of file Purity_module.cc.

TH1D* pdunedp::Purity::htbin_num[3][100]

private

Definition at line 228 of file Purity_module.cc.

TH1D* pdunedp::Purity::htbin_singlehits[3][100]

private

Definition at line 228 of file Purity_module.cc.

TH2D* pdunedp::Purity::hTrkTrajectory_0

private

Definition at line 229 of file Purity_module.cc.

TH2D* pdunedp::Purity::hTrkTrajectory_1

private

Definition at line 229 of file Purity_module.cc.

int pdunedp::Purity::SkipEvents =0

private

Definition at line 187 of file Purity_module.cc.

std::map<size_t, recob::Track > pdunedp::Purity::TrackList

private

Definition at line 221 of file Purity_module.cc.

---

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

• duneprototypes/duneprototypes/Protodune/dualphase/Purity_module.cc