|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
#ifndef DataFormats_CSCRecHit2D_H
#define DataFormats_CSCRecHit2D_H
/**
* \class CSCRecHit2D
* Describes a 2-dim reconstructed hit in one layer of an Endcap Muon CSC.
*
* \author Tim Cox et al.
*
*/
#include "DataFormats/Common/interface/RangeMap.h"
#include <DataFormats/TrackingRecHit/interface/RecHit2DLocalPos.h>
#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include <vector>
#include <map>
#include <iosfwd>
class CSCRecHit2D final : public RecHit2DLocalPos {
public:
typedef std::vector<int> ChannelContainer;
typedef edm::RangeMap<int, std::vector<float> > ADCContainer;
enum SharedInputType {
all = TrackingRecHit::all,
some = TrackingRecHit::some,
allWires,
someWires,
allStrips,
someStrips
};
static const unsigned int MAXSTRIPS = 3;
static const unsigned int MAXTIMEBINS = 4;
static const unsigned int N_ADC = MAXSTRIPS * MAXTIMEBINS;
CSCRecHit2D();
CSCRecHit2D(const CSCDetId& id,
const LocalPoint& pos,
const LocalError& err,
const ChannelContainer& channels,
const ADCContainer& adcs,
const ChannelContainer& wgroups,
float tpeak,
float posInStrip,
float errInStrip,
int quality,
short int badStrip = 0,
short int badWireGroup = 0,
int scaledWireTime = 0,
float energyDeposit = -995.);
~CSCRecHit2D() override;
/// RecHit2DLocalPos base class interface
CSCRecHit2D* clone() const override { return new CSCRecHit2D(*this); }
LocalPoint localPosition() const override { return theLocalPosition; }
LocalError localPositionError() const override { return theLocalError; }
CSCDetId cscDetId() const { return geographicalId(); }
/// Extracting strip channel numbers comprising the rechit - low
int channels(unsigned int i) const { return theStrips_[i]; }
unsigned int nStrips() const { return nStrips_; }
/// Extract the L1A phase bits from the StripChannelContainer - high
int channelsl1a(unsigned int i) const { return theL1APhaseBits_[i]; } /// L1A
/// Container of wire groups comprising the rechit
short int hitWire() const { return hitWire_; }
short int wgroupsBX() const { return theWGroupsBX_; }
unsigned int nWireGroups() const { return nWireGroups_; }
/// Map of strip ADCs for strips comprising the rechit
float adcs(unsigned int strip, unsigned int timebin) const { return theADCs_[strip * MAXTIMEBINS + timebin]; }
unsigned int nTimeBins() const { return nTimeBins_; }
/// Fitted peaking time
float tpeak() const { return theTpeak; }
/// The estimated position within the strip
float positionWithinStrip() const { return thePositionWithinStrip; };
/// The uncertainty of the estimated position within the strip
float errorWithinStrip() const { return theErrorWithinStrip; };
/// quality flag of the reconstruction
int quality() const { return theQuality; }
/// flags for involvement of 'bad' channels
short int badStrip() const { return theBadStrip; }
short int badWireGroup() const { return theBadWireGroup; }
// Calculated wire time in ns
float wireTime() const { return (float)theScaledWireTime / 100.; }
/// Energy deposited in the layer. Note: this value is dE. In order to
/// get the dE/dX, you will need to divide by the path length.
/// Specific failure values...
/// If the user has chosen not to use the gas gain correction ---> -998.
/// If the gas gain correction from the database is a bad value -> -997.
/// If it is an edge strip --------------------------------------> -996.
/// If gas-gain is OK, but the ADC vector is the wrong size ----> -999.
/// If the user has created the Rechit without the energy deposit> -995.
/// If the user has created the Rechit with no arguments --------> -994.
float energyDepositedInLayer() const { return theEnergyDeposit; }
/// Returns true if the two TrackingRecHits are using the same input information, false otherwise. In this case, looks at the geographical ID and channel numbers for strips and wires.
bool sharesInput(const TrackingRecHit* other, TrackingRecHit::SharedInputType what) const override;
/// Returns true if the two TrackingRecHits are using the same input information, false otherwise. In this case, looks at the geographical ID and channel numbers for strips and wires.
bool sharesInput(const TrackingRecHit* other, CSCRecHit2D::SharedInputType what) const;
/// Returns true if the two CSCRecHits are using the same input information, false otherwise. In this case, looks at the geographical ID and channel numbers for strips and wires.
bool sharesInput(const CSCRecHit2D* otherRecHit, CSCRecHit2D::SharedInputType what) const;
/// Print the content of the RecHit2D including L1A (for debugging)
void print() const;
private:
float theTpeak;
float thePositionWithinStrip;
float theErrorWithinStrip;
float theEnergyDeposit;
int theQuality;
int theScaledWireTime;
short int hitWire_;
short int theWGroupsBX_;
short int theBadStrip;
short int theBadWireGroup;
unsigned char nStrips_, nWireGroups_, nTimeBins_;
unsigned char theL1APhaseBits_[MAXSTRIPS];
unsigned char theStrips_[MAXSTRIPS];
float theADCs_[N_ADC];
LocalPoint theLocalPosition;
LocalError theLocalError;
};
/// Output operator for CSCRecHit2D
std::ostream& operator<<(std::ostream& os, const CSCRecHit2D& rh);
#endif
|