zcl::FastMMStudy Class Reference

Inheritance diagram for zcl::FastMMStudy:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	FastMMStudy (fhicl::ParameterSet const &p)
	FastMMStudy (FastMMStudy const &)=delete
	FastMMStudy (FastMMStudy &&)=delete
FastMMStudy &	operator= (FastMMStudy const &)=delete
FastMMStudy &	operator= (FastMMStudy &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
std::string	workerType () const
bool	modifiesEvent () const
void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
art::Handle< art::TriggerResults >	getTriggerResults (Event const &e) const
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)
base_engine_t &	createEngine (seed_t seed)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
detail::CachedProducts &	cachedProducts ()
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Private Member Functions
void	tset (std::string const &branch_name, double const &value)
void	tset2 (std::string const &branch_name, double const &value)
void	StdevCellsPerPlane (rb::Cluster const &hits, double &stdev_cells_per_plane_xview, double &stdev_cells_per_plane_yview) const
void	NumberOfCellsPerLength (rb::Cluster const &hits, double tracklength_of_xview, double tracklength_of_yview, double &number_of_cells_per_length_xview, double &number_of_cells_per_length_yview) const
float	Distance (unsigned const &xcell1, unsigned const &ycell1, unsigned const &plane1, unsigned const &xcell2, unsigned const &ycell2, unsigned const &plane2) const
bool	SurfAssign (rb::CellHit const &hit) const
int	GetSurfaceId (rb::CellHit const &hit) const
double	WeightedCenterCut (rb::Cluster const &hits, float PX_min, float PX_max, float PY_min, float PY_max, float CX_min, float CX_max, float CY_min, float CY_max, double &weighted_off_center_xx, double &weighted_off_center_xz, double &weighted_off_center_yy, double &weighted_off_center_yz) const
double	NumberOfHitsInOverlapPlanesCut (rb::Cluster const &hits, float PX_min, float PX_max, float PY_min, float PY_max) const
int	NumberOfSurfaceHits (rb::Cluster const &hits) const
void	LinFit (const std::vector< double > &x, const std::vector< double > &y, double *fitpar)

Private Attributes
std::string	_slicelabel
std::string	_sliceinstance
int	nxhits_MC
int	nyhits_MC
double	ti_MC
double	tf_MC
double	beta_MC
int	_xMin
int	_xMax
int	_yMin
int	_yMax
int	_zMin
int	_zMax
int	_deltx
int	_delty
int	_deltz
int	_sigmaP
int	_trueHits
int	_nSatHits
double	_stdevcellsperplane
double	_MaxHitsPerLength
double	_MaxOffCenter
int	_MinNumberSurfaceHits
int	PX_max
int	PY_max
int	PX_min
int	PY_min
int	CX_max
int	CX_min
int	CY_max
int	CY_min
long long	TX_max
long long	TY_max
long long	TX_min
long long	TY_min
double	branch_default_value_
TTree *	tree_
TTree *	tree2_
std::map< std::string, double >	t_
std::map< std::string, double >	t2_
art::ServiceHandle< art::TFileService >	tfs_

Detailed Description

Definition at line 61 of file FastMMStudy_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

zcl::FastMMStudy::FastMMStudy ( fhicl::ParameterSet const &  p )

explicit

Definition at line 116 of file FastMMStudy_module.cc.

  :
    EDAnalyzer(p),
    _slicelabel(p.get<std::string>("slice_label")),
    _sliceinstance(p.get<std::string>("slice_instance")),
    _xMin(p.get<int>("x_min")),
    _xMax(p.get<int>("x_max")),
    _yMin(p.get<int>("y_min")),
    _yMax(p.get<int>("y_max")),
    _zMin(p.get<int>("z_min")),
    _zMax(p.get<int>("z_max")),
    _deltx(p.get<int>("x_delt")),
    _delty(p.get<int>("y_delt")),
    _deltz(p.get<int>("z_delt")),
    _sigmaP(p.get<int>("sigma_P")),
    _stdevcellsperplane(p.get<double>("stdev_cells_per_plane")),
    _MaxHitsPerLength(p.get<double>("max_hits_per_length")),
    _MaxOffCenter(p.get<double>("max_off_center")),
    _MinNumberSurfaceHits(p.get<int>("min_number_surface_hits"))
{
}

zcl::FastMMStudy::FastMMStudy ( FastMMStudy const &  )

delete

zcl::FastMMStudy::FastMMStudy ( FastMMStudy &&  )

delete

Member Function Documentation

void zcl::FastMMStudy::analyze ( art::Event const &  e )

overridevirtual

I want to remove isolated hits from beginning

Implements art::EDAnalyzer.

Definition at line 139 of file FastMMStudy_module.cc.

References _deltx, _delty, _deltz, _nSatHits, _sigmaP, _slicelabel, _trueHits, _xMax, _xMin, _yMax, _yMin, _zMax, _zMin, rawdata::RawDigit::ADC(), rb::Cluster::AllCells(), sim::ParticleNavigator::begin(), beta_MC, compare_h5_caf::branch, branch_default_value_, rb::CellHit::Cell(), CX_max, CX_min, CY_max, CY_min, Distance(), sim::ParticleNavigator::end(), art::Event::event(), geom(), art::DataViewImpl::getByLabel(), make_syst_table_plots::h, hits(), MECModelEnuComparisons::i, geo::kX, LinFit(), simb::MCParticle::Mass(), std::max(), std::min(), rb::Cluster::NCell(), NumberOfCellsPerLength(), NumberOfHitsInOverlapPlanesCut(), NumberOfSurfaceHits(), nxhits_MC, nyhits_MC, simb::MCParticle::P(), cheat::BackTracker::ParticleNavigator(), cheat::BackTracker::ParticleToFLSHit(), simb::MCParticle::PdgCode(), rb::CellHit::Plane(), PX_max, PX_min, PY_max, PY_min, art::Event::run(), std::sqrt(), StdevCellsPerPlane(), art::Event::subRun(), t_, rawdata::RawDigit::TDC(), tf_MC, ti_MC, simb::MCParticle::TrackId(), tree2_, tree_, tset(), tset2(), TX_max, TX_min, TY_max, TY_min, rb::CellHit::View(), WeightedCenterCut(), and rb::Cluster::XCell().

                                                {
  //  std::unique_ptr<std::vector<novaddt::TriggerDecision>> 
  //    trigger_decisions(new std::vector<novaddt::TriggerDecision>);
  art::ServiceHandle<geo::Geometry> geom;
  art::ServiceHandle<calib::Calibrator> calib;
  art::Handle< std::vector<rb::Cluster> > slices;
  e.getByLabel(_slicelabel,  slices);
  //e.getByLabel(_slicelabel, _sliceinstance, slices);

  //Get MC info first:
  art::ServiceHandle<cheat::BackTracker> bt;
  const sim::ParticleNavigator& pnav = bt->ParticleNavigator();
  for (sim::ParticleNavigator::const_iterator i = pnav.begin(); i != pnav.end(); ++i) {
    nxhits_MC = 0;
    nyhits_MC = 0;
    const sim::Particle *p = (*i).second;
    if (p->PdgCode()!=42 ) continue;
    //  std::cout<<"Found monopole mctruth"<<std::endl;
    const std::vector<sim::FLSHit>& flshits = bt->ParticleToFLSHit(p->TrackId());
    unsigned ntruehits_Tot = flshits.size();
    if ( ntruehits_Tot==0 ) continue;
    //  std::cout<<"Monopole mctruth hits more than 0!"<<std::endl;
    double momentum2 = (p->P())*(p->P());
    double mass2     = (p->Mass())*(p->Mass());
    beta_MC = sqrt(momentum2 / (momentum2+mass2) );
    ti_MC = flshits[0].GetEntryT();
    tf_MC = flshits[0].GetExitT();
    size_t planeID_i = flshits[0].GetPlaneID();
    if (planeID_i%2) ++nxhits_MC;
    else ++nyhits_MC;
    for (unsigned h =1; h!= flshits.size(); ++h) {
      double ti = flshits[h].GetEntryT();
      double tf = flshits[h].GetExitT(); 
      if (ti < ti_MC) {
        ti_MC = ti;
      }
      if (tf > tf_MC) {
        tf_MC = tf;
      }
      if (flshits[h].GetPlaneID() % 2) ++nxhits_MC;
      else ++nyhits_MC;
    }
    tset2("betaMC",beta_MC);
    tset2("run_number", e.run());
    tset2("subrun_number", e.subRun());
    tset2("event_number", e.event());
    tree2_->Fill();
    //1 monopole per event, so you can break out from the loop:
    break;
  }

  ///I want to remove isolated hits from beginning
  //  std::cout<<"number of slices you are getting: "<<slices->size()<<std::endl;
  // std::cout<<"event start time: "<<event_time*15.625<<std::endl;
  // std::cout<<"TDCT0: "<<TDCT0*15.625<<std::endl;
  // std::cout<<"ti MC and tf MC: "<<ti_MC<<"  "<<tf_MC<<std::endl;
  for (unsigned i = 0; i!= slices->size(); ++i) {
    const rb::Cluster& hits = (*slices)[i];
    PX_max = 0;//P is for plane .enhao
    CX_max = 0;
    CY_max = 0;
    PY_max = 0;
    PX_min = 99999;
    CX_min = 99999;
    CY_min = 99999;
    PY_min = 99999;
    TX_max = hits.XCell(0)->TDC();
    TY_max = TX_max;
    TX_min = hits.XCell(0)->TDC();
    TY_min = TX_min;
    _trueHits = 0;
    _nSatHits   = 0;

    //int nb_hits = hits.NCell();
    //    std::cout<<"slice number = "<<i+1<<" number of hits: "<<nb_hits<<std::endl;
    int  entry_id = 0;
    bool PassThrough = false;
    bool Intrusion   = false;

    int nx = 0;
    int ny = 0;
    double sumADC=0;
    double stdevADC=0;
    double sq_sum_ADC=0;

    std::vector<double> x_hit;
    std::vector<double> y_hit;
    std::vector<double> zy_hit;
    std::vector<double> zx_hit;
    //const double peakToBaselineOffset = 96;
    for (auto const hit_ptr: hits.AllCells()) {
      const rb::CellHit hit = *hit_ptr;
      sumADC+=hit.ADC();
      sq_sum_ADC+=hit.ADC()*hit.ADC();
      if (hit.ADC()>3100) ++ _nSatHits;

      //X view hits
      if (hit.View()== geo::kX) {
        //std::cout<<"X hits"<<std::endl;
        ++nx;
        x_hit.push_back(hit.Cell());
        zx_hit.push_back(hit.Plane());
        if (PX_max < hit.Plane())   PX_max  = hit.Plane();
        if (PX_min > hit.Plane())   PX_min  = hit.Plane();
        if (CX_max < hit.Cell())    CX_max  = hit.Cell();
        if (CX_min > hit.Cell())    CX_min  = hit.Cell();
        if (TX_max < hit.TDC())     TX_max  = hit.TDC();
        if (TX_min > hit.TDC())     TX_min  = hit.TDC();  

        if (!Intrusion) {
          if (hit.Cell() < _xMin+_deltx) {
            entry_id  = 1;
            Intrusion = true;
          }
          else if (hit.Cell()> _xMax-_deltx) {
            entry_id  = 2;
            Intrusion = true;
          }
        }
        else if (!PassThrough) {
          if (hit.Cell() < _xMin+_deltx && entry_id != 1) PassThrough = true;
          else if (hit.Cell()> _xMax-_deltx && entry_id != 2) PassThrough = true;
        }
      } 
      //Y view hits
      else {
        //std::cout<<"Y hits"<<std::endl;
        ++ny;
        y_hit.push_back(hit.Cell());
        zy_hit.push_back(hit.Plane());
        if (PY_max < hit.Plane())   PY_max   = hit.Plane();
        if (PY_min > hit.Plane())   PY_min   = hit.Plane();
        if (CY_max < hit.Cell())    CY_max  = hit.Cell();
        if (CY_min > hit.Cell())    CY_min  = hit.Cell();
        if (TY_max < hit.TDC() )   TY_max   = hit.TDC();
        if (TY_min > hit.TDC() )   TY_min   = hit.TDC();

        if (!Intrusion) {
          if (hit.Cell() < _yMin+_delty) {
            entry_id  = 3;
            Intrusion = true;
          }
          else if (hit.Cell()> _yMax-_delty) {
            entry_id  = 4;
            Intrusion = true;
          }
        }
        else if (!PassThrough) {
          if (hit.Cell() < _yMin+_delty && entry_id != 3) PassThrough = true;
          else if (hit.Cell()> _yMax-_delty && entry_id != 4) PassThrough = true;
        }
      }
      //No matter which view:
      if (!Intrusion) {
        if (hit.Plane() < _zMin+_deltz) {
          entry_id  = 5;
          Intrusion = true;
        }
        else if (hit.Plane() > _zMax-_deltz) {
          entry_id  = 6;
          Intrusion = true;
        }
      }
      else if (!PassThrough) {
        if (hit.Plane() < _zMin+_deltz && entry_id != 5) PassThrough=true;
        else if (hit.Plane() > _zMax-_deltz && entry_id != 6) PassThrough=true;
      }//I think there is a bug if this hit is at corner,  it's cell value will set x intrusion true and entry id =1 and then it will pass through if it's plane val is 0 . by enhao
    }

    //Min-hits on both view:
    //       std::cout<<"xhits number : "<<nx<<"   yhits number : "<<ny<<std::endl;
    double preselection=0;
    if (_trueHits>1 || nx<3 || ny<3||!PassThrough||(std::min(PX_max,PY_max) < std::max(PX_min,PY_min)+_sigmaP)||(std::min(TX_max,TY_max) < std::max(TX_min,TY_min))) { //select cosmic I should add a handler in fcl file.
      continue;
    }
    // if (_trueHits<1 || nx<3 || ny<3||!PassThrough||(std::min(PX_max,PY_max) < std::max(PX_min,PY_min)+_sigmaP)||(std::min(TX_max,TY_max) < std::max(TX_min,TY_min))) {//select MC
    //tset("betaMC",beta_MC);
    //tset("trueHits",_trueHits);
    //tset("preselection",preselection);
    //  if(_trueHits>0) { 
    //    tree_->Fill();
    //std::cout<<"event id: "<<e.event()<<" trueHits: "<<_trueHits<<" nx: "<<nx <<" ny: "<< ny<<" Pass Through: "<<PassThrough<<" Pmax: "<<std::min(PX_max,PY_max)<<" Pmin: "<< std::max(PX_min,PY_min)<<" Tmax: "<<std::min(TX_max,TY_max) <<" Tmin:"<< std::max(TX_min,TY_min)<<std::endl;
    //  }
    //continue;
    //}

    //min number of surface hits test:
    //cut    if (NumberOfSurfaceHits(hits)<_MinNumberSurfaceHits) continue;
    float tracklength =Distance(CX_max, CY_max, PX_max, CX_min, CY_min, PX_min); 
    double tracklength_of_xview=Distance(CX_max, 0, PX_max, CX_min, 0, PX_min);
    double tracklength_of_yview=Distance(0, CY_max, PY_max, 0, CY_min, PY_min);
    double number_of_cells_per_length_xview=0;
    double number_of_cells_per_length_yview=0;
    NumberOfCellsPerLength(hits,tracklength_of_xview,tracklength_of_yview,number_of_cells_per_length_xview,number_of_cells_per_length_yview);
    double fitpar[3];
    LinFit(x_hit, zx_hit, fitpar);
    double r2x = fitpar[2];
    LinFit(y_hit, zy_hit, fitpar);
    double r2y = fitpar[2];
    //std::cout<<"hits/length is "<<hits.size()/tracklength<<std::endl;
    //cut    if (hits.size()/tracklength>_MaxHitsPerLength) continue;
    //I didn't add minhitsperlength because there might be noise hits make tracklength too long.

    //cut    if (StdevCellsPerPlane(hits)==false) continue;
    //cut    if(WeightedCenterCut(hits, PX_min, PX_max, PY_min, PY_max, CX_min,  CX_max,  CY_min,  CY_max )==false) continue;
    //    gROOT->SetBatch(kTRUE);
    stdevADC = std::sqrt( sq_sum_ADC*1.0/(nx+ny) - sumADC*1.0/(nx+ny)*sumADC*1.0/(nx+ny));

    // Fill Tree
    preselection=1;
    for (auto & branch : t_)
      branch.second = branch_default_value_;
    double weighted_off_center_xx=0, weighted_off_center_xz=0, weighted_off_center_yy=0, weighted_off_center_yz=0;
    double stdev_cells_per_plane_xview=0, stdev_cells_per_plane_yview=0;
    tset("run_number", e.run());
    tset("subrun_number", e.subRun());
    tset("event_number", e.event());
    tset("slice_number",i);
    //  tset("event_length", event_length);
    tset("number_of_surface_hits",NumberOfSurfaceHits(hits));
    tset("number_of_hits_per_length",hits.NCell()/tracklength);
    tset("number_of_cells_per_length_xview",number_of_cells_per_length_xview);
    tset("number_of_cells_per_length_yview",number_of_cells_per_length_yview);
    tset("track_length",tracklength);
    StdevCellsPerPlane(hits, stdev_cells_per_plane_xview, stdev_cells_per_plane_yview);
    tset("stdev_cells_per_plane_xview",stdev_cells_per_plane_xview);
    tset("stdev_cells_per_plane_yview",stdev_cells_per_plane_yview);
    tset("weighted_off_center",WeightedCenterCut(hits, PX_min, PX_max, PY_min, PY_max, CX_min,  CX_max,  CY_min,  CY_max, weighted_off_center_xx, weighted_off_center_xz, weighted_off_center_yy, weighted_off_center_yz ));
    tset("weighted_off_center_xx",weighted_off_center_xx);
    tset("weighted_off_center_xz",weighted_off_center_xz);
    tset("weighted_off_center_yy",weighted_off_center_yy);
    tset("weighted_off_center_yz",weighted_off_center_yz);
    tset("meanADC",sumADC*1.0/(nx+ny));
    tset("stdevADC",stdevADC);
    tset("nxhits",nx);
    tset("nyhits",ny);
    tset("percent_of_hits_in_overlap_planes",NumberOfHitsInOverlapPlanesCut(hits,PX_min,PX_max,PY_min,PY_max));
    tset("trueHits",_trueHits);
    tset("betaMC",beta_MC);
    tset("xhitsMC",nxhits_MC);
    tset("yhitsMC",nyhits_MC);
    tset("deltaTDC",std::min(TX_max,TY_max) - std::max(TX_min,TY_min));
    tset("deltaP",std::min(PX_max,PY_max) - std::max(PX_min,PY_min));
    tset("preselection",preselection);
    tset("nSatHits",_nSatHits);
    tset("r2x",r2x);
    tset("r2y",r2y);
    tset("dPX",PX_max-PX_min+1);
    tset("dPY",PY_max-PY_min+1);
    tset("dCX",CX_max-CX_min+1);
    tset("dCY",CY_max-CY_min+1);
    //   std::cout<<"event id: "<<e.event()<<" trueHits: "<<_trueHits<<" nx: "<<nx <<" ny: "<< ny<<" Pass Through: "<<PassThrough<<" Pmax: "<<std::min(PX_max,PY_max)<<" Pmin: "<< std::max(PX_min,PY_min)<<" Tmax: "<<std::min(TX_max,TY_max) <<" Tmin:"<< std::max(TX_min,TY_min)<<std::endl;
    tree_->Fill();
    //if(_trueHits>0) { 
    //  tree_->Fill();
    //}
  }
}

void zcl::FastMMStudy::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 398 of file FastMMStudy_module.cc.

References art::TFileDirectory::make(), t2_, t_, tfs_, tree2_, and tree_.

{
  tree_ = tfs_->make<TTree>("Event", "Monopole Event Tree");
  tree2_ = tfs_->make<TTree>("MC", "Monopole MC Tree");
  std::vector<std::string> branch_names =
  { "run_number", "subrun_number", "event_number", "slice_number","number_of_hits_per_length","number_of_cells_per_length_xview","number_of_cells_per_length_yview","track_length","stdev_cells_per_plane_xview","stdev_cells_per_plane_yview", "weighted_off_center", "weighted_off_center_xx", "weighted_off_center_xz", "weighted_off_center_yy", "weighted_off_center_yz", "percent_of_hits_in_overlap_planes","number_of_surface_hits","betaMC","xhitsMC","yhitsMC","trueHits","preselection","meanADC","stdevADC","nxhits","nyhits","deltaTDC","deltaP","nSatHits","r2x","r2y","dPX","dPY","dCX","dCY"};
  for (auto const& name : branch_names)
    tree_->Branch(name.c_str(), &t_[name], (name + "/D").c_str());

  std::vector<std::string> branch_names2 =
  { "run_number", "subrun_number", "event_number", "betaMC"};
  for (auto const& name : branch_names2)
    tree2_->Branch(name.c_str(), &t2_[name], (name + "/D").c_str());
}

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

  {
    return selectors_;
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 146 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 161 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 171 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed )

inherited

Referenced by evgen::CosmicPionGen::CosmicPionGen(), evgen::CosmicsGen::CosmicsGen(), fuzz::FuzzyKVertex::FuzzyKVertex(), g4n::G4Gen::G4Gen(), photrans::ImprovedTransport::ImprovedTransport(), art::MixFilter< T >::initEngine_(), evgen::IsoGen::IsoGen(), rsim::MakeNoiseSpectrumFile::MakeNoiseSpectrumFile(), photrans::PhotonSystRewriter::PhotonSystRewriter(), photrans::PhotonTransport::PhotonTransport(), rsim::ReadoutSim::ReadoutSim(), evgen::SingleGen::SingleGen(), and evgen::SupernovaGen::SupernovaGen().

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

CurrentProcessingContext const* art::EDAnalyzer::currentContext ( ) const

protectedinherited

Referenced by art::EDAnalyzer::workerType().

float zcl::FastMMStudy::Distance ( unsigned const &  xcell1, unsigned const &  ycell1, unsigned const &  plane1, unsigned const &  xcell2, unsigned const &  ycell2, unsigned const &  plane2  ) const

private

Definition at line 537 of file FastMMStudy_module.cc.

References std::sqrt().

Referenced by analyze().

                                                                                                                                                                                     {
  return std::sqrt((plane1-plane2)*(plane1-plane2)*2.82072+(xcell1-xcell2)*(xcell1-xcell2)+(ycell1-ycell2)*(ycell1-ycell2));
  //This part I should use Martin's smt::Constant::cell width, plane width
}

seed_t art::EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset, char const  key[] = "seed", seed_t const  implicit_seed = -1  )

inherited

Referenced by art::MixFilter< T >::initEngine_().

int zcl::FastMMStudy::GetSurfaceId ( rb::CellHit const &  hit ) const

private

Definition at line 553 of file FastMMStudy_module.cc.

References _deltx, _delty, _deltz, _xMax, _xMin, _yMax, _yMin, _zMax, _zMin, rb::CellHit::Cell(), geo::kX, rb::CellHit::Plane(), and rb::CellHit::View().

                                                           {
  int entry_id=0;
  if (hit.View() == geo::kX) {
    if (hit.Cell() < _xMin+_deltx) {
      entry_id  = 1;
    }
    else if (hit.Cell()> _xMax-_deltx) {
      entry_id  = 2;
    }
  } 
  //Y view hits
  else {
    if (hit.Cell() < _yMin+_delty) {
      entry_id  = 3;
    }
    else if (hit.Cell()> _yMax-_delty) {
      entry_id  = 4;
    }
  }
  //No matter which view:
  if (hit.Plane() < _zMin+_deltz) {
    entry_id  = 5;
  }
  else if (hit.Plane() > _zMax-_deltz) {
    entry_id  = 6;
  }
  return entry_id;

}

art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e ) const

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

  {
    return selectors_.getOneTriggerResults(e);
  }

void zcl::FastMMStudy::LinFit ( const std::vector< double > &  x, const std::vector< double > &  y, double *  fitpar  )

private

Definition at line 632 of file FastMMStudy_module.cc.

References a, b, DEFINE_ART_MODULE(), getGoodRuns4SAM::n, r(), std::sqrt(), submit_syst::x, xx, and submit_syst::y.

Referenced by analyze().

                                                                                                         {
  // http://en.wikipedia.org/wiki/Simple_linear_regression
  const auto n    = x_val.size();
  const auto x  = (std::accumulate(x_val.begin(), x_val.end(), 0.0))/n;                       // <x>
  const auto y  = (std::accumulate(y_val.begin(), y_val.end(), 0.0))/n;                       // <y>
  const auto xx = (std::inner_product(x_val.begin(), x_val.end(), x_val.begin(), 0.0))/n;     // <xx>
  const auto yy = (std::inner_product(y_val.begin(), y_val.end(), y_val.begin(), 0.0))/n;     // <yy>
  const auto xy = (std::inner_product(x_val.begin(), x_val.end(), y_val.begin(), 0.0))/n;     // <xy>

  const auto b    = (xy - x*y)/(xx - x*x);                                                    // slope
  const auto a    = y - b*x;                                                                  // intercept
  const auto r    = (xy - x*y)/sqrt((xx - x*x)*(yy - y*y));                                   // Rxy - coeffcient of determination
  fitpar[0] = a;
  fitpar[1] = b;
  fitpar[2] = r*r;

}

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it )

inherited

Definition at line 189 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 204 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it )

inherited

Definition at line 214 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

  {
    return false;
  }

static cet::exempt_ptr<Consumer> art::Consumer::non_module_context ( )

staticinherited

Referenced by art::RootOutput::beginSubRun(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EmptyEvent::readRun_(), and art::EmptyEvent::readSubRun_().

void zcl::FastMMStudy::NumberOfCellsPerLength ( rb::Cluster const &  hits, double  tracklength_of_xview, double  tracklength_of_yview, double &  number_of_cells_per_length_xview, double &  number_of_cells_per_length_yview  ) const

private

Definition at line 438 of file FastMMStudy_module.cc.

References rb::Cluster::AllCells(), rb::CellHit::Cell(), distance(), febshutoff_auto::end, MECModelEnuComparisons::i, it, geo::kX, geo::kY, rb::CellHit::Plane(), PX_max, PX_min, PY_max, PY_min, and rb::CellHit::View().

Referenced by analyze().

                                                                                                                                                                                                                  {
  //Maybe use hits instead of cells?
  std::vector<std::vector<int>> x_hits_plane(PX_max-PX_min+1,std::vector<int>(0));
  std::vector<std::vector<int>> y_hits_plane(PY_max-PY_min+1,std::vector<int>(0)); 
  for (auto const hit_ptr: hits.AllCells()) {
    const rb::CellHit hit = *hit_ptr;
    if (hit.View() == geo::kX) {
      //std::cout<<"X hits"<<std::endl;
      x_hits_plane[hit.Plane()-PX_min].push_back(hit.Cell()); 
    }
    if (hit.View() == geo::kY) {
      //std::cout<<"Y hits"<<std::endl;
      y_hits_plane[hit.Plane()-PY_min].push_back(hit.Cell()); 
    }
  }
  //following is to get unique number of cells hit per plane
  std::vector<int> x_cells_plane(PX_max-PX_min+1,0);
  std::vector<int> y_cells_plane(PY_max-PY_min+1,0); 
  for (unsigned i=0;i!=x_hits_plane.size();i++){
    std::sort(x_hits_plane[i].begin(),x_hits_plane[i].end());
    std::vector<int>::iterator it;
    it = std::unique (x_hits_plane[i].begin(),x_hits_plane[i].end());
    x_hits_plane[i].resize(std::distance(x_hits_plane[i].begin(),it));
    x_cells_plane[i]=x_hits_plane[i].size();
  }
  for (unsigned i=0;i!=y_hits_plane.size();i++){
    std::sort(y_hits_plane[i].begin(),y_hits_plane[i].end());
    std::vector<int>::iterator it;
    it = std::unique (y_hits_plane[i].begin(),y_hits_plane[i].end());
    y_hits_plane[i].resize(std::distance(y_hits_plane[i].begin(),it));
    y_cells_plane[i]=y_hits_plane[i].size();
  }
  //following is to get number of cells per track length in x view and y view
  int sum_of_cells=0;
  for (unsigned i=0;i!=x_cells_plane.size();i++){
    sum_of_cells+=x_cells_plane[i];
  }
  number_of_cells_per_length_xview=sum_of_cells/tracklength_of_xview;
  sum_of_cells=0;
  for (unsigned i=0;i!=y_cells_plane.size();i++){
    sum_of_cells+=y_cells_plane[i];
  }
  number_of_cells_per_length_yview=sum_of_cells/tracklength_of_yview;

}

double zcl::FastMMStudy::NumberOfHitsInOverlapPlanesCut ( rb::Cluster const &  hits, float  PX_min, float  PX_max, float  PY_min, float  PY_max  ) const

private

Definition at line 619 of file FastMMStudy_module.cc.

References rb::Cluster::AllCells(), cet::sqlite::max(), std::max(), min(), std::min(), rb::Cluster::NCell(), and rb::CellHit::Plane().

Referenced by analyze().

                                                                                                                                     {
  int max=std::min(PX_max,PY_max);
  int min=std::max(PX_min,PY_min);
  double number_hits_in_overlap_planes=0;
  for (auto const hit_ptr: hits.AllCells()) {
    const rb::CellHit hit = *hit_ptr;
    if(hit.Plane()<=max && hit.Plane()>=min)    number_hits_in_overlap_planes++;
  }

  return number_hits_in_overlap_planes/hits.NCell();

}

int zcl::FastMMStudy::NumberOfSurfaceHits ( rb::Cluster const &  hits ) const

private

Definition at line 583 of file FastMMStudy_module.cc.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), and SurfAssign().

Referenced by analyze().

                                                                  {
  int number_of_surface_hits=0;
  for(unsigned i=0;i!=hits.NCell();++i){
    if(SurfAssign(*(hits.Cell(i)))==true) number_of_surface_hits++;
  }
  return number_of_surface_hits;
}

FastMMStudy& zcl::FastMMStudy::operator= ( FastMMStudy const &  )

delete

FastMMStudy& zcl::FastMMStudy::operator= ( FastMMStudy &&  )

delete

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

  {
    return process_name_;
  }

void art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

References string.

  {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

  {
    return selector_config_id_;
  }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Referenced by art::RootOutput::endJob().

void zcl::FastMMStudy::StdevCellsPerPlane ( rb::Cluster const &  hits, double &  stdev_cells_per_plane_xview, double &  stdev_cells_per_plane_yview  ) const

private

Definition at line 483 of file FastMMStudy_module.cc.

References rb::Cluster::AllCells(), rb::CellHit::Cell(), distance(), febshutoff_auto::end, MECModelEnuComparisons::i, it, geo::kX, geo::kY, extractScale::mean, rb::CellHit::Plane(), PX_max, PX_min, PY_max, PY_min, std::sqrt(), HandyFuncs::stdev(), sum, and rb::CellHit::View().

Referenced by analyze().

                                                                                                                                             {
  //fluctuation of number of cells of each plane check by Enhao
  //I need to use cells not hits, because there might be many hits in same cell, and the cosmic shower has many different cells hits
  //Maybe use hits instead of cells?
  std::vector<std::vector<int>> x_hits_plane(PX_max-PX_min+1,std::vector<int>(0));
  std::vector<std::vector<int>> y_hits_plane(PY_max-PY_min+1,std::vector<int>(0)); 
  for (auto const hit_ptr: hits.AllCells()) {
    const rb::CellHit hit = *hit_ptr;
    if (hit.View() == geo::kX) {
      //std::cout<<"X hits"<<std::endl;
      x_hits_plane[hit.Plane()-PX_min].push_back(hit.Cell()); 
    }
    if (hit.View() == geo::kY) {
      //std::cout<<"Y hits"<<std::endl;
      y_hits_plane[hit.Plane()-PY_min].push_back(hit.Cell()); 
    }
  }
  double sum =0;
  double mean = 0;
  double sq_sum = 0;
  double stdev = 0;
  //following is to get unique number of cells hit per plane
  std::vector<int> x_cells_plane(PX_max-PX_min+1,0);
  std::vector<int> y_cells_plane(PY_max-PY_min+1,0); 
  for (unsigned i=0;i!=x_hits_plane.size();i++){
    std::sort(x_hits_plane[i].begin(),x_hits_plane[i].end());
    std::vector<int>::iterator it;
    it = std::unique (x_hits_plane[i].begin(),x_hits_plane[i].end());
    x_hits_plane[i].resize(std::distance(x_hits_plane[i].begin(),it));
    x_cells_plane[i]=x_hits_plane[i].size();
  }
  for (unsigned i=0;i!=y_hits_plane.size();i++){
    std::sort(y_hits_plane[i].begin(),y_hits_plane[i].end());
    std::vector<int>::iterator it;
    it = std::unique (y_hits_plane[i].begin(),y_hits_plane[i].end());
    y_hits_plane[i].resize(std::distance(y_hits_plane[i].begin(),it));
    y_cells_plane[i]=y_hits_plane[i].size();
  }
  sum = std::accumulate(x_cells_plane.begin(),x_cells_plane.end(),0.0);
  mean = sum / x_cells_plane.size();
  sq_sum = std::inner_product(x_cells_plane.begin(),x_cells_plane.end(),x_cells_plane.begin(),0.0);
  stdev = std::sqrt(sq_sum / x_cells_plane.size()-mean * mean);
  //        std::cout<<"x stdev is "<<stdev<<" "<<stdev/mean<<std::endl;
  //if (stdev/mean > _stdevcellsperplane) return false;
  stdev_cells_per_plane_xview=stdev/mean;
  sum = std::accumulate(y_cells_plane.begin(),y_cells_plane.end(),0.0);
  mean = sum / y_cells_plane.size();
  sq_sum = std::inner_product(y_cells_plane.begin(),y_cells_plane.end(),y_cells_plane.begin(),0.0);
  stdev = std::sqrt(sq_sum / y_cells_plane.size()-mean * mean);
  stdev_cells_per_plane_yview=stdev/mean;
  //        std::cout<<"y stdev is "<<stdev<<" "<<stdev/mean<<std::endl;
  //        std::cout<<"sum of x and y stdev is "<<sum_stdev<<std::endl;
  //if (stdev/mean > _stdevcellsperplane) return false;
}

bool zcl::FastMMStudy::SurfAssign ( rb::CellHit const &  hit ) const

private

Definition at line 542 of file FastMMStudy_module.cc.

References _deltx, _deltz, _xMax, _xMin, _zMax, _zMin, rb::CellHit::Cell(), and rb::CellHit::Plane().

Referenced by NumberOfSurfaceHits().

                                                         {
  //Check if this is close to surface
  if      ( hit.Cell()  > _xMax - _deltx   ) return true;
  else if ( hit.Cell()  < _xMin+_deltx     ) return true;
  else if ( hit.Plane() < _zMin+_deltz     ) return true;
  else if ( hit.Plane() > _zMax - _deltz   ) return true;

  return false;
}

void zcl::FastMMStudy::tset ( std::string const &  branch_name, double const &  value  )

private

Definition at line 413 of file FastMMStudy_module.cc.

References ana::assert(), om::cerr, allTimeWatchdog::endl, it, and t_.

Referenced by analyze().

{
  auto it = t_.find(branch_name);
  if (it == t_.end())
  {
    std::cerr << "\n\t Branch " << branch_name << " does not exist! \n"
      << std::endl;
    assert(false);
  }

  it->second = value;
}

void zcl::FastMMStudy::tset2 ( std::string const &  branch_name, double const &  value  )

private

Definition at line 426 of file FastMMStudy_module.cc.

References ana::assert(), om::cerr, allTimeWatchdog::endl, it, and t2_.

Referenced by analyze().

{
  auto it = t2_.find(branch_name);
  if (it == t2_.end())
  {
    std::cerr << "\n\t Branch " << branch_name << " does not exist! \n"
      << std::endl;
    assert(false);
  }

  it->second = value;
}

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::RootOutput::RootOutput().

  {
    return wantAllEvents_;
  }

bool art::EventObserverBase::wantEvent ( Event const &  e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

  {
    return selectors_.wantEvent(e);
  }

double zcl::FastMMStudy::WeightedCenterCut ( rb::Cluster const &  hits, float  PX_min, float  PX_max, float  PY_min, float  PY_max, float  CX_min, float  CX_max, float  CY_min, float  CY_max, double &  weighted_off_center_xx, double &  weighted_off_center_xz, double &  weighted_off_center_yy, double &  weighted_off_center_yz  ) const

private

Definition at line 590 of file FastMMStudy_module.cc.

References std::abs(), rawdata::RawDigit::ADC(), rb::Cluster::AllCells(), rb::CellHit::Cell(), CX_min, CY_min, geo::kX, geo::kY, std::max(), rb::CellHit::Plane(), PX_min, PY_min, and rb::CellHit::View().

Referenced by analyze().

                                                                                                                                                                                                                                                                                                            {
  float Weighted_PX_Center=0.5,Weighted_CX_Center=0.5,Weighted_PY_Center=0.5,Weighted_CY_Center=0.5,SumOfXHitsADC=0,SumOfYHitsADC=0;
  for (auto const hit_ptr: hits.AllCells()) {
    const rb::CellHit hit = *hit_ptr;
    if (hit.View() == geo::kX) {
      SumOfXHitsADC+=hit.ADC();
      Weighted_PX_Center+=hit.ADC()*(hit.Plane()-PX_min);
      Weighted_CX_Center+=hit.ADC()*(hit.Cell()-CX_min);
    }
    if (hit.View() == geo::kY) {
      SumOfYHitsADC+=hit.ADC();
      Weighted_PY_Center+=hit.ADC()*(hit.Plane()-PY_min);
      Weighted_CY_Center+=hit.ADC()*(hit.Cell()-CY_min);
    }
  }
  if(PX_max!=PX_min)  Weighted_PX_Center=Weighted_PX_Center/SumOfXHitsADC/(PX_max-PX_min);
  if(PY_max!=PY_min)  Weighted_PY_Center=Weighted_PY_Center/SumOfYHitsADC/(PY_max-PY_min);
  if(CX_max!=CX_min)  Weighted_CX_Center=Weighted_CX_Center/SumOfXHitsADC/(CX_max-CX_min);
  if(CY_max!=CY_min)  Weighted_CY_Center=Weighted_CY_Center/SumOfYHitsADC/(CY_max-CY_min);
  //  std::cout<<"weighted px, py, cx, cy center is "<<Weighted_PX_Center<<" "<<Weighted_PY_Center<<" "<<Weighted_CX_Center<<" "<<Weighted_CY_Center<<std::endl;
  // if(std::abs(Weighted_PX_Center-0.5)>_MaxOffCenter||std::abs(Weighted_PY_Center-0.5)>_MaxOffCenter||std::abs(Weighted_CX_Center-0.5)>_MaxOffCenter||std::abs(Weighted_CY_Center-0.5)>_MaxOffCenter) return false;
  //else  return true;
  weighted_off_center_xx=std::abs(Weighted_CX_Center-0.5);
  weighted_off_center_xz=std::abs(Weighted_PX_Center-0.5);
  weighted_off_center_yy=std::abs(Weighted_CY_Center-0.5); 
  weighted_off_center_yz=std::abs(Weighted_PY_Center-0.5);
  return std::max(std::max(std::abs(Weighted_PX_Center-0.5),std::abs(Weighted_PY_Center-0.5)),std::max(std::abs(Weighted_CX_Center-0.5),std::abs(Weighted_CY_Center-0.5)));
}

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

    {
      return "WorkerT<EDAnalyzer>";
    }

Member Data Documentation

int zcl::FastMMStudy::_deltx

private

Definition at line 85 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

int zcl::FastMMStudy::_delty

private

Definition at line 86 of file FastMMStudy_module.cc.

Referenced by analyze(), and GetSurfaceId().

int zcl::FastMMStudy::_deltz

private

Definition at line 87 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

double zcl::FastMMStudy::_MaxHitsPerLength

private

Definition at line 92 of file FastMMStudy_module.cc.

double zcl::FastMMStudy::_MaxOffCenter

private

Definition at line 93 of file FastMMStudy_module.cc.

int zcl::FastMMStudy::_MinNumberSurfaceHits

private

Definition at line 94 of file FastMMStudy_module.cc.

int zcl::FastMMStudy::_nSatHits

private

Definition at line 90 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::_sigmaP

private

Definition at line 88 of file FastMMStudy_module.cc.

Referenced by analyze().

std::string zcl::FastMMStudy::_sliceinstance

private

Definition at line 73 of file FastMMStudy_module.cc.

std::string zcl::FastMMStudy::_slicelabel

private

Definition at line 72 of file FastMMStudy_module.cc.

Referenced by analyze().

double zcl::FastMMStudy::_stdevcellsperplane

private

Definition at line 91 of file FastMMStudy_module.cc.

int zcl::FastMMStudy::_trueHits

private

Definition at line 89 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::_xMax

private

Definition at line 80 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

int zcl::FastMMStudy::_xMin

private

Definition at line 79 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

int zcl::FastMMStudy::_yMax

private

Definition at line 82 of file FastMMStudy_module.cc.

Referenced by analyze(), and GetSurfaceId().

int zcl::FastMMStudy::_yMin

private

Definition at line 81 of file FastMMStudy_module.cc.

Referenced by analyze(), and GetSurfaceId().

int zcl::FastMMStudy::_zMax

private

Definition at line 84 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

int zcl::FastMMStudy::_zMin

private

Definition at line 83 of file FastMMStudy_module.cc.

Referenced by analyze(), GetSurfaceId(), and SurfAssign().

double zcl::FastMMStudy::beta_MC

private

Definition at line 77 of file FastMMStudy_module.cc.

Referenced by analyze().

double zcl::FastMMStudy::branch_default_value_

private

Definition at line 97 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::CX_max

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::CX_min

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), and WeightedCenterCut().

int zcl::FastMMStudy::CY_max

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::CY_min

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), and WeightedCenterCut().

int zcl::FastMMStudy::nxhits_MC

private

Definition at line 75 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::nyhits_MC

private

Definition at line 75 of file FastMMStudy_module.cc.

Referenced by analyze().

int zcl::FastMMStudy::PX_max

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), NumberOfCellsPerLength(), and StdevCellsPerPlane().

int zcl::FastMMStudy::PX_min

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), NumberOfCellsPerLength(), StdevCellsPerPlane(), and WeightedCenterCut().

int zcl::FastMMStudy::PY_max

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), NumberOfCellsPerLength(), and StdevCellsPerPlane().

int zcl::FastMMStudy::PY_min

private

Definition at line 95 of file FastMMStudy_module.cc.

Referenced by analyze(), NumberOfCellsPerLength(), StdevCellsPerPlane(), and WeightedCenterCut().

std::map<std::string, double> zcl::FastMMStudy::t2_

private

Definition at line 101 of file FastMMStudy_module.cc.

Referenced by beginJob(), and tset2().

std::map<std::string, double> zcl::FastMMStudy::t_

private

Definition at line 100 of file FastMMStudy_module.cc.

Referenced by analyze(), beginJob(), and tset().

double zcl::FastMMStudy::tf_MC

private

Definition at line 76 of file FastMMStudy_module.cc.

Referenced by analyze().

art::ServiceHandle<art::TFileService> zcl::FastMMStudy::tfs_

private

Definition at line 102 of file FastMMStudy_module.cc.

Referenced by beginJob().

double zcl::FastMMStudy::ti_MC

private

Definition at line 76 of file FastMMStudy_module.cc.

Referenced by analyze().

TTree* zcl::FastMMStudy::tree2_

private

Definition at line 99 of file FastMMStudy_module.cc.

Referenced by analyze(), and beginJob().

TTree* zcl::FastMMStudy::tree_

private

Definition at line 98 of file FastMMStudy_module.cc.

Referenced by analyze(), and beginJob().

long long zcl::FastMMStudy::TX_max

private

Definition at line 96 of file FastMMStudy_module.cc.

Referenced by analyze().

long long zcl::FastMMStudy::TX_min

private

Definition at line 96 of file FastMMStudy_module.cc.

Referenced by analyze().

long long zcl::FastMMStudy::TY_max

private

Definition at line 96 of file FastMMStudy_module.cc.

Referenced by analyze().

long long zcl::FastMMStudy::TY_min

private

Definition at line 96 of file FastMMStudy_module.cc.

Referenced by analyze().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-07/ZukCluster/FastMMStudy_module.cc