LBNEQuickPiToNuVect Class Reference

#include <LBNEQuickPiToNu.hh>

Public Member Functions
void	YesDoThis ()
bool	doIt () const
void	open (const std::string &prefix)
void	close ()
void	addPion (int trNum1, int sign1, double e, G4ThreeVector p)
void	RingHangerPion (int trNum1, int sign, G4ThreeVector p, G4ThreeVector xPos)
void	blessPion (int trNum1, double e, G4ThreeVector p)
void	reset ()
void	evtOut ()

Static Public Member Functions
static LBNEQuickPiToNuVect *	Instance ()

Private Member Functions
	LBNEQuickPiToNuVect ()
	LBNEQuickPiToNuVect (LBNEQuickPiToNuVect const &)
LBNEQuickPiToNuVect &	operator= (LBNEQuickPiToNuVect const &)

Private Attributes
bool	doItForReal
std::vector< LBNEQuickPiToNu >	data
std::ofstream	fOut

Static Private Attributes
static LBNEQuickPiToNuVect *	m_pInstance =0

Detailed Description

Definition at line 49 of file LBNEQuickPiToNu.hh.

Constructor & Destructor Documentation

LBNEQuickPiToNuVect::LBNEQuickPiToNuVect ( )

private

Definition at line 29 of file LBNEQuickPiToNu.cc.

                                         {
 doItForReal = false;
}

LBNEQuickPiToNuVect::LBNEQuickPiToNuVect ( LBNEQuickPiToNuVect const &  )

private

Member Function Documentation

void LBNEQuickPiToNuVect::addPion ( int  trNum1, int  sign1, double  e, G4ThreeVector  p  )

inline

Definition at line 74 of file LBNEQuickPiToNu.hh.

                                                                           {
       LBNERunManager* pRunManager=(LBNERunManager*)LBNERunManager::GetRunManager();
       int evtno = pRunManager->GetCurrentEvent()->GetEventID();       
       LBNEQuickPiToNu a; a.evtNum = evtno; a.sign = sign1;
       a.trNum = trNum1; a.piEnergy=e;
        for (size_t k=0; k!=3; k++) a.piMomentumTarget[k]= p[k];
        data.push_back(a);
     }

void LBNEQuickPiToNuVect::blessPion ( int  trNum1, double  e, G4ThreeVector  p  )

inline

Definition at line 100 of file LBNEQuickPiToNu.hh.

                                                                  {
      bool foundPion = false;
      for (std::vector<LBNEQuickPiToNu>::iterator it=data.begin(); it!=data.end(); it++) {
        if ((it->trNum == -9999) || (it->trNum == -9998)) continue;
        if ((it->trNum != trNum1) && ((it->trNum-10000) != trNum1)) continue;
        for (size_t k=0; k!= 3; k++) {
           it->nuMomentum[k] = p[k];
         }
         it->nuEnergy = e;
         foundPion = true; return;
        if (!foundPion) { 
            std::ostringstream mStrStr;
            mStrStr << "Found Unknown pion.. Should not happen.. Stop here. Fatal " << std::endl;
            G4String mStr(mStrStr.str());
            G4Exception("LBNERunAction::BeginOfRunAction", " ", RunMustBeAborted, mStr.c_str()); 
        }
      }
    }

void LBNEQuickPiToNuVect::close ( )

inline

Definition at line 73 of file LBNEQuickPiToNu.hh.

{fOut.close();} 

bool LBNEQuickPiToNuVect::doIt ( ) const

inline

Definition at line 62 of file LBNEQuickPiToNu.hh.

{return doItForReal;}

void LBNEQuickPiToNuVect::evtOut ( )

inline

Definition at line 120 of file LBNEQuickPiToNu.hh.

                         {
      for (std::vector<LBNEQuickPiToNu>::const_iterator it=data.begin(); it!=data.end(); it++) {
        if ((it->trNum == -9999) || (it->trNum == -9998)) continue;
        if (std::abs(it->nuEnergy) < 0.001) continue; // gave no neutrino..
        fOut << " " << it->evtNum << " " << it->trNum << " " << it->sign << " " << it->piEnergy;
        fOut << " " << it->nuEnergy;
        fOut << " " << it->piMomentumTarget[0]/it->piMomentumTarget[2] << " " 
                    << it->piMomentumTarget[1]/it->piMomentumTarget[2];
        for (size_t k=0; k != 3; k++) fOut << " " << it->piPosAtHangerRing[k];
        const double rSq = it->piPosAtHangerRing[0]*it->piPosAtHangerRing[0] + 
                           it->piPosAtHangerRing[1]*it->piPosAtHangerRing[1];
        if (rSq > 1.0e-6)                  
           fOut << " " << it->piMomentumAtHangerRing[0]/it->piMomentumAtHangerRing[2] << " "
                       << it->piMomentumAtHangerRing[1]/it->piMomentumAtHangerRing[2];
        else fOut << " 0. 0. ";        
        fOut << " " << it->nuMomentum[0]/it->nuMomentum[2] << " " << it->nuMomentum[1]/it->nuMomentum[2];
        fOut << std::endl;       
      }
      this->reset();
    }

LBNEQuickPiToNuVect * LBNEQuickPiToNuVect::Instance ( )

static

Definition at line 21 of file LBNEQuickPiToNu.cc.

                                                   {
 
   if (m_pInstance== NULL) m_pInstance=new LBNEQuickPiToNuVect;

  return m_pInstance;
}

void LBNEQuickPiToNuVect::open ( const std::string &  prefix )

inline

Definition at line 63 of file LBNEQuickPiToNu.hh.

                                               {
       if (fOut.is_open()) {
          std::cerr << " From LBNEQuickPiToNuVect, file has been already opened, close it " << std::endl;
          fOut.close();
       }
       std::string fName(prefix); fName += std::string("QuickPitoNu_v1.txt");
       fOut.open(fName.c_str());
       fOut << " evtNum trNum sign piEnergy nuEnergy piPxTgt piPyTgt pixH2 piyH2 pizH2 piPxH2 piPyH2 ";
       fOut << " nuPx nuPy " << std::endl; // P stands for slope ,Px/Pz... 
     }

LBNEQuickPiToNuVect& LBNEQuickPiToNuVect::operator= ( LBNEQuickPiToNuVect const &  )

private

void LBNEQuickPiToNuVect::reset ( )

inline

Definition at line 118 of file LBNEQuickPiToNu.hh.

{ data.clear();}

void LBNEQuickPiToNuVect::RingHangerPion ( int  trNum1, int  sign, G4ThreeVector  p, G4ThreeVector  xPos  )

inline

Definition at line 82 of file LBNEQuickPiToNu.hh.

                                                                                            {
       bool foundPion = false;
       for (std::vector<LBNEQuickPiToNu>::iterator it=data.begin(); it!=data.end(); it++) {
         if ((it->trNum == -9999) || (it->trNum == -9998)) continue;
         if (it->trNum != trNum1) continue;
         for (size_t k=0; k!= 3; k++) {
           it->piMomentumAtHangerRing[k] = p[k];
           it->piPosAtHangerRing[k] = xPos[k];
         }
         foundPion = true; return;
       }
       if (!foundPion) { 
       const double e = std::sqrt((139.58*139.58) + p[0]*p[0] +  p[1]*p[1] + p[2]*p[2]);
        std::cerr << " Found new secondary pion, declare it with a higher track number " << std::endl;
        addPion(10000+trNum1, sign, e, p);
       }
     }

void LBNEQuickPiToNuVect::YesDoThis ( )

inline

Definition at line 61 of file LBNEQuickPiToNu.hh.

{ doItForReal = true;}

Member Data Documentation

std::vector<LBNEQuickPiToNu> LBNEQuickPiToNuVect::data

private

Definition at line 56 of file LBNEQuickPiToNu.hh.

bool LBNEQuickPiToNuVect::doItForReal

private

Definition at line 55 of file LBNEQuickPiToNu.hh.

std::ofstream LBNEQuickPiToNuVect::fOut

private

Definition at line 57 of file LBNEQuickPiToNu.hh.

LBNEQuickPiToNuVect * LBNEQuickPiToNuVect::m_pInstance =0

staticprivate

Definition at line 54 of file LBNEQuickPiToNu.hh.

---

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

• g4lbne/include/LBNEQuickPiToNu.hh
• g4lbne/src/LBNEQuickPiToNu.cc