| 3 |
|
|
| 4 |
|
#include <string> |
| 5 |
|
#include <iostream> |
| 6 |
+ |
#include <vector> |
| 7 |
+ |
#include <map> |
| 8 |
+ |
#include <cmath> |
| 9 |
|
|
| 10 |
|
#include "Rtypes.h" |
| 11 |
|
#include "TObject.h" |
| 23 |
|
,passGlobalL1_(false) |
| 24 |
|
,passGlobalHLT_(false) |
| 25 |
|
,trigHLT_(0) |
| 26 |
+ |
,csa07id_(-1) |
| 27 |
+ |
,csa07weight_(-1.) |
| 28 |
+ |
,csa07process_() |
| 29 |
|
,nBasicClusters_() |
| 30 |
|
,nSuperClusters_() |
| 31 |
|
,primaryVertex_(0) |
| 32 |
< |
,metCalo_() |
| 32 |
> |
,idParton1_(-1) |
| 33 |
> |
,xParton1_(-1.) |
| 34 |
> |
,idParton2_(-1) |
| 35 |
> |
,xParton2_(-1.) |
| 36 |
> |
,factorizationScale_(-1.) |
| 37 |
|
{;} |
| 38 |
< |
|
| 38 |
> |
|
| 39 |
|
~TRootEvent() {;} |
| 40 |
|
|
| 41 |
|
// Event number |
| 59 |
|
} |
| 60 |
|
} |
| 61 |
|
|
| 62 |
+ |
// CSA07 Process Id and Event Weight |
| 63 |
+ |
Int_t csa07id() const { return csa07id_; } |
| 64 |
+ |
Float_t csa07weight() const { return csa07weight_; } |
| 65 |
+ |
TString csa07process() const { return csa07process_; } |
| 66 |
+ |
|
| 67 |
|
// Nb of primary vertices |
| 68 |
|
unsigned nPrimaryVertices() const { return primaryVertex_.size(); } |
| 69 |
|
|
| 93 |
|
Double_t primaryVertex_y(unsigned i) const { return (primaryVertex_.size()>i ? primaryVertex_.at(i).x() : 0. ); } |
| 94 |
|
Double_t primaryVertex_z(unsigned i) const { return (primaryVertex_.size()>i ? primaryVertex_.at(i).x() : 0. ); } |
| 95 |
|
|
| 96 |
< |
// Missing Et as 3-vector |
| 97 |
< |
TVector3 metCaloVect() const { return metCalo_; } |
| 98 |
< |
// Missing Et magnitude |
| 99 |
< |
Double_t metCalo() const { return sqrt(metCalo_.Perp2()); } |
| 100 |
< |
|
| 96 |
> |
// PDF infos |
| 97 |
> |
// flavour first incoming parton |
| 98 |
> |
Int_t idParton1() const { return idParton1_; } |
| 99 |
> |
// energy fraction carried by first incoming parton |
| 100 |
> |
Float_t xParton1() const { return xParton1_; } |
| 101 |
> |
// flavour first incoming parton |
| 102 |
> |
Int_t idParton2() const { return idParton2_; } |
| 103 |
> |
// energy fraction carried by first incoming parton |
| 104 |
> |
Float_t xParton2() const { return xParton2_; } |
| 105 |
> |
// Factorization Scale Q |
| 106 |
> |
Float_t factorizationScale() const { return factorizationScale_; } |
| 107 |
> |
|
| 108 |
> |
|
| 109 |
|
|
| 110 |
|
void setNb(Int_t nb) { nb_ = nb; } |
| 111 |
|
void setGlobalL1(Int_t passGlobalL1) { passGlobalL1_ = passGlobalL1; } |
| 112 |
< |
void setGlobalHLT(Bool_t passGlobalHLT) { passGlobalHLT_ = passGlobalHLT_; } |
| 112 |
> |
void setGlobalHLT(Bool_t passGlobalHLT) { passGlobalHLT_ = passGlobalHLT; } |
| 113 |
|
void setTrigHLT(std::vector<Bool_t> trigHLT) |
| 114 |
|
{ |
| 115 |
|
trigHLT_.resize(trigHLT.size()); |
| 116 |
|
for (unsigned int i=0; i!=trigHLT.size(); ++i) trigHLT_[i]=trigHLT[i]; |
| 117 |
|
} |
| 118 |
|
|
| 119 |
< |
|
| 119 |
> |
void setCsa07id(Int_t csa07id) { csa07id_=csa07id; } |
| 120 |
> |
void setCsa07weight(Float_t csa07weight) { csa07weight_=csa07weight; } |
| 121 |
> |
void setCsa07process(TString csa07process) { csa07process_=csa07process; } |
| 122 |
> |
void setCsa07process(char* csa07process) { csa07process_=csa07process; } |
| 123 |
|
|
| 124 |
|
void setNBasicClusters(Int_t type,Int_t nBC) { nBasicClusters_[type]=nBC; } |
| 125 |
|
void setNSuperClusters(Int_t type,Int_t nSC) { nSuperClusters_[type]=nSC; } |
| 126 |
|
|
| 127 |
|
void addPrimaryVertex(TVector3 vertex) { primaryVertex_.push_back(vertex); } |
| 128 |
|
void clearPrimaryVertex(TVector3 vertex) { primaryVertex_.clear(); } |
| 129 |
< |
void setMetCaloVect(TVector3 metCaloVect) { metCalo_ = metCaloVect; } |
| 130 |
< |
|
| 129 |
> |
|
| 130 |
> |
void setIdParton1(Int_t idParton1) { idParton1_=idParton1; } |
| 131 |
> |
void setXParton1(Float_t xParton1) { xParton1_=xParton1; } |
| 132 |
> |
void setIdParton2(Int_t idParton2) { idParton2_=idParton2; } |
| 133 |
> |
void setXParton2(Float_t xParton2) { xParton2_=xParton2; } |
| 134 |
> |
void setFactorizationScale(Float_t factorizationScale) { factorizationScale_=factorizationScale; } |
| 135 |
> |
|
| 136 |
|
/* |
| 137 |
+ |
// FIXME |
| 138 |
|
friend std::ostream& operator<< (std::ostream& stream, const TRootEvent& event) { |
| 139 |
|
stream << "Event #"<< event.nb() <<" L1="<< event.trigL1() <<" HLT="<< event.trigHLT() |
| 140 |
< |
<< " Primary vertex x=" << event.primaryVertex_x() << " y=" << event.primaryVertex_y() << " z=" << event.primaryVertex_z() |
| 109 |
< |
<< " CaloMET=" << event.metCalo() << endl; |
| 140 |
> |
<< " Primary vertex x=" << event.primaryVertex_x() << " y=" << event.primaryVertex_y() << " z=" << event.primaryVertex_z(); |
| 141 |
|
return stream; |
| 142 |
|
}; |
| 143 |
|
*/ |
| 144 |
|
|
| 145 |
|
private: |
| 146 |
< |
|
| 146 |
> |
|
| 147 |
|
Int_t nb_; |
| 148 |
+ |
|
| 149 |
+ |
// Trigger Infos |
| 150 |
|
Bool_t passGlobalL1_; |
| 151 |
|
Bool_t passGlobalHLT_; |
| 152 |
|
std::vector<Bool_t> trigHLT_; |
| 153 |
+ |
|
| 154 |
+ |
// CSA07 Process ID and Weight |
| 155 |
+ |
Int_t csa07id_; |
| 156 |
+ |
Float_t csa07weight_; |
| 157 |
+ |
TString csa07process_; |
| 158 |
|
|
| 159 |
|
map<Int_t,Int_t> nBasicClusters_; |
| 160 |
|
map<Int_t,Int_t> nSuperClusters_; |
| 161 |
|
std::vector<TVector3> primaryVertex_; |
| 162 |
|
|
| 163 |
< |
TVector3 metCalo_; |
| 164 |
< |
|
| 165 |
< |
ClassDef (TRootEvent,2); |
| 163 |
> |
// PDF infos |
| 164 |
> |
Int_t idParton1_; |
| 165 |
> |
Float_t xParton1_; |
| 166 |
> |
Int_t idParton2_; |
| 167 |
> |
Float_t xParton2_; |
| 168 |
> |
Float_t factorizationScale_; |
| 169 |
> |
|
| 170 |
> |
ClassDef (TRootEvent,1); |
| 171 |
|
}; |
| 172 |
|
|
| 173 |
|
#endif |
