Recipe for creating a ZDC 1n peak in the 2015 PbPb data

To fix a wiring mixup, EM4- and HAD1- (the labels of the leafs) need to be interchanged: negEM4cFtsz->negHD1cFtsz, negHD1cFtsz->negEM4cFtsz.

Signal and background can be taken from timeslices 4,5 and 2,7, respectively. A slightly better choice is to align it channel by channel, according to the following table (-1 is "background", 1 is "signal", 0 is "not used")

  double fCmultiplier[NSIDE][NTYPE][NCH][NTS] =
{
  { //positive side
    { //HD
      {0, 0,-1, 1, 1, 0, 0,-1,0,0}, //channel 1
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 2
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 3
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 4
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 5
    },
    { //EM
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 1
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 2
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 3
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 4
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 5
    },
  },
  { //negative side
    { //HD
      {0, 0,-1, 1, 1, 0, 0,-1,0,0}, //channel 1
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 2
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 3
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 4
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 5
    },
    { //EM
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 1
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 2
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 3
      {0,-1, 0, 0, 1, 1, 0,-1,0,0}, //channel 4
      {0, 0,-1, 0, 1, 1, 0,-1,0,0}, //channel 5
    },
  }
};

Relative gains can be given as follows:

double fCgain[NSIDE][NTYPE][NCH] =
{
 {
   {0.500, 3.090, 2.423, 1.000, 0.000}, //HD+
   {1.000, 0.957, 1.000, 0.957, 1.000}, //EM+
 },
 {
   {0.417, 5.160, 4.071, 1.000, 0.000}, //HD-
   {0.330, 0.454, 1.875, 1.436, 1.000}, //EM-
 }
};

The relative gain of EM to HAD is 0.3 for the positive side and 0.5 for the negative side.

It should be checked (for HAD1, HAD1&2 or HAD1,2&3) that the signal maximum is in timeslice 4 or 5.

With this recipe, the single neutron peak becomes visible.

Example code for creating a single neutron peak in 2015 PbPb data can be downloaded here:

Here are some plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3:

2D plots of EM vs EMHADsum, for ZDC- and ZDC+, cent≥0 and cent≥100:

neg_centmin000_HD123check.png neg_centmin100_HD123check.png

pos_centmin000_HD123check.png pos_centmin100_HD123check.png

1D plots of EMHADsum:

neg1d_centmin000_HD123check.png neg1d_centmin100_HD123check.png

pos1d_centmin000_HD123check.png pos1d_centmin100_HD123check.png

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng neg1d_centmin000_HD123check.png r1 manage 18.2 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng neg1d_centmin100_HD123check.png r1 manage 18.5 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng neg_centmin000_HD123check.png r1 manage 68.6 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng neg_centmin100_HD123check.png r1 manage 55.6 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng pos1d_centmin000_HD123check.png r1 manage 18.3 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng pos1d_centmin100_HD123check.png r1 manage 18.5 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng pos_centmin000_HD123check.png r1 manage 57.4 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
PNGpng pos_centmin100_HD123check.png r1 manage 47.8 K 2017-06-01 - 12:11 MateCsanad Plots showing approximate single neutron peaks, using the signal maximum condition in HAD1,2&3
C source code filec singleneutron.C r1 manage 15.0 K 2017-06-01 - 11:28 MateCsanad Example code for creating a single neutron peak in 2015 PbPb data
Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2017-06-01 - MateCsanad
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback