68 fMassPeak = tfs->
make<TH1F>(
"massPeak",
";Reconstructed mass (MeV)", 125, 0, 250);
77 if(prongs->size() != 2)
return;
80 for(
int photIdx = 0; photIdx < 2; ++photIdx){
81 const rb::Prong& phot = (*prongs)[photIdx];
82 for(
unsigned int cellIdx = 0; cellIdx < phot.
NCell(); ++cellIdx){
91 const double dot = (*prongs)[0].Dir().Dot((*prongs)[1].
Dir());
93 const double m =
sqrt(2*E[0]*E[1]*(1-dot));
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
TutAnalyzer(const fhicl::ParameterSet &pset)
void reconfigure(const fhicl::ParameterSet &pset)
DEFINE_ART_MODULE(TestTMapFile)
Calibrated quantities relying on position in the orthogonal view. To generate a rb::CellHit from a rb...
void analyze(const art::Event &evt)
T get(std::string const &key) const
rb::RecoHit RecoHit(const art::Ptr< rb::CellHit > &chit) const
Return calibrated hit based on assumed W coordinate.
bool IsCalibrated() const
You MUST check here before accessing PECorr, MIP or GeV.
EDAnalyzer(Table< Config > const &config)
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
double dot(const std::vector< double > &x, const std::vector< double > &y)
T * make(ARGS...args) const
A Cluster with defined start position and direction.
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const