fTracks(),
hMEvents(NULL),
hTrackCent(NULL),
+ hTrigPt(NULL),
+ hTrackPt(NULL),
// fPhotonCorr(NULL),
// fPionCorr(NULL),
fDeltaAODBranchName("AliAODGammaConversion_gamma"),
MEHistograms->Add(hMEvents);
hTrackCent = new TH2I("hTrackCent", "N accepted tracks vs centrality",
- fAxisCent.GetNbins() > 2 ? 90 : 1, fAxisCent.GetBinLowEdge(1), fAxisCent.GetBinUpEdge(fAxisCent.GetNbins()),
- 750, 0, 1500);
+ fAxisCent.GetNbins() > 1 ? 10*(fAxisCent.GetXmax() - fAxisCent.GetXmin()) : 1,
+ fAxisCent.GetXmin(), fAxisCent.GetXmax(),
+ fAxisCent.GetNbins() > 1 ? 900 : 50, 0,
+ fAxisCent.GetNbins() > 1 ? 1800 : 50);
MEHistograms->Add(hTrackCent);
+ hTrigPt = new TH2F("hTrigPt", "trigger pt", 100, 0, 10, fAxisCent.GetNbins(), fAxisCent.GetXmin(), fAxisCent.GetXmax());
+ MEHistograms->Add(hTrigPt);
+ hTrackPt = new TH2F("hTrackPt", "track pt", 100, 0, 10, fAxisCent.GetNbins(), fAxisCent.GetXmin(), fAxisCent.GetXmax());
+ MEHistograms->Add(hTrackPt);
+
Int_t ntrackfilters[2] = {fTrackFilters[0].GetEntriesFast(), fTrackFilters[1].GetEntriesFast()};
fkTrackAxis = kTRUE;
fAxisTrackFilters.SetNameTitle("trackCuts", "trackCuts");
if(!photon) continue;
if(!fV0Filter || fV0Filter->PhotonIsSelected(photon, fInputEvent)) {
ggammas->Add(photon);
+ hTrigPt->Fill(photon->Pt(), centrality);
} else {
for(Int_t igf = 0; igf < fV0Filters[1].GetEntriesFast(); igf++) {
AliConversionCuts * gfilter = dynamic_cast<AliConversionCuts*>(fV0Filters[1].At(igf));
if(track->Pt() < aptlim[0] || track->Pt() > aptlim[1]) continue;
if(track->Eta() < aetalim[0] || track->Eta() > aetalim[1]) continue;
if(fTrackFilter->IsSelected(track)) {
+ hTrackPt->Fill(track->Pt(), centrality);
ttracks->Add(track);
} else {
///upside cuts
/// @file AliConversionTrackCuts.cxx
/// @author Svein Lindal
-/// @brief Base class for analysation of conversion particle - track correlations
+/// @brief Base class for analysation of conversion particle - track correlations - track cuts
#include "AliConversionTrackCuts.h"
fEsdTrackCutsExtra2(NULL),
fEvent(NULL),
fFilterBit(2048),
- fDCAZmax(-1),
- fDCAXYmax(-1),
+ fDCAZmax(3.2*3.2),
+ fDCAXYmax(2.4*2.4),
fOwnedTracks(),
fInitialized(kFALSE),
fhPhi(NULL),
const Int_t filterbit = fFilterBit;
if (filterbit == 128) {
- fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
- fEsdTrackCuts->SetMinNClustersTPC(70);
-
+ if(!fEsdTrackCuts) {
+ fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+ fEsdTrackCuts->SetMinNClustersTPC(70);
+ }
} else if (filterbit == 256) {
- // syst study
- fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
- fEsdTrackCuts->SetMinNClustersTPC(80);
- fEsdTrackCuts->SetMaxChi2PerClusterTPC(3);
- fEsdTrackCuts->SetMaxDCAToVertexZ(2.7);
- fEsdTrackCuts->SetMaxDCAToVertexXY(1.9);
-
+ if(!fEsdTrackCuts) {
+ // syst study
+ fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+ fEsdTrackCuts->SetMinNClustersTPC(80);
+ fEsdTrackCuts->SetMaxChi2PerClusterTPC(3);
+ fEsdTrackCuts->SetMaxDCAToVertexZ(2.7);
+ fEsdTrackCuts->SetMaxDCAToVertexXY(1.9);
+ }
} else if (filterbit == 512) {
- // syst study
- fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
- fEsdTrackCuts->SetMinNClustersTPC(60);
- fEsdTrackCuts->SetMaxChi2PerClusterTPC(5);
- fEsdTrackCuts->SetMaxDCAToVertexZ(3.7);
- fEsdTrackCuts->SetMaxDCAToVertexXY(2.9);
-
+ if(!fEsdTrackCuts) {
+ // syst study
+ fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+ fEsdTrackCuts->SetMinNClustersTPC(60);
+ fEsdTrackCuts->SetMaxChi2PerClusterTPC(5);
+ fEsdTrackCuts->SetMaxDCAToVertexZ(3.7);
+ fEsdTrackCuts->SetMaxDCAToVertexXY(2.9);
+ }
} else if (filterbit == 1024) {
if(!fEsdTrackCuts) {
fEsdTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
}
}
+
//______________________________________________________________________________
Bool_t AliConversionTrackCuts::AcceptTrack(AliESDtrack * track) {
//Check esd track
+ FillHistograms(kPreCut, track);
+
+
+ if( fFilterBit == 256) {
+
+ ///Standalone tpc tracks constrained
+ const AliExternalTrackParam * param = track->GetConstrainedParam();
+ if(param) {
+ AliESDtrack* esdTrack = new AliESDtrack(*track);
+ esdTrack->CopyFromVTrack(param);
+ track = esdTrack;
+ fOwnedTracks.Add(track);
+
+ if( !fEsdTrackCuts->IsSelected(track)) return kFALSE;
+
+ FillHistograms(1, track);
+
+ Double_t dca[2];
+ GetDCA(track, dca);
+
+ FillDCAHist(dca[1], dca[0], track);
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+
+ return kTRUE;
+ } else {
+ return kFALSE;
+ }
+
+ return kFALSE;
+ }
+
if(!fInitialized) {
DefineESDCuts();
- if(fDCAXYmax > 0) {
- if(fEsdTrackCuts) fEsdTrackCuts->SetMaxDCAToVertexXY(fDCAXYmax);
- }
- if(fDCAZmax > 0) {
- if(fEsdTrackCuts) fEsdTrackCuts->SetMaxDCAToVertexZ(fDCAZmax);
- }
+ // if(fDCAXYmax > 0) {
+ // if(fEsdTrackCuts) fEsdTrackCuts->SetMaxDCAToVertexXY(fDCAXYmax);
+ // }
+ // if(fDCAZmax > 0) {
+ // if(fEsdTrackCuts) fEsdTrackCuts->SetMaxDCAToVertexZ(fDCAZmax);
+ // }
fInitialized = kTRUE;
}
- FillHistograms(kPreCut, track);
+ Double_t dca[2];
+ GetDCA(track, dca);
- if( !fEsdTrackCuts->IsSelected(track)) return kFALSE;
-
-
///If only one track cuts then it has passed the cuts
if( !(fEsdTrackCutsExtra1 && fEsdTrackCutsExtra2)) {
FillHistograms(1, track);
+ FillDCAHist(dca[1], dca[0], track);
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
return kTRUE;
}
+
///If passing extra
if (fEsdTrackCutsExtra1 && fEsdTrackCutsExtra1->IsSelected(track)) {
- FillHistograms(2, track);
FillHistograms(1, track);
+ FillHistograms(2, track);
+
+ FillDCAHist(dca[1], dca[0], track);
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+
return kTRUE;
}
FillHistograms(3, track);
FillHistograms(1, track);
+
+ FillDCAHist(dca[1], dca[0], track);
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+
return kTRUE;
} else {
return kFALSE;
// ///Get impact parameters
// Double_t extCov[15];
// track->GetExternalCovariance(extCov);
- // Float_t b[2];
- // Float_t bCov[3];
- // track->GetImpactParameters(b,bCov);
- // if (bCov[0]<=0 || bCov[2]<=0) {
- // AliDebug(1, "Estimated b resolution lower or equal zero!");
- // bCov[0]=0; bCov[2]=0;
- // }
-
- // Float_t dcaToVertexXY = b[0];
- // Float_t dcaToVertexZ = b[1];
- // FillDCAHist(dcaToVertexZ, dcaToVertexXY, track);
// return kTRUE;
}
Bool_t AliConversionTrackCuts::AcceptTrack(AliAODTrack * track) {
//Check aod track
-
+
FillHistograms(kPreCut, track);
+
+ if (fFilterBit == 768) {
+ if(!track->IsHybridGlobalConstrainedGlobal()) return kFALSE;
+
+ if (!(track->GetStatus() & AliVTrack::kITSrefit)) {
+ return kFALSE;
+ }
+
+ //The cluster sharing cut can be done with:
+ Double_t frac = Double_t(track->GetTPCnclsS()) / Double_t(track->GetTPCncls());
+ if (frac > 0.4) return kFALSE;
+
+ ///Do dca xy cut!
+ FillHistograms(1, track);
+
+ ///DCA
+ Double_t dca[2] = { -999, -999};
+ //Bool_t dcaok =
+ GetDCA(track, dca);
+ FillDCAHist(dca[1], dca[0], track);
+
+
+ if(track->IsGlobalConstrained()) {
+ FillHistograms(3, track);
+ } else {
+ FillHistograms(2, track);
+ }
+
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+
+ return kTRUE;
+
+ ////////////////////////////////
+ //// Standalone
+ ////////////////////////////////
+ } else if(fFilterBit == 256) {
+ if(!track->IsTPCConstrained()) return kFALSE;
- if(!track->IsHybridGlobalConstrainedGlobal()) return kFALSE;
- if (!(track->GetStatus() & AliVTrack::kITSrefit)) {
- return kFALSE;
- }
-
+ ///DCA
+ Double_t dca[2] = { -999, -999};
+ GetDCA(track, dca);
- //The cluster sharing cut can be done with:
- Double_t frac = Double_t(track->GetTPCnclsS()) / Double_t(track->GetTPCncls());
- if (frac > 0.4) return kFALSE;
+ if( (dca[0]*dca[0]/fDCAXYmax + dca[1]*dca[1]/fDCAZmax) > 1 ) {
+ FillHistograms(3, track);
+ return kFALSE;
+ }
+ if(track->GetTPCncls() < 70) {
+ FillHistograms(4, track);
+ return kFALSE;
+ }
- ///Do dca xy cut!
- FillHistograms(1, track);
+ AliAODVertex * vtx = track->GetProdVertex();
+ if (vtx->GetType() == AliAODVertex::kKink ) {
+ FillHistograms(5, track);
+ return kFALSE;
+ }
- ///DCA
- Double_t dca[2] = { -999, -999};
- //Bool_t dcaok =
- GetDCA(track, dca);
- FillDCAHist(dca[1], dca[0], track);
-
+ if(track->Chi2perNDF() > 36) {
+ FillHistograms(6, track);
+ return kFALSE;
+ }
+ if(track->Chi2perNDF() > 26) {
+ FillHistograms(7, track);
+ return kFALSE;
+ }
+ if(track->Chi2perNDF() > 16) {
+ FillHistograms(8, track);
+ return kFALSE;
+ }
+ if(track->Chi2perNDF() > 4) {
+ FillHistograms(9, track);
+ return kFALSE;
+ }
+
+
+
+ FillDCAHist(dca[1], dca[0], track);
- if(track->IsGlobalConstrained()) {
- FillHistograms(3, track);
- } else {
FillHistograms(2, track);
+ if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+ return kTRUE;
+
}
-
- if(fhEtaPhi) fhEtaPhi->Fill(track->Eta(), track->Phi());
+ return kFALSE;
+}
+
+///______________________________________________________________________________
+Bool_t AliConversionTrackCuts::GetDCA(const AliESDtrack *track, Double_t dcaxyz[2]) {
+ ///Get track dca esd trck
+ Float_t dca[2];
+ Float_t bCov[3];
+ track->GetImpactParameters(dca,bCov);
+ if (bCov[0]<=0 || bCov[2]<=0) {
+ AliDebug(1, "Estimated b resolution lower or equal zero!");
+ bCov[0]=0; bCov[2]=0;
+ return kFALSE;
+ }
+
+ dcaxyz[0] = dca[0];
+ dcaxyz[1] = dca[1];
+
return kTRUE;
}
-///______________________________________________________________________________
+///_____________________________________________________________________________
Bool_t AliConversionTrackCuts::GetDCA(const AliAODTrack *track, Double_t dca[2]) {
+ ///Get track dca aod trck
if(track->TestBit(AliAODTrack::kIsDCA)){
dca[0]=track->DCA();
dca[1]=track->ZAtDCA();
Float_t xstart = etp.GetX();
if(xstart>3.) {
dca[0]=-999.;
- dca[1]=-999.;
+ dca[1]=-999.;
//printf("This method can be used only for propagation inside the beam pipe \n");
return kFALSE;
}
-
-
-
-
-
//_________________________________________________________________________________________________
void AliConversionTrackCuts::Print(const Option_t *) const
{