ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/TreeFiller/src/FillerEmbedWeight.cc
(Generate patch)

Comparing UserCode/MitProd/TreeFiller/src/FillerEmbedWeight.cc (file contents):
Revision 1.4 by pharris, Mon May 6 18:27:49 2013 UTC vs.
Revision 1.5 by pharris, Thu Jul 11 14:00:41 2013 UTC

# Line 22 | Line 22 | FillerEmbedWeight::FillerEmbedWeight(con
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   {
# Line 52 | Line 53 | void FillerEmbedWeight::FillDataBlock(co
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;
# Line 107 | Line 108 | void FillerEmbedWeight::FillDataBlock(co
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());  
# Line 182 | Line 184 | void FillerEmbedWeight::FillDataBlock(co
184    embedWeight->SetGenTau2VsGenTau1PtRec    (inGenTau2PtVsGenTau1PtRec);
185    embedWeight->SetGenTau2VsGenTau1EtaRec   (inGenTau2EtaVsGenTau1EtaRec);
186    embedWeight->SetDiTauMassVsGenDiTauPtRec (inDiTauMassVsGenDiTauPtRec);
187 <  embedWeight->SetWeight();
187 >  if(recHit_) embedWeight->SetWeight();
188    embedWeight_->Trim();
189   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines