Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

fix: decrease CTOF time bins from 40 ps to 25 ps #237

Merged
merged 5 commits into from
Sep 6, 2024
Merged

fix: decrease CTOF time bins from 40 ps to 25 ps #237

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
03e150d
fix: decrease CTOF time bins from 40 ps to 25 ps
c-dilks Aug 26, 2024
1a34a25
fix: data type
c-dilks Aug 27, 2024
df9274d
fix: broaden fit range
c-dilks Sep 3, 2024
32c4ed3
fix: tune fit range to about 1.2ns
c-dilks Sep 4, 2024
2253cd3
fix: revert to narrow fit range
c-dilks Sep 5, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion detectors/src/main/java/org/jlab/clas/timeline/fitter/CTOFFitter.groovy
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ class CTOFFitter {
double hAmp = h1.getBinContent(h1.getMaximumBin());
double hMu = h1.getAxis().getBinCenter(h1.getMaximumBin())
double hSigma = Math.min(HistoUtil.getHistoIQR(h1, 0.4) / 2.0, 0.2)
def rangeFactor = 1.5
def rangeFactor = 1.5 // NOTE: 5.0 gives us a range closer to the calibration suite (~1.2ns)
double rangeMin = (hMu - rangeFactor*hSigma);
double rangeMax = (hMu + rangeFactor*hSigma);
f1.setRange(rangeMin, rangeMax);
Expand Down
21 changes: 10 additions & 11 deletions monitoring/src/main/java/org/jlab/clas12/monitoring/central.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,8 @@ public class central {
public String outputDir;
public int rf_large_integer;
public boolean BackToBack;
public float STT, RFT, MinCTOF,MaxCTOF, minSTT, maxSTT;
public float STT, RFT, MinCTOF, MaxCTOF, minSTT, maxSTT;
public int NbinsCTOF;
public float[] CTOF_shft;
public double rfPeriod, rfoffset1, rfoffset2;
public int e_part_ind;
Expand Down Expand Up @@ -150,11 +151,9 @@ public central(int reqrunNum, String reqOutputDir, boolean reqTimeBased) {
if(runNum>0&&runNum<3210){
minSTT = 540;maxSTT=590;
}
MinCTOF=320;MaxCTOF=400;
if(runNum>0&&runNum<3210){
MinCTOF=-5;MaxCTOF=5;
}
MinCTOF=-5;MaxCTOF=5;
NbinsCTOF = 400; // 25 ps per bin
MinCTOF = -5;
MaxCTOF = 5;
H_CVT_t_STT = new H2F("H_CVT_t_STT","H_CVT_t_STT",100,0.,maxSTT,100,0.,maxSTT+(MinCTOF+MaxCTOF)/2);
H_CVT_t_STT.setTitle("CTOF vertex t vs Pion vertex t, neg. tracks");
H_CVT_t_STT.setTitleX("Pion vertex t (ns)");
Expand All @@ -165,22 +164,22 @@ public central(int reqrunNum, String reqOutputDir, boolean reqTimeBased) {
System.out.println(e);
}
System.out.println("CTOF range "+MinCTOF+" to "+MaxCTOF);
H_CVT_t_pad = new H2F("H_CVT_t_pad","H_CVT_t_pad",50,0.5,50.5,250,MinCTOF,MaxCTOF);
H_CVT_t_pad = new H2F("H_CVT_t_pad","H_CVT_t_pad",50,0.5,50.5,NbinsCTOF,MinCTOF,MaxCTOF);
H_CVT_t_pad.setTitle("CTOF vertex t - Pion vertex t vs pad, neg. tracks");
H_CVT_t_pad.setTitleX("pad");
H_CVT_t_pad.setTitleY("CTOF vertex t - Pion vertex t (ns)");
H_CVT_t = new H1F[49];
for(int p=0;p<49;p++){
H_CVT_t[p] = new H1F(String.format("H_CVT_t_p%d",p+1),String.format("H_CVT_t_p%d",p+1),250,MinCTOF,MaxCTOF);
H_CVT_t[p] = new H1F(String.format("H_CVT_t_p%d",p+1),String.format("H_CVT_t_p%d",p+1),NbinsCTOF,MinCTOF,MaxCTOF);
H_CVT_t[p].setTitle(String.format("pad %d CTOF vertex t - STT, neg. tracks",p+1));
H_CVT_t[p].setTitleX("Delta_t (ns)");
}
H_CVT_t[48].setTitle("all pad, CTOF vertex t - STT, neg. tracks");
H_CVT_t_pos = new H1F("H_CVT_t_pos","H_CVT_t_pos",250,MinCTOF,MaxCTOF);
H_CVT_t_pos = new H1F("H_CVT_t_pos","H_CVT_t_pos",NbinsCTOF,MinCTOF,MaxCTOF);
H_CVT_t_pos.setTitle("All CTOF pads, CTOF vertex t - Pion vertex t, pos. tracks");
H_CVT_t_pos.setTitleX("CTOF vertex t - Pion vertex t (ns)");

H_CVT_t_neg = new H1F("H_CVT_t_neg","H_CVT_t_neg",250,MinCTOF,MaxCTOF);
H_CVT_t_neg = new H1F("H_CVT_t_neg","H_CVT_t_neg",NbinsCTOF,MinCTOF,MaxCTOF);
H_CVT_t_neg.setTitle("All CTOF pads, CTOF vertex t - Pion vertex t, neg. tracks");
H_CVT_t_neg.setTitleX("CTOF vertex t - Pion vertex t (ns)");

Expand All @@ -190,7 +189,7 @@ public central(int reqrunNum, String reqOutputDir, boolean reqTimeBased) {
H_CTOF_neg_mass = new H1F("H_CTOF_neg_mass","H_CTOF_neg_mass",100,-0.5,2.0);
H_CTOF_neg_mass.setTitle("neg Mass^2");
H_CTOF_neg_mass.setTitleX("M^2 (GeV^2)");
H_CTOF_vt_pim = new H2F("H_CTOF_vt_pim","H_CTOF_vt_pim",100,-rfPeriod/2,rfPeriod/2,250,-5,5);
H_CTOF_vt_pim = new H2F("H_CTOF_vt_pim","H_CTOF_vt_pim",100,-rfPeriod/2,rfPeriod/2,NbinsCTOF,MinCTOF,MaxCTOF);
H_CTOF_vt_pim.setTitle("CTOF MIP (pi-) vertex time");
H_CTOF_vt_pim.setTitleX("vertex time - RFTime (ns)");
H_CTOF_vt_pim.setTitleY("vertex time - STTime (ns)");
Expand Down