PatBTagAnalyzer Class Reference

Inheritance diagram for PatBTagAnalyzer:

Public Member Functions
	PatBTagAnalyzer (const edm::ParameterSet &)
	~PatBTagAnalyzer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob () override

Private Attributes
std::string	bname [10]
double	BTagdisccut_
double	BTagdiscmax_
std::string	BTagdiscriminator_
bJetSelector	BTagger
PatBTagCommonHistos	BTagHistograms
std::string	BTagmethod_
BTagPerformance	BTagPerf [10]
std::string	BTagpurity_
bool	BTagverbose
std::string	cname [10]
std::string	discname [10]
std::map< std::string, TGraph * >	graphcontainer_
std::map< std::string, TGraphErrors * >	grapherrorscontainer_
std::map< std::string, TH2D * >	h2_
std::map< std::string, TH1D * >	histocontainer_
edm::EDGetTokenT< edm::View< pat::Jet > >	jetToken_
edm::ParameterSet	PatBjet_
std::map< int, std::string >	udsgname

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 22 of file PatBTagAnalyzer.cc.

Constructor & Destructor Documentation

PatBTagAnalyzer::PatBTagAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 71 of file PatBTagAnalyzer.cc.

                                                              :
  jetToken_(consumes<edm::View<pat::Jet> >(iConfig.getUntrackedParameter<edm::InputTag>("jetTag"))),
  PatBjet_(iConfig.getParameter< edm::ParameterSet >("BjetTag")),
  BTagpurity_(PatBjet_.getParameter<std::string>("purity")),
  BTagmethod_(PatBjet_.getUntrackedParameter<std::string>("tagger","TC2")),
  BTagdiscriminator_(PatBjet_.getParameter<std::string>("discriminator")),
  BTagverbose(PatBjet_.getUntrackedParameter<bool>("verbose",false)),
  BTagdisccut_(PatBjet_.getUntrackedParameter<double>("mindiscriminatorcut",5.0)),
  BTagdiscmax_(PatBjet_.getUntrackedParameter<double>("maxdiscriminatorcut",15.0)),
  BTagger(iConfig.getParameter< edm::ParameterSet >("BJetOperatingPoints")),
  BTagHistograms(iConfig)
{
   //now do what ever initialization is needed
}

PatBTagAnalyzer::~PatBTagAnalyzer ( )

override

Definition at line 87 of file PatBTagAnalyzer.cc.

{

   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void PatBTagAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 97 of file PatBTagAnalyzer.cc.

References funct::abs(), BTagPerformance::Add(), edm::View< T >::begin(), BTagdiscriminator_, BTagHistograms, BTagPerf, edm::View< T >::end(), PatBTagCommonHistos::Fill(), edm::Event::getByToken(), fwrapper::jets, and jetToken_.

{
   using namespace edm;

   // first: get all objects from the event.

   edm::Handle<edm::View<pat::Jet> > jetHandle;
   iEvent.getByToken(jetToken_,jetHandle);
   edm::View<pat::Jet> jets = *jetHandle; // get JETS

   // LOOP over all jets

   for(edm::View<pat::Jet>::const_iterator jet_iter = jets.begin(); jet_iter!=jets.end(); ++jet_iter){

       float bdiscriminator = jet_iter->bDiscriminator(BTagdiscriminator_);
       int flavor           = jet_iter->partonFlavour();
       //
       // Fill in for performance standard pt(uncorrected) >10 and abs(eta)<2.4
       if( jet_iter->correctedJet("raw").pt() > 10  &&
           fabs(jet_iter->eta()) < 2.4 ) {

           BTagPerf[0].Add(bdiscriminator,abs(flavor));

       }

       //Fill histograms
       BTagHistograms.Fill(jet_iter,"all");
       if (flavor ==  0  ) BTagHistograms.Fill(jet_iter,"no_flavor");
       if (flavor ==  5 || flavor ==  -5 ) BTagHistograms.Fill(jet_iter,"b");
       if (flavor ==  4 || flavor ==  -4 ) BTagHistograms.Fill(jet_iter,"c");
       if ((-4 < flavor && flavor < 4 && flavor != 0 )||(flavor == 21 || flavor == -21 ))
           BTagHistograms.Fill(jet_iter,"udsg");


   }//end loop over jets

}

void PatBTagAnalyzer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 136 of file PatBTagAnalyzer.cc.

References bname, BTagdisccut_, BTagdiscmax_, BTagHistograms, BTagmethod_, BTagPerf, cname, discname, BTagPerformance::GetN(), graphcontainer_, grapherrorscontainer_, histocontainer_, mps_fire::i, TFileService::make(), BTagPerformance::Set(), PatBTagCommonHistos::Set(), BTagPerformance::SetMaxDiscriminator(), BTagPerformance::SetMinDiscriminator(), PatBTagCommonHistos::Sumw2(), and udsgname.

{
  //
  // define some histograms using the framework tfileservice. Define the output file name in your .cfg.
  //
  edm::Service<TFileService> fs;

  TString suffix1="_test";

  //set performance variables collector
  for (int i=0; i < 10; i++){
    BTagPerf[i].Set(BTagmethod_);
    BTagPerf[i].SetMinDiscriminator(BTagdisccut_);
    BTagPerf[i].SetMaxDiscriminator(BTagdiscmax_);
  }

  histocontainer_["njets"]=fs->make<TH1D>("njets","jet multiplicity for jets with p_{T} > 50 GeV/c",10,0,10);
// Std. 30 pt uncorr cut for performance
  discname[0]="disc"+BTagmethod_+"_udsg";
  bname[0]   ="g"+BTagmethod_+"_b";
  cname[0]   ="g"+BTagmethod_+"_c";
  udsgname[0]="g"+BTagmethod_+"_udsg";

// 10 pt uncorr for performance + all,>0,>1,>2 tracks
  discname[1]="Uncor10_disc"+BTagmethod_+"_udsg";
  bname[1]   ="Uncor10_g"+BTagmethod_+"_b";
  cname[1]   ="Uncor10_g"+BTagmethod_+"_c";
  udsgname[1]="Uncor10_g"+BTagmethod_+"_udsg";
  discname[2]="Uncor10t0_disc"+BTagmethod_+"_udsg";
  bname[2]   ="Uncor10t0_g"+BTagmethod_+"_b";
  cname[2]   ="Uncor10t0_g"+BTagmethod_+"_c";
  udsgname[2]="Uncor10t0_g"+BTagmethod_+"_udsg";
  discname[3]="Uncor10t1_disc"+BTagmethod_+"_udsg";
  bname[3]   ="Uncor10t1_g"+BTagmethod_+"_b";
  cname[3]   ="Uncor10t1_g"+BTagmethod_+"_c";
  udsgname[3]="Uncor10t1_g"+BTagmethod_+"_udsg";
  discname[4]="Uncor10t2_disc"+BTagmethod_+"_udsg";
  bname[4]   ="Uncor10t2_g"+BTagmethod_+"_b";
  cname[4]   ="Uncor10t2_g"+BTagmethod_+"_c";
  udsgname[4]="Uncor10t2_g"+BTagmethod_+"_udsg";

// 30 pt corr for performance + all,>0,>1,>2 tracks
  discname[5]="Corr30_disc"+BTagmethod_+"_udsg";
  bname[5]   ="Corr30_g"+BTagmethod_+"_b";
  cname[5]   ="Corr30_g"+BTagmethod_+"_c";
  udsgname[5]="Corr30_g"+BTagmethod_+"_udsg";
  discname[6]="Corr30t0_disc"+BTagmethod_+"_udsg";
  bname[6]   ="Corr30t0_g"+BTagmethod_+"_b";
  cname[6]   ="Corr30t0_g"+BTagmethod_+"_c";
  udsgname[6]="Corr30t0_g"+BTagmethod_+"_udsg";
  discname[7]="Corr30t1_disc"+BTagmethod_+"_udsg";
  bname[7]   ="Corr30t1_g"+BTagmethod_+"_b";
  cname[7]   ="Corr30t1_g"+BTagmethod_+"_c";
  udsgname[7]="Corr30t1_g"+BTagmethod_+"_udsg";
  discname[8]="Corr30t2_disc"+BTagmethod_+"_udsg";
  bname[8]   ="Corr30t2_g"+BTagmethod_+"_b";
  cname[8]   ="Corr30t2_g"+BTagmethod_+"_c";
  udsgname[8]="Corr30t2_g"+BTagmethod_+"_udsg";

// check filter
  discname[9]="check_disc"+BTagmethod_+"_udsg";
  bname[9]   ="check_g"+BTagmethod_+"_b";
  cname[9]   ="check_g"+BTagmethod_+"_c";
  udsgname[9]="check_g"+BTagmethod_+"_udsg";

  for(int i=1; i<10;i++){
    graphcontainer_[discname[i]]      =fs->make<TGraph>(BTagPerf[i].GetN());       graphcontainer_[discname[i]]->SetName(discname[i].c_str());
    grapherrorscontainer_[bname[i]]   =fs->make<TGraphErrors>(BTagPerf[i].GetN()); grapherrorscontainer_[bname[i]]   ->SetName(bname[i].c_str());
    grapherrorscontainer_[cname[i]]   =fs->make<TGraphErrors>(BTagPerf[i].GetN()); grapherrorscontainer_[cname[i]]   ->SetName(cname[i].c_str());
    grapherrorscontainer_[udsgname[i]]=fs->make<TGraphErrors>(BTagPerf[i].GetN()); grapherrorscontainer_[udsgname[i]]->SetName(udsgname[i].c_str());
  }
  //Define histograms
  BTagHistograms.Set("all");
  BTagHistograms.Set("no_flavor");
  BTagHistograms.Set("b");
  BTagHistograms.Set("c");
  BTagHistograms.Set("udsg");

  // Set to save histogram errors
  BTagHistograms.Sumw2();

}

void PatBTagAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 221 of file PatBTagAnalyzer.cc.

References bname, BTagmethod_, BTagPerf, cname, DEFINE_FWK_MODULE, discname, BTagPerformance::Eval(), BTagPerformance::GetArray(), BTagPerformance::GetN(), graphcontainer_, grapherrorscontainer_, mps_fire::i, TFileDirectory::make(), TFileService::mkdir(), gen::n, and udsgname.

                        {
//ed++
   edm::Service<TFileService> fs;

// Save performance plots as Tgraphs


   for (int i=1;i<10;i++){
      BTagPerf[i].Eval();
      for (int n=0; n<BTagPerf[i].GetN();  n++ ){
         graphcontainer_[discname[i]]       ->SetPoint(n,BTagPerf[i].GetArray("udsg")[n],BTagPerf[i].GetArray("discriminator")[n]);
         grapherrorscontainer_[bname[i]]    ->SetPoint(n,BTagPerf[i].GetArray("b")[n],BTagPerf[i].GetArray("b")[n]);
         grapherrorscontainer_[cname[i]]    ->SetPoint(n,BTagPerf[i].GetArray("b")[n],BTagPerf[i].GetArray("c")[n]);
         grapherrorscontainer_[udsgname[i]] ->SetPoint(n,BTagPerf[i].GetArray("b")[n],BTagPerf[i].GetArray("udsg")[n]);
         grapherrorscontainer_[bname[i]]    ->SetPointError(n,BTagPerf[i].GetArray("bErr")[n],BTagPerf[i].GetArray("bErr")[n]);
         grapherrorscontainer_[cname[i]]    ->SetPointError(n,BTagPerf[i].GetArray("bErr")[n],BTagPerf[i].GetArray("cErr")[n]);
         grapherrorscontainer_[udsgname[i]] ->SetPointError(n,BTagPerf[i].GetArray("bErr")[n],BTagPerf[i].GetArray("udsgErr")[n]);
      }//end for over BTagPerf[i] elements
      graphcontainer_[discname[i]]     ->SetTitle("Jet udsg-mistagging");
      grapherrorscontainer_[bname[i]]   ->SetTitle("Jet b-efficiency");
      grapherrorscontainer_[cname[i]]   ->SetTitle("Jet c-mistagging");
      grapherrorscontainer_[udsgname[i]]->SetTitle("discriminator vs udsg-mistagging");
   }//end for over [i]


// Save default cut performance plot
   BTagPerf[0].Eval();

//  TFileDirectory TaggerDir = fs->mkdir(BTagmethod_);
//  TGraphErrors *BTagger_b    = new TGraphErrors(BTagTool.GetN(),
  TGraphErrors *BTagger_b    = fs->mkdir(BTagmethod_).make<TGraphErrors>(BTagPerf[0].GetN(),
                                            BTagPerf[0].GetArray("b").GetArray(),BTagPerf[0].GetArray("b").GetArray(),
                                            BTagPerf[0].GetArray("bErr").GetArray(),BTagPerf[0].GetArray("bErr").GetArray());

  TGraphErrors *BTagger_c    = new TGraphErrors(BTagPerf[0].GetN(),
                                            BTagPerf[0].GetArray("b").GetArray(),BTagPerf[0].GetArray("c").GetArray(),
                                            BTagPerf[0].GetArray("bErr").GetArray(),BTagPerf[0].GetArray("cErr").GetArray());

  TGraphErrors *BTagger_udsg = new TGraphErrors(BTagPerf[0].GetN(),
                                               BTagPerf[0].GetArray("b").GetArray(),BTagPerf[0].GetArray("udsg").GetArray(),
                                               BTagPerf[0].GetArray("bErr").GetArray(),BTagPerf[0].GetArray("udsgErr").GetArray());
  TGraph *discBTagger_udsg   = new TGraph(BTagPerf[0].GetN(),
                                   BTagPerf[0].GetArray("udsg").GetArray(),
                                   BTagPerf[0].GetArray("discriminator").GetArray());

  BTagger_b->SetName(bname[0].c_str());
  BTagger_c->SetName(cname[0].c_str());
  BTagger_udsg->SetName(udsgname[0].c_str());
  discBTagger_udsg->SetName(discname[0].c_str());

  BTagger_b->SetTitle("Jet b-efficiency");
  BTagger_c->SetTitle("Jet c-mistagging");
  BTagger_udsg->SetTitle("Jet udsg-mistagging");
  discBTagger_udsg->SetTitle("discriminator vs udsg-mistagging");


  for (int i=1;i<10;i++){
   graphcontainer_[discname[i]]      ->Write();
   grapherrorscontainer_[bname[i]]   ->Write();
   grapherrorscontainer_[cname[i]]   ->Write();
   grapherrorscontainer_[udsgname[i]]->Write();
  }

  BTagger_b->Write();
  BTagger_c->Write();
  BTagger_udsg->Write();
  discBTagger_udsg->Write();


}

Member Data Documentation

std::string PatBTagAnalyzer::bname[10]

private

Definition at line 49 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

double PatBTagAnalyzer::BTagdisccut_

private

Definition at line 44 of file PatBTagAnalyzer.cc.

Referenced by beginJob().

double PatBTagAnalyzer::BTagdiscmax_

private

Definition at line 45 of file PatBTagAnalyzer.cc.

Referenced by beginJob().

std::string PatBTagAnalyzer::BTagdiscriminator_

private

Definition at line 41 of file PatBTagAnalyzer.cc.

Referenced by analyze().

bJetSelector PatBTagAnalyzer::BTagger

private

Definition at line 67 of file PatBTagAnalyzer.cc.

PatBTagCommonHistos PatBTagAnalyzer::BTagHistograms

private

Definition at line 68 of file PatBTagAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::string PatBTagAnalyzer::BTagmethod_

private

Definition at line 40 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

BTagPerformance PatBTagAnalyzer::BTagPerf[10]

private

Definition at line 51 of file PatBTagAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

std::string PatBTagAnalyzer::BTagpurity_

private

Definition at line 39 of file PatBTagAnalyzer.cc.

bool PatBTagAnalyzer::BTagverbose

private

Definition at line 43 of file PatBTagAnalyzer.cc.

std::string PatBTagAnalyzer::cname[10]

private

Definition at line 50 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

std::string PatBTagAnalyzer::discname[10]

private

Definition at line 48 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

std::map<std::string,TGraph*> PatBTagAnalyzer::graphcontainer_

private

simple map to contain all graphs; graphs are booked in beginJob()

Definition at line 62 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

std::map<std::string,TGraphErrors*> PatBTagAnalyzer::grapherrorscontainer_

private

simple map to contain all graphs; graphs are booked in beginJob()

Definition at line 65 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().

std::map<std::string,TH2D*> PatBTagAnalyzer::h2_

private

simple map to contain 2D histograms; histograms are booked in beginJob()

Definition at line 59 of file PatBTagAnalyzer.cc.

std::map<std::string,TH1D*> PatBTagAnalyzer::histocontainer_

private

simple map to contain all histograms; histograms are booked in beginJob()

Definition at line 56 of file PatBTagAnalyzer.cc.

Referenced by beginJob().

edm::EDGetTokenT<edm::View<pat::Jet> > PatBTagAnalyzer::jetToken_

private

Definition at line 36 of file PatBTagAnalyzer.cc.

Referenced by analyze().

edm::ParameterSet PatBTagAnalyzer::PatBjet_

private

Definition at line 37 of file PatBTagAnalyzer.cc.

std::map<int,std::string> PatBTagAnalyzer::udsgname

private

Definition at line 52 of file PatBTagAnalyzer.cc.

Referenced by beginJob(), and endJob().