hESDTrueMotherGGInvMassPt(NULL),
hESDTrueConvGammaPt(NULL),
hESDTrueConvGammaR(NULL),
+ hESDTrueConvGammaRMC(NULL),
hESDTruePositronPt(NULL),
hESDTrueElectronPt(NULL),
hESDTrueSecConvGammaPt(NULL),
hESDTrueMotherGGInvMassPt(NULL),
hESDTrueConvGammaPt(NULL),
hESDTrueConvGammaR(NULL),
+ hESDTrueConvGammaRMC(NULL),
hESDTruePositronPt(NULL),
hESDTrueElectronPt(NULL),
hESDTrueSecConvGammaPt(NULL),
hESDEposEnegTruePhotonPsiPairDPhiPtCut = new TH2F*[fnCuts];
hESDEposEnegTrueJPsiInvMassPt = new TH2F*[fnCuts];
hESDTrueConvGammaR = new TH1F*[fnCuts];
+ hESDTrueConvGammaRMC = new TH1F*[fnCuts];
}
hESDTrueConvGammaR[iCut] = new TH1F("ESD_TrueConvGamma_R","ESD_TrueConvGamma_R",800,0,200);
fTrueList[iCut]->Add(hESDTrueConvGammaR[iCut]);
+
+ hESDTrueConvGammaRMC[iCut] = new TH1F("ESD_TrueConvGamma_R_MC","ESD_TrueConvGamma_R_MC",800,0,200);
+ fTrueList[iCut]->Add(hESDTrueConvGammaRMC[iCut]);
+
}
hESDTruePositronPt[iCut] = new TH1F("ESD_TruePositron_Pt","ESD_TruePositron_Pt",1000,0,25);
if( labelGamma < MCStack->GetNprimary() ){
if( fIsFromMBHeader ){
hESDTrueConvGammaPt[fiCut]->Fill(TruePhotonCandidate->Pt());
- if(fDoMesonQA) hESDTrueConvGammaR[fiCut]->Fill(TruePhotonCandidate->GetConversionRadius());
+ if(fDoMesonQA){
+ hESDTrueConvGammaR[fiCut]->Fill(TruePhotonCandidate->GetConversionRadius());
+ hESDTrueConvGammaRMC[fiCut]->Fill( posDaughter->R() );
+ }
}
} else {
if( fIsFromMBHeader){
if(trainConfig == 9 || trainConfig == 10 ){
cutnumberPhoton = "160084001001500000000"; //Offline V0 finder
- }
- else {
- cutnumberPhoton = "060084001001500000000"; //Online V0 finder
+ } else if ( trainConfig == 54 ){ // Warning!!! Those trains must run in stand-alone mode
+ cutnumberPhoton = "000700040000500000000";
+ } else {
+ cutnumberPhoton = "060084001000500000000"; //Online V0 finder //change
}
TString ElecCuts = "90005400000002000000"; //Electron Cuts
eventCutArray[0]="8000011"; photonCutArray[0] = "002093603007200000000"; ElecCutarray[0] = "90475400233202623700"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + New psi pair cut + New psi pair cut fPsiPairCut = 0.60 fDeltaPhiCutMin = 0.0 fDeltaPhiCutMax = 0.12
+} else if( trainConfig == 53 ) {
+
+ eventCutArray[0]="8000011"; photonCutArray[0] = "002093603007200000000"; ElecCutarray[0] = "90475400233202623700"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + New psi pair cut + New psi pair cut fPsiPairCut = 0.60 fDeltaPhiCutMin = 0.0 fDeltaPhiCutMax = 0.12
+
+} else if( trainConfig == 54 ) {
+
+ eventCutArray[0]="8000011"; photonCutArray[0] = "000700040000500000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 //Tracks 2011 + kAny + new psiPair Cut 0.60, 0.0 0.12 + 0 cm < Rconv < 180 cm //open cuts
+
}
if ( ( trainConfig >= 1 && trainConfig <= 6 ) || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 14 || trainConfig == 16 || trainConfig == 17 || trainConfig == 18 || trainConfig == 19 || trainConfig == 20 || trainConfig == 21 || trainConfig == 22 || trainConfig == 23 ||
trainConfig == 28 || trainConfig == 29 || trainConfig == 30 || trainConfig == 31 || trainConfig == 32 || trainConfig == 33 || trainConfig == 37 || trainConfig == 38 || trainConfig == 39 || trainConfig == 40 || trainConfig == 41 || trainConfig == 41 || trainConfig == 43 || trainConfig == 44 ||
- trainConfig == 45 || trainConfig == 46 || trainConfig == 47 || trainConfig == 48 || trainConfig == 49 || trainConfig == 50 || trainConfig == 51 || trainConfig == 52 ) {
+ trainConfig == 45 || trainConfig == 46 || trainConfig == 47 || trainConfig == 48 || trainConfig == 49 || trainConfig == 50 || trainConfig == 51 || trainConfig == 52 || trainConfig == 54 ) {
if (doWeighting){
if (generatorName.CompareTo("DPMJET")==0){