Bool_t AliAnalysisTaskGammaConvDalitzV1::Notify()
{
for(Int_t iCut = 0; iCut<fnCuts;iCut++){
+
if( !((AliConvEventCuts*)fCutEventArray->At(iCut))->GetDoEtaShift() ){
- /*if (((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift() != 0.){
- printf("Error: Gamma Conversion Dalitz Task %s :: Eta Shift not requested but set to %f, reset to 00. \n\n",
- (((AliConvEventCuts*)fCutEventArray->At(iCut))->GetCutNumber()).Data(),((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift());
- ((AliConvEventCuts*)fCutEventArray->At(iCut))->SetEtaShift(0.);
- ((AliDalitzElectronCuts*)fCutElectronArray->At(iCut))->SetEtaShift(0.);
-
-
- }*/
hEtaShift[iCut]->Fill(0.,0.);
continue; // No Eta Shift requested, continue
}
if( ((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift() == 0.0){ // Eta Shift requested but not set, get shift automatically
((AliConvEventCuts*)fCutEventArray->At(iCut))->GetCorrectEtaShiftFromPeriod(fV0Reader->GetPeriodName());
((AliConvEventCuts*)fCutEventArray->At(iCut))->DoEtaShift(kFALSE); // Eta Shift Set, make sure that it is called only once
- ((AliDalitzElectronCuts*)fCutElectronArray->At(iCut))->SetEtaShift( ((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift() );
- hEtaShift[iCut]->Fill(0.,(((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift()));
+ hEtaShift[iCut]->Fill(0.,(((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift()));
continue;
} else {
printf(" Gamma Conversion Dalitz Task %s :: Eta Shift Manually Set to %f \n\n",
(((AliConvEventCuts*)fCutEventArray->At(iCut))->GetCutNumber()).Data(),((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift());
((AliConvEventCuts*)fCutEventArray->At(iCut))->DoEtaShift(kFALSE); // Eta Shift Set, make sure that it is called only once
- ((AliDalitzElectronCuts*)fCutElectronArray->At(iCut))->SetEtaShift( ((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift() );
hEtaShift[iCut]->Fill(0.,(((AliConvEventCuts*)fCutEventArray->At(iCut))->GetEtaShift()));
}
}
fPIDResponse(NULL),
fesdTrackCuts(NULL),
fEtaCut(0.9),
- fEtaShift(0.0),
fDoEtaCut(kFALSE),
fPtMinCut(0.0),
fPtMaxCut(9999),
if( TMath::Abs( particle->GetPdgCode() ) != 11 ) return kFALSE;
if( fDoEtaCut ){
- if( particle->Eta() > (fEtaCut + fEtaShift) || particle->Eta() < (-fEtaCut + fEtaShift) )
+ if( particle->Eta() > fEtaCut || particle->Eta() < -fEtaCut )
return kFALSE;
}
}
if( fDoEtaCut ) {
- if( lTrack->Eta() > (fEtaCut + fEtaShift) || lTrack->Eta() < (-fEtaCut + fEtaShift) ) {
+ if( lTrack->Eta() > fEtaCut || lTrack->Eta() < -fEtaCut ) {
return kFALSE;
}
}
Bool_t SetEtaCut(Int_t etaCut);
Bool_t SetPtCut(Int_t ptCut);
Bool_t SetDCACut(Int_t dcaCut);
- void SetEtaShift(Double_t etaShift){fEtaShift = etaShift;}
Bool_t SetMinMomPiondEdxCut(Int_t piMinMomdedxSigmaCut);
Bool_t SetMaxMomPiondEdxCut(Int_t piMaxMomdedxSigmaCut);
Bool_t SetLowPRejectionCuts(Int_t LowPRejectionSigmaCut);
AliESDtrackCuts *fesdTrackCuts;
Double_t fEtaCut; //eta cutç
- Double_t fEtaShift;
Bool_t fDoEtaCut;
Double_t fPtMinCut;
Double_t fPtMaxCut;