// FileName: reco_MC_far_Chopped_R1_18_2.C // Created by: Joshua Boehm // Last Modified: May 17, 2006 // Correction to settings // Usage // cp /reco_MC_far_Chopped_R1_18_2.C reco_MC.C // loon -bq "reco_MC.C()" ///////////////////////////////////////////////////////// // Description of this // // This file will produce an MC sim of a Spill stream // // This is only meant to do MC FarDet processing // // for nue, numu, or nutau in R1-18-2. // // Processing occurs through a single path "Reco" // // Each output stream requires its own processing // // branch and so there are small offshoot branches // // Uses: CandSubShower Module // // CandChop Module // // Demuxer: Alt/Cambridge // // // // Output: // // Output Files: // // CandS.root - candidate output file // // ntupleStS.root - NtpSt ntuple output // // ntupleStS.sub.root - trimmed NtpSt // // // ///////////////////////////////////////////////////////// class JobC; int LoadLibraries(); int GetRunNumber(); void SetCalibrator(); void SetCandOutputWithBNtp(JobCModule &jcm, TString file); void SetCandOutputNoBNtp(JobCModule &jcm, TString file); void SetNtpStOutputWithBNtp(JobCModule &jcm, TString file); void SetNtpStOutputNoBNtp(JobCModule &jcm, TString file); void SetMSGLevels(JobC &jc); void reco_MC() { LoadLibraries(); JobC jc; //Create path - first the reroot job jc.Path.Create("Reco", "RerootToTruthModule::Get " "PhotonTransport::Get " "DetSim::Get " ); //And Now the Normal Reco chain jc.Path.Create("Reco1", "NoiseFilterModule::Ana " "RecordSetupModule::Get " "DigitListModule::Get " "DigitListModule::Reco " "ChopModule::Reco " "BiggestChopModule::Reco " "DeMuxDigitListModule::Reco " "StripSRListModule::Reco " "SliceSRListModule::Reco " "TrackSRListModule::Reco " "FitTrackSRListModule::Reco " "ClusterSRListModule::Reco " "SubShowerSRListModule::Reco " "ShowerSRListModule::Reco " "EventSRListModule::Reco " "RecordSetupModule::Reco " "Output::Put"); //Attach the normal reco chain after the reroot path jc.Path.Attach("Reco","Reco1"); //Input Parameters jc.Input.Set("Format=reroot"); jc.Input.Set("Streams=DaqSnarl"); //takes reroot data as input // Set up RerootExodus to build the event as MC (instead of Reroot) RerootExodus::SetVldSimFlag(SimFlag::kMC); jc.Path("Reco").SetAllFilters(1); //turn on reroot filters (Don't know what they are though) jc.Path("Reco1").SetAllFilters(0); //turn off the reconstruction filters ////////////////////////////////////// //Configure Reconstruction Software for Far Det MC "Spill" data . // For a full description of the different possible module options see the // Document: ConfigOptions.txt //Settings for PhotonTransport Module JobCModule& photon=jc.Path("Reco").Mod("PhotonTransport"); photon.Set("GeVPerMip=0.001737");//was 0.001730 //end of photon settings // This setting was introduced by Jeff Hartnell on 11/4/2005 jc.Path("Reco1").Mod("ClusterSRListModule").Set("ClusterListAlgConfig=FarBeam"); jc.Path("Reco1").Mod("TrackSRListModule").Set("TrackListAlgConfig=FarBeam"); jc.Path("Reco1").Mod("FitTrackSRListModule").Set("FitTrackListAlgConfig=FarBeam"); jc.Path("Reco1").Mod("EventSRListModule").Set("EventListAlgConfig=FarBeam"); //Set up the conditions for the CandSubShower Package jc.Path("Reco1").Mod("SubShowerSRListModule").Set("SubShowerSRListAlgConfig=default"); jc.Path("Reco1").Mod("ShowerSRListModule").Set("ShowerListAlgorithm=AlgShowerSSList"); jc.Path("Reco1").Mod("ShowerSRListModule").Set("ShowerListAlgConfig=default"); jc.Path("Reco1").Mod("ShowerSRListModule").Set("ListOut=CandShowerSRList"); //End CandSubShower Settings jc.Path("Reco1").Mod("DeMuxDigitListModule").Set("SwitchPersToTemp=1"); jc.Path("Reco1").Mod("DeMuxDigitListModule").Set("NameListOut=altdemux"); jc.Path("Reco1").Mod("DeMuxDigitListModule").Set("DeMuxDigitListAlgConfig=devel"); jc.Path("Reco1").Mod("StripSRListModule").Set("ListIn=altdemux"); //CandChop Package Settings jc.Path("Reco1").Mod("ChopModule").Set("ChopAlgorithm=AlgChopListFar"); jc.Path("Reco1").Mod("BiggestChopModule").Set("MinEnergy=0.0"); //No cut. jc.Path("Reco1").Mod("BiggestChopModule").Set("OmitLiEvents=1"); //kill LI jc.Path("Reco1").Mod("BiggestChopModule").Set("LiVetoWindowLow=-1.0e-6"); jc.Path("Reco1").Mod("BiggestChopModule").Set("LiVetoWindowLow=31.0e-6"); jc.Path("Reco1").Mod("SliceSRListModule").Set("SliceListAlgConfig=FarBeam"); // Settings to use CandChop instead of CandSlice (still required in the path) AlgFactory &af = AlgFactory::GetInstance(); AlgHandle alh = af.GetAlgHandle("AlgSliceSRList","FarBeam"); AlgConfig &alg = alh.GetAlgConfig(); alg.Set("PassAll",1); // Run NO Slicer ////////////////////////////////////// // Set the simulated date of the events as to allow for consistent calibration // RerootExodus::SetOverrideVldTimeStamp(); // Set Ugli to use only database. UgliLoanPool::SetAlwaysUseDbi(true); // End of Configuration settings // //////////////////////////////////////// //////////////////////////////////////////////////////////////// //Display the Configuration of modules jc.Path("Reco").Mod("RerootToTruthModule").Report(); jc.Path("Reco").Mod("PhotonTransport").Report(); jc.Path("Reco").Mod("DetSim").Report(); jc.Path("Reco1").Mod("DeMuxDigitListModule").Report(); jc.Path("Reco1").Mod("StripSRListModule").Report(); jc.Path("Reco1").Mod("ShowerSRListModule").Report(); jc.Path("Reco1").Mod("TrackSRListModule").Report(); jc.Path("Reco1").Mod("FitTrackSRListModule").Report(); ////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////// ///////////////////////////////////////////////// // Setting the Output Paths // // - This script produces the 3 output // // streams, cand, snts, sntp // ///////////////////////////////////////////////// //Set Candidate Output to CandA.root (all of the snarls included in output) //this will become the .cand. file SetCandOutputNoBNtp(jc.Path("Reco1").Mod("Output"), "CandS.root"); //Ntuple record has its own output file so needs its own output module jc.Path.Create("NtpSR", "NtpStModule::Get " "NtpSRModule::Reco " "NtpMCModule::Reco " //Necessary to write the MC info tree "NtpTHModule::Reco " //Necessary to write the MC th tree "Output::Put "); jc.Path("NtpSR").Mod("NtpMCModule").Set("UseStandard=1"); // use NtpSt jc.Path("NtpSR").Mod("NtpTHModule").Set("UseStandard=1"); // use NtpSt jc.Path("NtpSR").Mod("NtpSRModule").Set("UseStandard=1"); // use NtpSt SetNtpStOutputNoBNtp(jc.Path("NtpSR").Mod("Output"), "ntupleStS.root"); jc.Path.Attach("Reco1","NtpSR"); // Ntuple abridged record jc.Path.Create("NtpSRFilter", "NtpStFilterModule::Reco " "Output::Put "); SetNtpStOutputNoBNtp(jc.Path("NtpSRFilter").Mod("Output"), "ntupleStS.sub.root"); jc.Path.Attach("Reco1", "NtpSRFilter"); ////////////////////////////////////////////////////// /// End of Output Settings //// ////////////////////////////////////////////////////// SetMSGLevels(jc); //Other Run Options // jc.Path("All").Run(100); //Run over the first 100 snarls // jc.Input.Next(10); //Skip to Snarl 10 jc.Path("Reco").Run(); jc.Path("Reco").Report(); //Get Message Statistics jc.Msg.Stats(); //Output the message statistics to see what is generating output } ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// // Everything from this point forward are sub-functions // int LoadLibraries() { // Link dynamic libraries gSystem->Load("libDataUtil"); gSystem->Load("libRecoBase"); gSystem->Load("libNoiseFilter"); gSystem->Load("libFilterDigitSR"); gSystem->Load("libBField"); gSystem->Load("libNumericalMethods"); gSystem->Load("libSwimmer"); gSystem->Load("libDeMux"); gSystem->Load("libCandChop"); gSystem->Load("libCandStripSR"); gSystem->Load("libCandSliceSR"); gSystem->Load("libCandTrackSR"); gSystem->Load("libCandClusterSR"); gSystem->Load("libCandSubShowerSR"); gSystem->Load("libCandShowerSR"); gSystem->Load("libCandFitTrackSR"); gSystem->Load("libCandEventSR"); gSystem->Load("libVertexFinder"); gSystem->Load("libTimeCalibratorSR"); gSystem->Load("libAstroUtil"); gSystem->Load("libCandNtupleSR"); gSystem->Load("libMCNtuple"); gSystem->Load("libTruthHelperNtuple"); gSystem->Load("libStandardNtuple"); gSystem->Load("libPhotonTransport"); gSystem->Load("libDetSim"); gSystem->Load("libMCNtupleModule"); gSystem->Load("libCandNtupleSRModule"); gSystem->Load("libTruthHelperNtupleModule"); gSystem->Load("libStandardNtupleModule"); gSystem->Load("libCandShield"); gSystem->Load("libAltDeMux"); gSystem->Load("libFiltration"); } int GetRunNumber(){ //////////////////////////////////////////////////////// //Get the run number from the Input file name //Allows us to set the name of the dst file JobCEnv& jce = JobCEnv::Instance(); string filename = jce.GetFileName(0); string rnumstring = filename.substr(filename.find_last_of("_")-8,8); int run = atoi(rnumstring.c_str()); // string srnstring = filename.substr(filename.find_last_of("_")+1,4); return run; } void SetCandOutputWithBNtp(JobCModule &jcm, TString file){ jcm.Cmd("DefineStream Config ConfigRecord"); jcm.Cmd("DefineStream NtpBDLite NtpBDLiteRecord"); jcm.Set("Streams=NtpBDLite,Cand,Config,DaqSnarl"); jcm.Set("FileName="+file); } void SetCandOutputNoBNtp(JobCModule &jcm, TString file){ jcm.Cmd("DefineStream Config ConfigRecord"); jcm.Set("Streams=Cand,Config,DaqSnarl"); jcm.Set("FileName="+file); } void SetNtpStOutputWithBNtp(JobCModule &jcm, TString file){ jcm.Cmd("DefineStream NtpSt NtpStRecord"); jcm.Cmd("DefineStream NtpBDLite NtpBDLiteRecord"); jcm.Set("Streams=NtpBDLite,NtpSt"); jcm.Set("FileName="+file); } void SetNtpStOutputNoBNtp(JobCModule &jcm, TString file){ jcm.Cmd("DefineStream NtpSt NtpStRecord"); jcm.Set("Streams=NtpSt"); jcm.Set("FileName="+file); } void SetMSGLevels(JobC &jc){ jc.Msg.SetLevel("FitTrackSR","Warning"); jc.Msg.SetLevel("NoiseFilter","Fatal"); jc.Msg.SetLevel("Cand","Error"); jc.Msg.SetLevel("Calibrator","Fatal"); jc.Msg.SetLevel("Dbi","Fatal"); jc.Msg.SetLevel("SigCor Calibrator","Fatal"); jc.Msg.SetLevel("Time Calibrator","Fatal"); jc.Msg.SetLevel("MuonCalibrator","Fatal"); jc.Msg.SetLevel("PE Calibrator","Fatal"); jc.Msg.SetLevel("MapperCalibrator","Fatal"); jc.Msg.SetLevel("DetSim","Error"); jc.Msg.SetLevel("Reroot","Error"); jc.Msg.SetLevel("Photon","Warning"); jc.Msg.SetLevel("Per","Fatal"); jc.Msg.SetLevel("Io", "Fatal"); jc.Msg.SetLevel("DeMuxDigitListModule","Error"); jc.Msg.SetLevel("Plex","Error"); jc.Msg.SetLevel("AltDeMuxModule","Error"); jc.Msg.SetLevel("AltDeMux","Fatal"); jc.Msg.SetLevel("AlgAltDeMux","Fatal"); jc.Msg.SetLevel("AlgShowerSS", "Warning"); jc.Msg.SetLevel("SubShowerSR","Error"); jc.Msg.SetLevel("AlgFitTrackSR", "Warning"); }