[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskJetMassResponseDet.cxx

//
// Detector response jet mass analysis task.
//
// Author: M.Verweij

#include <TClonesArray.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <THnSparse.h>
#include <TList.h>
#include <TLorentzVector.h>
#include <TProfile.h>
#include <TChain.h>
#include <TSystem.h>
#include <TFile.h>
#include <TKey.h>

#include "AliVCluster.h"
#include "AliVTrack.h"
#include "AliEmcalJet.h"
#include "AliRhoParameter.h"
#include "AliLog.h"
#include "AliEmcalParticle.h"
#include "AliMCEvent.h"
#include "AliGenPythiaEventHeader.h"
#include "AliAODMCHeader.h"
#include "AliMCEvent.h"
#include "AliAnalysisManager.h"
#include "AliJetContainer.h"

#include "AliAODEvent.h"

#include "AliAnalysisTaskJetMassResponseDet.h"

ClassImp(AliAnalysisTaskJetMassResponseDet)

//________________________________________________________________________
AliAnalysisTaskJetMassResponseDet::AliAnalysisTaskJetMassResponseDet() : 
  AliAnalysisTaskEmcalJet("AliAnalysisTaskJetMassResponseDet", kTRUE),
  fContainerPart(0),
  fContainerDet(0),
  fJetMassType(kRaw),
  fh2PtVsMassJetPartAll(0),
  fh2PtVsMassJetPartMatch(0),
  fh2PtVsMassJetPartTagged(0),
  fh2PtVsMassJetPartTaggedMatch(0),
  fh2PtVsMassJetDetAll(0),
  fh2PtVsMassJetDetTagged(0),
  fh2EtaPhiMatchedDet(0),
  fh2EtaPhiMatchedPart(0),
  fhnMassResponse(0),
  fh1AreaPartAll(0),
  fh1AreaDetAll(0)
{
  // Default constructor.

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
AliAnalysisTaskJetMassResponseDet::AliAnalysisTaskJetMassResponseDet(const char *name) : 
  AliAnalysisTaskEmcalJet(name, kTRUE),  
  fContainerPart(0),
  fContainerDet(0),
  fJetMassType(kRaw),
  fh2PtVsMassJetPartAll(0),
  fh2PtVsMassJetPartMatch(0),
  fh2PtVsMassJetPartTagged(0),
  fh2PtVsMassJetPartTaggedMatch(0),
  fh2PtVsMassJetDetAll(0),
  fh2PtVsMassJetDetTagged(0),
  fh2EtaPhiMatchedDet(0),
  fh2EtaPhiMatchedPart(0),
  fhnMassResponse(0),
  fh1AreaPartAll(0),
  fh1AreaDetAll(0)
{
  // Standard constructor.

  SetMakeGeneralHistograms(kTRUE);
}

//________________________________________________________________________
AliAnalysisTaskJetMassResponseDet::~AliAnalysisTaskJetMassResponseDet()
{
  // Destructor.
}

//________________________________________________________________________
void AliAnalysisTaskJetMassResponseDet::UserCreateOutputObjects()
{
  // Create user output.

  AliAnalysisTaskEmcalJet::UserCreateOutputObjects();

  Bool_t oldStatus = TH1::AddDirectoryStatus();
  TH1::AddDirectory(kFALSE);

  const Int_t nBinsPt  = 200;
  const Double_t minPt = 0.;
  const Double_t maxPt = 200.;

  const Int_t nBinsM  = 100;
  const Double_t minM = 0.;
  const Double_t maxM = 50.;

  const Int_t nBinsConstEff  = 40;
  const Double_t minConstEff = 0.;
  const Double_t maxConstEff = 2.;

  // const Int_t nBinsConst = 26;
  // const Double_t minConst = -5.5;
  // const Double_t maxConst = 20.5;

  //Binning for THnSparse
  const Int_t nBinsSparse0 = 5;
  const Int_t nBins0[nBinsSparse0] = {nBinsM,nBinsM,nBinsPt,nBinsPt,nBinsConstEff};
  const Double_t xmin0[nBinsSparse0]  = { minM, minM, minPt, minPt, minConstEff};
  const Double_t xmax0[nBinsSparse0]  = { maxM, maxM, maxPt, maxPt, maxConstEff};

  //Create histograms
  TString histName = "";
  TString histTitle = "";
  
  histName = "fh2PtVsMassJetPartAll";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetPartAll = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetPartAll);

  histName = "fh2PtVsMassJetPartMatch";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetPartMatch = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetPartMatch);

  histName = "fh2PtVsMassJetPartTagged";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetPartTagged = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetPartTagged);

  histName = "fh2PtVsMassJetPartTaggedMatch";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetPartTaggedMatch = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetPartTaggedMatch);

  histName = "fh2PtVsMassJetDetAll";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetDetAll = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetDetAll);

  histName = "fh2PtVsMassJetDetTagged";
  histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
  fh2PtVsMassJetDetTagged = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
  fOutput->Add(fh2PtVsMassJetDetTagged);

  histName = "fh2EtaPhiMatchedDet";
  histTitle = TString::Format("%s;#eta;#varphi",histName.Data());
  fh2EtaPhiMatchedDet = new TH2F(histName.Data(),histTitle.Data(),100,-1.,1.,72,0.,TMath::TwoPi());
  fOutput->Add(fh2EtaPhiMatchedDet);

  histName = "fh2EtaPhiMatchedPart";
  histTitle = TString::Format("%s;#eta;#varphi",histName.Data());
  fh2EtaPhiMatchedPart = new TH2F(histName.Data(),histTitle.Data(),100,-1.,1.,72,0.,TMath::TwoPi());
  fOutput->Add(fh2EtaPhiMatchedPart);

  histName = "fhnMassResponse";
  histTitle = Form("%s;#it{M}_{det};#it{M}_{part};#it{p}_{T,det};#it{p}_{T,part};#it{N}_{const}^{det}/#it{N}_{const}^{part}",histName.Data());
  fhnMassResponse = new THnSparseF(histName.Data(),histTitle.Data(),nBinsSparse0,nBins0,xmin0,xmax0);
  fOutput->Add(fhnMassResponse);

  fh1AreaPartAll = new TH1D("fh1AreaPartAll","fh1AreaPartAll",100.,0.,1.);
  fOutput->Add(fh1AreaPartAll);

  fh1AreaDetAll = new TH1D("fh1AreaDetAll","fh1AreaDetAll",100.,0.,1.);
  fOutput->Add(fh1AreaDetAll);


  // =========== Switch on Sumw2 for all histos ===========
  for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
    TH1 *h1 = dynamic_cast<TH1*>(fOutput->At(i));
    if (h1){
      h1->Sumw2();
      continue;
    }
    THnSparse *hn = dynamic_cast<THnSparse*>(fOutput->At(i));
    if(hn)hn->Sumw2();
  }

  TH1::AddDirectory(oldStatus);

  PostData(1, fOutput); // Post data for ALL output slots > 0 here.
}

//________________________________________________________________________
Bool_t AliAnalysisTaskJetMassResponseDet::Run()
{
  // Run analysis code here, if needed. It will be executed before FillHistograms().
  return kTRUE;
}

//________________________________________________________________________
Bool_t AliAnalysisTaskJetMassResponseDet::FillHistograms()
{
  // Fill histograms.

  AliJetContainer *cPart = GetJetContainer(fContainerPart);
  AliJetContainer *cDet = GetJetContainer(fContainerDet);
  AliEmcalJet* jPart = NULL;
  AliEmcalJet* jDet = NULL;

  //loop on particle level jets
  Int_t nAccPart = 0;
  if(cPart) {
    cPart->ResetCurrentID();
    while((jPart = cPart->GetNextAcceptJet())) {
      fh2PtVsMassJetPartAll->Fill(jPart->Pt(),jPart->M());
      fh1AreaPartAll->Fill(jPart->Area());
      nAccPart++;
      jDet = jPart->ClosestJet();
      if(jDet) fh2PtVsMassJetPartMatch->Fill(jPart->Pt(),jPart->M());
      if(jPart->GetTagStatus()<1 || !jPart->GetTaggedJet())
	continue;
      fh2PtVsMassJetPartTagged->Fill(jPart->Pt(),jPart->M());
      if(jDet) fh2PtVsMassJetPartTaggedMatch->Fill(jPart->Pt(),jPart->M());
    }
  }
  
  //loop on detector level jets
  Int_t nAccDet = 0;
  if(cDet) {
    cDet->ResetCurrentID();
    while((jDet = cDet->GetNextAcceptJet())) {
      Double_t mjet = GetJetMass(jDet);     
      fh2PtVsMassJetDetAll->Fill(jDet->Pt(),mjet);
      fh1AreaDetAll->Fill(jDet->Area());
      nAccDet++;
      if(jDet->GetTagStatus()>=1 && jDet->GetTaggedJet())
	 fh2PtVsMassJetDetTagged->Fill(jDet->Pt(),mjet);
       
       //fill detector response
       jPart = jDet->ClosestJet();
       if(jPart) {
	 fh2EtaPhiMatchedDet->Fill(jDet->Eta(),jDet->Phi());
	 fh2EtaPhiMatchedPart->Fill(jPart->Eta(),jPart->Phi());

	 Int_t nConstPart = jPart->GetNumberOfConstituents();
	 Int_t nConstDet = jDet->GetNumberOfConstituents();
	 Double_t eff = -1.;
	 if(nConstPart>0) eff = (Double_t)nConstDet/((Double_t)nConstPart);
	 Double_t var[5] = {GetJetMass(jDet),jPart->M(),jDet->Pt(),jPart->Pt(),eff};
	 fhnMassResponse->Fill(var);
       }
    }
  }

  return kTRUE;
}

//________________________________________________________________________
Double_t AliAnalysisTaskJetMassResponseDet::GetJetMass(AliEmcalJet *jet) {
  //calc subtracted jet mass
  if(fJetMassType==kRaw)
    return jet->M();
  else if(fJetMassType==kDeriv)
    return jet->GetSecondOrderSubtracted();
  
  return 0;
}

//________________________________________________________________________
Bool_t AliAnalysisTaskJetMassResponseDet::RetrieveEventObjects() {
  //
  // retrieve event objects
  if (!AliAnalysisTaskEmcalJet::RetrieveEventObjects())
    return kFALSE;

  return kTRUE;
}

//_______________________________________________________________________
void AliAnalysisTaskJetMassResponseDet::Terminate(Option_t *) 
{
  // Called once at the end of the analysis.
}
Commit	Line	Data
252cb00e	1	//
	2	// Detector response jet mass analysis task.
	3	//
	4	// Author: M.Verweij
	5
	6	#include <TClonesArray.h>
	7	#include <TH1F.h>
	8	#include <TH2F.h>
	9	#include <TH3F.h>
	10	#include <THnSparse.h>
	11	#include <TList.h>
	12	#include <TLorentzVector.h>
	13	#include <TProfile.h>
	14	#include <TChain.h>
	15	#include <TSystem.h>
	16	#include <TFile.h>
	17	#include <TKey.h>
	18
	19	#include "AliVCluster.h"
	20	#include "AliVTrack.h"
	21	#include "AliEmcalJet.h"
	22	#include "AliRhoParameter.h"
	23	#include "AliLog.h"
	24	#include "AliEmcalParticle.h"
	25	#include "AliMCEvent.h"
	26	#include "AliGenPythiaEventHeader.h"
	27	#include "AliAODMCHeader.h"
	28	#include "AliMCEvent.h"
	29	#include "AliAnalysisManager.h"
	30	#include "AliJetContainer.h"
	31
	32	#include "AliAODEvent.h"
	33
	34	#include "AliAnalysisTaskJetMassResponseDet.h"
	35
	36	ClassImp(AliAnalysisTaskJetMassResponseDet)
	37
	38	//________________________________________________________________________
	39	AliAnalysisTaskJetMassResponseDet::AliAnalysisTaskJetMassResponseDet() :
	40	AliAnalysisTaskEmcalJet("AliAnalysisTaskJetMassResponseDet", kTRUE),
	41	fContainerPart(0),
	42	fContainerDet(0),
	43	fJetMassType(kRaw),
	44	fh2PtVsMassJetPartAll(0),
	45	fh2PtVsMassJetPartMatch(0),
	46	fh2PtVsMassJetPartTagged(0),
	47	fh2PtVsMassJetPartTaggedMatch(0),
	48	fh2PtVsMassJetDetAll(0),
	49	fh2PtVsMassJetDetTagged(0),
b3821b05	50	fh2EtaPhiMatchedDet(0),
b3821b05	51	fh2EtaPhiMatchedPart(0),
3c923384	52	fhnMassResponse(0),
	53	fh1AreaPartAll(0),
	54	fh1AreaDetAll(0)
252cb00e	55	{
	56	// Default constructor.
	57
	58	SetMakeGeneralHistograms(kTRUE);
	59	}
	60
	61	//________________________________________________________________________
	62	AliAnalysisTaskJetMassResponseDet::AliAnalysisTaskJetMassResponseDet(const char *name) :
	63	AliAnalysisTaskEmcalJet(name, kTRUE),
	64	fContainerPart(0),
	65	fContainerDet(0),
	66	fJetMassType(kRaw),
	67	fh2PtVsMassJetPartAll(0),
	68	fh2PtVsMassJetPartMatch(0),
	69	fh2PtVsMassJetPartTagged(0),
	70	fh2PtVsMassJetPartTaggedMatch(0),
	71	fh2PtVsMassJetDetAll(0),
	72	fh2PtVsMassJetDetTagged(0),
b3821b05	73	fh2EtaPhiMatchedDet(0),
b3821b05	74	fh2EtaPhiMatchedPart(0),
3c923384	75	fhnMassResponse(0),
	76	fh1AreaPartAll(0),
	77	fh1AreaDetAll(0)
252cb00e	78	{
	79	// Standard constructor.
	80
	81	SetMakeGeneralHistograms(kTRUE);
	82	}
	83
	84	//________________________________________________________________________
	85	AliAnalysisTaskJetMassResponseDet::~AliAnalysisTaskJetMassResponseDet()
	86	{
	87	// Destructor.
	88	}
	89
	90	//________________________________________________________________________
	91	void AliAnalysisTaskJetMassResponseDet::UserCreateOutputObjects()
	92	{
	93	// Create user output.
	94
	95	AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
	96
	97	Bool_t oldStatus = TH1::AddDirectoryStatus();
	98	TH1::AddDirectory(kFALSE);
	99
	100	const Int_t nBinsPt = 200;
	101	const Double_t minPt = 0.;
	102	const Double_t maxPt = 200.;
	103
	104	const Int_t nBinsM = 100;
	105	const Double_t minM = 0.;
	106	const Double_t maxM = 50.;
	107
51ff3835	108	const Int_t nBinsConstEff = 40;
b3821b05	109	const Double_t minConstEff = 0.;
51ff3835	110	const Double_t maxConstEff = 2.;
252cb00e	111
b3821b05	112	// const Int_t nBinsConst = 26;
	113	// const Double_t minConst = -5.5;
	114	// const Double_t maxConst = 20.5;
1c0626ce	115
252cb00e	116	//Binning for THnSparse
252cb00e	117	const Int_t nBinsSparse0 = 5;
b3821b05	118	const Int_t nBins0[nBinsSparse0] = {nBinsM,nBinsM,nBinsPt,nBinsPt,nBinsConstEff};
	119	const Double_t xmin0[nBinsSparse0] = { minM, minM, minPt, minPt, minConstEff};
	120	const Double_t xmax0[nBinsSparse0] = { maxM, maxM, maxPt, maxPt, maxConstEff};
252cb00e	121
	122	//Create histograms
	123	TString histName = "";
	124	TString histTitle = "";
	125
	126	histName = "fh2PtVsMassJetPartAll";
	127	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	128	fh2PtVsMassJetPartAll = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	129	fOutput->Add(fh2PtVsMassJetPartAll);
	130
	131	histName = "fh2PtVsMassJetPartMatch";
	132	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	133	fh2PtVsMassJetPartMatch = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	134	fOutput->Add(fh2PtVsMassJetPartMatch);
	135
	136	histName = "fh2PtVsMassJetPartTagged";
	137	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	138	fh2PtVsMassJetPartTagged = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	139	fOutput->Add(fh2PtVsMassJetPartTagged);
	140
	141	histName = "fh2PtVsMassJetPartTaggedMatch";
	142	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	143	fh2PtVsMassJetPartTaggedMatch = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	144	fOutput->Add(fh2PtVsMassJetPartTaggedMatch);
	145
	146	histName = "fh2PtVsMassJetDetAll";
	147	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	148	fh2PtVsMassJetDetAll = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	149	fOutput->Add(fh2PtVsMassJetDetAll);
	150
	151	histName = "fh2PtVsMassJetDetTagged";
	152	histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1}",histName.Data());
	153	fh2PtVsMassJetDetTagged = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM);
	154	fOutput->Add(fh2PtVsMassJetDetTagged);
	155
b3821b05	156	histName = "fh2EtaPhiMatchedDet";
	157	histTitle = TString::Format("%s;#eta;#varphi",histName.Data());
	158	fh2EtaPhiMatchedDet = new TH2F(histName.Data(),histTitle.Data(),100,-1.,1.,72,0.,TMath::TwoPi());
	159	fOutput->Add(fh2EtaPhiMatchedDet);
	160
	161	histName = "fh2EtaPhiMatchedPart";
	162	histTitle = TString::Format("%s;#eta;#varphi",histName.Data());
	163	fh2EtaPhiMatchedPart = new TH2F(histName.Data(),histTitle.Data(),100,-1.,1.,72,0.,TMath::TwoPi());
	164	fOutput->Add(fh2EtaPhiMatchedPart);
	165
252cb00e	166	histName = "fhnMassResponse";
51ff3835	167	histTitle = Form("%s;#it{M}_{det};#it{M}_{part};#it{p}_{T,det};#it{p}_{T,part};#it{N}_{const}^{det}/#it{N}_{const}^{part}",histName.Data());
252cb00e	168	fhnMassResponse = new THnSparseF(histName.Data(),histTitle.Data(),nBinsSparse0,nBins0,xmin0,xmax0);
	169	fOutput->Add(fhnMassResponse);
	170
3c923384	171	fh1AreaPartAll = new TH1D("fh1AreaPartAll","fh1AreaPartAll",100.,0.,1.);
	172	fOutput->Add(fh1AreaPartAll);
	173
	174	fh1AreaDetAll = new TH1D("fh1AreaDetAll","fh1AreaDetAll",100.,0.,1.);
	175	fOutput->Add(fh1AreaDetAll);
	176
252cb00e	177
	178	// =========== Switch on Sumw2 for all histos ===========
	179	for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
	180	TH1 h1 = dynamic_cast<TH1>(fOutput->At(i));
	181	if (h1){
	182	h1->Sumw2();
	183	continue;
	184	}
	185	THnSparse hn = dynamic_cast<THnSparse>(fOutput->At(i));
	186	if(hn)hn->Sumw2();
	187	}
	188
	189	TH1::AddDirectory(oldStatus);
	190
	191	PostData(1, fOutput); // Post data for ALL output slots > 0 here.
	192	}
	193
	194	//________________________________________________________________________
	195	Bool_t AliAnalysisTaskJetMassResponseDet::Run()
	196	{
	197	// Run analysis code here, if needed. It will be executed before FillHistograms().
252cb00e	198	return kTRUE;
	199	}
	200
	201	//________________________________________________________________________
	202	Bool_t AliAnalysisTaskJetMassResponseDet::FillHistograms()
	203	{
	204	// Fill histograms.
	205
1c0626ce	206	AliJetContainer *cPart = GetJetContainer(fContainerPart);
1c0626ce	207	AliJetContainer *cDet = GetJetContainer(fContainerDet);
252cb00e	208	AliEmcalJet* jPart = NULL;
	209	AliEmcalJet* jDet = NULL;
	210
	211	//loop on particle level jets
3c923384	212	Int_t nAccPart = 0;
1c0626ce	213	if(cPart) {
	214	cPart->ResetCurrentID();
	215	while((jPart = cPart->GetNextAcceptJet())) {
252cb00e	216	fh2PtVsMassJetPartAll->Fill(jPart->Pt(),jPart->M());
3c923384	217	fh1AreaPartAll->Fill(jPart->Area());
3c923384	218	nAccPart++;
252cb00e	219	jDet = jPart->ClosestJet();
	220	if(jDet) fh2PtVsMassJetPartMatch->Fill(jPart->Pt(),jPart->M());
	221	if(jPart->GetTagStatus()<1 \|\| !jPart->GetTaggedJet())
	222	continue;
	223	fh2PtVsMassJetPartTagged->Fill(jPart->Pt(),jPart->M());
	224	if(jDet) fh2PtVsMassJetPartTaggedMatch->Fill(jPart->Pt(),jPart->M());
	225	}
	226	}
	227
	228	//loop on detector level jets
3c923384	229	Int_t nAccDet = 0;
1c0626ce	230	if(cDet) {
	231	cDet->ResetCurrentID();
	232	while((jDet = cDet->GetNextAcceptJet())) {
252cb00e	233	Double_t mjet = GetJetMass(jDet);
252cb00e	234	fh2PtVsMassJetDetAll->Fill(jDet->Pt(),mjet);
3c923384	235	fh1AreaDetAll->Fill(jDet->Area());
	236	nAccDet++;
	237	if(jDet->GetTagStatus()>=1 && jDet->GetTaggedJet())
252cb00e	238	fh2PtVsMassJetDetTagged->Fill(jDet->Pt(),mjet);
	239
	240	//fill detector response
	241	jPart = jDet->ClosestJet();
	242	if(jPart) {
b3821b05	243	fh2EtaPhiMatchedDet->Fill(jDet->Eta(),jDet->Phi());
	244	fh2EtaPhiMatchedPart->Fill(jPart->Eta(),jPart->Phi());
	245
	246	Int_t nConstPart = jPart->GetNumberOfConstituents();
	247	Int_t nConstDet = jDet->GetNumberOfConstituents();
b3821b05	248	Double_t eff = -1.;
	249	if(nConstPart>0) eff = (Double_t)nConstDet/((Double_t)nConstPart);
	250	Double_t var[5] = {GetJetMass(jDet),jPart->M(),jDet->Pt(),jPart->Pt(),eff};
252cb00e	251	fhnMassResponse->Fill(var);
	252	}
	253	}
	254	}
3c923384	255
252cb00e	256	return kTRUE;
	257	}
	258
	259	//________________________________________________________________________
	260	Double_t AliAnalysisTaskJetMassResponseDet::GetJetMass(AliEmcalJet *jet) {
	261	//calc subtracted jet mass
	262	if(fJetMassType==kRaw)
	263	return jet->M();
	264	else if(fJetMassType==kDeriv)
	265	return jet->GetSecondOrderSubtracted();
	266
	267	return 0;
	268	}
	269
	270	//________________________________________________________________________
	271	Bool_t AliAnalysisTaskJetMassResponseDet::RetrieveEventObjects() {
	272	//
	273	// retrieve event objects
	274	if (!AliAnalysisTaskEmcalJet::RetrieveEventObjects())
	275	return kFALSE;
	276
	277	return kTRUE;
	278	}
	279
	280	//_______________________________________________________________________
	281	void AliAnalysisTaskJetMassResponseDet::Terminate(Option_t *)
	282	{
	283	// Called once at the end of the analysis.
	284	}
	285