Throw an error when --soloFeatures Velocyto is used with --soloType S…

…martSeq .
alexdobin · Aug 15, 2020 · 8537c16 · 8537c16
1 parent 1fe083b
commit 8537c16
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/source/ParametersSolo.cpp b/source/ParametersSolo.cpp
@@ -137,6 +137,12 @@ void ParametersSolo::initialize(Parameters *pPin)
         exitWithError(errOut.str(),std::cerr, pP->inOut->logMain, EXIT_CODE_PARAMETER, *pP);
     };
 
+    if (featureYes[SoloFeatureTypes::Velocyto] && type==SoloTypes::SmartSeq) {
+        string errOut = "EXITING because of fatal PARAMETERS error: --soloFeatures Velocyto is presently not compatible with --soloType SmartSeq .\n";
+        errOut       += "SOLUTION: re-run without --soloFeatures Velocyto .";
+        exitWithError(errOut, std::cerr, pP->inOut->logMain, EXIT_CODE_PARAMETER, *pP);
+    };
+
     if (type != SoloTypes::CB_samTagOut) {//gene quantification is needed
         if (featureYes[SoloFeatureTypes::Gene]) {
             pP->quant.gene.yes=true;