Andrea Dubla's update
authorrbailhac <Raphaelle.Bailhache@cern.ch>
Mon, 3 Feb 2014 13:05:28 +0000 (14:05 +0100)
committerrbailhac <Raphaelle.Bailhache@cern.ch>
Mon, 3 Feb 2014 13:05:28 +0000 (14:05 +0100)
PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.h
PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h

index aeafacd..bb793b7 100755 (executable)
@@ -80,55 +80,55 @@ private:
     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
@@ -147,18 +147,18 @@ private:
     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
@@ -166,13 +166,13 @@ private:
     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
index b23975b..c271e8c 100644 (file)
@@ -626,7 +626,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserExec(Option_t*)
                 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)
@@ -640,10 +640,12 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserExec(Option_t*)
                         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);
+
             }
         }
         //==========================================================================================================
@@ -673,7 +675,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserExec(Option_t*)
         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)
@@ -687,11 +689,11 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::UserExec(Option_t*)
                 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){
index 5408bd9..80da9bd 100644 (file)
@@ -98,7 +98,7 @@ private:
     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
@@ -111,20 +111,20 @@ private:
     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
@@ -136,21 +136,21 @@ private:
     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
@@ -164,15 +164,15 @@ private:
     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