[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliAnalysisTaskSEBkgLikeSignJPSI.cxx

/**************************************************************************
 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

///////////////////////////////////////////////////////////////////////////
//
// AliAnalysisTaskSE for reading both reconstructed JPSI -> ee candidates
// and like sign pairs and for drawing corresponding distributions
//
// Author: C.Di Giglio, carmelo.digiglio@ba.infn.it
///////////////////////////////////////////////////////////////////////////

#include <TSystem.h>
#include <TROOT.h>
#include <TClonesArray.h>
#include <TNtuple.h>
#include <TList.h>
#include <TH1F.h>

#include "AliAnalysisManager.h"
#include "AliAODHandler.h"
#include "AliAODEvent.h"
#include "AliAODVertex.h"
#include "AliAODTrack.h"
#include "AliAODRecoDecayHF2Prong.h"
#include "AliAnalysisVertexingHF.h"
#include "AliAnalysisTaskSE.h"
#include "AliAnalysisTaskSEBkgLikeSignJPSI.h"

ClassImp(AliAnalysisTaskSEBkgLikeSignJPSI)

//________________________________________________________________________
AliAnalysisTaskSEBkgLikeSignJPSI::AliAnalysisTaskSEBkgLikeSignJPSI():
AliAnalysisTaskSE(),
fOutput(0), 
fHistMassJPSI(0),
fHistMassLS(0),
fHistCtsJPSI(0),           
fHistCtsLS(0),
fHistCtsLSpos(0),
fHistCtsLSneg(0),
fHistCPtaJPSI(0),          
fHistCPtaLS(0),
fHistd0d0JPSI(0),          
fHistd0d0LS(0),
fHistDCAJPSI(0),           
fHistDCALS(0),
fVHF(0),
fTotPosPairs(0),
fTotNegPairs(0),
fLsNormalization(1.)
{
  //
  // Default constructor
  //
}

//________________________________________________________________________
AliAnalysisTaskSEBkgLikeSignJPSI::AliAnalysisTaskSEBkgLikeSignJPSI(const char *name):
AliAnalysisTaskSE(name),
fOutput(0),
fHistMassJPSI(0),
fHistMassLS(0),
fHistCtsJPSI(0),
fHistCtsLS(0),
fHistCtsLSpos(0),
fHistCtsLSneg(0),
fHistCPtaJPSI(0),
fHistCPtaLS(0),
fHistd0d0JPSI(0),
fHistd0d0LS(0),
fHistDCAJPSI(0),
fHistDCALS(0),
fVHF(0),
fTotPosPairs(0),
fTotNegPairs(0),
fLsNormalization(1.)
{
  //
  // Standard constructor
  //
  // Output slot #1 writes into a TList container
  DefineOutput(1,TList::Class());  //My private output
}

//________________________________________________________________________
AliAnalysisTaskSEBkgLikeSignJPSI::~AliAnalysisTaskSEBkgLikeSignJPSI()
{
  // Destructor
  if (fOutput) {
    delete fOutput;
    fOutput = 0;
  }

  if (fVHF) {
    delete fVHF;
    fVHF = 0;
  }

}
//________________________________________________________________________
void AliAnalysisTaskSEBkgLikeSignJPSI::Init()
{
  // Initialization

  if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI::Init() \n");

  gROOT->LoadMacro("ConfigVertexingHF.C");

  fVHF = (AliAnalysisVertexingHF*)gROOT->ProcessLine("ConfigVertexingHF()");
  fVHF->PrintStatus();

  return;
}

//________________________________________________________________________
void AliAnalysisTaskSEBkgLikeSignJPSI::UserCreateOutputObjects()
{
  // Create the output container
  //
  if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI::UserCreateOutputObjects() \n");

  // Several histograms are more conveniently managed in a TList
  fOutput = new TList();
  fOutput->SetOwner();

  fHistMassJPSI = new TH1F("fHistMassJPSI", "J/#Psi invariant mass; M [GeV]; Entries",200,2.8,3.25);
  fHistMassJPSI->Sumw2();
  fHistMassJPSI->SetMinimum(0);
  fOutput->Add(fHistMassJPSI);

  fHistMassLS = new TH1F("fHistMassLS", "Like sign pairs invariant mass; M [GeV]; Entries",200,2.8,3.25);
  fHistMassLS->Sumw2();
  fHistMassLS->SetMinimum(0);
  fOutput->Add(fHistMassLS);

  fHistCtsJPSI = new TH1F("fHistCtsJPSI", "J/#Psi cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
  fHistCtsJPSI->Sumw2();
  fHistCtsJPSI->SetMinimum(0);
  fOutput->Add(fHistCtsJPSI);

  fHistCtsLS = new TH1F("fHistCtsLS", "Like sign pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
  fHistCtsLS->Sumw2();
  fHistCtsLS->SetMinimum(0);
  fOutput->Add(fHistCtsLS);

  fHistCtsLSpos = new TH1F("fHistCtsLSpos", "Like sign ++ pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
  fHistCtsLSpos->Sumw2();
  fHistCtsLSpos->SetMinimum(0);
  fOutput->Add(fHistCtsLSpos);

  fHistCtsLSneg = new TH1F("fHistCtsLSneg", "Like sign -- pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
  fHistCtsLSneg->Sumw2();
  fHistCtsLSneg->SetMinimum(0);
  fOutput->Add(fHistCtsLSneg);

  fHistCPtaJPSI = new TH1F("fHistCPtaJPSI", "J/#Psi cosine of pointing angle; Cos#Theta_{point}; Entries",200,-1.,1.);
  fHistCPtaJPSI->Sumw2();
  fHistCPtaJPSI->SetMinimum(0);
  fOutput->Add(fHistCPtaJPSI);

  fHistCPtaLS = new TH1F("fHistCPtaLS", "Like sign pairs cosine of pointing angle; Cos#Theta_{point}; Entries",200,-1.,1.);
  fHistCPtaLS->Sumw2();
  fHistCPtaLS->SetMinimum(0);
  fOutput->Add(fHistCPtaLS);

  fHistd0d0JPSI = new TH1F("fHistd0d0JPSI", "J/#Psi product of impact parameters; d0xd0 [#mu m^{2}]; Entries",200,-100000.,100000.);
  fHistd0d0JPSI->Sumw2(); 
  fHistd0d0JPSI->SetMinimum(0);
  fOutput->Add(fHistd0d0JPSI);

  fHistd0d0LS = new TH1F("fHistd0d0LS", "Like sign pairs product of impact parameters; d0xd0 [#mu m^{2}]; Entries",200,-100000.,100000.);
  fHistd0d0LS->Sumw2();
  fHistd0d0LS->SetMinimum(0);
  fOutput->Add(fHistd0d0LS);

  fHistDCAJPSI = new TH1F("fHistDCAJPSI", "J/#Psi distance of closest approach; dca [10^{2}#mu m]; Entries",100,0.,5.);
  fHistDCAJPSI->Sumw2(); 
  fHistDCAJPSI->SetMinimum(0);
  fOutput->Add(fHistDCAJPSI);

  fHistDCALS = new TH1F("fHistDCALS", "Like sign pairs distance of closest approach; dca [10^{2}#mu m]; Entries",100,0.,5.);
  fHistDCALS->Sumw2(); 
  fHistDCALS->SetMinimum(0);
  fOutput->Add(fHistDCALS);

  return;
}

//________________________________________________________________________
void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t */*option*/)
{
  // Execute analysis for current event:
  // heavy flavor candidates association to MC truth
  
  AliAODEvent *aod = dynamic_cast<AliAODEvent*> (InputEvent());

  TClonesArray *arrayJPSItoEle = 0;
  TClonesArray *arrayLikeSign = 0;

  if(!aod && AODEvent() && IsStandardAOD()) {
    // In case there is an AOD handler writing a standard AOD, use the AOD 
    // event in memory rather than the input (ESD) event.    
    aod = dynamic_cast<AliAODEvent*> (AODEvent());
    // in this case the braches in the deltaAOD (AliAOD.VertexingHF.root)
    // have to taken from the AOD event hold by the AliAODExtension
    AliAODHandler* aodHandler = (AliAODHandler*) 
      ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
    if(aodHandler->GetExtensions()) {
      AliAODExtension *ext = (AliAODExtension*)aodHandler->GetExtensions()->FindObject("AliAOD.VertexingHF.root");
      AliAODEvent *aodFromExt = ext->GetAOD();
      // load Jpsi candidates   
      arrayJPSItoEle=(TClonesArray*)aodFromExt->GetList()->FindObject("JPSItoEle");
      // load like sign candidates
      arrayLikeSign=(TClonesArray*)aodFromExt->GetList()->FindObject("LikeSign2Prong");
    }
  } else {
    // load Jpsi candidates                                                   
    arrayJPSItoEle=(TClonesArray*)aod->GetList()->FindObject("JPSItoEle");
    // load like sign candidates
    arrayLikeSign=(TClonesArray*)aod->GetList()->FindObject("LikeSign2Prong");
  }


  if(!arrayJPSItoEle) {
    printf("AliAnalysisTaskSEBkgLikeSignJPSI::UserExec: JPSItoEle branch not found!\n");
    return;
  }
  if(!arrayLikeSign) {
    printf("AliAnalysisTaskSEBkgLikeSignJPSI::UserExec: LikeSign2Prong branch not found!\n");
    return;
  }

  // fix for temporary bug in ESDfilter 
  // the AODs with null vertex pointer didn't pass the PhysSel
  if(!aod->GetPrimaryVertex()) return;

  // AOD primary vertex
  AliAODVertex *vtx1 = (AliAODVertex*)aod->GetPrimaryVertex();

  // make trkIDtoEntry register (temporary)
  Int_t trkIDtoEntry[100000];
  for(Int_t it=0;it<aod->GetNumberOfTracks();it++) {
    AliAODTrack *track = aod->GetTrack(it);
    trkIDtoEntry[track->GetID()]=it;
  }

  // loop over Like sign candidates
  Int_t nPosPairs=0,nNegPairs=0;
  Int_t nLikeSign = arrayLikeSign->GetEntriesFast();
  if(fDebug>1) printf("+++\n+++Number of like sign pairs ---> %d \n+++\n", nLikeSign);

  for(Int_t iLikeSign = 0; iLikeSign < nLikeSign; iLikeSign++) {
    AliAODRecoDecayHF2Prong *d = (AliAODRecoDecayHF2Prong*)arrayLikeSign->UncheckedAt(iLikeSign);
    Bool_t unsetvtx=kFALSE;
    if(!d->GetOwnPrimaryVtx()) {
        d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
        unsetvtx=kTRUE;
    }
    Int_t okBtoJPSIls=0;
    if(d->SelectBtoJPSI(fVHF->GetBtoJPSICuts(),okBtoJPSIls)) {
       AliAODTrack *trk0 = (AliAODTrack*)d->GetDaughter(0);
       fHistMassLS->Fill(d->InvMassJPSIee());
       fHistCPtaLS->Fill(d->CosPointingAngle());
       fHistd0d0LS->Fill(1e8*d->Prodd0d0());
       fHistDCALS->Fill(100*d->GetDCA());
       //PostData(1,fOutput);
       if(!trk0) {
          trk0=aod->GetTrack(trkIDtoEntry[d->GetProngID(0)]);
          printf("references to standard AOD not available \n");
       }
       if((trk0->Charge())==1) {
          nPosPairs++;
          fHistCtsLS->Fill(d->CosThetaStar(0,443,11,11));
          fHistCtsLSpos->Fill(d->CosThetaStar(0,443,11,11));
          //PostData(1,fOutput);
        } else {
          nNegPairs++;
          fHistCtsLS->Fill(d->CosThetaStarJPSI());
          fHistCtsLSneg->Fill(d->CosThetaStarJPSI());
          //PostData(1,fOutput);
        }
       PostData(1,fOutput);
    }
    if(unsetvtx) d->UnsetOwnPrimaryVtx();
  }

  if(fDebug>1) printf("------------ N. of positive pairs in Event ----- %d \n", nPosPairs);
  if(fDebug>1) printf("------------ N. of negative pairs in Event ----- %d \n", nNegPairs);

  fTotPosPairs += nPosPairs;
  fTotNegPairs += nNegPairs;

  // loop over JPSI candidates
  Int_t nBtoJpsiToEle = arrayJPSItoEle->GetEntriesFast();
  if(fDebug>1) printf("Number of like JPSI -> ee candidates ---> %d \n", nBtoJpsiToEle);

  for (Int_t iBtoJpsiToEle = 0; iBtoJpsiToEle < nBtoJpsiToEle; iBtoJpsiToEle++) {
    AliAODRecoDecayHF2Prong *d = (AliAODRecoDecayHF2Prong*)arrayJPSItoEle->UncheckedAt(iBtoJpsiToEle);
    Bool_t unsetvtx=kFALSE;
    if(!d->GetOwnPrimaryVtx()) {
      d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
      unsetvtx=kTRUE;
    }
    Int_t okBtoJPSI=0;
    if(d->SelectBtoJPSI(fVHF->GetBtoJPSICuts(),okBtoJPSI)) {
      fHistMassJPSI->Fill(d->InvMassJPSIee());
      fHistCtsJPSI->Fill(d->CosThetaStarJPSI());
      fHistd0d0JPSI->Fill(1e8*d->Prodd0d0());
      fHistCPtaJPSI->Fill(d->CosPointingAngle());
      fHistDCAJPSI->Fill(100*d->GetDCA());
      PostData(1,fOutput);
    }
    if(unsetvtx) d->UnsetOwnPrimaryVtx();
  }

  return;
}

//________________________________________________________________________
void AliAnalysisTaskSEBkgLikeSignJPSI::Terminate(Option_t */*option*/)
{
  // Terminate analysis
  //
  if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI: Terminate() \n");

  fOutput = dynamic_cast<TList*> (GetOutputData(1));
  if (!fOutput) {     
    printf("ERROR: fOutput not available\n");
    return;
  }

  fLsNormalization = 2.*TMath::Sqrt(fTotPosPairs*fTotNegPairs); 

  fHistMassJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistMassJPSI"));
  fHistMassLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistMassLS"));
  fHistCtsJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsJPSI"));
  fHistCtsLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLS"));
  fHistCtsLSpos = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLSpos"));
  fHistCtsLSneg = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLSneg"));
  fHistCPtaJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCPtaJPSI"));
  fHistCPtaLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCPtaLS"));
  fHistd0d0JPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistd0d0JPSI"));
  fHistd0d0LS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistd0d0LS"));
  fHistDCAJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCAJPSI"));
  fHistDCALS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCALS"));

  if(fLsNormalization>0.) {
    fHistMassLS->Scale((1/fLsNormalization)*fHistMassLS->GetEntries());
    fHistCtsLS->Scale((1/fLsNormalization)*fHistCtsLS->GetEntries());
    fHistCtsLSpos->Scale((1/fLsNormalization)*fHistCtsLSpos->GetEntries());
    fHistCtsLSneg->Scale((1/fLsNormalization)*fHistCtsLSneg->GetEntries());
    fHistCPtaLS->Scale((1/fLsNormalization)*fHistCPtaLS->GetEntries());
    fHistd0d0LS->Scale((1/fLsNormalization)*fHistd0d0LS->GetEntries());
    fHistDCALS->Scale((1/fLsNormalization)*fHistDCALS->GetEntries());
  }

  return;
}
Commit	Line	Data
adef5d4e	1	/**************************************************************************
	2	* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	///////////////////////////////////////////////////////////////////////////
	17	//
	18	// AliAnalysisTaskSE for reading both reconstructed JPSI -> ee candidates
	19	// and like sign pairs and for drawing corresponding distributions
	20	//
	21	// Author: C.Di Giglio, carmelo.digiglio@ba.infn.it
	22	///////////////////////////////////////////////////////////////////////////
	23
	24	#include <TSystem.h>
	25	#include <TROOT.h>
	26	#include <TClonesArray.h>
cbfc2f04	27	#include <TNtuple.h>
	28	#include <TList.h>
	29	#include <TH1F.h>
adef5d4e	30
b557eb43	31	#include "AliAnalysisManager.h"
b557eb43	32	#include "AliAODHandler.h"
adef5d4e	33	#include "AliAODEvent.h"
	34	#include "AliAODVertex.h"
	35	#include "AliAODTrack.h"
	36	#include "AliAODRecoDecayHF2Prong.h"
cbfc2f04	37	#include "AliAnalysisVertexingHF.h"
adef5d4e	38	#include "AliAnalysisTaskSE.h"
	39	#include "AliAnalysisTaskSEBkgLikeSignJPSI.h"
	40
	41	ClassImp(AliAnalysisTaskSEBkgLikeSignJPSI)
	42
	43	//________________________________________________________________________
	44	AliAnalysisTaskSEBkgLikeSignJPSI::AliAnalysisTaskSEBkgLikeSignJPSI():
	45	AliAnalysisTaskSE(),
	46	fOutput(0),
	47	fHistMassJPSI(0),
	48	fHistMassLS(0),
	49	fHistCtsJPSI(0),
	50	fHistCtsLS(0),
	51	fHistCtsLSpos(0),
	52	fHistCtsLSneg(0),
	53	fHistCPtaJPSI(0),
	54	fHistCPtaLS(0),
	55	fHistd0d0JPSI(0),
	56	fHistd0d0LS(0),
	57	fHistDCAJPSI(0),
	58	fHistDCALS(0),
	59	fVHF(0),
	60	fTotPosPairs(0),
	61	fTotNegPairs(0),
	62	fLsNormalization(1.)
	63	{
	64	//
	65	// Default constructor
	66	//
adef5d4e	67	}
	68
	69	//________________________________________________________________________
	70	AliAnalysisTaskSEBkgLikeSignJPSI::AliAnalysisTaskSEBkgLikeSignJPSI(const char *name):
	71	AliAnalysisTaskSE(name),
	72	fOutput(0),
	73	fHistMassJPSI(0),
	74	fHistMassLS(0),
	75	fHistCtsJPSI(0),
	76	fHistCtsLS(0),
	77	fHistCtsLSpos(0),
	78	fHistCtsLSneg(0),
	79	fHistCPtaJPSI(0),
	80	fHistCPtaLS(0),
	81	fHistd0d0JPSI(0),
	82	fHistd0d0LS(0),
	83	fHistDCAJPSI(0),
	84	fHistDCALS(0),
	85	fVHF(0),
	86	fTotPosPairs(0),
	87	fTotNegPairs(0),
	88	fLsNormalization(1.)
	89	{
	90	//
9d6d35b0	91	// Standard constructor
adef5d4e	92	//
	93	// Output slot #1 writes into a TList container
	94	DefineOutput(1,TList::Class()); //My private output
	95	}
	96
	97	//________________________________________________________________________
	98	AliAnalysisTaskSEBkgLikeSignJPSI::~AliAnalysisTaskSEBkgLikeSignJPSI()
	99	{
	100	// Destructor
	101	if (fOutput) {
	102	delete fOutput;
	103	fOutput = 0;
	104	}
	105
	106	if (fVHF) {
	107	delete fVHF;
	108	fVHF = 0;
	109	}
	110
	111	}
	112	//________________________________________________________________________
	113	void AliAnalysisTaskSEBkgLikeSignJPSI::Init()
	114	{
	115	// Initialization
	116
	117	if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI::Init() \n");
	118
	119	gROOT->LoadMacro("ConfigVertexingHF.C");
	120
	121	fVHF = (AliAnalysisVertexingHF*)gROOT->ProcessLine("ConfigVertexingHF()");
	122	fVHF->PrintStatus();
	123
	124	return;
adef5d4e	125	}
	126
	127	//________________________________________________________________________
	128	void AliAnalysisTaskSEBkgLikeSignJPSI::UserCreateOutputObjects()
	129	{
	130	// Create the output container
	131	//
	132	if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI::UserCreateOutputObjects() \n");
	133
	134	// Several histograms are more conveniently managed in a TList
	135	fOutput = new TList();
	136	fOutput->SetOwner();
	137
	138	fHistMassJPSI = new TH1F("fHistMassJPSI", "J/#Psi invariant mass; M [GeV]; Entries",200,2.8,3.25);
	139	fHistMassJPSI->Sumw2();
	140	fHistMassJPSI->SetMinimum(0);
	141	fOutput->Add(fHistMassJPSI);
	142
	143	fHistMassLS = new TH1F("fHistMassLS", "Like sign pairs invariant mass; M [GeV]; Entries",200,2.8,3.25);
	144	fHistMassLS->Sumw2();
	145	fHistMassLS->SetMinimum(0);
	146	fOutput->Add(fHistMassLS);
	147
	148	fHistCtsJPSI = new TH1F("fHistCtsJPSI", "J/#Psi cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
	149	fHistCtsJPSI->Sumw2();
	150	fHistCtsJPSI->SetMinimum(0);
	151	fOutput->Add(fHistCtsJPSI);
	152
	153	fHistCtsLS = new TH1F("fHistCtsLS", "Like sign pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
	154	fHistCtsLS->Sumw2();
	155	fHistCtsLS->SetMinimum(0);
	156	fOutput->Add(fHistCtsLS);
	157
	158	fHistCtsLSpos = new TH1F("fHistCtsLSpos", "Like sign ++ pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
	159	fHistCtsLSpos->Sumw2();
	160	fHistCtsLSpos->SetMinimum(0);
	161	fOutput->Add(fHistCtsLSpos);
	162
	163	fHistCtsLSneg = new TH1F("fHistCtsLSneg", "Like sign -- pairs cosine of decay angle; Cos#Theta*; Entries",200,-1.,1.);
	164	fHistCtsLSneg->Sumw2();
	165	fHistCtsLSneg->SetMinimum(0);
	166	fOutput->Add(fHistCtsLSneg);
	167
	168	fHistCPtaJPSI = new TH1F("fHistCPtaJPSI", "J/#Psi cosine of pointing angle; Cos#Theta_{point}; Entries",200,-1.,1.);
	169	fHistCPtaJPSI->Sumw2();
	170	fHistCPtaJPSI->SetMinimum(0);
	171	fOutput->Add(fHistCPtaJPSI);
	172
	173	fHistCPtaLS = new TH1F("fHistCPtaLS", "Like sign pairs cosine of pointing angle; Cos#Theta_{point}; Entries",200,-1.,1.);
	174	fHistCPtaLS->Sumw2();
	175	fHistCPtaLS->SetMinimum(0);
	176	fOutput->Add(fHistCPtaLS);
	177
	178	fHistd0d0JPSI = new TH1F("fHistd0d0JPSI", "J/#Psi product of impact parameters; d0xd0 [#mu m^{2}]; Entries",200,-100000.,100000.);
	179	fHistd0d0JPSI->Sumw2();
	180	fHistd0d0JPSI->SetMinimum(0);
	181	fOutput->Add(fHistd0d0JPSI);
	182
	183	fHistd0d0LS = new TH1F("fHistd0d0LS", "Like sign pairs product of impact parameters; d0xd0 [#mu m^{2}]; Entries",200,-100000.,100000.);
	184	fHistd0d0LS->Sumw2();
	185	fHistd0d0LS->SetMinimum(0);
	186	fOutput->Add(fHistd0d0LS);
	187
	188	fHistDCAJPSI = new TH1F("fHistDCAJPSI", "J/#Psi distance of closest approach; dca [10^{2}#mu m]; Entries",100,0.,5.);
189	fHistDCAJPSI->Sumw2();
190	fHistDCAJPSI->SetMinimum(0);
191	fOutput->Add(fHistDCAJPSI);
192
193	fHistDCALS = new TH1F("fHistDCALS", "Like sign pairs distance of closest approach; dca [10^{2}#mu m]; Entries",100,0.,5.);
194	fHistDCALS->Sumw2();
195	fHistDCALS->SetMinimum(0);
196	fOutput->Add(fHistDCALS);
197
198	return;
199	}
200
201	//________________________________________________________________________
202	void AliAnalysisTaskSEBkgLikeSignJPSI::UserExec(Option_t /option*/)
203	{
204	// Execute analysis for current event:
205	// heavy flavor candidates association to MC truth
206
207	AliAODEvent aod = dynamic_cast<AliAODEvent> (InputEvent());
208
b557eb43	209	TClonesArray *arrayJPSItoEle = 0;
	210	TClonesArray *arrayLikeSign = 0;
	211
	212	if(!aod && AODEvent() && IsStandardAOD()) {
	213	// In case there is an AOD handler writing a standard AOD, use the AOD
	214	// event in memory rather than the input (ESD) event.
	215	aod = dynamic_cast<AliAODEvent*> (AODEvent());
	216	// in this case the braches in the deltaAOD (AliAOD.VertexingHF.root)
	217	// have to taken from the AOD event hold by the AliAODExtension
	218	AliAODHandler* aodHandler = (AliAODHandler*)
	219	((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
	220	if(aodHandler->GetExtensions()) {
	221	AliAODExtension ext = (AliAODExtension)aodHandler->GetExtensions()->FindObject("AliAOD.VertexingHF.root");
	222	AliAODEvent *aodFromExt = ext->GetAOD();
	223	// load Jpsi candidates
	224	arrayJPSItoEle=(TClonesArray*)aodFromExt->GetList()->FindObject("JPSItoEle");
	225	// load like sign candidates
	226	arrayLikeSign=(TClonesArray*)aodFromExt->GetList()->FindObject("LikeSign2Prong");
	227	}
	228	} else {
	229	// load Jpsi candidates
	230	arrayJPSItoEle=(TClonesArray*)aod->GetList()->FindObject("JPSItoEle");
	231	// load like sign candidates
	232	arrayLikeSign=(TClonesArray*)aod->GetList()->FindObject("LikeSign2Prong");
adef5d4e	233	}
adef5d4e	234
b557eb43	235
adef5d4e	236	if(!arrayJPSItoEle) {
	237	printf("AliAnalysisTaskSEBkgLikeSignJPSI::UserExec: JPSItoEle branch not found!\n");
	238	return;
	239	}
adef5d4e	240	if(!arrayLikeSign) {
	241	printf("AliAnalysisTaskSEBkgLikeSignJPSI::UserExec: LikeSign2Prong branch not found!\n");
	242	return;
	243	}
	244
7c23877d	245	// fix for temporary bug in ESDfilter
	246	// the AODs with null vertex pointer didn't pass the PhysSel
	247	if(!aod->GetPrimaryVertex()) return;
	248
adef5d4e	249	// AOD primary vertex
	250	AliAODVertex vtx1 = (AliAODVertex)aod->GetPrimaryVertex();
	251
	252	// make trkIDtoEntry register (temporary)
	253	Int_t trkIDtoEntry[100000];
	254	for(Int_t it=0;it<aod->GetNumberOfTracks();it++) {
	255	AliAODTrack *track = aod->GetTrack(it);
	256	trkIDtoEntry[track->GetID()]=it;
	257	}
	258
	259	// loop over Like sign candidates
	260	Int_t nPosPairs=0,nNegPairs=0;
	261	Int_t nLikeSign = arrayLikeSign->GetEntriesFast();
3b1598da	262	if(fDebug>1) printf("+++\n+++Number of like sign pairs ---> %d \n+++\n", nLikeSign);
adef5d4e	263
	264	for(Int_t iLikeSign = 0; iLikeSign < nLikeSign; iLikeSign++) {
	265	AliAODRecoDecayHF2Prong d = (AliAODRecoDecayHF2Prong)arrayLikeSign->UncheckedAt(iLikeSign);
	266	Bool_t unsetvtx=kFALSE;
	267	if(!d->GetOwnPrimaryVtx()) {
	268	d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
	269	unsetvtx=kTRUE;
	270	}
	271	Int_t okBtoJPSIls=0;
	272	if(d->SelectBtoJPSI(fVHF->GetBtoJPSICuts(),okBtoJPSIls)) {
1fd36d27	273	AliAODTrack trk0 = (AliAODTrack)d->GetDaughter(0);
adef5d4e	274	fHistMassLS->Fill(d->InvMassJPSIee());
	275	fHistCPtaLS->Fill(d->CosPointingAngle());
	276	fHistd0d0LS->Fill(1e8*d->Prodd0d0());
adef5d4e	277	fHistDCALS->Fill(100*d->GetDCA());
1fd36d27	278	//PostData(1,fOutput);
adef5d4e	279	if(!trk0) {
	280	trk0=aod->GetTrack(trkIDtoEntry[d->GetProngID(0)]);
	281	printf("references to standard AOD not available \n");
	282	}
	283	if((trk0->Charge())==1) {
	284	nPosPairs++;
1fd36d27	285	fHistCtsLS->Fill(d->CosThetaStar(0,443,11,11));
	286	fHistCtsLSpos->Fill(d->CosThetaStar(0,443,11,11));
	287	//PostData(1,fOutput);
adef5d4e	288	} else {
adef5d4e	289	nNegPairs++;
1fd36d27	290	fHistCtsLS->Fill(d->CosThetaStarJPSI());
adef5d4e	291	fHistCtsLSneg->Fill(d->CosThetaStarJPSI());
1fd36d27	292	//PostData(1,fOutput);
adef5d4e	293	}
cbfc2f04	294	PostData(1,fOutput);
cbfc2f04	295	}
adef5d4e	296	if(unsetvtx) d->UnsetOwnPrimaryVtx();
	297	}
	298
3b1598da	299	if(fDebug>1) printf("------------ N. of positive pairs in Event ----- %d \n", nPosPairs);
3b1598da	300	if(fDebug>1) printf("------------ N. of negative pairs in Event ----- %d \n", nNegPairs);
adef5d4e	301
	302	fTotPosPairs += nPosPairs;
	303	fTotNegPairs += nNegPairs;
	304
	305	// loop over JPSI candidates
	306	Int_t nBtoJpsiToEle = arrayJPSItoEle->GetEntriesFast();
3b1598da	307	if(fDebug>1) printf("Number of like JPSI -> ee candidates ---> %d \n", nBtoJpsiToEle);
adef5d4e	308
	309	for (Int_t iBtoJpsiToEle = 0; iBtoJpsiToEle < nBtoJpsiToEle; iBtoJpsiToEle++) {
	310	AliAODRecoDecayHF2Prong d = (AliAODRecoDecayHF2Prong)arrayJPSItoEle->UncheckedAt(iBtoJpsiToEle);
	311	Bool_t unsetvtx=kFALSE;
	312	if(!d->GetOwnPrimaryVtx()) {
	313	d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
	314	unsetvtx=kTRUE;
	315	}
	316	Int_t okBtoJPSI=0;
	317	if(d->SelectBtoJPSI(fVHF->GetBtoJPSICuts(),okBtoJPSI)) {
	318	fHistMassJPSI->Fill(d->InvMassJPSIee());
	319	fHistCtsJPSI->Fill(d->CosThetaStarJPSI());
	320	fHistd0d0JPSI->Fill(1e8*d->Prodd0d0());
	321	fHistCPtaJPSI->Fill(d->CosPointingAngle());
	322	fHistDCAJPSI->Fill(100*d->GetDCA());
	323	PostData(1,fOutput);
	324	}
cbfc2f04	325	if(unsetvtx) d->UnsetOwnPrimaryVtx();
adef5d4e	326	}
	327
	328	return;
	329	}
	330
	331	//________________________________________________________________________
	332	void AliAnalysisTaskSEBkgLikeSignJPSI::Terminate(Option_t /option*/)
	333	{
	334	// Terminate analysis
	335	//
	336	if(fDebug > 1) printf("AnalysisTaskSEBkgLikeSignJPSI: Terminate() \n");
	337
	338	fOutput = dynamic_cast<TList*> (GetOutputData(1));
	339	if (!fOutput) {
	340	printf("ERROR: fOutput not available\n");
	341	return;
	342	}
	343
	344	fLsNormalization = 2.TMath::Sqrt(fTotPosPairsfTotNegPairs);
	345
	346	fHistMassJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistMassJPSI"));
	347	fHistMassLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistMassLS"));
	348	fHistCtsJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsJPSI"));
	349	fHistCtsLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLS"));
	350	fHistCtsLSpos = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLSpos"));
	351	fHistCtsLSneg = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCtsLSneg"));
	352	fHistCPtaJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCPtaJPSI"));
	353	fHistCPtaLS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistCPtaLS"));
	354	fHistd0d0JPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistd0d0JPSI"));
	355	fHistd0d0LS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistd0d0LS"));
	356	fHistDCAJPSI = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCAJPSI"));
	357	fHistDCALS = dynamic_cast<TH1F*>(fOutput->FindObject("fHistDCALS"));
	358
7fcc8b15	359	if(fLsNormalization>0.) {
	360	fHistMassLS->Scale((1/fLsNormalization)*fHistMassLS->GetEntries());
	361	fHistCtsLS->Scale((1/fLsNormalization)*fHistCtsLS->GetEntries());
	362	fHistCtsLSpos->Scale((1/fLsNormalization)*fHistCtsLSpos->GetEntries());
	363	fHistCtsLSneg->Scale((1/fLsNormalization)*fHistCtsLSneg->GetEntries());
	364	fHistCPtaLS->Scale((1/fLsNormalization)*fHistCPtaLS->GetEntries());
	365	fHistd0d0LS->Scale((1/fLsNormalization)*fHistd0d0LS->GetEntries());
	366	fHistDCALS->Scale((1/fLsNormalization)*fHistDCALS->GetEntries());
	367	}
adef5d4e	368
	369	return;
	370	}