MissingEtTruth Class Reference

#include <MissingEtTruth.h>

Inheritance diagram for MissingEtTruth:

Collaboration diagram for MissingEtTruth:

Public Types
enum	TruthIndex {   Int = 0 , NonInt = 1 , IntCentral = 2 , IntFwd = 3 ,   IntOutCover = 4 , Muons = 5 , Size = 6 , NotValid = 999 }
enum	Source {   Calo = 0 , Calib = 1 , Truth = 2 , Muon = 3 ,   CaloRaw = 4 , Final = 5 , Cryo = 6 , Topo = 7 ,   Ref = 8 , DeadMaterial = 9 , Track = 10 , ObjElectron = 100 ,   ObjMuon = 101 , ObjJet = 102 , ObjIdTrk = 103 , ObjMiniJet = 104 ,   ObjRest = 105 , ObjFinal = 106 , Unknown = 999 }

Public Member Functions
	MissingEtTruth ()
	MissingEtTruth (MissingEtRegions *aRegion)
	MissingEtTruth (MissingET::Source aSource, MissingEtRegions *aRegion)
virtual	~MissingEtTruth ()
void	addPart (HepMC::ConstGenParticlePtr aPart, double etaFull)
void	setExTruth (TruthIndex aTruth, double theEx)
void	setEyTruth (TruthIndex aTruth, double theEy)
void	setEtSumTruth (TruthIndex aTruth, double theEtSum)
void	setExTruthVec (std::vector< double > &&exCaloVec)
void	setEyTruthVec (std::vector< double > &&exCaloVec)
void	setEtSumTruthVec (std::vector< double > &&etSumCaloVec)
double	exTruth (TruthIndex aTruth) const
double	eyTruth (TruthIndex aTruth) const
double	etSumTruth (TruthIndex aTruth) const
const std::vector< double > &	exTruthVec () const
const std::vector< double > &	eyTruthVec () const
const std::vector< double > &	etSumTruthVec () const
virtual void	setEx (double theEx)
virtual void	setEy (double theEy)
virtual void	setEtSum (double theSum)
virtual void	add (double theEx, double theEy, double theEt, double theEta, double weight)
virtual double	etx () const
virtual double	ety () const
virtual double	sumet () const
virtual double	et () const
virtual double	phi () const
virtual const MissingEtRegions *	getRegions () const
virtual Source	getSource () const

Protected Member Functions
void	setup ()

Protected Attributes
std::vector< double >	m_exTruth
std::vector< double >	m_eyTruth
std::vector< double >	m_etSumTruth
std::unique_ptr< MissingEtRegions >	m_regions
Source	m_source
double	m_ex
double	m_ey
double	m_etSum

Detailed Description

Definition at line 27 of file MissingEtTruth.h.

Member Enumeration Documentation

Source

enum MissingET::Source

inherited

Enumerator
Calo
Calib
Truth
Muon
CaloRaw
Final
Cryo
Topo
Ref
DeadMaterial
Track
ObjElectron
ObjMuon
ObjJet
ObjIdTrk
ObjMiniJet
ObjRest
ObjFinal
Unknown

Definition at line 28 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

    { 
      Calo        = 0,
      Calib       = 1,
      Truth       = 2,
      Muon        = 3, 
      CaloRaw     = 4,
      Final       = 5,
      Cryo        = 6,
      Topo        = 7,
      Ref         = 8,
      DeadMaterial= 9,
      Track       = 10,
      ObjElectron = 100,
      ObjMuon     = 101,
      ObjJet      = 102,
      ObjIdTrk    = 103,
      ObjMiniJet  = 104,
      ObjRest     = 105,
      ObjFinal    = 106,
      Unknown     = 999
    };

TruthIndex

enum MissingEtTruth::TruthIndex

Enumerator
Int
NonInt
IntCentral
IntFwd
IntOutCover
Muons
Size
NotValid

Definition at line 31 of file MissingEtTruth.h.

   { 
     Int        = 0,  // all interacting particles till abs(eta) <= 5
     NonInt     = 1,  // non interacting particles
     IntCentral = 2,  // interacting particles in central region : 
                      // till eta of the full sim
     IntFwd     = 3,  // interacting particles in Forward region :
                      // with eta_full_sim < abs(eta) <= 5 
     IntOutCover= 4,  // interacting particles with abs(eta) > 5
     Muons      = 5,  // truth muons ( all )              
     Size       = 6,                  
     NotValid   = 999 
   };         

Constructor & Destructor Documentation

MissingEtTruth() [1/3]

MissingEtTruth::MissingEtTruth

(

)

Definition at line 24 of file MissingEtTruth.cxx.

                               : MissingET(MissingET::Truth)
{
  this->setup();
}

MissingEtTruth() [2/3]

MissingEtTruth::MissingEtTruth

(

MissingEtRegions *

aRegion

)

Definition at line 29 of file MissingEtTruth.cxx.

  : MissingET(MissingET::Truth,aRegion)
{
  this->setup();
}

MissingEtTruth() [3/3]

MissingEtTruth::MissingEtTruth	(	MissingET::Source	aSource,
		MissingEtRegions *	aRegion )

Definition at line 35 of file MissingEtTruth.cxx.

  : MissingET(aSource,aRegion)
{ 
  this->setup();
}

~MissingEtTruth()

MissingEtTruth::~MissingEtTruth ( )

virtual

Definition at line 42 of file MissingEtTruth.cxx.

{ 
  this->setup();
}

Member Function Documentation

add()

void MissingET::add ( double theEx, double theEy, double theEt, double theEta, double weight )

virtualinherited

Definition at line 127 of file MissingET.cxx.

{        
  // Note the minus sign! We want to calculate missing ET 
  m_ex -= theEx * weight;
  m_ey -= theEy * weight;
  m_etSum += theEt * weight;
    
  // regions
  if ( m_regions != 0 ) m_regions->addReg(theEx, theEy, theEt, theEta, weight);
} 

addPart()

void MissingEtTruth::addPart	(	HepMC::ConstGenParticlePtr	aPart,
		double	etaFull )

Definition at line 47 of file MissingEtTruth.cxx.

{
  
  double weight = 1.;
  double exp  = (Part)->momentum().px();
  double eyp  = (Part)->momentum().py();
  double etp  = (Part)->momentum().perp();
  double etap = (Part)->momentum().pseudoRapidity();
  int idp    = (Part)->pdg_id();
 
  TruthIndex thePart = NotValid;
 
  if ( ! MC::isSimInteracting(Part) )  // if IsGenNonInteracting returns true
  {  
    thePart = NonInt;    // all NON interacting particles (neutrinos + SUSY particles)
    m_exTruth[NonInt]    += exp;    
    m_eyTruth[NonInt]    += eyp;    
    m_etSumTruth[NonInt] += etp; 
    
  } else {                 
    if (abs(idp) != 13)  // all interacting particles and NO muons 
    {
      if( fabs(etap) <= 5.) {
        // update the global sums and eta regions in eta <= 5
        this->add(exp, eyp, etp, etap, weight);
         m_exTruth[Int]    -= exp;    
         m_eyTruth[Int]    -= eyp;    
         m_etSumTruth[Int] += etp;   
    //thePart = Int;
      }                          
      if( fabs(etap) <= etaFull)                    thePart = IntCentral;     
      if( fabs(etap) > etaFull && fabs(etap) <= 5.) thePart = IntFwd; 
      if( fabs(etap) > 5.)                          thePart = IntOutCover;       
    } else {
      //if (etp > 6000. && fabs(etap) <= 2.5)   //  select  "measurable" muons
      //{
        thePart = Muons; //  all muons
      //}
    }
    
   if (thePart != NotValid) { // Can be == if etap is a nan !! Fixed by Scott Snyder
     m_exTruth[(size_t)thePart] -= exp;    
     m_eyTruth[(size_t)thePart] -= eyp;    
     m_etSumTruth[(size_t)thePart] += etp;
   }
  
  }  
  
     
    
}

et()

double MissingET::et ( ) const

virtualinherited

Definition at line 153 of file MissingET.cxx.

{
  return sqrt(m_ex*m_ex + m_ey*m_ey);
}

etSumTruth()

double MissingEtTruth::etSumTruth

(

MissingEtTruth::TruthIndex

theTruth

)

const

Definition at line 157 of file MissingEtTruth.cxx.

{
  return m_etSumTruth[theTruth] ;
}

etSumTruthVec()

const std::vector< double > & MissingEtTruth::etSumTruthVec

(

)

const

Definition at line 175 of file MissingEtTruth.cxx.

{
  return m_etSumTruth;
}

etx()

double MissingET::etx ( ) const

virtualinherited

Definition at line 141 of file MissingET.cxx.

{ 
  return m_ex; 
}

ety()

double MissingET::ety ( ) const

virtualinherited

Definition at line 145 of file MissingET.cxx.

{ 
  return m_ey; 
}

exTruth()

double MissingEtTruth::exTruth

(

MissingEtTruth::TruthIndex

theTruth

)

const

Definition at line 147 of file MissingEtTruth.cxx.

{
  return m_exTruth[theTruth] ;
}

exTruthVec()

const std::vector< double > & MissingEtTruth::exTruthVec

(

)

const

Definition at line 163 of file MissingEtTruth.cxx.

{
  return m_exTruth;
}

eyTruth()

double MissingEtTruth::eyTruth

(

MissingEtTruth::TruthIndex

theTruth

)

const

Definition at line 152 of file MissingEtTruth.cxx.

{
  return m_eyTruth[theTruth] ;
}

eyTruthVec()

const std::vector< double > & MissingEtTruth::eyTruthVec

(

)

const

Definition at line 169 of file MissingEtTruth.cxx.

{
  return m_eyTruth;
}

getRegions()

const MissingEtRegions * MissingET::getRegions ( ) const

virtualinherited

Definition at line 164 of file MissingET.cxx.

{
  return m_regions.get();
}

getSource()

MissingET::Source MissingET::getSource ( ) const

virtualinherited

Definition at line 171 of file MissingET.cxx.

{
  return m_source;
}

phi()

double MissingET::phi ( ) const

virtualinherited

Definition at line 157 of file MissingET.cxx.

{
 
  return atan2(m_ey,m_ex);
}

setEtSum()

void MissingET::setEtSum ( double theSum )

virtualinherited

Definition at line 104 of file MissingET.cxx.

{
  m_etSum = theEtSum;
}

setEtSumTruth()

void MissingEtTruth::setEtSumTruth	(	MissingEtTruth::TruthIndex	theTruth,
		double	theEtSum )

Definition at line 116 of file MissingEtTruth.cxx.

{
  if ( theTruth < Size ) m_etSumTruth[(size_t)theTruth] = theEtSum;
  
}

setEtSumTruthVec()

void MissingEtTruth::setEtSumTruthVec

(

std::vector< double > &&

etSumCaloVec

)

Definition at line 138 of file MissingEtTruth.cxx.

{
  assert (etSumTruthVec.size() == Size);
  m_etSumTruth = std::move (etSumTruthVec);
}

setEx()

void MissingET::setEx ( double theEx )

virtualinherited

Definition at line 96 of file MissingET.cxx.

{
  m_ex = theEx;
}

setExTruth()

void MissingEtTruth::setExTruth	(	MissingEtTruth::TruthIndex	theTruth,
		double	theEx )

Definition at line 102 of file MissingEtTruth.cxx.

{
  if ( theTruth < Size ) m_exTruth[(size_t)theTruth] = theEx;
}

setExTruthVec()

void MissingEtTruth::setExTruthVec

(

std::vector< double > &&

exCaloVec

)

Definition at line 124 of file MissingEtTruth.cxx.

{
  assert (exTruthVec.size() == Size);
  m_exTruth = std::move (exTruthVec);
}

setEy()

void MissingET::setEy ( double theEy )

virtualinherited

Definition at line 100 of file MissingET.cxx.

{
  m_ey = theEy;
}

setEyTruth()

void MissingEtTruth::setEyTruth	(	MissingEtTruth::TruthIndex	theTruth,
		double	theEy )

Definition at line 109 of file MissingEtTruth.cxx.

{
  if ( theTruth < Size ) m_eyTruth[(size_t)theTruth] = theEy;
}

setEyTruthVec()

void MissingEtTruth::setEyTruthVec

(

std::vector< double > &&

exCaloVec

)

Definition at line 131 of file MissingEtTruth.cxx.

{
  assert (eyTruthVec.size() == Size);
  m_eyTruth = std::move (eyTruthVec);
}

setup()

void MissingEtTruth::setup ( )

protected

Definition at line 181 of file MissingEtTruth.cxx.

{
  // initialize data members
  
  m_exTruth.resize(MissingEtTruth::Size,0.);
  m_eyTruth.resize(MissingEtTruth::Size,0.);
  m_etSumTruth.resize(MissingEtTruth::Size,0.);
}

sumet()

double MissingET::sumet ( ) const

virtualinherited

Definition at line 149 of file MissingET.cxx.

{ 
  return m_etSum; 
}

Member Data Documentation

m_etSum

double MissingET::m_etSum

protectedinherited

Definition at line 115 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

m_etSumTruth

std::vector<double> MissingEtTruth::m_etSumTruth

protected

Definition at line 79 of file MissingEtTruth.h.

m_ex

double MissingET::m_ex

protectedinherited

Definition at line 113 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

m_exTruth

std::vector<double> MissingEtTruth::m_exTruth

protected

Definition at line 77 of file MissingEtTruth.h.

m_ey

double MissingET::m_ey

protectedinherited

Definition at line 114 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

m_eyTruth

std::vector<double> MissingEtTruth::m_eyTruth

protected

Definition at line 78 of file MissingEtTruth.h.

m_regions

std::unique_ptr<MissingEtRegions> MissingET::m_regions

protectedinherited

Definition at line 108 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

m_source

Source MissingET::m_source

protectedinherited

Definition at line 110 of file Reconstruction/MissingETEvent/MissingETEvent/MissingET.h.

---

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

• MissingEtTruth.h
• MissingEtTruth.cxx