VHbbDataFormats/interface/VHbbEventAuxInfo.h

#ifndef VHbbEventAuxInfo__H 
#define VHbbEventAuxInfo__H 

#include <TLorentzVector.h>
#include <TVector2.h>
#include <vector>

class VHbbEventAuxInfo{
 public:


  class PrimaryVertexInfo {
  public:
    PrimaryVertexInfo() : nVertices (-99){}
    TVector3 firstPVInPT2,firstPVInProb;
    int nVertices;
  };

  class PUInfo{
  public:
    PUInfo(): rho(-99) {}
  public:
    float rho;
  };


  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:
  PUInfo puInfo;
  
 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;       

};
#endif
Revision:	1.8
Committed:	Fri Sep 9 08:05:27 2011 UTC (13 years, 7 months ago) by tboccali
Content type:	text/plain
Branch:	MAIN
Changes since 1.7:	+1 -66 lines
Log Message:	many small changes
#	User	Rev	Content
1	tboccali	1.1	#ifndef VHbbEventAuxInfo__H
2			#define VHbbEventAuxInfo__H
3
4			#include <TLorentzVector.h>
5			#include <TVector2.h>
6			#include <vector>
7
8			class VHbbEventAuxInfo{
9			public:
10
11
12			class PrimaryVertexInfo {
13			public:
14	tboccali	1.5	PrimaryVertexInfo() : nVertices (-99){}
15	tboccali	1.1	TVector3 firstPVInPT2,firstPVInProb;
16	tboccali	1.5	int nVertices;
17	tboccali	1.1	};
18
19			class PUInfo{
20			public:
21			PUInfo(): rho(-99) {}
22			public:
23			float rho;
24			};
25
26
27	tboccali	1.2	class ParticleMCInfo {
28			public:
29	tboccali	1.4	ParticleMCInfo(): status(-99), momid(-99), gmomid(-99), charge(-99){}
30	tboccali	1.2	public:
31			int status;
32			int momid;
33			int gmomid;
34	tboccali	1.4	float charge;
35	tboccali	1.3	TLorentzVector p4;
36	tboccali	1.2	// int ndau;
37			std::vector<int> dauid;
38			std::vector<TLorentzVector> dauFourMomentum;
39			};
40
41	arizzi	1.7	double genBBDeltaR() const
42	arizzi	1.6	{
43	arizzi	1.7	if(mcB.size() > 0 && mcBbar.size() > 0)
44	arizzi	1.6	return mcB[0].p4.DeltaR(mcBbar[0].p4);
45	arizzi	1.7	else return -99;
46			}
47			double genCCDeltaR() const
48			{
49			if(mcC.size() >=2)
50			return mcC[0].p4.DeltaR(mcC[1].p4);
51			else return -99;
52	arizzi	1.6	}
53
54	tboccali	1.1	public:
55			PUInfo puInfo;
56
57	tboccali	1.8	PrimaryVertexInfo pvInfo;
58	tboccali	1.2	std::vector<ParticleMCInfo> mcH;
59			std::vector<ParticleMCInfo> mcW;
60			std::vector<ParticleMCInfo> mcZ;
61	tboccali	1.1
62	tboccali	1.2	std::vector<ParticleMCInfo> mcBbar;
63			std::vector<ParticleMCInfo> mcB;
64			std::vector<ParticleMCInfo> mcC;
65	tboccali	1.1
66			};
67			#endif