<commit_before>{
   GeFiCa::PointContactRZ *detector2 = new GeFiCa::PointContactRZ(346,506);
   detector2->RUpperBound=3.45;
   detector2->RLowerBound=-3.45;
   detector2->ZUpperBound=5.05;
   detector2->PointBegin=-0.14;
   detector2->PointEnd=0.1400;

   //TF2 *im=new TF2("f","-0.19175e10-0.025e10*y");
   TF2 *im=new TF2("f","-0.318e10+0.025e10*y");
   //TF1 *im1=new TF1("f","-0.318e10+0.025e10*x",0,6.9);

   detector2->MaxIterations=1e6;
   detector2->Csor=1.994;
   detector2->V0=2500*GeFiCa::volt;
   detector2->V1=0*GeFiCa::volt;

   //TF1 *im=new TF1("","pol1",-0.318e10,0.025e10)
   detector2->Impurity="-0.318e10+0.025e10*y";//-0.01e10/GeFiCa::cm3);
   //detector2->SetImpurity(0e10/GeFiCa::cm3);
   
   //detector2->CalculateField(GeFiCa::kSOR2);
   //detector2->SaveField("point2dSOR2.root");
   //detector2->LoadField("point2dSOR23.root");
   
/*
   // calculate fields
   GeFiCa::Planar1D *detector = new GeFiCa::Planar1D(505);
   detector->UpperBound=5.05;
   detector->V1=2500*GeFiCa::volt;
   detector->V0=0*GeFiCa::volt;
   detector->SetImpurity(-0.01e10/GeFiCa::cm3);
   detector->CalculateField(GeFiCa::kAnalytic);
   detector->SaveField("planar1dTrue.root");
*/

  
   TCanvas * cvs=new TCanvas();
   gStyle->SetOptTitle(kTRUE);
   gStyle->SetPadTopMargin(0.02);
   gStyle->SetPadRightMargin(0.01);
   gStyle->SetPadLeftMargin(0.0999999999);
   gStyle->SetLabelFont(22,"XY");
   gStyle->SetLabelSize(0.05,"XY");
   gStyle->SetTitleSize(0.05,"XY");
   gStyle->SetTitleFont(22,"XY");
   gStyle->SetLegendFont(22);
  /* 
   TChain *ta = new TChain("t");
   ta->Add("planar1dTrue.root");
   ta->Draw("e1:c1");
   TGraph *ga = new TGraph(ta->GetSelectedRows(), ta->GetV2(), ta->GetV1());
*/

   // generate graphics
   TChain *tn = new TChain("t");
   //tn->Add("point2dSOR2.root");
   //tn->Draw("c2*10:p","c1<0.00&&c1>-0.05");
  // TGraph *gn = new TGraph(tn->GetSelectedRows(), tn->GetV2(), tn->GetV1());
 /* 
   TChain *ta = new TChain("t");
   ta->Add("point2dSOR2.root");
   ta->Draw("c2*10:e2","c1<0.01&&c1>-0.01","colz");
   TGraph *ga = new TGraph(ta->GetSelectedRows(), ta->GetV2(), ta->GetV1());
  
*/
   
  TTree *t = new TTree("t","t");
  t->ReadFile("/home/byron/mjd_siggen/fields/p1/evq.nob", "r:z:v");
  t->AddFriend("t2=t","point2dSOR2.root");
  t->Draw("z:(t2.p-v)","r>0.1&r<0.3","");
  //t->Draw("z:r:(t2.p-v)","","colz");
  

   TGraph *gn = new TGraph(t->GetSelectedRows(), t->GetV2(), t->GetV1());
/*
   // make final plot
   gn->SetMarkerColor(kBlue);
   gn->SetMarkerStyle(8);
   gn->SetMarkerSize(0.3);
   //ga->SetLineColor(kRed);
   //gn->GetXaxis()->SetTitle("Thickness [cm]");
   //gn->GetYaxis()->SetTitle("Potential [V]");
   //gn->SetTitle("");
   
   gn->Draw("ap");
   //ga->Draw("l");<commit_msg>3time grid<commit_after>{
   GeFiCa::PointContactRZ *detector2 = new GeFiCa::PointContactRZ(1036,506);
   detector2->RUpperBound=3.45;
   detector2->RLowerBound=-3.45;
   detector2->ZUpperBound=5.05;
   detector2->PointBegin=-0.13;
   detector2->PointEnd=0.1300;

   //TF2 *im=new TF2("f","-0.19175e10-0.025e10*y");
   TF2 *im=new TF2("f","-0.318e10+0.025e10*y");
   //TF1 *im1=new TF1("f","-0.318e10+0.025e10*x",0,6.9);

   detector2->MaxIterations=1e6;
   detector2->Precision=1e-8;
   detector2->Csor=1.996;
   detector2->V0=2500*GeFiCa::volt;
   detector2->V1=0*GeFiCa::volt;

   //TF1 *im=new TF1("","pol1",-0.318e10,0.025e10)
   detector2->Impurity="-0.318e10+0.025e10*y";//-0.01e10/GeFiCa::cm3);
   //detector2->SetImpurity(0e10/GeFiCa::cm3);
   
   detector2->CalculateField(GeFiCa::kSOR2);
   detector2->SaveField("point2dSOR2k2.root");
   //detector2->LoadField("point2dSOR23.root");
   
/*
   // calculate fields
   GeFiCa::Planar1D *detector = new GeFiCa::Planar1D(505);
   detector->UpperBound=5.05;
   detector->V1=2500*GeFiCa::volt;
   detector->V0=0*GeFiCa::volt;
   detector->SetImpurity(-0.01e10/GeFiCa::cm3);
   detector->CalculateField(GeFiCa::kAnalytic);
   detector->SaveField("planar1dTrue.root");
*/

  
   TCanvas * cvs=new TCanvas();
   gStyle->SetOptTitle(kTRUE);
   gStyle->SetPadTopMargin(0.02);
   gStyle->SetPadRightMargin(0.01);
   gStyle->SetPadLeftMargin(0.0999999999);
   gStyle->SetLabelFont(22,"XY");
   gStyle->SetLabelSize(0.05,"XY");
   gStyle->SetTitleSize(0.05,"XY");
   gStyle->SetTitleFont(22,"XY");
   gStyle->SetLegendFont(22);
  /* 
   TChain *ta = new TChain("t");
   ta->Add("planar1dTrue.root");
   ta->Draw("e1:c1");
   TGraph *ga = new TGraph(ta->GetSelectedRows(), ta->GetV2(), ta->GetV1());
*/

   // generate graphics
   TChain *tn = new TChain("t");
   //tn->Add("point2dSOR2.root");
   //tn->Draw("c2*10:p","c1<0.00&&c1>-0.05");
  // TGraph *gn = new TGraph(tn->GetSelectedRows(), tn->GetV2(), tn->GetV1());
  
   TChain *ta = new TChain("t");
   ta->Add("point2dSOR2k.root");
   //ta->Draw("c2*10:c1:p","","colz");
//   TGraph *ga = new TGraph(ta->GetSelectedRows(), ta->GetV2(), ta->GetV1());
  

   
  TTree *t = new TTree("t","t");
  t->ReadFile("/home/byron/mjd_siggen/fields/p1/evqp.nob", "r:z:v");
  t->AddFriend("t2=t","point2dSOR2k.root");
  //t->Draw("z:(t2.p-v)","z!=1&r!=1&z<34.5&r>0.1&r<0.3","");
  t->Draw("z:r:(t2.p-v)","z!=1&r!=1&z<34.5","colz");
  

 //  TGraph *gn = new TGraph(t->GetSelectedRows(), t->GetV2(), t->GetV1());
/*
   // make final plot
   gn->SetMarkerColor(kBlue);
   gn->SetMarkerStyle(8);
   gn->SetMarkerSize(0.3);
   //ga->SetLineColor(kRed);
   //gn->GetXaxis()->SetTitle("Thickness [cm]");
   //gn->GetYaxis()->SetTitle("Potential [V]");
   //gn->SetTitle("");
   
   gn->Draw("ap");
   //ga->Draw("l");
{'fraction_non_alphanumeric': 0.24400792935664084, 'fraction_numerical': 0.0893854748603352, 'mean_word_length': 4.0592525068368275, 'pattern_counts': {'":': 0, '<': 10, '<?xml version=': 0, '>': 133, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 2, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '3588216', 'n_tokens_mistral': 2656, 'n_tokens_neox': 2474, 'n_words': 245}