Signal2Noise Class Reference

Inheritance diagram for Signal2Noise:

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

Private Member Functions
void	reset ()

Private Attributes
std::string	fRawDigitLabel
std::string	fRawInstanceLabel
std::string	fWireProducerLabel
std::string	fWireInstanceLabel
std::string	fHitModuleLabel
std::string	fTrackModuleLabel
std::string	fCalorimetryModuleLabel
bool	fSaveWaveForm
std::vector< int >	fSelectedWires
TTree *	fEventTree
int	event
int	run
int	subrun
double	evttime
int	year_month_date
int	hour_min_sec
int	ntrks
int	trkid [kMaxTracks]
float	trkstart [kMaxTracks][3]
float	trkend [kMaxTracks][3]
float	trklen [kMaxTracks]
float	trkthetaxz [kMaxTracks]
float	trkthetayz [kMaxTracks]
float	trkstartcosxyz [kMaxTracks][3]
float	trkendcosxyz [kMaxTracks][3]
int	ntrkhits [kMaxTracks][3]
float	trkdqdx [kMaxTracks][3][kMaxHits]
float	trkx [kMaxTracks][3][kMaxHits]
float	trkt [kMaxTracks][3][kMaxHits]
double	trkhitx [kMaxHits][3][kMaxHits]
double	trkhity [kMaxHits][3][kMaxHits]
double	trkhitz [kMaxHits][3][kMaxHits]
int	wireid [kMaxTracks][kMaxHits]
int	chid [kMaxTracks][kMaxHits]
int	tpcid [kMaxTracks][kMaxHits]
float	hit_plane [kMaxTracks][kMaxHits]
double	cosgma [kMaxTracks][kMaxHits]
float	amp [kMaxTracks][kMaxHits]
int	tamp [kMaxTracks][kMaxHits]
float	ped [kMaxTracks][kMaxHits]
float	noiserms [kMaxTracks][kMaxHits]
float	noisermsfit [kMaxTracks][kMaxHits]
int	fNticks
int	fNticksReadout
float	fSampleRate
TH1F *	fWaveForm [10]
float	fMaxNoise = 40.
TH1F *	fWaveFormHist [10]

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 85 of file Signal2Noise_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 176 of file Signal2Noise_module.cc.

  : EDAnalyzer{p}  // ,
  // More initializers here.
{
  // Call appropriate consumes<>() for any products to be retrieved by this module.
  fRawDigitLabel          = p.get<std::string>("RawDigitLabel");
  fRawInstanceLabel       = p.get<std::string>("RawInstanceLabel", "daq");
  fWireProducerLabel      = p.get<std::string>("WireProducerLabel");
  fWireInstanceLabel      = p.get<std::string>("WireInstanceLabel", "dataprep");
  fHitModuleLabel         = p.get<std::string>("HitModuleLabel");
  fTrackModuleLabel       = p.get<std::string>("TrackModuleLabel");
  fCalorimetryModuleLabel = p.get<std::string>("CalorimetryModuleLabel");
  fSaveWaveForm        = p.get<bool>("SaveWaveForm");
  fSelectedWires          = p.get<std::vector<int>>("SelectedWires");

  if (fRawDigitLabel.empty() && fWireProducerLabel.empty()) {
    throw cet::exception("AdcThresholdRoiFinder") << "Both RawDigitLabel and WireProducerLabel are empty";
  }

  if ((!fRawDigitLabel.empty()) && (!fWireProducerLabel.empty())){
    throw cet::exception("AdcThresholdRoiFinder") << "Only one of RawDigitLabel and WireProducerLabel should be set";
  }

  // DetectorPropertiesService
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataForJob();
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataForJob(clockData);
  fNticks = detProp.NumberTimeSamples(); // number of clock ticks per event
  fNticksReadout = detProp.ReadOutWindowSize(); // number of clock ticks per readout window
  fSampleRate = sampling_rate(clockData); // period of the TPC readout electronics clock
  cout << "Numer of clock ticks per event: " << fNticks << endl;
  cout << "Numer of clock ticks per readout window: " << fNticksReadout << endl;
  cout << "Sampling rate: " << fSampleRate << endl;
}

Signal2Noise::Signal2Noise ( Signal2Noise const &  )

delete

Signal2Noise::Signal2Noise ( Signal2Noise &&  )

delete

Member Function Documentation

void Signal2Noise::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 211 of file Signal2Noise_module.cc.

{
  // Implementation of required member function here.
  reset();

  // event info
  run = e.run();
  subrun = e.subRun();
  event = e.id().event();
  art::Timestamp ts = e.time();
  TTimeStamp tts(ts.timeHigh(), ts.timeLow());
  evttime = tts.AsDouble();
  
  //cout << "tts " << tts << endl;
  //cout << "evttime " << evttime << endl;
  
  UInt_t year=0;
  UInt_t month=0;
  UInt_t day=0;

  year_month_date=tts.GetDate(kTRUE,0,&year,&month,&day);

  UInt_t hour=0;
  UInt_t min=0;
  UInt_t sec=0;

  hour_min_sec=tts.GetTime(kTRUE,0,&hour,&min,&sec);

  // channel status
  lariov::ChannelStatusProvider const& channelStatus
    = art::ServiceHandle<lariov::ChannelStatusService const>()->GetProvider();

  // get RawDigit
  std::vector< art::Ptr<raw::RawDigit> > rawdigitlist;
  art::InputTag itag1(fRawDigitLabel, fRawInstanceLabel);
  auto rawdigitListHandle = e.getHandle< std::vector<raw::RawDigit> >(itag1);
  if (rawdigitListHandle) {
    art::fill_ptr_vector(rawdigitlist, rawdigitListHandle);
  }

  //cout << "rawdigitlist.size():  " << rawdigitlist.size() << endl;

  // get Wire
  std::vector<art::Ptr<recob::Wire> > wirelist;
  art::InputTag itag2(fWireProducerLabel, "dataprep");
  auto wireListHandle = e.getHandle< std::vector<recob::Wire> >(itag2);
  if (wireListHandle) {
    art::fill_ptr_vector(wirelist, wireListHandle);
  }

  //cout << "wirelist.size():  " << wirelist.size() << endl;

  // hit
  std::vector< art::Ptr<recob::Hit> > hitlist;
  auto hitListHandle = e.getHandle< std::vector<recob::Hit> >(fHitModuleLabel);
  if (hitListHandle) {
    art::fill_ptr_vector(hitlist, hitListHandle);
  }

  // track
  std::vector< art::Ptr<recob::Track> > tracklist;
  auto trackListHandle = e.getHandle< std::vector<recob::Track> >(fTrackModuleLabel);
  if (trackListHandle) {
    art::fill_ptr_vector(tracklist, trackListHandle);
  }
 
  art::ServiceHandle<geo::Geometry> geom;

  art::FindManyP<recob::Hit> fmtrkhit(trackListHandle, e, fTrackModuleLabel);
  art::FindManyP<anab::Calorimetry> fmtrkcalo(trackListHandle, e, fCalorimetryModuleLabel);
  //art::FindManyP<raw::RawDigit> fmhitrawdigit(hitListHandle, e, fHitModuleLabel);
  //art::FindManyP<recob::Wire> fmhitwire(hitListHandle, e, fHitModuleLabel);
  //art::FindManyP<raw::RawDigit> fmwirerawdigit(wireListHandle, e, "digitwire");
  art::FindManyP<recob::Hit, recob::TrackHitMeta> fmhittrkmeta(trackListHandle, e, fTrackModuleLabel);

  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataFor(e);

  // waveform: save several waveforms for check
  int nwaveform = 0; // only save 10 waveforms
  int nwaveform_plane_0 = 0; // only save 3 waveforms
  int nwaveform_plane_1 = 0; // only save 3 waveforms
  int nwaveform_plane_2 = 0; // only save 3 waveforms

  ntrks = 0;
  for (const auto& trk : tracklist) {

    trkid[ntrks] = trk->ID();
    trklen[ntrks] = trk->Length();

    // TVector3
    auto start = trk->Vertex();
    auto end = trk->End();
    auto start_dir = trk->VertexDirection();
    auto end_dir = trk->EndDirection();

    trkthetaxz[ntrks] = std::atan2(start_dir.X(), start_dir.Z());
    trkthetayz[ntrks]= std::atan2(start_dir.Y(), start_dir.Z());
    
    trkstart[ntrks][0] = start.X();
    trkstart[ntrks][1] = start.Y();
    trkstart[ntrks][2] = start.Z();

    trkend[ntrks][0] = end.X();
    trkend[ntrks][1] = end.Y();
    trkend[ntrks][2] = end.Z();

    trkstartcosxyz[ntrks][0] = start_dir.X();
    trkstartcosxyz[ntrks][1] = start_dir.Y();
    trkstartcosxyz[ntrks][2] = start_dir.Z();

    trkendcosxyz[ntrks][0] = end_dir.X();
    trkendcosxyz[ntrks][1] = end_dir.Y();
    trkendcosxyz[ntrks][2] = end_dir.Z();
    
    // calometry
    if (fmtrkcalo.isValid()) {
      std::vector<art::Ptr<anab::Calorimetry>> calos = fmtrkcalo.at(ntrks);
      for (size_t icalo=0; icalo<calos.size(); icalo++) {
        if (!calos[icalo]) continue;
        if (!calos[icalo]->PlaneID().isValid) continue;
        int planenum = calos[icalo]->PlaneID().Plane;
        if (planenum<0 || planenum>2) continue;

        const size_t NHits = calos[icalo] -> dEdx().size();
        ntrkhits[ntrks][planenum] = NHits;

        double minx = 1e9;
        for (size_t iHit=0; iHit<NHits; ++iHit) {
          cout << "plane: "<< planenum << "; pitch: " << (calos[icalo]->TrkPitchVec())[iHit] << endl;
          if ((calos[icalo]->TrkPitchVec())[iHit]>1) continue;
          const auto& TrkPos = (calos[icalo] -> XYZ())[iHit];
          if (TrkPos.X()<minx)
            minx = TrkPos.X();
        }// loop NHits

        for(size_t iHit = 0; iHit < NHits; ++iHit) {
          if ((calos[icalo]->TrkPitchVec())[iHit]>1) continue;
          const auto& TrkPos1 = (calos[icalo] -> XYZ())[iHit];
          double x = TrkPos1.X()-minx; //subtract the minx to get correct t0
          double XDriftVelocity = detProp.DriftVelocity()*1e-3; //cm/ns
          double t = x/(XDriftVelocity*1000); //change the velocity units to cm/ns to cm/us
          trkx[ntrks][planenum][iHit] = x;
          trkt[ntrks][planenum][iHit] = t;
          trkdqdx[ntrks][planenum][iHit] = (calos[icalo]->dQdx())[iHit];
        } // loop over NHits iHit
      } // end loop over icalo
    } //  end if fmtrkcalo

    // hits associated with trk
    std::vector<art::Ptr<recob::Hit>> allhits = fmtrkhit.at(ntrks); // use either trk.key() or ntrks 
    
    for (size_t ihit=0; ihit<allhits.size(); ihit++) {
      // wire plane information
      unsigned int wireplane = allhits[ihit]->WireID().Plane;
      if (wireplane <0 || wireplane>2) continue;
      unsigned int wire = allhits[ihit]->WireID().Wire;
      unsigned int tpc = allhits[ihit]->WireID().TPC;
      unsigned int channel = allhits[ihit]->Channel();
       
      if (channelStatus.IsBad(channel)) continue;

      // hit position: not all hits are associated with space points, using neighboring space points to interpolate
      double xyz[3] = {-9999.0, -9999.0, -9999.0};
      bool fBadhit = false;

      if (fmhittrkmeta.isValid()) {
        auto vhit = fmhittrkmeta.at(ntrks);
        auto vmeta = fmhittrkmeta.data(ntrks);

        for (size_t ii=0; ii<vhit.size(); ii++) {
          if (vhit[ii].key() == allhits[ihit].key()) {
            
            // nb.  LArPandoraTrackCreation_module.cc fills the max of a signed int in an unsigned int
            // to indicate an invalid index

            if (vmeta[ii]->Index() >= (unsigned int) std::numeric_limits<int32_t>::max()) {
              fBadhit = true;
              continue;
            }

            if (vmeta[ii]->Index() >= trk->NumberTrajectoryPoints()) {
              throw cet::exception("Calorimetry_module.cc") << "Requested track trajectory index "<<vmeta[ii]->Index()<<" exceeds the total number of trajectory points "<<trk->NumberTrajectoryPoints()<<" for track index "<<ntrks<<". Something is wrong with the track reconstruction. Please contact tjyang@fnal.gov!!";
            }

            if (!trk->HasValidPoint(vmeta[ii]->Index())) {
              fBadhit = true;
              continue;
            }
            
            auto loc = trk->LocationAtPoint(vmeta[ii]->Index());
            xyz[0] = loc.X();
            xyz[1] = loc.Y();
            xyz[2] = loc.Z();
            
            break;
          } // if (vhit[ii].key() == allhits[ihit].key())
        
        } // end of for ii

      } // if fmhittrkmeta.isValid()

      if (fBadhit) continue;
      
      trkhitx[ntrks][wireplane][ihit] = xyz[0];
      trkhity[ntrks][wireplane][ihit] = xyz[1];
      trkhitz[ntrks][wireplane][ihit] = xyz[2];

      wireid[ntrks][ihit] = wire;
      chid[ntrks][ihit] = channel;
      tpcid[ntrks][ihit] = tpc;
      hit_plane[ntrks][ihit] = wireplane;

      // calculate track angle w.r.t. wire
      double angleToVert = geom->WireAngleToVertical(geom->View(allhits[ihit]->WireID()), allhits[ihit]->WireID().TPC, allhits[ihit]->WireID().Cryostat)-0.5*::util::pi<>();
      
      //cout << "tpc: " << tpc << "; plane: " << wireplane << ";  wire: " << wire <<  "channel: " << channel << "; WireangleToVert: " << angleToVert << "; x: " << xyz[0] << endl;

      const auto& dir = trk->DirectionAtPoint(0);
      // angleToVert: return the angle w.r.t y+ axis, anti-closewise
      // dir: 3d track direction: u = (x,y,z);
      // vector that perpendicular to wires in yz plane v = (0, sin(angleToVert), cos(angleToVert))   
      // cos gamma = u.Dot(v)/(u.mag()*v.mag()) here, u.mag()=v.mag()=1
      double tmp_cosgamma = abs(sin(angleToVert)*dir.Y() + std::cos(angleToVert)*dir.Z());
      cosgma[ntrks][ihit] = tmp_cosgamma;

      //cout << "track direction: " << trkstartcosxyz[ntrks][0] << ", " << trkstartcosxyz[ntrks][1] << ", " << trkstartcosxyz[ntrks][2]  << endl;
      //cout << "angleToVert: " << angleToVert << endl;
      //cout << "dir: " << dir.X() << ", " << dir.Y() << ", " << dir.Z() << endl;

      /*
      // check wire direction on each plane
      cout << geom->Plane(wireplane).Wire(wire).ThetaZ(true) << endl;
      double wirestart[3];
      double wireend[3];
      geom->Plane(wireplane).Wire(wire).GetStart(wirestart);
      geom->Plane(wireplane).Wire(wire).GetEnd(wireend);
      cout << "wirestart: (" << wirestart[0] << ", " << wirestart[1] << ", "<<  wirestart[2] << ")" << endl;
      cout << "wireend: (" << wireend[0] << ", " << wireend[1] << ", "<<  wireend[2] << ")" << endl;
      */
      
      int datasize = fNticks;
      
      std::vector<float> adcvec(datasize);

      // loop over wires
      // use either rawdigitlist or wirelist (one is empty, the other is not) to find the associated ADCVec with same channel of the hit
      if (!rawdigitlist.empty()) {
        int key_rawdigit = -1;
        for (size_t ird=0; ird<rawdigitlist.size(); ++ird) {
          if (rawdigitlist[ird]->Channel() == channel) {
            key_rawdigit = ird;
            break;
          }
        }

        if (key_rawdigit == -1) continue; // in case of poor bad channel configuration
        int datasize_tmp = rawdigitlist[key_rawdigit]->Samples();
        if (datasize_tmp != datasize) continue; // in case of poor bad channel configuration

        // to use a compressed RawDigit, one has to create a new buffer, fill and use it
        std::vector<short> rawadc(datasize); // create a buffer
        raw::Uncompress(rawdigitlist[key_rawdigit]->ADCs(), rawadc, rawdigitlist[key_rawdigit]->Compression());

        // pedestal
        ped[ntrks][ihit] = rawdigitlist[key_rawdigit]->GetPedestal(); // Pedestal level (ADC counts)

        for (size_t jj=0; jj<rawadc.size(); jj++) {
          adcvec[jj] = rawadc[jj] - ped[ntrks][ihit];
        }
      } // if (!rawdigitlist.empty())
      else if (!wirelist.empty()) {
        int key_wire = -1;
        for (size_t iw=0; iw<wirelist.size(); ++iw) {
          if (wirelist[iw]->Channel() == channel) {
            key_wire = iw;
            break;
          }
        }

        if (key_wire == -1) continue;
        const auto & signal = wirelist[key_wire]->Signal();
        if (int(signal.size()) != datasize) continue;

        for (size_t jj=0; jj<signal.size(); jj++) {
          adcvec[jj] = signal[jj];
        }
      } // if (!wirelist.empty())

      // ROI from the reconstructed hits
      int t0 = allhits[ihit]->PeakTime() - 5*(allhits[ihit]->RMS());
      if (t0<0) t0 = 0;
      int t1 = allhits[ihit]->PeakTime() + 5*(allhits[ihit]->RMS());
      if (t1>= datasize) t1 = datasize - 1;
      //cout << "t0: " << t0 << " ; t1: " << t1 << endl;
      
      // maximum pulse height of waveform
      float temp_max_pulseheight = -9999.;
      //float temp_min_pulseheight = 9999.;
      int temp_t_max_pulseheight = -1; // time in unit of ticks
      for (int itime=t0; itime <=t1; itime++) {
        if (adcvec[itime] > temp_max_pulseheight) {
          temp_max_pulseheight = adcvec[itime];
          temp_t_max_pulseheight = itime;
        }
        
        //if (adcvec[itime] < temp_min_pulseheight) {
        //  temp_min_pulseheight = adcvec[itime];
        //}
      }
      //if (temp_max_pulseheight < 0) {
      //  std::cout << "amp: " << temp_max_pulseheight << "; plane: " << wireplane << std::endl;
      //  cout << "amp min: " << temp_min_pulseheight << endl;
      //  cout << "t0: " << t0 << " ; t1: " << t1 << endl;
      //}

      amp[ntrks][ihit] = temp_max_pulseheight;
      tamp[ntrks][ihit] = temp_t_max_pulseheight;

      // noise rms calculation: ideally, this should be done for all wires, not only wires that have hits
      // method 1: calculate rms directly
      int start_ped = 0; 
      int end_ped = datasize-1; 
      float temp_sum = 0.;
      int temp_number = 0;
      for (int iped=start_ped; iped<=end_ped; iped++) {
        if (iped > t0 && iped < t1) continue; // ideally we should use this to skip ROI region
        if (abs(adcvec[iped]) > fMaxNoise) continue; // skip ROI with a threshold, protection for multiple hits on a wire
        temp_sum += adcvec[iped]*adcvec[iped];
        temp_number++;
      }
      noiserms[ntrks][ihit] = sqrt(temp_sum/temp_number);
     
      // method 2: fit noise histogram with a gaus
      TH1F *h1_noise = new TH1F(TString::Format("noise_trk%d_hit%d",ntrks, (int)ihit), TString::Format("noise_trk%d_hit%d",ntrks, (int)ihit), (int)fMaxNoise, -fMaxNoise, fMaxNoise);

      // fill the readout datasize for each wire with hit: signal is included but would not affect the noise rms since signals are far way from the noise peak. One may also exclude signals by using ROI threshold cuts
      for (int jj=0; jj<datasize; jj++) {
        if (jj > t0 && jj < t1) continue; // ideally we should use this to skip ROI region
        if (abs(adcvec[jj]) > fMaxNoise) continue; // skip ROI with a threshold, protection for multiple hits on a wire
        h1_noise->Fill(adcvec[jj]);
      }
      TF1 *f1_noise = new TF1("f1_noise", "gaus" , -fMaxNoise, fMaxNoise);
      double par[3];
      h1_noise->Fit(f1_noise, "WWQ");
      f1_noise->GetParameters(&par[0]);
      noisermsfit[ntrks][ihit] = par[2]; // sigma from gaus fit

      if (fSaveWaveForm && nwaveform<10 && nwaveform_plane_0<4 && nwaveform_plane_1<4 && nwaveform_plane_2<5) {
        if (wireplane==0) nwaveform_plane_0++;
        if (wireplane==1) nwaveform_plane_1++;
        if (wireplane==2) nwaveform_plane_2++;

        fWaveForm[nwaveform]->SetNameTitle(Form("plane_%d_AdcChannel_%d", wireplane,  channel), Form("AdcChannel%d", channel));

        for (int jj=0; jj<datasize; jj++) {
          fWaveForm[nwaveform]->SetBinContent(jj+1, adcvec[jj]);
        }//fWaveForm

        fWaveFormHist[nwaveform]->SetNameTitle(Form("Noise_%d_AdcChannel_%d", wireplane,  channel), Form("NhistChannel%d", channel));
        for (int tt=1; tt<=h1_noise->GetNbinsX(); tt++){
          fWaveFormHist[nwaveform]->SetBinContent(tt, h1_noise->GetBinContent(tt));
        }//fWaveFormHist
        nwaveform++;
      }
      
      delete h1_noise;
      delete f1_noise;

    } // end of for ihit
    
    ++ ntrks;
  } // end of for trk
  fEventTree->Fill();
}

void Signal2Noise::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 587 of file Signal2Noise_module.cc.

                           {
  art::ServiceHandle<art::TFileService> tfs;
  fEventTree = tfs->make<TTree>("Event", "Event");
  fEventTree->Branch("event", &event, "event/I");
  fEventTree->Branch("run", &run, "run/I");
  fEventTree->Branch("subrun", &subrun, "subrun/I");
  fEventTree->Branch("evttime", &evttime, "evttime/D");
  fEventTree->Branch("year_month_date", &year_month_date, "year_month_date/I");
  fEventTree->Branch("hour_min_sec", &hour_min_sec, "hour_min_sec/I");
  
  // track
  fEventTree->Branch("ntrks", &ntrks, "ntrks/I");
  fEventTree->Branch("trkid", trkid, "trkid[ntrks]/I");
  fEventTree->Branch("trkstart", trkstart, "trkstart[ntrks][3]/F");
  fEventTree->Branch("trkend", trkend, "trkend[ntrks][3]/F");
  fEventTree->Branch("trklen", trklen, "trklen[ntrks]/F");
  fEventTree->Branch("trkthetaxz", trkthetaxz, "trkthetaxz[ntrks]/F");
  fEventTree->Branch("trkthetayz", trkthetayz, "trkthetayz[ntrks]/F");
  fEventTree->Branch("trkstartcosxyz", trkstartcosxyz, "trkstartcosxyz[ntrks][3]/F");
  fEventTree->Branch("trkendcosxyz", trkendcosxyz, "trkendcosxyz[ntrks][3]/F");
 
  fEventTree->Branch("ntrkhits", ntrkhits, "ntrkhits[ntrks][3]/I");
  fEventTree->Branch("trkdqdx", trkdqdx, "trkdqdx[ntrks][3][1000]/F");
  //fEventTree->Branch("trkdedx", trkdedx, "trkdedx[ntrks][3][1000]/F");
  fEventTree->Branch("trkx", trkx, "trkx[ntrks][3][1000]/F");
  fEventTree->Branch("trkt", trkt, "trkt[ntrks][3][1000]/F");
  


  // hit
  fEventTree->Branch("trkhitx", trkhitx, "trkhitx[ntrks][3][1000]/D");
  fEventTree->Branch("trkhity", trkhity, "trkhity[ntrks][3][1000]/D");
  fEventTree->Branch("trkhitz", trkhitz, "trkhitz[ntrks][3][1000]/D");
 
  fEventTree->Branch("wireid", wireid, "wireid[ntrks][1000]/I");
  fEventTree->Branch("chid", wireid, "chid[ntrks][1000]/I");
  fEventTree->Branch("tpcid", wireid, "tpcid[ntrks][1000]/I");

  fEventTree->Branch("hit_plane", hit_plane, "hit_plane[ntrks][1000]/F");
  fEventTree->Branch("ped", ped, "ped[ntrks][1000]/F");
  fEventTree->Branch("amp", amp, "amp[ntrks][1000]/F");
  fEventTree->Branch("tamp", tamp, "tamp[ntrks][1000]/I");
  fEventTree->Branch("cosgma", cosgma, "cosgma[ntrks][1000]/D");

  fEventTree->Branch("noiserms", noiserms, "noiserms[ntrks][1000]/F");
  fEventTree->Branch("noisermsfit", noisermsfit, "noisermsfit[ntrks][1000]/F");


  // waveform
  if (fSaveWaveForm) {
    for (int i=0; i<10; i++) {
      fWaveForm[i] = tfs->make<TH1F>(Form("waveform_%d",i), "wire waveform", fNticksReadout, 0, fNticksReadout);
      fWaveForm[i]->SetStats(0);
      fWaveForm[i]->GetXaxis()->SetTitle("Time [ticks]");
      fWaveForm[i]->GetYaxis()->SetTitle("ADC");
      fWaveFormHist[i] = tfs->make<TH1F>(Form("Noise_%d",i), "noise", (int)fMaxNoise, -fMaxNoise, fMaxNoise);
    }
  }
}

Signal2Noise& Signal2Noise::operator= ( Signal2Noise const &  )

delete

Signal2Noise& Signal2Noise::operator= ( Signal2Noise &&  )

delete

void Signal2Noise::reset ( )

private

Definition at line 647 of file Signal2Noise_module.cc.

                        {
  run = -99999;
  subrun = -99999;
  event = -99999;
  evttime = -99999;
  year_month_date = -99999;
  hour_min_sec = -99999;

  
  ntrks = -99999;
  for (size_t i=0; i<kMaxTracks; ++i) {
    trkid[i] = -1;
    trklen[i] = -1.0;
    trkthetaxz[i] = -9999.0;
    trkthetayz[i] = -9999.0;
    for (int j=0; j<3; j++) {
      trkstart[i][j] = -9999.0;
      trkend[i][j] = -9999.0;
      trkstartcosxyz[i][j] = -9999.0;
      trkendcosxyz[i][j] = -9999.0;
      
      ntrkhits[i][j] = -9999;

      for (int k=0; k<kMaxHits; k++) {
        trkhitx[i][j][k] = -9999.0;
        trkhity[i][j][k] = -9999.0;
        trkhitz[i][j][k] = -9999.0;

        trkdqdx[i][j][k] = -9999.0;
        //trkdedx[i][j][k] = -9999.0;
        trkx[i][j][k] = -9999.0;
        trkt[i][j][k] = -9999.0;
      }
    }

    for (int k=0; k<kMaxHits; k++) {
      wireid[i][k] = -99;
      chid[i][k] = -99;
      tpcid[i][k] = -99;

      hit_plane[i][k] = -1;
      ped[i][k] = -9999.0;
      amp[i][k] = -9999.0;
      tamp[i][k] = -1;
      
      noiserms[i][k] = -999.;
      noisermsfit[i][k] = -999.;

      cosgma[i][k] = -99.;
    
    }

  }


}

Member Data Documentation

float Signal2Noise::amp[kMaxTracks][kMaxHits]

private

Definition at line 159 of file Signal2Noise_module.cc.

int Signal2Noise::chid[kMaxTracks][kMaxHits]

private

Definition at line 153 of file Signal2Noise_module.cc.

double Signal2Noise::cosgma[kMaxTracks][kMaxHits]

private

Definition at line 157 of file Signal2Noise_module.cc.

int Signal2Noise::event

private

Definition at line 122 of file Signal2Noise_module.cc.

double Signal2Noise::evttime

private

Definition at line 125 of file Signal2Noise_module.cc.

std::string Signal2Noise::fCalorimetryModuleLabel

private

Definition at line 111 of file Signal2Noise_module.cc.

TTree* Signal2Noise::fEventTree

private

Definition at line 119 of file Signal2Noise_module.cc.

std::string Signal2Noise::fHitModuleLabel

private

Definition at line 109 of file Signal2Noise_module.cc.

float Signal2Noise::fMaxNoise = 40.

private

Definition at line 171 of file Signal2Noise_module.cc.

int Signal2Noise::fNticks

private

Definition at line 167 of file Signal2Noise_module.cc.

int Signal2Noise::fNticksReadout

private

Definition at line 168 of file Signal2Noise_module.cc.

std::string Signal2Noise::fRawDigitLabel

private

Definition at line 105 of file Signal2Noise_module.cc.

std::string Signal2Noise::fRawInstanceLabel

private

Definition at line 106 of file Signal2Noise_module.cc.

float Signal2Noise::fSampleRate

private

Definition at line 169 of file Signal2Noise_module.cc.

bool Signal2Noise::fSaveWaveForm

private

Definition at line 112 of file Signal2Noise_module.cc.

std::vector<int> Signal2Noise::fSelectedWires

private

Definition at line 113 of file Signal2Noise_module.cc.

std::string Signal2Noise::fTrackModuleLabel

private

Definition at line 110 of file Signal2Noise_module.cc.

TH1F* Signal2Noise::fWaveForm[10]

private

Definition at line 170 of file Signal2Noise_module.cc.

TH1F* Signal2Noise::fWaveFormHist[10]

private

Definition at line 172 of file Signal2Noise_module.cc.

std::string Signal2Noise::fWireInstanceLabel

private

Definition at line 108 of file Signal2Noise_module.cc.

std::string Signal2Noise::fWireProducerLabel

private

Definition at line 107 of file Signal2Noise_module.cc.

float Signal2Noise::hit_plane[kMaxTracks][kMaxHits]

private

Definition at line 155 of file Signal2Noise_module.cc.

int Signal2Noise::hour_min_sec

private

Definition at line 127 of file Signal2Noise_module.cc.

float Signal2Noise::noiserms[kMaxTracks][kMaxHits]

private

Definition at line 163 of file Signal2Noise_module.cc.

float Signal2Noise::noisermsfit[kMaxTracks][kMaxHits]

private

Definition at line 164 of file Signal2Noise_module.cc.

int Signal2Noise::ntrkhits[kMaxTracks][3]

private

Definition at line 141 of file Signal2Noise_module.cc.

int Signal2Noise::ntrks

private

Definition at line 130 of file Signal2Noise_module.cc.

float Signal2Noise::ped[kMaxTracks][kMaxHits]

private

Definition at line 161 of file Signal2Noise_module.cc.

int Signal2Noise::run

private

Definition at line 123 of file Signal2Noise_module.cc.

int Signal2Noise::subrun

private

Definition at line 124 of file Signal2Noise_module.cc.

int Signal2Noise::tamp[kMaxTracks][kMaxHits]

private

Definition at line 160 of file Signal2Noise_module.cc.

int Signal2Noise::tpcid[kMaxTracks][kMaxHits]

private

Definition at line 154 of file Signal2Noise_module.cc.

float Signal2Noise::trkdqdx[kMaxTracks][3][kMaxHits]

private

Definition at line 142 of file Signal2Noise_module.cc.

float Signal2Noise::trkend[kMaxTracks][3]

private

Definition at line 133 of file Signal2Noise_module.cc.

float Signal2Noise::trkendcosxyz[kMaxTracks][3]

private

Definition at line 138 of file Signal2Noise_module.cc.

double Signal2Noise::trkhitx[kMaxHits][3][kMaxHits]

private

Definition at line 148 of file Signal2Noise_module.cc.

double Signal2Noise::trkhity[kMaxHits][3][kMaxHits]

private

Definition at line 149 of file Signal2Noise_module.cc.

double Signal2Noise::trkhitz[kMaxHits][3][kMaxHits]

private

Definition at line 150 of file Signal2Noise_module.cc.

int Signal2Noise::trkid[kMaxTracks]

private

Definition at line 131 of file Signal2Noise_module.cc.

float Signal2Noise::trklen[kMaxTracks]

private

Definition at line 134 of file Signal2Noise_module.cc.

float Signal2Noise::trkstart[kMaxTracks][3]

private

Definition at line 132 of file Signal2Noise_module.cc.

float Signal2Noise::trkstartcosxyz[kMaxTracks][3]

private

Definition at line 137 of file Signal2Noise_module.cc.

float Signal2Noise::trkt[kMaxTracks][3][kMaxHits]

private

Definition at line 145 of file Signal2Noise_module.cc.

float Signal2Noise::trkthetaxz[kMaxTracks]

private

Definition at line 135 of file Signal2Noise_module.cc.

float Signal2Noise::trkthetayz[kMaxTracks]

private

Definition at line 136 of file Signal2Noise_module.cc.

float Signal2Noise::trkx[kMaxTracks][3][kMaxHits]

private

Definition at line 144 of file Signal2Noise_module.cc.

int Signal2Noise::wireid[kMaxTracks][kMaxHits]

private

Definition at line 152 of file Signal2Noise_module.cc.

int Signal2Noise::year_month_date

private

Definition at line 126 of file Signal2Noise_module.cc.

---

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

• duneprototypes/duneprototypes/Iceberg/Signal2Noise/Signal2Noise_module.cc