SMMADC_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       SMMADC
// Module Type: analyzer
// File:        SMMADC_module.cc
//
// Generated at Thu Sep 28 15:31:46 2013 by Zukai Wang using artmod
// from art v1_02_06.
////////////////////////////////////////////////////////////////////////

#include "RecoBase/CellHit.h"
#include "Simulation/Particle.h"
#include "Simulation/ParticleNavigator.h"
#include "MCCheater/BackTracker.h"
#include "Simulation/FLSHit.h"
#include "Simulation/FLSHitList.h"

#include "art/Framework/Core/EDAnalyzer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"
#include "art/Framework/Services/Optional/TFileService.h"

#include <fstream>
#include <string>
#include <vector>
#include "TTree.h"

class TTree;

namespace ddt {
  class SMMADC;
}


class ddt::SMMADC : public art::EDAnalyzer {
public:
  explicit SMMADC(fhicl::ParameterSet const & p);
  virtual ~SMMADC();

  void beginJob();
  void analyze(art::Event const & e) override;

private:
  std::string fCellHitInput;
  TTree*      SMtree;

  float beta,plc,edep;
  int ADC;
};


ddt::SMMADC::SMMADC(fhicl::ParameterSet const & p) :
  EDAnalyzer(p),
  fCellHitInput(p.get<std::string>("CellHitInput"))
{
}

ddt::SMMADC::~SMMADC()
{
}

void ddt::SMMADC::beginJob()
{
 
  art::ServiceHandle<art::TFileService> tfs;
  SMtree = tfs->make<TTree>("SMHitsTree","");

  SMtree->Branch("beta",&beta,"beta/F");
  SMtree->Branch("ADC",&ADC,"ADC/I");
  SMtree->Branch("plc",&plc,"plc/F");
  SMtree->Branch("edep",&edep,"edep/F");
}

void ddt::SMMADC::analyze(art::Event const & e)
{
  art::Handle<std::vector<rb::CellHit> > hitcol;
  e.getByLabel(fCellHitInput, hitcol);

  art::ServiceHandle<cheat::BackTracker> bt;

  const sim::ParticleNavigator& pnav = bt->ParticleNavigator();
  for (sim::ParticleNavigator::const_iterator i = pnav.begin(); i != pnav.end(); ++i){
    const sim::Particle *p = (*i).second;
    if (p->PdgCode() == 42) {
      beta = sqrt((p->P())*(p->P())/((p->P())*(p->P())+ (p->Mass())*(p->Mass())) );
      break;
    }
  }


  for (unsigned int i = 0; i < hitcol->size(); ++i) {
    art::Ptr<rb::CellHit> hit(hitcol, i);
    if (!bt->IsNoise(hit)) {
      ADC = hit->ADC();
      const std::vector<sim::FLSHit >& flshits =  bt->HitToFLSHit(hit);
      plc = flshits[0].GetTotalPathLength();
      edep = flshits[0].GetEdep();

      SMtree->Fill();
    }
  }  

  return;
}



DEFINE_ART_MODULE(ddt::SMMADC)