| 1 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
| 2 |
sixie |
1.7 |
// $Id: EvtSelData.h,v 1.6 2010/01/07 17:09:02 loizides Exp $
|
| 3 |
loizides |
1.1 |
//
|
| 4 |
|
|
// EvtSelData
|
| 5 |
|
|
//
|
| 6 |
|
|
// Class to store event properties on which one may want to select clean events.
|
| 7 |
|
|
// This may be needed for early collision data only.
|
| 8 |
|
|
//
|
| 9 |
|
|
// Authors: C.Loizides
|
| 10 |
|
|
//--------------------------------------------------------------------------------------------------
|
| 11 |
|
|
|
| 12 |
|
|
#ifndef MITANA_DATATREE_EVENTSELDATA_H
|
| 13 |
|
|
#define MITANA_DATATREE_EVENTSELDATA_H
|
| 14 |
|
|
|
| 15 |
|
|
#include "MitAna/DataTree/interface/DataBase.h"
|
| 16 |
|
|
|
| 17 |
|
|
namespace mithep
|
| 18 |
|
|
{
|
| 19 |
|
|
class EvtSelData : public DataBase
|
| 20 |
|
|
{
|
| 21 |
|
|
public:
|
| 22 |
|
|
EvtSelData() : fHcalNeg(0), fHcalPos(0),
|
| 23 |
|
|
fHfNeg(0), fHfPos(0), fHfNegTime(0), fHfPosTime(0),
|
| 24 |
|
|
fCaNeg(0), fCaPos(0), fCaNegTime(0), fCaPosTime(0),
|
| 25 |
|
|
fZdcNeg(0), fZdcPos(0), fZdcNegTime(0), fZdcPosTime(0),
|
| 26 |
edwenger |
1.5 |
fPxbHits(0), fPxHits(0), fClusVtxQual(0), fClusVtxDiff(0),
|
| 27 |
loizides |
1.6 |
fNHfNegHits(0), fNHfPosHits(0), fNHfTowersP(0), fNHfTowersN(0),
|
| 28 |
|
|
fSumEsubEpPos(0), fSumEaddEpPos(0), fSumEsubEpNeg(0),
|
| 29 |
|
|
fSumEaddEpNeg(0), fSumHfEsubEpPos(0), fSumHfEaddEpPos(0),
|
| 30 |
|
|
fSumHfEsubEpNeg(0), fSumHfEaddEpNeg(0), fHPTrkFrac(0) {}
|
| 31 |
loizides |
1.1 |
EvtSelData(Double_t HcalNeg, Double_t HcalPos,
|
| 32 |
|
|
Double_t HfNeg, Double_t HfPos, Double_t HfNegTime, Double_t HfPosTime,
|
| 33 |
|
|
Double_t CaNeg, Double_t CaPos, Double_t CaNegTime, Double_t CaPosTime,
|
| 34 |
|
|
Double_t ZdcNeg, Double_t ZdcPos, Double_t ZdcNegTime, Double_t ZdcPosTime,
|
| 35 |
edwenger |
1.5 |
Int_t PxbHits, Int_t PxHits, Double_t ClusVtxQual, Double_t ClusVtxDiff,
|
| 36 |
loizides |
1.6 |
Int_t nHfNegHits, Int_t nHfPosHits, Int_t nHfTowersP, Int_t nHfTowersN,
|
| 37 |
|
|
Double_t sumEsubEpPos, Double_t sumEaddEpPos, Double_t sumEsubEpNeg,
|
| 38 |
|
|
Double_t sumEaddEpNeg, Double_t sumHfEsubEpPos, Double_t sumHfEaddEpPos,
|
| 39 |
|
|
Double_t sumHfEsubEpNeg, Double_t sumHfEaddEpNeg, Double_t HPTrkFrac);
|
| 40 |
loizides |
1.1 |
~EvtSelData() {}
|
| 41 |
|
|
|
| 42 |
loizides |
1.2 |
Double_t HcalNeg() const { return fHcalNeg; }
|
| 43 |
|
|
Double_t HcalPos() const { return fHcalPos; }
|
| 44 |
loizides |
1.1 |
Double_t CastorNeg() const { return fCaNeg; }
|
| 45 |
|
|
Double_t CastorPos() const { return fCaPos; }
|
| 46 |
|
|
Double_t CastorNegTime() const { return fCaNegTime; }
|
| 47 |
|
|
Double_t CastorPosTime() const { return fCaPosTime; }
|
| 48 |
loizides |
1.4 |
Double_t ClusVtxQual() const { return fClusVtxQual;}
|
| 49 |
|
|
Double_t ClusVtxDiff() const { return fClusVtxDiff;}
|
| 50 |
edwenger |
1.5 |
Double_t HPTrkFrac() const { return fHPTrkFrac; }
|
| 51 |
loizides |
1.1 |
Double_t HfNeg() const { return fHfNeg; }
|
| 52 |
|
|
Double_t HfPos() const { return fHfPos; }
|
| 53 |
|
|
Double_t HfNegTime() const { return fHfNegTime; }
|
| 54 |
|
|
Double_t HfPosTime() const { return fHfPosTime; }
|
| 55 |
|
|
Int_t NpixBarrel() const { return fPxbHits; }
|
| 56 |
edwenger |
1.3 |
Int_t Npix() const { return fPxHits; }
|
| 57 |
sixie |
1.7 |
Int_t NHfTowersP() const { return fNHfTowersP; }
|
| 58 |
|
|
Int_t NHfTowersN() const { return fNHfTowersN; }
|
| 59 |
loizides |
1.1 |
Double_t ZdcNeg() const { return fZdcNeg; }
|
| 60 |
|
|
Double_t ZdcPos() const { return fZdcPos; }
|
| 61 |
|
|
Double_t ZdcNegTime() const { return fZdcNegTime; }
|
| 62 |
|
|
Double_t ZdcPosTime() const { return fZdcPosTime; }
|
| 63 |
|
|
EObjType ObjType() const { return kEvtSelData; }
|
| 64 |
|
|
void Set (Double_t HcalNeg, Double_t HcalPos,
|
| 65 |
|
|
Double_t HfNeg, Double_t HfPos, Double_t HfNegTime, Double_t HfPosTime,
|
| 66 |
|
|
Double_t CaNeg, Double_t CaPos, Double_t CaNegTime, Double_t CaPosTime,
|
| 67 |
|
|
Double_t ZdcNeg, Double_t ZdcPos, Double_t ZdcNegTime, Double_t ZdcPosTime,
|
| 68 |
edwenger |
1.5 |
Int_t PxbHits, Int_t PxHits, Double_t ClusVtxQual, Double_t ClusVtxDiff,
|
| 69 |
loizides |
1.6 |
Int_t nHfNegHits, Int_t nHfPosHits, Int_t nHfTowersP, Int_t nHfTowersN,
|
| 70 |
|
|
Double_t sumEsubEpPos, Double_t sumEaddEpPos, Double_t sumEsubEpNeg,
|
| 71 |
|
|
Double_t sumEaddEpNeg, Double_t sumHfEsubEpPos, Double_t sumHfEaddEpPos,
|
| 72 |
|
|
Double_t sumHfEsubEpNeg, Double_t sumHfEaddEpNeg,
|
| 73 |
edwenger |
1.5 |
Double_t HPTrkFrac);
|
| 74 |
loizides |
1.1 |
|
| 75 |
|
|
protected:
|
| 76 |
loizides |
1.6 |
Double_t fHcalNeg; //energy HCAL negative side
|
| 77 |
|
|
Double_t fHcalPos; //energy HCAL positive side
|
| 78 |
|
|
Double_t fHfNeg; //energy HF negative side
|
| 79 |
|
|
Double_t fHfPos; //energy HF positive side
|
| 80 |
|
|
Double_t fHfNegTime; //energy weighted HF time on negative side
|
| 81 |
|
|
Double_t fHfPosTime; //energy weighted HF time on positive side
|
| 82 |
|
|
Double_t fCaNeg; //energy CASTOR negative side
|
| 83 |
|
|
Double_t fCaPos; //energy CASTOR positive side
|
| 84 |
|
|
Double_t fCaNegTime; //energy weighted CASTOR time on negative side
|
| 85 |
|
|
Double_t fCaPosTime; //energy weighted CASTOR time on positive side
|
| 86 |
|
|
Double_t fZdcNeg; //energy ZDC negative side
|
| 87 |
|
|
Double_t fZdcPos; //energy ZDC positive side
|
| 88 |
|
|
Double_t fZdcNegTime; //energy weighted ZDC time on negative side
|
| 89 |
|
|
Double_t fZdcPosTime; //energy weighted ZDC time on positive side
|
| 90 |
|
|
Int_t fPxbHits; //number of pixel rechits in the three barrel layers
|
| 91 |
|
|
Int_t fPxHits; //number of pixel rechits in all barrel and forward layers
|
| 92 |
|
|
Double_t fClusVtxQual; //incompatibility of pixel cluster shapes with vertex (ratio)
|
| 93 |
|
|
Double_t fClusVtxDiff; //incompatibility of pixel cluster shapes with vertex (diff)
|
| 94 |
|
|
Int_t fNHfNegHits; //hf neg hits above threshold
|
| 95 |
|
|
Int_t fNHfPosHits; //hf pos hits above threshold
|
| 96 |
|
|
Int_t fNHfTowersP; //hf neg calo towers above threshold
|
| 97 |
|
|
Int_t fNHfTowersN; //hf pos calo towers above threshold
|
| 98 |
|
|
Double_t fSumEsubEpPos; //sum E sub Ep for pos calo towers
|
| 99 |
|
|
Double_t fSumEaddEpPos; //sum E add Ep for pos calo towers
|
| 100 |
|
|
Double_t fSumEsubEpNeg; //sum E sub Ep for neg calo towers
|
| 101 |
|
|
Double_t fSumEaddEpNeg; //sum E add Ep for neg calo towers
|
| 102 |
|
|
Double_t fSumHfEsubEpPos; //sum E sub Ep for pos hf calo towers
|
| 103 |
|
|
Double_t fSumHfEaddEpPos; //sum E add Ep for pos hf calo towers
|
| 104 |
|
|
Double_t fSumHfEsubEpNeg; //sum E sub Ep for neg hf calo towers
|
| 105 |
|
|
Double_t fSumHfEaddEpNeg; //sum E add Ep for neg hf calo towers
|
| 106 |
|
|
Double_t fHPTrkFrac; //fraction of high-purity tracks out of all with "loose" cuts
|
| 107 |
loizides |
1.1 |
|
| 108 |
edwenger |
1.5 |
ClassDef(EvtSelData, 3) // Event selection data
|
| 109 |
loizides |
1.1 |
};
|
| 110 |
|
|
}
|
| 111 |
|
|
#endif
|
