diff --git a/Makefile b/Makefile index 44bb1c5..012c3ed 100644 --- a/Makefile +++ b/Makefile @@ -11,7 +11,7 @@ FFLAGS = -Wall -O2 CXXFLAGS += -std=c++11 LDFLAGS += -std=c++11 -FJCONFIG = /Users/mverweij/soft/fastjet330-install/bin/fastjet-config +FJCONFIG = /Users/mverweij/soft/fastjet332-install/bin/fastjet-config INCLUDE += `$(FJCONFIG) --cxxflags` LIBRARIES += `$(FJCONFIG) --libs --plugins` -lfastjetcontribfragile @@ -31,14 +31,14 @@ COMMONSRC = F77SRC = COMMONOBJ = -PROGSRC = runConversionQPYTHIA.cc runCreatePythiaEvents.cc runCreatePythiaEventsPartonLevel.cc runCreateThermalEvents.cc runCSVariations.cc runFromFile.cc runJetPerformance.cc runJetProfile.cc runJetProfileJewelSub.cc runJewelSub.cc runSDGenVarious.cc runSDGenVariousJewelSub.cc runSharedLayerSubtraction.cc runSimpleJetAnalysis.cc runtest.cc -PROGOBJ = runConversionQPYTHIA.o runCreatePythiaEvents.o runCreatePythiaEventsPartonLevel.o runCreateThermalEvents.o runCSVariations.o runFromFile.o runJetPerformance.o runJetProfile.o runJetProfileJewelSub.o runJewelSub.o runSDGenVarious.o runSDGenVariousJewelSub.o runSharedLayerSubtraction.o runSimpleJetAnalysis.o runtest.o +PROGSRC = runConversionQPYTHIA.cc runCreatePythiaEvents.cc runCreatePythiaEventsPartonLevel.cc runCreateThermalEvents.cc runCreateVinciaEvents.cc runCSVariations.cc runFromFile.cc runJetPerformance.cc runJetProfile.cc runJetProfileJewelSub.cc runJewelSub.cc runSDGenVarious.cc runSDGenVariousJewelSub.cc runSharedLayerSubtraction.cc runSimpleJetAnalysis.cc runtest.cc +PROGOBJ = runConversionQPYTHIA.o runCreatePythiaEvents.o runCreatePythiaEventsPartonLevel.o runCreateThermalEvents.o runCreateVinciaEvents.o runCSVariations.o runFromFile.o runJetPerformance.o runJetProfile.o runJetProfileJewelSub.o runJewelSub.o runSDGenVarious.o runSDGenVariousJewelSub.o runSharedLayerSubtraction.o runSimpleJetAnalysis.o runtest.o INCLUDE += LIBRARIES += -LPU14 -lPU14 -lz -all: runConversionQPYTHIA runCreatePythiaEvents runCreatePythiaEventsPartonLevel runCreateThermalEvents runCSVariations runFromFile runJetPerformance runJetProfile runJetProfileJewelSub runJewelSub runSDGenVarious runSDGenVariousJewelSub runSharedLayerSubtraction runSimpleJetAnalysis runtest +all: runConversionQPYTHIA runCreatePythiaEvents runCreatePythiaEventsPartonLevel runCreateThermalEvents runCreateVinciaEvents runCSVariations runFromFile runJetPerformance runJetProfile runJetProfileJewelSub runJewelSub runSDGenVarious runSDGenVariousJewelSub runSharedLayerSubtraction runSimpleJetAnalysis runtest runConversionQPYTHIA: runConversionQPYTHIA.o $(COMMONOBJ) @@ -53,6 +53,9 @@ runCreatePythiaEventsPartonLevel: runCreatePythiaEventsPartonLevel.o $(COMMONOB runCreateThermalEvents: runCreateThermalEvents.o $(COMMONOBJ) $(CXX) $(LDFLAGS) -o $@ $@.o $(COMMONOBJ) $(LIBRARIES) +runCreateVinciaEvents: runCreateVinciaEvents.o $(COMMONOBJ) + $(CXX) $(LDFLAGS) -o $@ $@.o $(COMMONOBJ) $(LIBRARIES) + runCSVariations: runCSVariations.o $(COMMONOBJ) $(CXX) $(LDFLAGS) -o $@ $@.o $(COMMONOBJ) $(LIBRARIES) @@ -94,7 +97,7 @@ clean: rm -vf $(COMMONOBJ) $(PROGOBJ) realclean: clean - rm -vf runConversionQPYTHIA runCreatePythiaEvents runCreatePythiaEventsPartonLevel runCreateThermalEvents runCSVariations runFromFile runJetPerformance runJetProfile runJetProfileJewelSub runJewelSub runSDGenVarious runSDGenVariousJewelSub runSharedLayerSubtraction runSimpleJetAnalysis runtest + rm -vf runConversionQPYTHIA runCreatePythiaEvents runCreatePythiaEventsPartonLevel runCreateThermalEvents runCreateVinciaEvents runCSVariations runFromFile runJetPerformance runJetProfile runJetProfileJewelSub runJewelSub runSDGenVarious runSDGenVariousJewelSub runSharedLayerSubtraction runSimpleJetAnalysis runtest .cc.o: $< $(CXX) $(CXXFLAGS) $(INCLUDE) -c $< -o $@ @@ -124,6 +127,9 @@ runCreatePythiaEventsPartonLevel.o: include/extraInfo.hh include/extraInfo.hh runCreatePythiaEventsPartonLevel.o: PU14/CmdLine.hh runCreateThermalEvents.o: include/ProgressBar.h include/thermalEvent.hh runCreateThermalEvents.o: include/extraInfo.hh PU14/CmdLine.hh +runCreateVinciaEvents.o: include/ProgressBar.h include/vinciaEvent.hh +runCreateVinciaEvents.o: include/extraInfo.hh include/extraInfo.hh +runCreateVinciaEvents.o: PU14/CmdLine.hh runCSVariations.o: include/ProgressBar.h PU14/EventMixer.hh PU14/CmdLine.hh runCSVariations.o: PU14/EventSource.hh PU14/CmdLine.hh PU14/PU14.hh runCSVariations.o: PU14/HepPID/ParticleIDMethods.hh include/jetCollection.hh diff --git a/PU14/Makefile b/PU14/Makefile index e9d3228..48d5418 100644 --- a/PU14/Makefile +++ b/PU14/Makefile @@ -11,7 +11,7 @@ FFLAGS = -Wall -O2 CXXFLAGS += -std=c++11 LDFLAGS += -std=c++11 -FJCONFIG = /Users/mverweij/soft/fastjet330-install/bin/fastjet-config +FJCONFIG = /Users/mverweij/soft/fastjet332-install/bin/fastjet-config INCLUDE += `$(FJCONFIG) --cxxflags` LIBRARIES += `$(FJCONFIG) --libs --plugins` -lfastjetcontribfragile LIBRARIES += -L/usr/local/Cellar/gsl/2.5/lib -lgsl -lgslcblas diff --git a/README_Simple.md b/README_Simple.md index 95e64ae..351ba5b 100644 --- a/README_Simple.md +++ b/README_Simple.md @@ -26,18 +26,18 @@ cd .. ### Install fastjet ```sh -curl -O http://fastjet.fr/repo/fastjet-3.3.0.tar.gz -tar zxvf fastjet-3.3.0.tar.gz -cd fastjet-3.3.0/ +curl -O http://fastjet.fr/repo/fastjet-3.3.2.tar.gz +tar zxvf fastjet-3.3.2.tar.gz +cd fastjet-3.3.2/ -./configure --prefix=$PWD/../fastjet-install +./configure --prefix=$PWD/../fastjet332-install make make check make install -FASTJET=$PWD/../fastjet-install +FASTJET=$PWD/../fastjet332-install cd .. -export FJ_CONTRIB_VER=1.026 +export FJ_CONTRIB_VER=1.041 curl -Lo source.tar.gz http://fastjet.hepforge.org/contrib/downloads/fjcontrib-"$FJ_CONTRIB_VER".tar.gz tar xzf source.tar.gz cd fjcontrib-"$FJ_CONTRIB_VER" diff --git a/include/softDropGroomer.hh b/include/softDropGroomer.hh index 9367a06..24c63c7 100644 --- a/include/softDropGroomer.hh +++ b/include/softDropGroomer.hh @@ -12,6 +12,7 @@ #include "fastjet/PseudoJet.hh" #include "fastjet/ClusterSequenceArea.hh" #include "fastjet/contrib/SoftDrop.hh" +#include "fastjet/contrib/Recluster.hh" #include "TMath.h" #include "TH2.h" @@ -60,14 +61,6 @@ private : std::vector injectedtrack_z_; std::vector injectedtrack_theta_; std::vector injectedtrack_pt_; - std::vector jetProfile0_; - std::vector jetProfile1_; - std::vector jetProfile2_; - std::vector jetProfile3_; - std::vector jetProfile4_; - std::vector jetProfile5_; - std::vector jetProfile6_; - std::vector jetProfile7_; std::vector> constituents_; std::vector> constituents1_; @@ -100,7 +93,6 @@ public : std::vector getInjectedz() const; std::vector getInjectedtheta() const; std::vector getInjectedpt() const; - std::vector getJetProfile(Int_t i) const; std::vector getRecur_N() const; std::vector getSubJetLeadingTrackPt() const; @@ -108,10 +100,6 @@ public : std::vector doGrooming(jetCollection &c); std::vector doGrooming(std::vector v); std::vector doGrooming(); - - TH2F* calculateProfile(jetCollection &c); - TH2F* calculateProfile(std::vector v); - TH2F* calculateProfile(); std::vector> getConstituents() {return constituents_;} std::vector> getConstituents1() {return constituents1_;} @@ -263,19 +251,6 @@ std::vector softDropGroomer::getInjectedpt() const return injectedtrack_pt_; } -std::vector softDropGroomer::getJetProfile(Int_t i) const -{ - if(i==0) return jetProfile0_; - if(i==1) return jetProfile1_; - if(i==2) return jetProfile2_; - if(i==3) return jetProfile3_; - if(i==4) return jetProfile4_; - if(i==5) return jetProfile5_; - if(i==6) return jetProfile6_; - if(i==7) return jetProfile7_; - -} - std::vector softDropGroomer::doGrooming(jetCollection &c) { return doGrooming(c.getJet()); @@ -287,48 +262,6 @@ std::vector softDropGroomer::doGrooming(std::vector v) -{ - setInputJets(v); - return calculateProfile(); -} - -TH2F* softDropGroomer::calculateProfile() -{ - TH2F *h2PtJetRho = new TH2F("h2PtJetRhotemp","",150,0,150,20,0,0.5); - int nJets = fjInputs_.size(); - double min_track_pt = 1; - double max_track_r = 0.5; - for(int i = 0;i tracks_all = fjInputs_[i].constituents(); - std::vector tracks, tracks_neutral; - SelectorIsCharged().sift(tracks_all, tracks, tracks_neutral); - Double_t jetPt = fjInputs_[i].perp(); - Double_t jetEta = fjInputs_[i].pseudorapidity(); - Double_t jetPhi = fjInputs_[i].phi(); - Double_t trackPt,trackEta,trackPhi,trackDelR,delEta,delPhi; - for(int j=0;j min_track_pt && trackDelR < max_track_r) { - h2PtJetRho->Fill(jetPt,trackDelR,eventWeight*(trackPt/jetPt)); - } - } - - } - return h2PtJetRho; -} - Double_t softDropGroomer::leading_track_pt(fastjet::PseudoJet jet) { std::vector< fastjet::PseudoJet > constituents = jet.constituents(); @@ -350,23 +283,18 @@ std::vector softDropGroomer::doGrooming() drBranches_.reserve(fjInputs_.size()); dr12_.reserve(fjInputs_.size()); sjmass_.reserve(fjInputs_.size()); - logdr12_.reserve(fjInputs_.size()); - logzgdr12_.reserve(fjInputs_.size()); - sjleadingtrack_.reserve(fjInputs_.size()); - jetProfile0_.reserve(fjInputs_.size()); - jetProfile1_.reserve(fjInputs_.size()); - jetProfile2_.reserve(fjInputs_.size()); - jetProfile3_.reserve(fjInputs_.size()); - jetProfile4_.reserve(fjInputs_.size()); - jetProfile5_.reserve(fjInputs_.size()); - jetProfile6_.reserve(fjInputs_.size()); - jetProfile7_.reserve(fjInputs_.size()); + //NEW STUFF that is not properly treated + // logdr12_.reserve(fjInputs_.size()); + // logzgdr12_.reserve(fjInputs_.size()); + // sjleadingtrack_.reserve(fjInputs_.size()); + //-- + for(fastjet::PseudoJet& jet : fjInputs_) { std::vector particles, ghosts; fastjet::SelectorIsPureGhost().sift(jet.constituents(), ghosts, particles); - fastjet::JetDefinition jet_def(fastjet::cambridge_algorithm,999.); + fastjet::JetDefinition jet_def(fastjet::cambridge_algorithm,fastjet::JetDefinition::max_allowable_R); fastjet::ClusterSequence cs(particles, jet_def); std::vector tempJets = fastjet::sorted_by_pt(cs.inclusive_jets()); @@ -374,64 +302,43 @@ std::vector softDropGroomer::doGrooming() fjOutputs_.push_back(fastjet::PseudoJet(0.,0.,0.,0.)); zg_.push_back(-1.); drBranches_.push_back(-1.); + dr12_.push_back(-1.); + // sjmass_.push_back(-1.); continue; } - Double_t jetProfile0=0,jetProfile1=0,jetProfile2=0,jetProfile3=0,jetProfile4=0,jetProfile5=0,jetProfile6=0,jetProfile7=0; - std::vector tracks = tempJets[0].constituents(); - Double_t jetEta = tempJets[0].pseudorapidity(); - Double_t jetPhi = tempJets[0].phi(); - Double_t trackEta,trackPhi,trackDelR,delEta,delPhi; - for(int i=0;i0.05 && trackDelR<0.1) jetProfile1+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.1 && trackDelR<0.15) jetProfile2+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.15 && trackDelR<0.2) jetProfile3+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.2 && trackDelR<0.25) jetProfile4+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.25 && trackDelR<0.3) jetProfile5+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.3 && trackDelR<0.35) jetProfile6+=tracks[i].perp()/tempJets[0].perp(); - if(trackDelR>0.35 && trackDelR<0.4) jetProfile7+=tracks[i].perp()/tempJets[0].perp(); - } - jetProfile0_.push_back(jetProfile0/0.05); - jetProfile1_.push_back(jetProfile1/0.05); - jetProfile2_.push_back(jetProfile2/0.05); - jetProfile3_.push_back(jetProfile3/0.05); - jetProfile4_.push_back(jetProfile4/0.05); - jetProfile5_.push_back(jetProfile5/0.05); - jetProfile6_.push_back(jetProfile6/0.05); - jetProfile7_.push_back(jetProfile7/0.05); - - fastjet::contrib::SoftDrop * sd = new fastjet::contrib::SoftDrop(beta_, zcut_, r0_ ); sd->set_verbose_structure(true); - fastjet::contrib::Recluster *reclusterer; - if(fReclusterAlgo == 2) reclusterer = new fastjet::contrib::Recluster(fastjet::kt_algorithm,1,true); - if(fReclusterAlgo == 1) reclusterer = new fastjet::contrib::Recluster(fastjet::antikt_algorithm,1,true); - if(fReclusterAlgo == 0) reclusterer = new fastjet::contrib::Recluster(fastjet::cambridge_algorithm,1,true); - if(fReclusterAlgo == 3) reclusterer = new fastjet::contrib::Recluster(fastjet::genkt_algorithm,1,0.5); - sd->set_reclustering(true,reclusterer); - RecursiveParents(tempJets[0],bflagAdded,fAdditionalTracks,kRecursiveAlgo); + // fastjet::JetAlgorithm jetalgo(fastjet::cambridge_algorithm); + // if(fReclusterAlgo == 1) jetalgo=fastjet::antikt_algorithm; + // else if(fReclusterAlgo == 2) jetalgo=fastjet::kt_algorithm; + // else if(fReclusterAlgo == 3) jetalgo=fastjet::genkt_algorithm; + + // fastjet::Recluster recluster(jetalgo,1,fastjet::Recluster::keep_only_hardest); + // // if(fReclusterAlgo == 3) recluster = fastjet::Recluster(jetalgo,1,0.5); + // sd->set_reclustering(true,&recluster); + + // RecursiveParents(tempJets[0],bflagAdded,fAdditionalTracks,kRecursiveAlgo); //do we have to do this here? fastjet::PseudoJet transformedJet = tempJets[0]; if ( transformedJet == 0 ) { fjOutputs_.push_back(fastjet::PseudoJet(0.,0.,0.,0.)); zg_.push_back(-1.); drBranches_.push_back(0.); + dr12_.push_back(0.); + sjmass_.push_back(0.); if(sd) { delete sd; sd = 0;} + //if(reclusterer) { delete reclusterer; sd = 0;} continue; } transformedJet = (*sd)(transformedJet); double zg = transformedJet.structure_of().symmetry(); int ndrop = transformedJet.structure_of().dropped_count(); - //std::vector dropped_symmetry = transformedJet.structure_of().dropped_symmetry(); - + zg_.push_back(zg); + drBranches_.push_back(ndrop); + //get distance between the two subjets std::vector subjets; if ( transformedJet.has_pieces() ) { @@ -439,36 +346,36 @@ std::vector softDropGroomer::doGrooming() fastjet::PseudoJet subjet1 = subjets[0]; fastjet::PseudoJet subjet2 = subjets[1]; dr12_.push_back(subjet1.delta_R(subjet2)); - if(subjet1.perp() > subjet2.perp()){ - sjmass_.push_back(subjet2.m()/subjet2.perp()); - sjleadingtrack_.push_back(leading_track_pt(subjet2)); - } - else{ - sjmass_.push_back(subjet1.m()/subjet1.perp()); - sjleadingtrack_.push_back(leading_track_pt(subjet1)); - } - - if(subjet1.delta_R(subjet2)!=0){ - logdr12_.push_back(TMath::Log(1/subjet1.delta_R(subjet2))); - logzgdr12_.push_back(TMath::Log(zg*subjet1.delta_R(subjet2))); - } - else{ - logdr12_.push_back(0.); - logzgdr12_.push_back(0.); - } + + // if(subjet1.perp() > subjet2.perp()){ + // sjmass_.push_back(subjet2.m()/subjet2.perp()); + // sjleadingtrack_.push_back(leading_track_pt(subjet2)); + // } + // else{ + // sjmass_.push_back(subjet1.m()/subjet1.perp()); + // sjleadingtrack_.push_back(leading_track_pt(subjet1)); + // } + + // if(subjet1.delta_R(subjet2)!=0){ + // logdr12_.push_back(TMath::Log(1/subjet1.delta_R(subjet2))); + // logzgdr12_.push_back(TMath::Log(zg*subjet1.delta_R(subjet2))); + // } + // else{ + // logdr12_.push_back(0.); + // logzgdr12_.push_back(0.); + // } } else { dr12_.push_back(-1.); - sjmass_.push_back(-1.); - sjleadingtrack_.push_back(-1.); - logdr12_.push_back(0.); - logzgdr12_.push_back(0.); + // sjmass_.push_back(-1.); + // sjleadingtrack_.push_back(-1.); + // logdr12_.push_back(0.); + // logzgdr12_.push_back(0.); } - fjOutputs_.push_back( transformedJet ); //put CA reclusterd jet after softDrop into vector - zg_.push_back(zg); - drBranches_.push_back(ndrop); + fjOutputs_.push_back( transformedJet ); //put CA reclustered jet after softDrop into vector if(sd) { delete sd; sd = 0;} + // if(reclusterer) { delete reclusterer; sd = 0;} } return fjOutputs_; } @@ -478,7 +385,7 @@ void softDropGroomer::RecursiveParents(fastjet::PseudoJet fJet, Bool_t bflagAdde std::vector fInputVectors; fInputVectors.clear(); - double xflagalgo=0; + //double xflagalgo=0; Float_t kpval=0; fInputVectors = fJet.constituents(); if(bflagAdded && fJet.m()!=0){ @@ -487,14 +394,14 @@ void softDropGroomer::RecursiveParents(fastjet::PseudoJet fJet, Bool_t bflagAdde //add tracks to the jet prior to the reclusterer in case of iterative mapping of splittings //MyJet.reset(fJet.four_mom()); - Double_t omegac=0.5*fqhat*fxlength*fxlength/0.2; - Double_t thetac=TMath::Sqrt(12*0.2/(fqhat*TMath::Power(fxlength,3))); + //Double_t omegac=0.5*fqhat*fxlength*fxlength/0.2; + //Double_t thetac=TMath::Sqrt(12*0.2/(fqhat*TMath::Power(fxlength,3))); Double_t xQs=TMath::Sqrt(fqhat*fxlength); for(Int_t i=0;i(0), fastjet::BestFJ30 ); + fastjet::JetDefinition fJetDef(jetalgo, 1., kpval, static_cast(0), fastjet::BestFJ30 ); fastjet::ClusterSequence fClustSeqSA(fInputVectors, fJetDef); std::vector fOutputJets; @@ -575,19 +482,9 @@ fastjet::JetDefinition fJetDef(jetalgo, 1., kpval, static_cast softDropGroomer::doGroomingWithJewelSub(jetCollection &c, std::vector particlesDummy) diff --git a/runJetProfile.cc b/runJetProfile.cc index 2f4ed51..37a45a0 100644 --- a/runJetProfile.cc +++ b/runJetProfile.cc @@ -173,9 +173,9 @@ int main (int argc, char ** argv) { jetCollectionSigSD.addVector("logztheta", sdgSig.getLogZgDR12()); jetCollectionSigSD.addVector("sjMass", sdgSig.getSubJetMass()); jetCollectionSigSD.addVector("leadingtrack_pt", sdgSig.getSubJetLeadingTrackPt()); - for(int j=0;j<8;j++){ - jetCollectionSigSD.addVector(Form("jetProfile%d",j), sdgSig.getJetProfile(j)); - } + // for(int j=0;j<8;j++){ + // jetCollectionSigSD.addVector(Form("jetProfile%d",j), sdgSig.getJetProfile(j)); + // } jetCollection jetCollectionSigSD_Recur(groomedJets_Sig); DesignRecurColl(&jetCollectionSigSD_Recur,sdgSig); @@ -212,9 +212,9 @@ int main (int argc, char ** argv) { jetCollectionCSSD.addVector(Form("log_%dztheta",ics), sdgCS.getLogZgDR12()); jetCollectionCSSD.addVector(Form("sjMass_CS%d",ics), sdgCS.getSubJetMass()); jetCollectionCSSD.addVector(Form("leadingtrack_pt_CS%d",ics), sdgCS.getSubJetLeadingTrackPt()); - for(int j=0;j<8;j++){ - jetCollectionCSSD.addVector(Form("jetProfile%d",j), sdgCS.getJetProfile(j)); - } + // for(int j=0;j<8;j++){ + // jetCollectionCSSD.addVector(Form("jetProfile%d",j), sdgCS.getJetProfile(j)); + // } jetCollection jetCollectionCSSD_Recur(groomedJets_CS); DesignRecurColl(&jetCollectionCSSD_Recur,sdgCS,ics); jetCollection jetCollectionCSSD_Recur_Toy(groomedJets_CS_Toy); diff --git a/runJetProfileJewelSub.cc b/runJetProfileJewelSub.cc index fd62b9b..a1980ec 100644 --- a/runJetProfileJewelSub.cc +++ b/runJetProfileJewelSub.cc @@ -161,8 +161,8 @@ int main (int argc, char ** argv) { sdgSig_Toy.setMediumParameters(1,3); std::vector groomedJets_Sig = sdgSig.doGroomingWithJewelSub(jetCollectionSig, particlesDummy); std::vector groomedJets_Sig_Toy = sdgSig_Toy.doGroomingWithJewelSub(jetCollectionSig, particlesDummy); - TH2F *temph2PtJetRho = sdgSig.calculateProfile(jetCollectionSig); - h2PtJetRho->Add(temph2PtJetRho,1); + // TH2F *temph2PtJetRho = sdgSig.calculateProfile(jetCollectionSig); + // h2PtJetRho->Add(temph2PtJetRho,1); jetCollection jetCollectionSigSD(groomedJets_Sig); jetCollectionSigSD.addVector("zgSigSD", sdgSig.getZgs()); jetCollectionSigSD.addVector("ndropSigSD", sdgSig.getNDroppedSubjets()); @@ -171,9 +171,9 @@ int main (int argc, char ** argv) { jetCollectionSigSD.addVector("logztheta", sdgSig.getLogZgDR12()); jetCollectionSigSD.addVector("sjMass", sdgSig.getSubJetMass()); jetCollectionSigSD.addVector("leadingtrack_pt", sdgSig.getSubJetLeadingTrackPt()); - for(int j=0;j<8;j++){ - jetCollectionSigSD.addVector(Form("jetProfile%d",j), sdgSig.getJetProfile(j)); - } + // for(int j=0;j<8;j++){ + // jetCollectionSigSD.addVector(Form("jetProfile%d",j), sdgSig.getJetProfile(j)); + // } jetCollection jetCollectionSigSD_Recur(groomedJets_Sig); jetCollectionSigSD_Recur.addVector("recur_jetpt", sdgSig.getRecur_JetPt()); @@ -220,9 +220,9 @@ int main (int argc, char ** argv) { jetCollectionCSSD.addVector(Form("log_%dztheta",ics), sdgCS.getLogZgDR12()); jetCollectionCSSD.addVector(Form("sjMass_CS%d",ics), sdgCS.getSubJetMass()); jetCollectionCSSD.addVector(Form("leadingtrack_pt_CS%d",ics), sdgCS.getSubJetLeadingTrackPt()); - for(int j=0;j<8;j++){ - jetCollectionCSSD.addVector(Form("jetProfile%d",j), sdgCS.getJetProfile(j)); - } + // for(int j=0;j<8;j++){ + // jetCollectionCSSD.addVector(Form("jetProfile%d",j), sdgCS.getJetProfile(j)); + // } jetCollection jetCollectionCSSD_Recur(groomedJets_CS); jetCollectionCSSD_Recur.addVector(Form("recur%d_jetpt",ics), sdgCS.getRecur_JetPt()); jetCollectionCSSD_Recur.addVector(Form("recur%d_logdr12",ics), sdgCS.getRecur_LogDR12());