CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
HLTMCtruth Class Reference

#include <HLTMCtruth.h>

Public Member Functions

void analyze (const edm::Handle< reco::CandidateView > &mctruth, const double &pthat, const edm::Handle< std::vector< SimTrack > > &simTracks, const edm::Handle< std::vector< SimVertex > > &simVertices, TTree *tree)
 
 HLTMCtruth ()
 
void setup (const edm::ParameterSet &pSet, TTree *tree)
 

Private Attributes

bool _Debug
 
bool _Monte
 
float * mceta
 
float * mcphi
 
int * mcpid
 
float * mcpt
 
int * mcstatus
 
float * mcvx
 
float * mcvy
 
float * mcvz
 
int nab
 
int nbb
 
int nel3
 
int nmcpart
 
int nmu3
 
int nwenu
 
int nwmunu
 
int nzee
 
int nzmumu
 
float ptEleMax
 
float pthatf
 
float ptMuMax
 

Detailed Description

$Date: November 2006 $Revision:

Author
P. Bargassa - Rice U.

Definition at line 29 of file HLTMCtruth.h.

Constructor & Destructor Documentation

HLTMCtruth::HLTMCtruth ( )

Definition at line 15 of file HLTMCtruth.cc.

References _Debug, and _Monte.

15  {
16 
17  //set parameter defaults
18  _Monte=false;
19  _Debug=false;
20 }
HLTMCtruth::_Debug
bool _Debug
Definition: HLTMCtruth.h:51
HLTMCtruth::_Monte
bool _Monte
Definition: HLTMCtruth.h:51

Member Function Documentation

void HLTMCtruth::analyze ( const edm::Handle< reco::CandidateView > &  mctruth,
const double &  pthat,
const edm::Handle< std::vector< SimTrack > > &  simTracks,
const edm::Handle< std::vector< SimVertex > > &  simVertices,
TTree *  tree 
)

Analyze the Data

Definition at line 69 of file HLTMCtruth.cc.

References _Monte, funct::abs(), ztail::d, reco::Candidate::daughter(), eta, reco::Candidate::eta(), i, edm::HandleBase::isValid(), j, mceta, mcphi, mcpid, mcpt, mcstatus, mcvx, mcvy, mcvz, nab, nbb, nel3, nmcpart, nmu3, reco::Candidate::numberOfDaughters(), nwenu, nwmunu, nzee, nzmumu, AlCaHLTBitMon_ParallelJobs::p, reco::Candidate::pdgId(), reco::Candidate::phi(), EnergyCorrector::pt, reco::Candidate::pt(), ptEleMax, pthatf, ptMuMax, alignCSCRings::r, mathSSE::sqrt(), reco::Candidate::status(), reco::Candidate::vx(), reco::Candidate::vy(), reco::Candidate::vz(), x, y, and z.

Referenced by HLTAnalyzer::analyze().

73  {
74 
75  //std::cout << " Beginning HLTMCtruth " << std::endl;
76 
77  if (_Monte) {
78  int nmc = 0;
79  int mu3 = 0;
80  int el3 = 0;
81  int mab = 0;
82  int mbb = 0;
83  int wel = 0;
84  int wmu = 0;
85  int zee = 0;
86  int zmumu = 0;
87 
88  ptEleMax = -999.0;
89  ptMuMax = -999.0;
90  pthatf = pthat;
91 
92  if((simTracks.isValid())&&(simVertices.isValid())){
93  for (unsigned int j=0; j<simTracks->size(); j++) {
94  int pdgid = simTracks->at(j).type();
95  if (abs(pdgid)!=13) continue;
96  double pt = simTracks->at(j).momentum().pt();
97  if (pt<2.5) continue;
98  double eta = simTracks->at(j).momentum().eta();
99  if (abs(eta)>2.5) continue;
100  if (simTracks->at(j).noVertex()) continue;
101  int vertIndex = simTracks->at(j).vertIndex();
102  double x = simVertices->at(vertIndex).position().x();
103  double y = simVertices->at(vertIndex).position().y();
104  double r = sqrt(x*x+y*y);
105  if (r>150.) continue; // I think units are cm here
106  double z = simVertices->at(vertIndex).position().z();
107  if (abs(z)>300.) continue; // I think units are cm here
108  mu3 += 1;
109  break;
110  }
111  }
112 
113  if (mctruth.isValid()){
114 
115  for (size_t i = 0; i < mctruth->size(); ++ i) {
116  const reco::Candidate & p = (*mctruth)[i];
117 
118  mcpid[nmc] = p.pdgId();
119  mcstatus[nmc] = p.status();
120  mcpt[nmc] = p.pt();
121  mceta[nmc] = p.eta();
122  mcphi[nmc] = p.phi();
123  mcvx[nmc] = p.vx();
124  mcvy[nmc] = p.vy();
125  mcvz[nmc] = p.vz();
126 
127  if ((mcpid[nmc]==24)||(mcpid[nmc]==-24)) { // Checking W -> e/mu nu
128  size_t idg = p.numberOfDaughters();
129  for (size_t j=0; j != idg; ++j){
130  const reco::Candidate & d = *p.daughter(j);
131  if ((d.pdgId()==11)||(d.pdgId()==-11)){wel += 1;}
132  if ((d.pdgId()==13)||(d.pdgId()==-13)){wmu += 1;}
133 // if ( (abs(d.pdgId())!=24) && ((mcpid[nmc])*(d.pdgId())>0) )
134 // {std::cout << "Wrong sign between mother-W and daughter !" << std::endl;}
135  }
136  }
137  if (mcpid[nmc]==23) { // Checking Z -> 2 e/mu
138  size_t idg = p.numberOfDaughters();
139  for (size_t j=0; j != idg; ++j){
140  const reco::Candidate & d = *p.daughter(j);
141  if (d.pdgId()==11){zee += 1;}
142  if (d.pdgId()==-11){zee += 2;}
143  if (d.pdgId()==13){zmumu += 1;}
144  if (d.pdgId()==-13){zmumu += 2;}
145  }
146  }
147 
148  // Set-up flags, based on Pythia-generator information, for avoiding double-counting events when
149  // using both pp->{e,mu}X AND QCD samples
150 // if (((mcpid[nmc]==13)||(mcpid[nmc]==-13))&&(mcpt[nmc]>2.5)) {mu3 += 1;} // Flag for muons with pT > 2.5 GeV/c
151  if (((mcpid[nmc]==11)||(mcpid[nmc]==-11))&&(mcpt[nmc]>2.5)) {el3 += 1;} // Flag for electrons with pT > 2.5 GeV/c
152 
153  if (mcpid[nmc]==-5) {mab += 1;} // Flag for bbar
154  if (mcpid[nmc]==5) {mbb += 1;} // Flag for b
155 
156  if ((mcpid[nmc]==13)||(mcpid[nmc]==-13))
157  {if (p.pt()>ptMuMax) {ptMuMax=p.pt();} } // Save max pt of generated Muons
158  if ((mcpid[nmc]==11)||(mcpid[nmc]==-11))
159  {if (p.pt() > ptEleMax) ptEleMax=p.pt();} // Save max pt of generated Electrons
160 
161  nmc++;
162  }
163 
164  }
165  // else {std::cout << "%HLTMCtruth -- No MC truth information" << std::endl;}
166 
167  nmcpart = nmc;
168  nmu3 = mu3;
169  nel3 = el3;
170  nbb = mbb;
171  nab = mab;
172  nwenu = wel;
173  nwmunu = wmu;
174  if((zee%3)==0){nzee = zee/3;}
175 // else {std::cout << "Z does not decay in e+ e- !" << std::endl;}
176  if ((zmumu%3)==0){nzmumu = zmumu/3;}
177 // else {std::cout << "Z does not decay in mu+ mu- !" << std::endl;}
178 
179  }
180 
181 }
i
int i
Definition: DBlmapReader.cc:9
HLTMCtruth::mcpid
int * mcpid
Definition: HLTMCtruth.h:46
reco::Candidate::daughter
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
HLTMCtruth::mcpt
float * mcpt
Definition: HLTMCtruth.h:45
reco::Candidate::pt
virtual double pt() const =0
transverse momentum
HLTMCtruth::pthatf
float pthatf
Definition: HLTMCtruth.h:48
reco::Candidate::status
virtual int status() const =0
status word
reco::Candidate::vx
virtual double vx() const =0
x coordinate of vertex position
reco::Candidate::vy
virtual double vy() const =0
y coordinate of vertex position
HLTMCtruth::mceta
float * mceta
Definition: HLTMCtruth.h:45
HLTMCtruth::nmcpart
int nmcpart
Definition: HLTMCtruth.h:47
ztail.d
tuple d
Definition: ztail.py:151
HLTMCtruth::nzmumu
int nzmumu
Definition: HLTMCtruth.h:47
reco::Candidate::numberOfDaughters
virtual size_type numberOfDaughters() const =0
number of daughters
HLTMCtruth::ptEleMax
float ptEleMax
Definition: HLTMCtruth.h:49
HLTMCtruth::ptMuMax
float ptMuMax
Definition: HLTMCtruth.h:49
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
HLTMCtruth::mcvy
float * mcvy
Definition: HLTMCtruth.h:45
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
j
int j
Definition: DBlmapReader.cc:9
HLTMCtruth::_Monte
bool _Monte
Definition: HLTMCtruth.h:51
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
HLTMCtruth::mcvx
float * mcvx
Definition: HLTMCtruth.h:45
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
HLTMCtruth::mcphi
float * mcphi
Definition: HLTMCtruth.h:45
reco::Candidate::vz
virtual double vz() const =0
z coordinate of vertex position
HLTMCtruth::mcvz
float * mcvz
Definition: HLTMCtruth.h:45
HLTMCtruth::nwmunu
int nwmunu
Definition: HLTMCtruth.h:47
HLTMCtruth::mcstatus
int * mcstatus
Definition: HLTMCtruth.h:46
reco::Candidate::phi
virtual double phi() const =0
momentum azimuthal angle
reco::Candidate::eta
virtual double eta() const =0
momentum pseudorapidity
void HLTMCtruth::setup ( const edm::ParameterSet &  pSet,
TTree *  tree 
)

Definition at line 23 of file HLTMCtruth.cc.

References _Debug, _Monte, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), mceta, mcphi, mcpid, mcpt, mcstatus, mcvx, mcvy, mcvz, nab, nbb, nel3, nmcpart, nmu3, nwenu, nwmunu, nzee, nzmumu, ptEleMax, pthatf, and ptMuMax.

Referenced by HLTAnalyzer::HLTAnalyzer().

23  {
24 
25  edm::ParameterSet myMCParams = pSet.getParameter<edm::ParameterSet>("RunParameters") ;
26  std::vector<std::string> parameterNames = myMCParams.getParameterNames() ;
27 
28  for ( std::vector<std::string>::iterator iParam = parameterNames.begin();
29  iParam != parameterNames.end(); iParam++ ){
30  if ( (*iParam) == "Monte" ) _Monte = myMCParams.getParameter<bool>( *iParam );
31  else if ( (*iParam) == "Debug" ) _Debug = myMCParams.getParameter<bool>( *iParam );
32  }
33 
34  const int kMaxMcTruth = 10000;
35  mcpid = new int[kMaxMcTruth];
36  mcstatus = new int[kMaxMcTruth];
37  mcvx = new float[kMaxMcTruth];
38  mcvy = new float[kMaxMcTruth];
39  mcvz = new float[kMaxMcTruth];
40  mcpt = new float[kMaxMcTruth];
41  mceta = new float[kMaxMcTruth];
42  mcphi = new float[kMaxMcTruth];
43 
44  // MCtruth-specific branches of the tree
45  HltTree->Branch("NMCpart",&nmcpart,"NMCpart/I");
46  HltTree->Branch("MCpid",mcpid,"MCpid[NMCpart]/I");
47  HltTree->Branch("MCstatus",mcstatus,"MCstatus[NMCpart]/I");
48  HltTree->Branch("MCvtxX",mcvx,"MCvtxX[NMCpart]/F");
49  HltTree->Branch("MCvtxY",mcvy,"MCvtxY[NMCpart]/F");
50  HltTree->Branch("MCvtxZ",mcvz,"MCvtxZ[NMCpart]/F");
51  HltTree->Branch("MCpt",mcpt,"MCpt[NMCpart]/F");
52  HltTree->Branch("MCeta",mceta,"MCeta[NMCpart]/F");
53  HltTree->Branch("MCphi",mcphi,"MCphi[NMCpart]/F");
54  HltTree->Branch("MCPtHat",&pthatf,"MCPtHat/F");
55  HltTree->Branch("MCmu3",&nmu3,"MCmu3/I");
56  HltTree->Branch("MCel3",&nel3,"MCel3/I");
57  HltTree->Branch("MCbb",&nbb,"MCbb/I");
58  HltTree->Branch("MCab",&nab,"MCab/I");
59  HltTree->Branch("MCWenu",&nwenu,"MCWenu/I");
60  HltTree->Branch("MCWmunu",&nwmunu,"MCmunu/I");
61  HltTree->Branch("MCZee",&nzee,"MCZee/I");
62  HltTree->Branch("MCZmumu",&nzmumu,"MCZmumu/I");
63  HltTree->Branch("MCptEleMax",&ptEleMax,"MCptEleMax/F");
64  HltTree->Branch("MCptMuMax",&ptMuMax,"MCptMuMax/F");
65 
66 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
HLTMCtruth::mcpid
int * mcpid
Definition: HLTMCtruth.h:46
HLTMCtruth::mcpt
float * mcpt
Definition: HLTMCtruth.h:45
HLTMCtruth::pthatf
float pthatf
Definition: HLTMCtruth.h:48
HLTMCtruth::_Debug
bool _Debug
Definition: HLTMCtruth.h:51
HLTMCtruth::mceta
float * mceta
Definition: HLTMCtruth.h:45
HLTMCtruth::nmcpart
int nmcpart
Definition: HLTMCtruth.h:47
HLTMCtruth::nzmumu
int nzmumu
Definition: HLTMCtruth.h:47
HLTMCtruth::ptEleMax
float ptEleMax
Definition: HLTMCtruth.h:49
HLTMCtruth::ptMuMax
float ptMuMax
Definition: HLTMCtruth.h:49
HLTMCtruth::mcvy
float * mcvy
Definition: HLTMCtruth.h:45
HLTMCtruth::_Monte
bool _Monte
Definition: HLTMCtruth.h:51
edm::ParameterSet::getParameterNames
std::vector< std::string > getParameterNames() const
Definition: ParameterSet.cc:748
HLTMCtruth::mcvx
float * mcvx
Definition: HLTMCtruth.h:45
HLTMCtruth::mcphi
float * mcphi
Definition: HLTMCtruth.h:45
HLTMCtruth::mcvz
float * mcvz
Definition: HLTMCtruth.h:45
HLTMCtruth::nwmunu
int nwmunu
Definition: HLTMCtruth.h:47
HLTMCtruth::mcstatus
int * mcstatus
Definition: HLTMCtruth.h:46

Member Data Documentation

bool HLTMCtruth::_Debug
private

Definition at line 51 of file HLTMCtruth.h.

Referenced by HLTMCtruth(), and setup().

bool HLTMCtruth::_Monte
private

Definition at line 51 of file HLTMCtruth.h.

Referenced by analyze(), HLTMCtruth(), and setup().

float * HLTMCtruth::mceta
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float * HLTMCtruth::mcphi
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int* HLTMCtruth::mcpid
private

Definition at line 46 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float * HLTMCtruth::mcpt
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int * HLTMCtruth::mcstatus
private

Definition at line 46 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float* HLTMCtruth::mcvx
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float * HLTMCtruth::mcvy
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float * HLTMCtruth::mcvz
private

Definition at line 45 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nab
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nbb
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nel3
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nmcpart
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nmu3
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nwenu
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nwmunu
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nzee
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

int HLTMCtruth::nzmumu
private

Definition at line 47 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float HLTMCtruth::ptEleMax
private

Definition at line 49 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float HLTMCtruth::pthatf
private

Definition at line 48 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

float HLTMCtruth::ptMuMax
private

Definition at line 49 of file HLTMCtruth.h.

Referenced by analyze(), and setup().

Generated for CMSSW Reference Manual by doxygen 1.8.5