void lhapdfdraw(Int_t nfiles=1, Int_t process=3)
{
  Int_t IDs[3];
  if (process == 1) {IDs[0]=11; IDs[1]=12; IDs[2]=13;} //ff->jets
  if (process == 2) {IDs[0]=28; IDs[1]=28; IDs[2]=28;} //fg->gets
  if (process == 3) {IDs[0]=53; IDs[1]=68; IDs[2]=68;} //gg->jets
  
  TH1F *hEtaJets3 = new TH1F("hEtaJets3", "hEtaJets3", 50, -5.5, 5.5);
  TH1F *hEtaJets5 = new TH1F("hEtaJets5", "hEtaJets5", 50, -5.5, 5.5);
  TH1F *hEtaJets7 = new TH1F("hEtaJets7", "hEtaJets7", 50, -5.5, 5.5);

  TH1F *hPtJets3 = new TH1F("hPtJets3", "hPtJets3", 50, 0, 150);
  TH1F *hPtJets5 = new TH1F("hPtJets5", "hPtJets5", 50, 0, 150);
  TH1F *hPtJets7 = new TH1F("hPtJets7", "hPtJets7", 50, 0, 150);

  TH1F *hPtJets3_eta1 = new TH1F("hPtJets3_eta1", "hPtJets3_eta1", 50, 0, 150);
  TH1F *hPtJets3_eta2 = new TH1F("hPtJets3_eta2", "hPtJets3_eta2", 50, 0, 150);
  TH1F *hPtJets3_eta3 = new TH1F("hPtJets3_eta3", "hPtJets3_eta3", 50, 0, 150);
  TH1F *hPtJets5_eta1 = new TH1F("hPtJets5_eta1", "hPtJets5_eta1", 50, 0, 150);
  TH1F *hPtJets5_eta2 = new TH1F("hPtJets5_eta2", "hPtJets5_eta2", 50, 0, 150);
  TH1F *hPtJets5_eta3 = new TH1F("hPtJets5_eta3", "hPtJets5_eta3", 50, 0, 150);
  TH1F *hPtJets7_eta1 = new TH1F("hPtJets7_eta1", "hPtJets7_eta1", 50, 0, 150);
  TH1F *hPtJets7_eta2 = new TH1F("hPtJets7_eta2", "hPtJets7_eta2", 50, 0, 150);
  TH1F *hPtJets7_eta3 = new TH1F("hPtJets7_eta3", "hPtJets7_eta3", 50, 0, 150);

  TH1F *x1_3 = new TH1F("x1_3", "x1_3", 50, -6, 0);
  TH1F *x1_5 = new TH1F("x1_5", "x1_5", 50, -6, 0);
  TH1F *x1_7 = new TH1F("x1_7", "x1_7", 50, -6, 0);
  TH1F *x2_3 = new TH1F("x2_3", "x2_3", 50, -6, 0);
  TH1F *x2_5 = new TH1F("x2_5", "x2_5", 50, -6, 0);
  TH1F *x2_7 = new TH1F("x2_7", "x2_7", 50, -6, 0);
  TH1F *xf_3 = new TH1F("xf_3", "xf_3", 50, -6, 0);
  TH1F *xf_5 = new TH1F("xf_5", "xf_5", 50, -6, 0);
  TH1F *xf_7 = new TH1F("xf_7", "xf_7", 50, -6, 0);
  TH1F *XX_3 = new TH1F("XX_3", "XX_3", 50, -6, 0);
  TH1F *XX_5 = new TH1F("XX_5", "XX_5", 50, -6, 0);
  TH1F *XX_7 = new TH1F("XX_7", "XX_7", 50, -6, 0);
  TH1F *xx_3 = new TH1F("xx_3", "xx_3", 50, -6, 0);
  TH1F *xx_5 = new TH1F("xx_5", "xx_5", 50, -6, 0);
  TH1F *xx_7 = new TH1F("xx_7", "xx_7", 50, -6, 0);

  TH1F *Q2_3 = new TH1F("Q2_3", "Q2_3", 50, 0, 500);
  TH1F *Q2_5 = new TH1F("Q2_5", "Q2_5", 50, 0, 500);
  TH1F *Q2_7 = new TH1F("Q2_7", "Q2_7", 50, 0, 500);


  TH2F *x1_Q2_3 = new TH2F("x1_Q2_3", "x1_Q2_3", 50, -6,0, 50, 0,500 );
  TH2F *x2_Q2_3 = new TH2F("x2_Q2_3", "x2_Q2_3", 50, -6,0, 50, 0,500 );
  TH2F *x1_Q2_5 = new TH2F("x1_Q2_5", "x1_Q2_5", 50, -6,0, 50, 0,500 );
  TH2F *x2_Q2_5 = new TH2F("x2_Q2_5", "x2_Q2_5", 50, -6,0, 50, 0,500 );
  TH2F *x1_Q2_7 = new TH2F("x1_Q2_7", "x1_Q2_7", 50, -6,0, 50, 0,500 );
  TH2F *x2_Q2_7 = new TH2F("x2_Q2_7", "x2_Q2_7", 50, -6,0, 50, 0,500 );

  TH2F *XX_pt_3 = new TH2F("XX_pt_3", "x vs Pt of jets", 100, -6,0, 50, 0,120 );
  TH2F *xx_pt_3 = new TH2F("xx_pt_3", "x vs Pt of jets", 100, -6,0, 50, 0,120 );
  TH2F *XX_pt_5 = new TH2F("XX_pt_5", "x vs Pt of jets", 100, -6,0, 50, 0,120 );
  TH2F *xx_pt_5 = new TH2F("xx_pt_5", "x vs Pt of jets", 100, -6,0, 50, 0,120 );
  TH2F *XX_pt_7 = new TH2F("XX_pt_7", "x vs Pt of jets", 100, -6,0, 50, 0,120 );
  TH2F *xx_pt_7 = new TH2F("xx_pt_7", "x vs Pt of jets", 100, -6,0, 50, 0,120 );

  TH2F *XX_xx_1_1_3 = new TH2F("XX_xx_1_1_3", "X1-x2 at |eta|<1.3 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_2_1_3 = new TH2F("XX_xx_2_1_3", "X1-x2 at 1.3<|eta|<3.0 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_3_1_3 = new TH2F("XX_xx_3_1_3", "X1-x2 at 3.0<|eta|<5.0 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_1_2_3 = new TH2F("XX_xx_1_2_3", "X1-x2 at |eta|<1.3 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_2_2_3 = new TH2F("XX_xx_2_2_3", "X1-x2 at 1.3<|eta|<3.0 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_3_2_3 = new TH2F("XX_xx_3_2_3", "X1-x2 at 3.0<|eta|<5.0 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_1_3_3 = new TH2F("XX_xx_1_3_3", "X1-x2 at |eta|<1.3 and 60<|pt|; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_2_3_3 = new TH2F("XX_xx_2_3_3", "X1-x2 at 1.3<|eta|<3.0 and 60<|pt|; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_3_3_3 = new TH2F("XX_xx_3_3_3", "X1-x2 at 3.0<|eta|<5.0 and 60<|pt|; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );

  TH2F *XX_xx_1_1_5 = new TH2F("XX_xx_1_1_5", "X1-x2 at |eta|<1.3 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_2_1_5 = new TH2F("XX_xx_2_1_5", "X1-x2 at 1.3<|eta|<3.0 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_3_1_5 = new TH2F("XX_xx_3_1_5", "X1-x2 at 3.0<|eta|<5.0 and 10<|pt|<30; log10(X1) Large; log10(x2) small", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_1_2_5 = new TH2F("XX_xx_1_2_5", "X1-x2 at |eta|<1.3 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_2_2_5 = new TH2F("XX_xx_2_2_5", "X1-x2 at 1.3<|eta|<3.0 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_3_2_5 = new TH2F("XX_xx_3_2_5", "X1-x2 at 3.0<|eta|<5.0 and 30<|pt|<60; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_1_3_5 = new TH2F("XX_xx_1_3_5", "X1-x2 at |eta|<1.3 and 60<|pt|; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_2_3_5 = new TH2F("XX_xx_2_3_5", "X1-x2 at 1.3<|eta|<3.0 and 60<|pt|; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_3_3_5 = new TH2F("XX_xx_3_3_5", "X1-x2 at 3.0<|eta|<5.0 and 60<|pt|; log10(X1) Large; log10(x2) small", 50, -4,0, 50, -6,0 );

  TH2F *XX_xx_1_1_7 = new TH2F("XX_xx_1_1_7", "XX_xx_1_1_7", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_2_1_7 = new TH2F("XX_xx_2_1_7", "XX_xx_2_1_7", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_3_1_7 = new TH2F("XX_xx_3_1_7", "XX_xx_3_1_7", 100, -4,0, 100, -6,0 );
  TH2F *XX_xx_1_2_7 = new TH2F("XX_xx_1_2_7", "XX_xx_1_2_7", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_2_2_7 = new TH2F("XX_xx_2_2_7", "XX_xx_2_2_7", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_3_2_7 = new TH2F("XX_xx_3_2_7", "XX_xx_3_2_7", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_1_3_7 = new TH2F("XX_xx_1_3_7", "XX_xx_1_3_7", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_2_3_7 = new TH2F("XX_xx_2_3_7", "XX_xx_2_3_7", 50, -4,0, 50, -6,0 );
  TH2F *XX_xx_3_3_7 = new TH2F("XX_xx_3_3_7", "XX_xx_3_3_7", 50, -4,0, 50, -6,0 );

  TH2F *jets_Q2_3 = new TH2F("jets_Q2_3", "mass j_{1} + j_{2} vs Q^{2}", 50, 0,500, 50, 0,500);
  TH2F *jets_Q2_5 = new TH2F("jets_Q2_5", "mass j_{1} + j_{2} vs Q^{2}", 50, 0,500, 50, 0,500);
  TH2F *jets_Q2_7 = new TH2F("jets_Q2_7", "mass j_{1} + j_{2} vs Q^{2}", 50, 0,500, 50, 0,500);


  //-----------------------------------------------------------------------
  
  Int_t njet2, IDproc;
  Float_t pxj2[50], pyj2[50], pzj2[50], ptj2[50], etaj2[50];
  Float_t ptJet[50], MJETS;
  Float_t x1, x2, xf, Q2, XX, xx;  
  Float_t sigma[4], coef, cs[3];
  Float_t LUMI = 1000.;                               //Luminosity needed nb^-1
  Int_t nevt;					//number of events
  // Int_t nfiles;					//number of files
  Int_t iFlag;
  //-----------------------------------------------------------------------
  
  // nfiles = 3;
  nevt = nfiles*200000;
  
  for (Int_t file_num = 1; file_num <= nfiles; file_num++)
    {
      cout  <<" file_num "<< file_num << endl;
      TFile *f1_1 =  TFile::Open(Form("rfio:/castor/cern.ch/user/a/andrey/lhapdf/pp14TeV_lha_%d.root", file_num));
      TTree *h107 = (TTree*)f1_1->Get("h101");
      TTree *h105 = (TTree*)f1_1->Get("h102");
      TTree *h103 = (TTree*)f1_1->Get("h103");
      TTree *h111 = (TTree*)f1_1->Get("h111");
      
      h111->SetBranchAddress("cs",&cs);
      
      h103->SetBranchAddress("njet2",&njet2);
      h103->SetBranchAddress("etaj2",&etaj2);
      h103->SetBranchAddress("pxj2",&pxj2);
      h103->SetBranchAddress("pyj2",&pyj2);
      h103->SetBranchAddress("pzj2",&pzj2);
      h103->SetBranchAddress("ptj2",&ptj2);
      h103->SetBranchAddress("x1",&x1);
      h103->SetBranchAddress("x2",&x2);
      h103->SetBranchAddress("xf",&xf);
      h103->SetBranchAddress("Q2",&Q2);
      h103->SetBranchAddress("Idproc",&IDproc);
      
      h111->GetEntry(0);   
      sigma[1] = cs[2]*1000000; 
      sigma[2] = cs[1]*1000000; 
      sigma[3] = cs[0]*1000000; 

      cout<<" sigma[1] "<< sigma[1] << endl;
      cout<<" sigma[2] "<< sigma[2] << endl;
      cout<<" sigma[3] "<< sigma[3] << endl;

      coef = sigma[1]*LUMI/nevt;
      cout  << "coef for PDF 'set 1'  ->> "<<coef <<endl;
      
      Int_t nentries = (Int_t)h103->GetEntries();
      for (Int_t i=0; i<nentries; i++)
	{
	  h103->GetEntry(i);
	  
	  if(IDproc == IDs[0] || IDproc == IDs[1] || IDproc == IDs[2])
	    {
	      
	      x1_3 -> Fill(log10(x1));
	      x2_3 -> Fill(log10(x2));
	      xf_3 -> Fill(log10(fabs(xf)));
	      
	      XX = x1 > x2  ? x1 : x2;
	      xx = x1 > x2  ? x2 : x1;
	      //cout <<"x1:"<<x1<<"  x2:"<<x2<<"  XX:"<<XX<<"  xx:"<<xx<<endl;
	      XX_3 -> Fill(log10(XX));
	      xx_3 -> Fill(log10(xx));
	      
	      Q2_3 -> Fill(Q2);
	      x1_Q2_3 -> Fill(log10(x1), Q2); 
	      x2_Q2_3 -> Fill(log10(x2), Q2); 
	      
	    for (Int_t k=0; k<njet2; k++)
	    {
	    	for (Int_t l=k; l<njet2; l++)
	    	{
	              if(k!=l && fabs(etaj2[k])>3.0  && fabs(etaj2[l])>3.0 && etaj2[k]*etaj2[l]<0)
			{
	                  MJETS = ptj2[k]*ptj2[l];
        	          jets_Q2_3 -> Fill(MJETS, Q2);
//		  cout<<"ALT: "<<njet2<<" pt1 "<<ptj2[k]<<" pt2 "<<ptj2[l]<<" eta1 "<<etaj2[k]<<" eta2 "<<etaj2[l]<<k<<l<<endl;
			k=njet2;
			break;
			}
	    	}	    
	    }
	  
	      
	      iFlag=0;
	      for (Int_t j=0; j<njet2; j++)
		{
		  if (fabs(etaj2[j]) > 3.0){
	  		  XX_pt_3 -> Fill(log10(XX), ptj2[j]); 
			  xx_pt_3 -> Fill(log10(xx), ptj2[j]);}
		  
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])<1.3) {XX_xx_1_1_3 -> Fill(log10(XX), log10(xx)) ;}
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_1_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_1_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])<1.3) {XX_xx_1_2_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_2_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_2_3 -> Fill(log10(XX), log10(xx));}
	      
	      if(ptj2[j]>60 && fabs(etaj2[j])<1.3) {XX_xx_1_3_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_3_3 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_3_3 -> Fill(log10(XX), log10(xx));}
	      
	      hPtJets3 -> Fill(ptj2[j], coef);
	      hEtaJets3 -> Fill(etaj2[j], coef);
	      
	      if (fabs(etaj2[j]) <= 1.3) {hPtJets3_eta1 -> Fill(ptj2[j], coef);}
	      if (1.3 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 3.0) {hPtJets3_eta2 -> Fill(ptj2[j], coef);}
	      if (3.0 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 5.0) {hPtJets3_eta3 -> Fill(ptj2[j], coef);}
	    }
	  
	}
    }
  
  h105->SetBranchAddress("njet2",&njet2);
  h105->SetBranchAddress("etaj2",&etaj2);
  h105->SetBranchAddress("pxj2",&pxj2);
  h105->SetBranchAddress("pyj2",&pyj2);
  h105->SetBranchAddress("pzj2",&pzj2);
  h105->SetBranchAddress("ptj2",&ptj2);
  h105->SetBranchAddress("x1",&x1);
  h105->SetBranchAddress("x2",&x2);
  h105->SetBranchAddress("xf",&xf);
  h105->SetBranchAddress("Q2",&Q2);
  h105->SetBranchAddress("Idproc",&IDproc);
  
  
  coef = sigma[2]*LUMI/nevt;
  cout  << "coef for PDF 'set 2'  ->> "<<coef <<endl;
  
  Int_t nentries = (Int_t)h105->GetEntries();
  for (Int_t i=0; i<nentries; i++)
    {
      h105->GetEntry(i);
      
      if(IDproc == IDs[0] || IDproc == IDs[1] || IDproc == IDs[2])
	{
	  x1_5 -> Fill(log10(x1));
	  x2_5 -> Fill(log10(x2));
	  xf_5 -> Fill(log10(fabs(xf)));
          XX = x1 > x2  ? x1 : x2;
          xx = x1 > x2  ? x2 : x1;
	  XX_5 -> Fill(log10(XX));
	  xx_5 -> Fill(log10(xx));
	  Q2_5 -> Fill(Q2);
	  x1_Q2_5 -> Fill(log10(x1), Q2); 
	  x2_Q2_5 -> Fill(log10(x2), Q2); 
	  

	    for (Int_t k=0; k<njet2; k++)
	    {
	    	for (Int_t l=k; l<njet2; l++)
	    	{
	              if(k!=l && fabs(etaj2[k])>3.0  && fabs(etaj2[l])>3.0 && etaj2[k]*etaj2[l]<0)
			{
	                  MJETS = ptj2[k]*ptj2[l];
        	          jets_Q2_5 -> Fill(MJETS, Q2);
//		  cout<<"ALT: "<<njet2<<" pt1 "<<ptj2[k]<<" pt2 "<<ptj2[l]<<" eta1 "<<etaj2[k]<<" eta2 "<<etaj2[l]<<k<<l<<endl;
			k=njet2;
			break;
			}
	    	}	    
	    }
	  
	  iFlag=0;
	  for (Int_t j=0; j<njet2; j++)
	    {
	      if (fabs(etaj2[j]) > 3.0){
		XX_pt_5 -> Fill(log10(XX), ptj2[j]); 
		xx_pt_5 -> Fill(log10(xx), ptj2[j]);}
	      
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])<1.3) {XX_xx_1_1_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_1_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_1_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])<1.3) {XX_xx_1_2_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_2_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_2_5 -> Fill(log10(XX), log10(xx));}
	      
	      if(ptj2[j]>60 && fabs(etaj2[j])<1.3) {XX_xx_1_3_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_3_5 -> Fill(log10(XX), log10(xx));}
	      if(ptj2[j]>60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_3_5 -> Fill(log10(XX), log10(xx));}
	      
	      hPtJets5 -> Fill(ptj2[j], coef);
	      hEtaJets5 -> Fill(etaj2[j], coef);
	      
	      if (fabs(etaj2[j]) <= 1.3) {hPtJets5_eta1 -> Fill(ptj2[j], coef);}
	      if (1.3 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 3.0) {hPtJets5_eta2 -> Fill(ptj2[j], coef);}
	      if (3.0 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 5.0) {hPtJets5_eta3 -> Fill(ptj2[j], coef);}
	    }
	  
	}
    }
h107->SetBranchAddress("njet2",&njet2);
h107->SetBranchAddress("etaj2",&etaj2);
h107->SetBranchAddress("pxj2",&pxj2);
h107->SetBranchAddress("pyj2",&pyj2);
h107->SetBranchAddress("pzj2",&pzj2);
h107->SetBranchAddress("ptj2",&ptj2);
h107->SetBranchAddress("x1",&x1);
h107->SetBranchAddress("x2",&x2);
h107->SetBranchAddress("xf",&xf);
h107->SetBranchAddress("Q2",&Q2);
h107->SetBranchAddress("Idproc",&IDproc);


coef = sigma[3]*LUMI/nevt;
cout  << "coef for PDF 'set 3'  ->> "<<coef <<endl;
Int_t nentries = (Int_t)h107->GetEntries();
for (Int_t i=0; i<nentries; i++)
{
  h107->GetEntry(i);
  if(IDproc == IDs[0] || IDproc == IDs[1] || IDproc == IDs[2])
    {
      x1_7 -> Fill(log10(x1));
      x2_7 -> Fill(log10(x2));
      xf_7 -> Fill(log10(fabs(xf)));
      XX = x1 > x2  ? x1 : x2;
      xx = x1 > x2  ? x2 : x1;
      XX_7 -> Fill(log10(XX));
      xx_7 -> Fill(log10(xx));
      Q2_7 -> Fill(Q2);
      x1_Q2_7 -> Fill(log10(x1), Q2); 
      x2_Q2_7 -> Fill(log10(x2), Q2); 
      
	    for (Int_t k=0; k<njet2; k++)
	    {
	    	for (Int_t l=k; l<njet2; l++)
	    	{
	              if(k!=l && fabs(etaj2[k])>3.0  && fabs(etaj2[l])>3.0 && etaj2[k]*etaj2[l]<0)
			{
	                  MJETS = ptj2[k]*ptj2[l];
        	          jets_Q2_7 -> Fill(MJETS, Q2);
//		  cout<<"ALT: "<<njet2<<" pt1 "<<ptj2[k]<<" pt2 "<<ptj2[l]<<" eta1 "<<etaj2[k]<<" eta2 "<<etaj2[l]<<k<<l<<endl;
			k=njet2;
			break;
			}
	    	}	    
	    }
	  
      iFlag=0;
	  for (Int_t j=0; j<njet2; j++)
	    {
		  if (fabs(etaj2[j]) > 3.0){
		   XX_pt_7 -> Fill(log10(XX), ptj2[j]); 
		   xx_pt_7 -> Fill(log10(xx), ptj2[j]);}

	  if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])<1.3) {XX_xx_1_1_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_1_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>10 && ptj2[j]<30 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_1_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])<1.3) {XX_xx_1_2_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_2_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>30 && ptj2[j]<60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_2_7 -> Fill(log10(XX), log10(xx));}
	  
	  if(ptj2[j]>60 && fabs(etaj2[j])<1.3) {XX_xx_1_3_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>60 && fabs(etaj2[j])>1.3 && etaj2[j]<3.0) {XX_xx_2_3_7 -> Fill(log10(XX), log10(xx));}
	  if(ptj2[j]>60 && fabs(etaj2[j])>3.0 && etaj2[j]<5.0) {XX_xx_3_3_7 -> Fill(log10(XX), log10(xx));}
	  
	  hPtJets7 -> Fill(ptj2[j], coef);
	  hEtaJets7 -> Fill(etaj2[j], coef);
	  if (fabs(etaj2[j]) <= 1.3) {hPtJets7_eta1 -> Fill(ptj2[j], coef);}
	  if (1.3 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 3.0) {hPtJets7_eta2 -> Fill(ptj2[j], coef);}
	  if (3.0 < fabs(etaj2[j]) && fabs(etaj2[j]) <= 5.0) {hPtJets7_eta3 -> Fill(ptj2[j], coef);}
           }
  
       }
   }
}

if (process == 1) {TFile ftemp("temp21_ff.root", "recreate");}
if (process == 2) {TFile ftemp("temp22_fg.root", "recreate");}
if (process == 3) {TFile ftemp("temp23_gg.root", "recreate");}
hPtJets3 -> Sumw2();
hPtJets3_eta1 -> Sumw2();
hPtJets3_eta2 -> Sumw2();
hPtJets3_eta3 -> Sumw2();
hPtJets3 -> Write();
hPtJets3_eta1 -> Write();
hPtJets3_eta2 -> Write();
hPtJets3_eta3 -> Write();

hPtJets5 -> Sumw2();
hPtJets5_eta1 -> Sumw2();
hPtJets5_eta2 -> Sumw2();
hPtJets5_eta3 -> Sumw2();
hPtJets5 -> Write();
hPtJets5_eta1 -> Write();
hPtJets5_eta2 -> Write();
hPtJets5_eta3 -> Write();

hPtJets7 -> Sumw2();
hPtJets7_eta1 -> Sumw2();
hPtJets7_eta2 -> Sumw2();
hPtJets7_eta3 -> Sumw2();
hPtJets7 -> Write();
hPtJets7_eta1 -> Write();
hPtJets7_eta2 -> Write();
hPtJets7_eta3 -> Write();

hEtaJets3 -> Sumw2();
hEtaJets5 -> Sumw2();
hEtaJets7 -> Sumw2();
hEtaJets3 -> Write();
hEtaJets5 -> Write();
hEtaJets7 -> Write();


x1_3 -> Write();
x1_5 -> Write();
x1_7 -> Write();
x2_3 -> Write();
x2_5 -> Write();
x2_7 -> Write();
xf_3 -> Write();
xf_5 -> Write();
xf_7 -> Write();

XX_3 -> Write();
XX_5 -> Write();
XX_7 -> Write();
xx_3 -> Write();
xx_5 -> Write();
xx_7 -> Write();

Q2_3 -> Write();
Q2_5 -> Write();
Q2_7 -> Write();
x1_Q2_3 -> Write(); 
x2_Q2_3 -> Write(); 
x1_Q2_5 -> Write(); 
x2_Q2_5 -> Write(); 
x1_Q2_7 -> Write(); 
x2_Q2_7 -> Write(); 

XX_pt_3 -> Write(); 
xx_pt_3 -> Write(); 
XX_pt_5 -> Write(); 
xx_pt_5 -> Write(); 
XX_pt_7 -> Write(); 
xx_pt_7 -> Write(); 

XX_xx_1_1_3 -> Write(); 
XX_xx_2_1_3 -> Write(); 
XX_xx_3_1_3 -> Write(); 
XX_xx_1_2_3 -> Write(); 
XX_xx_2_2_3 -> Write(); 
XX_xx_3_2_3 -> Write(); 
XX_xx_1_3_3 -> Write(); 
XX_xx_2_3_3 -> Write(); 
XX_xx_3_3_3 -> Write(); 

XX_xx_1_1_5 -> Write(); 
XX_xx_2_1_5 -> Write(); 
XX_xx_3_1_5 -> Write(); 
XX_xx_1_2_5 -> Write(); 
XX_xx_2_2_5 -> Write(); 
XX_xx_3_2_5 -> Write(); 
XX_xx_1_3_5 -> Write(); 
XX_xx_2_3_5 -> Write(); 
XX_xx_3_3_5 -> Write(); 

XX_xx_1_1_7 -> Write(); 
XX_xx_2_1_7 -> Write(); 
XX_xx_3_1_7 -> Write(); 
XX_xx_1_2_7 -> Write(); 
XX_xx_2_2_7 -> Write(); 
XX_xx_3_2_7 -> Write(); 
XX_xx_1_3_7 -> Write(); 
XX_xx_2_3_7 -> Write(); 
XX_xx_3_3_7 -> Write(); 

jets_Q2_3 -> Write();
jets_Q2_5 -> Write();
jets_Q2_7 -> Write();
}