Prophecy4fMerger Class Reference

#include <Prophecy4fMerger.h>

Collaboration diagram for Prophecy4fMerger:

Public Member Functions
	Prophecy4fMerger ()
virtual	~Prophecy4fMerger ()
int	alulb4 (double *ps, double *pi, double *pf)
int	alulf4 (double *ps, double *pi, double *pf)
int	alulob (double *ps, double *pi, double *pf)
int	alulof (double *ps, double *pi, double *pf)
double	alupcm (double em0, double em1, double em2)
int	rescms (double *p, double *p1, double *p2, double m1, double m2)
void	setIO (const std::string &powheg, const std::string &prophecy4e, const std::string &prophecy4mu, const std::string &prophecy2e2mu, const std::string &outlhe, bool debug)
void	setRandomSeed (unsigned long long seed)
void	merge ()
bool	isPHevent (const TLorentzVector &higgs, const TLorentzVector &sum_daugh_rest_init)
double	setParticleMass (int id) const
bool	fileExists (const std::string &filename)
void	print (const std::string &field)
void	print (const std::string &field, int value)

Protected Attributes
std::string	m_inPowheg
std::string	m_inProphecy4e
std::string	m_inProphecy4mu
std::string	m_inProphecy2e2mu
std::string	m_outLHE
bool	m_phEvent
bool	m_debug
TRandom3	m_rand

Static Protected Attributes
static constexpr double	m_deltaM = 3.
static constexpr double	m_muonMass = 0.10566
static constexpr double	m_electronMass = 0.00051
static constexpr double	m_tauMass = 1.777
static const long	m_electronID = 11
static const long	m_muonID = 13
static const long	m_tauID = 15
static const long	m_zID = 23
static const long	m_higgsID = 25
static const long	m_photonID = 22
static const long	m_neutrinoEl = 12
static const long	m_neutrinoMu = 14
static const long	m_neutrinoTau = 16

Detailed Description

Definition at line 31 of file Prophecy4fMerger.h.

Constructor & Destructor Documentation

Prophecy4fMerger()

Prophecy4fMerger::Prophecy4fMerger

(

)

Definition at line 20 of file Prophecy4fMerger.cxx.

                                   : m_phEvent(false), m_debug(false), m_rand(0) { 
}

~Prophecy4fMerger()

Prophecy4fMerger::~Prophecy4fMerger ( )

virtual

Definition at line 23 of file Prophecy4fMerger.cxx.

                                   {
}

Member Function Documentation

alulb4()

int Prophecy4fMerger::alulb4	(	double *	ps,
		double *	pi,
		double *	pf
	)

Definition at line 32 of file Prophecy4fMerger.cxx.

                                                              {
 
    // transform pi rest to lab frame of ps - output pf
    // For description, see LOREN4 doc in CERNLIB
 
    if(ps[3] == ps[4]){
        pf[0] = pi[0];
        pf[1] = pi[1];
        pf[2] = pi[2];
        pf[3] = pi[3];
    }
    else{
        const double pf4 = (pi[0] * ps[0] + pi[1] * ps[1] + pi[2] * ps[2] + pi[3] * ps[3]) / ps[4];
        const double fn = (pf4 + pi[3]) / (ps[3] + ps[4]);
        pf[0] = pi[0] + fn * ps[0];
        pf[1] = pi[1] + fn * ps[1];
        pf[2] = pi[2] + fn * ps[2];
        pf[3] = pf4;
    }
    return 0;
  
}

alulf4()

int Prophecy4fMerger::alulf4	(	double *	ps,
		double *	pi,
		double *	pf
	)

Definition at line 55 of file Prophecy4fMerger.cxx.

                                                              {
 
    // transform pi lab to rest frame of ps - output pf
    // For description, see LOREN4 doc in CERNLIB
  
    if(ps[3] == ps[4]){
        pf[0] = pi[0];
        pf[1] = pi[1];
        pf[2] = pi[2];
        pf[3] = pi[3];
    }
    else{
        const double pf4 = (pi[3] * ps[3] - pi[2] * ps[2] - pi[1] * ps[1] - pi[0] * ps[0]) / ps[4];
        const double fn = (pf4 + pi[3]) / (ps[3] + ps[4]);
        pf[0] = pi[0] - fn * ps[0]; 
        pf[1] = pi[1] - fn * ps[1]; 
        pf[2] = pi[2] - fn * ps[2]; 
        pf[3] = pf4;
    }
    return 0;
  
}

alulob()

int Prophecy4fMerger::alulob	(	double *	ps,
		double *	pi,
		double *	pf
	)

Definition at line 78 of file Prophecy4fMerger.cxx.

                                                              {
  
    //extern int alulb4(double *, double *, double *);
  
    alulb4(ps, pi, pf);
    pf[4] = pi[4];
    return 0;
  
} 

alulof()

int Prophecy4fMerger::alulof	(	double *	ps,
		double *	pi,
		double *	pf
	)

Definition at line 88 of file Prophecy4fMerger.cxx.

                                                              {
  
    //extern int alulf4(double *, double *,  double *);
  
    alulf4(ps, pi, pf);
    pf[4] = pi[4];
  
    return 0;
  
} 

alupcm()

double Prophecy4fMerger::alupcm	(	double	em0,
		double	em1,
		double	em2
	)

Definition at line 99 of file Prophecy4fMerger.cxx.

                                                                 {
 
    // calculate the momentum of particles 1 and 2 in the rest frame of 0 given the three masses
 
    double ret_val;
 
    //ems = abs(em1 + em2);
    //emd = abs(em1 - em2);
    double ems = (em1 + em2);
    if(ems < 0){
        ems=-(em1 + em2);
    }
    double emd = (em1 - em2);
    if(emd < 0){
        emd = -(em1 - em2);
    }
  
    if(em0 < ems || em0 < emd){
        ret_val = -1.f;
    }
    else if(em0 == ems || em0 == emd){
        ret_val = 0.f;
    }
    else{
        ret_val = sqrt((em0 + emd) * (em0 - emd) * (em0 + ems) * (em0 - ems)) * .5f / em0;
    }
    return ret_val;
  
}

fileExists()

bool Prophecy4fMerger::fileExists

(

const std::string &

filename

)

Definition at line 521 of file Prophecy4fMerger.cxx.

                                                          {
  
    std::ifstream ifile(filename.c_str());
    if( ifile.good() ){
        ifile.close();
        return true;
    }
    ifile.close();
    return false;
 
}

isPHevent()

bool Prophecy4fMerger::isPHevent	(	const TLorentzVector &	higgs,
		const TLorentzVector &	sum_daugh_rest_init
	)

Definition at line 486 of file Prophecy4fMerger.cxx.

                                                                           {
 
    if( std::abs(higgs.M()-sum_daugh_rest_init.M())>m_deltaM ){
        ANA_MSG_INFO("POWHEG event with Higgs off-mass shell: " << higgs.M() << " GeV");
        return true;
    }
    return false;
  
}

merge()

void Prophecy4fMerger::merge

(

)

Definition at line 192 of file Prophecy4fMerger.cxx.

                            {
 
    print(" Opening Powheg LHE file ... " + m_inPowheg);
    if( !fileExists(m_inPowheg) ){
        ANA_MSG_ERROR("Input Powheg LHE not Found!"
                  << " Aborting ... ");
        return;
    }
    LHEF::Reader readH( m_inPowheg.c_str() );
  
    print(" Opening Prophecy4f LHE file for 4e ... " + m_inProphecy4e);
    if( !fileExists(m_inProphecy4e) ){
        ANA_MSG_INFO( "Input Prophecy4f for 4e LHE not Found! " << m_inProphecy4e
                  << " Aborting ... " );
        return;
    }
    LHEF::Reader read4e( m_inProphecy4e.c_str() );
  
    print(" Opening Prophecy4f LHE file for 4mu ... " + m_inProphecy4mu);
    if( !fileExists(m_inProphecy4mu) ){
        ANA_MSG_ERROR("Input Prophecy4f for 4mu LHE not Found! " << m_inProphecy4mu
                  << " Aborting ... ");
        return;
    }
    LHEF::Reader read4mu( m_inProphecy4mu.c_str() );
  
    print(" Opening Prophecy4f for 2e2mu LHE file ... " + m_inProphecy2e2mu);
    if( !fileExists(m_inProphecy2e2mu) ){
        ANA_MSG_ERROR("Input Prophecy4f for 2e2mu LHE not Found! " << m_inProphecy2e2mu
                  << " Aborting ... ");
        return;
    }
    LHEF::Reader read2e2mu( m_inProphecy2e2mu.c_str() );
  
    print(" Opening Out LHE file ... " + m_outLHE);
    LHEF::Writer writeLHE( m_outLHE.c_str() );
  
    writeLHE.headerBlock() << readH.headerBlock;
    writeLHE.initComments() << readH.initComments;
    writeLHE.heprup = readH.heprup;
    writeLHE.init();
 
    long neve = 0;
    LHEF::Reader* read4f = 0;
  
  
    // Read in each Higgs event and one Prophecy4f event, merge them and write out the merged event
    // randomly chose between 4mu, 4e, 2mu2e and 2e2mu. Note that 2mu2e is always taken as 2e2mu,
    // since it is generated this way.
  
    while( readH.readEvent() ){
        // randomly choose between 4mu, 4e, 2mu2e and 2e2mu
        int ID_l[4] = {-999, -999, -999, -999};
        int decayModeType = 4.0 * m_rand.Rndm();   
     
        if (decayModeType == 0){
            //4e
                        read4f = &read4e;
            for(int i = 0; i < 4; i++) ID_l[i] = m_electronID;
        }
        else if (decayModeType == 1 || decayModeType == 2){
            //2e2mu
                        read4f = &read2e2mu;
            for(int i = 0; i < 2; i++) ID_l[i] = m_electronID;
            for(int i = 2; i < 4; i++) ID_l[i] = m_muonID;
        }
        else {
            //4mu
                        read4f = &read4mu;
            for(int i = 0; i < 4; i++) ID_l[i] = m_muonID;
        }
    
        if( !read4f->readEvent() ){
            ANA_MSG_INFO("Still Powheg events but no more Prophecy4f ");
            ANA_MSG_INFO("events! Writing out LHE file ... Events processed so far " << neve );
            break;
        }
    
        ++neve;
        if(neve%1000==0){
            print("Events processed so far ", neve);
        }
        if( readH.outsideBlock.length() )
            ANA_MSG_INFO( readH.outsideBlock );
    
        int nup_org=readH.hepeup.NUP;
 
        for(int nup=0; nup<nup_org; nup++){
            if(readH.hepeup.IDUP[nup]==m_higgsID)
                readH.hepeup.ISTUP[nup]=2; // set Higgs status code to 2
        }
        writeLHE.hepeup = readH.hepeup;
 
        double P4_Z[5][2] = {{0}}; // momentum and mass of Zi i=1,2 set to 0 P4_Z[5] = (px, py, pz, E, mZ)
        double P4_l[5][5];
    
        double Pph[5] = {-999., -999., -999., -999., -999.};
        double daughter2[5][5];
    
        for(int i=0; i<nup_org; i++){
            if(readH.hepeup.IDUP[i]==m_higgsID) {
                TLorentzVector higgs( readH.hepeup.PUP[i][0],
                                      readH.hepeup.PUP[i][1],
                                      readH.hepeup.PUP[i][2],
                                      readH.hepeup.PUP[i][3] );
    
                TVector3 boostvec = higgs.BoostVector();
                TLorentzRotation boost(-boostvec);
                TLorentzRotation backboost(boostvec);
    
                /* boost to Higgs rest frame */
                TLorentzVector r_higgs(higgs);
                r_higgs*=boost;
    
                /* Here add daughter particles to Higgs in rest
                   boost back with backboost */
                writeLHE.hepeup.resize(nup_org+read4f->hepeup.NUP);
    
                TLorentzVector sum_daughter;
                TLorentzVector sum_daughter_rest;
                TLorentzVector sum_daughter_rest_init;
    
                for(int k=0; k<read4f->hepeup.NUP; k++){
                    TLorentzVector daughter(read4f->hepeup.PUP[k][0],
                                            read4f->hepeup.PUP[k][1],
                                            read4f->hepeup.PUP[k][2],
                                            read4f->hepeup.PUP[k][3]);
                    sum_daughter_rest_init+=daughter;
                }
    
                for(int j=0; j<read4f->hepeup.NUP; j++){
                    TLorentzVector daughter(read4f->hepeup.PUP[j][0],
                                            read4f->hepeup.PUP[j][1],
                                            read4f->hepeup.PUP[j][2],
                                            read4f->hepeup.PUP[j][3]);
                    sum_daughter_rest+=daughter;
      
                    daughter *= backboost;
                    daughter *= higgs.M()/sum_daughter_rest_init.M();
      
                    /* POWHEG event with Higgs off-mass shell */
                    m_phEvent=isPHevent(higgs,sum_daughter_rest_init);
                    if( m_phEvent ){
                        break;
                    }
      
                    sum_daughter+=daughter;
      
                    if( std::abs(read4f->hepeup.IDUP[j])==m_electronID ||
                        std::abs(read4f->hepeup.IDUP[j])==m_muonID      ){
        
                        //ID_l[j]    = read4f->hepeup.IDUP[j];
                        P4_l[0][j] = daughter.Px();
                        P4_l[1][j] = daughter.Py();
                        P4_l[2][j] = daughter.Pz();
                        P4_l[3][j] = daughter.E();
        
                    }       
                    else if( std::abs(read4f->hepeup.IDUP[j])==m_photonID ){
                        Pph[0] = daughter.Px();
                        Pph[1] = daughter.Py();
                        Pph[2] = daughter.Pz();
                        Pph[3] = daughter.E();
                        Pph[4] =0.0;
                    }
                }           
    
                for(int n=0; n<4; n++){
                    if( std::abs(ID_l[n])==m_electronID ||
                        std::abs(ID_l[n])==m_muonID  ) {
        
                        // Z1 -> l1 + l2 and Z2 -> l3 + l4
                        P4_Z[n][0] = P4_l[n][0] + P4_l[n][1];
                        P4_Z[n][1] = P4_l[n][2] + P4_l[n][3];
                    }
                }
    
                double mass1, mass2;
                P4_Z[4][0] = sqrt(P4_Z[3][0] * P4_Z[3][0] - P4_Z[0][0] * P4_Z[0][0] - P4_Z[1][0] * P4_Z[1][0] - P4_Z[2][0] * P4_Z[2][0]);
                P4_Z[4][1] = sqrt(P4_Z[3][1] * P4_Z[3][1] - P4_Z[0][1] * P4_Z[0][1] - P4_Z[1][1] * P4_Z[1][1] - P4_Z[2][1] * P4_Z[2][1]);
                mass1 = setParticleMass(ID_l[0]);
                mass2 = setParticleMass(ID_l[2]);
    
                double p1[5], p2[5], pZ1[5], p3[5], p4[5], pZ2[5];
                for(int u=0; u<5; u++){
      
                    pZ1[u] = P4_Z[u][0];
                    p1[u]  = P4_l[u][0];
                    p2[u]  = P4_l[u][1];
                    pZ2[u] = P4_Z[u][1];
                    p3[u]  = P4_l[u][2];
                    p4[u]  = P4_l[u][3];
      
                }
    
                p1[4] = 0.0;
                p2[4] = 0.0;
                p3[4] = 0.0;
                p4[4] = 0.0;
    
                // Give masses to the leptons using Alpgen/FHerwig routines
                // Note: prophecy leptons are generated massless
                rescms(pZ1, p1, p2, mass1, mass1);
                rescms(pZ2, p3, p4, mass2, mass2);
 
                for(int d =0; d<5; d++){
                    daughter2[0][d] = p1[d];
                    daughter2[1][d] = p2[d];
                    daughter2[2][d] = p3[d];
                    daughter2[3][d] = p4[d];
                    daughter2[4][d] = Pph[d];
                }
 
                // print out children
                bool doDebug = false;
                doDebug = true;
                if (doDebug) {
 
                    ANA_MSG_DEBUG("Event: " << neve);
                    ANA_MSG_DEBUG(" Higgs and children, m, pt, eta, phi" );
                    ANA_MSG_DEBUG("Higgs " << higgs.M() << ", " << higgs.Pt() << "," << higgs.Eta() << "," << higgs.Phi());
 
                    
                    TLorentzVector higgsFromChildren;
                    for(int d =0; d < read4f->hepeup.NUP; d++){
                        TLorentzVector child;
                        child.SetXYZM( daughter2[d][0],
                                       daughter2[d][1],
                                       daughter2[d][2],
                                       daughter2[d][4] );
 
                        ANA_MSG_DEBUG("child " << d << " " << child.M() << ", " << child.Pt() << ", " << child.Eta() << ", " << child.Phi());
                        higgsFromChildren += child;
                    }
                    ANA_MSG_DEBUG("Higgs mass " << higgs.M() << ", from children " << higgsFromChildren.M() << ", diff " << higgs.M() - higgsFromChildren.M());
                }
                break;
            }
        }
 
        /* Write LHE */
        if(m_phEvent){
            writeLHE.hepeup = readH.hepeup;
        }
        else{
      
            for(int jp=0; jp<read4f->hepeup.NUP; jp++){
                if (jp < 4){
                    writeLHE.hepeup.IDUP[nup_org+jp] = (jp % 2 == 0)?  ID_l[jp] : -ID_l[jp];
 
                }
                else{
                    writeLHE.hepeup.IDUP[nup_org+jp] = read4f->hepeup.IDUP[jp];       
                }  
        
                writeLHE.hepeup.ISTUP[nup_org+jp]=read4f->hepeup.ISTUP[jp];
                if( read4f->hepeup.ICOLUP[jp].first==0 &&
                    read4f->hepeup.ICOLUP[jp].second==0  ) {
                    writeLHE.hepeup.ICOLUP[nup_org+jp].first=0;
                    writeLHE.hepeup.ICOLUP[nup_org+jp].second=0;
                }
                else {
                    writeLHE.hepeup.ICOLUP[nup_org+jp].first=read4f->hepeup.ICOLUP[jp].first+100;
                    writeLHE.hepeup.ICOLUP[nup_org+jp].second=read4f->hepeup.ICOLUP[jp].second+100;
                }  
                writeLHE.hepeup.MOTHUP[nup_org+jp].first=3;//i-1;
                writeLHE.hepeup.MOTHUP[nup_org+jp].second=3;//i-1;
                writeLHE.hepeup.PUP[nup_org+jp][0]=daughter2[jp][0];
                writeLHE.hepeup.PUP[nup_org+jp][1]=daughter2[jp][1];
                writeLHE.hepeup.PUP[nup_org+jp][2]=daughter2[jp][2];
                writeLHE.hepeup.PUP[nup_org+jp][3]=daughter2[jp][3];
                writeLHE.hepeup.PUP[nup_org+jp][4]=daughter2[jp][4];
                writeLHE.hepeup.VTIMUP[nup_org+jp]=read4f->hepeup.VTIMUP[jp];
                writeLHE.hepeup.SPINUP[nup_org+jp]=read4f->hepeup.SPINUP[jp];
            }
        }
 
        // set powheg weights to negative, if prophecy weight is -1
 
        // Get Prophecy4f weight
            writeLHE.hepeup.XWGTUP *= read4f->hepeup.XWGTUP;
            double prophecyWeight = read4f->hepeup.XWGTUP;
        if ( writeLHE.hepeup.weights_compressed.size() > 0 ) {
            for ( int i = 1, N = writeLHE.hepeup.weights_compressed.size(); i < N; ++i ) {
                writeLHE.hepeup.weights_compressed[i] *= prophecyWeight;
    
            }
        }
        
        writeLHE.eventComments() << readH.eventComments;
        writeLHE.writeEvent();
    }
}

print() [1/2]

void Prophecy4fMerger::print

(

const std::string &

field

)

Definition at line 533 of file Prophecy4fMerger.cxx.

                                                  {
    if(m_debug){
        ANA_MSG_DEBUG(field);
    }
}

print() [2/2]

void Prophecy4fMerger::print	(	const std::string &	field,
		int	value
	)

Definition at line 539 of file Prophecy4fMerger.cxx.

                                       {
 
    if(m_debug){
        ANA_MSG_DEBUG("DEBUG:: " << value << " " << field);
    }
  
}

rescms()

int Prophecy4fMerger::rescms	(	double *	p,
		double *	p1,
		double *	p2,
		double	m1,
		double	m2
	)

Definition at line 129 of file Prophecy4fMerger.cxx.

                                                                                   {
 
    // p1 and p2 are the four-vector decay products of p. Here we change the masses of p1 and p2 to be
    // m1 and m2. Note that for Prophecy, the masses in p1 and p2 are zero.
  
    double po1[5], po2[5], pcm, pcmo;
      
    const double m = p[4];
    const double mo1 = p1[4];
    const double mo2 = p2[4];
  
    // transform p1 to po1 to be in rest frame of p, sampe for p2
    alulof(p, p1, po1);
    alulof(p, p2, po2);
  
    // Calculate the momentum of p1/p2 in the rest frame of p (po1/po2) for the original mass and new
    // masses, and then rescale the rest frame momenta po1 and po2. (Note: eo1 = eo2 in cms of p, so
    // po1 and po2 in cms of p change according to e^2 - m^2 as m changes.)
    if (std::max(mo1,mo2) < 1e-6){
        pcmo = m / 2.;
    }
    else{
        pcmo = alupcm(m, mo1, mo2);
    }
    pcm = alupcm(m, m1, m2);
  
    //rescale the cms momenta, po1 and po2, to account for the new masses used 
    for (int il = 0; il < 4; il++) {
        po1[il] = pcm / pcmo * po1[il];
        po2[il] = pcm / pcmo * po2[il];
    }
  
    //set energy and mass
    po1[3] = sqrt(pcm*pcm + m1*m1);
    po2[3] = sqrt(pcm*pcm + m2*m2);
    po1[4] = m1;
    po2[4] = m2;
  
    //boost back
    alulob(p, po1, p1);
    alulob(p, po2, p2); 
  
    return 0;
  
}

setIO()

void Prophecy4fMerger::setIO	(	const std::string &	powheg,
		const std::string &	prophecy4e,
		const std::string &	prophecy4mu,
		const std::string &	prophecy2e2mu,
		const std::string &	outlhe,
		bool	debug
	)

Definition at line 175 of file Prophecy4fMerger.cxx.

{
    m_inPowheg        = powheg;
    m_inProphecy4e    = prophecy4e;
    m_inProphecy4mu   = prophecy4mu;
    m_inProphecy2e2mu = prophecy2e2mu;
    m_outLHE          = outlhe;
    m_phEvent         = false;
    m_debug           = debug;
  
}

setParticleMass()

double Prophecy4fMerger::setParticleMass

(

int

id

)

const

Definition at line 498 of file Prophecy4fMerger.cxx.

                                                     {
  
    double mass = 0.0;
 
    if( std::abs(id) == m_neutrinoEl ||
        std::abs(id) == m_neutrinoMu ||
        std::abs(id) == m_neutrinoTau ){
        mass = 0.0;
    }
    else if( std::abs(id) == m_electronID ){
        mass = m_electronMass;
    }
    else if( std::abs(id) == m_muonID ){
        mass = m_muonMass;
    }
    else if( std::abs(id) == m_tauID ){
        mass = m_tauMass;
    }
  
    return mass;
  
}

setRandomSeed()

void Prophecy4fMerger::setRandomSeed

(

unsigned long long

seed

)

Definition at line 26 of file Prophecy4fMerger.cxx.

{
    m_rand.SetSeed(seed);
    ANA_MSG_INFO("Prophecy4fMerger::setRandomSeed: set random seed to " << seed << ".");
}

Member Data Documentation

m_debug

bool Prophecy4fMerger::m_debug

protected

Definition at line 108 of file Prophecy4fMerger.h.

m_deltaM

constexpr double Prophecy4fMerger::m_deltaM = 3.

staticconstexprprotected

Definition at line 110 of file Prophecy4fMerger.h.

m_electronID

const long Prophecy4fMerger::m_electronID = 11

staticprotected

Definition at line 114 of file Prophecy4fMerger.h.

m_electronMass

constexpr double Prophecy4fMerger::m_electronMass = 0.00051

staticconstexprprotected

Definition at line 112 of file Prophecy4fMerger.h.

m_higgsID

const long Prophecy4fMerger::m_higgsID = 25

staticprotected

Definition at line 118 of file Prophecy4fMerger.h.

m_inPowheg

std::string Prophecy4fMerger::m_inPowheg

protected

Definition at line 102 of file Prophecy4fMerger.h.

m_inProphecy2e2mu

std::string Prophecy4fMerger::m_inProphecy2e2mu

protected

Definition at line 105 of file Prophecy4fMerger.h.

m_inProphecy4e

std::string Prophecy4fMerger::m_inProphecy4e

protected

Definition at line 103 of file Prophecy4fMerger.h.

m_inProphecy4mu

std::string Prophecy4fMerger::m_inProphecy4mu

protected

Definition at line 104 of file Prophecy4fMerger.h.

m_muonID

const long Prophecy4fMerger::m_muonID = 13

staticprotected

Definition at line 115 of file Prophecy4fMerger.h.

m_muonMass

constexpr double Prophecy4fMerger::m_muonMass = 0.10566

staticconstexprprotected

Definition at line 111 of file Prophecy4fMerger.h.

m_neutrinoEl

const long Prophecy4fMerger::m_neutrinoEl = 12

staticprotected

Definition at line 120 of file Prophecy4fMerger.h.

m_neutrinoMu

const long Prophecy4fMerger::m_neutrinoMu = 14

staticprotected

Definition at line 121 of file Prophecy4fMerger.h.

m_neutrinoTau

const long Prophecy4fMerger::m_neutrinoTau = 16

staticprotected

Definition at line 122 of file Prophecy4fMerger.h.

m_outLHE

std::string Prophecy4fMerger::m_outLHE

protected

Definition at line 106 of file Prophecy4fMerger.h.

m_phEvent

bool Prophecy4fMerger::m_phEvent

protected

Definition at line 107 of file Prophecy4fMerger.h.

m_photonID

const long Prophecy4fMerger::m_photonID = 22

staticprotected

Definition at line 119 of file Prophecy4fMerger.h.

m_rand

TRandom3 Prophecy4fMerger::m_rand

protected

Definition at line 109 of file Prophecy4fMerger.h.

m_tauID

const long Prophecy4fMerger::m_tauID = 15

staticprotected

Definition at line 116 of file Prophecy4fMerger.h.

m_tauMass

constexpr double Prophecy4fMerger::m_tauMass = 1.777

staticconstexprprotected

Definition at line 113 of file Prophecy4fMerger.h.

m_zID

const long Prophecy4fMerger::m_zID = 23

staticprotected

Definition at line 117 of file Prophecy4fMerger.h.

---

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

• Prophecy4fMerger.h
• Prophecy4fMerger.cxx