const char* AliConversionMesonCuts::fgkCutNames[AliConversionMesonCuts::kNCuts] = {
-"MesonKind",
-"BackgroundScheme",
-"NumberOfBGEvents",
-"DegreesForRotationMethod",
-"RapidityMesonCut",
-"RCut",
-"AlphaMesonCut",
-"Chi2MesonCut",
-"SharedElectronCuts",
-"RejectToCloseV0s",
-"UseMCPSmearing",
+ "MesonKind",
+ "BackgroundScheme",
+ "NumberOfBGEvents",
+ "DegreesForRotationMethod",
+ "RapidityMesonCut",
+ "RCut",
+ "AlphaMesonCut",
+ "Chi2MesonCut",
+ "SharedElectronCuts",
+ "RejectToCloseV0s",
+ "UseMCPSmearing",
};
//________________________________________________________________________
-AliConversionMesonCuts::AliConversionMesonCuts(const char *name,const char *title) : AliAnalysisCuts(name,title),
- fHistograms(NULL),
- fMaxR(200),
- fChi2CutMeson(1000),
- fAlphaMinCutMeson(0),
- fAlphaCutMeson(1),
- fRapidityCutMeson(1),
- fUseRotationMethodInBG(kFALSE),
- fdoBGProbability(kFALSE),
- fUseTrackMultiplicityForBG(kFALSE),
- fnDegreeRotationPMForBG(0),
- fNumberOfBGEvents(0),
- fOpeningAngle(0.005),
- fDoToCloseV0sCut(kFALSE),
- fminV0Dist(200.),
- fDoSharedElecCut(kFALSE),
- fUseMCPSmearing(kFALSE),
- fPBremSmearing(0),
- fPSigSmearing(0),
- fPSigSmearingCte(0),
- fBrem(NULL),
- fRandom(0),
- fElectronLabelArray(NULL),
- fCutString(NULL),
- hMesonCuts(NULL),
- hMesonBGCuts(NULL)
+AliConversionMesonCuts::AliConversionMesonCuts(const char *name,const char *title) :
+ AliAnalysisCuts(name,title),
+ fHistograms(NULL),
+ fMesonKind(0),
+ fMaxR(200),
+ fChi2CutMeson(1000),
+ fAlphaMinCutMeson(0),
+ fAlphaCutMeson(1),
+ fRapidityCutMeson(1),
+ fUseRotationMethodInBG(kFALSE),
+ fDoBG(kTRUE),
+ fdoBGProbability(kFALSE),
+ fUseTrackMultiplicityForBG(kFALSE),
+ fnDegreeRotationPMForBG(0),
+ fNumberOfBGEvents(0),
+ fOpeningAngle(0.005),
+ fDoToCloseV0sCut(kFALSE),
+ fminV0Dist(200.),
+ fDoSharedElecCut(kFALSE),
+ fUseMCPSmearing(kFALSE),
+ fPBremSmearing(0),
+ fPSigSmearing(0),
+ fPSigSmearingCte(0),
+ fBrem(NULL),
+ fRandom(0),
+ fElectronLabelArray(NULL),
+ fBackgroundHandler(0),
+ fCutString(NULL),
+ hMesonCuts(NULL),
+ hMesonBGCuts(NULL)
+
{
- for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=0;}
- fCutString=new TObjString((GetCutNumber()).Data());
- fElectronLabelArray = new Int_t[500];
- if (fBrem == NULL){
- fBrem = new TF1("fBrem","pow(-log(x),[0]/log(2.0)-1.0)/TMath::Gamma([0]/log(2.0))",0.00001,0.999999999);
- // tests done with 1.0e-14
- fBrem->SetParameter(0,fPBremSmearing);
- fBrem->SetNpx(100000);
- }
+ for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=0;}
+ fCutString=new TObjString((GetCutNumber()).Data());
+ fElectronLabelArray = new Int_t[500];
+ if (fBrem == NULL){
+ fBrem = new TF1("fBrem","pow(-log(x),[0]/log(2.0)-1.0)/TMath::Gamma([0]/log(2.0))",0.00001,0.999999999);
+ // tests done with 1.0e-14
+ fBrem->SetParameter(0,fPBremSmearing);
+ fBrem->SetNpx(100000);
+ }
}
//________________________________________________________________________
AliConversionMesonCuts::~AliConversionMesonCuts() {
- // Destructor
- //Deleting fHistograms leads to seg fault it it's added to output collection of a task
- // if(fHistograms)
- // delete fHistograms;
- // fHistograms = NULL;
+ // Destructor
+ //Deleting fHistograms leads to seg fault it it's added to output collection of a task
+ // if(fHistograms)
+ // delete fHistograms;
+ // fHistograms = NULL;
if(fCutString != NULL){
delete fCutString;
fCutString = NULL;
}
//________________________________________________________________________
-void AliConversionMesonCuts::InitCutHistograms(TString name, Bool_t preCut){
-
- // Initialize Cut Histograms for QA (only initialized and filled if function is called)
-
- if(fHistograms != NULL){
- delete fHistograms;
- fHistograms=NULL;
- }
- if(fHistograms==NULL){
- fHistograms=new TList();
- if(name=="")fHistograms->SetName(Form("ConvMesonCuts_%s",GetCutNumber().Data()));
- else fHistograms->SetName(Form("%s_%s",name.Data(),GetCutNumber().Data()));
- }
-
- // Meson Cuts
- hMesonCuts=new TH1F(Form("MesonCuts %s",GetCutNumber().Data()),"MesonCuts",10,-0.5,9.5);
- hMesonCuts->GetXaxis()->SetBinLabel(1,"in");
- hMesonCuts->GetXaxis()->SetBinLabel(2,"undef rapidity");
- hMesonCuts->GetXaxis()->SetBinLabel(3,"rapidity cut");
- hMesonCuts->GetXaxis()->SetBinLabel(4,"opening angle");
- hMesonCuts->GetXaxis()->SetBinLabel(5,"alpha max");
- hMesonCuts->GetXaxis()->SetBinLabel(6,"alpha min");
- hMesonCuts->GetXaxis()->SetBinLabel(7,"out");
- fHistograms->Add(hMesonCuts);
-
- hMesonBGCuts=new TH1F(Form("MesonBGCuts %s",GetCutNumber().Data()),"MesonBGCuts",10,-0.5,9.5);
- hMesonBGCuts->GetXaxis()->SetBinLabel(1,"in");
- hMesonBGCuts->GetXaxis()->SetBinLabel(2,"undef rapidity");
- hMesonBGCuts->GetXaxis()->SetBinLabel(3,"rapidity cut");
- hMesonBGCuts->GetXaxis()->SetBinLabel(4,"opening angle");
- hMesonBGCuts->GetXaxis()->SetBinLabel(5,"alpha max");
- hMesonBGCuts->GetXaxis()->SetBinLabel(6,"alpha min");
- hMesonBGCuts->GetXaxis()->SetBinLabel(7,"out");
- fHistograms->Add(hMesonBGCuts);
+void AliConversionMesonCuts::InitCutHistograms(TString name){
+
+ // Initialize Cut Histograms for QA (only initialized and filled if function is called)
+
+ if(fHistograms != NULL){
+ delete fHistograms;
+ fHistograms=NULL;
+ }
+ if(fHistograms==NULL){
+ fHistograms=new TList();
+ if(name=="")fHistograms->SetName(Form("ConvMesonCuts_%s",GetCutNumber().Data()));
+ else fHistograms->SetName(Form("%s_%s",name.Data(),GetCutNumber().Data()));
+ }
+
+ // Meson Cuts
+ hMesonCuts=new TH1F(Form("MesonCuts %s",GetCutNumber().Data()),"MesonCuts",10,-0.5,9.5);
+ hMesonCuts->GetXaxis()->SetBinLabel(1,"in");
+ hMesonCuts->GetXaxis()->SetBinLabel(2,"undef rapidity");
+ hMesonCuts->GetXaxis()->SetBinLabel(3,"rapidity cut");
+ hMesonCuts->GetXaxis()->SetBinLabel(4,"opening angle");
+ hMesonCuts->GetXaxis()->SetBinLabel(5,"alpha max");
+ hMesonCuts->GetXaxis()->SetBinLabel(6,"alpha min");
+ hMesonCuts->GetXaxis()->SetBinLabel(7,"out");
+ fHistograms->Add(hMesonCuts);
+
+ hMesonBGCuts=new TH1F(Form("MesonBGCuts %s",GetCutNumber().Data()),"MesonBGCuts",10,-0.5,9.5);
+ hMesonBGCuts->GetXaxis()->SetBinLabel(1,"in");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(2,"undef rapidity");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(3,"rapidity cut");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(4,"opening angle");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(5,"alpha max");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(6,"alpha min");
+ hMesonBGCuts->GetXaxis()->SetBinLabel(7,"out");
+ fHistograms->Add(hMesonBGCuts);
}
//________________________________________________________________________
-Bool_t AliConversionMesonCuts::MesonIsSelectedMC(TParticle *fMCMother,AliStack *fMCStack,Bool_t bMCDaughtersInAcceptance){
- // Returns true for all pions within acceptance cuts for decay into 2 photons
- // If bMCDaughtersInAcceptance is selected, it requires in addition that both daughter photons are within acceptance cuts
-
- if(!fMCStack)return kFALSE;
-
- if(fMCMother->GetPdgCode()==111 || fMCMother->GetPdgCode()==221){
- if(fMCMother->R()>fMaxR) return kFALSE; // cuts on distance from collision point
-
- Double_t rapidity = 10.;
- if(fMCMother->Energy() - fMCMother->Pz() == 0 || fMCMother->Energy() + fMCMother->Pz() == 0){
- rapidity=8.;
- } else{
- rapidity = 0.5*(TMath::Log((fMCMother->Energy()+fMCMother->Pz()) / (fMCMother->Energy()-fMCMother->Pz())));
- }
-
- // Rapidity Cut
- if(TMath::Abs(rapidity)>fRapidityCutMeson)return kFALSE;
-
- // Select only -> 2y decay channel
- if(fMCMother->GetNDaughters()!=2)return kFALSE;
-
- for(Int_t i=0;i<2;i++){
- TParticle *MDaughter=fMCStack->Particle(fMCMother->GetDaughter(i));
- // Is Daughter a Photon?
- if(MDaughter->GetPdgCode()!=22)return kFALSE;
- // Is Photon in Acceptance?
- // if(bMCDaughtersInAcceptance){
- // if(!PhotonIsSelectedMC(MDaughter,fMCStack)){return kFALSE;}
- // }
- }
- return kTRUE;
- }
- return kFALSE;
+Bool_t AliConversionMesonCuts::MesonIsSelectedMC(TParticle *fMCMother,AliStack *fMCStack){
+ // Returns true for all pions within acceptance cuts for decay into 2 photons
+ // If bMCDaughtersInAcceptance is selected, it requires in addition that both daughter photons are within acceptance cuts
+
+ if(!fMCStack)return kFALSE;
+
+ if(fMCMother->GetPdgCode()==111 || fMCMother->GetPdgCode()==221){
+ if(fMCMother->R()>fMaxR) return kFALSE; // cuts on distance from collision point
+
+ Double_t rapidity = 10.;
+ if(fMCMother->Energy() - fMCMother->Pz() == 0 || fMCMother->Energy() + fMCMother->Pz() == 0){
+ rapidity=8.;
+ } else{
+ rapidity = 0.5*(TMath::Log((fMCMother->Energy()+fMCMother->Pz()) / (fMCMother->Energy()-fMCMother->Pz())));
+ }
+
+ // Rapidity Cut
+ if(abs(rapidity)>fRapidityCutMeson)return kFALSE;
+
+ // Select only -> 2y decay channel
+ if(fMCMother->GetNDaughters()!=2)return kFALSE;
+
+ for(Int_t i=0;i<2;i++){
+ TParticle *MDaughter=fMCStack->Particle(fMCMother->GetDaughter(i));
+ // Is Daughter a Photon?
+ if(MDaughter->GetPdgCode()!=22)return kFALSE;
+ // Is Photon in Acceptance?
+ // if(bMCDaughtersInAcceptance){
+ // if(!PhotonIsSelectedMC(MDaughter,fMCStack)){return kFALSE;}
+ // }
+ }
+ return kTRUE;
+ }
+ return kFALSE;
}
//________________________________________________________________________
Bool_t AliConversionMesonCuts::MesonIsSelectedMCDalitz(TParticle *fMCMother,AliStack *fMCStack){
- // Returns true for all pions within acceptance cuts for decay into 2 photons
- // If bMCDaughtersInAcceptance is selected, it requires in addition that both daughter photons are within acceptance cuts
+ // Returns true for all pions within acceptance cuts for decay into 2 photons
+ // If bMCDaughtersInAcceptance is selected, it requires in addition that both daughter photons are within acceptance cuts
- if(!fMCStack)return kFALSE;
+ if(!fMCStack)return kFALSE;
- if(fMCMother->GetPdgCode()==111 || fMCMother->GetPdgCode()==221){
+ if(fMCMother->GetPdgCode()==111 || fMCMother->GetPdgCode()==221){
- if(fMCMother->R()>fMaxR) return kFALSE; // cuts on distance from collision point
-
- Double_t rapidity = 10.;
- if(fMCMother->Energy() - fMCMother->Pz() == 0 || fMCMother->Energy() + fMCMother->Pz() == 0){
- rapidity=8.;
- }
- else{
- rapidity = 0.5*(TMath::Log((fMCMother->Energy()+fMCMother->Pz()) / (fMCMother->Energy()-fMCMother->Pz())));
- }
+ if(fMCMother->R()>fMaxR) return kFALSE; // cuts on distance from collision point
+
+ Double_t rapidity = 10.;
+ if(fMCMother->Energy() - fMCMother->Pz() == 0 || fMCMother->Energy() + fMCMother->Pz() == 0){
+ rapidity=8.;
+ }
+ else{
+ rapidity = 0.5*(TMath::Log((fMCMother->Energy()+fMCMother->Pz()) / (fMCMother->Energy()-fMCMother->Pz())));
+ }
- // Rapidity Cut
- if(TMath::Abs(rapidity)>fRapidityCutMeson)return kFALSE;
-
- // Select only -> Dalitz decay channel
- if(fMCMother->GetNDaughters()!=3)return kFALSE;
-
- Int_t daughterPDGs[3] = {0,0,0};
- Int_t index = 0;
-
-// iParticle->GetFirstDaughter(); idxPi0 <= iParticle->GetLastDaughter()
-
- for(Int_t i=fMCMother->GetFirstDaughter(); i<= fMCMother->GetLastDaughter();i++){
- TParticle *MDaughter=fMCStack->Particle(i);
- // Is Daughter a Photon or an electron?
- daughterPDGs[index]=MDaughter->GetPdgCode();
- index++;
- }
- for (Int_t j=0;j<2;j++){
-
- for (Int_t i=0;i<2;i++){
- if (daughterPDGs[i] > daughterPDGs[i+1]){
- Int_t interMed = daughterPDGs[i] ;
- daughterPDGs[i] = daughterPDGs[i+1];
- daughterPDGs[i+1] = interMed;
- }
- }
- }
- if (daughterPDGs[0] == -11 && daughterPDGs[1] == 11 && daughterPDGs[2] == 22) return kTRUE;
- }
- return kFALSE;
+ // Rapidity Cut
+ if(abs(rapidity)>fRapidityCutMeson)return kFALSE;
+
+ // Select only -> Dalitz decay channel
+ if(fMCMother->GetNDaughters()!=3)return kFALSE;
+
+ Int_t daughterPDGs[3] = {0,0,0};
+ Int_t index = 0;
+
+ // iParticle->GetFirstDaughter(); idxPi0 <= iParticle->GetLastDaughter()
+
+ for(Int_t i=fMCMother->GetFirstDaughter(); i<= fMCMother->GetLastDaughter();i++){
+ TParticle *MDaughter=fMCStack->Particle(i);
+ // Is Daughter a Photon or an electron?
+ daughterPDGs[index]=MDaughter->GetPdgCode();
+ index++;
+ }
+ for (Int_t j=0;j<2;j++){
+
+ for (Int_t i=0;i<2;i++){
+ if (daughterPDGs[i] > daughterPDGs[i+1]){
+ Int_t interMed = daughterPDGs[i] ;
+ daughterPDGs[i] = daughterPDGs[i+1];
+ daughterPDGs[i+1] = interMed;
+ }
+ }
+ }
+ if (daughterPDGs[0] == -11 && daughterPDGs[1] == 11 && daughterPDGs[2] == 22) return kTRUE;
+ }
+ return kFALSE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::MesonIsSelected(AliAODConversionMother *pi0,Bool_t IsSignal)
{
- // Selection of reconstructed Meson candidates
- // Use flag IsSignal in order to fill Fill different
- // histograms for Signal and Background
- TH1 *hist=0x0;
-
-
-
- if(IsSignal){hist=hMesonCuts;}
- else{hist=hMesonBGCuts;}
-
- Int_t cutIndex=0;
-
- if(hist)hist->Fill(cutIndex);
- cutIndex++;
-
- // Undefined Rapidity -> Floating Point exception
- if((pi0->E()+pi0->Pz())/(pi0->E()-pi0->Pz())<=0){
- if(hist)hist->Fill(cutIndex);
- cutIndex++;
- return kFALSE;
- }
- else{
- // PseudoRapidity Cut --> But we cut on Rapidity !!!
- cutIndex++;
- if(TMath::Abs(pi0->Rapidity())>fRapidityCutMeson){
- //if(TMath::Abs(pi0->PseudoRapidity())>fRapidityCutMeson){
- if(hist)hist->Fill(cutIndex);
- return kFALSE;
- }
- }
- cutIndex++;
-
- // Opening Angle Cut
- //fOpeningAngle=2*TMath::ATan(0.134/pi0->P());// physical minimum opening angle
- if( pi0->GetOpeningAngle() < fOpeningAngle){
- if(hist)hist->Fill(cutIndex);
- return kFALSE;
- }
- cutIndex++;
-
- // Alpha Max Cut
- if(pi0->GetAlpha()>fAlphaCutMeson){
- if(hist)hist->Fill(cutIndex);
- return kFALSE;
- }
- cutIndex++;
-
- // Alpha Min Cut
- if(pi0->GetAlpha()<fAlphaMinCutMeson){
- if(hist)hist->Fill(cutIndex);
- return kFALSE;
- }
- cutIndex++;
-
- if(hist)hist->Fill(cutIndex);
- return kTRUE;
+ // Selection of reconstructed Meson candidates
+ // Use flag IsSignal in order to fill Fill different
+ // histograms for Signal and Background
+ TH1 *hist=0x0;
+
+ if(IsSignal){hist=hMesonCuts;}
+ else{hist=hMesonBGCuts;}
+
+ Int_t cutIndex=0;
+
+ if(hist)hist->Fill(cutIndex);
+ cutIndex++;
+
+ // Undefined Rapidity -> Floating Point exception
+ if((pi0->E()+pi0->Pz())/(pi0->E()-pi0->Pz())<=0){
+ if(hist)hist->Fill(cutIndex);
+ cutIndex++;
+ return kFALSE;
+ }
+ else{
+ // PseudoRapidity Cut --> But we cut on Rapidity !!!
+ cutIndex++;
+ if(abs(pi0->Rapidity())>fRapidityCutMeson){
+ if(hist)hist->Fill(cutIndex);
+ return kFALSE;
+ }
+ }
+ cutIndex++;
+
+ // Opening Angle Cut
+ //fOpeningAngle=2*TMath::ATan(0.134/pi0->P());// physical minimum opening angle
+ if( pi0->GetOpeningAngle() < fOpeningAngle){
+ if(hist)hist->Fill(cutIndex);
+ return kFALSE;
+ }
+ cutIndex++;
+
+ // Alpha Max Cut
+ if(pi0->GetAlpha()>fAlphaCutMeson){
+ if(hist)hist->Fill(cutIndex);
+ return kFALSE;
+ }
+ cutIndex++;
+
+ // Alpha Min Cut
+ if(pi0->GetAlpha()<fAlphaMinCutMeson){
+ if(hist)hist->Fill(cutIndex);
+ return kFALSE;
+ }
+ cutIndex++;
+
+ if(hist)hist->Fill(cutIndex);
+ return kTRUE;
}
///________________________________________________________________________
///________________________________________________________________________
-Bool_t AliConversionMesonCuts::UpdateCutString(cutIds cutID, Int_t value) {
-///Update the cut string (if it has been created yet)
-
- if(fCutString && fCutString->GetString().Length() == kNCuts) {
-// cout << "Updating cut id in spot number " << cutID << " to " << value << endl;
- fCutString->SetString(GetCutNumber());
- } else {
-// cout << "fCutString not yet initialized, will not be updated" << endl;
- return kFALSE;
- }
-// cout << fCutString->GetString().Data() << endl;
- return kTRUE;
+Bool_t AliConversionMesonCuts::UpdateCutString() {
+ ///Update the cut string (if it has been created yet)
+
+ if(fCutString && fCutString->GetString().Length() == kNCuts) {
+ // cout << "Updating cut id in spot number " << cutID << " to " << value << endl;
+ fCutString->SetString(GetCutNumber());
+ } else {
+ // cout << "fCutString not yet initialized, will not be updated" << endl;
+ return kFALSE;
+ }
+ // cout << fCutString->GetString().Data() << endl;
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::InitializeCutsFromCutString(const TString analysisCutSelection ) {
- // Initialize Cuts from a given Cut string
- AliInfo(Form("Set Meson Cutnumber: %s",analysisCutSelection.Data()));
- if(analysisCutSelection.Length()!=kNCuts) {
- AliError(Form("Cut selection has the wrong length! size is %d, number of cuts is %d", analysisCutSelection.Length(), kNCuts));
- return kFALSE;
- }
- if(!analysisCutSelection.IsDigit()){
- AliError("Cut selection contains characters");
- return kFALSE;
- }
+ // Initialize Cuts from a given Cut string
+ AliInfo(Form("Set Meson Cutnumber: %s",analysisCutSelection.Data()));
+ if(analysisCutSelection.Length()!=kNCuts) {
+ AliError(Form("Cut selection has the wrong length! size is %d, number of cuts is %d", analysisCutSelection.Length(), kNCuts));
+ return kFALSE;
+ }
+ if(!analysisCutSelection.IsDigit()){
+ AliError("Cut selection contains characters");
+ return kFALSE;
+ }
- const char *cutSelection = analysisCutSelection.Data();
- #define ASSIGNARRAY(i) fCuts[i] = cutSelection[i] - '0'
- for(Int_t ii=0;ii<kNCuts;ii++){
- ASSIGNARRAY(ii);
- }
-
- // Set Individual Cuts
- for(Int_t ii=0;ii<kNCuts;ii++){
+ const char *cutSelection = analysisCutSelection.Data();
+#define ASSIGNARRAY(i) fCuts[i] = cutSelection[i] - '0'
+ for(Int_t ii=0;ii<kNCuts;ii++){
+ ASSIGNARRAY(ii);
+ }
+
+ // Set Individual Cuts
+ for(Int_t ii=0;ii<kNCuts;ii++){
if(!SetCut(cutIds(ii),fCuts[ii]))return kFALSE;
- }
+ }
- //PrintCuts();
- return kTRUE;
+ //PrintCuts();
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetCut(cutIds cutID, const Int_t value) {
- ///Set individual cut ID
-
- //cout << "Updating cut " << fgkCutNames[cutID] << " (" << cutID << ") to " << value << endl;
- switch (cutID) {
- case kMesonKind:
- if( SetMesonKind(value)) {
- fCuts[kMesonKind] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
-
- case kchi2MesonCut:
- if( SetChi2MesonCut(value)) {
- fCuts[kchi2MesonCut] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kalphaMesonCut:
- if( SetAlphaMesonCut(value)) {
- fCuts[kalphaMesonCut] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kRCut:
- if( SetRCut(value)) {
- fCuts[kRCut] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kRapidityMesonCut:
- if( SetRapidityMesonCut(value)) {
- fCuts[kRapidityMesonCut] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kBackgroundScheme:
- if( SetBackgroundScheme(value)) {
- fCuts[kBackgroundScheme] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kDegreesForRotationMethod:
- if( SetNDegreesForRotationMethod(value)) {
- fCuts[kDegreesForRotationMethod] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kNumberOfBGEvents:
- if( SetNumberOfBGEvents(value)) {
- fCuts[kNumberOfBGEvents] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kuseMCPSmearing:
- if( SetMCPSmearing(value)) {
- fCuts[kuseMCPSmearing] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
-
-
- case kElecShare:
- if( SetSharedElectronCut(value)) {
- fCuts[kElecShare] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kToCloseV0s:
- if( SetToCloseV0sCut(value)) {
- fCuts[kToCloseV0s] = value;
- UpdateCutString(cutID, value);
- return kTRUE;
- } else return kFALSE;
-
- case kNCuts:
- cout << "Error:: Cut id out of range"<< endl;
- return kFALSE;
- }
+ ///Set individual cut ID
+
+ //cout << "Updating cut " << fgkCutNames[cutID] << " (" << cutID << ") to " << value << endl;
+ switch (cutID) {
+ case kMesonKind:
+ if( SetMesonKind(value)) {
+ fCuts[kMesonKind] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+ case kchi2MesonCut:
+ if( SetChi2MesonCut(value)) {
+ fCuts[kchi2MesonCut] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+ case kalphaMesonCut:
+ if( SetAlphaMesonCut(value)) {
+ fCuts[kalphaMesonCut] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+ case kRCut:
+ if( SetRCut(value)) {
+ fCuts[kRCut] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kRapidityMesonCut:
+ if( SetRapidityMesonCut(value)) {
+ fCuts[kRapidityMesonCut] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kBackgroundScheme:
+ if( SetBackgroundScheme(value)) {
+ fCuts[kBackgroundScheme] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kDegreesForRotationMethod:
+ if( SetNDegreesForRotationMethod(value)) {
+ fCuts[kDegreesForRotationMethod] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kNumberOfBGEvents:
+ if( SetNumberOfBGEvents(value)) {
+ fCuts[kNumberOfBGEvents] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kuseMCPSmearing:
+ if( SetMCPSmearing(value)) {
+ fCuts[kuseMCPSmearing] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+ case kElecShare:
+ if( SetSharedElectronCut(value)) {
+ fCuts[kElecShare] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+ case kToCloseV0s:
+ if( SetToCloseV0sCut(value)) {
+ fCuts[kToCloseV0s] = value;
+ UpdateCutString();
+ return kTRUE;
+ } else return kFALSE;
+
+ case kNCuts:
+ cout << "Error:: Cut id out of range"<< endl;
+ return kFALSE;
+ }
- cout << "Error:: Cut id " << cutID << " not recognized "<< endl;
- return kFALSE;
+ cout << "Error:: Cut id " << cutID << " not recognized "<< endl;
+ return kFALSE;
}
///________________________________________________________________________
void AliConversionMesonCuts::PrintCuts() {
- // Print out current Cut Selection
- for(Int_t ic = 0; ic < kNCuts; ic++) {
- printf("%-30s : %d \n", fgkCutNames[ic], fCuts[ic]);
- }
+ // Print out current Cut Selection
+ for(Int_t ic = 0; ic < kNCuts; ic++) {
+ printf("%-30s : %d \n", fgkCutNames[ic], fCuts[ic]);
+ }
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetMesonKind(Int_t mesonKind){
- // Set Cut
- switch(mesonKind){
- case 0:
- fMesonKind=0;
- break;
- case 1:
- fMesonKind=1;;
- break;
- default:
- cout<<"Warning: Meson kind not defined"<<mesonKind<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(mesonKind){
+ case 0:
+ fMesonKind=0;
+ break;
+ case 1:
+ fMesonKind=1;;
+ break;
+ default:
+ cout<<"Warning: Meson kind not defined"<<mesonKind<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetRCut(Int_t rCut){
- // Set Cut
- switch(rCut){
- case 0:
- fMaxR = 180.;
- break;
- case 1:
- fMaxR = 180.;
- break;
- case 2:
- fMaxR = 180.;
- break;
- case 3:
- fMaxR = 70.;
- break;
- case 4:
- fMaxR = 70.;
- break;
- case 5:
- fMaxR = 180.;
- break;
- // High purity cuts for PbPb
- case 9:
- fMaxR = 180.;
- break;
- default:
- cout<<"Warning: rCut not defined"<<rCut<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(rCut){
+ case 0:
+ fMaxR = 180.;
+ break;
+ case 1:
+ fMaxR = 180.;
+ break;
+ case 2:
+ fMaxR = 180.;
+ break;
+ case 3:
+ fMaxR = 70.;
+ break;
+ case 4:
+ fMaxR = 70.;
+ break;
+ case 5:
+ fMaxR = 180.;
+ break;
+ // High purity cuts for PbPb
+ case 9:
+ fMaxR = 180.;
+ break;
+ default:
+ cout<<"Warning: rCut not defined"<<rCut<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetChi2MesonCut(Int_t chi2MesonCut){
- // Set Cut
- switch(chi2MesonCut){
- case 0: // 100.
- fChi2CutMeson = 100.;
- break;
- case 1: // 50.
- fChi2CutMeson = 50.;
- break;
- case 2: // 30.
- fChi2CutMeson = 30.;
- break;
- case 3:
- fChi2CutMeson = 200.;
- break;
- case 4:
- fChi2CutMeson = 500.;
- break;
- case 5:
- fChi2CutMeson = 1000.;
- break;
- default:
- cout<<"Warning: Chi2MesonCut not defined "<<chi2MesonCut<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(chi2MesonCut){
+ case 0: // 100.
+ fChi2CutMeson = 100.;
+ break;
+ case 1: // 50.
+ fChi2CutMeson = 50.;
+ break;
+ case 2: // 30.
+ fChi2CutMeson = 30.;
+ break;
+ case 3:
+ fChi2CutMeson = 200.;
+ break;
+ case 4:
+ fChi2CutMeson = 500.;
+ break;
+ case 5:
+ fChi2CutMeson = 1000.;
+ break;
+ default:
+ cout<<"Warning: Chi2MesonCut not defined "<<chi2MesonCut<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetAlphaMesonCut(Int_t alphaMesonCut)
{ // Set Cut
- switch(alphaMesonCut){
- case 0: // 0- 0.7
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.7;
- break;
- case 1: // 0-0.5
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.5;
- break;
- case 2: // 0.5-1
- fAlphaMinCutMeson = 0.5;
- fAlphaCutMeson = 1.;
- break;
- case 3: // 0.0-1
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 1.;
- break;
- case 4: // 0-0.65
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.65;
- break;
- case 5: // 0-0.75
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.75;
- break;
- case 6: // 0-0.8
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.8;
- break;
- case 7: // 0.0-0.85
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.85;
- break;
- case 8: // 0.0-0.6
- fAlphaMinCutMeson = 0.0;
- fAlphaCutMeson = 0.6;
- break;
- default:
- cout<<"Warning: AlphaMesonCut not defined "<<alphaMesonCut<<endl;
- return kFALSE;
- }
- return kTRUE;
+ switch(alphaMesonCut){
+ case 0: // 0- 0.7
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.7;
+ break;
+ case 1: // 0-0.5
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.5;
+ break;
+ case 2: // 0.5-1
+ fAlphaMinCutMeson = 0.5;
+ fAlphaCutMeson = 1.;
+ break;
+ case 3: // 0.0-1
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 1.;
+ break;
+ case 4: // 0-0.65
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.65;
+ break;
+ case 5: // 0-0.75
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.75;
+ break;
+ case 6: // 0-0.8
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.8;
+ break;
+ case 7: // 0.0-0.85
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.85;
+ break;
+ case 8: // 0.0-0.6
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.6;
+ break;
+ case 9: // 0.0-0.3
+ fAlphaMinCutMeson = 0.0;
+ fAlphaCutMeson = 0.3;
+ break;
+ default:
+ cout<<"Warning: AlphaMesonCut not defined "<<alphaMesonCut<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetRapidityMesonCut(Int_t RapidityMesonCut){
- // Set Cut
- switch(RapidityMesonCut){
- case 0: //
- fRapidityCutMeson = 0.9;
- break;
- case 1: //
- fRapidityCutMeson = 0.8;
- break;
- case 2: //
- fRapidityCutMeson = 0.7;
- break;
- case 3: //
- fRapidityCutMeson = 0.6;
- break;
- case 4: //
- fRapidityCutMeson = 0.5;
- break;
-
- default:
- cout<<"Warning: RapidityMesonCut not defined "<<RapidityMesonCut<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(RapidityMesonCut){
+ case 0: //
+ fRapidityCutMeson = 0.9;
+ break;
+ case 1: //
+ fRapidityCutMeson = 0.8;
+ break;
+ case 2: //
+ fRapidityCutMeson = 0.7;
+ break;
+ case 3: //
+ fRapidityCutMeson = 0.6;
+ break;
+ case 4: //
+ fRapidityCutMeson = 0.5;
+ break;
+ case 5: //
+ fRapidityCutMeson = 0.85;
+ break;
+ case 6: //
+ fRapidityCutMeson = 0.75;
+ break;
+
+ default:
+ cout<<"Warning: RapidityMesonCut not defined "<<RapidityMesonCut<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetBackgroundScheme(Int_t BackgroundScheme){
- // Set Cut
- switch(BackgroundScheme){
- case 0: //Rotation
- fUseRotationMethodInBG=kTRUE;
- fdoBGProbability=kFALSE;
- break;
- case 1: // mixed event with V0 multiplicity
- fUseRotationMethodInBG=kFALSE;
- fUseTrackMultiplicityForBG=kFALSE;
- fdoBGProbability=kFALSE;
- break;
- case 2: // mixed event with track multiplicity
- fUseRotationMethodInBG=kFALSE;
- fUseTrackMultiplicityForBG=kTRUE;
- fdoBGProbability=kFALSE;
- break;
- case 3: //Rotation
- fUseRotationMethodInBG=kTRUE;
- fdoBGProbability=kTRUE;
- break;
- default:
- cout<<"Warning: BackgroundScheme not defined "<<BackgroundScheme<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(BackgroundScheme){
+ case 0: //Rotation
+ fUseRotationMethodInBG=kTRUE;
+ fdoBGProbability=kFALSE;
+ break;
+ case 1: // mixed event with V0 multiplicity
+ fUseRotationMethodInBG=kFALSE;
+ fUseTrackMultiplicityForBG=kFALSE;
+ fdoBGProbability=kFALSE;
+ break;
+ case 2: // mixed event with track multiplicity
+ fUseRotationMethodInBG=kFALSE;
+ fUseTrackMultiplicityForBG=kTRUE;
+ fdoBGProbability=kFALSE;
+ break;
+ case 3: //Rotation
+ fUseRotationMethodInBG=kTRUE;
+ fdoBGProbability=kTRUE;
+ break;
+ case 4: //No BG calculation
+ cout << "no BG calculation should be done" << endl;
+ fUseRotationMethodInBG=kFALSE;
+ fdoBGProbability=kFALSE;
+ fDoBG=kFALSE;
+ break;
+ case 5: //Rotation
+ fUseRotationMethodInBG=kTRUE;
+ fdoBGProbability=kFALSE;
+ fBackgroundHandler = 1;
+ break;
+ case 6: // mixed event with V0 multiplicity
+ fUseRotationMethodInBG=kFALSE;
+ fUseTrackMultiplicityForBG=kFALSE;
+ fdoBGProbability=kFALSE;
+ fBackgroundHandler = 1;
+ break;
+ case 7: // mixed event with track multiplicity
+ fUseRotationMethodInBG=kFALSE;
+ fUseTrackMultiplicityForBG=kTRUE;
+ fdoBGProbability=kFALSE;
+ fBackgroundHandler = 1;
+ break;
+ case 8: //Rotation
+ fUseRotationMethodInBG=kTRUE;
+ fdoBGProbability=kTRUE;
+ fBackgroundHandler = 1;
+ break;
+ default:
+ cout<<"Warning: BackgroundScheme not defined "<<BackgroundScheme<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetNDegreesForRotationMethod(Int_t DegreesForRotationMethod){
- // Set Cut
- switch(DegreesForRotationMethod){
- case 0:
- fnDegreeRotationPMForBG = 5;
- break;
- case 1:
- fnDegreeRotationPMForBG = 10;
- break;
- case 2:
- fnDegreeRotationPMForBG = 15;
- break;
- case 3:
- fnDegreeRotationPMForBG = 20;
- break;
- default:
- cout<<"Warning: DegreesForRotationMethod not defined "<<DegreesForRotationMethod<<endl;
- return kFALSE;
- }
- fCuts[kDegreesForRotationMethod]=DegreesForRotationMethod;
- return kTRUE;
+ // Set Cut
+ switch(DegreesForRotationMethod){
+ case 0:
+ fnDegreeRotationPMForBG = 5;
+ break;
+ case 1:
+ fnDegreeRotationPMForBG = 10;
+ break;
+ case 2:
+ fnDegreeRotationPMForBG = 15;
+ break;
+ case 3:
+ fnDegreeRotationPMForBG = 20;
+ break;
+ default:
+ cout<<"Warning: DegreesForRotationMethod not defined "<<DegreesForRotationMethod<<endl;
+ return kFALSE;
+ }
+ fCuts[kDegreesForRotationMethod]=DegreesForRotationMethod;
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetNumberOfBGEvents(Int_t NumberOfBGEvents){
- // Set Cut
- switch(NumberOfBGEvents){
- case 0:
- fNumberOfBGEvents = 5;
- break;
- case 1:
- fNumberOfBGEvents = 10;
- break;
- case 2:
- fNumberOfBGEvents = 15;
- break;
- case 3:
- fNumberOfBGEvents = 20;
- break;
- case 4:
- fNumberOfBGEvents = 2;
- break;
- case 5:
- fNumberOfBGEvents = 50;
- break;
- case 6:
- fNumberOfBGEvents = 80;
- break;
- case 7:
- fNumberOfBGEvents = 100;
- break;
- default:
- cout<<"Warning: NumberOfBGEvents not defined "<<NumberOfBGEvents<<endl;
- return kFALSE;
- }
- return kTRUE;
+ // Set Cut
+ switch(NumberOfBGEvents){
+ case 0:
+ fNumberOfBGEvents = 5;
+ break;
+ case 1:
+ fNumberOfBGEvents = 10;
+ break;
+ case 2:
+ fNumberOfBGEvents = 15;
+ break;
+ case 3:
+ fNumberOfBGEvents = 20;
+ break;
+ case 4:
+ fNumberOfBGEvents = 2;
+ break;
+ case 5:
+ fNumberOfBGEvents = 50;
+ break;
+ case 6:
+ fNumberOfBGEvents = 80;
+ break;
+ case 7:
+ fNumberOfBGEvents = 100;
+ break;
+ default:
+ cout<<"Warning: NumberOfBGEvents not defined "<<NumberOfBGEvents<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
-
-// ///________________________________________________________________________
-// Bool_t AliConversionMesonCuts::SetPsiPairCut(Int_t psiCut) {
-//
-// switch(psiCut) {
-// case 0:
-// fPsiPairCut = 10000; //
-// break;
-// case 1:
-// fPsiPairCut = 0.1; //
-// break;
-// case 2:
-// fPsiPairCut = 0.05; // Standard
-// break;
-// case 3:
-// fPsiPairCut = 0.035; //
-// break;
-// case 4:
-// fPsiPairCut = 0.15; //
-// break;
-// case 5:
-// fPsiPairCut = 0.2; //
-// break;
-// case 6:
-// fPsiPairCut = 0.03; //
-// break;
-// case 7:
-// fPsiPairCut = 0.025; //
-// break;
-// case 8:
-// fPsiPairCut = 0.01; //
-// break;
-// default:
-// cout<<"Warning: PsiPairCut not defined "<<fPsiPairCut<<endl;
-// return kFALSE;
-// }
-//
-// return kTRUE;
-// }
-
-
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetSharedElectronCut(Int_t sharedElec) {
- switch(sharedElec){
- case 0:
- fDoSharedElecCut = kFALSE;
- break;
- case 1:
- fDoSharedElecCut = kTRUE;
- break;
- default:
- cout<<"Warning: Shared Electron Cut not defined "<<sharedElec<<endl;
- return kFALSE;
- }
+ switch(sharedElec){
+ case 0:
+ fDoSharedElecCut = kFALSE;
+ break;
+ case 1:
+ fDoSharedElecCut = kTRUE;
+ break;
+ default:
+ cout<<"Warning: Shared Electron Cut not defined "<<sharedElec<<endl;
+ return kFALSE;
+ }
- return kTRUE;
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetToCloseV0sCut(Int_t toClose) {
- switch(toClose){
- case 0:
- fDoToCloseV0sCut = kFALSE;
- fminV0Dist = 250;
- break;
- case 1:
- fDoToCloseV0sCut = kTRUE;
- fminV0Dist = 1;
- break;
- case 2:
- fDoToCloseV0sCut = kTRUE;
- fminV0Dist = 2;
- break;
- case 3:
- fDoToCloseV0sCut = kTRUE;
- fminV0Dist = 3;
- break;
- default:
- cout<<"Warning: Shared Electron Cut not defined "<<toClose<<endl;
- return kFALSE;
- }
- return kTRUE;
+ switch(toClose){
+ case 0:
+ fDoToCloseV0sCut = kFALSE;
+ fminV0Dist = 250;
+ break;
+ case 1:
+ fDoToCloseV0sCut = kTRUE;
+ fminV0Dist = 1;
+ break;
+ case 2:
+ fDoToCloseV0sCut = kTRUE;
+ fminV0Dist = 2;
+ break;
+ case 3:
+ fDoToCloseV0sCut = kTRUE;
+ fminV0Dist = 3;
+ break;
+ default:
+ cout<<"Warning: Shared Electron Cut not defined "<<toClose<<endl;
+ return kFALSE;
+ }
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::SetMCPSmearing(Int_t useMCPSmearing)
{// Set Cut
- switch(useMCPSmearing){
- case 0:
- fUseMCPSmearing=0;
- fPBremSmearing=1.;
- fPSigSmearing=0.;
- fPSigSmearingCte=0.;
- break;
- case 1:
- fUseMCPSmearing=1;
- fPBremSmearing=1.0e-14;
- fPSigSmearing=0.;
- fPSigSmearingCte=0.;
- break;
- case 2:
- fUseMCPSmearing=1;
- fPBremSmearing=1.0e-15;
- fPSigSmearing=0.0;
- fPSigSmearingCte=0.;
- break;
- case 3:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.003;
- fPSigSmearingCte=0.002;
- break;
- case 4:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.003;
- fPSigSmearingCte=0.007;
- break;
- case 5:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.003;
- fPSigSmearingCte=0.016;
- break;
- case 6:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.007;
- fPSigSmearingCte=0.016;
- break;
- case 7:
- fUseMCPSmearing=1;
- fPBremSmearing=1.0e-16;
- fPSigSmearing=0.0;
- fPSigSmearingCte=0.;
- break;
- case 8:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.007;
- fPSigSmearingCte=0.014;
- break;
- case 9:
- fUseMCPSmearing=1;
- fPBremSmearing=1.;
- fPSigSmearing=0.007;
- fPSigSmearingCte=0.011;
- break;
-
- default:
- AliError("Warning: UseMCPSmearing not defined");
- return kFALSE;
- }
- return kTRUE;
+ switch(useMCPSmearing){
+ case 0:
+ fUseMCPSmearing=0;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.;
+ fPSigSmearingCte=0.;
+ break;
+ case 1:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.0e-14;
+ fPSigSmearing=0.;
+ fPSigSmearingCte=0.;
+ break;
+ case 2:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.0e-15;
+ fPSigSmearing=0.0;
+ fPSigSmearingCte=0.;
+ break;
+ case 3:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.003;
+ fPSigSmearingCte=0.002;
+ break;
+ case 4:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.003;
+ fPSigSmearingCte=0.007;
+ break;
+ case 5:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.003;
+ fPSigSmearingCte=0.016;
+ break;
+ case 6:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.007;
+ fPSigSmearingCte=0.016;
+ break;
+ case 7:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.0e-16;
+ fPSigSmearing=0.0;
+ fPSigSmearingCte=0.;
+ break;
+ case 8:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.007;
+ fPSigSmearingCte=0.014;
+ break;
+ case 9:
+ fUseMCPSmearing=1;
+ fPBremSmearing=1.;
+ fPSigSmearing=0.007;
+ fPSigSmearingCte=0.011;
+ break;
+
+ default:
+ AliError("Warning: UseMCPSmearing not defined");
+ return kFALSE;
+ }
+ return kTRUE;
}
-
-// Bool_t AliConversionMesonCuts::PsiPairCut(const AliConversionPhotonBase * photon) const {
-//
-// if(photon->GetPsiPair() > fPsiPairCut){
-// return kFALSE;}
-// else{return kTRUE;}
-//
-// }
-
///________________________________________________________________________
TString AliConversionMesonCuts::GetCutNumber(){
- // returns TString with current cut number
- TString a(kNCuts);
- for(Int_t ii=0;ii<kNCuts;ii++){
- a.Append(Form("%d",fCuts[ii]));
- }
- return a;
+ // returns TString with current cut number
+ TString a(kNCuts);
+ for(Int_t ii=0;ii<kNCuts;ii++){
+ a.Append(Form("%d",fCuts[ii]));
+ }
+ return a;
}
///________________________________________________________________________
void AliConversionMesonCuts::FillElectonLabelArray(AliAODConversionPhoton* photon, Int_t nV0){
- Int_t posLabel = photon->GetTrackLabelPositive();
- Int_t negLabel = photon->GetTrackLabelNegative();
+ Int_t posLabel = photon->GetTrackLabelPositive();
+ Int_t negLabel = photon->GetTrackLabelNegative();
- fElectronLabelArray[nV0*2] = posLabel;
- fElectronLabelArray[(nV0*2)+1] = negLabel;
+ fElectronLabelArray[nV0*2] = posLabel;
+ fElectronLabelArray[(nV0*2)+1] = negLabel;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::RejectSharedElectronV0s(AliAODConversionPhoton* photon, Int_t nV0, Int_t nV0s){
- Int_t posLabel = photon->GetTrackLabelPositive();
- Int_t negLabel = photon->GetTrackLabelNegative();
+ Int_t posLabel = photon->GetTrackLabelPositive();
+ Int_t negLabel = photon->GetTrackLabelNegative();
- for(Int_t i = 0; i<nV0s*2;i++){
- if(i==nV0*2) continue;
- if(i==(nV0*2)+1) continue;
- if(fElectronLabelArray[i] == posLabel){
- return kFALSE;}
- if(fElectronLabelArray[i] == negLabel){
- return kFALSE;}
- }
+ for(Int_t i = 0; i<nV0s*2;i++){
+ if(i==nV0*2) continue;
+ if(i==(nV0*2)+1) continue;
+ if(fElectronLabelArray[i] == posLabel){
+ return kFALSE;}
+ if(fElectronLabelArray[i] == negLabel){
+ return kFALSE;}
+ }
- return kTRUE;
+ return kTRUE;
}
///________________________________________________________________________
Bool_t AliConversionMesonCuts::RejectToCloseV0s(AliAODConversionPhoton* photon, TList *photons, Int_t nV0){
- Double_t posX = photon->GetConversionX();
- Double_t posY = photon->GetConversionY();
- Double_t posZ = photon->GetConversionZ();
-
- for(Int_t i = 0;i<photons->GetEntries();i++){
- if(nV0 == i) continue;
- AliAODConversionPhoton *photonComp = (AliAODConversionPhoton*) photons->At(i);
- Double_t posCompX = photonComp->GetConversionX();
- Double_t posCompY = photonComp->GetConversionY();
- Double_t posCompZ = photonComp->GetConversionZ();
+ Double_t posX = photon->GetConversionX();
+ Double_t posY = photon->GetConversionY();
+ Double_t posZ = photon->GetConversionZ();
+
+ for(Int_t i = 0;i<photons->GetEntries();i++){
+ if(nV0 == i) continue;
+ AliAODConversionPhoton *photonComp = (AliAODConversionPhoton*) photons->At(i);
+ Double_t posCompX = photonComp->GetConversionX();
+ Double_t posCompY = photonComp->GetConversionY();
+ Double_t posCompZ = photonComp->GetConversionZ();
- Double_t dist = pow((posX - posCompX),2)+pow((posY - posCompY),2)+pow((posZ - posCompZ),2);
+ Double_t dist = pow((posX - posCompX),2)+pow((posY - posCompY),2)+pow((posZ - posCompZ),2);
- if(dist < fminV0Dist*fminV0Dist){
- if(photon->GetChi2perNDF() < photonComp->GetChi2perNDF()) return kTRUE;
- else {
- return kFALSE;}
- }
+ if(dist < fminV0Dist*fminV0Dist){
+ if(photon->GetChi2perNDF() < photonComp->GetChi2perNDF()) return kTRUE;
+ else {
+ return kFALSE;}
+ }
- }
- return kTRUE;
+ }
+ return kTRUE;
}
///________________________________________________________________________