Get the analysis.
1763 Int_t JetLoPass = 0;
1764 Int_t JetHiPass = 0;
1766 if (triggerResults.
isValid()){
1769 const unsigned int nTrig(triggerNames.
size());
1770 for (
unsigned int i=0;
i<nTrig;++
i)
1784 std::cout <<
">>> Trigger Lo = " << JetLoPass
1785 <<
" Hi = " << JetHiPass
1794 if (!vertexHandle.
isValid()) {
1795 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1796 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1799 if ( vertexHandle.
isValid() ){
1801 numPV = vertexCollection.size();
1811 if (!gtReadoutRecord.
isValid()) {
1812 LogInfo(
"JetAnalyzer") <<
"JetAnalyzer: Could not find GT readout record" << std::endl;
1813 if (
verbose_)
std::cout <<
"JetAnalyzer: Could not find GT readout record product" << std::endl;
1817 bool techTriggerResultBxF =
false;
1818 bool techTriggerResultBx0 =
false;
1820 if (!gtReadoutRecord.
isValid()) {
1821 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find GT readout record" << std::endl;
1822 if (
verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find GT readout record product" << std::endl;
1832 if(techTriggerResultBx0!=0){
1861 bool pass_Z_selection=
false;
1875 for (
unsigned int i=0;
i<Muons->size();
i++){
1876 bool pass_muon_id=
false;
1877 bool pass_muon_iso=
false;
1878 double dxy=fabs((*Muons)[
i].muonBestTrack()->
dxy());
1879 double dz=fabs((*Muons)[
i].muonBestTrack()->
dz());
1881 dxy=fabs((*Muons)[
i].muonBestTrack()->
dxy((*vertexHandle)[0].
position()));
1882 dz=fabs((*Muons)[
i].muonBestTrack()->
dz((*vertexHandle)[0].
position()));
1884 if((*Muons)[
i].
pt()>20 && fabs((*Muons)[
i].
eta())<2.3){
1885 if((*Muons)[
i].isGlobalMuon() && (*Muons)[
i].isPFMuon() && (*Muons)[
i].
globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && (*Muons)[
i].numberOfMatchedStations() > 1 && dxy < 0.2 && (*Muons)[
i].numberOfMatchedStations() > 1 && dz<0.5 && (*Muons)[
i].
innerTrack()->hitPattern().numberOfValidPixelHits() > 0 && (*Muons)[
i].
innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5){
1889 float muonIsoPFdb = ((*Muons)[
i].pfIsolationR04().sumChargedHadronPt
1891 if(muonIsoPFdb<0.12){
1895 if(pass_muon_id && pass_muon_iso){
1896 if((*Muons)[
i].
pt()>pt0){
1897 mu_index1=mu_index0;
1900 pt0=(*Muons)[
i].pt();
1901 }
else if ((*Muons)[
i].
pt()>pt1){
1903 pt1=(*Muons)[
i].pt();
1908 if(mu_index0>=0 && mu_index1>=0){
1909 if((*Muons)[mu_index0].
charge()*(*Muons)[mu_index1].
charge()<0){
1910 zCand=(*Muons)[mu_index0].polarP4()+(*Muons)[mu_index1].polarP4();
1911 if(fabs(zCand.M()-91.)<20 && zCand.Pt()>30){
1912 pass_Z_selection=
true;
1929 bool jetCollectionIsValid =
false;
1936 if (
jetCleaningFlag_ && (!jetCollectionIsValid || !bPrimaryVertex || !dcsDecision))
return;
1938 unsigned int collSize=-1;
1948 if (
isPFJet_) collSize=pfJets->size();
1953 std::vector<Jet> corJets;
1964 bool cleaned_first_jet=
false;
1965 bool cleaned_second_jet=
false;
1968 int ind1_mu_vetoed=-1;
1969 double pt1_mu_vetoed=-1;
1970 int ind2_mu_vetoed=-1;
1971 double pt2_mu_vetoed=-1;
1972 bool cleaned_first_jet_mu_vetoed=
false;
1973 bool cleaned_second_jet_mu_vetoed=
false;
1982 bool pass_correction_flag=
false;
1986 pass_correction_flag=
true;
1990 pass_correction_flag=
true;
1993 for (
unsigned int ijet=0; ijet<collSize; ijet++) {
1996 bool pass_uncorrected=
false;
1997 bool pass_corrected=
false;
1999 correctedJet=(*caloJets)[ijet];
2005 correctedJet=(*pfJets)[ijet];
2008 correctedJet=(*patJets)[ijet];
2011 pass_uncorrected=
true;
2014 pass_uncorrected=
true;
2018 scale = jetCorr->
correction((*caloJets)[ijet]);
2021 scale = jetCorr->
correction((*pfJets)[ijet]);
2026 pass_corrected=
true;
2032 recoJets.push_back(correctedJet);
2034 bool jetpassid=
true;
2035 bool Thiscleaned=
true;
2041 reco::JetID jetID = (*jetID_ValueMap_Handle)[calojetref];
2045 Thiscleaned=jetpassid;
2047 if(Thiscleaned && pass_corrected){
2055 Thiscleaned=jetpassid;
2057 if(Thiscleaned && pass_uncorrected){
2064 if(Thiscleaned && pass_corrected){
2069 if (fabs((*caloJets)[ijet].
eta()) <= 1.3) {
2072 }
else if(fabs((*caloJets)[ijet].
eta()) <3.0){
2093 if ( !(*patJets)[ijet].hasPFSpecific() )
2099 Thiscleaned = jetpassid;
2100 JetIDWPU = jetpassid;
2102 if(Thiscleaned && pass_uncorrected){
2108 if(Thiscleaned && pass_corrected){
2111 if(fabs(correctedJet.
eta()) <= 1.3) {
2112 if(correctedJet.
pt()<=50.){
2115 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2118 if(correctedJet.
pt()>140.){
2125 }
else if(fabs(correctedJet.
eta()) <= 3) {
2126 if(correctedJet.
pt()<=50.){
2129 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2132 if(correctedJet.
pt()>140.){
2139 }
else if(fabs(correctedJet.
eta()) <= 5) {
2140 if(correctedJet.
pt()<=50.){
2143 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2146 if(correctedJet.
pt()>140.){
2159 int puidmvaflag=-10;
2160 int puidcutflag=-10;
2161 puidmva=(*puJetIdMva)[pfjetref];
2162 puidcut=(*puJetId)[pfjetref];
2163 puidmvaflag=(*puJetIdFlagMva)[pfjetref];
2164 puidcutflag=(*puJetIdFlag)[pfjetref];
2167 Thiscleaned = jetpassid;
2170 if(Thiscleaned && pass_uncorrected){
2176 if(Thiscleaned && pass_corrected){
2219 if (correctedJet.
pt()<= 50) {
2224 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2229 if (correctedJet.
pt()>140.) {
2234 if (fabs(correctedJet.
eta()) <= 1.3) {
2236 if (correctedJet.
pt()<=50.) {
2256 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2276 if (correctedJet.
pt()>140.) {
2299 }
else if(fabs(correctedJet.
eta()) <= 3) {
2301 if (correctedJet.
pt()<=50.) {
2321 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2341 if (correctedJet.
pt()>140.) {
2368 if (correctedJet.
pt()<=50.) {
2382 if (correctedJet.
pt()>50. && correctedJet.
pt()<=140.) {
2396 if (correctedJet.
pt()>140.) {
2440 if(techTriggerResultBx0 && techTriggerResultBxF){
2442 if(fabs(correctedJet.
eta()) <= 1.3) {
2447 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2452 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2457 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2461 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2494 if(techTriggerResultBx0 && !techTriggerResultBxF){
2496 if(fabs(correctedJet.
eta()) <= 1.3) {
2501 }
else if (correctedJet.
eta() > -3.0 && correctedJet.
eta() <= -1.3) {
2506 }
else if (correctedJet.
eta() >= 1.3 && correctedJet.
eta() < 3.0) {
2511 }
else if (correctedJet.
eta() > -5.0 && correctedJet.
eta() <= -3.0) {
2515 }
else if (correctedJet.
eta() >= 3.0 && correctedJet.
eta() < 5.0) {
2544 if(pass_Z_selection){
2545 if(
deltaR((*Muons)[mu_index0].
eta(),(*Muons)[mu_index0].
phi(),correctedJet.
eta(),correctedJet.
phi())>0.2 &&
deltaR((*Muons)[mu_index1].eta(),(*Muons)[mu_index1].phi(),correctedJet.
eta(),correctedJet.
phi())>0.2 ){
2546 if(correctedJet.
pt()>pt1_mu_vetoed){
2547 pt2_mu_vetoed=pt1_mu_vetoed;
2548 ind2_mu_vetoed=ind1_mu_vetoed;
2549 cleaned_second_jet_mu_vetoed=cleaned_first_jet_mu_vetoed;
2550 pt1_mu_vetoed=correctedJet.
pt();
2551 ind1_mu_vetoed=ijet;
2552 cleaned_first_jet_mu_vetoed=JetIDWPU;
2553 }
else if(correctedJet.
pt()>pt2_mu_vetoed){
2554 pt2_mu_vetoed=correctedJet.
pt();
2555 ind2_mu_vetoed=ijet;
2556 cleaned_second_jet_mu_vetoed=JetIDWPU;
2561 if(correctedJet.
pt()>
pt1){
2567 cleaned_second_jet=cleaned_first_jet;
2568 pt1=correctedJet.
pt();
2570 cleaned_first_jet=JetIDWPU;
2571 }
else if(correctedJet.
pt()>
pt2){
2575 pt2=correctedJet.
pt();
2577 cleaned_second_jet=JetIDWPU;
2578 }
else if(correctedJet.
pt()>pt3){
2579 pt3=correctedJet.
pt();
2583 if(!pass_corrected){
2592 if(fabs(correctedJet.
eta())<3.0){
2598 if(fabs(correctedJet.
eta())<3.0){
2604 if(Thiscleaned && pass_corrected){
2631 if((*patJets)[ijet].hasSubjets(
"SoftDropPuppi")){
2634 if((*patJets)[ijet].hasSubjets(
"SoftDropPuppi") && !(*patJets)[ijet].subjets(
"SoftDropPuppi").empty()){
2639 if((*patJets)[ijet].subjets(
"SoftDropPuppi").size()>1){
2652 if((*patJets)[ijet].hasSubjets(
"SoftDropPuppi")){
2655 if((*patJets)[ijet].hasSubjets(
"SoftDropPuppi") && !(*patJets)[ijet].subjets(
"SoftDropPuppi").empty()){
2660 if((*patJets)[ijet].subjets(
"SoftDropPuppi").size()>1){
2693 if ((*patJets)[ijet].hasUserFloat(
"patPuppiJetSpecificProducer:puppiMultiplicity"))
2694 nConstituents = (*patJets)[ijet].userFloat(
"patPuppiJetSpecificProducer:puppiMultiplicity");
2706 if (fabs(correctedJet.
eta()) <= 1.3) {
2713 }
else if (fabs(correctedJet.
eta()) <= 3) {
2741 if(cleaned_first_jet){
2746 if(numofjets>1 && cleaned_second_jet) {
2747 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2749 dphi=2*acos(-1.)-dphi;
2762 if(
jetCleaningFlag_ && recoJets.size()>1 && cleaned_first_jet && cleaned_second_jet ){
2769 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2771 dphi=2*acos(-1.)-dphi;
2779 reco::JetID jetID1 = (*jetID_ValueMap_Handle)[calojetref1];
2781 reco::JetID jetID2 = (*jetID_ValueMap_Handle)[calojetref2];