ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc
(Generate patch)

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.cc (file contents):
Revision 1.93 by wulsin, Wed Jul 3 07:42:46 2013 UTC vs.
Revision 1.94 by jbrinson, Wed Jul 3 09:30:38 2013 UTC

# Line 109 | Line 109 | OSUAnalysis::OSUAnalysis (const edm::Par
109      if(tempInputCollection == "jet-photon pairs") tempInputCollection = "photon-jet pairs";
110      if(tempInputCollection == "jet-muon pairs") tempInputCollection = "muon-jet pairs";
111      if(tempInputCollection == "jet-met pairs")  tempInputCollection = "met-jet pairs";
112 +    if(tempInputCollection == "track-jet pairs")  tempInputCollection = "track-jet pairs";
113      if(tempInputCollection == "event-track pairs")   tempInputCollection = "track-event pairs";
114      if(tempInputCollection == "secondary muon-muon pairs")   tempInputCollection = "muon-secondary muon pairs";
115      if(tempInputCollection == "secondary jet-muon pairs")   tempInputCollection = "muon-secondary jet pairs";
# Line 292 | Line 293 | OSUAnalysis::OSUAnalysis (const edm::Par
293        if(tempInputCollection == "jet-electron pairs") tempInputCollection = "electron-jet pairs";
294        if(tempInputCollection == "jet-muon pairs") tempInputCollection = "muon-jet pairs";
295        if(tempInputCollection == "jet-met pairs")  tempInputCollection = "met-jet pairs";
296 +      if(tempInputCollection == "track-jet pairs")  tempInputCollection = "track-jet pairs";
297        if(tempInputCollection == "jet-photon pairs") tempInputCollection = "photon-jet pairs";
298        if(tempInputCollection == "event-track pairs")   tempInputCollection = "track-event pairs";
299        if(tempInputCollection == "secondary muon-muon pairs")   tempInputCollection = "muon-secondary muon pairs";
# Line 576 | Line 578 | OSUAnalysis::OSUAnalysis (const edm::Par
578          else if(currentObject == "muon-jet pairs")            currentObject = "muonJetPairs";
579          else if(currentObject == "photon-jet pairs")          currentObject = "photonJetPairs";
580          else if(currentObject == "met-jet pairs")             currentObject = "metJetPairs";
581 +        else if(currentObject == "track-jet pairs")           currentObject = "trackJetPairs";
582          else if(currentObject == "muon-secondary jet pairs")  currentObject = "muonSecondaryJetPairs";
583          else if(currentObject == "photon-secondary jet pairs")  currentObject = "photonSecondaryJetPairs";
584          else if(currentObject == "electron-secondary jet pairs")  currentObject = "electronSecondaryJetPairs";
# Line 833 | Line 836 | OSUAnalysis::produce (edm::Event &event,
836                                                                    cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
837                                                                    cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
838                                                                    "muon-jet pairs");
839 <        else if(currentObject == "met-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),jets.product(), \
839 >        else if(currentObject == "met-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,mets.product(),jets.product(), \
840                                                                    cumulativeFlags.at("mets").at(flagsForPairCutsIndex), \
841                                                                    cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
842                                                                    "met-jet pairs");
843 +        else if(currentObject == "track-jet pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,tracks.product(),jets.product(), \
844 +                                                                 cumulativeFlags.at("tracks").at(flagsForPairCutsIndex), \
845 +                                                                 cumulativeFlags.at("jets").at(flagsForPairCutsIndex), \
846 +                                                                 "track-jet pairs");
847          else if(currentObject == "muon-photon pairs") setObjectFlags(currentCut,currentCutIndex,individualFlags,cumulativeFlags,muons.product(),photons.product(), \
848                                                                    cumulativeFlags.at("muons").at(flagsForPairCutsIndex), \
849                                                                    cumulativeFlags.at("photons").at(flagsForPairCutsIndex), \
# Line 1032 | Line 1039 | OSUAnalysis::produce (edm::Event &event,
1039              else if(currentHistogram.inputCollection == "muon-jet pairs") fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
1040                                                                                            cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
1041                                                                                            cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
1042 <            else if(currentHistogram.inputCollection == "met-jet pairs")  fill1DHistogram(histo,currentHistogram, muons.product(),jets.product(), \
1042 >            else if(currentHistogram.inputCollection == "met-jet pairs")  fill1DHistogram(histo,currentHistogram, mets.product(),jets.product(), \
1043                                                                                            cumulativeFlags.at("mets").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
1044                                                                                            cumulativeFlags.at("met-jet pairs").at(currentDir),scaleFactor);
1045 +            else if(currentHistogram.inputCollection == "track-jet pairs")  fill1DHistogram(histo,currentHistogram,tracks.product(),jets.product(), \
1046 +                                                                                          cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("jets").at(currentDir),
1047 +                                                                                          cumulativeFlags.at("track-jet pairs").at(currentDir),scaleFactor);
1048              else if(currentHistogram.inputCollection == "muon-photon pairs") fill1DHistogram(histo,currentHistogram, muons.product(),photons.product(), \
1049                                                                                            cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("photons").at(currentDir),
1050                                                                                            cumulativeFlags.at("muon-photon pairs").at(currentDir),scaleFactor);
# Line 1121 | Line 1131 | OSUAnalysis::produce (edm::Event &event,
1131              else if(currentHistogram.inputCollection == "muon-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
1132                                                                                            cumulativeFlags.at("muons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
1133                                                                                            cumulativeFlags.at("muon-jet pairs").at(currentDir),scaleFactor);
1134 <            else if(currentHistogram.inputCollection == "met-jet pairs") fill2DHistogram(histo,currentHistogram,muons.product(),jets.product(), \
1134 >            else if(currentHistogram.inputCollection == "met-jet pairs") fill2DHistogram(histo,currentHistogram,mets.product(),jets.product(), \
1135                                                                                           cumulativeFlags.at("mets").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
1136                                                                                           cumulativeFlags.at("met-jet pairs").at(currentDir),scaleFactor);
1137 +            else if(currentHistogram.inputCollection == "track-jet pairs") fill2DHistogram(histo,currentHistogram,tracks.product(),jets.product(),
1138 +                                                                                         cumulativeFlags.at("tracks").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
1139 +                                                                                         cumulativeFlags.at("track-jet pairs").at(currentDir),scaleFactor);
1140              else if(currentHistogram.inputCollection == "photon-jet pairs") fill2DHistogram(histo,currentHistogram,photons.product(),jets.product(), \
1141                                                                                            cumulativeFlags.at("photons").at(currentDir),cumulativeFlags.at("jets").at(currentDir), \
1142                                                                                            cumulativeFlags.at("photon-jet pairs").at(currentDir),scaleFactor);
# Line 1185 | Line 1198 | OSUAnalysis::produce (edm::Event &event,
1198            else if(currentObject == "electron-jet pairs")                 objectToPlot = "electronJetPairs";
1199            else if(currentObject == "muon-jet pairs")                     objectToPlot = "muonJetPairs";
1200            else if(currentObject == "met-jet pairs")                      objectToPlot = "metJetPairs";
1201 +          else if(currentObject == "track-jet pairs")                    objectToPlot = "trackJetPairs";
1202            else if(currentObject == "muon-photon pairs")                  objectToPlot = "muonPhotonPairs";
1203            else if(currentObject == "photon-jet pairs")                   objectToPlot = "photonJetPairs";
1204            else if(currentObject == "jet-jet pairs")                      objectToPlot = "dijetPairs";
# Line 3484 | Line 3498 | OSUAnalysis::valueLookup (const BNphoton
3498    return value;
3499   }
3500  
3501 + // track-jet pair valueLookup
3502 + double
3503 + OSUAnalysis::valueLookup (const BNtrack* object1, const BNjet* object2, string variable, string function, string &stringValue){
3504 +
3505 +  double value = 0.0;
3506 +
3507 +  if(variable == "deltaPhi") value = fabs(deltaPhi(object1->phi,object2->phi));
3508 +  else if(variable == "deltaR") value = deltaR(object1->eta,object1->phi,object2->eta,object2->phi);
3509 +
3510 +  else{cout << "WARNING: invalid variable '" << variable << "'\n"; value = -999;}
3511 +  value = applyFunction(function, value);
3512 +
3513 +  return value;
3514 +
3515 + }
3516 +
3517 +
3518 +
3519   // met-jet pair valueLookup
3520   double
3521   OSUAnalysis::valueLookup (const BNmet* object1, const BNjet* object2, string variable, string function, string &stringValue){

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines