22 |
|
|
23 |
|
genInfo_(Conf().getUntrackedParameter<bool> ("useGenInfo","True")), |
24 |
|
recHit_ (Conf().getUntrackedParameter<bool> ("useRecHit","True")), |
25 |
+ |
useMuRad_ (Conf().getUntrackedParameter<bool> ("useMuonRad","True")), |
26 |
|
mitName_(Conf().getUntrackedParameter<string>("mitName",Names::gkEmbedWeightBrn)), |
27 |
|
embedWeight_(new mithep::EmbedWeightArr) |
28 |
|
{ |
53 |
|
// Fill missing energy from edm collection into our collection. |
54 |
|
embedWeight_->Delete(); |
55 |
|
|
56 |
< |
double inEmbedWeightValue = 0; |
57 |
< |
double inSpinnerWeightValue = 0; |
58 |
< |
double inSpinnerFlipWeightValue = 0; |
59 |
< |
double inSpinnerPlusWeightValue = 0; |
60 |
< |
double inSpinnerMinusWeightValue = 0; |
61 |
< |
double inMuEffWeightValue = 0; |
62 |
< |
double inMuEffWeightUpValue = 0; |
63 |
< |
double inMuEffWeightDownValue = 0; |
64 |
< |
double inMuRadWeightValue = 0; |
65 |
< |
double inMuRadWeightUpValue = 0; |
66 |
< |
double inMuRadWeightDownValue = 0; |
67 |
< |
double inGenTau2PtVsGenTau1PtGenValue = 0; |
68 |
< |
double inGenTau2EtaVsGenTau1EtaGenValue = 0; |
69 |
< |
double inDiTauMassVsGenDiTauPtGenValue = 0; |
70 |
< |
double inGenTau2PtVsGenTau1PtRecValue = 0; |
71 |
< |
double inGenTau2EtaVsGenTau1EtaRecValue= 0; |
72 |
< |
double inDiTauMassVsGenDiTauPtRecValue = 0; |
56 |
> |
double inEmbedWeightValue = 1; |
57 |
> |
double inSpinnerWeightValue = 1; |
58 |
> |
double inSpinnerFlipWeightValue = 1; |
59 |
> |
double inSpinnerPlusWeightValue = 1; |
60 |
> |
double inSpinnerMinusWeightValue = 1; |
61 |
> |
double inMuEffWeightValue = 1; |
62 |
> |
double inMuEffWeightUpValue = 1; |
63 |
> |
double inMuEffWeightDownValue = 1; |
64 |
> |
double inMuRadWeightValue = 1; |
65 |
> |
double inMuRadWeightUpValue = 1; |
66 |
> |
double inMuRadWeightDownValue = 1; |
67 |
> |
double inGenTau2PtVsGenTau1PtGenValue = 1; |
68 |
> |
double inGenTau2EtaVsGenTau1EtaGenValue = 1; |
69 |
> |
double inDiTauMassVsGenDiTauPtGenValue = 1; |
70 |
> |
double inGenTau2PtVsGenTau1PtRecValue = 1; |
71 |
> |
double inGenTau2EtaVsGenTau1EtaRecValue= 1; |
72 |
> |
double inDiTauMassVsGenDiTauPtRecValue = 1; |
73 |
|
|
74 |
|
if(!genInfo_) { |
75 |
|
Handle<double> hEmbedWeight; |
108 |
|
Handle<double> hMuEffWeightDown; |
109 |
|
event.getByLabel(edm::InputTag(edmMuEffName_,"weightDown","EmbeddedRECO"),hMuEffWeightDown); |
110 |
|
inMuEffWeightDownValue = *(hMuEffWeightDown.product()); |
111 |
< |
|
112 |
< |
Handle<double> hMuRadWeight; |
113 |
< |
event.getByLabel(edm::InputTag(edmMuRadName_,"weight","EmbeddedRECO"),hMuRadWeight); |
114 |
< |
inMuRadWeightValue = *(hMuRadWeight.product()); |
115 |
< |
|
116 |
< |
Handle<double> hMuRadWeightUp; |
117 |
< |
event.getByLabel(edm::InputTag(edmMuRadName_,"weightUp","EmbeddedRECO"),hMuRadWeightUp); |
118 |
< |
inMuRadWeightUpValue = *(hMuRadWeightUp.product()); |
119 |
< |
|
120 |
< |
Handle<double> hMuRadWeightDown; |
121 |
< |
event.getByLabel(edm::InputTag(edmMuRadName_,"weightDown","EmbeddedRECO"),hMuRadWeightDown); |
122 |
< |
inMuRadWeightDownValue = *(hMuRadWeightDown.product()); |
123 |
< |
|
111 |
> |
|
112 |
> |
if(useMuRad_) { |
113 |
> |
Handle<double> hMuRadWeight; |
114 |
> |
event.getByLabel(edm::InputTag(edmMuRadName_,"weight","EmbeddedRECO"),hMuRadWeight); |
115 |
> |
inMuRadWeightValue = *(hMuRadWeight.product()); |
116 |
> |
|
117 |
> |
Handle<double> hMuRadWeightUp; |
118 |
> |
event.getByLabel(edm::InputTag(edmMuRadName_,"weightUp","EmbeddedRECO"),hMuRadWeightUp); |
119 |
> |
inMuRadWeightUpValue = *(hMuRadWeightUp.product()); |
120 |
> |
|
121 |
> |
Handle<double> hMuRadWeightDown; |
122 |
> |
event.getByLabel(edm::InputTag(edmMuRadName_,"weightDown","EmbeddedRECO"),hMuRadWeightDown); |
123 |
> |
inMuRadWeightDownValue = *(hMuRadWeightDown.product()); |
124 |
> |
} |
125 |
|
Handle<double> hGenTau2PtVsGenTau1PtGen; |
126 |
|
event.getByLabel(edm::InputTag(edmKineReweightGenName_,"genTau2PtVsGenTau1Pt"),hGenTau2PtVsGenTau1PtGen); |
127 |
|
inGenTau2PtVsGenTau1PtGenValue = *(hGenTau2PtVsGenTau1PtGen.product()); |
184 |
|
embedWeight->SetGenTau2VsGenTau1PtRec (inGenTau2PtVsGenTau1PtRec); |
185 |
|
embedWeight->SetGenTau2VsGenTau1EtaRec (inGenTau2EtaVsGenTau1EtaRec); |
186 |
|
embedWeight->SetDiTauMassVsGenDiTauPtRec (inDiTauMassVsGenDiTauPtRec); |
187 |
< |
embedWeight->SetWeight(); |
187 |
> |
if(recHit_) embedWeight->SetWeight(); |
188 |
|
embedWeight_->Trim(); |
189 |
|
} |