* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-const int c_array_size = 23;
+const int c_array_size = 29; // RRnewTOF c_array_size was increased to 27 in order to put a TOF PID cut on electrons
class AliAnalysisDataContainer;
class AliGammaConversionHistograms;
//Svein
Bool_t kGCRunGammaJetTask = kFALSE;
/** ---------------------------------- define cuts here ------------------------------------*/
-TString kGCAnalysisCutSelectionId="90035620401003321022000"; // do not change here, use -set-cut-selection in argument instead
+TString kGCAnalysisCutSelectionId="90035620401003321022000000090"; // do not change here, use -set-cut-selection in argument instead
Int_t kGCNEventsForBGCalculation=20;
Bool_t kGCUseMultiplicity = kFALSE;
Int_t kGCUseMultiplicityBin=0;
+Int_t kGCIsHeavyIon = 0;
+Int_t kGCUseCentrality = 0;
+Int_t kGCUseCentralityBin = 0;
+Int_t kGCUseCorrectedTPCClsInfo = 0;
+Int_t kGCUseMCPSmearing=0;
+
+Double_t kGCPBremSmearing=1.;
+Double_t kGCPSigSmearing=0.;
+Double_t kGCPSigSmearingCte=0.;
+
/** ---------------------------------- define pi0 dalitz cuts here ------------------------------------*/
Bool_t kGCRunStandalone = kTRUE;
Bool_t kGCplotMCPi0Energy = kFALSE;
Bool_t kGCplotMCPi0Mass = kTRUE;
Bool_t kGCplotMCPi0Alpha = kTRUE;
+Bool_t kGCplotMCEtaAlpha = kTRUE;
Bool_t kGCplotMCPi0OpeningAngle = kTRUE;
Bool_t kGCplotMCPi0R = kTRUE;
Bool_t kGCplotMCPi0ZR = kFALSE;
Bool_t kGCplotESDENTPCClusters = kTRUE;
Bool_t kGCplotESDENITSClusters = kTRUE;
Bool_t kGCplotESDENTPCClustersToFP = kTRUE;
+Bool_t kGCplotESDENTPCClustersToFR = kTRUE;
Bool_t kGCplotESDETPCchi2 = kTRUE;
Bool_t kGCplotESDPEnergy = kFALSE;
Bool_t kGCplotESDPNTPCClusters = kTRUE;
Bool_t kGCplotESDPNITSClusters = kTRUE;
Bool_t kGCplotESDPNTPCClustersToFP = kTRUE;
+Bool_t kGCplotESDPNTPCClustersToFR = kTRUE;
Bool_t kGCplotESDPTPCchi2 = kTRUE;
Bool_t kGCplotESDConvGammaEnergy = kFALSE;
Bool_t kGCplotESDConvGammaPtvsEta = kTRUE;
Bool_t kGCplotESDConvGammaPtvsChi2 = kTRUE;
Bool_t kGCplotESDConvGammaEtavsChi2 = kTRUE;
+Bool_t kGCplotESDConvGammaPtvsTOF = kTRUE;
Bool_t kGCplotESDTrueDalitzContaminationR = kTRUE;
+Bool_t kGCplotESDTruePi0DalitzContaminationR = kTRUE;
+Bool_t kGCplotESDTrueEtaDalitzContaminationR = kTRUE;
+Bool_t kGCplotESDTrueCombinatorialContaminationR = kTRUE;
+Bool_t kGCplotESDTrueCombinatorialElecContaminationR = kTRUE;
+Bool_t kGCplotESDTrueHadronicContaminationR = kTRUE;
+Bool_t kGCplotESDTrueCombinatorialContaminationPt = kTRUE;
+Bool_t kGCplotESDTrueCombinatorialElecContaminationPt = kTRUE;
+Bool_t kGCplotESDTrueHadronicContaminationPt = kTRUE;
+
+
+Bool_t kGCplotESDTrueBackground = kTRUE;
Bool_t kGCplotESDTrueConvGammaEnergy = kFALSE;
Bool_t kGCplotESDTrueConvGammaPt = kTRUE;
Bool_t kGCplotESDTrueConvGammaEta = kTRUE;
Bool_t kGCplotESDTrueConversionMCZR = kFALSE;
Bool_t kGCplotESDTrueConversionMCXY = kFALSE;
+Bool_t kGCplotESDNoCutAllV0Pt = kTRUE;
Bool_t kGCplotESDNoCutConvGammaEnergy = kFALSE;
Bool_t kGCplotESDNoCutConvGammaPt = kTRUE;
Bool_t kGCplotESDNoCutConvGammaEta = kTRUE;
Bool_t kGCplotESDCutZ = kTRUE;
Bool_t kGCplotESDCutMinClsTPC = kTRUE;
Bool_t kGCplotESDCutMinClsTPCToF = kTRUE;
+Bool_t kGCplotESDCutPhotonAsymmetry = kTRUE;
Bool_t kGCplotESDGoodV0s = kTRUE;
Bool_t kGCplotESDAllV0s = kTRUE;
Bool_t kGCplotESDAllV0sCurrentFinder = kTRUE;
//EventQuality-plot
-Int_t kGCnXBinsEvtQ= 8;
+Int_t kGCnXBinsEvtQ= 9;
Double_t kGCfirstXBinEvtQ=-1.5;
-Double_t kGClastXBinEvtQ=6.5;
+Double_t kGClastXBinEvtQ=7.5;
//R-plots
Int_t kGCnXBinsR = 400;
Double_t kGCfirstXBinPt = 0.;
Double_t kGClastXBinPt = 50.;
+//TOF-plots RRnewTOF start /////////////////////////////////////////
+Int_t kGCnXBinsTOFsignal = 500;
+Double_t kGCfirstXBinTOFsignal = -10000.;
+Double_t kGClastXBinTOFsignal = 40000.;
+// end RRnewTOF ////////////////////////////////////////////////////
+
//Eta-plots
Int_t kGCnXBinsEta = 40;
Double_t kGCfirstXBinEta = -2.;
Double_t kGClastXBinNITSClusters = 6.5;
//TPCcluster to Findable-plots
-Int_t kGCnYBinsClsToF = 120;
+Int_t kGCnYBinsClsToF = 200;
Double_t kGCfirstYBinClsToF = 0.;
-Double_t kGClastYBinClsToF = 1.2;
+Double_t kGClastYBinClsToF = 2.0;
//TPCchi2 -plots
Int_t kGCnXBinsTPCchi2 = 100;
Double_t kGCPIDnSigmaAbovePionLine=0;
Double_t kGCPIDMinPnSigmaAbovePionLine=1.;
Double_t kGCPIDMaxPnSigmaAbovePionLine=3.;
+Double_t kGCPIDnSigmaAbovePionLineHighPt=0;
+
+Bool_t kGCuseTOFpid = kFALSE; // RRnewTOF start //////////
+Double_t kGCtofPIDnSigmaBelowElectronLine=-100;
+Double_t kGCtofPIDnSigmaAboveElectronLine=100; // RRnewTOF end
/**------- Flag to apply rejection at LowP of Kaons, protons , pions------------*/
Bool_t kGCdoPionRejectionLowP=kTRUE;
Bool_t kGCdoQtGammaSelection=kTRUE;
Double_t kGCQtMax=100.;
+Bool_t kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+Double_t kGCHighPtQtMax=100.; // RRnew
+Double_t kGCPtBorderForQt=100.; // RRnew
+/**----Flag to apply cut on the photon asymmetry -----*/
+Bool_t kGCdoPhotonAsymmetryCut= kTRUE;
+Double_t kGCMinPPhotonAsymmetryCut=100.;
+Double_t kGCMinPhotonAsymmetry=0.;
Bool_t scanArguments(TString arguments){
cout<<"Number Of files to analyze: "<<kGCnumberOfFilesToAnalyze<<endl;
build();//build (if necessary) and load the libraries needed
-
+ LoadLibraries();
gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C"); // load the CreateChain macro
}
v0Reader->SetPIDnSigmaAboveElectronLine(kGCPIDnSigmaAboveElectronLine);
v0Reader->SetPIDnSigmaBelowElectronLine(kGCPIDnSigmaBelowElectronLine);
v0Reader->SetPIDnSigmaAbovePionLine(kGCPIDnSigmaAbovePionLine);
+ v0Reader->SetPIDnSigmaAbovePionLineHighPt(kGCPIDnSigmaAbovePionLineHighPt);
v0Reader->SetPIDMinPnSigmaAbovePionLine(kGCPIDMinPnSigmaAbovePionLine);
v0Reader->SetPIDMaxPnSigmaAbovePionLine(kGCPIDMaxPnSigmaAbovePionLine);
v0Reader->SetOnFlyFlag(kGCUseOnFlyV0Finder);
v0Reader->SetCalculateBackground(kGCcalculateBackground);
+ v0Reader->SetDoTOFsigmaCut(kGCuseTOFpid); // RRnewTOF
+ v0Reader->SetTofPIDnSigmaAboveElectronLine(kGCtofPIDnSigmaAboveElectronLine); // RRnewTOF
+ v0Reader->SetTofPIDnSigmaBelowElectronLine(kGCtofPIDnSigmaBelowElectronLine); // RRnewTOF
+
// for the rejection at LowP based on Dedx signal
v0Reader->SetDoKaonRejectionLowP(kGCdoKaonRejectionLowP);
v0Reader->SetPIDMinPProtonRejectionLowP(kGCPIDMinPProtonRejectionLowP);
v0Reader->SetPIDMinPPionRejectionLowP(kGCPIDMinPPionRejectionLowP);
v0Reader->SetDoQtGammaSelection(kGCdoQtGammaSelection);
+ v0Reader->SetDoHighPtQtGammaSelection(kGCdoHighPtQtGammaSelection); // RRnew
v0Reader->SetQtMax(kGCQtMax);
+ v0Reader->SetHighPtQtMax(kGCHighPtQtMax); // RRnew
+ v0Reader->SetPtBorderForQt(kGCPtBorderForQt); // RRnew
+
+ v0Reader->SetDoPhotonAsymmetryCut(kGCdoPhotonAsymmetryCut);
+ v0Reader->SetMinPPhotonAsymmetryCut(kGCMinPPhotonAsymmetryCut);
+ v0Reader->SetMinPhotonAsymmetry(kGCMinPhotonAsymmetry);
+
+
+ kGCNEventsForBGCalculation= kGCnumberOfRotationEventsForBG;
+ cout<< "number of Events used for mixing::"<<kGCNEventsForBGCalculation<<endl;
v0Reader->SetNEventsForBG(kGCNEventsForBGCalculation);
// Create the GammaConversionTask
else{
cout<<"Using mixed event for bg."<<endl;
}
+ if(kGCUseTrackMultiplicityForBG){
+ cout<<"Using track multiplicity for bck"<<endl;
+ }else{
+ cout<<"Using V0 multiplicity for bck"<<endl;
+ }
gammaconversion->SetDoRotation(kGCUseRotationMethodInBG);
cout<<"Using :"<<kGCnumberOfRotationEventsForBG<<" rotations in bg calculation"<<endl;
gammaconversion->SetNumberOfRotationsBG(kGCnumberOfRotationEventsForBG);
gammaconversion->SetUseMultiplicity(kGCUseMultiplicity);
gammaconversion->SetUseMultiplicityBin(kGCUseMultiplicityBin);
+ v0Reader->SetIsHeavyIon(kGCIsHeavyIon);
+ v0Reader->SetUseCorrectedTPCClsInfo(kGCUseCorrectedTPCClsInfo);
+ gammaconversion->SetUseCentrality(kGCUseCentrality);
+ if(kGCUseCentrality){
+ gammaconversion->SetUseCentralityBin(kGCUseCentralityBin);
+ }
+ v0Reader->SetUseMCPSmearing(kGCUseMCPSmearing);
+ v0Reader->SetPBremSmearing(kGCPBremSmearing);
+ v0Reader->SetPSigSmearing(kGCPSigSmearing);
+ v0Reader->SetPSigSmearingCte(kGCPSigSmearingCte);
+
+
// for CF
gammaconversion->SetCFManager(man);
gammaconversion->SetDoCF(kGCrunCF);
// Define Output Event Handler and add
if(kGCWriteAOD){
gammaconversion->SetForceAOD(kGCForceAOD);
+ gammaconversion->SetAODBranchName(Form("GammaConv_%s", kGCAnalysisCutSelectionId.Data()));
- if( kGCrunOnTrain ) {
+ if( kGCrunOnTrain ) {
+
AliAODHandler * aodHandler = dynamic_cast<AliAODHandler*>(mgr->GetOutputEventHandler());
if(!aodHandler) {
::Error("This task requires an AOD handler");
- return -1;
+ return NULL;
}
-
+
gammaconversion->SetDeltaAODFileName(kGCDeltaAODFilename);
if(kGCDeltaAODFilename.Length() > 0) {
} else {
if(kGCDeltaAODFilename.Length() == 0 ) {
cout << "Error:: Need a file name for the AOD"<<endl;
- return;
+ return NULL;
}
AliAODHandler* aodHandler = new AliAODHandler();
aodHandler->SetOutputFileName(kGCDeltaAODFilename);
aodHandler->SetCreateNonStandardAOD();
mgr->SetOutputEventHandler(aodHandler);
}
+ } else {
+ gammaconversion->SetCreateAOD(kFALSE);
}
// Connect I/O to the task
mgr->ConnectInput (gammaconversion, 0, cinput1);
if(mgr->GetCommonOutputContainer())
mgr->ConnectOutput(gammaconversion, 0, mgr->GetCommonOutputContainer());
-
+
mgr->ConnectOutput(gammaconversion, 1, coutput2);
mgr->ConnectOutput(gammaconversion, 2, coutput3);
-
+
if(kGCRunGammaJetTask) {
AliAnalysisTaskGammaJet * gammaJetTask = new AliAnalysisTaskGammaJet("GammaJetTask");
if(kGCrunOnTrain) {
return gammaconversion;
}
+void LoadLibraries() {
+
+ TStopwatch timer;
+ timer.Start();
+ gSystem->Load("libTree.so");
+ gSystem->Load("libGeom");
+
+ gSystem->Load("libSTEERBase.so");
+ gSystem->Load("libVMC.so");
+ gSystem->Load("libESD.so");
+ gSystem->Load("libAOD.so");
+ gSystem->Load("libANALYSIS.so");
+ gSystem->Load("libANALYSISalice.so");
+ gSystem->Load("libCORRFW.so");
+ gSystem->Load("libPWG4GammaConv.so");
+
+ // gSystem->ChangeDirectory(pwd.Data());
+
+}
void build() {
TStopwatch timer;
histograms->AddHistogram("MC_DecayOmegaGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
histograms->AddHistogram("MC_DecayK0sGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
histograms->AddHistogram("MC_DecayEtapGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("MC_DecayPhiGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
histograms->AddHistogram("MC_DecayAllGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
}
if(kGCplotESDENTPCClusters == kTRUE){ histograms->AddHistogram("ESD_E_nTPCClusters" ,"" , kGCnXBinsNTPCClusters, kGCfirstXBinNTPCClusters, kGClastXBinNTPCClusters, "", "");}
if(kGCplotESDENITSClusters == kTRUE){ histograms->AddHistogram("ESD_E_nITSClusters" ,"" , kGCnXBinsNITSClusters, kGCfirstXBinNITSClusters, kGClastXBinNITSClusters, "", "");}
if(kGCplotESDENTPCClustersToFP== kTRUE){ histograms->AddHistogram("ESD_E_nTPCClustersToFP" ,"" ,kGCnXBinsP, kGCfirstXBinP, kGClastXBinP,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+ if(kGCplotESDENTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_E_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+ if(kGCplotESDENTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_TrueConversion_E_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
if(kGCplotESDETPCchi2 == kTRUE){ histograms->AddHistogram("ESD_E_TPCchi2" ,"" , kGCnXBinsTPCchi2, kGCfirstXBinTPCchi2, kGClastXBinTPCchi2, "", "");}
if(kGCplotESDPNTPCClusters == kTRUE){ histograms->AddHistogram("ESD_P_nTPCClusters" ,"" , kGCnXBinsNTPCClusters, kGCfirstXBinNTPCClusters, kGClastXBinNTPCClusters, "", "");}
if(kGCplotESDPNITSClusters == kTRUE){ histograms->AddHistogram("ESD_P_nITSClusters" ,"" , kGCnXBinsNITSClusters, kGCfirstXBinNITSClusters, kGClastXBinNITSClusters, "", "");}
if(kGCplotESDPNTPCClustersToFP== kTRUE){ histograms->AddHistogram("ESD_P_nTPCClustersToFP" ,"" ,kGCnXBinsP, kGCfirstXBinP, kGClastXBinP,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+ if(kGCplotESDPNTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_P_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+ if(kGCplotESDPNTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_TrueConversion_P_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
if(kGCplotESDPTPCchi2 == kTRUE){ histograms->AddHistogram("ESD_P_TPCchi2" ,"" , kGCnXBinsTPCchi2, kGCfirstXBinTPCchi2, kGClastXBinTPCchi2, "", "");}
if(kGCplotESDConvGammaEnergy == kTRUE){ histograms->AddHistogram("ESD_ConvGamma_Energy" ,"" , kGCnXBinsEnergy, kGCfirstXBinEnergy, kGClastXBinEnergy, "", "");}
if(kGCplotESDConvGammaPtvsEta == kTRUE){ histograms->AddHistogram("ESD_ConvGamma_Pt_Eta","", kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta,"","" );}
if(kGCplotESDConvGammaPtvsChi2 == kTRUE){ histograms->AddHistogram("ESD_ConvGamma_Pt_Chi2" ,"" ,kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, kGCnXBinsGammaChi2, kGCfirstXBinGammaChi2, kGClastXBinGammaChi2, "", "");}
if(kGCplotESDConvGammaEtavsChi2 == kTRUE){ histograms->AddHistogram("ESD_ConvGamma_Eta_Chi2" ,"" ,kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, kGCnXBinsGammaChi2, kGCfirstXBinGammaChi2, kGClastXBinGammaChi2, "", "");}
-
+ histograms->AddHistogram("ESD_ConvGamma_EandP_P_dT" ,"" , kGCnXBinsP, kGCfirstXBinP, kGClastXBinP, kGCnXBinsTOFsignal, kGCfirstXBinTOFsignal, kGClastXBinTOFsignal, "", "",0); // RRnewTOF
if(kGCplotESDConversionR == kTRUE){ histograms->AddHistogram("ESD_Conversion_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
if(kGCplotESDTrueDalitzContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueDalitzContamination_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
-
+ if(kGCplotESDTruePi0DalitzContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueConvDalitzPi0_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
+ if(kGCplotESDTrueEtaDalitzContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueConvDalitzEta_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
+ if(kGCplotESDTrueCombinatorialContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueConvCombinatorial_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
+ if(kGCplotESDTrueCombinatorialElecContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueConvCombinatorialElec_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
+ if(kGCplotESDTrueHadronicContaminationR == kTRUE){ histograms->AddHistogram("ESD_TrueConvHadronicBck_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
+ if(kGCplotESDTrueCombinatorialContaminationPt == kTRUE){ histograms->AddHistogram("ESD_TrueConvCombinatorial_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotESDTrueCombinatorialElecContaminationPt == kTRUE){ histograms->AddHistogram("ESD_TrueConvCombinatorialElec_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotESDTrueHadronicContaminationPt == kTRUE){ histograms->AddHistogram("ESD_TrueConvHadronicBck_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+
+
+ if(kGCplotESDTrueBackground){
+ histograms->AddHistogram("ESD_TrueConvCombinatorialDaughter_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsPt,kGCfirstXBinPt , kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvHadronicBckDaughter_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsPt,kGCfirstXBinPt , kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPiDaughter_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsPt,kGCfirstXBinPt , kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPiPDaughter_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsPt,kGCfirstXBinPt , kGClastXBinPt, "", "");
+
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPi_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPi_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPiP_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialPiP_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialElecPi_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvCombinatorialElecPi_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvMeson_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvMeson_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvLambda_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConvLambda_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ }
+
+ histograms->AddHistogram("ESD_TrueConvDalitzPi0_SinglePos_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvDalitzPi0_SingleNeg_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvDalitzPi0_SinglePos_kFirst_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConvDalitzPi0_SingleNeg_kFirst_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");
+ histograms->AddHistogram("ESD_TrueConversion_SinglePos_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,"", "");
+ histograms->AddHistogram("ESD_TrueConversion_SingleNeg_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConversion_SinglePos_kFirst_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+ histograms->AddHistogram("ESD_TrueConversion_SingleNeg_kFirst_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+
+
if(kGCplotESDTrueConvGammaEnergy == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Energy" ,"" , kGCnXBinsEnergy, kGCfirstXBinEnergy, kGClastXBinEnergy, "", "");}
if(kGCplotESDTrueConvGammaPt == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotESDTrueConvGammaEta == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Eta" ,"" , kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, "", "");}
if(kGCplotESDTrueConvGammaNDF == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_NDF" ,"" , kGCnXBinsGammaNDF, kGCfirstXBinGammaNDF, kGClastXBinGammaNDF, "", "");}
if(kGCplotESDTrueConvGammaRapid == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Rapid" ,"" , kGCnXBinsRapid, kGCfirstXBinRapid, kGClastXBinRapid, "", "");}
if(kGCplotESDTrueConvGammaPtvsEta == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Pt_Eta" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt,kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, "", "");}
+
if(kGCplotESDTrueConvGammaPtvsChi2 == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Pt_Chi2" ,"" ,kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, kGCnXBinsGammaChi2, kGCfirstXBinGammaChi2, kGClastXBinGammaChi2, "", "");}
if(kGCplotESDTrueConvGammaEtavsChi2 == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_Eta_Chi2" ,"" ,kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, kGCnXBinsGammaChi2, kGCfirstXBinGammaChi2, kGClastXBinGammaChi2, "", "");}
if(kGCplotESDTrueConvGammaMCPtEta == kTRUE){ histograms->AddHistogram("ESD_TrueConvGamma_MC_Pt_Eta" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, "", "");}
if(kGCplotESDTrueConversionMCZR == kTRUE){ histograms->AddHistogram("ESD_TrueConversion_MC_ZR" ,"" , kGCnXBinsZR, kGCfirstXBinZR, kGClastXBinZR, kGCnYBinsZR, kGCfirstYBinZR, kGClastYBinZR, "", "");}
if(kGCplotESDTrueConversionMCXY == kTRUE){ histograms->AddHistogram("ESD_TrueConversion_MC_XY" ,"" , kGCnXBinsXY, kGCfirstXBinXY, kGClastXBinXY, kGCnYBinsXY, kGCfirstYBinXY, kGClastYBinXY, "", "");}
-
-
-
+ histograms->AddHistogram("ESD_TrueConvGamma_EandP_P_dT" ,"" , kGCnXBinsP, kGCfirstXBinP, kGClastXBinP, kGCnXBinsTOFsignal, kGCfirstXBinTOFsignal, kGClastXBinTOFsignal, "", "",0); // RRnewTOF
+ histograms->AddHistogram("ESD_TrueConvCombinatorial_DaughtersNotElec_P_dT" ,"" , kGCnXBinsP, kGCfirstXBinP, kGClastXBinP, kGCnXBinsTOFsignal, kGCfirstXBinTOFsignal, kGClastXBinTOFsignal, "", "",0); // RRnewTOF
+ histograms->AddHistogram("ESD_TrueConvHadronicBck_Daughters_P_dT" ,"" , kGCnXBinsP, kGCfirstXBinP, kGClastXBinP, kGCnXBinsTOFsignal, kGCfirstXBinTOFsignal, kGClastXBinTOFsignal, "", "",0); // RRnewTOF
+
+ histograms->AddHistogram("ESD_NoCutConvGamma_EandP_P_dT" ,"" , kGCnXBinsP, kGCfirstXBinP, kGClastXBinP, kGCnXBinsTOFsignal, kGCfirstXBinTOFsignal, kGClastXBinTOFsignal, "", "",0); // RRnewTOF
+ if(kGCplotESDNoCutAllV0Pt == kTRUE){ histograms->AddHistogram("ESD_NoCutAllV0_Pt" ,"" , kGCnXBinsPt,kGCfirstXBinPt , kGClastXBinPt, "", "");}
if(kGCplotESDNoCutConvGammaEnergy == kTRUE){ histograms->AddHistogram("ESD_NoCutConvGamma_Energy" ,"" , kGCnXBinsEnergy, kGCfirstXBinEnergy, kGClastXBinEnergy, "", "");}
if(kGCplotESDNoCutConvGammaPt == kTRUE){ histograms->AddHistogram("ESD_NoCutConvGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotESDNoCutConvGammaEta == kTRUE){ histograms->AddHistogram("ESD_NoCutConvGamma_Eta" ,"" , kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, "", "");}
if(kGCplotESDCutPionRejectionLowP==kTRUE){histograms->AddHistogram("ESD_CutPionRejectionLowP_InvMass" ,"dedx PionRejection LowP" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutKaonRejectionLowP==kTRUE){histograms->AddHistogram("ESD_CutKaonRejectionLowP_InvMass" ,"dedx KaonRejection LowP" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutQtGammaSelection==kTRUE){histograms->AddHistogram("ESD_CutQt_InvMass","ESD_CutQt_InvMass",kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
+ histograms->AddHistogram("ESD_CutTOFsigmaElec_InvMass", "ESD_CutTOFsigmaElec_InvMass",kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"",""); // RRnewTOF
if(kGCplotESDCutProtonRejectionLowP==kTRUE){histograms->AddHistogram("ESD_CutProtonRejectionLowP_InvMass" ,"dedx ProtonRejection LowP" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutR == kTRUE){histograms->AddHistogram("ESD_CutR_InvMass" ,"Above RMax" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutZ == kTRUE){histograms->AddHistogram("ESD_CutZ_InvMass" ,"Out of reconstruction area" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutMinClsTPC == kTRUE){histograms->AddHistogram("ESD_CutMinNClsTPC_InvMass" ,"Out of reconstruction area" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDCutMinClsTPCToF == kTRUE){histograms->AddHistogram("ESD_CutMinNClsTPCToF_InvMass" ,"Out of reconstruction area" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
+ if(kGCplotESDCutPhotonAsymmetry== kTRUE){histograms->AddHistogram("ESD_CutPhotonAsymmetry_InvMass" ,"Out of reconstruction area" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDGoodV0s == kTRUE){histograms->AddHistogram("ESD_GoodV0s_InvMass" ,"Good V0s" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotESDAllV0s == kTRUE){histograms->AddHistogram("ESD_AllV0s_InvMass" ,"All V0s" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}
if(kGCplotPi0Spectra == kTRUE){
histograms->AddHistogram("ESD_Mother_alfa","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinAlpha, kGClastXBinAlpha,"#alpha","Counts");
histograms->AddHistogram("ESD_Mother_alfa_Pi0","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinAlpha, kGClastXBinAlpha,"#alpha","Counts");
-
+ histograms->AddHistogram("ESD_Mother_alfa_Eta","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinAlpha, kGClastXBinAlpha,"#alpha","Counts");
+ histograms->AddHistogram("ESD_Background_alfa_Pi0","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinAlpha, kGClastXBinAlpha,"#alpha","Counts");
+ histograms->AddHistogram("ESD_Background_alfa_Eta","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinAlpha, kGClastXBinAlpha,"#alpha","Counts");
// histograms->AddHistogram("ESD_Mother_InvMass_vs_Pt" ,"Invariant Mass vs Pt" , kGCnXBinsSpectra, kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Pt [GeV]");
histograms->AddHistogram("ESD_Mother_InvMass_vs_Pt" ,"Invariant Mass vs Pt" , kGCnXBinsSpectra, kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Pt [GeV]");
histograms->AddHistogram("ESD_Mother_InvMass_vs_Pt_alpha" ,"Invariant Mass vs Pt" , kGCnXBinsSpectra, kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Pt [GeV]");
}
// if(kGCdoNeutralMesonV0MCCheck == kTRUE){
- histograms->AddHistogram("ESD_TrueBckGG_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
- histograms->AddHistogram("ESD_TrueBckCont_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
- histograms->AddHistogram("ESD_TruePi0Sec_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
-
+ histograms->AddHistogram("ESD_TrueBckGG_InvMass_vs_Pt","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
+ histograms->AddHistogram("ESD_TrueBckCont_InvMass_vs_Pt","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
+ histograms->AddHistogram("ESD_TruePi0Sec_InvMass_vs_Pt","Invariant mass vs Pt",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
+ histograms->AddHistogram("ESD_TruePi0DalitzCont_InvMass_vs_Pt","Invariant mass vs Pt",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
histograms->AddHistogram("ESD_TruePi0_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
histograms->AddHistogram("ESD_TruePi0_InvMass_1212","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
histograms->AddHistogram("ESD_TruePi0_InvMass_0912","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
if(kGCplotMCPi0Eta == kTRUE){ histograms->AddHistogram("MC_Pi0_Eta" ,"" , kGCnXBinsEta, kGCfirstXBinEta, kGClastXBinEta, "", "");}
if(kGCplotMCPi0Rapid == kTRUE){ histograms->AddHistogram("MC_Pi0_Rapid" ,"" , kGCnXBinsRapid, kGCfirstXBinRapid, kGClastXBinRapid, "", "");}
if(kGCplotMCPi0PtvsRapid == kTRUE){ histograms->AddHistogram("MC_Pi0_Pt_vs_Rapid" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, kGCnXBinsRapid, kGCfirstXBinRapid, kGClastXBinRapid, "", "");}
+ if(kGCplotMCPi0PtvsRapid == kTRUE){ histograms->AddHistogram("MC_Pi0_Pt_vs_Rapid_allDaughters" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, kGCnXBinsRapid, kGCfirstXBinRapid, kGClastXBinRapid, "", "");} // RR primary Pi0 debug
if(kGCplotMCPi0Phi == kTRUE){ histograms->AddHistogram("MC_Pi0_Phi" ,"" , kGCnXBinsPhi, kGCfirstXBinPhi, kGClastXBinPhi, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_SD_EvtQ1_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_SD_EvtQ2_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_SD_EvtQ3_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_SD_EvtQ4_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_SD_EvtQ5_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_DD_EvtQ1_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_DD_EvtQ2_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_DD_EvtQ3_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_DD_EvtQ4_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_DD_EvtQ5_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_ND_EvtQ1_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_ND_EvtQ2_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_ND_EvtQ3_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
-
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_ND_EvtQ4_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
+ if(kGCplotMCPi0Pt == kTRUE){ histograms->AddHistogram("MC_ND_EvtQ5_Pi0_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0PtFiducial == kTRUE){ histograms->AddHistogram("MC_Pi0_Pt_Fiducial" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0PtWithinAcceptanceFiducial == kTRUE){ histograms->AddHistogram("MC_Pi0_Pt_withinAcceptance_Fiducial" ,"" , kGCnXBinsPt,kGCfirstXBinPt, kGClastXBinPt, "", "");}
if(kGCplotMCPi0Energy == kTRUE){ histograms->AddHistogram("MC_Pi0_Energy" ,"" , kGCnXBinsEnergy, kGCfirstXBinEnergy, kGClastXBinEnergy, "", "");}
if(kGCplotMCPi0Mass == kTRUE){ histograms->AddHistogram("MC_Pi0_Mass" ,"" , kGCnXBinsPi0Mass, kGCfirstXBinPi0Mass, kGClastXBinPi0Mass, "", "");}
if(kGCplotMCPi0Alpha == kTRUE){ histograms->AddHistogram("MC_Pi0_alpha" ,"" , kGCnXBinsPi0Mass, kGCfirstXBinPi0Alpha, kGClastXBinPi0Alpha, "", "");}
-
+ if(kGCplotMCEtaAlpha == kTRUE){ histograms->AddHistogram("MC_Eta_alpha" ,"" , kGCnXBinsPi0Mass, kGCfirstXBinPi0Alpha, kGClastXBinPi0Alpha, "", "");}
if(kGCplotMCPi0OpeningAngle == kTRUE){ histograms->AddHistogram("MC_Pi0_GammaDaughter_OpeningAngle" ,"" , kGCnXBinsOpeningAngle, kGCfirstXBinOpeningAngle, kGClastXBinOpeningAngle, "", "");}
if(kGCplotMCPi0R == kTRUE){ histograms->AddHistogram("MC_Pi0_R" ,"" , kGCnXBinsR, kGCfirstXBinR, kGClastXBinR, "", "");}
if(kGCplotMCPi0ZR == kTRUE){ histograms->AddHistogram("MC_Pi0_ZR" ,"" , kGCnXBinsZR, kGCfirstXBinZR, kGClastXBinZR, kGCnYBinsZR, kGCfirstYBinZR, kGClastYBinZR, "", "");}
//--------------------------------------------------- 2 gamma Background -------------------------------------------------------
if(kGCcalculateBackground==kTRUE){
+ histograms->AddHistogram("ESD_GlobalPrimaryVtxZ","Z primary vertex Global",300, -15., 15.,"Z[cm]","counts");
+ histograms->AddHistogram("ESD_SPDPrimaryVtxZ","Z primary vertex SPD",300, -15., 15.,"Z[cm]","counts");
histograms->AddHistogram("ESD_Z_distribution" ,"Z primary vertex" , 2000, -30, 30,"Z[cm]","counts");
histograms->AddHistogram("ESD_multiplicity_distribution" ,"multiplicity distribution" , 200, 0, 200,"counts","Multiplicity");
histograms->AddHistogram("ESD_ZvsMultiplicity" ,"Z vs Multiplicity" , 1000, -10, 10,200,0,200,"Z[cm]","Multiplicity");
Int_t removePileUp=array[20];
Int_t selectV0AND=array[21];
Int_t multiplicityBin=array[22];
-
+ Int_t isHeavyIon=array[23];
+ Int_t useCentrality=array[24];
+ Int_t centralityBin=array[25];
+ Int_t TOFelectronPID=array[26]; // RRnewTOF
+ Int_t useMCPSmearing=array[27];
+ Int_t doPhotonAsymmetryCut=array[28];
+
+ cout<<"doPhotonAsymmetryCut::"<<doPhotonAsymmetryCut<<endl;
+ cout<<"useMCPSmearing::"<<useMCPSmearing<<endl;
+ cout<<"TOFelectronPID: "<<TOFelectronPID<<endl; // RRnewTOF
+ cout<<"CentralityBin::"<< centralityBin <<endl;
+ cout<<"Use Centrality::"<< useCentrality <<endl;
+ cout<<"Heavy Ion::"<< isHeavyIon<<endl;
cout<<"Multiplicity Bin::"<< multiplicityBin<<endl;
cout<<"Select V0AND::"<< selectV0AND<<endl;
cout<<"Remove PileUp::"<< removePileUp<<endl;
cout<<"v0FinderType: "<<v0FinderType <<endl;
cout<<"goodId: "<<goodId <<endl;
+
if(goodId !=9){
cout<<"Analysis Cut Selection too short or does not start with 9"<<endl;
return iResult;
switch(pidedxSigmaCut){
case 0: // -10
kGCPIDnSigmaAbovePionLine=-10;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
break;
case 1: // 0
kGCPIDnSigmaAbovePionLine=0;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
break;
case 2: // 1
kGCPIDnSigmaAbovePionLine=1;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
break;
case 3: // 1
kGCPIDnSigmaAbovePionLine=-1;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
break;
case 4: // 1
kGCPIDnSigmaAbovePionLine=-1.5;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
break;
case 5: // 1
kGCPIDnSigmaAbovePionLine=2.;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
+ break;
+ case 6: // 1
+ kGCPIDnSigmaAbovePionLine=2.;
+ kGCPIDnSigmaAbovePionLineHighPt=0.5;
+ break;
+ case 7: // 1
+ kGCPIDnSigmaAbovePionLine=3.5;
+ kGCPIDnSigmaAbovePionLineHighPt=-10;
+ break;
+ case 8: // 1
+ kGCPIDnSigmaAbovePionLine=2.;
+ kGCPIDnSigmaAbovePionLineHighPt=1.;
break;
-
default:
return iResult;
}
case 5: // 0.3 GeV
kGCPIDMinPnSigmaAbovePionLine=0.3;
break;
- case 6: // 0.3 GeV
+ case 6: // 0.25 GeV // RRnew label was wrong in comment
kGCPIDMinPnSigmaAbovePionLine=0.25;
break;
case 7: // 0.4 GeV
break;
case 4: // 60% of findable clusters
kGCminClsTPCCutToF= 0.6;
+ kGCUseCorrectedTPCClsInfo=0;
break;
case 5: // 0% of findable clusters
kGCminClsTPCCutToF= 0.0;
+ kGCUseCorrectedTPCClsInfo=1;
break;
case 6: // 0% of findable clusters
kGCminClsTPCCutToF= 0.7;
+ kGCUseCorrectedTPCClsInfo=0;
break;
case 7: // 0% of findable clusters
kGCminClsTPCCutToF= 0.35;
+ kGCUseCorrectedTPCClsInfo=0;
+ break;
+ case 8:
+ kGCminClsTPCCutToF= 0.35;
+ kGCUseCorrectedTPCClsInfo=1;
+ break;
+ case 9:
+ kGCminClsTPCCutToF= 0.6;
+ kGCUseCorrectedTPCClsInfo=1;
break;
-
default:
return iResult;
}
kGCPIDnSigmaAtLowPAroundProtonLine=0.;
kGCPIDnSigmaAtLowPAroundPionLine=2.;
break;
-
default:
return iResult;
}
switch(QtMaxCut){
case 0: //
kGCQtMax=1.;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
break;
case 1:
kGCQtMax=0.1;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
break;
case 2:
kGCQtMax=0.07;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
break;
case 3:
kGCQtMax=0.05;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
break;
case 4:
kGCQtMax=0.03;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
+ break;
+ case 5: // RR try to improve (get rid of) low InvMass peak in PbPb
+ kGCQtMax=0.02;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
+ break; // end RR ///////////////////////////////////////////////
+ case 6: // RRnew start: pT dependent qT cut
+ kGCQtMax=0.02;
+ kGCdoHighPtQtGammaSelection=kTRUE;
+ kGCHighPtQtMax=0.06;
+ kGCPtBorderForQt=2.5;
+ break; // RRnew end ////////////////////////////////////////////
+ case 7:
+ kGCQtMax=0.15;
+ kGCdoHighPtQtGammaSelection=kFALSE; // RRnew
+ kGCHighPtQtMax=100.; // RRnew
+ kGCPtBorderForQt=100.; // RRnew
break;
default:
return iResult;
switch(BackgroundScheme){
case 0: //Rotation
kGCUseRotationMethodInBG=kTRUE;
+ kGCdoBGProbability=kFALSE;
break;
- case 1: // mixed event
+ case 1: // mixed event with track multiplicity
kGCUseRotationMethodInBG=kFALSE;
+ kGCUseTrackMultiplicityForBG=kFALSE;
+ kGCdoBGProbability=kFALSE;
+ break;
+ case 2: // mixed event with V0 multiplicity
+ kGCUseRotationMethodInBG=kFALSE;
+ kGCUseTrackMultiplicityForBG=kTRUE;
+ kGCdoBGProbability=kFALSE;
+ break;
+ case 3: //Rotation
+ kGCUseRotationMethodInBG=kTRUE;
+ kGCdoBGProbability=kTRUE;
break;
-
default:
return iResult;
}
case 3:
kGCnumberOfRotationEventsForBG = 20;
break;
-
+ case 4:
+ kGCnumberOfRotationEventsForBG = 2;
+ break;
+ case 5:
+ kGCnumberOfRotationEventsForBG = 50;
+ break;
+ case 6:
+ kGCnumberOfRotationEventsForBG = 80;
+ break;
+ case 7:
+ kGCnumberOfRotationEventsForBG = 100;
+ break;
default:
return iResult;
}
default:
return iResult;
}
- iResult=1;
- return iResult;
+ switch(isHeavyIon){
+ case 0:
+ kGCIsHeavyIon=0;
+ break;
+ case 1:
+ kGCIsHeavyIon=1;
+ break;
+ default:
+ return iResult;
+ }
+
+ switch(useCentrality){
+ case 0:
+ kGCUseCentrality=0;
+ break;
+ case 1:
+ kGCUseCentrality=1;// GetCentralityClass10("V0M")
+ break;
+ case 2:
+ kGCUseCentrality=2; // GetCentralityClass10("CL1")
+ break;
+ case 3://////////////// RRnew start //////////////////////////////////////////////////////////////////////////////////////////////////////
+ kGCUseCentrality=3;// GetCentralityClass10("V0M") but with merged bins: see AliAnalysisTaskGammaConversion.cxx
+ break;
+ case 4:
+ kGCUseCentrality=4;// GetCentralityClass10("CL1") but with merged bins: see AliAnalysisTaskGammaConversion.cxx
+ break;/////////////// RRnew end ////////////////////////////////////////////////////////////////////////////////////////////////////////
+ default:
+ return iResult;
+ }
+
+ switch(centralityBin){
+ case 0:
+ kGCUseCentralityBin=0;
+ break;
+ case 1:
+ kGCUseCentralityBin=1;
+ break;
+ case 2:
+ kGCUseCentralityBin=2;
+ break;
+ case 3:
+ kGCUseCentralityBin=3;
+ break;
+ case 4:
+ kGCUseCentralityBin=4;
+ break;
+ case 5:
+ kGCUseCentralityBin=5;
+ break;
+ case 6:
+ kGCUseCentralityBin=6;
+ break;
+ case 7:
+ kGCUseCentralityBin=7;
+ break;
+ case 8:
+ kGCUseCentralityBin=8;
+ break;
+ case 9:
+ kGCUseCentralityBin=9;
+ break;
+ default:
+ return iResult;
+ }
+
+ switch(TOFelectronPID){ // RRnewTOF start //////////////////////////////////////////////////////////////////////////
+ case 0: // no cut
+ kGCuseTOFpid = kFALSE;
+ kGCtofPIDnSigmaBelowElectronLine=-100;
+ kGCtofPIDnSigmaAboveElectronLine=100;
+ break;
+ case 1: // -7,7
+ kGCuseTOFpid = kTRUE;
+ kGCtofPIDnSigmaBelowElectronLine=-7;
+ kGCtofPIDnSigmaAboveElectronLine=7;
+ break;
+ case 2: // -5,5
+ kGCuseTOFpid = kTRUE;
+ kGCtofPIDnSigmaBelowElectronLine=-5;
+ kGCtofPIDnSigmaAboveElectronLine=5;
+ break;
+ case 3: // -3,5
+ kGCuseTOFpid = kTRUE;
+ kGCtofPIDnSigmaBelowElectronLine=-3;
+ kGCtofPIDnSigmaAboveElectronLine=5;
+ break;
+ case 4: // -2,3
+ kGCuseTOFpid = kTRUE;
+ kGCtofPIDnSigmaBelowElectronLine=-2;
+ kGCtofPIDnSigmaAboveElectronLine=3;
+ break;
+ default:
+ return iResult;
+ } //////////////////////// RRnewTOF end //////////////////////////////////////////////////////////////////////////
+
+ switch(useMCPSmearing){
+ case 0:
+ kGCUseMCPSmearing=0;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.;
+ kGCPSigSmearingCte=0.;
+ break;
+ case 1:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.0e-14;
+ kGCPSigSmearing=0.;
+ kGCPSigSmearingCte=0.;
+ break;
+ case 2:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.0e-15;
+ kGCPSigSmearing=0.0;
+ kGCPSigSmearingCte=0.;
+ break;
+ case 3:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.003;
+ kGCPSigSmearingCte=0.002;
+ break;
+ case 4:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.003;
+ kGCPSigSmearingCte=0.007;
+ break;
+ case 5:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.003;
+ kGCPSigSmearingCte=0.016;
+ break;
+ case 6:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.007;
+ kGCPSigSmearingCte=0.016;
+ break;
+ case 7:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.0e-16;
+ kGCPSigSmearing=0.0;
+ kGCPSigSmearingCte=0.;
+ break;
+ case 8:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.007;
+ kGCPSigSmearingCte=0.014;
+ break;
+ case 9:
+ kGCUseMCPSmearing=1;
+ kGCPBremSmearing=1.;
+ kGCPSigSmearing=0.007;
+ kGCPSigSmearingCte=0.011;
+ break;
+
+ default:
+ return iResult;
+ }
+ switch(doPhotonAsymmetryCut){
+ case 0:
+ kGCdoPhotonAsymmetryCut=0;
+ kGCMinPPhotonAsymmetryCut=100.;
+ kGCMinPhotonAsymmetry=0.;
+ break;
+ case 1:
+ kGCdoPhotonAsymmetryCut=1;
+ kGCMinPPhotonAsymmetryCut=3.5;
+ kGCMinPhotonAsymmetry=0.04;
+ break;
+ case 2:
+ kGCdoPhotonAsymmetryCut=1;
+ kGCMinPPhotonAsymmetryCut=3.5;
+ kGCMinPhotonAsymmetry=0.06;
+ break;
+ default:
+ return iResult;
+ }
+
+ iResult = 1;
+ return iResult;
}
ASSIGNARRAY(20);
ASSIGNARRAY(21);
ASSIGNARRAY(22);
+ ASSIGNARRAY(23);
+ ASSIGNARRAY(24);
+ ASSIGNARRAY(25);
+ ASSIGNARRAY(26); // RRnewTOF
+ ASSIGNARRAY(27);
+ ASSIGNARRAY(28);
}
}