// AOD track cuts
- for(Int_t iTrackBit = 0; iTrackBit < gBitMax; iTrackBit++){
- fHistTrackStats->Fill(gCentrality,iTrackBit,aodTrack->TestFilterBit(1<<iTrackBit));
-
- if(aodTrack->TestFilterBit(1<<iTrackBit)){
- fHistKinematics[iCharge][iTrackBit]->Fill(vEta,vPhi,vPt);
- fHistDCAconstrained[iCharge][iTrackBit]->Fill(vDCAconstrainedxy,vDCAconstrainedz);
- fHistDCAglobal[iCharge][iTrackBit]->Fill(vDCAglobalx,vDCAglobaly,vDCAglobalz);
- fHistChiClus[iCharge][iTrackBit]->Fill(vChi2,vClus);
- }
-
- }//bit loop
+ if(iCharge > -1){
+ for(Int_t iTrackBit = 0; iTrackBit < gBitMax; iTrackBit++){
+ fHistTrackStats->Fill(gCentrality,iTrackBit,aodTrack->TestFilterBit(1<<iTrackBit));
+
+ if(aodTrack->TestFilterBit(1<<iTrackBit)){
+ fHistKinematics[iCharge][iTrackBit]->Fill(vEta,vPhi,vPt);
+ fHistDCAconstrained[iCharge][iTrackBit]->Fill(vDCAconstrainedxy,vDCAconstrainedz);
+ fHistDCAglobal[iCharge][iTrackBit]->Fill(vDCAglobalx,vDCAglobaly,vDCAglobalz);
+ fHistChiClus[iCharge][iTrackBit]->Fill(vChi2,vClus);
+ }
+ }//bit loop
+ }//charge positive or negative
}//track loop
-
+
return;
}
//________________________________________________________________________
AliEbyEMultFluctuationTask::AliEbyEMultFluctuationTask(const char *name)
-: AliAnalysisTaskSE(name),fAOD(0), fAODVertex(0),fHistNchPt(0),fHistNchEta(0),fHistNchEtaCent(0),fHistNchPhi(0),fHistDCAxy(0),fHistDCAz(0),fHistnclus(0),fHistchi2ndf(0),fHistchi2ndfvscs(0),fHistVz(0),fHistMultV0A(0),fHistMultV0C(0),fHistMultV0total(0),My_ntuple(0),fOutputList(0),fCentralityEstimator("V0M"),fCentralityCounter(0),fEventCounter(0),histcounter(0)
+: AliAnalysisTaskSE(name),fAOD(0), fAODVertex(0),fHistNchPt(0),fHistNchEta(0),fHistNchEtaCent(0),fHistNchPhi(0),fHistDCAxy(0),fHistDCAz(0),fHistnclus(0),fHistchi2ndf(0),fHistchi2ndfvscs(0),fHistVz(0),fHistMultV0A(0),fHistMultV0C(0),fHistMultV0total(0),My_ntuple(0),fOutputList(0),fCentralityEstimator("V0M"),fCentralityBins20(kFALSE),fCentralityCounter(0),fEventCounter(0),histcounter(0)
{
//________________________________________________________________________
Bool_t AliEbyEMultFluctuationTask :: SelectEvent(AliAODVertex* vertex)
{
-if(vertex)
-
- if(vertex->GetNContributors() < 0) return kFALSE;
-
+ if(vertex){
+
+ if(vertex->GetNContributors() < 0) return kFALSE;
+
+
Double_t lvx = vertex->GetX();
Double_t lvy = vertex->GetY();
Double_t lvz = vertex->GetZ();
if(TMath::Abs(lvx) > 0.3) return kFALSE;
- if(TMath::Abs(lvy) > 0.3) return kFALSE;
- if(TMath::Abs(lvz) > 10) return kFALSE;
- if(vertex->GetType()==AliAODVertex::kPileupSPD) return kFALSE;
-
- fEventCounter->Fill(5);
- fHistVz->Fill(lvz);
-return kTRUE;
+ if(TMath::Abs(lvy) > 0.3) return kFALSE;
+ if(TMath::Abs(lvz) > 10) return kFALSE;
+ if(vertex->GetType()==AliAODVertex::kPileupSPD) return kFALSE;
+
+ fEventCounter->Fill(5);
+ fHistVz->Fill(lvz);
+ return kTRUE;
+ }
+ else{
+ return kFALSE;
+ }
}
//_____________________________________________________________________________
Int_t AliEbyEMultFluctuationTask :: SelectTrack(AliAODTrack* track)
// Get the event vertex.
fAODVertex = fAOD->GetPrimaryVertex();
-
- if(TMath::Abs((fAODVertex->GetZ())-((fAOD->GetPrimaryVertexSPD())->GetZ())) > 0.5) return;
- if (!fAODVertex) {
-
- return;
+ if (!fAODVertex) {
+ return;
}
+
+ if(TMath::Abs((fAODVertex->GetZ())-((fAOD->GetPrimaryVertexSPD())->GetZ())) > 0.5) return;
if (!SelectEvent(fAODVertex)) return;
class AliEbyEMultFluctuationTask : public AliAnalysisTaskSE {
public:
- AliEbyEMultFluctuationTask() : AliAnalysisTaskSE(), fAOD(0), fAODVertex(0),fHistNchPt(0),fHistNchEta(0),fHistNchEtaCent(0),fHistNchPhi(0),fHistDCAxy(0),fHistDCAz(0),fHistnclus(0),fHistchi2ndf(0),fHistchi2ndfvscs(0),fHistVz(0),fHistMultV0A(0),fHistMultV0C(0),fHistMultV0total(0),My_ntuple(0),fOutputList(0),fCentralityEstimator("V0M"),fCentralityCounter(0),fEventCounter(0),histcounter(0){
+ AliEbyEMultFluctuationTask() : AliAnalysisTaskSE(), fAOD(0), fAODVertex(0),fHistNchPt(0),fHistNchEta(0),fHistNchEtaCent(0),fHistNchPhi(0),fHistDCAxy(0),fHistDCAz(0),fHistnclus(0),fHistchi2ndf(0),fHistchi2ndfvscs(0),fHistVz(0),fHistMultV0A(0),fHistMultV0C(0),fHistMultV0total(0),My_ntuple(0),fOutputList(0),fCentralityEstimator("V0M"),fCentralityBins20(kFALSE),fCentralityCounter(0),fEventCounter(0),histcounter(0){
for(Int_t ibin=0;ibin<91;ibin++)
{
fMult[ibin]=NULL;