VHbbDataFormats/interface/VHbbEventAuxInfo.h

#ifndef VHbbEventAuxInfo__H 
#define VHbbEventAuxInfo__H 

#include <TLorentzVector.h>
#include <TVector2.h>
#include "Math/SMatrix.h"
#include <vector>
#include<map>
class VHbbEventAuxInfo{
 public:


  class PrimaryVertexInfo {
  public:
    PrimaryVertexInfo() : nVertices (-99){}
    TVector3 firstPVInPT2,firstPVInProb;
    ROOT::Math::SMatrix<double, 3, 3, ROOT::Math::MatRepSym<double, 3> > efirstPVInPT2,efirstPVInProb;
    int nVertices;
  };

  class PUInfo{
  public:
    PUInfo(): rho(-99), rho25(-99),rho25Iso(-99),rhoNeutral(-99),truePU(9999) {}
  public:
    float rho, rho25,rhoNeutral,rho25Iso;
    std::map<int,unsigned int>  pus;
    unsigned int truePU;
  };


  class ParticleMCInfo {
  public:
    ParticleMCInfo(): status(-99), momid(-99), gmomid(-99), charge(-99){}
  public:
    int status;
    int momid;
    int gmomid;
    float charge;
    TLorentzVector p4;
    //    int ndau;
    std::vector<int> dauid;
    std::vector<TLorentzVector> dauFourMomentum;    
  };

 double genBBDeltaR() const
 {
   if(mcB.size() > 0 && mcBbar.size() > 0)
   return mcB[0].p4.DeltaR(mcBbar[0].p4);
   else return -99;
 }
 double genCCDeltaR() const
 {
   if(mcC.size() >=2)
   return mcC[0].p4.DeltaR(mcC[1].p4);
   else return -99;
 }
 
 public:

    int run;
    int event; 
    int lumi; 
    int bunch; 
    
  PUInfo puInfo;
  float weightMCProd;  
 PrimaryVertexInfo pvInfo;
  std::vector<ParticleMCInfo> mcH;
  std::vector<ParticleMCInfo> mcW;
  std::vector<ParticleMCInfo> mcZ;
  
  std::vector<ParticleMCInfo> mcBbar;
  std::vector<ParticleMCInfo> mcB;
  std::vector<ParticleMCInfo> mcC;       

  std::vector<ParticleMCInfo> mcTau;
};
#endif
Revision:	1.16.2.1
Committed:	Wed Nov 28 02:48:02 2012 UTC (12 years, 5 months ago) by ntran
Content type:	text/plain
Branch:	hbbsubstructDevPostHCP
CVS Tags:	hbbsubstructDev_11, hbbsubstructDev_10, hbbsubstructDev_9, hbbsubstructDev_8, hbbsubstructDev_7
Changes since 1.16:	+6 -0 lines
Log Message:	move substructure tools to a new area
#	User	Rev	Content
1	tboccali	1.1	#ifndef VHbbEventAuxInfo__H
2			#define VHbbEventAuxInfo__H
3
4			#include <TLorentzVector.h>
5			#include <TVector2.h>
6	tboccali	1.12	#include "Math/SMatrix.h"
7	tboccali	1.1	#include <vector>
8	tboccali	1.9	#include<map>
9	tboccali	1.1	class VHbbEventAuxInfo{
10			public:
11
12
13			class PrimaryVertexInfo {
14			public:
15	tboccali	1.5	PrimaryVertexInfo() : nVertices (-99){}
16	tboccali	1.1	TVector3 firstPVInPT2,firstPVInProb;
17	tboccali	1.12	ROOT::Math::SMatrix<double, 3, 3, ROOT::Math::MatRepSym<double, 3> > efirstPVInPT2,efirstPVInProb;
18	tboccali	1.5	int nVertices;
19	tboccali	1.1	};
20
21			class PUInfo{
22			public:
23	arizzi	1.16	PUInfo(): rho(-99), rho25(-99),rho25Iso(-99),rhoNeutral(-99),truePU(9999) {}
24	tboccali	1.1	public:
25	arizzi	1.16	float rho, rho25,rhoNeutral,rho25Iso;
26	tboccali	1.9	std::map<int,unsigned int> pus;
27	arizzi	1.15	unsigned int truePU;
28	tboccali	1.1	};
29
30
31	tboccali	1.2	class ParticleMCInfo {
32			public:
33	tboccali	1.4	ParticleMCInfo(): status(-99), momid(-99), gmomid(-99), charge(-99){}
34	tboccali	1.2	public:
35			int status;
36			int momid;
37			int gmomid;
38	tboccali	1.4	float charge;
39	tboccali	1.3	TLorentzVector p4;
40	tboccali	1.2	// int ndau;
41			std::vector<int> dauid;
42			std::vector<TLorentzVector> dauFourMomentum;
43			};
44
45	arizzi	1.7	double genBBDeltaR() const
46	arizzi	1.6	{
47	arizzi	1.7	if(mcB.size() > 0 && mcBbar.size() > 0)
48	arizzi	1.6	return mcB[0].p4.DeltaR(mcBbar[0].p4);
49	arizzi	1.7	else return -99;
50			}
51			double genCCDeltaR() const
52			{
53			if(mcC.size() >=2)
54			return mcC[0].p4.DeltaR(mcC[1].p4);
55			else return -99;
56	arizzi	1.6	}
57
58	tboccali	1.1	public:
59	ntran	1.16.2.1
60			int run;
61			int event;
62			int lumi;
63			int bunch;
64
65	tboccali	1.1	PUInfo puInfo;
66	arizzi	1.10	float weightMCProd;
67	tboccali	1.8	PrimaryVertexInfo pvInfo;
68	tboccali	1.2	std::vector<ParticleMCInfo> mcH;
69			std::vector<ParticleMCInfo> mcW;
70			std::vector<ParticleMCInfo> mcZ;
71	tboccali	1.1
72	tboccali	1.2	std::vector<ParticleMCInfo> mcBbar;
73			std::vector<ParticleMCInfo> mcB;
74			std::vector<ParticleMCInfo> mcC;
75	tboccali	1.1
76	sethzenz	1.13	std::vector<ParticleMCInfo> mcTau;
77	tboccali	1.1	};
78			#endif