Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);
Bool_t fDebug; //! enable debug mode
- AliAODEvent *fAOD; //AOD object
- TList *fOutputList; //output list
- AliHFEcuts *fCuts; //Cut Collection
- Bool_t fIdentifiedAsOutInz; //Out Of Range in z
- Bool_t fPassTheEventCut; //Pass The Event Cut
- AliCFManager *fCFM; //!Correction Framework Manager
- AliHFEpid *fPID; //PID
- AliHFEpidQAmanager *fPIDqa; //! PID QA manager
+ AliAODEvent *fAOD; //AOD object
+ TList *fOutputList;//output list
+ AliHFEcuts *fCuts; //Cut Collection
+ Bool_t fIdentifiedAsOutInz;//Out Of Range in z
+ Bool_t fPassTheEventCut;//Pass The Event Cut
+ AliCFManager *fCFM;//!Correction Framework Manager
+ AliHFEpid *fPID;//PID
+ AliHFEpidQAmanager *fPIDqa; //! PID QA manager
AliFlowTrackCuts *fCutsRP; // track cuts for reference particles
AliFlowTrackCuts *fNullCuts; // dummy cuts for flow event tracks
- AliFlowEvent *fFlowEvent; //! flow events Inclusive e
+ AliFlowEvent *fFlowEvent; // flow events Inclusive e
const char *fkCentralityMethod; // method used to determine centrality (V0 by default)
Double_t fCentrality; // event centrality for QA
Double_t fCentralityMin; // lower bound of cenrality bin
Double_t fCentralityMax; // upper bound of centrality bin
- Double_t fInvmassCut; //invariant mass cut value
- Double_t fpTCutmin; //pt cut value min
- Double_t fpTCutmax; //pt cut value max
- Int_t fTrigger; //invariant mass cut value
+ Double_t fInvmassCut;//invariant mass cut value
+ Double_t fpTCutmin;//pt cut value min
+ Double_t fpTCutmax; //pt cut value max
+ Int_t fTrigger; //invariant mass cut value
TH1F *fPhi; //! QA plot of azimuthal distribution of tracks used for event plane estimation
TH1F *fEta; //! QA plot of eta distribution of tracks used for event plane estimation
TH1F *fVZEROA; //! QA plot vzeroa multiplicity (all tracks in event)
TH1F *fVZEROC; //! QA plot vzeroc multiplicity (all tracks in event)
TH1F *fTPCM; //! QA plot TPC multiplicity (tracks used for event plane estimation)
- TH1F *fNoEvents; //no of events
- TH1F *fInclusiveElecPt; // Inclusive elec pt
- TH2F *fTPCnsigma; //TPC n sigma vs p
- TH2F *fTPCnsigmaAft; //TPC n sigma vs p after HFE pid
- TH2F *fITSnsigmaAft; //ITS n sigma vs p after HFE pid
- TH2F *fTPCnsigmaVSptAft; //lodviow
- TH2F *fTOFns; //track TOFnSigma
- TH2F *fTOFnsAft; //track TOFnSigma after PID
- TH2F *fTOFBetaAft;
+ TH1F *fNoEvents; //!no of events
+ TH1F *fInclusiveElecPt; //!Inclusive elec pt
+ TH2F *fTPCnsigma;//!TPC n sigma vs p
+ TH2F *fTPCnsigmaAft;//!TPC n sigma vs p after HFE pid
+ TH2F *fITSnsigmaAft;//!ITS n sigma vs p after HFE pid
+ TH2F *fTPCnsigmaVSptAft;//!lodviow
+ TH2F *fTOFns;//!track TOFnSigma
+ TH2F *fTOFnsAft;//!track TOFnSigma after PID
+ TH2F *fTOFBetaAft;//!track TOFnSigma after PID2
TH1F *fCentralityPass; // ! QA histogram of events that pass centrality cut
TH1F *fCentralityNoPass; //! QA histogram of events that do not pass centrality cut
- TH1F *fInvmassLS1; //LS Invmass for all rec par
- TH1F *fInvmassULS1;//ULS Invmass for all rec par
- TH1F *fPhotoElecPt; //photonic elec pt
- TH1F *fSemiInclElecPt; //Semi inclusive ele pt
- TH1F *fULSElecPt; //ULS elec Pt
- TH1F *fLSElecPt;// LS elec pt
+ TH1F *fInvmassLS1; //!LS Invmass for all rec par
+ TH1F *fInvmassULS1;//!ULS Invmass for all rec par
+ TH1F *fPhotoElecPt;//!photonic elec pt
+ TH1F *fSemiInclElecPt;//!Semi inclusive ele pt
+ TH1F *fULSElecPt; //! ULS elec Pt
+ TH1F *fLSElecPt;//! LS elec pt
TH2F *fMultCorAfterCuts; //! QA profile global and tpc multiplicity after outlier cut
TH2F *fMultvsCentr; //! QA profile of centralty vs multiplicity
TProfile *fSubEventDPhiv2;
- TH1D *EPVzA;//ep v0a
- TH1D *EPVzC;//ep v0c
- TH1D *EPTPC;//ep tpc
+ TH1D *EPVzA;//!ep v0a
+ TH1D *EPVzC;//!ep v0c
+ TH1D *EPTPC;//!ep tpc
THnSparseF *fV2Phi;//! v2 analysis of EP-V0
- TH1D *fvertex;//vertex
+ TH1D *fvertex;//!vertex
TH2F *fMultCorBeforeCuts; //! QA profile global and tpc multiplicity after outlier cut
THnSparseF *fQAPid; //! v2 analysis of EP-V0
Double_t fminITSnsigmaLowpT; //ID cuts its
Int_t fVz; //vertex range
Double_t fOpeningAngleCut; //openingAngle cut value
Bool_t fOP_angle; //to shitch on and off the op_angle cut
- TH1F *fOpeningAngleLS; //opening angle for LS pairs
- TH1F *fOpeningAngleULS; //opening angle for ULS pairs
+ TH1F *fOpeningAngleLS; //!opening angle for LS pairs
+ TH1F *fOpeningAngleULS; //!opening angle for ULS pairs
AliSelectNonHFE *fNonHFE;//new elienos stuff
Bool_t fDCA;//selection PHelectron
- TH2F *fITSnsigma; //TPC n sigma vs p
- TH2F *fTPCnsigmaAftITSTOF; //TPC n sigma vs p
+ TH2F *fITSnsigma;//!TPC n sigma vs p
+ TH2F *fTPCnsigmaAftITSTOF; //!TPC n sigma vs p
TH2F *fTPCnsigmaAftTOF; //TPC n sigma vs p
- TH2F *fITSnsigmaAftTOF;//jbd
- TH2F *fITSvsTOF; //TPC n sigma vs p
- TH2F *fTPCvsITS; //TPC n sigma vs p
- TH2F *fTPCvsITSafterTOF; //TPC n sigma vs p
- TH2F *fTPCvsTOF;//jbd
+ TH2F *fITSnsigmaAftTOF;//!jbd
+ TH2F *fITSvsTOF;//!TPC n sigma vs p
+ TH2F *fTPCvsITS;//!TPC n sigma vs p
+ TH2F *fTPCvsITSafterTOF;//!TPC n sigma vs p
+ TH2F *fTPCvsTOF;//!jbd
Bool_t fMultCut;//for mult correlationcut
TH2F *fMultCorAfterCentrBeforeCuts; //! QA profile global and tpc multiplicity after outlier cut
TH2F *fMultCorAfterVZTRKComp;//!after cent comp
TH1F *fCentralityAfterVZTRK;//!cent chneck2
TH1F *fCentralityAfterCorrCut;//!cent chneck2
TH2F *fMultCorAfterCorrCut;//!after cent comp
- TH1D *EPVz;//v0cep
- TH1D *EPTPCp;//tpcep
- TH1D *EPTPCn;//!tpcep
+ TH1D *EPVz;//!v0cep
+ TH1D *EPTPCp;//!tpcep
+ TH1D *EPTPCn;//!tpcen
TProfile *fSubEventDPhiv2new;
THnSparseF *fV2Phivzerotot;//! v2 analysis of EP-V0
- TH1F *fHistCentrDistr;// isto for Centr Flat
- TH1F *fCentralityNoPassForFlattening; //! QA histogram of events that do not pass centrality cut for flattening
+ TH1F *fHistCentrDistr;//-> isto for Centr Flat
+ TH1F *fCentralityNoPassForFlattening; //!QA histogram of events that do not pass centrality cut for flattening
AliAnalysisTaskFlowITSTPCTOFQCSP(const AliAnalysisTaskFlowITSTPCTOFQCSP&); // not implemented
AliFlowTrack *sTrackCont = new AliFlowTrack();
sTrackCont->Set(track);
sTrackCont->SetID(track->GetID());
- sTrackCont->SetForRPSelection(kFALSE);
+ sTrackCont->SetForRPSelection(kTRUE);
sTrackCont->SetForPOISelection(kTRUE);
sTrackCont->SetMass(2637);
for(int iRPs=0; iRPs!=fFlowEventCont->NumberOfTracks(); ++iRPs)
if(fDebug) printf(" was in RP set");
// cout << sTrack->GetID() <<" == " << iRP->GetID() << " was in RP set" <<endl;
iRPCont->SetForRPSelection(kFALSE);
- fFlowEventCont->SetNumberOfRPs(fFlowEventCont->GetNumberOfRPs() - 1);
+ // fFlowEventCont->SetNumberOfRPs(fFlowEventCont->GetNumberOfRPs() - 1);
}
} //end of for loop on RPs
fFlowEventCont->InsertTrack(((AliFlowTrack*) sTrackCont));
+ fFlowEventCont->SetNumberOfPOIs(fFlowEventCont->GetNumberOfPOIs()+1);
+
}
}
//==========================================================================================================
AliFlowTrack *sTrack = new AliFlowTrack();
sTrack->Set(track);
sTrack->SetID(track->GetID());
- sTrack->SetForRPSelection(kFALSE);
+ sTrack->SetForRPSelection(kTRUE);
sTrack->SetForPOISelection(kTRUE);
sTrack->SetMass(263732);
for(int iRPs=0; iRPs!=fFlowEvent->NumberOfTracks(); ++iRPs)
if(fDebug) printf(" was in RP set");
// cout << sTrack->GetID() <<" == " << iRP->GetID() << " was in RP set" <<endl;
iRP->SetForRPSelection(kFALSE);
- fFlowEvent->SetNumberOfRPs(fFlowEvent->GetNumberOfRPs() - 1);
+ // fFlowEvent->SetNumberOfRPs(fFlowEvent->GetNumberOfRPs() - 1);
}
} //end of for loop on RPs
fFlowEvent->InsertTrack(((AliFlowTrack*) sTrack));
-
+ fFlowEvent->SetNumberOfPOIs(fFlowEvent->GetNumberOfPOIs()+1);
if(fDCA){
AliHFEpidQAmanager *fPIDqa; //! PID QA manager
AliFlowTrackCuts *fCutsRP; // track cuts for reference particles
AliFlowTrackCuts *fNullCuts; // dummy cuts for flow event tracks
- AliFlowEvent *fFlowEvent; //! flow events Inclusive e
+ AliFlowEvent *fFlowEvent; // flow events Inclusive e
const char *fkCentralityMethod; // method used to determine centrality (V0 by default)
Double_t fCentrality; // event centrality for QA
Double_t fCentralityMin; // lower bound of cenrality bin
TH1F *fVZEROA; //! QA plot vzeroa multiplicity (all tracks in event)
TH1F *fVZEROC; //! QA plot vzeroc multiplicity (all tracks in event)
TH1F *fTPCM; //! QA plot TPC multiplicity (tracks used for event plane estimation)
- TH1F *fNoEvents; //no of events
- TH2F *fTrkEovPBef; //track E/p before HFE pid
+ TH1F *fNoEvents; //!no of events
+ TH2F *fTrkEovPBef; //!track E/p before HFE pid
// TH2F *fdEdxBef; //track dEdx vs p before HFE pid
- TH1F *fInclusiveElecPt; // Inclusive elec pt
- TH2F *fTPCnsigma; //TPC n sigma vs p
- TH2F *fTPCnsigmaAft; //TPC n sigma vs p after HFE pid
+ TH1F *fInclusiveElecPt; //! Inclusive elec pt
+ TH2F *fTPCnsigma; //!TPC n sigma vs p
+ TH2F *fTPCnsigmaAft; //!TPC n sigma vs p after HFE pid
TH1F *fCentralityPass; // ! QA histogram of events that pass centrality cut
TH1F *fCentralityNoPass; //! QA histogram of events that do not pass centrality cut
- TH1F *fInvmassLS1; //LS Invmass for all rec par
- TH1F *fInvmassULS1;//ULS Invmass for all rec par
- TH1F *fPhotoElecPt; //photonic elec pt
- TH1F *fSemiInclElecPt; //Semi inclusive ele pt
- TH1F *fULSElecPt; //ULS elec Pt
- TH1F *fLSElecPt;// LS elec pt
+ TH1F *fInvmassLS1; //!LS Invmass for all rec par
+ TH1F *fInvmassULS1;//!ULS Invmass for all rec par
+ TH1F *fPhotoElecPt; //!photonic elec pt
+ TH1F *fSemiInclElecPt; //!Semi inclusive ele pt
+ TH1F *fULSElecPt; //!ULS elec Pt
+ TH1F *fLSElecPt;//! LS elec pt
Double_t fminTPC; //ID cuts tpc
Double_t fmaxTPC; //ID cuts tpc
Double_t fminEovP; //ID cuts eovp
Double_t fDispersion;//ID cuts SS
TH2F *fMultCorAfterCuts; //! QA profile global and tpc multiplicity after outlier cut
TH2F *fMultvsCentr; //! QA profile of centralty vs multiplicity
- TProfile *fSubEventDPhiv2;
- TH1D *EPVzA;//v0aep
- TH1D *EPVzC;//v0cep
- TH1D *EPTPC;//tpcep
+ TProfile *fSubEventDPhiv2;//!lll
+ TH1D *EPVzA;//!v0aep
+ TH1D *EPVzC;//!v0cep
+ TH1D *EPTPC;//!tpcep
THnSparseF *fV2Phi;//! v2 analysis of EP-V0
THnSparseD *fSparseElectronHadron;//! Trk matching sparse for v1 clusterizer
- TH1D *fvertex;//huge ThNsparse
+ TH1D *fvertex;//!huge ThNsparse
TH2F *fMultCorBeforeCuts; //! QA profile global and tpc multiplicity after outlier cut
Bool_t fSideBandsFlow;//flow from side bands for contamination
Bool_t fPhiminusPsi;//flow from phi minus psi method
- AliFlowEvent *fFlowEventCont; //! flow events for elect Contamination
+ AliFlowEvent *fFlowEventCont; // flow events for elect Contamination
Bool_t fpurity; //for purity evaluation
THnSparseD *fSparseElectronpurity;//! Trk matching sparse for v1 clusterizer
- TH1F *fOpeningAngleLS; //opening angle for LS pairs
- TH1F *fOpeningAngleULS; //opening angle for ULS pairs
+ TH1F *fOpeningAngleLS; //!opening angle for LS pairs
+ TH1F *fOpeningAngleULS; //!opening angle for ULS pairs
AliSelectNonHFE *fNonHFE;//new elienos stuff
Bool_t fDCA;//selection PHelectron
Double_t fOpeningAngleCut; //openingAngle cut value
TH1F *fCentralityAfterVZTRK;//!cent chneck2
TH1F *fCentralityAfterCorrCut;//!cent chneck2
TH2F *fMultCorAfterCorrCut;//!after cent comp
- TH1D *EPVz;//v0cep
- TH1D *EPTPCp;//tpcep
+ TH1D *EPVz;//!v0cep
+ TH1D *EPTPCp;//!tpcep
TH1D *EPTPCn;//!tpcep
TProfile *fSubEventDPhiv2new;
THnSparseF *fV2Phivzerotot;//! v2 analysis of EP-V0
- TH1F *fHistCentrDistr;// isto for Centr Flat
+ TH1F *fHistCentrDistr;//-> isto for Centr Flat
TH1F *fCentralityNoPassForFlattening; //! QA histogram of events that do not pass centrality cut for flattening
- TH1F *fInvmassLS1highpt; //LS Invmass for all rec par high pt
- TH1F *fInvmassULS1highpt;//ULS Invmass for all rec par high pt
+ TH1F *fInvmassLS1highpt; //!LS Invmass for all rec par high pt
+ TH1F *fInvmassULS1highpt;//!ULS Invmass for all rec par high pt