sbnd::NuAna Class Reference

A module to extract the info from the Monte Carlo generator GENIE. More...

Inheritance diagram for sbnd::NuAna:

Public Member Functions
	NuAna (fhicl::ParameterSet const &pset)
virtual	~NuAna ()
void	beginJob ()
void	analyze (const art::Event &evt)
void	beginSubRun (const art::SubRun &subrun)
void	reset ()

Private Attributes
NuAnaAlg	fNuAnaAlg
std::string	fMode
double	fBaseline
bool	fFullOscTrue
std::string	fGenieModuleLabel
std::string	fLarg4ModuleLabel
bool	fFluxReweight
bool	fXSecReweight
std::vector< std::string >	fWeights
std::vector< float >	fWeightRangeSigma
unsigned int	fRandSeed
int	fNWeights
unsigned int	FinalRandSeed
std::vector< std::vector< float > >	reweightingSigmas
TTree *	fTreeTot
std::vector< std::vector< float > >	eventWeights
std::vector< float >	neutMom
std::vector< std::vector< float > >	leptonPos
std::vector< std::vector< float > >	leptonMom
double	Elep
double	thetaLep
double	phiLep
int	NPi0FinalState
int	NGamma
int	NChargedPions
bool	foundAllPhotons
std::vector< std::vector< float > >	p1PhotonConversionPos
std::vector< std::vector< float > >	p1PhotonConversionMom
std::vector< std::vector< float > >	p2PhotonConversionPos
std::vector< std::vector< float > >	p2PhotonConversionMom
std::vector< std::vector< float > >	miscPhotonConversionPos
std::vector< std::vector< float > >	miscPhotonConversionMom
std::vector< std::vector < std::vector< float > > >	chargedPionPos
std::vector< std::vector < std::vector< float > > >	chargedPionMom
std::vector< int >	chargedPionSign
std::vector< std::vector< float > >	pionPos
std::vector< std::vector< float > >	pionMom
std::vector< int >	GeniePDG
std::vector< std::vector< float > >	GenieMomentum
std::vector< std::string >	GenieProc
std::vector< std::vector< float > >	genieReweights
int	iflux
int	nuchan
int	inno
int	parid
int	ndecay
int	isCC
int	mode
double	enugen
double	nuleng
std::vector< float >	vertex
std::vector< float >	neutVertexInWindow
std::vector< float >	ParentVertex
std::vector< float >	nuParentMomAtDecay
std::vector< float >	nuParentMomAtProd
std::vector< float >	nuParentMomTargetExit
int	ptype
int	tptype
double	POT

Detailed Description

A module to extract the info from the Monte Carlo generator GENIE.

Definition at line 50 of file NuAna_module.cc.

Constructor & Destructor Documentation

sbnd::NuAna::NuAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 180 of file NuAna_module.cc.

    : EDAnalyzer(pset)
    , fMode             (pset.get< std::string >              ("Mode"))
    , fBaseline         (pset.get< double >                   ("Baseline"))
    , fFullOscTrue      (pset.get< bool >                     ("FullOsc"))
    , fGenieModuleLabel (pset.get< std::string >              ("GenieModuleLabel"))
    , fLarg4ModuleLabel (pset.get< std::string >              ("LArG4ModuleLabel"))
    , fFluxReweight     (pset.get< bool >                     ("FluxReweight"))
    , fXSecReweight     (pset.get< bool >                     ("XSecReweight"))
    , fWeights          (pset.get< std::vector<std::string> > ("Weights"))
    , fRandSeed         (pset.get< unsigned int >             ("RandSeed"))
    , fNWeights         (pset.get< int >                      ("NWeights"))
  {

    // This function sets up the ttrees
    // get access to the TFile service
    art::ServiceHandle<art::TFileService> tfs;

    // PoTTree  = tfs->make<TTree>("POT", "POT");

    std::ofstream libraryFile;
    libraryFile.open("loadLibs.C");
    libraryFile << "#include <vector>\n";
    libraryFile << "#ifdef __MAKECINT__\n";
    libraryFile << "#pragma link C++ class std::vector<std::vector<float> >+;\n";
    libraryFile << "#pragma link C++ class std::vector<std::vector<std::vector<float> > >+;\n";
    libraryFile << "#endif\n";
    libraryFile.close();
    gROOT -> ProcessLine(".L loadLibs.C+");

    fTreeTot = tfs->make<TTree>("EventsTot", "Event info for ALL types");

    fTreeTot->Branch("POT",&POT,"POT/D");

    // Neutrino/event variables:
    fTreeTot->Branch("iflux",    &iflux,    "iflux/I");
    fTreeTot->Branch("nuchan",   &nuchan,   "nuchan/I");
    fTreeTot->Branch("inno",     &inno,     "inno/I");
    fTreeTot->Branch("enugen",   &enugen,   "enugen/D");
    fTreeTot->Branch("nuleng",   &nuleng,   "nuleng/D");
    fTreeTot->Branch("isCC",     &isCC,     "isCC/I");
    fTreeTot->Branch("mode",     &mode,     "mode/I");
    fTreeTot->Branch("ThetaLep", &thetaLep, "ThetaLep/D");
    fTreeTot->Branch("PhiLep",   &phiLep,   "PhiLep/D");
    fTreeTot->Branch("Elep",     &Elep,     "Elep/D");
    fTreeTot->Branch("neutMom", "neutMom", &neutMom, 32000, 0);
    fTreeTot->Branch("vertex",  "vertex",  &vertex, 32000, 0);

    // Genie Variables
    fTreeTot->Branch("GeniePDG",       &GeniePDG);
    fTreeTot->Branch("GenieMomentum", "std::vector<std::vector<float> >",  &GenieMomentum);
    fTreeTot->Branch("GenieProc",      &GenieProc);
          // "var","std::vector<std::vector<float> >",&vec

    // Flux variables:
    fTreeTot->Branch("ptype", &ptype, "ptype/I");
    fTreeTot->Branch("tptype",&tptype,"tptype/I");
    fTreeTot->Branch("ndecay",&ndecay,"ndecay/I");
    fTreeTot->Branch("neutVertexInWindow", "neutVertexInWindow", &neutVertexInWindow, 32000, 0);
    fTreeTot->Branch("ParentVertex", "ParentVertex", &ParentVertex, 32000, 0);
    fTreeTot->Branch("nuParentMomAtDecay", "nuParentMomAtDecay", &nuParentMomAtDecay, 32000, 0);
    fTreeTot->Branch("nuParentMomAtProd",  "nuParentMomAtProd",  &nuParentMomAtProd, 32000, 0);
    fTreeTot->Branch("nuParentMomTargetExit", "nuParentMomTargetExit", &nuParentMomTargetExit, 32000, 0);

    // larg4 info
    fTreeTot->Branch("leptonPos","std::vector< std::vector<float> > ", &leptonPos, 32000, 0);
    fTreeTot->Branch("leptonMom","std::vector< std::vector<float> > ", &leptonMom, 32000, 0);
    fTreeTot->Branch("NPi0FinalState", &NPi0FinalState, "NPi0FinalState/I");
    fTreeTot->Branch("NGamma",         &NGamma,         "NGamma/I");
    fTreeTot->Branch("FoundPhotons",   &foundAllPhotons,"FoundAllPhotons/B");
    fTreeTot->Branch("p1PhotonConversionPos","std::vector< std::vector<float> > ", &p1PhotonConversionPos, 32000, 0);
    fTreeTot->Branch("p1PhotonConversionMom","std::vector< std::vector<float> > ", &p1PhotonConversionMom, 32000, 0);
    fTreeTot->Branch("p2PhotonConversionPos","std::vector< std::vector<float> > ", &p2PhotonConversionPos, 32000, 0);
    fTreeTot->Branch("p2PhotonConversionMom","std::vector< std::vector<float> > ", &p2PhotonConversionMom, 32000, 0);
    fTreeTot->Branch("miscPhotonConversionPos","std::vector< std::vector<float> > ", &miscPhotonConversionPos, 32000, 0);
    fTreeTot->Branch("miscPhotonConversionMom","std::vector< std::vector<float> > ", &miscPhotonConversionMom, 32000, 0);
    fTreeTot->Branch("PionPos","std::vector< std::vector<float> >", &pionPos, 32000, 0);
    fTreeTot->Branch("PionMom","std::vector< std::vector<float> >", &pionMom, 32000, 0);
    fTreeTot->Branch("ChargedPionPos","std::vector<std::vector< std::vector<float> > >",&chargedPionPos, 32000,0);
    fTreeTot->Branch("ChargedPionMom","std::vector<std::vector< std::vector<float> > >",&chargedPionMom, 32000,0);
    fTreeTot->Branch("ChargedPionSign","ChargedPionSign",&chargedPionSign, 32000,0);


    if (fXSecReweight){
      fTreeTot->Branch("MultiWeight","MultiWeight",&eventWeights,32000,0);
      fTreeTot->Branch("FinalRandSeed", &FinalRandSeed, "FinalRandSeed/I");
    }
    if (fFluxReweight){
      fTreeTot->Branch("MultiWeight","MultiWeight",&eventWeights,32000,0);
    }

  }

virtual sbnd::NuAna::~NuAna ( )

inlinevirtual

Definition at line 53 of file NuAna_module.cc.

{};

Member Function Documentation

void sbnd::NuAna::analyze

(

const art::Event &

evt

)

Definition at line 379 of file NuAna_module.cc.

                                        {

    this -> reset();

    //get the MC generator information out of the event
    //these are all handles to mc information.
    art::Handle< std::vector<simb::MCTruth> > mclistGENIE;
    art::Handle< std::vector<simb::MCParticle> > mclistLARG4;
    art::Handle< std::vector<simb::MCFlux> > mcflux;
    art::Handle< std::vector<simb::GTruth> > mcgtruth;

    // Start by getting the POT for this neutrino's file:
    // art::Handle< sumdata::POTSummary > potHandle;
    // evt.getByLabel(fGenieModuleLabel, potHandle);
    // POT = potHandle->totpot;


    //actually go and get the stuff
    evt.getByLabel(fGenieModuleLabel,mclistGENIE);
    evt.getByLabel(fGenieModuleLabel,mcflux);
    evt.getByLabel(fGenieModuleLabel,mcgtruth);
    if (!fFullOscTrue)
        evt.getByLabel(fLarg4ModuleLabel,mclistLARG4);

    // contains the mctruth object from genie
    art::Ptr<simb::MCTruth> mc(mclistGENIE,0);

    // contains the mcflux object
    art::Ptr<simb::MCFlux > flux(mcflux,0);

    // contains the gtruth object
    art::Ptr<simb::GTruth > gtruth(mcgtruth,0);



    // Contains the neutrino info
    simb::MCNeutrino neutrino = mc -> GetNeutrino();

    // std::cout << "The interaction info is: \n"
    //           << "  gtruth->ftgtPDG................." << gtruth->ftgtPDG << "\n"
    //           << "  gtruth->ftgtZ..................." << gtruth->ftgtZ << "\n"
    //           << "  gtruth->ftgtA..................." << gtruth->ftgtA << "\n"
    //           << "  gtruth->fGint..................." << gtruth->fGint << "\n"
    //           << "  gtruth->fGscatter..............." << gtruth->fGscatter << "\n"
    //           << "  gtruth->fweight................." << gtruth->fweight << "\n"
    //           << "  neutrino.Mode()................." << neutrino.Mode() << "\n"
    //           << "  neutrino.InteractionType()......" << neutrino.InteractionType() << "\n"
    //           << "  neutrino.CCNC()................." << neutrino.CCNC() << "\n"
    //           << "  neutrino.Target()..............." << neutrino.Target() << "\n"
    //           << std::endl;

    // Now start packing up the variables to fill the tree
    // In general, have the algorithm do this:


    // get the basic neutrino info:
    fNuAnaAlg.packNeutrinoInfo( &neutrino,
                                nuchan,
                                inno,
                                enugen,
                                isCC,
                                mode,
                                thetaLep,
                                phiLep,
                                Elep,
                                neutMom,
                                vertex);


    // Pack up the flux info:
    fNuAnaAlg.packFluxInfo(     flux,
                                ptype, tptype, ndecay,
                                neutVertexInWindow,
                                ParentVertex,
                                nuParentMomAtDecay,
                                nuParentMomAtProd,
                                nuParentMomTargetExit);

    // nuleng needs special attention:
    TVector3 length(vertex[0] - ParentVertex[0],
                    vertex[1] - ParentVertex[1],
                    100*fBaseline + vertex[2] - ParentVertex[2]);
    nuleng = length.Mag();

    // Pack up the genie info:
    fNuAnaAlg.packGenieInfo(    mc,
                                GeniePDG,
                                GenieMomentum,
                                GenieProc,
                                NPi0FinalState,
                                NGamma,
                                NChargedPions);


    // pack up the larg4 photon info:
    if(!fFullOscTrue)
        fNuAnaAlg.packLarg4Info(mclistLARG4, isCC, NPi0FinalState,
                                NGamma, NChargedPions,
                                leptonPos,
                                leptonMom,
                                p1PhotonConversionPos,
                                p1PhotonConversionMom,
                                p2PhotonConversionPos,
                                p2PhotonConversionMom,
                                miscPhotonConversionPos,
                                miscPhotonConversionMom,
                                pionPos,
                                pionMom,
                                chargedPionPos,
                                chargedPionMom,
                                chargedPionSign);






    // If needed, set up the weights.

    if (fFluxReweight){
        fNuAnaAlg.packFluxWeight(flux, eventWeights);
    }

    // Find a new weight for this event:
    // std::vector<std::vector<float> > weights;
    if (fXSecReweight){
      fNuAnaAlg.calcWeight(mc, gtruth,eventWeights);
      // std::cout << "Printing weights for this event:\n";
      // for (auto s : fWeights) std::cout << s << "\t";
      //    std::cout << "Total\n";
      //    for (unsigned int row = 0; row < eventWeights.front().size(); row ++){
      //      for (unsigned int column = 0; column < eventWeights.size(); column ++){
      //        std::cout << eventWeights[column][row] << "\t";
      //      }
      //    std::cout << std::endl;
      //    }

    }

    // Fill the ttree:
    fTreeTot->Fill();


    return;
  }

void sbnd::NuAna::beginJob

(

)

Definition at line 274 of file NuAna_module.cc.

                      {

    std::cout << "\nThe mode is " << fMode << std::endl;
    std::cout << "This set of events is ";
    if (fFullOscTrue) std::cout << "fullosc" << std::endl;
    else std::cout << "not fullosc" << std::endl;
    std::cout << "The baseline for this detector is "
              << fBaseline << "m." << std::endl << std::endl;

    reset();
    // gROOT->ProcessLine(".L loadDictionaries.C+");

    art::ServiceHandle<geo::Geometry> geom;
    // configure the geometry in the worker function:
    fNuAnaAlg.configureGeometry(geom);

    if (fXSecReweight){
      std::vector<reweight> reweights;
      fNuAnaAlg.parseWeights(fWeights, reweights);
      FinalRandSeed = fNuAnaAlg.prepareSigmas(fNWeights,
                                              fRandSeed, reweightingSigmas);
      fNuAnaAlg.configureReWeight(reweights, reweightingSigmas);
    }

    return;
  }

void sbnd::NuAna::beginSubRun

(

const art::SubRun &

subrun

)

Definition at line 302 of file NuAna_module.cc.

                                                {

    // Go through POTSummary objects
    art::Handle< sumdata::POTSummary > potHandle;
    subrun.getByLabel(fGenieModuleLabel, potHandle);
    const sumdata::POTSummary& potSum = (*potHandle);
    POT = potSum.totpot;
    // PoTTree->Fill();
    return;
  }

void sbnd::NuAna::reset

(

)

Definition at line 312 of file NuAna_module.cc.

                   {

    // This method takes ANYTHING that goes into the ntuple and sets it to default.

    // Start by making sure the appropriate vectors are cleared and emptied:
    eventWeights.clear();
    leptonPos.clear();
    leptonMom.clear();
    p1PhotonConversionPos.clear();
    p1PhotonConversionMom.clear();
    p2PhotonConversionPos.clear();
    p2PhotonConversionMom.clear();
    miscPhotonConversionPos.clear();
    miscPhotonConversionMom.clear();
    pionPos.clear();
    pionMom.clear();
    chargedPionPos.clear();
    chargedPionMom.clear();
    chargedPionSign.clear();

    GeniePDG.clear();
    GenieMomentum.clear();
    GenieProc.clear();

    NPi0FinalState  = 0;
    NChargedPions   = 0;
    NGamma          = 0;
    foundAllPhotons = true;

    ptype  = 0;
    tptype = 0;

    // POT = 0;

    // Clear out the root vectors
    vertex.clear();
    ParentVertex.clear();
    nuParentMomAtDecay.clear();
    nuParentMomAtProd.clear();
    nuParentMomTargetExit.clear();
    neutMom.clear();

    // iflux    = -999;             // represents the sample, 0 = nu, 1 = nu_fosc, 2 = nubar, 3 = nubar_fosc
    nuchan   = -999;             // Type of neutrino interaction
    inno     = -999;             // Type of neutrino, PDG code.  Not sure why it's inno...
    parid    = -999;             // ID of the parent, for flux reweighing
    ndecay   = -999;             // Type of decay, for flux reweighing
    isCC     = -999;             // isCC event? isCC == 1 means CC, isCC == 0 means NC
    mode     = -999;             // beam mode
    enugen   = -999;             // Energy of the neutrino (both)
    nuleng   = -999;             // Length the neutrino traveled.
    Elep     = -999;
    thetaLep = -999;
    phiLep   = -999;
    // POT      = -999;             //POT of the whole file.

    // eventReweight.clear();
    // eventReweight.resize(7);
    // for (auto & weight : eventReweight)
    // {
    //   weight.clear();
    //   weight.resize(1000);
    // }

    return;
  }

Member Data Documentation

std::vector< std::vector < std::vector<float> > > sbnd::NuAna::chargedPionMom

private

Definition at line 132 of file NuAna_module.cc.

std::vector< std::vector < std::vector<float> > > sbnd::NuAna::chargedPionPos

private

Definition at line 131 of file NuAna_module.cc.

std::vector< int > sbnd::NuAna::chargedPionSign

private

Definition at line 133 of file NuAna_module.cc.

double sbnd::NuAna::Elep

private

Definition at line 107 of file NuAna_module.cc.

double sbnd::NuAna::enugen

private

Definition at line 154 of file NuAna_module.cc.

std::vector<std::vector<float> > sbnd::NuAna::eventWeights

private

Definition at line 98 of file NuAna_module.cc.

double sbnd::NuAna::fBaseline

private

Definition at line 71 of file NuAna_module.cc.

bool sbnd::NuAna::fFluxReweight

private

Definition at line 79 of file NuAna_module.cc.

bool sbnd::NuAna::fFullOscTrue

private

Definition at line 72 of file NuAna_module.cc.

std::string sbnd::NuAna::fGenieModuleLabel

private

Definition at line 75 of file NuAna_module.cc.

unsigned int sbnd::NuAna::FinalRandSeed

private

Definition at line 87 of file NuAna_module.cc.

std::string sbnd::NuAna::fLarg4ModuleLabel

private

Definition at line 76 of file NuAna_module.cc.

std::string sbnd::NuAna::fMode

private

Definition at line 70 of file NuAna_module.cc.

NuAnaAlg sbnd::NuAna::fNuAnaAlg

private

Definition at line 65 of file NuAna_module.cc.

int sbnd::NuAna::fNWeights

private

Definition at line 85 of file NuAna_module.cc.

bool sbnd::NuAna::foundAllPhotons

private

Definition at line 115 of file NuAna_module.cc.

unsigned int sbnd::NuAna::fRandSeed

private

Definition at line 84 of file NuAna_module.cc.

TTree* sbnd::NuAna::fTreeTot

private

Definition at line 95 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::fWeightRangeSigma

private

Definition at line 83 of file NuAna_module.cc.

std::vector<std::string> sbnd::NuAna::fWeights

private

Definition at line 82 of file NuAna_module.cc.

bool sbnd::NuAna::fXSecReweight

private

Definition at line 81 of file NuAna_module.cc.

std::vector<std::vector<float> > sbnd::NuAna::GenieMomentum

private

Definition at line 142 of file NuAna_module.cc.

std::vector<int> sbnd::NuAna::GeniePDG

private

Definition at line 141 of file NuAna_module.cc.

std::vector<std::string> sbnd::NuAna::GenieProc

private

Definition at line 143 of file NuAna_module.cc.

std::vector< std::vector< float > > sbnd::NuAna::genieReweights

private

Definition at line 145 of file NuAna_module.cc.

int sbnd::NuAna::iflux

private

Definition at line 147 of file NuAna_module.cc.

int sbnd::NuAna::inno

private

Definition at line 149 of file NuAna_module.cc.

int sbnd::NuAna::isCC

private

Definition at line 152 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::leptonMom

private

Definition at line 106 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::leptonPos

private

Definition at line 105 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::miscPhotonConversionMom

private

Definition at line 129 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::miscPhotonConversionPos

private

Definition at line 128 of file NuAna_module.cc.

int sbnd::NuAna::mode

private

Definition at line 153 of file NuAna_module.cc.

int sbnd::NuAna::NChargedPions

private

Definition at line 113 of file NuAna_module.cc.

int sbnd::NuAna::ndecay

private

Definition at line 151 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::neutMom

private

Definition at line 100 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::neutVertexInWindow

private

Definition at line 157 of file NuAna_module.cc.

int sbnd::NuAna::NGamma

private

Definition at line 112 of file NuAna_module.cc.

int sbnd::NuAna::NPi0FinalState

private

Definition at line 111 of file NuAna_module.cc.

int sbnd::NuAna::nuchan

private

Definition at line 148 of file NuAna_module.cc.

double sbnd::NuAna::nuleng

private

Definition at line 155 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::nuParentMomAtDecay

private

Definition at line 159 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::nuParentMomAtProd

private

Definition at line 160 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::nuParentMomTargetExit

private

Definition at line 161 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::p1PhotonConversionMom

private

Definition at line 123 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::p1PhotonConversionPos

private

Definition at line 122 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::p2PhotonConversionMom

private

Definition at line 125 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::p2PhotonConversionPos

private

Definition at line 124 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::ParentVertex

private

Definition at line 158 of file NuAna_module.cc.

int sbnd::NuAna::parid

private

Definition at line 150 of file NuAna_module.cc.

double sbnd::NuAna::phiLep

private

Definition at line 108 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::pionMom

private

Definition at line 138 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::pionPos

private

Definition at line 137 of file NuAna_module.cc.

double sbnd::NuAna::POT

private

Definition at line 163 of file NuAna_module.cc.

int sbnd::NuAna::ptype

private

Definition at line 162 of file NuAna_module.cc.

std::vector< std::vector<float> > sbnd::NuAna::reweightingSigmas

private

Definition at line 89 of file NuAna_module.cc.

double sbnd::NuAna::thetaLep

private

Definition at line 108 of file NuAna_module.cc.

int sbnd::NuAna::tptype

private

Definition at line 162 of file NuAna_module.cc.

std::vector<float> sbnd::NuAna::vertex

private

Definition at line 156 of file NuAna_module.cc.

---

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

• NuAna_module.cc