From ae73f4c98438212f45f4dcba7437424d65ae3e67 Mon Sep 17 00:00:00 2001 From: kleinb Date: Mon, 9 Jul 2012 12:31:02 +0000 Subject: [PATCH] take only good quality tracks as trigger --- PWGJE/AliAnalysisTaskJetCore.cxx | 21 ++++++++++++++------- PWGJE/AliAnalysisTaskJetCore.h | 2 ++ 2 files changed, 16 insertions(+), 7 deletions(-) diff --git a/PWGJE/AliAnalysisTaskJetCore.cxx b/PWGJE/AliAnalysisTaskJetCore.cxx index c44acb2f8e2..e8bc954e592 100644 --- a/PWGJE/AliAnalysisTaskJetCore.cxx +++ b/PWGJE/AliAnalysisTaskJetCore.cxx @@ -71,6 +71,7 @@ fVtxZMax(10.), fEvtClassMin(0), fEvtClassMax(4), fFilterMask(0), +fFilterType(2), fRadioFrac(0.2), fMinDist(0.1), fCentMin(0.), @@ -182,6 +183,7 @@ fVtxZMax(10.), fEvtClassMin(0), fEvtClassMax(4), fFilterMask(0), +fFilterType(2), fRadioFrac(0.2), fMinDist(0.1), fCentMin(0.), @@ -364,8 +366,8 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fh2JetCoreMethod1C60 = new TH2F("JetCoreMethod1C60","",150, 0., 150.,100, 0., 1.5); fh2JetCoreMethod2C60 = new TH2F("JetCoreMethod2C60","",150, 0., 150.,100, 0., 1.5);} - fh3JetTrackC3060=new TH3F("JetTrackC3060","",50,0,50,150,0.,150.,35,-0.5,3.5); - fh3JetTrackC20=new TH3F("JetTrackC20","",50,0,50,150,0.,150.,35,-0.5,3.5); + fh3JetTrackC3060=new TH3F("JetTrackC3060","",50,0,50,150,0.,150.,35,0.,3.5); + fh3JetTrackC20=new TH3F("JetTrackC20","",50,0,50,150,0.,150.,35,0.,3.5); if(fAngStructCloseTracks>0){ fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,10.); fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,10.); @@ -391,8 +393,8 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fh3JetDensity=new TH3F("Jet density vs mutliplicity A>0.4","",100,0.,4000.,100,0.,5.,10,0.,50.); fh3JetDensityA4=new TH3F("Jet density vs multiplicity A>0.4","",100,0.,4000.,100,0.,5.,10,0.,50.); - fh2RPJets=new TH2F("RPJet","",fNRPBins,0.,3.15,150,0.,150.); - fh2RPT=new TH2F("RPTrigger","",fNRPBins,0.,3.15,150,0.,150.); + fh2RPJets=new TH2F("RPJet","",fNRPBins,0.,fNRPBins,150,0.,150.); + fh2RPT=new TH2F("RPTrigger","",fNRPBins,0.,fNRPBins,150,0.,150.); fh3spectriggeredC20RP = new TH3F("Triggered spectrumC20RP","",3,0.,3.,140,-80.,200.,10,0.,50.); fh3spectriggeredC20 = new TH3F("Triggered spectrumC20","",100,0.,1.,140,-80.,200.,10,0.,50.); fh3spectriggeredC3060 = new TH3F("Triggered spectrumC3060","",100,0.,1.,140,-80.,200.,10,0.,50.); @@ -453,7 +455,7 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fOutputList->Add(fh3JetDensityA4); fOutputList->Add(fh2RPJets); fOutputList->Add(fh2RPT); - + fOutputList->Add(fh3spectriggeredC20RP); fOutputList->Add(fh3spectriggeredC20); fOutputList->Add(fh3spectriggeredC3060); @@ -654,7 +656,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) return;} fh2Ntriggers->Fill(centValue,partback->Pt()); - Int_t phiBinT = GetPhiBin(partback->Phi()-fRPAngle); + Int_t phiBinT = GetPhiBin(RelativePhi(partback->Phi(),fRPAngle)); if(centValue<20.) fh2RPT->Fill(phiBinT,partback->Pt()); Double_t accep=2.*TMath::Pi()*1.8; @@ -666,7 +668,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) phibig = jetbig->Phi(); ptbig = jetbig->Pt(); if(ptbig==0) continue; - Int_t phiBin = GetPhiBin(phibig-fRPAngle); + Int_t phiBin = GetPhiBin(RelativePhi(phibig,fRPAngle)); areabig = jetbig->EffectiveAreaCharged(); Double_t ptcorr=ptbig-rho*areabig; if((etabigfJetEtaMax)) continue; @@ -981,7 +983,12 @@ Int_t AliAnalysisTaskJetCore::GetListOfTracks(TList *list){ for(int it = 0;it < aod->GetNumberOfTracks();++it){ AliAODTrack *tr = aod->GetTrack(it); + Bool_t bGood = false; + if(fFilterType == 0)bGood = true; + else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal(); + else if(fFilterType == 2)bGood = tr->IsHybridGlobalConstrainedGlobal(); if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue; + if(bGood==false) continue; if(TMath::Abs(tr->Eta())>0.9)continue; if(tr->Pt()<0.15)continue; list->Add(tr); diff --git a/PWGJE/AliAnalysisTaskJetCore.h b/PWGJE/AliAnalysisTaskJetCore.h index 1fed31418fe..44192218d81 100644 --- a/PWGJE/AliAnalysisTaskJetCore.h +++ b/PWGJE/AliAnalysisTaskJetCore.h @@ -68,6 +68,7 @@ public: virtual void SetEvtClassMin(Int_t evtClass) { fEvtClassMin = evtClass; } virtual void SetEvtClassMax(Int_t evtClass) { fEvtClassMax = evtClass; } virtual void SetFilterMask(UInt_t i){fFilterMask = i;} + virtual void SetFilterType(Int_t iType){fFilterType=iType;} virtual void SetRadioFrac(Float_t radiofrac) { fRadioFrac = radiofrac; } virtual void SetMinDist(Float_t minDist) { fMinDist = minDist; } virtual void SetCentMin(Float_t cent) { fCentMin = cent; } @@ -119,6 +120,7 @@ private: Int_t fEvtClassMin; // lower bound on event class Int_t fEvtClassMax; // upper bound on event class UInt_t fFilterMask; // filter bit for slecected tracks + Int_t fFilterType; Float_t fRadioFrac; //!size of the concentric cone Float_t fMinDist; Float_t fCentMin; // lower bound on centrality -- 2.39.3