Add multiplicity histogram for filter-bit-4 tracks
authorzconesa <zaida.conesa.del.valle@cern.ch>
Thu, 6 Mar 2014 15:34:43 +0000 (00:34 +0900)
committerzconesa <zaida.conesa.del.valle@cern.ch>
Thu, 6 Mar 2014 15:35:17 +0000 (00:35 +0900)
PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx

index 40a95da..a444ce0 100644 (file)
@@ -818,6 +818,8 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
 
       hname="hMult";
       TH1F* hMult=new TH1F(hname.Data(),"Multiplicity;multiplicity;Entries",10000,-0.5,9999.5);
+      hname="hMultFBit4";
+      TH1F* hMultFBit4=new TH1F(hname.Data(),"Multiplicity (global+tracklet) with filter bit 4;multiplicity;Entries",10000,-0.5,9999.5);
       hname="hMultComb05";
       TH1F* hMultC05=new TH1F(hname.Data(),"Multiplicity (global+tracklet) in |#eta|<0.5;multiplicity;Entries",10000,-0.5,9999.5);
       hname="hMultComb08";
@@ -825,6 +827,7 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
 
       fOutputTrack->Add(hNtracklets);
       fOutputTrack->Add(hMult);
+      fOutputTrack->Add(hMultFBit4);
       fOutputTrack->Add(hMultC05);
       fOutputTrack->Add(hMultC08);
     }
@@ -1800,6 +1803,7 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
 
   if(fOnOff[0] || fOnOff[1]){
     //loop on tracks in the event
+    Int_t ntracksFBit4=0;
     for (Int_t k=0;k<ntracks;k++){
       AliAODTrack* track=aod->GetTrack(k);
 
@@ -1809,6 +1813,7 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
       if(!track->PropagateToDCA(vtx1,aod->GetMagneticField(),99999.,d0z0,covd0z0)) continue;
       if(track->TestFilterMask(AliAODTrack::kTrkGlobalNoDCA)){
        ((TH1F*)fOutputTrack->FindObject("hd0TracksFilterBit4"))->Fill(d0z0[0]);
+       ntracksFBit4++;
       }
       ULong_t trStatus=track->GetStatus();
       if(trStatus&AliESDtrack::kITSrefit){
@@ -2072,6 +2077,7 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
       if (fReadMC) ((TH1F*)fOutputTrack->FindObject("hdistrFakeTr"))->Fill(isFakeTrack);
       ((TH1F*)fOutputTrack->FindObject("hdistrGoodTr"))->Fill(isGoodTrack);
       ((TH1F*)fOutputTrack->FindObject("hdistrSelTr"))->Fill(isSelTrack);
+      ((TH1F*)fOutputTrack->FindObject("hMultFBit4"))->Fill(ntracksFBit4);
     }
 
     if(!isSimpleMode){