--- tof/flight/ToFNuclei/src/ToFNuclei.cpp 2008/10/13 13:41:39 1.1 +++ tof/flight/ToFNuclei/src/ToFNuclei.cpp 2009/03/02 11:42:51 1.2 @@ -1,3 +1,12 @@ +//***************************************************** +//24-feb-2009: introduced "beta_in" and "Set_Beta": +// beta_in is set to 100 by default. If it is not overwritten +// by "Set_Beta" (ToFNuclei.h), the code is used as before. +// If beta_in was set by the user using "Set_Beta", this beta +// is taken for the calculations. +//26-feb-2009: Inserted S115B calibartions from Lucia +// +//***************************************************** #include //-------------------------------------- @@ -14,6 +23,8 @@ // L2 = l2p; // + beta_in = 100.; + // if ( !L2->IsORB() ) printf(" WARNING: OrbitalInfo Tree is needed, the plugin could not work properly without it \n"); // OBT = 0; @@ -30,10 +41,12 @@ ncount_trk = 0.; charge_trk = 0.; memset (charge_layer_trk,0,6*sizeof(Float_t)); // + memset (charge_PMT_trk,0,48*sizeof(Float_t)); // memset (layerinfo_trk,0,6*sizeof(Int_t)); // ncount_std = 0.; charge_std = 0.; memset (charge_layer_std,0,6*sizeof(Float_t)); // + memset (charge_PMT_std,0,48*sizeof(Float_t)); // memset (layerinfo_std,0,6*sizeof(Int_t)); // }; @@ -83,7 +96,9 @@ newentry = true; }; // - if ( !newentry ) return; + if ( !newentry ) { + return; + } // tr = ntr; // @@ -96,7 +111,8 @@ Clear(); - // + // cout<<"beta_in = "<GetOrbitalInfo()->absTime; @@ -108,100 +124,235 @@ //if (L2->GetTrkLevel2()->GetNTracks()!=1) return; if (L2->GetTrkLevel2()->GetNTracks()==1) { +/* Float_t PMTsat[48] = { 3162.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37, 3130.05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03, 3094.09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53, 3097.11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98, 3081.37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66, 3106.33, 3094.85, 3150.85, 3118.8, 3096.24, 3118.47,3111.36, 3117.11 -} -; +}; +*/ +// new values from Napoli dec 2008 +Float_t PMTsat[48] = { +3176.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27, +3146.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14, +3106.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62, +3112.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26, +3095.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18, +3119.36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 }; + +for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.; // safety margin + + +Float_t C0t[48][9] = { +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92, +0, 1, 2, 3, 4, 5, 6, 8, 92} ; -Float_t C0t[48][9] = {0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, -0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92, 0, 1, 2, 3, 4, 5, 6, 8, 92 }; Float_t A0t[48][9] = { -0, 0.90472, 5.19581, 7.99982, 11, 13, 16, 18, 2000, 0, 0.7774, 5.15121, 7.99976, 11, 13, 16, 18, 2000, 0, 1.18414, 5.1952, 8, 11, 13, 16, 17.9999, 2000, 0, 1.0189, 5.61869, 8, 11, 13.0006, 15.9998, 18.0014, 2000, 0, 1.08725, 4.4728, 7.99979, 11, -13.0002, 16, 18, 2000, 0, 0.946462, 4.69554, 7.99994, 11, 13.0001, 16, 18, 2000, 0, 1.13176, 4.57483, 7.99999, 11, 13.0002, 15.9999, 18.0002, 2000, 0, 1.00872, 4.99316, 8, 11, 13, 16, 18, 2000, 0, 1.19891, 5.01355, 8, 11, 13.0003, 16, 18, 2000, 0, 0.5, -2, 4.5, 5.50003, 6.50028, 8.00004, 8.99989, 2000, 0, 1.16961, 4.59028, 8, 11, 13.0001, 16, 18, 2000, 0, 0.989643, 5.21958, 7.99998, 11, 13, 16, -17.9999, 2000, 0, 0.952029, 5.19598, 7.99998, 11, 13.0005, 16.0005, 18, 2000, 0, 0.736833, 4.95596, 8, 11, 13, 16, 18, 2000, 0, 1.34007, 3.74275, 8, 11, 13, 16, 18, 2000, 0, 1.47101, 4.31154, 8, 11, 13, 16, 18, -2000, 0, 1.3754, 4.83141, 8.00002, 11, 13, 16, 18, 2000, 0, 1.47712, 5.79783, 8, 11, 13, 16, 18, 2000, 0, 1.19127, 5.35796, 7.99996, 11, 13, -16, 18, 2000, 0, 1.25745, 5.98391, 8, 11, 13.0004, 16.0001, 18.0001, 2000, 0, 1.18024, 5.40107, 7.99999, 11, 13, 16, 18.0001, 2000, 0, 0.961883, -5.10846, 7.99999, 11, 13, 16, 18, 2000, 0, 1.23192, 6.21098, 8, 11, 13.0005, 16, 18.0003, 2000, 0, 1.15846, 5.23573, 7.99995, -11, 13, 16, 18.0002, 2000, 0, 1.26163, 5.5425, 8.00001, 11, 13.0002, 16.0006, 18, 2000, 0, 1.12754, 5.57064, 7.99999, 11, 13, 16, 18.0002, -2000, 0, 1.22768, 6.06283, 8, 11, 13.0002, 16, 18, 2000, 0, 1.83569, 5.87584, 8.00006, 11, 13.0005, 16, 18, -2000, 0, 1.61434, 5.56634, 7, 7.99999, 8.99999, 14.999, 16, 2000, 0, 1.61287, 5.46416, 7.99989, 10, 12, 14.9997, 16, 2000, 0, 1.49053, -6.21581, 9.00001, 12, 14, 16, 22, 2000, 0, 1.64579, 5.08498, 7.99989, 10, 12, 15, 16, 2000, 0, 1.83474, 6.29746, 9.00001, 12, 14.0001, 16.0001, 22, -2000, 0, 1.80066, 6.19725, 9, 12, 14, 16, 22, 2000, 0, 1.9025, 5.83736, 8, 10, 12, 14.9999, 16, 2000, 0, 1.90163, 5.62801, 7.99999, 10, 12, 15, 16, 2000, 0, 2.5518, -5.97325, 9.00012, 12, 15, 18, 22, 2000, 0, 2.53476, 5.75312, 9.00007, 12, 15, 18, 22, 2000, 0, 2.67784, 6.9232, 9.00005, 12, 15, 18, 22, 2000, 0, 2.62769, 7.22552, 9.00003, -12.0001, 15, 18, 22, 2000, 0, 2.58088, 6.30413, 9.0001, 12, 15, 18, 22, 2000, 0, 2.48219, 7.06085, 9.00005, 12, 15, 18, 22, 2000, 0, 3.23157, 5.20668, 9.00007, 12, 15, -18, 22, 2000, 0, 3.38244, 6.45172, 9.00006, 12.0001, 15, 18.0001, 22, 2000, 0, 3.14595, 5.29036, 9.00008, 12, 15, 18, 22.0001, 2000, 0, 3.34692, 6.17283, 9.00041, 12, -15, 18.0001, 22, 2000, 0, 3.16915, 4.25973, 9.00015, 12, 15, 18, 22, 2000, 0, 3.25809, 5.58971, 9.00068, 12, 15, 18, -22.0001, 2000 } ; +0, 0.90472, 5.19581, 7.99982, 11, 13, 16, 18, 2000, +0, 0.7774, 5.15121, 7.99976, 11, 13, 16, 18, 2000, +0, 1.18414, 5.1952, 8, 11, 13, 16, 17.9999, 2000, +0, 1.0189, 5.61869, 8, 11, 13.0006, 15.9998, 18.0014, 2000, +0, 1.08725, 4.4728, 7.99979, 11, 13.0002, 16, 18, 2000, +0, 0.946462, 4.69554, 7.99994, 11, 13.0001, 16, 18, 2000, +0, 1.13176, 4.57483, 7.99999, 11, 13.0002, 15.9999, 18.0002, 2000, +0, 1.00872, 4.99316, 8, 11, 13, 16, 18, 2000, +0, 1.19891, 5.01355, 8, 11, 13.0003, 16, 18, 2000, +0, 0.5, 2.284, 3.12, 8.795, 6.257, 7.841, 6.1743, 2000, +0, 1.16961, 4.59028, 8, 11, 13.0001, 16, 18, 2000, +0, 0.989643, 5.21958, 7.99998, 11, 13, 16, 17.9999, 2000, +0, 0.952029, 5.19598, 7.99998, 11, 13.0005, 16.0005, 18, 2000, +0, 0.736833, 4.95596, 8, 11, 13, 16, 18, 2000, +0, 1.34007, 3.74275, 8, 11, 13, 16, 18, 2000, +0, 1.47101, 4.31154, 8, 11, 13, 16, 18, 2000, +0, 1.3754, 4.83141, 8.00002, 11, 13, 16, 18, 2000, +0, 1.47712, 5.79783, 8, 11, 13, 16, 18, 2000, +0, 1.19127, 5.35796, 7.99996, 11, 13, 16, 18, 2000, +0, 1.25745, 5.98391, 8, 11, 13.0004, 16.0001, 18.0001, 2000, +0, 1.18024, 5.40107, 7.99999, 11, 13, 16, 18.0001, 2000, +0, 0.961883, 5.10846, 7.99999, 11, 13, 16, 18, 2000, +0, 1.23192, 6.21098, 8, 11, 13.0005, 16, 18.0003, 2000, +0, 1.15846, 5.23573, 7.99995, 11, 13, 16, 18.0002, 2000, +0, 1.26163, 5.5425, 8.00001, 11, 13.0002, 16.0006, 18, 2000, +0, 1.12754, 5.57064, 7.99999, 11, 13, 16, 18.0002, 2000, +0, 1.22768, 6.06283, 8, 11, 13.0002, 16, 18, 2000, +0, 1.83569, 5.87584, 8.00006, 11, 13.0005, 16, 18, 2000, +0, 1.61434, 5.56634, 7, 7.99999, 8.99999, 14.999, 16, 2000, +0, 1.61287, 5.46416, 7.99989, 10, 12, 14.9997, 16, 2000, +0, 1.49053, 6.21581, 9.00001, 12, 14, 16, 22, 2000, +0, 1.64579, 5.08498, 7.99989, 10, 12, 15, 16, 2000, +0, 1.83474, 6.29746, 9.00001, 12, 14.0001, 16.0001, 22, 2000, +0, 1.80066, 6.19725, 9, 12, 14, 16, 22, 2000, +0, 1.9025, 5.83736, 8, 10, 12, 14.9999, 16, 2000, +0, 1.90163, 5.62801, 7.99999, 10, 12, 15, 16, 2000, +0, 2.5518, 5.97325, 9.00012, 12, 15, 18, 22, 2000, +0, 2.53476, 5.75312, 9.00007, 12, 15, 18, 22, 2000, +0, 2.67784, 6.9232, 9.00005, 12, 15, 18, 22, 2000, +0, 2.62769, 7.22552, 9.00003, 12.0001, 15, 18, 22, 2000, +0, 2.58088, 6.30413, 9.0001, 12, 15, 18, 22, 2000, +0, 2.48219, 7.06085, 9.00005, 12, 15, 18, 22, 2000, +0, 3.23157, 5.20668, 9.00007, 12, 15, 18, 22, 2000, +0, 3.38244, 6.45172, 9.00006, 12.0001, 15, 18.0001, 22, 2000, +0, 3.14595, 5.29036, 9.00008, 12, 15, 18, 22.0001, 2000, +0, 3.34692, 6.17283, 9.00041, 12, 15, 18.0001, 22, 2000, +0, 3.16915, 4.25973, 9.00015, 12, 15, 18, 22, 2000, +0, 3.25809, 5.58971, 9.00068, 12, 15, 18, 22.0001, 2000} ; + Float_t A1t[48][9] = { -0, -0.148721, -3.39495, -4.31459, -0.929148, 3.54132, -0.781274, 1.31968, 0, 0, 0.300986, -2.80331, -4.23227, 0.298916, 4.10416, -5.18516, 7.03041, 0, -0, -0.590177, -3.12529, -0.44614, 0.620419, 5.19609, 9.88925, 14.9819, 0, 0, -0.173358, -3.98938, -0.616916, 0.405404, -6.25045, 10.5087, 24.7055, 0, 0, -0.32707, -1.1686, -3.46701, -0.284242, 4.09514, -0.561233, -0.166844, 0, 0, 0.0440221, -1.5147, -2.67239, --1.08848, 5.05632, 3.92546, -0.0525454, 0, 0, -0.331836, -1.25371, -1.59317, 0.654738, 4.41673, 5.99879, -14.9057, 0, 0, 0.00587183, -2.17841, -1.48496, 0.477997, 6.24062, 6.19841, 11.4798, 0, 0, -0.515989, -2.1195, -1.44376, 0.145546, 5.41926, -5.03587, 10.7126, 0, 0, 0, 0, -0.5, 0, 4.35377, 6.97709, 15.0047, 0, 0, -0.554732, -1.6306, -1.54677, -0.687092, 3.97612, 4.05449, 3.44324, 0, 0, -0.0314164, -2.83515, -1.29946, 2.31597, 7.15645, 6.95353, 12.8573, 0, 0, -0.218413, -3.04095, -1.50327, -0.291698, 5.27119, 9.06051, 18.4417, 0, 0, 0.410005, -2.45452, -0.851601, 0.206665, 6.44906, 7.53132, -10.5636, 0, 0, -1.29063, 0.463767, -0.720394, 1.27269, 3.66394, 4.52102, 7.33444, 0, 0, -1.54073, -1.60146, -1.91721, -0.640132, 3.58881, 0.561348, -7.33444, 0, 0, -0.558527, -2.43481, -1.64492, 1.51694, 6.05449, 4.67951, 5.40004, 0, 0, -1.44736, -5.265, -1.28878, -0.238869, 4.13668, 4.58264, -3.28758, 0, 0, -0.596315, -3.35555, -2.03624, -0.148223, 6.16332, 4.71164, 15.6531, 0, 0, -0.740826, -4.79843, 0.345734, -0.210722, 6.32702, 9.76698, 19.9151, 0, 0, -0.541817, --3.30557, -1.76057, -0.216039, 5.79923, 8.00931, 16.4061, 0, 0, -0.0503124, -2.66635, -1.70196, -1.07446, 3.46058, 1.45824, --2.21932, 0, 0, -0.495955, -5.78012, -0.174017, 0.675206, 6.30264, 8.30058, 20.3899, 0, 0, -0.315843, -2.8796, -3.04323, -0.681945, 5.49248, 7.42757, -17.1667, 0, 0, -0.643179, -3.5364, -0.245549, 0.694213, 5.79727, 10.2911, 9.61651, 0, 0, -0.315433, -3.56856, -1.25894, 2.05501, 6.85829, 8.74895, -14.9955, 0, 0, -0.478207, -4.82517, 0.688668, 1.40183, 6.78743, 8.02806, 2.65821, 0, 0, -1.63072, -5.05269, -0.582789, 1.22014, 6.80023, 8.8931, 5.15517, 0, 0, --1.50826, -3.68958, 0.420744, 3.25124, 4.28046, -16.9712, -12, 0, 0, -1.52021, -3.37459, -2.41086, 0.547517, -0.56013, -8.59935, -3, 0, 0, -0.998702, -4.7187, -2.33974, --0.0974348, 3.10138, 2.3293, -5.71135, 0, 0, -1.66188, -2.39145, -2.55052, 0.29005, -0.16697, -3.97018, -1.5, 0, 0, --2.04431, -4.95214, -1.51889, 1.92081, 5.92476, 8.94417, -0.172422, 0, 0, -1.94155, -4.95644, -2.47063, 0.0233995, 2.57015, 4.84173, --0.555407, 0, 0, -2.22045, -4.21875, -1.23892, 0.911188, -0.705415, -11.5648, -7, 0, 0, -2.25418, -3.71067, --1.87307, 0.560423, 0.206679, -4.27293, -3, 0, 0, -3.90179, -5.14354, -4.08454, -3.13545, -3.57868, -5.13198, -1.0019, 0, 0, -3.78789, -4.50313, --5.24307, -3.24931, -0.319336, -1.407, -3.77821, 0, 0, -4.16337, -7.12861, -2.29639, -0.786276, 2.58494, 4.16162, 7.22747, 0, 0, -4.04189, -7.82707, --1.34893, -1.05519, 4.33166, 4.87264, 12.3328, 0, 0, -3.85946, -5.71064, -4.43564, -1.53025, 1.42505, -0.0695653, 4.33418, 0, 0, -3.65794, -8.19828, --3.32738, -1.25785, -1.13124, -6.79745, -8.37113, 0, 0, -5.571, -2.49596, -2.8623, -2.74219, -1.43215, -0.359863, 0.25062, 0, 0, -6.00932, -5.60393, --2.33967, -1.53792, 3.67721, 5.49137, 8.34507, 0, 0, -5.36442, -3.19657, -2.9975, -1.73078, 0.0558693, -0.179221, 9.45231, 0, 0, -5.91858, -5.04606, -2.88549, -0.953526, 1.17816, 5.43427, 0.871257, 0, -0, -5.53832, -0.498886, -3.85789, -2.61557, -2.45863, -1.44554, -7.21927, 0, 0, -5.70059, -3.94476, -3.23927, 0.15573, 1.85295, 3.96162, 8.70362, 0 } ; +0, -0.148721, -3.39495, -4.31459, -0.929148, 3.54132, -0.781274, 1.31968, 0, +0, 0.300986, -2.80331, -4.23227, 0.298916, 4.10416, -5.18516, 7.03041, 0, +0, -0.590177, -3.12529, -0.44614, 0.620419, 5.19609, 9.88925, 14.9819, 0, +0, -0.173358, -3.98938, -0.616916, 0.405404, 6.25045, 10.5087, 24.7055, 0, +0, -0.32707, -1.1686, -3.46701, -0.284242, 4.09514, -0.561233, -0.166844, 0, +0, 0.0440221, -1.5147, -2.67239, -1.08848, 5.05632, 3.92546, -0.0525454, 0, +0, -0.331836, -1.25371, -1.59317, 0.654738, 4.41673, 5.99879, 14.9057, 0, +0, 0.00587183, -2.17841, -1.48496, 0.477997, 6.24062, 6.19841, 11.4798, 0, +0, -0.515989, -2.1195, -1.44376, 0.145546, 5.41926, 5.03587, 10.7126, 0, +0, 0, 0.0682, 3.2814, -6.1989, 8.05589, 9.75967, 23.797, 0, +0, -0.554732, -1.6306, -1.54677, 0.687092, 3.97612, 4.05449, 3.44324, 0, +0, -0.0314164, -2.83515, -1.29946, 2.31597, 7.15645, 6.95353, 12.8573, 0, +0, -0.218413, -3.04095, -1.50327, 0.291698, 5.27119, 9.06051, 18.4417, 0, +0, 0.410005, -2.45452, -0.851601, 0.206665, 6.44906, 7.53132, 10.5636, 0, +0, -1.29063, 0.463767, -0.720394, 1.27269, 3.66394, 4.52102, 7.33444, 0, +0, -1.54073, -1.60146, -1.91721, -0.640132, 3.58881, 0.561348, 7.33444, 0, +0, -0.558527, -2.43481, -1.64492, 1.51694, 6.05449, 4.67951, 5.40004, 0, +0, -1.44736, -5.265, -1.28878, -0.238869, 4.13668, 4.58264, 3.28758, 0, +0, -0.596315, -3.35555, -2.03624, -0.148223, 6.16332, 4.71164, 15.6531, 0, +0, -0.740826, -4.79843, 0.345734, -0.210722, 6.32702, 9.76698, 19.9151, 0, +0, -0.541817, -3.30557, -1.76057, -0.216039, 5.79923, 8.00931, 16.4061, 0, +0, -0.0503124, -2.66635, -1.70196, -1.07446, 3.46058, 1.45824, -2.21932, 0, +0, -0.495955, -5.78012, -0.174017, 0.675206, 6.30264, 8.30058, 20.3899, 0, +0, -0.315843, -2.8796, -3.04323, -0.681945, 5.49248, 7.42757, 17.1667, 0, +0, -0.643179, -3.5364, -0.245549, 0.694213, 5.79727, 10.2911, 9.61651, 0, +0, -0.315433, -3.56856, -1.25894, 2.05501, 6.85829, 8.74895, 14.9955, 0, +0, -0.478207, -4.82517, 0.688668, 1.40183, 6.78743, 8.02806, 2.65821, 0, +0, -1.63072, -5.05269, -0.582789, 1.22014, 6.80023, 8.8931, 5.15517, 0, +0, -1.50826, -3.68958, 0.420744, 3.25124, 4.28046, -16.9712, -12, 0, +0, -1.52021, -3.37459, -2.41086, 0.547517, -0.56013, -8.59935, -3, 0, +0, -0.998702, -4.7187, -2.33974, -0.0974348, 3.10138, 2.3293, -5.71135, 0, +0, -1.66188, -2.39145, -2.55052, 0.29005, -0.16697, -3.97018, -1.5, 0, +0, -2.04431, -4.95214, -1.51889, 1.92081, 5.92476, 8.94417, -0.172422, 0, +0, -1.94155, -4.95644, -2.47063, 0.0233995, 2.57015, 4.84173, -0.555407, 0, +0, -2.22045, -4.21875, -1.23892, 0.911188, -0.705415, -11.5648, -7, 0, +0, -2.25418, -3.71067, -1.87307, 0.560423, 0.206679, -4.27293, -3, 0, +0, -3.90179, -5.14354, -4.08454, -3.13545, -3.57868, -5.13198, -1.0019, 0, +0, -3.78789, -4.50313, -5.24307, -3.24931, -0.319336, -1.407, -3.77821, 0, +0, -4.16337, -7.12861, -2.29639, -0.786276, 2.58494, 4.16162, 7.22747, 0, +0, -4.04189, -7.82707, -1.34893, -1.05519, 4.33166, 4.87264, 12.3328, 0, +0, -3.85946, -5.71064, -4.43564, -1.53025, 1.42505, -0.0695653, 4.33418, 0, +0, -3.65794, -8.19828, -3.32738, -1.25785, -1.13124, -6.79745, -8.37113, 0, +0, -5.571, -2.49596, -2.8623, -2.74219, -1.43215, -0.359863, 0.25062, 0, +0, -6.00932, -5.60393, -2.33967, -1.53792, 3.67721, 5.49137, 8.34507, 0, +0, -5.36442, -3.19657, -2.9975, -1.73078, 0.0558693, -0.179221, 9.45231, 0, +0, -5.91858, -5.04606, -2.88549, -0.953526, 1.17816, 5.43427, 0.871257, 0, +0, -5.53832, -0.498886, -3.85789, -2.61557, -2.45863, -1.44554, -7.21927, 0, +0, -5.70059, -3.94476, -3.23927, 0.15573, 1.85295, 3.96162, 8.70362, 0 }; + Float_t A2t[48][9] = { -0, 0.447567, 2.19777, 4.56362, 2.7416, 2.05411, 10.11, 19.0312, 0, 0, 0.11417, 1.62292, 5.00684, 3, -3.00006, 17.3401, 17.8502, 0, 0, 0.613665, 1.84729, 0.717669, 3, 2.53701, 1.23897, 10.2529, 0, 0, 0.318447, 2.30647, 0.936633, 3, 3, 2.95973, 3, 0, 0, -0.398201, 0.701095, 3.91236, 3, 3, 11.5693, 17.9749, 0, 0, 0.202853, 0.806934, 3.11779, 3.57657, 1.59933, -6.71492, 23.6721, 0, 0, 0.378081, 0.60131, 1.91722, 1.97014, 3, 5.08036, 9.64776, 0, 0, 0.185602, 1.13243, 1.95726, 2.25608, 1.31382, -6.06097, 15.9104, 0, 0, 0.519026, 1.16139, 1.92471, 3, 3, 6.24901, 15.0465, 0, 0, 0, 0, 0, 0, 3, 1.24657, -0.76179, 0, 0, 0.584808, 1.01299, 1.55521, 1.3038, 3, 7.23646, 20.977, 0, 0, 0.252569, 1.60163, 1.83587, 0.630466, 0.777058, 6.0142, 15.4738, 0, 0, -0.434331, 1.84922, 1.83168, 3, 3, 3, 6.44375, 0, 0, 0.0330463, 1.50249, 1.15961, 3, 0.810556, 4.20175, -16.0172, 0, 0, 1.15272, -0.461141, 0.6358, 0.566375, 3.20343, 5.26698, 14.6336, 0, 0, 1.24677, 1.12736, 1.76074, 3, 3.77723, 10.7814, 14.6336, 0, -0, 0.551726, 1.59633, 1.18411, 0.119382, 0.843864, 6.67491, 20.4289, 0, 0, 1.17008, 3.53223, 1.7298, 3, -3.51043, 7.10543, 22.3253, 0, 0, 0.612142, 1.96606, 2.64693, 3.25016, 1.25888, 7.04821, 11.7283, 0, 0, 0.69148, 2.84536, 0.249509, 4.29848, 3, -3.58229, 9.01681, 0, 0, 0.565678, 1.96366, 2.32387, 3.15843, 1.66386, 3.07223, 7.67847, 0, 0, 0.286252, -1.61749, 2.08437, 3.96745, 3.46972, 9.29207, 23.9444, 0, 0, 0.427677, 3.58328, 0.303785, 3, 3, 4.66547, 6.77271, 0, 0, 0.367136, 1.62323, -3.53394, 3.22174, 1.19152, 2.67962, 7.37085, 0, 0, 0.571177, 2.03267, 0.714761, 2.31548, 3, 3, 23.2035, 0, -0, 0.410871, 2.02075, 2.01824, 1.2684, 1.18522, 3.21976, 12.8378, 0, 0, 0.463877, 2.86112, 0.170916, 3, 3, 5.97413, 28.4085, 0, 0, 1.30088, -3.26557, 2, 3, 3, 4.77584, 22.8794, 0, 0, 1.02941, 2.12558, 0.872765, 2.23727, 4.7981, 24.6338, 24.6338, 0, 0, 1.05783, 1.89525, 2.85289, 3.0898, -7.83407, 17.8466, 20, 0, 0, 0.87851, 2.49044, 1.91387, 2.41541, 3.8418, 9.41288, 24.0375, 0, 0, 1.18188, 1.273, 2.9275, 3.06704, 7.12642, 12.7629, 15, -0, 0, 1.53682, 2.57011, 1.06205, 0.506326, 1.68314, 3.58498, 23.1595, 0, 0, 1.46302, 2.68515, 1.99381, 2.17816, 4.41356, 6.53634, -20.7136, 0, 0, 1.4731, 2.3654, 1.53025, 2.96327, 8.70277, 21.6253, 25, 0, 0, 1.53808, 2.07089, 2.22662, 3.14766, 7.57131, 14.2471, 20, 0, 0, -2.54185, 2.94079, 2.66362, 3.52374, 6.80629, 11.3325, 14.5186, 0, 0, 2.43664, 2.49436, 4.17446, 3.77162, -3.31898, 7.43896, 18.2548, 0, 0, 2.66465, 3.9858, 1.0855, 1.91916, 1.82719, 3.62289, 11.3872, 0, 0, 2.59071, 4.39003, 0.0757547, 3, 0.158868, -3.55454, 5.71387, 0, 0, 2.48014, 3.16071, 3.1975, 1.7259, 1.53532, 6.33812, 11.2325, 0, 0, 2.35236, -4.94998, 1.78836, 1.9115, 4.74233, 13.6973, 21.965, 0, 0, 3.52215, 0.812087, 1.1768, 2.95239, 4.3202, 6.04614, 13.6556, 0, 0, 3.79446, 2.80606, -0.88371, 3, 0.614106, 2.52782, 12.4441, 0, 0, 3.41363, 1.57296, 1.59021, 2.47593, 3.41146, 6.76304, 5.59471, 0, 0, 3.72529, 2.50889, 2, 3, 4.06815, 2.78676, -19.6427, 0, 0, 3.54713, -0.296079, 2.10851, 2.42037, 5.37573, 6.78065, 22.298, 0, 0, 3.61517, 1.96909, 2, 0.249865, 2.46639, 3.42483, 8.1652, 0 } ; +0, 0.447567, 2.19777, 4.56362, 2.7416, 2.05411, 10.11, 19.0312, 0, +0, 0.11417, 1.62292, 5.00684, 3, 3.00006, 17.3401, 17.8502, 0, +0, 0.613665, 1.84729, 0.717669, 3, 2.53701, 1.23897, 10.2529, 0, +0, 0.318447, 2.30647, 0.936633, 3, 3, 2.95973, 3, 0, +0, 0.398201, 0.701095, 3.91236, 3, 3, 11.5693, 17.9749, 0, +0, 0.202853, 0.806934, 3.11779, 3.57657, 1.59933, 6.71492, 23.6721, 0, +0, 0.378081, 0.60131, 1.91722, 1.97014, 3, 5.08036, 9.64776, 0, +0, 0.185602, 1.13243, 1.95726, 2.25608, 1.31382, 6.06097, 15.9104, 0, +0, 0.519026, 1.16139, 1.92471, 3, 3, 6.24901, 15.0465, 0, +0, 0, 0, -1.972, 4.7828, -3.1154, -2.11319, -4.7497, 0, +0, 0.584808, 1.01299, 1.55521, 1.3038, 3, 7.23646, 20.977, 0, +0, 0.252569, 1.60163, 1.83587, 0.630466, 0.777058, 6.0142, 15.4738, 0, +0, 0.434331, 1.84922, 1.83168, 3, 3, 3, 6.44375, 0, +0, 0.0330463, 1.50249, 1.15961, 3, 0.810556, 4.20175, 16.0172, 0, +0, 1.15272, -0.461141, 0.6358, 0.566375, 3.20343, 5.26698, 14.6336, 0, +0, 1.24677, 1.12736, 1.76074, 3, 3.77723, 10.7814, 14.6336, 0, +0, 0.551726, 1.59633, 1.18411, 0.119382, 0.843864, 6.67491, 20.4289, 0, +0, 1.17008, 3.53223, 1.7298, 3, 3.51043, 7.10543, 22.3253, 0, +0, 0.612142, 1.96606, 2.64693, 3.25016, 1.25888, 7.04821, 11.7283, 0, +0, 0.69148, 2.84536, 0.249509, 4.29848, 3, 3.58229, 9.01681, 0, +0, 0.565678, 1.96366, 2.32387, 3.15843, 1.66386, 3.07223, 7.67847, 0, +0, 0.286252, 1.61749, 2.08437, 3.96745, 3.46972, 9.29207, 23.9444, 0, +0, 0.427677, 3.58328, 0.303785, 3, 3, 4.66547, 6.77271, 0, +0, 0.367136, 1.62323, 3.53394, 3.22174, 1.19152, 2.67962, 7.37085, 0, +0, 0.571177, 2.03267, 0.714761, 2.31548, 3, 3, 23.2035, 0, +0, 0.410871, 2.02075, 2.01824, 1.2684, 1.18522, 3.21976, 12.8378, 0, +0, 0.463877, 2.86112, 0.170916, 3, 3, 5.97413, 28.4085, 0, +0, 1.30088, 3.26557, 2, 3, 3, 4.77584, 22.8794, 0, +0, 1.02941, 2.12558, 0.872765, 2.23727, 4.7981, 24.6338, 24.6338, 0, +0, 1.05783, 1.89525, 2.85289, 3.0898, 7.83407, 17.8466, 20, 0, +0, 0.87851, 2.49044, 1.91387, 2.41541, 3.8418, 9.41288, 24.0375, 0, +0, 1.18188, 1.273, 2.9275, 3.06704, 7.12642, 12.7629, 15, 0, +0, 1.53682, 2.57011, 1.06205, 0.506326, 1.68314, 3.58498, 23.1595, 0, +0, 1.46302, 2.68515, 1.99381, 2.17816, 4.41356, 6.53634, 20.7136, 0, +0, 1.4731, 2.3654, 1.53025, 2.96327, 8.70277, 21.6253, 25, 0, +0, 1.53808, 2.07089, 2.22662, 3.14766, 7.57131, 14.2471, 20, 0, +0, 2.54185, 2.94079, 2.66362, 3.52374, 6.80629, 11.3325, 14.5186, 0, +0, 2.43664, 2.49436, 4.17446, 3.77162, 3.31898, 7.43896, 18.2548, 0, +0, 2.66465, 3.9858, 1.0855, 1.91916, 1.82719, 3.62289, 11.3872, 0, +0, 2.59071, 4.39003, 0.0757547, 3, 0.158868, 3.55454, 5.71387, 0, +0, 2.48014, 3.16071, 3.1975, 1.7259, 1.53532, 6.33812, 11.2325, 0, +0, 2.35236, 4.94998, 1.78836, 1.9115, 4.74233, 13.6973, 21.965, 0, +0, 3.52215, 0.812087, 1.1768, 2.95239, 4.3202, 6.04614, 13.6556, 0, +0, 3.79446, 2.80606, 0.88371, 3, 0.614106, 2.52782, 12.4441, 0, +0, 3.41363, 1.57296, 1.59021, 2.47593, 3.41146, 6.76304, 5.59471, 0, +0, 3.72529, 2.50889, 2, 3, 4.06815, 2.78676, 19.6427, 0, +0, 3.54713, -0.296079, 2.10851, 2.42037, 5.37573, 6.78065, 22.298, 0, +0, 3.61517, 1.96909, 2, 0.249865, 2.46639, 3.42483, 8.1652, 0} ; + // S115B failure: 1158700000 // 2006-09-19 21:06:40 UTC -if (time > 1158700000.) { // first workaround: just use calibration of PMT S115A +if (time < 1158700000.) { // first workaround: just use calibration of PMT S115A for ( Int_t jj=0; jj<9; jj++ ){ A0t[9][jj] = A0t[8][jj] ; A1t[9][jj] = A1t[8][jj] ; @@ -223,24 +374,28 @@ Int_t iz; Float_t x1[9],y1[9]; Int_t n1 = 9; -Float_t Layercharge[6] = {1000., 1000., 1000., 1000., 1000., 1000. }; + +Float_t Layercharge[6] = {1005., 1005., 1005., 1005., 1005., 1005. }; Int_t Layerinfo[6] = {0, 0, 0, 0, 0, 0 }; +Float_t adc[48],adc_c[48],PMTcharge_arr[48]; + //------------------------------------------------------------------ PamTrack *track = L2->GetTrack(0); Float_t beta_mean = fabs(track->GetToFTrack()->beta[12]); + +// New!! If "beta_in" != 100. then overwrite beta_mean with beta_in +if (beta_in != 100.) beta_mean=beta_in; + Float_t def = track->GetTrkTrack()->al[4]; Float_t rig = 1./def; Int_t ncount=0; Float_t charge_chi=1000.; -if ((beta_mean > 0.05) && (beta_mean < 1.5)) { - - //------------------------------------------------------------- //------------------ ToF Charge ------------------------------ //------------------------------------------------------------- @@ -248,12 +403,24 @@ // fill adc and tdc Int_t ich,jj,kk; - Float_t adc[48],adc_c[48],PMTcharge_arr[48]; for (ich=0; ich<48; ich++) adc[ich] = 4095.; for (ich=0; ich<48; ich++) adc_c[ich] = 1000.; for (ich=0; ich<48; ich++) PMTcharge_arr[ich] = 1000.; + + +if ((beta_mean < 0.05) || (beta_mean > 2.0)) { + for (Int_t jj=0; jj<48;jj++) charge_PMT_trk[jj] = 1001.; + for (Int_t jj=0; jj<6;jj++) charge_layer_trk[jj] = 1001.; + for (Int_t jj=0; jj<6;jj++) layerinfo_trk[jj] = 1001.; + charge_trk = 1001.; + ncount_trk = 0; + } + + +if ((beta_mean > 0.05) && (beta_mean < 2.0)) { + Float_t betahelp = pow(beta_mean, 1.2); //======================================================================= @@ -274,10 +441,12 @@ Int_t adcflag = track->GetToFTrack()->adcflag[ipmt]; L2->GetToFLevel2()->GetPMTIndex(pmtadc,jj,kk); if (adcflag==0) adc_c[pmtadc] = dEdx; - if (adc[pmtadc] > PMTsat[pmtadc]) adc_c[pmtadc] = 1000.; - cout<PMTsat[pmtadc]) && (adc[pmtadc]<4095)) adc_c[pmtadc] = 1000.; if (adc[pmtadc] < PMTsat[pmtadc]) { //---------------------------------------------------------------------------------------- @@ -296,7 +465,9 @@ spl2->Delete(); PMTcharge_arr[pmtadc] = charge; } // pmtadc != ... + //---------------------------------------------------------------------------------------- + if ((pmtadc==1) || (pmtadc==9) || (pmtadc==28) || (pmtadc==29) || (pmtadc==34)) { for ( Int_t jj=0; jj<8; jj++ ){ yl = A0t[pmtadc][jj]+A1t[pmtadc][jj]*beta_mean+A2t[pmtadc][jj]*beta_mean*beta_mean; @@ -312,6 +483,7 @@ } // if pmtadc ==.... //---------------------------------------------------------------------------------------- } // adc < sat + } // ipmt... //==================================================================== @@ -339,7 +511,8 @@ for (Int_t ilay=0; ilay<6; ilay ++) { Int_t jj = PaddleIdOfTrack[ilay] ; -Float_t xhelp = 1000.; +//Float_t xhelp = 1000.; +Float_t xhelp = 1005.; if (jj == -1) Layerinfo[ilay] = -1; @@ -353,15 +526,11 @@ if (adc[ih+1] == 4095) Layerinfo[ilay] = Layerinfo[ilay] + 1; if ((adc[ih] > PMTsat[ih]) && (adc[ih] < 4095)) Layerinfo[ilay] = Layerinfo[ilay] + 10; if ((adc[ih+1] > PMTsat[ih+1]) && (adc[ih+1] < 4095)) Layerinfo[ilay] = Layerinfo[ilay] + 10; -// cout< "< 1158700000.) { // first workaround: just use calibration of PMT S115A +if (time < 1158700000.) { // first workaround: just use calibration of PMT S115A for ( Int_t jj=0; jj<9; jj++ ){ A0t[9][jj] = A0t[8][jj] ; A1t[9][jj] = A1t[8][jj] ; @@ -641,28 +826,43 @@ Float_t x1[9],y1[9]; Int_t n1 = 9; -Float_t Layercharge[6] = {1000., 1000., 1000., 1000., 1000., 1000. }; +Float_t Layercharge[6] = {1005., 1005., 1005., 1005., 1005., 1005. }; + Int_t Layerinfo[6] = {0, 0, 0, 0, 0, 0 }; +Float_t adc[48],adc_c[48],PMTcharge_arr[48]; + +Float_t beta_mean = fabs(track->beta[12]); - Float_t beta_mean = fabs(track->beta[12]); +// New!! If "beta_in" != 100. then overwrite beta_mean with beta_in +if (beta_in != 100.) beta_mean=beta_in; Int_t ncount=0; Float_t charge_chi=1000.; - if ((beta_mean > 0.05) && (beta_mean < 1.5)) { //------------------------------------------------------------- // fill adc and tdc Int_t ich,jj,kk; - Float_t adc[48],adc_c[48],PMTcharge_arr[48]; for (ich=0; ich<48; ich++) adc[ich] = 4095.; for (ich=0; ich<48; ich++) adc_c[ich] = 1000.; for (ich=0; ich<48; ich++) PMTcharge_arr[ich] = 1000.; + +if ((beta_mean < 0.05) || (beta_mean > 2.0)) { + for (Int_t jj=0; jj<48;jj++) charge_PMT_std[jj] = 1001.; + for (Int_t jj=0; jj<6;jj++) charge_layer_std[jj] = 1001.; + for (Int_t jj=0; jj<6;jj++) layerinfo_std[jj] = 1001.; + charge_std = 1001.; + ncount_std = 0; + } + + + if ((beta_mean > 0.05) && (beta_mean < 2.0)) { + Float_t betahelp = pow(beta_mean, 1.2); //======================================================================= @@ -685,7 +885,11 @@ Int_t adcflag = track->adcflag[ipmt]; L2->GetToFLevel2()->GetPMTIndex(pmtadc,jj,kk); if (adcflag==0) adc_c[pmtadc] = dEdx; - if (adc[pmtadc] > PMTsat[pmtadc]) adc_c[pmtadc] = 1000.; + + if (adc[pmtadc] == 4095) adc_c[pmtadc] = 0.; + if (adc[pmtadc] == 4095) PMTcharge_arr[pmtadc] = 0.; + + if ((adc[pmtadc]>PMTsat[pmtadc]) && (adc[pmtadc]<4095)) adc_c[pmtadc] = 1000.; if (adc[pmtadc] < PMTsat[pmtadc]) { @@ -768,7 +972,8 @@ for (Int_t ilay=0; ilay<6; ilay ++) { Int_t jj = PaddleIdOfTrack[ilay] ; -Float_t xhelp = 1000.; +//Float_t xhelp = 1000.; +Float_t xhelp = 1005.; if (jj == -1) Layerinfo[ilay] = -1; @@ -801,7 +1006,8 @@ for (Int_t ii=0; ii<6; ii++) { Int_t ic = 0; Float_t corr[4] = {1., 0., 0., 0.}; -for (Int_t iii=1; iii<9; iii++) if (fabs(Layercharge[ii] - iii)<0.50) ic = iii; +//for (Int_t iii=1; iii<9; iii++) if (fabs(Layercharge[ii] - iii)<0.50) ic = iii; +for (Int_t iii=1; iii<10; iii++) if (fabs(Layercharge[ii] - iii)<0.50) ic = iii; if ((ic>0)&&(ic<=5)) for (Int_t iii=0; iii<4; iii++) corr[iii] = corr_a[ic-1][iii]; if (ic >5) for (Int_t iii=0; iii<4; iii++) corr[iii] = corr_a[5][iii]; @@ -810,15 +1016,20 @@ if ((ic>0)&&(ic<=5)) c_corr = c_corr/(1.*ic) ; if (ic>5) c_corr = c_corr/6. ; +//cout< "; Layercharge[ii] = Layercharge[ii]/c_corr; +//cout<0) charge_mean = charge_mean/(1.*icount) ; +if (icount>0) { + +//for (Int_t ii=0; ii<6; ii++) cout< charge_v; + Float_t charge_median = 0.; + + //////// get the median charge //////////// + for(Int_t ii=0;ii<6;ii++) { + if (Layercharge[ii] < 1000) charge_v.push_back(Layercharge[ii]); + } + + std::sort(charge_v.begin(),charge_v.end()); + charge_median = charge_v[(int)(charge_v.size()/2+charge_v.size()%2)-1]; + +//cout<0 //--- Calculate charge for the second time, use residuals of the single //--- measurements to get a chi2 value @@ -873,13 +1106,17 @@ // charge_tof = charge_chi; - } // 0.05 < fabs( beta) < 1.5 +// changed 30-01-2009 WM + for (Int_t jj=0; jj<48;jj++) charge_PMT_std[jj] = PMTcharge_arr[jj]; for (Int_t jj=0; jj<6;jj++) charge_layer_std[jj] = Layercharge[jj]; for (Int_t jj=0; jj<6;jj++) layerinfo_std[jj] = Layerinfo[jj]; charge_std = charge_chi; ncount_std = ncount; + } // 0.05 < fabs( beta) < 2.0 + + } // standalone track found }