DY : DYGto2LG-1Jets_MLL-10to50_PTG-XXtoXX_amcatnloFXFX-pythia8

[sihyunjeon]

No description provided.

[sihyunjeon]

can you check this @agrohsje ? cuts.f and run_card.dat ... hope i didn't make any mistakes

[sihyunjeon]

can you also follow the logic in the link below fortran makes me crazy?

alliso i am not sure what this is doing

GridpackFiles/Cards/MadGraph5_aMCatNLO/DY/DYGto2LG-1Jets_MLL-10to50_PTG-100to200_amcatnloFXFX-pythia8/DYGto2LG-1Jets_MLL-10to50_PTG-100to200_amcatnloFXFX-pythia8_cuts.f

Lines 278 to 335 in dd890aa

	do i=1,nexternal
	if (istatus(i).eq.1 .and. ipdg(i).eq.22) then
	is_a_ph(i)=.true.
	else
	is_a_ph(i)=.false.
	endif
	enddo
	if (ptgmin.ne.0d0) then
	nph=0
	do j=nincoming+1,nexternal
	if (is_a_ph(j)) then
	nph=nph+1
	do i=0,3
	pgamma(i,nph)=p(i,j)
	enddo
	endif
	enddo
	if(nph.eq.0)goto 444
	if(isoEM)then
	nem=nph
	do k=1,nem
	do i=0,3
	pem(i,k)=pgamma(i,k)
	enddo
	enddo
	do j=nincoming+1,nexternal
	if (is_a_lp(j).or.is_a_lm(j)) then
	nem=nem+1
	do i=0,3
	pem(i,nem)=p(i,j)
	enddo
	endif
	enddo
	endif
	alliso=.true.
	j=0
	do while(j.lt.nph.and.alliso)
	c Loop over all photons
	j=j+1
	ptg=pt(pgamma(0,j))
	if(ptg.lt.ptgmin)then
	passcuts_user=.false.
	return
	endif
	if(ptg.gt.230)then
	passcuts_user=.false.
	return
	endif
	if (etagamma.gt.0d0) then
	if (abs(eta(pgamma(0,j))).gt.etagamma) then
	passcuts_user=.false.
	return
	endif
	endif

[sihyunjeon]

@agrohsje reminder

[agrohsje]

One thought: Should we make the cuts.f snipped more like a patch so we are less version dependent? The code allow user patches, so everything should be in place to do it. I am also happy to keep going now. But just in the future we might consider. What do you think?