Functions | |
| int | FindFileType (ANtpHeaderInfo *header) |
| void | ConstructBadRuns () |
| bool | IsNDRunGood (int, int) |
| int | ConvertBeamRunToIndex (BeamType::BeamType_t beam, int runType) |
| void | ConvertIndexToBeamRun (int index, BeamType::BeamType_t &beam, int &runType) |
| int | FindRunType (ANtpHeaderInfo *header) |
| int | FindRunType (TString fileName) |
| TString | WhichCCPIDFiles (Detector::Detector_t detType, ReleaseType::Release_t mcType, BeamType::BeamType_t beamType) |
| TString | WhichMdaCoeffHistDef (Bool_t isCoeff, Detector::Detector_t detType, ReleaseType::Release_t mcType, BeamType::BeamType_t beamType) |
| int | RecoEnergyScaleValueToKey (double track, double shower) |
| void | RecoEnergyKeyToScaleValue (int key, double &track, double &shower) |
| void | MakeDataQualityPlots (std::vector< TH1D * > &dq, TString type) |
| void | Make2DDataQualityPlots (std::vector< TH2D * > &dq2D, TString type) |
Variables | |
| const double | kBaseLineFar = 735. |
| const double | kBaseLineNear = 1. |
| const int | kHighMultClean = 0 |
| const int | kLowMultClean = 1 |
| const int | kNoClean = 2 |
| const int | kBRStopper = 0 |
| const int | kPittStopper = 1 |
| const int | kCCStopper = 2 |
| const int | kBRFiducialVolume = 0 |
| const int | kOxFiducialVolume = 1 |
| const int | kCCFiducialVolume = 2 |
| const int | kTOFiducialVolume = 3 |
| const int | kRockMuonFiducialVolume = 4 |
| const int | kBREventAlgorithm = 0 |
| const int | kTOEventAlgorithm = 1 |
| const int | kCCEventAlgorithm = 2 |
| const int | kUnkEventAlgorithm = 3 |
| const int | kOxEventAlgorithm = 4 |
| const int | kBRSnarlAlgorithm = 0 |
| const int | kCCSnarlAlgorithm = 1 |
| const int | kStandardSpillSelection = 0 |
| const int | kCCSpillSelection = 1 |
| const int | kCCFinalSelection = 0 |
| const int | kNCFinalSelection = 1 |
| const int | kCCCuts = 0 |
| const int | kNCCuts = 1 |
| const int | kTOCuts = 2 |
| const int | kOxCuts = 3 |
| const int | kNCCCFidCuts = 4 |
| const int | kNC = 0 |
| const int | kCC = 1 |
| const int | kUnknown = 2 |
| const int | kMCTotal = 3 |
| const int | kTau = 4 |
| const int | kUnknownFile = 0 |
| const int | kBeamFile = 1 |
| const int | kTauFile = 2 |
| const int | kElectronFile = 3 |
| const int | kMockFile = 4 |
| const int | kRunAll = 0 |
| const int | kRunI = 1 |
| const int | kRunII = 2 |
| const int | kRunIII = 3 |
| const int | kMaxRun = kRunII |
| const int | kBadRuns [38] |
| const TString | kRunNames [kMaxRun+1] |
| const int | kNuMuToNuMu = 0 |
| const int | kNuMuToNuTau = 1 |
| const int | kNuMuToNuE = 2 |
| const int | kNuMuToNuS = 3 |
| const int | kNuEToNuE = 4 |
| const int | k3Flavor = 0 |
| const int | k4Flavor = 1 |
| const int | k4FlavorDelta43Is0 = 2 |
| const int | k4FlavorDelta41Is0 = 3 |
| const int | k4FlavorAngles = 4 |
| const int | k4FlavorAnglesDelta43Big = 5 |
| const int | k4FlavorAnglesDelta41Is0 = 6 |
| const int | k4FlavorAnglesDelta43Is0 = 7 |
| const int | kData = 1 |
| const int | kMC = 0 |
| const int | kNCCCExtractionTRann = 0 |
| const int | kNCCCExtractionTO = 1 |
| const int | kNCCCExtractionAS = 2 |
| const int | kNCCCExtractionNS = 3 |
| const int | kNCCCExtractionRO = 4 |
| const int | kNCCCExtractionKA = 5 |
| const int | kNCCCExtractionKAD = 6 |
| const int | kNCCCExtractionRPann = 7 |
| const int | kNCCCExtractionRPann2 = 8 |
| const int | kNCCCExtractionDP = 9 |
| const int | kNCCCExtractionADM = 10 |
| const int | kNCCCExtractionTR = 11 |
| const int | kNCCCExtractionPL = 12 |
| const int | kNCCCExtractionTOm = 13 |
| const int | kNCCCExtractionUnk = 14 |
| const int | kNumExtractions = 15 |
| const TString | kExtractionNames [15] |
| const int | kExtrapolationUnk = 0 |
| const int | kExtrapolationMQ = 1 |
| const int | kExtrapolationNS = 2 |
| const int | kExtrapolationDP = 3 |
| const int | kExtrapolationRS = 4 |
| const int | kExtrapolationJK = 5 |
| const int | kExtrapolationJK_fit = 6 |
| const int | kNumExtrapolations = 7 |
| const TString | kExtrapolationNames [7] |
| const int | kNumBeamNames = 14 |
| const TString | kBeamNames [14] |
| const int | kNumRangeLowerLimitsADM = 13 |
| const double | kRangeLowerLimitsADM [13] |
| const int | kTruthEnergy = 0 |
| const int | kTruthY = 1 |
| const int | kRecoEnergy = 2 |
| const int | kRecoY = 3 |
| const int | kNumEfficiencyAndPurityBaseNames = 4 |
| const TString | kEfficiencyAndPurityBaseNames [4] |
| const int | kUS3Sqr = 0 |
| const int | kUMu3Sqr = 1 |
| const int | kDeltaMSqr = 2 |
| const int | kNumUMu3SqrBins = 100 |
| const double | kUMu3SqrStart = 1.e-3 |
| const double | kUMu3SqrEnd = 1. |
| const double | kDeltaUMu3Sqr = (kUMu3SqrEnd-kUMu3SqrStart)/(1.*kNumUMu3SqrBins) |
| const int | kNumUS3SqrBins = 100 |
| const double | kUS3SqrStart = 1.e-3 |
| const double | kUS3SqrEnd = 1. |
| const double | kDeltaUS3Sqr = (kUS3SqrEnd-kUS3SqrStart)/(1.*kNumUS3SqrBins) |
| const int | kNumDeltaMSqrBins = 300 |
| const double | kDeltaMSqrStart = 1.5 |
| const double | kDeltaMSqrEnd = 4.5 |
| const double | kDeltaDeltaMSqr = (kDeltaMSqrEnd-kDeltaMSqrStart)/(1.*kNumDeltaMSqrBins) |
| const double | k127 = 1.267 |
| const int | kNumEnergyBinsNear = 120 |
| const int | kNumEnergyBinsFar = 23 |
| const double | kEnergyBinWidthNear = 1. |
| const double | kEnergyBinWidthFar = 1. |
| const double | kMaxEnergy = 120. |
| const double | kEnergyBinsNear [kNumEnergyBinsNear+1] |
| const double | kEnergyBinsFar [kNumEnergyBinsFar+1] |
| const double | kRecoEnergyChange = 0.005 |
| const int | kNumTrackChanges = 33 |
| const double | kInitialTrackValue = 0.92 |
| const int | kNumShowerChanges = 81 |
| const double | kInitialShowerValue = 0.8 |
| const int | kma_qe = 0 |
| const int | kma_res = 1 |
| const int | kqel_fa0 = 2 |
| const int | kqel_eta = 3 |
| const int | kres_omega = 4 |
| const int | kres_z = 5 |
| const int | kma_coh = 6 |
| const int | kcoh_r0 = 7 |
| const int | kcoh_rei = 8 |
| const int | kkno_a1 = 9 |
| const int | kkno_a2 = 10 |
| const int | kkno_a3 = 11 |
| const int | kkno_a4 = 12 |
| const int | kkno_b = 13 |
| const int | kkno_r112 = 14 |
| const int | kkno_r122 = 15 |
| const int | kkno_r132 = 16 |
| const int | kkno_r142 = 17 |
| const int | kkno_r113 = 18 |
| const int | kkno_r123 = 19 |
| const int | kkno_r133 = 20 |
| const int | kkno_r143 = 21 |
| const int | kkno_r212 = 22 |
| const int | kkno_r222 = 23 |
| const int | kkno_r232 = 24 |
| const int | kkno_r242 = 25 |
| const int | kkno_r213 = 26 |
| const int | kkno_r223 = 27 |
| const int | kkno_r233 = 28 |
| const int | kkno_r243 = 29 |
| const int | kDISFACT = 30 |
| const int | kCCMA = 31 |
| const int | kkno_r112122 = 32 |
| const int | kkno_r113123 = 33 |
| const int | kkno_r212222 = 34 |
| const int | kkno_r213223 = 35 |
| const int | kTrackEnergy = 36 |
| const int | kShowerEnergy = 37 |
| const int | kFiducialVolume = 38 |
| const int | kNormalization = 39 |
| const int | kAbsoluteHadronicCalibration = 40 |
| const int | kShowerEnergyOffset = 41 |
| const int | kNCBackground = 42 |
| const int | kPIDCut = 43 |
| const int | kRelativeHadronicCalibration = 44 |
| const int | kCCBackground = 45 |
| const int | kLowCompleteness = 46 |
| const int | kSKZP = 47 |
| const int | kNCCrossSection = 48 |
| const int | kNuBarCrossSection = 49 |
| const int | kNCFarCleanNoise = 50 |
| const int | kNCFarCleanCR = 51 |
| const int | kNCNearClean = 52 |
| const int | kNCCleanRunDiff = 53 |
| const int | kNCRunDiff = 54 |
| const int | kNumNeugenParameters = 36 |
| const int | kNumCrossSectionParameters = 2 |
| const int | kNumRecoParameters = 16 |
| const int | kNumSKZPParameters = 1 |
| const int | kNumParameters = kNumNeugenParameters+kNumCrossSectionParameters+kNumRecoParameters+kNumSKZPParameters |
| const TString | kParameterNames [kNumParameters] |
| const double | kParameterDefaults [kNumParameters] |
| const double | kParameterSigmas [kNumParameters] |
| const double | kParameterDefaultsCarrot [kNumParameters] |
| const double | kParameterSigmasCarrot [kNumParameters] |
| const int | kDQVtxX = 0 |
| const int | kDQVtxY = 1 |
| const int | kDQVtxZ = 2 |
| const int | kDQEventLength = 3 |
| const int | kDQNumTracks = 4 |
| const int | kDQTrackExtension = 5 |
| const int | kDQEventsVsTime = 6 |
| const int | kDQPOTVsTime = 7 |
| const int | kDQShowerE = 8 |
| const int | kDQTrackMomentum = 9 |
| const int | kDQRadialVtx = 10 |
| const int | kDQRadialSqrVtx = 11 |
| const int | kDQEventsPerSnarl = 12 |
| const int | kDQPHFraction = 13 |
| const int | kDQDeltaTSpill = 14 |
| const int | kDQAvStripsPerPlane = 15 |
| const int | kDQTransverseRMS = 16 |
| const int | kDQTrackDCosZVtx = 17 |
| const int | kDQTrackDCosYVtx = 18 |
| const int | kDQTrackEndZ = 19 |
| const int | kDQTrackEndY = 20 |
| const int | kDQEndMetersToEdge = 21 |
| const int | kDQTotalStrips = 22 |
| const int | kDQPulseHeight = 23 |
| const int | kDQEnergy = 24 |
| const int | kDQEnergyQ1 = 25 |
| const int | kDQEnergyQ2 = 26 |
| const int | kDQEnergyQ3 = 27 |
| const int | kDQEnergyQ4 = 28 |
| const int | kDQEnergyZ1 = 29 |
| const int | kDQEnergyZ2 = 30 |
| const int | kDQEnergyZ3 = 31 |
| const int | kDQEnergyA1 = 32 |
| const int | kDQEnergyA2 = 33 |
| const int | kDQNumDists = 34 |
| const TString | kDQNames [34] |
| const TString | kDQTitles [34] |
| const int | kDQ2DEventsPerPOT = 0 |
| const int | kDQ2DEventsPerPOTLT05GeV = 1 |
| const int | kDQ2DEventsPerPOTLT1GeV = 2 |
| const int | kDQ2DEventsPerPOTGT1GeV = 3 |
| const int | kDQ2DNumDists = 4 |
| const TString | kDQ2DNames [4] |
| const TString | kDQ2DTitles [4] |
|
|
Definition at line 392 of file NCType.cxx. References mBadRuns, and mBadRunsIter. Referenced by NCAnalysisModule::NCAnalysisModule(), and NCExtrapolationModule::NCExtrapolationModule(). 00392 {
00393 vector<int> temp;
00394
00395 for( uint i=0; i < sizeof(kBadRuns)/sizeof(*kBadRuns); ++i ){
00396 mBadRuns[ kBadRuns[i] ] = temp;
00397 }
00398
00399 mBadRuns[9805].push_back(17);
00400 mBadRuns[9805].push_back(18);
00401
00402 mBadRunsIter = mBadRuns.begin();
00403
00404 return;
00405 }
|
|
||||||||||||
|
Definition at line 222 of file NCType.cxx. References BeamType::ToZarko(). Referenced by NCExtrapolation::AddEvent(), NCExtrapolationModule::AddEventsToExtrapolations(), NCExtrapolationModule::EndJob(), and NCExtrapolationModule::SetPOTValues(). 00224 {
00225 return (runType + 100*BeamType::ToZarko(beam));
00226 }
|
|
||||||||||||||||
|
Definition at line 229 of file NCType.cxx. References BeamType::FromZarko(). Referenced by NCExtrapolation::NCExtrapolation(), and NCExtrapolationModule::SetPOTValues(). 00232 {
00233 int beamNum = TMath::FloorNint(index/100.);
00234 beam = BeamType::FromZarko(beamNum);
00235 runType = index - beamNum*100;
00236
00237 return;
00238 }
|
|
|
Definition at line 289 of file NCType.cxx. References ANtpHeaderInfo::run. Referenced by NCExtrapolationModule::AddEventsToExtrapolations(). 00290 {
00291 //check if the event is from a tau file. the mc run numbers have the convention
00292 //<ikey><itgt><iflv><S><RRRR>
00293 //<ikey> = 1 near, 2 far
00294 //<itgt> = 1 in detector only, 2 in rock only, 3 detector + rock, 4 in fiducial only
00295 //<iflv> = 0 beam, 1 nue, 3 nutau, 4 numu, nue inverted, 5 far oscillated mock data
00296 //<S> = split set
00297 //<RRRR> = run number
00298
00299 int fileType = NCType::kBeamFile;
00300
00301 if(header->run > 21300000
00302 && header->run < 21400000) fileType = NCType::kTauFile;
00303
00304 //dont know what these runs are just now
00305 else if(header->run > 21400000
00306 && header->run < 21510000) fileType = NCType::kElectronFile;
00307
00308 else if(header->run > 21910000
00309 && header->run < 21910004) fileType = NCType::kMockFile;
00310
00311 return fileType;
00312 }
|
|
|
Definition at line 264 of file NCType.cxx. 00265 {
00266 //run II conditions
00267 if(fileName.Contains("2006-05")
00268 || fileName.Contains("2006-06")
00269 || fileName.Contains("2006-07")
00270 || fileName.Contains("2006-08")
00271 || fileName.Contains("2006-09")
00272 || fileName.Contains("2006-10")
00273 || fileName.Contains("2006-11")
00274 || fileName.Contains("2006-12")
00275 || fileName.Contains("2007-01")
00276 || fileName.Contains("2007-02")
00277 || fileName.Contains("2007-03")
00278 || fileName.Contains("2007-04")
00279 || fileName.Contains("2007-05")
00280 || fileName.Contains("2007-06")
00281 || fileName.Contains("2007-07")
00282 || fileName.Contains("2007-08")
00283 ) return NCType::kRunII;
00284
00285 return NCType::kRunI;
00286 }
|
|
|
Definition at line 241 of file NCType.cxx. References ANtpHeaderInfo::dataType, ANtpHeaderInfo::detector, and ANtpHeaderInfo::run. Referenced by NCExtrapolationModule::AddEventsToExtrapolations(), NCExtrapolationModule::AddFilesToChain(), NCExtrapolationModule::CombineDataQualityPlots(), NCAnalysisModule::FillDataQualityPlots(), NCExtrapolationModule::FillDataQualityPlotsSpecial(), and NCExtrapolationModule::SetPOTValues(). 00242 {
00243 //run I corresponds to anything before run 10160 in the near detector and
00244 //35724 in the far detector
00245 if(header->detector == (int)Detector::kFar
00246 && header->dataType == (int)SimFlag::kData){
00247 if(header->run > 35723){
00248 return NCType::kRunII;
00249 }
00250 }
00251 else if(header->detector == (int)Detector::kNear
00252 && header->dataType == (int)SimFlag::kData){
00253 if(header->run > 10159){
00254 return NCType::kRunII;
00255 }
00256 }
00257
00258 return NCType::kRunI;
00259 }
|
|
||||||||||||
|
Definition at line 409 of file NCType.cxx. References mBadRuns, mBadRunsIter, and run(). Referenced by NCAnalysisModule::ExtractNCCC(), NCExtrapolationModule::FinalEventCheck(), and NCAnalysisUtils::FinalEventCheck(). 00409 {
00410
00411 if( run > (--mBadRuns.end())->first ){ return true;}// run is beyond NCType::kBadRuns[] index
00412 else{
00413
00414 mBadRunsIter = mBadRuns.find( run );
00415
00416 if( mBadRunsIter == mBadRuns.end() ){ return true;}// run is not in NCType::kBadRuns[] list
00417 else if( mBadRunsIter->second.size() == 0 ){ return false;}// Whole run is bad
00418 else{
00419
00420 for( uint i=0; i < mBadRunsIter->second.size(); ++i){
00421 if( subrun == mBadRunsIter->second.at(i) ){ return false;}// subrun is bad
00422 }
00423
00424 return true;
00425
00426 }
00427 }
00428
00429 return true;
00430 }
|
|
||||||||||||
|
Definition at line 367 of file NCType.cxx. Referenced by NCAnalysisModule::NCAnalysisModule(), and NCExtrapolationModule::NCExtrapolationModule(). 00369 {
00370
00371 int binsx[4] = { 40, 40, 40, 40};
00372 int binsy[4] = { 50, 50, 50, 50};
00373 double startx[4] = {0., 0., 0., 0.};
00374 double starty[4] = {0., 0., 0., 0.};
00375 double endx[4] = {40., 40., 40., 40.};
00376 double endy[4] = {50., 50., 50., 50.};
00377
00378 for(int k = 0; k < NCType::kDQ2DNumDists; ++k){
00379 dq2D.push_back(new TH2D(NCType::kDQ2DNames[k]+type, NCType::kDQ2DTitles[k],
00380 binsx[k], startx[k], endx[k],
00381 binsy[k], starty[k], endy[k]));
00382 dq2D[k]->GetXaxis()->CenterTitle();
00383 dq2D[k]->GetYaxis()->CenterTitle();
00384 }//end loop over histograms
00385
00386 return;
00387 }
|
|
||||||||||||
|
Definition at line 315 of file NCType.cxx. Referenced by NCAnalysisModule::NCAnalysisModule(), and NCExtrapolationModule::NCExtrapolationModule(). 00317 {
00318
00319 VldTimeStamp start_ts(2005,5,1,0,0,0);
00320 VldTimeStamp end_ts(2007,5,1,0,0,0);
00321
00322 int bins[34] = { 20, 20, 32, 50, 10,
00323 50, 24, 24, 50, 100,
00324 25, 25, 10, 25, 100,
00325 10, 20, 50, 50, 32,
00326 20, 10, 100, 100, 100,
00327 100, 100, 100, 100, 100,
00328 100, 100, 100, 100};
00329 double start[34] = {-4., -4., 0., 0., 0.,
00330 -30., 0., 0., 0., -50.,
00331 0., 0., 0., 0., -50.,
00332 0., 0., -1., -1., 0.,
00333 -4., 0., 0., 0., 0.,
00334 0., 0., 0., 0., 0.,
00335 0., 0., 0., 0.};
00336 double end[34] = { 4., 4., 32., 500., 10.,
00337 20., 24., 24., 25., 50.,
00338 5., 25., 10., 1., 50.,
00339 5., 1., 1., 1., 32.,
00340 4., 4., 1000., 10., 100.,
00341 100., 100., 100., 100., 100.,
00342 100., 100., 100., 100.};
00343 TString months[25] = {"May '05", "Jun '05", "Jul '05", "Aug '05",
00344 "Sep '05", "Oct '05", "Nov '05", "Dec '05",
00345 "Jan '06", "Feb '06", "Mar '06", "Apr '06",
00346 "May '06", "Jun '06", "Jul '06", "Aug '06",
00347 "Sep '06", "Oct '06", "Nov '06", "Dec '06",
00348 "Jan '07", "Feb '07", "Mar '07", "Apr '07"};
00349
00350 for(int k = 0; k < NCType::kDQNumDists; ++k){
00351 dq.push_back(new TH1D(NCType::kDQNames[k]+type, "", bins[k], start[k], end[k]));
00352 dq[k]->SetXTitle(NCType::kDQTitles[k]);
00353 dq[k]->GetXaxis()->CenterTitle();
00354 dq[k]->SetYTitle("Events");
00355 dq[k]->GetYaxis()->CenterTitle();
00356 if(k == NCType::kDQEventsVsTime || k == NCType::kDQPOTVsTime){
00357 if(k == NCType::kDQPOTVsTime) dq[k]->SetYTitle("10^{12} POT");
00358 for(int l = 0; l < bins[k]; ++l)
00359 dq[k]->GetXaxis()->SetBinLabel(l+1, months[l]);
00360 }//end set labels for events vs time histogram
00361 }//end loop over histograms
00362
00363 return;
00364 }
|
|
||||||||||||||||
|
Definition at line 212 of file NCType.cxx. 00213 {
00214 double keyD = key*1.e-4;
00215 track = TMath::Floor(keyD)*1.e-3;
00216 shower = 10.*(keyD-track*1.e3);
00217
00218 return;
00219 }
|
|
||||||||||||
|
Definition at line 199 of file NCType.cxx. 00200 {
00201
00202 int key = 0;
00203 int trackI = TMath::Nint(track*1.e7);
00204 int showerI = TMath::Nint(shower*1.e3);
00205
00206 key = trackI + showerI;
00207
00208 return key;
00209 }
|
|
||||||||||||||||
|
Definition at line 138 of file NCType.cxx. References ReleaseType::IsCarrot(), ReleaseType::IsDaikon(), and MSG. Referenced by NCExtractionTOm::DoExtraction(), NCExtractionTO::DoExtraction(), and NCExtractionDP::FillPDFs(). 00140 {
00141 if (mcType==ReleaseType::kUnknown) {
00142 MSG("NCType",Msg::kWarning)
00143 << "Can't figure out MC version, defaulting to carrot.\n";
00144 mcType = ReleaseType::kCarrot;
00145 }
00146
00147 if(detType == Detector::kFar) {
00148 if(ReleaseType::IsCarrot(mcType)){
00149 switch (beamType){
00150 case BeamType::kL010z185i: return "dp_pdf_far_le.root"; break;
00151 case BeamType::kUnknown: return "dp_pdf_far_le.root"; break;
00152 default: return ""; break;
00153 }
00154 } else if(ReleaseType::IsDaikon(mcType)) {
00155 switch (beamType){
00156 case BeamType::kL010z185i: return "dp_pdf_far_le_cedar_daikon.root"; break;
00157 case BeamType::kUnknown: return "dp_pdf_far_le_cedar_daikon.root"; break;
00158 default: return ""; break;
00159 }
00160 } else return "";
00161 } else if(detType == Detector::kNear) {
00162 if(ReleaseType::IsCarrot(mcType)){
00163 switch (beamType){
00164 case BeamType::kL000z200i: return "dp_pdf_near_le.root"; break;
00165 case BeamType::kL010z185i: return "dp_pdf_near_le.root"; break;
00166 case BeamType::kL050z200i: return "dp_pdf_near_le.root"; break;
00167 case BeamType::kL100z200i: return "dp_pdf_near_me.root"; break;
00168 case BeamType::kL200z200i: return "dp_pdf_near_le.root"; break;
00169 case BeamType::kL250z200i: return "dp_pdf_near_he.root"; break;
00170 case BeamType::kL010z000i: return "dp_pdf_near_le.root"; break;
00171 case BeamType::kL010z170i: return "dp_pdf_near_le.root"; break;
00172 case BeamType::kL010z200i: return "dp_pdf_near_le.root"; break;
00173 case BeamType::kL010z185i_lowintensity: return "dp_pdf_near_le.root"; break;
00174 case BeamType::kL150z200i: return "dp_pdf_near_le.root"; break;
00175 case BeamType::kUnknown: return "dp_pdf_near_le.root"; break;
00176 default: return ""; break;
00177 }
00178 }else if(ReleaseType::IsDaikon(mcType)){
00179 switch (beamType){
00180 case BeamType::kL000z200i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00181 case BeamType::kL010z185i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00182 case BeamType::kL050z200i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00183 case BeamType::kL100z200i: return "dp_pdf_near_L100z200i_cedar_daikon.root"; break;
00184 case BeamType::kL200z200i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00185 case BeamType::kL250z200i: return "dp_pdf_near_L250z200i_cedar_daikon.root"; break;
00186 case BeamType::kL010z000i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00187 case BeamType::kL010z170i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00188 case BeamType::kL010z200i: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00189 case BeamType::kL010z185i_lowintensity: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00190 case BeamType::kL150z200i: return "dp_pdf_near_L150z200i_cedar_daikon.root"; break;
00191 case BeamType::kUnknown: return "dp_pdf_near_L010z185i_cedar_daikon.root"; break;
00192 default: return ""; break;
00193 }
00194 } else return "";
00195 } else return "";
00196 }
|
|
||||||||||||||||||||
|
Definition at line 20 of file NCType.cxx. References ReleaseType::IsCarrot(), ReleaseType::IsDaikon(), and MSG. Referenced by NCExtractionAS::ReadCalibInfoFromFile(). 00023 {
00024 if (mcType==ReleaseType::kUnknown) {
00025 MSG("NCType",Msg::kWarning)
00026 << "Can't figure out MC version, defaulting to carrot.\n";
00027 mcType = ReleaseType::kCarrot;
00028 }
00029
00030 if(detType == Detector::kFar) {
00031 if(ReleaseType::IsCarrot(mcType)){
00032 switch (beamType){
00033 case BeamType::kL010z185i:
00034 if (isCoeff) return "Mda_Coeff_Far_L010z185i_BirchCarrot.dat";
00035 else return "Mda_HistDef_Far_L010z185i_BirchCarrot.dat";
00036 case BeamType::kUnknown:
00037 if (isCoeff) return "Mda_Coeff_Far_L010z185i_BirchCarrot.dat";
00038 else return "Mda_HistDef_Far_L010z185i_BirchCarrot.dat";
00039 default: return "";
00040 }
00041 } else if(ReleaseType::IsDaikon(mcType)) {
00042 switch (beamType){
00043 case BeamType::kL010z185i:
00044 if (isCoeff) return "Mda_Coeff_Far_L010z185i_CedarDaikon.dat";
00045 else return "Mda_HistDef_Far_L010z185i_CedarDaikon.dat";
00046 case BeamType::kUnknown:
00047 if (isCoeff) return "Mda_Coeff_Far_L010z185i_CedarDaikon.dat";
00048 else return "Mda_HistDef_Far_L010z185i_CedarDaikon.dat";
00049 default: return "";
00050 }
00051 } else return "";
00052 } else if(detType == Detector::kNear) {
00053 if(ReleaseType::IsCarrot(mcType)){
00054 switch (beamType){
00055 case BeamType::kL000z200i:
00056 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00057 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00058 case BeamType::kL010z185i:
00059 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00060 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00061 case BeamType::kL050z200i:
00062 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00063 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00064 case BeamType::kL100z200i:
00065 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00066 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00067 case BeamType::kL200z200i:
00068 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00069 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00070 case BeamType::kL250z200i:
00071 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00072 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00073 case BeamType::kL010z000i:
00074 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00075 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00076 case BeamType::kL010z170i:
00077 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00078 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00079 case BeamType::kL010z200i:
00080 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00081 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00082 case BeamType::kL010z185i_lowintensity:
00083 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00084 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00085 case BeamType::kL150z200i:
00086 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00087 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00088 case BeamType::kUnknown:
00089 if (isCoeff) return "Mda_Coeff_Near_L010z185i_BirchCarrot.dat";
00090 else return "Mda_HistDef_Near_L010z185i_BirchCarrot.dat";
00091 default: return "";
00092 }
00093 }else if(ReleaseType::IsDaikon(mcType)){
00094 switch (beamType){
00095 case BeamType::kL000z200i:
00096 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00097 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00098 case BeamType::kL010z185i:
00099 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00100 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00101 case BeamType::kL050z200i:
00102 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00103 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00104 case BeamType::kL100z200i:
00105 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00106 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00107 case BeamType::kL200z200i:
00108 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00109 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00110 case BeamType::kL250z200i:
00111 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00112 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00113 case BeamType::kL010z000i:
00114 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00115 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00116 case BeamType::kL010z170i:
00117 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00118 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00119 case BeamType::kL010z200i:
00120 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00121 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00122 case BeamType::kL010z185i_lowintensity:
00123 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00124 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00125 case BeamType::kL150z200i:
00126 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00127 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00128 case BeamType::kUnknown:
00129 if (isCoeff) return "Mda_Coeff_Near_L010z185i_CedarDaikon.dat";
00130 else return "Mda_HistDef_Near_L010z185i_CedarDaikon.dat";
00131 default: return "";
00132 }
00133 } else return "";
00134 } else return "";
00135 }
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: { 7629, 7683, 7697, 7716, 7937, 7942,
7981, 7982, 7988, 8165, 8168, 8169,
8214, 8308, 8385, 8423, 8426, 8645,
9805, 9809, 9831, 9867, 10271, 10896,
11231, 11318, 11386, 11598, 11778, 11782,
11785, 11881, 11893, 11977, 12200, 12540,
12543, 12587}
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {"Do Not Use",
"LE10_185kA",
"ME",
"HE",
"pME",
"pHE",
"LE10_LowIntensity",
"LE10_170kA",
"LE10_200kA",
"LE10_NoHorn",
"LE10_All",
"LE10_All185",
"LE10_Far",
"All"}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {"eventsPerPOT",
"eventsPerPOTLT05GeV",
"eventsPerPOTLT1GeV",
"eventsPerPOTGT1GeV"}
|
|
|
|
|
|
Initial value: {";POT/spill (#times10^{12});Events/Spill",
";POT/spill (#times10^{12});Events/Spill",
";POT/spill (#times10^{12});Events/Spill",
";POT/spill (#times10^{12});Events/Spill"}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {"truthEnergy",
"truthY",
"recoEnergy",
"recoY"}
|
|
|
|
|
|
Initial value: { 0.0, 1.0, 2.0, 3.0, 4.0,
5.0, 6.0, 7.0, 8.0, 9.0,
10.0, 11.0, 12.0, 13.0, 14.0,
15.0, 16.0, 17.0, 18.0, 19.0,
20.0, 25.0, 30.0, kMaxEnergy}
|
|
|
Initial value: { 0.0, 1.0, 2.0, 3.0, 4.0,
5.0, 6.0, 7.0, 8.0, 9.0,
10.0, 11.0, 12.0, 13.0, 14.0,
15.0, 16.0, 17.0, 18.0, 19.0,
20.0, 21.0, 22.0, 23.0, 24.0,
25.0, 26.0, 27.0, 28.0, 29.0,
30.0, 31.0, 32.0, 33.0, 34.0,
35.0, 36.0, 37.0, 38.0, 39.0,
40.0, 41.0, 42.0, 43.0, 44.0,
45.0, 46.0, 47.0, 48.0, 49.0,
50.0, 51.0, 52.0, 53.0, 54.0,
55.0, 56.0, 57.0, 58.0, 59.0,
60.0, 61.0, 62.0, 63.0, 64.0,
65.0, 66.0, 67.0, 68.0, 69.0,
70.0, 71.0, 72.0, 73.0, 74.0,
75.0, 76.0, 77.0, 78.0, 79.0,
80.0, 81.0, 82.0, 83.0, 84.0,
85.0, 86.0, 87.0, 88.0, 89.0,
90.0, 91.0, 92.0, 93.0, 94.0,
95.0, 96.0, 97.0, 98.0, 99.0,
100.0, 101.0, 102.0, 103.0, 104.0,
105.0, 106.0, 107.0, 108.0, 109.0,
110.0, 111.0, 112.0, 113.0, 114.0,
115.0, 116.0, 117.0, 118.0, 119.0,
kMaxEnergy}
|
|
|
|
|
|
|
|
|
Initial value: {"TRann.",
"TO.",
"AS.",
"NS.",
"RO.",
"KA.",
"KAD.",
"RPann.",
"RPann2.",
"DP.",
"ADM.",
"TR.",
"PL.",
"TOm.",
"Unk."}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {"Unknown",
"Marquardt",
"BeamMatrix",
"NDFit",
"NearFit",
"Far_Near",
"Far_NearFIT"}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {0.990, 1.120, -1.267, 0.120, 1.050,
0.762, 1.000, 1.000, 0.300, 0.400,
-0.20, 0.020, 0.800, 1.420,
0.100, 0.300, 0.300, 0.100, 1.000,
1.000, 1.000, 1.000,
0.100, 0.300, 0.300, 0.100, 1.000,
1.000, 1.000, 1.000,
1.000, 1.000, 1.000, 1.00, 1.000,
1.000,
0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.,
0.,
0.,
0.000, 0.000,
0.000,0.000,0.000,0.000,0.000}
Definition at line 447 of file NCType.h. Referenced by NCExtrapolationMQ::PrepareNearDetector(). |
|
|
Initial value: {1.032, 1.032, -1.26, 0.120, 1.050,
0.75, 1.00, 1.00, 0.300, 0.500,
0.0, 0.200, 0.200, 1.30,
0.200, 0.200, 0.200, 0.200, 1.000,
1.000, 1.000, 1.000,
0.200, 0.200, 0.200, 0.200, 1.000,
1.000, 1.000, 1.000,
1.000, 1.000, 1.000, 1.00, 1.000,
1.000,
0., 0., 0., 0., 0.,
0., 0., 0., 0., 0.,
0.,
0.,
0.000, 1.000,
0.000,0.000,0.000,0.000,0.000}
|
|
|
|
|
|
Initial value: {0.1485, 0.168, 0.005, 0.007, 0.10,
0.1250, 0.3000, 0.300, 0.001, 0.02,
0.0200, 0.0200, 0.020, 0.040,
0.1000, 0.1000, 0.100, 0.100, 0.20,
0.2000, 0.2000, 0.200,
0.1000, 0.1000, 0.100, 0.100, 0.20,
0.2000, 0.2000, 0.200,
0.200, 0.15, 0.33, 0.2000, 0.33, 0.2000,
0.020, 0.11,
0.10, 0.04,
0.110, 0.150,
0.25, 0.02,
0.03, 0.150,
0.10,
1.,
1., 0.04,
1.,1.,1.,1.,1.}
Definition at line 472 of file NCType.h. Referenced by NCExtrapolationMQ::PrepareNearDetector(). |
|
|
Initial value: {0.1032, 0.1032, 0.005, 0.007, 0.10,
0.1250, 0.3000, 0.300, 0.001, 0.02,
0.0200, 0.0200, 0.020, 0.040,
0.1000, 0.1000, 0.100, 0.100, 0.20,
0.2000, 0.2000, 0.200,
0.1000, 0.1000, 0.100, 0.100, 0.20,
0.2000, 0.2000, 0.200,
0.200, 0.15, 0.33, 0.2000, 0.33, 0.2000,
0.020, 0.11, 0.10, 0.03, 0.056,
0.150, 0.50, 0.02, 0.02, 0.050,
0.219,
1.,
1.0, 0.04,
1.,1.,1.,1.,1.}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {0.0, 1.0, 1.5, 2.0, 2.5,
3.0, 3.5, 4.0, 5.0, 7.0,
10., 15, 30.}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Initial value: {"RunAll",
"RunI",
"RunII"}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.3.9.1